Conformational analysis of six- and twelve-membered ring compounds by molecular dynamics

A molecular dynamics (MD)-based conformational analysis has been performed on a number of cycloalkanes in order to demonstrate the reliability and generality of MD as a tool for conformational analysis. MD simulations on cyclohexane and a series of methyl-substituted cyclohexanes were performed at temperatures between 400 and 1200 K. Depending on the simulation temperature, different types of interconversions (twist-boat–twist-boat, twist- boat–chair and chair–chair) could be observed, and the MD simulations demonstrated the expected correlation between simulation temperature and ring inversion barriers. A series of methyl-substituted 1,3- dioxanes were investigated at 1000 K, and the number of chair–chair interconversions could be quantitatively correlated to the experimentally determined ring inversion barrier. Similarly, the distribution of sampled minimum-energy conformations correlated with the energy-derived Boltzmann distribution. The macrocyclic ring system cyclododecane was subjected to an MD simulation at 1000 K and 71 different conformations could be sampled. These conformations were compared with the results of previously reported conformational analyses using stochastic search methods, and the MD method provided 19 out of the 20 most stable conformations found in the MM2 force field. Finally, the general performance of the MD method for conformational analysis is discussed.     

Conformational analysis and molecular dynamics simulation of cellobiose and larger cellooligomers

The conformational behavior of cellobiose (D -glc-ß(1→4)-D -glc), cellotetraose, and cellooctaose was studied by a combination of energy minimization and molecular dynamics simulations in vacuo at 400 K. These diand oligosaccharide models have considerable flexibility and exhibit a variety of different motions in glycosidic and exocyclic torsions. The glycosidic ?, ψ torsions moved frequently between two local minima on the cellobiose energy surface in the region of known crystal structures. Transitions of the hydroxymethyl side chain were observed between gt,gg, and tg conformations accompanied by changes in intramolecular hydrogen bonding patterns. A reasonable fit to the experimental optical rotation and nuclear magnetic resonance vicinal coupling data of cellobiose in solution required a distribution of its conformations. The oligomers, although generally extended, assumed a more coiled or twisted shape than is observed in the crystalline state of cellulose and exhibited considerable backbone motion due to local ring rotations about the glycosidic bonds. Long-lived transitions to structures having torsion angles 180° from the major minima (ring flips) introduced kinks and bends into the tetramer and octamer. While the glucose rings of the structures remained primarily in the ⁴C₁ conformation, twist and boat structures were also observed in the tetramer and octamer structures. Reducing the simulation temperature to 300 K eliminated some of the transitions seen at 400 K. © 1993 John Wiley & Sons, Inc.     

Conformational analysis of methylphenidate: comparison of molecular orbital and molecular mechanics methods

Summary Methylphenidate (MP) binds to the cocaine binding site on the dopamine transporter and inhibits reuptake of dopamine, but does not appear to have the same abuse potential as cocaine. This study, part of a comprehensive effort to identify a drug treatment for cocaine abuse, investigates the effect of choice of calculation technique and of solvent model on the conformational potential energy surface (PES) of MP and a rigid methylphenidate (RMP) analogue which exhibits the same dopamine transporter binding affinity as MP. Conformational analysis was carried out by the AM1 and AM1/SM5.4 semiempirical molecular orbital methods, a molecular mechanics method (Tripos force field with the dielectric set equal to that of vacuum or water) and the HF/6-31G* molecular orbital method in vacuum phase. Although all three methods differ somewhat in the local details of the PES, the general trends are the same for neutral and protonated MP. In vacuum phase, protonation has a distinctive effect in decreasing the regions of space available to the local conformational minima. Solvent has little effect on the PES of the neutral molecule and tends to stabilize the protonated species. The random search (RS) conformational analysis technique using the Tripos force field was found to be capable of locating the minima found by the molecular orbital methods using systematic grid search. This suggests that the RS/Tripos force field/vacuum phase protocol is a reasonable choice for locating the local minima of MP. However, the Tripos force field gave significantly larger phenyl ring rotational barriers than the molecular orbital methods for MP and RMP. For both the neutral and protonated cases, all three methods found the phenyl ring rotational barriers for the RMP conformers/invertamers (denoted as cte, tte, and cta) to be: cte, tte> MP > cta. Solvation has negligible effect on the phenyl ring rotational barrier of RMP. The B3LYP/6-31G* density functional method was used to calculate the phenyl ring rotational barrier for neutral MP and gave results very similar to those of the HF/6-31G* method.     

Conformational analysis of hydroxymatairesinol in aqueous solution with molecular dynamics

Molecular dynamics simulations were performed on the naturally occuring lignan hydroxymatairesinol (HMR) using the GROMACS software. The aim of this study was to explore the conformational behavior of HMR in aqueous solution adopting the TIP4P model. The topology of HMR was constructed by hand and HMR was modeled with the OPLS‐AA force field implemented in GROMACS. The five torsional angles in HMR were properly analyzed during the simulations. Correlations through certain patterns were observed between the angles. The determining property for the conformation preferred in aqueous solution was found to be the dipole moment and not the lowest energy in gas phase. The solvation effects on HMR was also studied by quantum chemical calculations applying the COnductorlike Screening MOdel (COSMO), the results of which were compared with results from a previous study using the Polarized Continuum Model (PCM). In the present work, COSMO was found to give more credible relative energies than PCM. © 2009 Wiley Periodicals, Inc., J Comput Chem, 2009     

Sialyldisaccharide conformations: a molecular dynamics perspective

Sialyldisaccharides are significant terminal components of glycoconjugates and their negative charge and conformation are extensively utilized in molecular recognition processes. The conformation and flexibility of four biologically important sialyldisaccharides [Neu5Acα(2-3)Gal, Neu5Acα(2-6)Gal, Neu5Acα(2-8)Neu5Ac and Neu5Acα(2-9)Neu5Ac] are studied using Molecular Dynamics simulations of 20 ns duration to deduce the conformational preferences of the sialyldisaccharides and the interactions which stabilize the conformations. This study clearly describes the possible conformational models of sialyldisaccharides deduced from 20 ns Molecular Dynamics simulations and our results confirm the role of water in the structural stabilization of sialyldisaccharides. An extensive analysis on the sialyldisaccharide structures available in PDB also confirms the conformational regions found by experiments are detected in MD simulations of 20 ns duration. The three dimensional structural coordinates for all the MD derived sialyldisaccharide conformations are deposited in the 3DSDSCAR database and these conformational models will be useful for glycobiologists and biotechnologists to understand the biological functions of sialic acid containing glycoconjugates.     

Cluster analysis of molecular conformations

We describe a method for locating clusters of geometrically similar conformers in ensembles of chemical conformations. We first calculate the pairwise interconformational distance matrix in either torsional or Cartesian space and then use an agglomerative, single-link clustering method to define a hierarchy of clusterings in the same space. Especially good clusterings are distinguished by high values of the separation ratio: the ratio of the shortest intercluster distance to the characteristic threshold distance defining the clustering. We also discuss other statistics. The method has been embodied in a program called XCluster, which can display the distance matrix, the hierarchy of clusterings, and the clustering statistics in a variety of formats. XCluster can also write out the clustered conformations for subsequent or simultaneous viewing with a molecular visualization program. We demonstrate the sorts of insight that this approach affords with examples obtained from conformational search and molecular dynamics procedures. © 1994 by John Wiley & Sons, Inc.     

Conformational analysis of colchicine and isocolchicine by molecular mechanics

The structures of isocolchicine ( ) and colchicine ( ) have been calculated using the MMX routine. The low energy conformations for isocolchicine and colchicine fit well with x-ray crystallographic data. The B ring atropisomer of isocolchicine, which can be spectroscopically observed, is calculated to be <1 kcal/mole higher in energy than . The boat-boat inversion conformer of colchicine, which has been predicted to be important in the binding of to tubulin, is also calculated. The B ring geometry of this isomer does not differ to the extent previously predicted.     

Conformational analysis of a tetrasaccharide based on NMR spectroscopy and molecular dynamics simulations

The conformational preference of the human milk oligosaccharide lacto-N-neotetraose, beta-d-Galp-(1 --> 4)-beta-d-GlcpNAc-(1 --> 3)-beta-d-Galp-(1 --> 4)-d-Glcp, has been analyzed using (1)H,(1)H T-ROESY and (1)H,(13)C trans-glycosidic J coupling experiments in isotropic solution and (1)H,(13)C residual dipolar couplings (RDCs) obtained in lyotropic liquid crystalline media. Molecular dynamics simulations of the tetrasaccharide with explicit water as the solvent revealed that two conformational states are significantly populated at the psi glycosidic torsion angle, defined by C(anomeric)-O-C-H, of the (1 --> 3)-linkage. Calculation of order parameters, related to the molecular shape, were based on the inertia tensor and fitting of experimental RDCs to different conformational states showed that psi(+) > 0 degrees is the major and psi(-) < 0 degrees is the minor conformation in solution, in complete agreement with a two-state analysis based on the T-ROESY data. Attention was also given to the effect of salt (200 mM NaCl) in the anisotropic medium, which was a ternary mixture of n-octyl-penta(ethylene glycol), n-octanol, and D(2)O.     

Analysis of the RGD sequence in protein structures: comparison to the conformations of the RGDW and DRGDW peptides determined by molecular dynamics simulations

Geometric properties of the RGD sequence in a data set of protein crystal and NMR structures deposited in the Protein Data Bank were examined to identify structural characteristics that are related to cell adhesion activity. Interatomic distances and dihedral angles are examined. These geometric measures are then used in an analysis of the conformations of the RGDW and _DRGDW peptides obtained from molecular dynamics simulations (Stote RH, et al. (2000) J Phys ChemB 104:1624). This analysis leads to the suggestion that differences in the accessible conformations contribute to the difference in biological activity between the RGDW and the _DRGDW peptides. Received: 15 August 2000 / Accepted: 4 October 2000 / Published online: 21 March 2001     

Conformational analysis of gossypol and its derivatives by molecular mechanics

Conformations and inversion pathways leading to racemization of all the tautomers of gossypol, gossypolone, anhydrogossypol, and a diethylamine Schiff's base of gossypol were investigated with MM3(2000). All forms have hindered rotation because of clashes between the methyl carbon atom and oxygen-containing moieties ortho to the bond linking the two naphthalene rings. Inversion energies generally agree with available experimental data. Gossypol preferentially inverts in its dihemiacetal tautomeric form through the cis pathway (where similar groups clash). Gossypolone inverts more easily than gossypol, and preferentially through the trans pathway (where dissimilar groups clash) when one of its outer rings has an enol-keto group and the other has an aldehyde group. Anhydrogossypol racemizes through the cis pathway. The bridge bond and the ortho exo-cyclic bonds in all the structures bend from planarity, and the inner naphthalene rings pucker to accommodate the inversion. For gossypol, the transition is achieved through greater bending of the exo-cyclic bonds (up to 12°) and less distortion of the inner benzyl rings (q≤0.34 Å), (up to 12.7°) . For gossypolone the transition occurs with greater distortion of the inner benzyl rings (q≤0.63 Å) and less out-of-plane bending (up to 8.4°). By isolating individual clashes, their contribution to the overall barrier can be analyzed, as shown for the dialdehyde tautomer of gossypol.     

Conformational analysis of dehydrodidemnin B (aplidine) by NMR spectroscopy and molecular mechanics/dynamics calculations.

Dehydrodidemnin B (DDB or aplidine), a potent antitumoral natural product currently in phase II clinical trials, exists as an approximately 1:1 mixture of two slowly interconverting conformations. These are sufficiently long-lived so as to allow their resolution by HPLC. NMR spectroscopy shows that this phenomenon is a consequence of restricted rotation about the Pyr-Pro(8) terminal amide bond of the molecule's side chain. The same technique also indicates that the overall three-dimensional structures of both the cis and trans isomers of DDB are similar despite the conformational change. Molecular dynamics simulations with different implicit and explicit solvent models show that the ensembles of three-dimensional structures produced are indeed similar for both the cis and trans isomers. These studies also show that hydrogen bonding patterns in both isomers are alike and that each one is stabilized by a hydrogen bond between the pyruvyl unit at the terminus of the molecule's side chain and the Thr(6) residue situated at the junction betwen the macrocycle and the molecule's side chain. Nevertheless, each conformational isomer forms this hydrogen bond using a different pyruvyl carbonyl group: CO(2) in the case of the cis isomer and CO(1) in the case of the trans isomer.     

Conformational and dynamic effects of pendant-PEK chains on rigid-rod polymers by molecular dynamics simulation

Molecular dynamics simulations were performed with models of a single-component ‘hairy-rod’ molecular composite composed of flexible meta poly(aryl ether ketone) (mPEK) chains (the ‘hairs’) grafted to a poly(p-phenylene benzobisthiazole) (PBZT) rigid-rod backbone. The molecular-composite concept, conceived to circumvent solubility problems and improve compressive strength of PBZT, relies on an even distribution of rods in a coil-like matrix. Two-molecule simulations show that the pendant mPEK chains associate with the rods causing large bends in the rod backbone. Simulations of bulk systems imply that the mPEK chains greatly reduce the correlation of rod orientation without altering rod spacing. Radial distribution and orientation correlation functions as well as correlation volumes are computed.     

Conformational analysis of oxazolidines by molecular and quantum mechanics and NMR spectroscopy

A combined experimental (¹HNMR) and theoretical (molecular mechanics and PCILO) study of (4S,5S,6R)-4-carbome-thoxyethylenyl-N-carbobenzyloxy-5-methyl-6-phenyl oxazolidine has been carried out with the purpose of contributing to a better understanding of the steric and electronic factors responsible for the observed diastereoisomeric selectivity in the dihydroxylation of oxazolidines leading to optically active α, β dialkoxyaldehydes.     

Comparison between self-guided Langevin dynamics and molecular dynamics simulations for structure refinement of protein loop conformations

This article presents a comparative analysis of two replica‐exchange simulation methods for the structure refinement of protein loop conformations, starting from low‐resolution predictions. The methods are self‐guided Langevin dynamics (SGLD) and molecular dynamics (MD) with a Nosé–Hoover thermostat. We investigated a small dataset of 8‐ and 12‐residue loops, with the shorter loops placed initially from a coarse‐grained lattice model and the longer loops from an enumeration assembly method (the Loopy program). The CHARMM22 + CMAP force field with a generalized Born implicit solvent model (molecular‐surface parameterized GBSW2) was used to explore conformational space. We also assessed two empirical scoring methods to detect nativelike conformations from decoys: the all‐atom distance‐scaled ideal‐gas reference state (DFIRE‐AA) statistical potential and the Rosetta energy function. Among the eight‐residue loop targets, SGLD out performed MD in all cases, with a median of 0.48 Å reduction in global root‐mean‐square deviation (RMSD) of the loop backbone coordinates from the native structure. Among the more challenging 12‐residue loop targets, SGLD improved the prediction accuracy over MD by a median of 1.31 Å, representing a substantial improvement. The overall median RMSD for SGLD simulations of 12‐residue loops was 0.91 Å, yielding refinement of a median 2.70 Å from initial loop placement. Results from DFIRE‐AA and the Rosetta model applied to rescoring conformations failed to improve the overall detection calculated from the CHARMM force field. We illustrate the advantage of SGLD over the MD simulation model by presenting potential‐energy landscapes for several loop predictions. Our results demonstrate that SGLD significantly outperforms traditional MD in the generation and populating of nativelike loop conformations and that the CHARMM force field performs comparably to other empirical force fields in identifying these conformations from the resulting ensembles. Published 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

Conformational analysis of ochratoxin a by NMR spectroscopy and computational molecular modeling

Two-dimensional NMR spectroscopy has been used for a complete assignment of the proton and carbon-13 spectra of the metabolite from Aspergillus ochraceus, ochratoxin A. In addition, phase-sensitive nuclear Overhauser effect spectrometry experiments and computational molecular modeling (MM2 and MMFF force field programs) have been employed to examine the conformational properties of ochratoxin A in chloroform solutions. Particular attention has been given to intramolecular hydrogen-bonding formation involving the phenolic group on dihydroisocoumarin, which may be responsible for the toxic mechanism of ochratoxin A.     

Calculation of alcohol conformations by molecular mechanics

The geometries and energies of simple alcohols were calculated with a molecular mechanics force field. The force field requires the application of the charge interaction model with charges calculated by the CNDO/2 method, the importance of electrostatic interactions for the equilibrium of rotamers about the C-O bond exceeds that of van der Waals interactions. The calculated rotamer populations are discussed with regard to the value of ¹H NMR coupling constants ³JHCOH and other experimental data.     

Application of torsion angle molecular dynamics for efficient sampling of protein conformations

We investigate the application of torsion angle molecular dynamics (TAMD) to augment conformational sampling of peptides and proteins. Interesting conformational changes in proteins mainly involve torsional degrees of freedom. Carrying out molecular dynamics in torsion space does not only explicitly sample the most relevant degrees of freedom, but also allows larger integration time steps with elimination of the bond and angle degrees of freedom. However, the covalent geometry needs to be fixed during internal coordinate dynamics, which can introduce severe distortions to the underlying potential surface in the extensively parameterized modern Cartesian-based protein force fields. A "projection" approach (Katritch et al. J Comput Chem 2003, 24, 254-265) is extended to construct an accurate internal coordinate force field (ICFF) from a source Cartesian force field. Torsion crossterm corrections constructed from local molecular fragments, together with softened van der Waals and electrostatic interactions, are used to recover the potential surface and incorporate implicit bond and angle flexibility. MD simulations of dipeptide models demonstrate that full flexibility in both the backbone phi/psi and side chain chi1 angles are virtually restored. The efficacy of TAMD in enhancing conformational sampling is then further examined by folding simulations of small peptides and refinement experiments of protein NMR structures. The results show that an increase of several fold in conformational sampling efficiency can be reliably achieved. The current study also reveals some complicated intrinsic properties of internal coordinate dynamics, beyond energy conservation, that can limit the maximum size of the integration time step and thus the achievable gain in sampling efficiency.     

Polarization effects on peptide conformations at water–membrane interface by molecular dynamics simulation

The electrostatic image method was applied to investigate the conformation of peptides characterized by different hydrophobicities in a water–membrane interface model. The interface was represented by a surface of discontinuity between two media with different dielectric constants, taking into account the difference between the polarizabilities of the aqueous medium and the hydrocarbon one. The method consists of a substitution of the real problem, which involves the charges and the induced polarization at the surface of discontinuity, by a simpler problem formed with charges and their images. The electric field due to the polarization induced at the surface by charge q was calculated using a hypothetical charge q′ (image of q), symmetrically located on the opposite side of the surface. The value of q′ was determined using the appropriate electrostatic boundary conditions at the surface. By means of this procedure, the effect of the interface can be introduced easily in the usual force field. We included this extension in the computational package that we are developing for molecular dynamics simulations (Thor ). The peptides studied included hydrophilic tetraaspartic acid (Asp–Asp–Asp–Asp), tetralysine (Lys–Lys–Lys–Lys), hydrophobic tetrapeptide (His–Phe–Arg–Trp), an amphiphilic fragment of β-endorphin, and the signal sequence of the E. coli λ-receptor. The simulation results are in agreement with known experimental data regarding the behavior of peptides at the water–membrane interface. An analysis of the conformational dynamics of the signal sequence peptide at the interface was performed over the course of a few nanoseconds. ©1999 John Wiley & Sons, Inc. J Comput Chem 20: 971–982, 1999