Ab initio molecular dynamics simulations of excited hydrogen halides and methyl halides

Using excited-state ab initio molecular dynamics simulations employing the complete-active-space self-consistent-field approach, we study the mechanism of photodissociation in terms of time evolution of structure, kinetic energy, charges and potential energy for the first excited state of hydrogen halides and methyl halides. Although the hydrogen halides and methyl halides are similar in the photodissociation mechanism, their dynamics are slightly different. The presence of the methyl group causes delay in photodissociation as compared to hydrogen halides.

Structure, dynamics, and interactions of jacalin. Insights from molecular dynamics simulations examined in conjunction with results of X-ray studies

Molecular dynamics simulations have been carried out on all the jacalin-carbohydrate complexes of known structure, models of unliganded molecules derived from the complexes and also models of relevant complexes where X-ray structures are not available. Results of the simulations and the available crystal structures involving jacalin permit delineation of the relatively rigid and flexible regions of the molecule and the dynamical variability of the hydrogen bonds involved in stabilizing the structure. Local flexibility appears to be related to solvent accessibility. Hydrogen bonds involving side chains and water bridges involving buried water molecules appear to be important in the stabilization of loop structures. The lectin-carbohydrate interactions observed in crystal structures, the average parameters pertaining to them derived from simulations, energetic contribution of the stacking residue estimated from quantum mechanical calculations, and the scatter of the locations of carbohydrate and carbohydrate-binding residues are consistent with the known thermodynamic parameters of jacalin-carbohydrate interactions. The simulations, along with X-ray results, provide a fuller picture of carbohydrate binding by jacalin than provided by crystallographic analysis alone. The simulations confirm that in the unliganded structures water molecules tend to occupy the positions occupied by carbohydrate oxygens in the lectin-carbohydrate complexes. Population distributions in simulations of the free lectin, the ligands, and the complexes indicate a combination of conformational selection and induced fit. Proteins 2009; 77:760-777.

Molecular Characterization of Sclerospora graminicola, the Incitant of Pearl Millet Downy Mildew Using ISSR Markers

The DNA polymorphism among 22 isolates of Sclerospora graminicola, the causal agent of downy mildew disease of pearl millet was assessed using 20 inter simple sequence repeats (ISSR) primers. The objective of the study was to examine the effectiveness of using ISSR markers for unravelling the extent and pattern of genetic diversity in 22 S. graminicola isolates collected from different host cultivars in different states of India. The 19 functional ISSR primers generated 410 polymorphic bands and revealed 89% polymorphism and were able to distinguish all the 22 isolates. Polymorphic bands used to construct an unweighted pair group method of averages (UPGMA) dendrogram based on Jaccard's co-efficient of similarity and principal coordinate analysis resulted in the formation of four major clusters of 22 isolates. The standardized Nei genetic distance among the 22 isolates ranged from 0.0050 to 0.0206. The UPGMA clustering using the standardized genetic distance matrix resulted in the identification of four clusters of the 22 isolates with bootstrap values ranging from 15 to 100. The 3D-scale data supported the UPGMA results, which resulted into four clusters amounting to 70% variation among each other. However, comparing the two methods show that sub clustering by dendrogram and multi dimensional scaling plot is slightly different. All the S. graminicola isolates had distinct ISSR genotypes and cluster analysis origin. The results of ISSR fingerprints revealed significant level of genetic diversity among the isolates and that ISSR markers could be a powerful tool for fingerprinting and diversity analysis in fungal pathogens.

Tryst with a molecular and supramolecular oddity: an anti-Furst-Plattner epoxide opening product as a designer additive for accessing an elusive polymorph

Additive induced polymorphism of a conformationally locked tetraacetate 3 in presence of its diastereomer 4 is described. The ester 3 was specially crafted on a trans-decalin backbone to relegate the O-H center dot center dot center dot O H-bond donors to the molecular interior and have the peripheral H-bond acceptors in 1,3-syndiaxial relationship. The supramolecular assembly of 3 was destined to evolve along two mutually exclusive pathways, namely one, which employs intermolecular O-H center dot center dot center dot O H-bonds (pathway 1) and the other that sacrifices these for intramolecular O-H center dot center dot center dot O H-bonds and settles for a crystal packing dictated by weak intermolecular interactions alone (pathway 2). Exploiting the similarity between the self-assemblies of 4 and the two recently reported dimorphs of 3, the ester 3 has been stimulated to follow the elusive non-hierarchical pathway 2 through preferential inhibition of pathway 1. Interestingly, the inhibitor 4 was obtained serendipitously en route 3 via an apparent breakdown of Furst-Plattner rule.

A comparative molecular dynamics study of diffusion of n-decane and 3-methyl pentane in Y zeolite

Molecular dynamics simulations are reported on the structure and dynamics of n-decane and 3-methylpentane in zeolite NaY. We have calculated several properties such as the center of mass-center of mass rdf, the end-end distance distribution, bond angle distribution and dihedral angle distribution. We have also analysed trajectory to obtain diffusivity and velocity autocorrelation function (VACF). Surprisingly, the diffusivity of 3-methylpentane which is having larger cross-section perpendicular to the long molecular axis is higher than n-decane at 300 K. Activation energies have been obtained from simulations performed at 200 K, 300 K, 350 K, 400 K and 450 K in the NVE ensemble. These results can be understood in terms of the previously known levitation effect. Arrhenious plot has higher value of slope for n-decane (5 center dot 9 kJ/mol) than 3-methylpentane (3 center dot 7 kJ/mol) in agreement with the prediction of levitation effect.

Self-Assembly of a Pd-6-Molecular Double-Square and a Cu-3-Trigonalbipyramidal Cage via a New Tripodal Flexible Ligand

A new tripodal flexible ligand (L) containing pyrazolyl functionality has been prepared and successfully used to obtain a pd(6) (1) molecular double-square and a cu(3) trigonalbipyramidal cage (2), where complex 1 represents the first example of a double-square obtained using a flexible tripodal ligand.

Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.

Conformational study of valinomycin: a molecular dynamics approach

Valinomycin is a highly flexible cyclic dodecadepsipeptide that transports ions across membranes. Such a flexibility in the conformation is required for its biological function since it has to encounter a variety of environments and liganding state. Exploration of conformational space of this molecule is therefore important and is one of the objectives of the present study that has been carried out by means of high temperature Molecular Dynamics. Further, the stability of the known bracelet-like structure of the uncomplexed valinomycin and the inherent flexibility around this structure has been investigated. The uncomplexed form of valinomycin has been simulated at 75–100 K for 1 ns in order to elucidate the average conformational properties. An alanine-analog of valinomycin has been simulated under identical conditions in order to evaluate the effect of sidechain on the conformational properties, The studies confirm the effect of sidechain on conformational equilibrium.

The Role of UPF0157 in the Folding of M-tuberculosis Dephosphocoenzyme A Kinase and the Regulation of the Latter by CTP

Background: Targeting the biosynthetic pathway of Coenzyme A (CoA) for drug development will compromise multiple cellular functions of the tubercular pathogen simultaneously. Structural divergence in the organization of the penultimate and final enzymes of CoA biosynthesis in the host and pathogen and the differences in their regulation mark out the final enzyme, dephosphocoenzyme A kinase (CoaE) as a potential drug target. Methodology/Principal Findings: We report here a complete biochemical and biophysical characterization of the M. tuberculosis CoaE, an enzyme essential for the pathogen's survival, elucidating for the first time the interactions of a dephosphocoenzyme A kinase with its substrates, dephosphocoenzyme A and ATP; its product, CoA and an intrinsic yet novel inhibitor, CTP, which helps modulate the enzyme's kinetic capabilities providing interesting insights into the regulation of CoaE activity. We show that the mycobacterial enzyme is almost 21 times more catalytically proficient than its counterparts in other prokaryotes. ITC measurements illustrate that the enzyme follows an ordered mechanism of substrate addition with DCoA as the leading substrate and ATP following in tow. Kinetic and ITC experiments demonstrate that though CTP binds strongly to the enzyme, it is unable to participate in DCoA phosphorylation. We report that CTP actually inhibits the enzyme by decreasing its Vmax. Not surprisingly, a structural homology search for the modeled mycobacterial CoaE picks up cytidylmonophosphate kinases, deoxycytidine kinases, and cytidylate kinases as close homologs. Docking of DCoA and CTP to CoaE shows that both ligands bind at the same site, their interactions being stabilized by 26 and 28 hydrogen bonds respectively. We have also assigned a role for the universal Unknown Protein Family 0157 (UPF0157) domain in the mycobacterial CoaE in the proper folding of the full length enzyme. Conclusions/Significance: In view of the evidence presented, it is imperative to assign a greater role to the last enzyme of Coenzyme A biosynthesis in metabolite flow regulation through this critical biosynthetic pathway.

Cellular and molecular regulation of mammalian blastocyst hatching

In mammals including humans, failure in blastocyst hatching and implantation leads to early embryonic loss and infertility. Prior to implantation, the blastocyst must hatch out of its acellular glycoprotein coat, the zona pellucida (ZP). The phenomenon of blastocyst hatching is believed to be regulated by (i) dynamic cellular components such as actin-based trophectodermal projections (TEPs), and (ii) a variety of autocrine and paracrine molecules such as growth factors, cytokines and proteases. The spatio-temporal regulation of zona lysis by blastocyst-derived cellular and molecular signaling factors is being keenly investigated. Our studies show that hamster blastocyst hatching is acelerated by growth factors such as heparin binding-epidermal growth factor and leukemia inhibitory factor and that embryo-derived, cysteine proteases including cathepsins are responsible for blastocyst hatching. Additionally, we believe that cyclooxygenase-generated prostaglandins, estradiol-17 beta mediated estrogen receptor-alpha signaling and possibly NF kappa B could be involved in peri-hatching development. Moreover, we show that TEPs are intimately involved with lysing ZP and that the TEPs potentially enrich and harbor hatching-enabling factors. These observations provide new insights into our understanding of the key cellular and molecular regulators involved in the phenomenon of mammalian blastocyst hatching, which is essential for the establishment of early pregnancy.

Ligand Specificity of Group I Biotin Protein Ligase of Mycobacterium tuberculosis

Background: Fatty acids are indispensable constituents of mycolic acids that impart toughness & permeability barrier to the cell envelope of M. tuberculosis. Biotin is an essential co-factor for acetyl-CoA carboxylase (ACC) the enzyme involved in the synthesis of malonyl-CoA, a committed precursor, needed for fatty acid synthesis. Biotin carboxyl carrier protein (BCCP) provides the co-factor for catalytic activity of ACC. Methodology/Principal Findings: BPL/BirA (Biotin Protein Ligase), and its substrate, biotin carboxyl carrier protein (BCCP) of Mycobacterium tuberculosis (Mt) were cloned and expressed in E. coli BL21. In contrast to EcBirA and PhBPL, the similar to 29.5 kDa MtBPL exists as a monomer in native, biotin and bio-5'AMP liganded forms. This was confirmed by molecular weigt profiling by gel filtration on Superdex S-200 and Dynamic Light Scattering (DLS). Computational docking of biotin and bio-5'AMP to MtBPL show that adenylation alters the contact residues for biotin. MtBPL forms 11 H-bonds with biotin, relative to 35 with bio-5'AMP. Docking simulations also suggest that bio-5'AMP hydrogen bonds to the conserved `GRGRRG' sequence but not biotin. The enzyme catalyzed transfer of biotin to BCCP was confirmed by incorporation of radioactive biotin and by Avidin blot. The K-m for BCCP was similar to 5.2 mu M and similar to 420 nM for biotin. MtBPL has low affinity (K-b = 1.06 x 10(-6) M) for biotin relative to EcBirA but their K-m are almost comparable suggesting that while the major function of MtBPL is biotinylation of BCCP, tight binding of biotin/bio-5'AMP by EcBirA is channeled for its repressor activity. Conclusions/Significance: These studies thus open up avenues for understanding the unique features of MtBPL and the role it plays in biotin utilization in M. tuberculosis.

Quantum Interference in Coherent Molecular Conductance

Coherent electronic transport through individual molecules is crucially sensitive to quantum interference. We investigate the zero-bias and zero-temperature conductance through pi-conjugated annulene molecules weakly coupled to two leads for different source-drain configurations, finding an important reduction for certain transmission channels and for particular geometries as a consequence of destructive quantum interference between states with definite momenta. When translational symmetry is broken by an external perturbation we find an abrupt increase of the conductance through those channels. Previous studies concentrated on the effect at the Fermi energy, where this effect is very small. By analyzing the effect of symmetry breaking on the main transmission channels we find a much larger response thus leading to the possibility of a larger switching of the conductance through single molecules.

Molecular characterisation of ABC transporter type FtsE and FtsX proteins of Mycobacterium tuberculosis.

Elicitation of drug resistance and various survival strategies inside host macrophages have been the hallmarks of Mycobacterium tuberculosis as a successful pathogen. ATP Binding Cassette (ABC) transporter type proteins are known to be involved in the efflux of drugs in bacterial and mammalian systems. FtsE, an ABC transporter type protein, in association with the integral membrane protein FtsX, is involved in the assembly of potassium ion transport proteins and probably of cell division proteins as well, both of which being relevant to tubercle bacillus. In this study, we cloned ftsE gene of M. tuberculosis, overexpressed and purified. The recombinant MtFtsE-6xHis protein and the native MtFtsE protein were found localized on the membrane of E. coli and M. tuberculosis cells, respectively. MtFtsE-6xHis protein showed ATP binding in vitro, for which the K42 residue in the Walker A motif was found essential. While MtFtsE-6xHis protein could partially complement growth defect of E. coli ftsE temperature-sensitive strain MFT1181, co-expression of MtFtsE and MtFtsX efficiently complemented the growth defect, indicating that the MtFtsE and MtFtsX proteins might be performing an associated function. MtFtsE and MtFtsX-6xHis proteins were found to exist as a complex on the membrane of E. coli cells co-expressing the two proteins.

The crystal and molecular structure of benzyloxycarbonyl-a-amino-isobutyryl-L-prolyl-methylamide. The observation of an X2-Pro3 Type III b-turn

The crystal and molecular structure of N-benzyloxycarbonyl-a-aminoisobutyryl-L-prolyl methylamide, the amino terminal dipeptide fragment of alamethicin, has been determined using direct methods. The compound crystallizes in the orthorhombic system with the space group P212-21. Cell dimensions are a = 7.705 A, b = 11.365 A, and c = 21.904 A. The structure has been refined using conventional procedures to a final R factor of 0.054. The molecular structure possesses a 4 - 1 intramolecular N-H--0 hydrogen bond formed between the CO group of the urethane moiety and the NH group of the methylamide function. The peptide backbone adopts the type 111 P-turn conformation, with 42 = -51.0°, +* = -39.7",&j = -65.0', $3 = -25.4'. An unusual feature is the occurrence of the proline residue at position 3 of the P-turn. The observed structure supports the view that Aib residues initiate the formation of type 111 @-turn conformations. The pyrrolidine ring is puckered in Cy-exo fashion.

An Insight to the Dynamics of Conserved Water Molecular Triad in IMPDH II (Human): Recognition of Cofactor and Substrate to Catalytic Arg 322

Inosine 5' monophosphate dehydrogenase (IMPDH II) is a key enzyme involved in the de novo biosynthesis pathway of purine nucleotides and is also considered to be an excellent target for cancer inhibitor design. The conserve R 322 residue (in human) is thought to play some role in the recognition of inhibitor and cofactor through the catalytic D 364 and N 303. The 15 ns simulation and the water dynamics of the three different PDB structures (1B3O, 1NF7, and 1NFB) of human IMPDH by CHARMM force field have clearly indicated the involvement of three conserved water molecules (W-L, W-M, and W-C) in the recognition of catalytic residues (R 322, D 364, and N 303) to inhibitor and cofactor. Both the guanidine nitrogen atoms (NH1 and NH 2) of the R 322 have anchored the di- and mono-nucleotide (cofactor and inhibitor) binding domains via the conserved W-C and W-L water molecules. Another conserved water molecule W-M seems to bridge the two domains including the R 322 and also the W-C and W-L through seven centers H-bonding coordination. The conserved water molecular triad (W-C - W-M - W-L) in the protein complex may thought to play some important role in the recognition of inhibitor and cofactor to the protein through R 322 residue.