Study of the influence of temperature on the dynamics of the catalytic cleft in 1,3-1,4-β-glucanase by molecular dynamics simulations

The dependence of some molecular motions in the enzyme 1,3-1,4-β-glucanase from Bacillus licheniformis on temperature changes and the role of the calcium ion in them were explored. For this purpose, two molecular dynamics simulated trajectories along 4 ns at low (300 K) and high (325 K) temperatures were generated by the GROMOS96 package. Several structural and thermodynamic parameters were calculated, including entropy values, solvation energies, and essential dynamics (ED). In addition, thermoinactivation experiments to study the influence of the calcium ion and some residues on the activity were conducted. The results showed the release of the calcium ion, which, in turn, significantly affected the movements of loops 1, 2, and 3, as shown by essential dynamics. These movements differ at low and high temperatures and affect dramatically the activity of the enzyme, as observed by thermoinactivation studies. The first two authors contributed equally to this work     

Reverse-docking study of the TADDOL-catalyzed asymmetric hetero-Diels–Alder reaction

The reverse-docking of a TADDOL catalyst to rigid transition-state (TS) representations of an asymmetric hetero-Diels–Alder reaction is described. The resulting docking poses represent a simplified geometric model of the TS for the catalyzed reaction. The conformational space of the catalyst in proximity to the catalyst-free TS models is sampled stochastically and the energetically favored poses are subjected to a clustering procedure to highlight structural attributes compatible with organocatalysis. Each pose is scored and ranked based on its molecular-mechanics docking energy. The reverse-docking procedure reveals a clear energetic trend in favor of the experimentally preferred product enantiomers. Analysis of the best poses suggests a geometric model that is consistent with principles of molecular recognition, catalysis, and experimental data.      

Conformational stability and three-dimensional model of the δ-opioid pharmacophore for the extended antiparallel dimer structure of Met-enkephalin in water

The conformational stability of the extended antiparallel dimer structure of Met-enkephalin in water was analyzed by examining the hydration structure of enkephalin using molecular dynamics simulations. The result shows that, despite of the hydrophicility of the terminal atoms in the pentapeptide, the main contributor for the stability of the dimer in water is the four intermolecular hydrogen bonds between the Gly² and Phe⁴ groups. The three-dimensional model of the δ-opioid pharmacophore for this dimer structure was also established. Such a model was demonstrated to match the δ-opioid pharmacophore query derived from the non-peptides SIOM, TAN-67, and OMI perfectly. This result thus strongly supports the assumption that the dimer structure of Met-enkephalin is a possible δ-receptor binding conformation. Figure Schematic model of the extended antiparallel dimer structure of Met-enkephalin     

The effect of salt concentration on DNA conformation transition: a molecular-dynamics study

We performed three 3-ns molecular dynamics simulations of d(CGCGAATTCGCG)₂ using the AMBER 8 package to determine the effect of salt concentration on DNA conformational transitions. All the simulations were started with A-DNA, with different salt concentrations, and converged with B-DNA with similar conformational parameters. However, the dynamic processes of the three MD simulations were very different. We found that the conformation transition was slow in the solution with higher salt concentration. To determine the cause of this retardation, we performed three additional 1.5-ns simulations starting with B-DNA and with the salt concentrations corresponding to the simulations mentioned above. However, astonishingly, there was no delayed conformation evolution found in any of the three simulations. Thus, our simulation conclusion is that higher salt concentrations slows the A → B conformation transition, but have no effect on the final stable structure. Figure A-DNA and B-DNA. (a) is the canonical A-DNA, and (b) is the canonical B-DNA. Looking from the central major groove     

Cold-active enzymes studied by comparative molecular dynamics simulation

Enzymes from cold-adapted species are significantly more active at low temperatures, even those close to zero Celsius, but the rationale of this adaptation is complex and relatively poorly understood. It is commonly stated that there is a relationship between the flexibility of an enzyme and its catalytic activity at low temperature. This paper gives the results of a study using molecular dynamics simulations performed for five pairs of enzymes, each pair comprising a cold-active enzyme plus its mesophilic or thermophilic counterpart. The enzyme pairs included α-amylase, citrate synthase, malate dehydrogenase, alkaline protease and xylanase. Numerous sites with elevated flexibility were observed in all enzymes; however, differences in flexibilities were not striking. Nevertheless, amino acid residues common in both enzymes of a pair (not present in insertions of a structure alignment) are generally more flexible in the cold-active enzymes. The further application of principle component analysis to the protein dynamics revealed that there are differences in the rate and/or extent of opening and closing of the active sites. The results indicate that protein dynamics play an important role in catalytic processes where structural rearrangements, such as those required for active site access by substrate, are involved. They also support the notion that cold adaptation may have evolved by selective changes in regions of enzyme structure rather than in global change to the whole protein. Figure Collective motions in C^α atoms of the active site of cold-active xylanase Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

A theoretical investigation of cytotoxic activity of celastroid triterpenoids

Quantum chemical calculations at the B3LYP/6-31G* level of theory have been carried out on 20 celastroid triterpenoids to obtain a set of molecular electronic properties and to correlate these with cytotoxic activities. The cytotoxic activities of these compounds can be roughly correlated with electronic effects related to nucleophilic addition to C(6) of the compounds: The energies of the frontier molecular orbitals (E _HOMO and E _LUMO), the HOMO-LUMO energy gap, the dipole moment, the charge on C(6), and the electrophilicity on C(6). Figure LUMO of Pristimerin.     

Molecular dynamics studies of a hexameric purine nucleoside phosphorylase

Purine nucleoside phosphorylase (PNP) (EC.2.4.2.1) is an enzyme that catalyzes the cleavage of N-ribosidic bonds of the purine ribonucleosides and 2-deoxyribonucleosides in the presence of inorganic orthophosphate as a second substrate. This enzyme is involved in purine-salvage pathway and has been proposed as a promising target for design and development of antimalarial and antibacterial drugs. Recent elucidation of the three-dimensional structure of PNP by X-ray protein crystallography left open the possibility of structure-based virtual screening initiatives in combination with molecular dynamics simulations focused on identification of potential new antimalarial drugs. Most of the previously published molecular dynamics simulations of PNP were carried out on human PNP, a trimeric PNP. The present article describes for the first time molecular dynamics simulations of hexameric PNP from Plasmodium falciparum (PfPNP). Two systems were simulated in the present work, PfPNP in ligand free form, and in complex with immucillin and sulfate. Based on the dynamical behavior of both systems the main results related to structural stability and protein-drug interactions are discussed.      

Molecular modeling of the three-dimensional structure of GLP-1R and its interactions with several agonists

Glucagon-like peptide-1 receptor (GLP-1R) is a promising molecular target for developing drugs treating type 2 diabetes. We have predicted the complete three-dimensional structure of GLP-1R and the binding modes of several GLP-1R agonists, including GLP-1, Boc5, and Cpd1, through a combination of homology modeling, molecular docking, and long-time molecular dynamics simulation on a lipid bilayer. Our model can reasonably interpret the results of a number of mutation experiments regarding GLP-1R as well as the successful modification to GLP-1 by Liraglutide. Our model is also validated by a recently revealed crystal structure of the extracellular domain of GLP-1R. An activation mechanism of GLP-1R agonists is proposed based on the principal component analysis and normal mode analysis on our predicted GLP-1R structure. Before the complete structure of GLP-1R is determined through experimental means, our model may serve as a valuable reference for characterizing the interactions between GLP-1R and its agonists. Figure Comparison of our predicted model of rGLP-1R (left) with the recently revealed crystal structure of hGLP-1R (right)     

A disulfide linked model of the complement protein C8γ complexed with C8α indel peptide

In a recent study C8γ (complement protein) with Cys40Ala substitution and a C8α derived peptide with Cys164Ala substitution were co-crystallized and their binding mode was revealed. Computer modeling and molecular dynamics simulations were performed in order to check the hypothesis that the residues Ala164 of C8α and Ala40 of C8γ occupied the right position if cysteine residues were in their place for disulfide bonding. Substitution of these two alanine residues with cysteine along with disulfide bond creation via molecular modeling and subsequent molecular dynamics simulation of the complex corroborated the hypothesis, which was also confirmed from recent crystallographic data. Average RMSD between backbone atoms of the indel peptide during the MD trajectory in comparison with the corresponding sequence of crystal structure of the C8α/γ complex was found only 0.085 nm. Figure Modeling the C*y/α comlexation. Ribbon representation of the C8y complexed with C8α indel peptide initial (green/cyan) X-ray structure and the final MD conformation (magenta/orange) after imposing the disulfide link. Average RMSD between backbone atoms of the indel peptide during MD trajectory in comparison with the corresponding sequence of crystal structure of the C8α/y complex was found only 0.085nm. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evolutionarily evolved discriminators in the 3-TPR domain of the Toc64 family involved in protein translocation at the outer membrane of chloroplasts and mitochondria

Transport of polypeptides across membranes is a general and essential cellular process utilised by molecular machines. At least one component of these complexes contains a domain composed of three tetratricopeptide repeat (3-TPR) motifs. We have focussed on the receptor Toc64 to elucidate the evolved functional specifications of its 3-TPR domain. Toc64 is a component of the Toc core complex and functionally replaces Tom70 at the outer membrane of mitochondria in plants. Its 3-TPR domain recognises the conserved C-terminus of precursor-bound chaperones. We built homology models of the 3-TPR domain of chloroplastic Toc64 from different species and of the mitochondrial isoform from Arabidopsis. Guided by modelling, we identified residues essential for functional discrimination of the differently located isoforms to be located almost exclusively on the convex surface of the 3-TPR domain. The only exception is at568Ser/ps557Met, which is positioned in the ligand-binding groove. The functional implications of the homology models are discussed. Figure Motion contained within the 2nd eigenvector of the Calpha covariance matrix of the 3-TPR domain of atToc64-V indicated by a porcupine plot Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Computational models for the shuttling motion of the macrocycle in rotaxane-based molecular switches

The shuttling motion of a macrocycle in rotaxane-based molecular switching devices has been studied by computational density functional methods. In the test case, energy profiles corresponding to the dethreading process of different types of guest molecules in a cyclobis(paraquat-p-phenylene) host verified the experimental preference of the tetrathiafulvalene recognition site over the dioxynaphthalene site in a Stoddart-Heath type molecular device. Furthermore, modification of the redox state of either the macrocycle or the guest molecule resulted in considerable changes in the computational energy profiles, which can be utilized in explaining the behavior of the host-guest system. In order to study the effect of chemical oxidation/reduction in the guest molecule, we have investigated a prototypical shaft including two octahedral ruthenium complexes linked by a conjugated C₁₄ carbon chain, where the shuttling motion can be triggered by changing the electronic environment of the active complexes with ligand exchange reactions. The computational results also indicated effective communication between the macrocycle and the conjugated carbon chain, therefore showing the importance of non-covalent host-guest interactions in the control of the motion. Figure A computational model for a chemically controlled molecular switch     

Interaction of copper organometallic precursors with barrier layers of Ti, Ta and W and their nitrides: a first-principles molecular dynamics study

Processes for the deposition of copper films on transition metal barrier layers by means CVD using organometallic precursors are often found to lead to poor adhesion characteristics of the grown film. By means of first-principles molecular dynamics simulations, we show that the source of the problem is the strong reactivity of the surfaces toward the precursors, which decompose spontaneously upon contact with the surface leading to contamination of the interface. Our simulations consider Ti, Ta, and W as barrier layers, and Cu(hfac)-(tmvs) as precursor. In contrast, we show that surfaces of these metals properly passivated with nitrogen, in such a way that only N atoms are exposed on the surface, are much less active and do not lead to decomposition of the precursor. We propose this passivation procedure as a practical solution to the adhesion problem. Figure CupraSelect on the WN (100) surface     

Study of a ligand complexed with Cdk2/Cdk4 by computer simulation

Cyclin-dependent kinases (Cdks) play important roles in the regulation of the cell cycle. Their inhibitors have entered clinical trials to treat cancer. Very recently, Davis et al. (Nat Struct Biol 9:745–749, 2002) have found a ligand NU6102, which has a high affinity with cyclin-dependent kinase 2 (K _i =6 nM) but a low affinity with cyclin-dependent kinase 4 (K _i =1,600 nM). To understand the selectivity, we use homology modeling, molecular docking, molecular dynamics and free-energy calculations to analyze the interactions. A rational 3D model of the Cdk4–NU6102 complex is built. Asp86 is a key residue that recognizes NU6102 more effectively with Cdk2 rather than Cdk4. Good binding free energies are obtained. Energetic analysis reveals that van der Waals interaction and nonpolar contributions to solvent are favorable in the formation of complexes and the sulfonamide group of the ligand plays a crucial role for binding selectivity between Cdk2 and Cdk4. Figure Two-dimensional representative for the interacting model of NU6102 complexed with the Cdk4 from a predicted structure by LIGPLOT.      

Homology modeling of a novel epoxide hydrolase (EH) from Aspergillus niger SQ-6: structure-activity relationship in expoxides inhibiting EH activity

The 3D structure of a novel epoxide hydrolase from Aspergillus niger SQ-6 (sqEH) was constructed by using homology modeling and molecular dynamics simulations. Based on the 3D model, Asp191, His369 and Glu343 were predicted as catalytic triad. The putative active pocket is a hydrophobic environment and is rich in some important non—polar residues (Pro318, Trp282, Pro319, Pro317 and Phe242). Using three sets of epoxide inhibitors for docking study, the interaction energies of sqEH with each inhibitor are consistent with their inhibitory effects in previous experiments. Moreover, a critical water molecule which closes to the His369 was identified to be an ideal position for the hydrolysis step of the reaction. Two tyrosine residues (Tyr249 and Tyr312) are able to form hydrogen bonds with the epoxide oxygen atom to maintain the initial binding and positioning of the substrate in the active pocket. These docked complex models can well interpret the substrate specificity of sqEH, which could be relevant for the structural—based design of specific epoxide inhibitors. Figure       

Oxidation mechanism in the metabolism of (S)-N-[1-(3-morpholin-4-ylphenyl)ethyl]-3-phenylacrylamide on oxyferryl active site in CYP3A4 Cytochrome: DFT modeling

The metabolism mechanism of (S)-N-[1-(3-morpholin-4ylphenyl)ethyl]-3-phenylacrylamide, mediated by CYP3A4 Cytochrome has been investigated by density functional QM calculations aided with molecular mechanics/molecular dynamics simulations. Two different orientations of phenyl ring for substrate approach toward oxyferryl center, imposing two subsequent rearrangement pathways have been investigated. Starting from σ-complex in perpendicular orientation enzymatic mechanism involves consecutive proton shuttle intermediate, which further leads to the formation of alcohol and ketone. Parallel conformation leads solely to ketone product by 1,2 hydride shift. Although parallel and perpendicular σ-complexes are energetically equivalent both for the gas phase or PCM solvent model, molecular dynamics studies in full CYP3A4 environment show that perpendicular conformation of the σ-complex should be privileged, stabilized by hydrophobic interactions of phenylacrylamide chain. After assessing probability of the two conformations we postulate that the alcohol, accessible with the lowest energy barriers should be the major metabolite for studied substrate and CYP3A4 enzyme. Figure Orientation of phenyl ring towards porphyrin plane selected by substrate interaction with enzymatic cavity channels enzymatic reaction     

Rotational model for nematic phase stability of fluorinated phenylbicyclohexane liquid crystals

A mechanical molecular rotation model for liquid crystal (LC) systems is employed to evaluate phase transition temperature of fluorinated phenylbicyclohexane isomeric LC compounds. Results show that when a fluorine atom is substituted along the molecular long axis, an LC molecule acquires high rotational speed and its rotation becomes stable, thereby resulting in a better thermal stability of the nematic phase. A novel explanation is proposed for the behavior of the nematic-isotropic phase of the LC system when a heavy atom is substituted along the molecular long axis. Figure Molecular conformation of fluorinated bicyclohexylphenyl compounds. . The fluorine atoms are substituted in different positions 2, 3, 4, and 5 of the phenyl ring, respectively. The axis expresses molecular long rotation axis.     

A hydrophobic similarity analysis of solvation effects on nucleic acid bases

We investigate the changes in the solvation properties of the natural nucleic acid bases due to the formation of the canonical Watson–Crick hydrogen-bonded complexes. To this end, the changes in the free energy of solvation of the bases induced upon hydrogen-bonded dimerization are analyzed by means of the hydrophobic similarity index, which relies on the atomic contributions to the free energy of solvation determined by the partitioning method implemented in the framework of the MST continuum model. Such an index is also used to examine the hydrophobic similarity between the canonical nucleic acid bases and a series of highly apolar analogues, which have been designed as potential candidates to expand the genetic alphabet. The ability of these analogues to be incorporated into modified DNA duplexes can be related to the large reduction in the hydrophilicity of the natural bases upon formation of the canonical hydrogen-bonded dimers. The results illustrate the suitability of the hydrophobic similarity index to rationalize the role played by solvation in molecular recognition. Proceedings of “Modeling Interactions in Biomolecules II”, Prague, September 5th–9th, 2005.     

Possible dynamic anchor points in a benzoxazinone derivative–human oxytocin receptor system — a molecular docking and dynamics calculation

In this study, we performed a molecular docking and dynamics simulation for a benzoxazinone–human oxytocin receptor system to determine the possible hydrophobic and electrostatic interaction points in the dynamic complex. After the homology modeling, the ligand was docked into the putative active using AutoDock 3.05. After the application of energetic and structural filters, the complexes obtained were further refined with a simulated annealing protocol (AMBER8) to remove steric clashes. Three complexes were selected for subjection to the molecular dynamics simulation (5 ns), and the results on the occurrence of average anchor points showed a stable complex between the benzoxazinone derivative and the receptor. The complex could be used as a good starting point for further analysis with site-directed mutagenesis, or further computational research. Figure The location of the ligands (complex B – blue; complex E – red; and complex F – green) in the transmembrane regions (TM1 – red; TM2 – blue; TM3 – yellow; TM4 – purple; TM5 – orange; TM6 – cyan; TM7 – pink) of the hOTR. For clarity, the EC and IC loops are not shown Electronic Supplementary Material Supplementary material is available at     

Optimization of cutting schemes for the evaluation of molecular electrostatic potentials in proteins via Moving-Domain QM/MM

This work presents new developments of the moving-domain QM/MM (MoD-QM/MM) method for modeling protein electrostatic potentials. The underlying goal of the method is to map the electronic density of a specific protein configuration into a point-charge distribution. Important modifications of the general strategy of the MoD-QM/MM method involve new partitioning and fitting schemes and the incorporation of dynamic effects via a single-step free energy perturbation approach (FEP). Selection of moderately sized QM domains partitioned between and C (from C=O), with incorporation of delocalization of electrons over neighboring domains, results in a marked improvement of the calculated molecular electrostatic potential (MEP). More importantly, we show that the evaluation of the electrostatic potential can be carried out on a dynamic framework by evaluating the free energy difference between a non-polarized MEP and a polarized MEP. A simplified form of the potassium ion channel protein Gramicidin-A from Bacillus brevis is used as the model system for the calculation of MEP. Figure Schematic representation of the Moving Domain QM/MM method