Conformational properties of penicillins: quantum chemical calculations and molecular dynamics simulations of benzylpenicillin

Herein, we present theoretical results on the conformational properties of benzylpenicillin, which are characterized by means of quantum chemical calculations (MP2/6-31G* and B3LYP/6-31G*) and classical molecular dynamics simulations (5 ns) both in the gas phase and in aqueous solution. In the gas phase, the benzylpenicillin conformer in which the thiazolidine ring has the carboxylate group oriented axially is the most favored one. Both intramolecular CH. O and dispersion interactions contribute to stabilize the axial conformer with respect to the equatorial one. In aqueous solution, a molecular dynamics simulation predicts a relative population of the axial:equatorial conformers of 0.70:0.30 in consonance with NMR experimental data. Overall, the quantum chemical calculations as well as the simulations give insight into substituent effects, the conformational dynamics of benzylpenicillin, the frequency of ring-puckering motions, and the correlation of side chain and ring-puckering motions.     

Orientation and conformational preference of leucine-enkephalin at the surface of a hydrated dimyristoylphosphatidylcholine bilayer: NMR and MD simulation

The morphogenic opiate pentapeptide leucine-enkephalin (lenk) in a hydrated dimyristoylphosphatidylcholine (DMPC) bilayer is studied using NMR spectroscopy and molecular dynamics simulation. Contrary to the frequent assumption that the peptide attains a single fixed conformation in the presence of membranes, we find that the lenk molecule is flexible, switching between specific bent conformations. The constraints to the orientation of the aromatic rings that are identified by the NMR experiment are found by the MD simulation to be related to the depth of the peptide in the bilayer. The motion of the N-H vectors of the peptide bonds with respect to the magnetic field direction as observed by MD largely explain the magnitude of the observed residual dipolar coupling (RDC), which are much reduced over the static (15)N-(1)H coupling. The measured RDCs are nevertheless significantly larger than the predicted ones, possibly due the absence of long-time motions in the simulations. The conformational behavior of lenk at the DMPC surface is compared to that in the aqueous solution, both in the neutral and in the zwitterionic forms.     

Comparison of aqueous molecular dynamics with NMR relaxation and residual dipolar couplings favors internal motion in a mannose oligosaccharide.

An investigation has been performed to assess how aqueous dynamical simulations of flexible molecules can be compared against NMR data. The methodology compares state-of-the-art NMR data (residual dipolar coupling, NOESY, and (13)C relaxation) to molecular dynamics simulations in water over several nanoseconds. In contrast to many previous applications of residual dipolar coupling in structure investigations of biomolecules, the approach described here uses molecular dynamics simulations to provide a dynamic representation of the molecule. A mannose pentasaccharide, alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->3)-alpha-D-Manp-(1-->2)-D-Manp, was chosen as the model compound for this study. The presence of alpha-linked mannan is common to many glycopeptides, and therefore an understanding of the structure and the dynamics of this molecule is of both chemical and biological importance. This paper sets out to address the following questions. (1) Are the single structures which have been used to interpret residual dipolar couplings a useful representation of this molecule? (2) If dynamic flexibility is included in a representation of the molecule, can relaxation and residual dipolar coupling data then be simultaneously satisfied? (3) Do aqueous molecular dynamics simulations provide a reasonable representation of the dynamics present in the molecule and its interaction with water? In summary, two aqueous molecular dynamics simulations, each of 20 ns, were computed. They were started from two distant conformations and both converged to one flexible ensemble. The measured residual dipolar couplings were in agreement with predictions made by averaging the whole ensemble and from a specific single structure selected from the ensemble. However, the inclusion of internal motion was necessary to rationalize the relaxation data. Therefore, it is proposed that although residual dipolar couplings can be interpreted as a single-structure, this may not be a correct interpretation of molecular conformation in light of other experimental data. Second, the methodology described here shows that the ensembles from aqueous molecular dynamics can be effectively tested against experimental data sets. In the simulation, significant conformational motion was observed at each of the linkages, and no evidence for intramolecular hydrogen bonds at either alpha(1-->2) or alpha(1-->3) linkages was found. This is in contrast to simulations of other linkages, such as beta(1-->4), which are often predicted to maintain intramolecular hydrogen bonds and are coincidentally predicted to have less conformational freedom in solution.     

Kinetic barriers in the isomerization of substituted ureas: implications for computer-aided drug design

Urea derivatives are ubiquitously found in many chemical disciplines. N,N′-substituted ureas may show different conformational preferences depending on their substitution pattern. The high energetic barrier for isomerization of the cis and trans state poses additional challenges on computational simulation techniques aiming at a reproduction of the biological properties of urea derivatives. Herein, we investigate energetics of urea conformations and their interconversion using a broad spectrum of methodologies ranging from data mining, via quantum chemistry to molecular dynamics simulation and free energy calculations. We find that the inversion of urea conformations is inherently slow and beyond the time scale of typical simulation protocols. Therefore, extra care needs to be taken by computational chemists to work with appropriate model systems. We find that both knowledge-driven approaches as well as physics-based methods may guide molecular modelers towards accurate starting structures for expensive calculations to ensure that conformations of urea derivatives are modeled as adequately as possible.     

A method for conformational sampling of loops in proteins based on adiabatic decoupling and temperature or force scaling

A method for conformational Boltzmann sampling of loops in proteins in aqueous solution is presented that is based on adiabatic decoupling molecular dynamics (MD) simulation with temperature or force scaling. To illustrate the enhanced sampling, the loop from residues 33 to 43 in the bovine protein ribonuclease A is adiabatically decoupled from the rest of the protein and the solvent with a mass scaling factor s(m) =1000 and the sampling is enhanced with a scaling of the temperature using s(T) =2 or of the force using s(V) =0.667. Over 5 ns of simulation the secondary structure of the protein remains unaltered while a combined dihedral-angle conformational cluster analysis shows an increase of conformations outside the first most populated cluster of loop conformations for adiabatic decoupling MD with temperature scaling using s(T) =2 or force scaling using s(V) =0.667 compared to the standard MD simulation. The atom-positional root-mean-square fluctuations of the C(α) atoms of the loop show an increase in the movement of the loop as well, indicating that adiabatic decoupling MD with upscaling of the temperature or downscaling of the force is a promising method for conformational Boltzmann sampling.     

Conformational behaviour of the antineoplastic peptide dolastatin-10 and of two mutated derivatives

Summary The three-dimensional structure of dolastatin-10, an extremely potent cytostatic and antineoplastic peptide extracted from the mollusc Dolabella auricularia, has not yet been fully characterized in an experimental way. By means of a systematic conformational search of the natural peptide and of two mutated analogs, carried out both in vacuo and in aqueous solution, the present work allows to obtain insights into the conformational preferences of this remarkable compound. In addition, the ability to form intra- and intermolecular H-bonds as a function both of the sequence and of the conformation is discussed. The search for the best molecular conformations has been carried out using a molecular mechanics approach, based on the CVFF potential. Dolastatin-10 contains some unusual amino acids for which no experimental structural data are available. In order to check the reliability of the CVFF potential in predicting structures of such nonconventional amino acids, geometry optimizations have been carried out using the ab initio Hartree-Fock procedure. The CVFF parameterization is found to be adequate also for nonconventional amino acids.     

Molecular mechanics conformational analysis of tylosin

The conformations of the 16-membered macrolide antibiotic tylosin were studied with molecular mechanics (AMBER* force field) including modelling of the effect of the solvent on the conformational preferences (GB/SA). A Monte Carlo conformational search procedure was used for finding the most probable low-energy conformations. The present study provides complementary data to recently reported analysis of the conformations of tylosin based on NMR techniques. A search for the low-energy conformations of protynolide, a 16-membered lactone containing the same aglycone as tylosin, was also carried out, and the results were compared with the observed conformation in the crystal as well as with the most probable conformations of the macrocyclic ring of tylosin. The dependence of the results on force field was also studied by utilizing the MM3 force field. Some particular conformations were computed with the semiempirical molecular orbital methods AM1 and PM3.     

Application of self-organizing maps in conformational analysis of lipids.

The characteristics of lipid assemblies are important for the functions of biological membranes. This has led to an increasing utilization of molecular dynamics simulations for the elucidation of the structural features of biomembranes. We have applied the self-organizing map (SOM) to the analysis of the complex conformational data from a 1-ns molecular dynamics simulation of PLPC phospholipids in a membrane assembly. Mapping of 1.44 million molecular conformations to a two-dimensional array of neurons revealed, without human intervention, the main conformational features in hours. Both the whole molecule and the characteristics of the unsaturated fatty acid chains were analyzed. All major structural features were easily distinguished, such as the orientational variability of the headgroup, the mainly trans state dihedral angles of the sn-1 chain, and both straight and bent conformations of the unsaturated sn-2 chain. Furthermore, presentation of the trajectory of an individual lipid molecule on the map provides information on conformational dynamics. The present results suggest that the SOM method provides a powerful tool for routinely gaining rapid insight to the main molecular conformations as well as to the conformational dynamics of any simulated molecular assembly without the requirement of a priori knowledge.     

A comparison of distance geometry and molecular dynamics simulation techniques for conformational analysis of β-cyclodextrin

Distance geometry and molecular dynamics simulation techniques were compared in their ability to search the conformational potential energy surface of β-cyclodextrin. Structures generated by the DISGEO program were minimized using three different atomic point charge sets. Some of these structures were used as starting points for molecular dynamics simulation in vacuo at 298K. The distance geometry results showed that the global features of the conformational potential energy surface were generally independent of the point charge set. The distance geometry technique was able to find structures of lower energy than those obtained by direct minimization of the X-ray or neutron diffraction structures. However, the molecular dynamics simulation technique was consistently able to find structures of lower energy than those generated by distance geometry. Root mean square fit of the trajectory structures to the starting structure showed that the simulation allowed the molecule to explore regions of the potential energy surface in the neighborhood of the starting structure. Both the distance geometry and molecular dynamics simulation techniques showed that β-cyclodextrin can adopt a wide range of conformations in the gas phase and that these conformations are much less symmetrical than the crystalline structure.     

Characterization of structural transitions from aqueous solution to a lipid phase for α-MSH

Fluorescence spectroscopy results show that the α-melanocyte-stimulating hormone peptide (α-MSH) interacts with acidic lipid vesicles. Detectable structural changes are concomitant with the passage of a tryptophan residue from aqueous to lipidic media. The observed multiexponential decay of fluorescence, rationalized as originating from three rotameric populations of the tryptophan residue, has been used together with a matrix algorithm to find the most probable conformational families of α-MSH in water and lipid environments. A model is discussed in which the same conformational families occur in various phases, although with different probabilities. A conformational family in which χ₁ of the Trp9 side chain is in the trans-rotameric conformation is shown to have structural features highly appropriate to interact with negatively charged biological membranes, which are also in accordance with previous molecular dynamics simulations and with structures engineered in α-MSH analogs that show an increased potency in biological essays. The gauche minus and gauche plus side-chain conformations of Trp9, on the other hand, yield conformations more likely to predominate in aqueous solution. NMR spectroscopy measurements of α-MSH analogs indicate the existence in aqueous solution of a β strand in the vicinity of Trp9. A similar structural feature was found in the present conformational analysis for the gauche minus and gauche plus side-chain rotamers of Trp9. © 1995 John Wiley & Sons, Inc.     

RNA unrestrained molecular dynamics ensemble improves agreement with experimental NMR data compared to single static structure: a test case

Nuclear magnetic resonance (NMR) provides structural and dynamic information reflecting an average, often non-linear, of multiple solution-state conformations. Therefore, a single optimized structure derived from NMR refinement may be misleading if the NMR data actually result from averaging of distinct conformers. It is hypothesized that a conformational ensemble generated by a valid molecular dynamics (MD) simulation should be able to improve agreement with the NMR data set compared with the single optimized starting structure. Using a model system consisting of two sequence-related self-complementary ribonucleotide octamers for which NMR data was available, 0.3 ns particle mesh Ewald MD simulations were performed in the AMBER force field in the presence of explicit water and counterions. Agreement of the averaged properties of the molecular dynamics ensembles with NMR data such as homonuclear proton nuclear Overhauser effect (NOE)-based distance constraints, homonuclear proton and heteronuclear ¹H–³¹P coupling constant (J) data, and qualitative NMR information on hydrogen bond occupancy, was systematically assessed. Despite the short length of the simulation, the ensemble generated from it agreed with the NMR experimental constraints more completely than the single optimized NMR structure. This suggests that short unrestrained MD simulations may be of utility in interpreting NMR results. As expected, a 0.5 ns simulation utilizing a distance dependent dielectric did not improve agreement with the NMR data, consistent with its inferior exploration of conformational space as assessed by 2-D RMSD plots. Thus, ability to rapidly improve agreement with NMR constraints may be a sensitive diagnostic of the MD methods themselves.     

单嘧磺隆晶体-活性构象转换的分子动力学模拟

应用分子动力学模拟方法对单嘧磺隆在水、正辛醇和正辛烷3种不同溶剂中的构象行为、单嘧磺隆与3种溶剂之间的相互作用能及氢键相互作用进行了计算研究. 计算结果表明, 在3种不同的溶剂中, 单嘧磺隆的优势构象不同; 其构象转换过程, 特别是转换成活性构象的过程主要发生在水溶液中; 与溶剂分子间的相互作用是分子构象行为的决定因素; 单嘧磺隆的脲桥部分可以和含氢键接受体的溶剂形成氢键, 分子间与分子内氢键的竞争可能是从晶体构象转换成活性构象的主要驱动力.     

Computing the (1)H NMR spectrum of a bulk ionic liquid from snapshots of car-parrinello molecular dynamics simulations.

We have investigated the performance of several computational protocols in predicting the NMR spectrum of a molecular ion in a complex liquid phase such as an ionic liquid. To do this, we computed the proton NMR chemical shifts of the 1-ethyl-3-methylimidazolium cation [emim](+) in [emim][Cl]. Environmental effects on the imidazolium ring proton chemical shifts are quite significant and must be taken into account explicitly. Calculations performed on the isolated imidazolium cation as well as on the [emim][Cl] ion pair grossly fail to reproduce the correct spacing between proton signals. In contrast, calculations performed on clusters extracted from the trajectory of a Car-Parrinello molecular dynamics simulation yield very good results.     

Molecular dynamics simulation of the A-DNA to B-DNA transition in aqueous RbCl solution

Unrestrained molecular dynamics (MD) simulations have been carried out to characterize the stability of DNA conformations and the dynamics of A-DNA→B-DNA conformational transitions in aqueous RbCl solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbCl concentration on the dynamics of the A→B conformational transition in the DNA duplex d(CGCGAATTCGCG)2 . Canonical Aand B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was monitored for 3.0 ns and the distances between the C5′ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simulations. In addition, increasing the RbCl concentration in aqueous solution hindered the A→B conformational transition and the transition in aqueous RbCl solution was faster than that in aqueous NaCl solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A→B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The rubidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help explain the effect of RbCl concentration on the dynamics of the A→B conformational transition.     

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 search by potential energy annealing: Algorithm and application to cyclosporin A

Summary A major problem in modelling (biological) macromolecules is the search for low-energy conformations. The complexity of a conformational search problem increases exponentially with the number of degrees of freedom which means that a systematic search can only be performed for very small structures. Here we introduce a new method (PEACS) which has a far better performance than conventional search methods.To show the advantages of PEACS we applied it to the refinement of Cyclosporin A and compared the results with normal molecular dynamics (MD) refinement. The structures obtained with PEACS were lower in energy and agreed with the NMR parameters much better than those obtained with MD. From the results it is further clear that PEACS samples a much larger part of the available conformational space than MD does.     

Interpreting protein structural dynamics from NMR chemical shifts

In this investigation, semiempirical NMR chemical shift prediction methods are used to evaluate the dynamically averaged values of backbone chemical shifts obtained from unbiased molecular dynamics (MD) simulations of proteins. MD-averaged chemical shift predictions generally improve agreement with experimental values when compared to predictions made from static X-ray structures. Improved chemical shift predictions result from population-weighted sampling of multiple conformational states and from sampling smaller fluctuations within conformational basins. Improved chemical shift predictions also result from discrete changes to conformations observed in X-ray structures, which may result from crystal contacts, and are not always reflective of conformational dynamics in solution. Chemical shifts are sensitive reporters of fluctuations in backbone and side chain torsional angles, and averaged (1)H chemical shifts are particularly sensitive reporters of fluctuations in aromatic ring positions and geometries of hydrogen bonds. In addition, poor predictions of MD-averaged chemical shifts can identify spurious conformations and motions observed in MD simulations that may result from force field deficiencies or insufficient sampling and can also suggest subsets of conformational space that are more consistent with experimental data. These results suggest that the analysis of dynamically averaged NMR chemical shifts from MD simulations can serve as a powerful approach for characterizing protein motions in atomistic detail.     

An application of coupled reference interaction site model/molecular dynamics to the conformational analysis of the alanine dipeptide

We present an application of our recently proposed coupled reference interaction site model (RISM) molecular dynamics (MD) solvation free energy methodology [Freedman and Truong, Chem. Phys. Lett. 381, 362 (2003); J. Chem. Phys. 121, 2187 (2004)] to study the conformational stability of alanine dipeptide in aqueous solution. In this methodology, radial distribution functions obtained from a single MD simulation are substituted into a RISM expression for solvation free energy. Consequently, iterative solution of the RISM equation is not needed. The relative solvation free energies of seven different conformations of the alanine dipeptide in aqueous solution are calculated. Results from the coupled RISM/MD methodology are in good agreement with those from earlier simulations using the accurate free energy perturbation approach, showing that the alphaR conformation is most stabilized by solution. This study establishes a framework for applying this coupled RISM/MD method to larger biological systems.     

Dependency of ligand free energy landscapes on charge parameters and solvent models

We performed replica-exchange molecular dynamics (REMD) simulations of six ligands to examine the dependency of their free energy landscapes on charge parameters and solvent models. Six different charge parameter sets for each ligand were first generated by RESP and AM1-BCC methods using three different conformations independently. RESP charges showed some conformational dependency. On the other hand, AM1-BCC charges did not show conformational dependency and well reproduced the overall trend of RESP charges. The free energy landscapes obtained from the REMD simulations of ligands in vacuum, Generalized-Born (GB), and TIP3P solutions were then analyzed. We found that even small charge differences can produce qualitatively different landscapes in vacuum condition, but the differences tend to be much smaller under GB and TIP3P conditions. The simulations in the GB model well reproduced the landscapes in the TIP3P model using only a fraction of the computational cost. The protein-bound ligand conformations were rarely the global minimum states, but similar conformations were found to exist in aqueous solution without proteins in regions close to the global minimum, local minimum or intermediate states.