Search PubMed⌕ Search

SEARCH · Search PubMed

Results for “MD simulation”

Search indexed PubMed citations on genomics, clinical trials, systematic reviews and public health. Explore titles, authors and supplied subject terms, then open the PubMed record.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6Linked to original sources
PubMed2002

Bridge function and other structural properties of core-softened model fluids from molecular dynamics simulations.

Molecular dynamics (MD) simulations in three dimensions have been performed for a fluid with particles interacting with a continuous version of the Stell-Hemmer core-softened potential that in two dimensions has been known to reproduce most of the static and dynamic anomalies of liquid water. The pair distribution function obtained from the MD simulation is extrapolated with the help of integral equation theory with a suitable closure relation and the bridge function is extracted. A strong dependence of the bridge function on the system size, i.e., the total number of particles (N) used in the simulation box is observed, which leads to spurious values of the structure factor at long wavelengths. A simple self-consistent correction scheme for the finite size effect has been adopted to correct the bridge function and this scheme produces the correct bridge function even for a small system size. The effects of temperature, number density, and potential parameters on the pair distribution functions and extracted bridge functions are extensively studied.

Journal Article↗
PubMed1993

A model for the dynemicin-A cleavage of DNA using molecular dynamics simulation.

Molecular dynamics (MD) simulations to model possible reaction pathways of the dynemicin-A-DNA cleavage mechanism were performed. Two base-pairs sequences, ApCpTpApCpTp-3'/TpGpApTpGpAp-5' and CpApCpGpGpGp-3'/GpTpGpCpCpCp-5', were considered in the calculations. A model based on a prior study of intercalation of dynemicin-A and posterior activation of the drug was assumed in this study. The minimum energy minor groove intercalation complexes for dynemicin-A were used as starting structures in the MD simulations for the reactive intermediate species involved in the postulated action mechanism. The dynemicin-A diol derivative product of the opening of the epoxy ring was used as a "steric mimic" ligand for the DNA-reactive diaryl intermediate. The calculated changes in the geometry of the intercalation complex, due to the opening of the epoxy ring, correspond to the approach of the postulated intermolecular reaction centers in the intercalation states that are responsible for the highest observed DNA cleavage frequency observed. Conversely, unfavorable reaction geometries were found for the intercalation modes corresponding to the lowest observed DNA cutting frequencies.

Anthraquinones↗
PubMed2006

LiTFSI structure and transport in ethylene carbonate from molecular dynamics simulations.

Molecular dynamics (MD) simulations using a many-body polarizable force field were performed on ethylene carbonate (EC) doped with lithium bistrifluoromethanesulfonamide (LiTFSI) salt as a function of temperature and salt concentration. At 313 K Li+ was coordinated by 2.7-3.2 EC carbonyl oxygen atoms and 0.67-1.05 TFSI- oxygen atoms at EC:Li = 10 and 20 salt concentrations. In completely dissociated electrolytes, however, Li+ was solvated by approximately 3.8 carbonyl oxygen atoms from EC on average. The probability of ions to participate ion aggregates decreased exponentially with an increase in the size of the aggregate. Ion and solvent self-diffusion coefficients and conductivity predicted by MD simulations were in good agreement with experiments. Approximately half of the charge was transported by charged ion aggregates with the other half carried by free (uncomplexed by counterion) ions. Investigation of the Li+ transport mechanism revealed that contribution from the Li+ diffusion together with its coordination shell to the total Li+ transport is similar to the contribution arising from Li+ exchanging solvent molecules in its first coordination shell with solvents from the outer shells.

Journal Article↗
PubMed1991

Mass-weighted molecular dynamics simulation of cyclic polypeptides.

A modified molecular dynamics (MD) method in which atomic masses are weighted was developed previously for studying the conformational flexibility of neuroregulating tetrapeptide Phe-Met-Arg-Phe-amide (FMRF-amide). The method has now been applied to longer and constrained molecules, namely a disulfide-linked cyclic hexapeptide, c[CYFQNC], and its linear and "pseudo-cyclic" analogues. The sampling of dehedral conformational space of teh linear hexapeptide in mass-weighted MD simulations was found to be improved significantly over conventional MD simulations, as in the case of the shorter FMRF-amide molecule studied previously. In the cyclic hexapeptide, the internal constraint of the molecule due to the intramolecular disulfide bond (hence the absence of free terminals in the molecule) does not adversely affect the significant improvement of conformational sampling in mass-weighted MD simulations over normal MD simulations. The pseudo-cyclic polypeptide is identical to the linear CYFQNC molecule in amino acid sequence (i.e., side chains of the cysteine residues are reduced), but the positions of its two terminal heavy atoms were held fixed in space such that the molecule has a nearly cyclic conformation. For this molecule, the mass-weighted MD simulation generated a wide range of polypeptide backbone conformations covering the internal dihedral degrees of freedom; moreover, the physical space of the pseudo-cyclic structure was also sampled in a complete revolution of the entire molecular fragment about the two fixed termini during the simulation. These characteristics suggest that mass-weighted MD can also be an extremely useful method for conformational analyses of constrained molecules and, in particular, for modeling loops on protein surfaces.

Amino Acid Sequence↗
PubMed2000

Cholesterol effects on the phosphatidylcholine bilayer polar region: a molecular simulation study.

A molecular dynamics (MD) simulation of a fully hydrated, liquid-crystalline dimyristoylphosphatidylcholine (DMPC)-Chol bilayer membrane containing approximately 22 mol% Chol was carried out for 4.3 ns. The bilayer reached thermal equilibrium after 2.3 ns of MD simulation. A 2.0-ns trajectory generated during 2.3-4.3 ns of MD simulation was used for analyses to determine the effects of Chol on the membrane/water interfacial region. In this region, 70% of Chol molecules are linked to DMPC molecules via short-distance interactions, where the Chol hydroxyl group (OH-Chol) is 1) charge paired to methyl groups of the DMPC choline moiety ( approximately 34%), via the hydroxyl oxygen atom (Och); 2) water bridged to carbonyl ( approximately 19%) and nonester phosphate ( approximately 14%) oxygen atoms, via both Och and the hydroxyl hydrogen atom (Hch); and 3) directly hydrogen (H) bonded to carbonyl ( approximately 11%) and nonester phosphate ( approximately 5%) oxygen atoms, via Hch ( approximately 17% of DMPC-Chol links are multiple). DMPC's gamma-chain carbonyl oxygen atom is involved in 44% of water bridges and 51% of direct H bonds formed between DMPC and Chol. On average, a Chol molecule forms 0.9 links with DMPC molecules, while a DMPC molecule forms 2.2 and 0.3 links with DMPC and Chol molecules, respectively. OH-Chol makes hydrogen bonds with 1.1 water molecules, preferentially via Hch. The average number of water molecules H bonded to the DMPC headgroup is increased by 7% in the presence of Chol. These results indicate that inclusion of Chol decreases interlipid links and increases hydration in the polar region of the membrane.

Cholesterol↗
PubMed2005

The behaviour of beta-carotene in the phosphatidylcholine bilayer as revealed by a molecular simulation study.

A molecular dynamics (MD) simulation of the fully hydrated bilayer made of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) and containing beta-carotene (beta-Car) molecules was carried out as a complementary approach to experimental techniques to investigate the orientation of beta-Car in the lipid membrane as well as its influence on the bilayer properties. The bilayer reached thermal equilibrium after 1200 ps of MD simulation and the productive run was carried out for 2800 ps. The results indicate that the carotene rings are located in the region occupied by the carbonyl groups of the POPC gamma-chain with no trace of penetration towards the centre of the bilayer. Carotene exhibits an ordering effect on both the beta- and the gamma-chain. While the fully saturated gamma-chain is affected evenly along, the order of the mono-unsaturated beta-chain is modified mainly below the double bond. In general, a high value of the order parameter and the chain tilt in the range from 11.4 degrees to 26.7 degrees were observed for the beta-Car molecules. However, for chain segment adjacent to methyl groups the value of the order parameter is low and the tilt angle is close to 75 degrees . Moreover, the probability of trans conformation being generally close to 1.0 along the beta-Car chain is reduced for these segments. Our MD simulation study suggests two pools of the preferential orientation of beta-Car: a slightly bent structure corresponding to a small chain tilt angle and a rather stretched structure that corresponds to a higher chain tilt. The results are discussed in the light of experimental findings.

Computer Simulation↗
PubMed2002

Molecular dynamics applied to nucleic acids.

In this Account, we focus on molecular dynamics (MD) simulations involving fully solvated nucleic acids. Historically, MD simulations were first applied to proteins and several years later to nucleic acids. The first MD simulations of DNA were carried out in vacuo, but nowadays fully solvated systems are common practice. Recently, technical improvements have made it possible to conduct accurate MD simulations of highly charged nucleic acids. The state-of-the-art of MD simulations and a number of applications on various nucleic acid systems are discussed.

Computer Simulation↗
PubMed2006

Dynamical motions of lipids and a finite size effect in simulations of bilayers.

Molecular dynamics (MD) simulations of dipalmitoylphosphatidylcholine bilayers composed of 72 and 288 lipids are used to examine system size dependence on dynamical properties associated with the particle mesh Ewald (PME) treatment of electrostatic interactions. The lateral diffusion constant Dl is 2.92 x 10(-7) and 0.95 x 10(-7) cm2/s for 72 and 288 lipids, respectively. This dramatic finite size effect originates from the correlation length of lipid diffusion, which extends to next-nearest neighbors in the 288 lipid system. Consequently, diffusional events in smaller systems can propagate across the boundaries of the periodic box. The internal dynamics of lipids calculated from the PME simulations are independent of the system size. Specifically, reorientational correlation functions for the slowly relaxing phosphorus-glycerol hydrogen, phosphorus-nitrogen vectors, and more rapidly relaxing CH vectors in the aliphatic chains are equivalent for the 72 and 288 lipid simulations. A third MD simulation of a bilayer with 72 lipids using spherical force-shift electrostatic cutoffs resulted in interdigitated chains, thereby rendering this cutoff method inappropriate.

1,2-Dipalmitoylphosphatidylcholine↗
PubMed1994

Molecular dynamics simulations of phospholipid bilayers.

Molecular dynamics (MD) simulations at 37 degrees C have been performed on three phospholipid bilayer systems composed of the lipids DLPE, DOPE, and DOPC. The model used included 24 explicit lipid molecules and explicit waters of solvation in the polar head group regions, together with constant-pressure periodic boundary conditions in three dimensions. Using this model, a MD simulation samples part of an infinite planar lipid bilayer. The lipid dynamics and packing behavior were characterized. Furthermore, using the results of the simulations, a number of diverse properties including bilayer structural parameters, hydrocarbon chain order parameters, dihedral conformations, electron density profile, hydration per lipid, and water distribution along the bilayer normal were calculated. Many of these properties are available for the three lipid systems chosen, making them well suited for evaluating the model and protocols used in these simulations by direct comparisons with experimental data. The calculated MD behavior, chain disorder, and lipid packing parameter, i.e. the ratio of the effective areas of hydrocarbon tails and head group per lipid (a(t)/ah), correctly predict the aggregation preferences of the three lipids observed experimentally at 37 degrees C, namely: a gel bilayer for DLPE, a hexagonal tube for DOPE, and a liquid crystalline bilayer for DOPC. In addition, the model and conditions used in the MD simulations led to good agreement of the calculated properties of the bilayers with available experimental results, demonstrating the reliability of the simulations. The effects of the cis unsaturation in the hydrocarbon chains of DOPE and DOPC, compared to the fully saturated one in DLPE, as well as the effects of the different polar head groups of PC and PE with the same unsaturated chains on the lipid packing and bilayer structure have been investigated. The results of these studies indicate the ability of MD methods to provide molecular-level insights into the structure and dynamics of lipid assemblies.

Algorithms↗
PubMed1999

Structural and dynamic differences of the estrogen receptor DNA-binding domain, binding as a dimer and as a monomer to DNA: molecular dynamics simulation studies.

Molecular dynamics (MD) simulations of the estrogen receptor DNA-binding domain (ERDBD) as a dimer in complex with its DNA response element (ERE) show a significant difference in both structure and dynamics, compared to a MD simulation of monomeric ERDBD bound to its half-site response element (EREH). The C-terminal zinc binding domain (ZnII), including a region (helix II) which is in a helical conformation in ERE-(ERDBD)2, is considerably more flexible in EREH-ERDBD than in the dimeric complex. In EREH-ERDBD, all helical hydrogen bonds in helix II are broken and the entire ZnII region is detached from a hydrogen bonding network that in ERE-(ERDBD)2 connects to other parts of the protein as well as to the DNA. The regions that become flexible in EREH-ERDBD are identical to the regions where the NMR solution structure of free ERDBD is poorly ordered. This strongly suggests that dimerisation of ERDBD is required for ordering of the ZnII region and that monomeric binding to DNA is not sufficient for the ordering. This contrasts to the glucocorticoid receptor DNA-binding domain (GRDBD) which has essentially the same mobility (uniform and limited), regardless of whether it is free as a monomer in solution, bound as a monomer to its half-site response element or in a dimeric complex with the full response element. The hydrogen bonding network that connects ZnII with other parts of the protein and to DNA is almost identical in ERDBD and GRDBD. However, in GRDBD there is also a serine (in the N-terminal zinc coordinating region) with a central role in this network, connecting to the ZnII region. This serine is replaced by a glycine in ERDBD and we suggest that this substitution is sufficient for destabilisation of the network, thus leading to a more flexible ZnII region, which becomes ordered first upon forming a complex with another ERDBD and DNA.

Amino Acid Sequence↗
PubMed2001

Effects of phospholipid unsaturation on the membrane/water interface: a molecular simulation study.

Molecular dynamics (MD) simulations of fully hydrated bilayers in the liquid-crystalline state made of 1-palmitoyl-2-oleoyl-phosphatidylcholine (POPC) or 1-palmitoyl-2-elaidoyl-phosphatidylcholine (PEPC) were carried out to investigate the effect of the incorporation of a double bond in the phosphatidylcholine (PC) beta-chain (cis or trans) on the membrane/water interface. The bilayers reached thermal equilibrium after 3 and 1 ns of MD simulations, respectively, and productive runs were carried out for 3 ns for each bilayer. As reference systems, the 1,2-dimyristoyl-phosphatidylcholine (DMPC) bilayer (M. Pasenkiewicz-Gierula, Y. Takaoka, H. Miyagawa, K. Kitamura, and A. Kusumi, 1999, Biophys. J. 76:1228-1240) and DMPC-cholesterol (Chol) bilayer containing 22 mol % Chol (M. Pasenkiewicz-Gierula, T. Róg, K. Kitamura, A. and Kusumi, 2000, Biophys. J. 78:1376-1389) were used. The study shows that at the interface of POPC, PEPC, and DMPC-Chol bilayers, average numbers of PC-water and PC-PC interactions are similar and, respectively, greater and smaller than in the DMPC bilayer. The average area/PC in mono-unsaturated bilayers is approximately 4 A(2) larger than in the DMPC bilayer; nevertheless, a strong correlation was found between a single molecular area (SMA) of a PC and the number of interactions this PC makes; i.e., PCs (either saturated or unsaturated) with the same SMA form similar numbers of intermolecular links. The numbers and corresponding SMAs are distributed about averages pertinent to each bilayer. No significant difference between cis and trans bonds was found.

Binding Sites↗
PubMed2005

Anionic phospholipid interactions with the potassium channel KcsA: simulation studies.

Molecular dynamics (MD) simulations have been used to unmask details of specific interactions of anionic phospholipids with intersubunit binding sites on the surface of the bacterial potassium channel KcsA. Crystallographic data on a diacyl glycerol fragment at this site were used to model phosphatidylethanolamine (PE), or phosphatidylglycerol (PG), or phosphatidic acid (PA) at the intersubunit binding sites. Each of these models of a KcsA-lipid complex was embedded in phosphatidyl choline bilayer and explored in a 20 ns MD simulation. H-bond analysis revealed that in terms of lipid-protein interactions PA > PG >> PE and revealed how anionic lipids (PG and PA) bind to a site provided by two key arginine residues (R(64) and R(89)) at the interface between adjacent subunits. A 27 ns simulation was performed in which KcsA (without any lipids initially modeled at the R(64)/R(89) sites) was embedded in a PE/PG bilayer. There was a progressive specific increase over the course of the simulation in the number of H-bonds of PG with KcsA. Furthermore, two specific PG binding events at R(64)/R(89) sites were observed. The phosphate oxygen atoms of bound PG formed H-bonds to the guanidinium group of R(89), whereas the terminal glycerol H-bonded to R(64). Overall, this study suggests that simulations can help identify and characterize sites for specific lipid interactions on a membrane protein surface.

Anions↗
PubMed2006

X-ray diffraction "fingerprinting" of DNA structure in solution for quantitative evaluation of molecular dynamics simulation.

Solution state x-ray diffraction fingerprinting is demonstrated as a method for experimentally assessing the accuracy of molecular dynamics (MD) simulations. Fourier transforms of coordinate data from MD simulations are used to produce reciprocal space "fingerprints" of atomic pair distance correlations that are characteristic of the ensemble and are the direct numerical analogues of experimental solution x-ray diffraction (SXD). SXD experiments and MD simulations were carried out to test the ability of experiment and simulation to resolve sequence-dependent modifications in helix conformation for B-form DNA. SXD experiments demonstrated that solution-state poly(AT) and poly(A)-poly(T) duplex DNA sequences exist in ensembles close to canonical B-form and B'-form structures, respectively. In contrast, MD simulations analyzed in terms of SXD fingerprints are shown to deviate from experiment, most significantly for poly(A)-poly(T) duplex DNA. Compared with experiment, MD simulation shortcomings were found to include both mismatches in simulated conformer structures and number population within the ensembles. This work demonstrates an experimental approach for quantitatively evaluating MD simulations and other coordinate models to simulate biopolymer structure in solution and suggests opportunities to use solution diffraction data as experimental benchmarks for developing supramolecular force fields optimized for a range of in situ applications.

DNA↗
PubMed2003

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials.

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.

Amino Acid Sequence↗
PubMed2005

Car-Parrinello molecular dynamics simulation of Fe 3+ (aq).

The optimized geometry and energetic properties of Fe(D2O)n 3+ clusters, with n = 4 and 6, have been studied with density-functional theory calculations and the BLYP functional, and the hydration of a single Fe 3+ ion in a periodic box with 32 water molecules at room temperature has been studied with Car-Parrinello molecular dynamics and the same functional. We have compared the results from the CPMD simulation with classical MD simulations, using a flexible SPC-based water model and the same number of water molecules, to evaluate the relative strengths and weaknesses of the two MD methods. The classical MD simulations and the CPMD simulations both give Fe-water distances in good agreement with experiment, but for the intramolecular vibrations, the classical MD yields considerably better absolute frequencies and ion-induced frequency shifts. On the other hand, the CPMD method performs considerably better than the classical MD in describing the intramolecular geometry of the water molecule in the first hydration shell and the average first shell...second shell hydrogen-bond distance. Differences between the two methods are also found with respect to the second-shell water orientations. The effect of the small box size (32 vs 512 water molecules) was evaluated by comparing results from classical simulations using different box sizes; non-negligible effects are found for the ion-water distance and the tilt angles of the water molecules in the second hydration shell and for the O-D stretching vibrational frequencies of the water molecules in the first hydration shell.

Journal Article↗
PubMed2006

Molecular dynamics simulation study of the structural characteristics of water molecules confined in functionalized carbon nanotubes.

Molecular dynamics (MD) simulations were performed to study the structural properties of water molecules confined in functionalized carbon nanotubes (CNTs). Four CNTs, two armchair-type (6, 6), (7, 7) and two zigzag-type (10, 0), (12, 0) CNTs, representing different helicities and different diameters, were chosen and functionalized at their open ends by the hydrophilic -COOH and the hydrophobic -CH3 groups. The structural properties of water molecules inside the functionalized CNTs, including the orientation distributions of dipole moment and O-H bonds, the length of the single-file water chain, and the average number of hydrogen bonds, were analyzed during a process of simulations. MD simulation results in this work showed that the -CH3 functional groups exert little special effects on the structural properties of water molecules. It is mainly due to the relatively small size of the -CH3 group and its hydrophobic nature, which is consistent with hydrophobic CNTs. For CNTs functionalized by -COOH groups, the configurations of -COOH groups, incurvature or excurvature, determine whether water molecules can enter the CNTs. The incurvature or excurvature configurations of -COOH groups are the results of synergy effects of the CNTs' helicity and diameter and control the flow direction of water molecules in CNTs.

Journal Article↗
PubMed1998

Essential dynamics of DNA containing a cis.syn cyclobutane thymine dimer lesion.

Conformational properties of a UV-damaged DNA decamer containing a cis.syn cyclobutane thymine dimer (PD) have been investigated by molecular dynamics (MD) simulations. Results from MD simulations of the damaged decamer DNA show a kink of approximately 21.7 degrees at the PD damaged site and a disruption of H bonding between the 5'-thymine of the PD and its complementary adenine. However, no extra-helical flipping of the 3'-adenine complementary to the PD was observed. Comparison to two undamaged DNA decamers, one with the same sequence and the other with an AT replacing the TT sequence, indicates that these properties are specific to the damaged DNA. Essential dynamics (ED) derived from the MD trajectories of the three DNAs show that the backbone phosphate between the two adenines complementary to the PD of the damaged DNA has considerably larger mobility than the rest of the molecule and occurs only in the damaged DNA. As observed in the crystal structure of T4 endonuclease V in a complex with the damaged DNA, the interaction of the enzyme with the damaged DNA can lead to bending along the flexible joint and to induction of adenine flipping into an extra-helical position. Such motions may play an important role in damage recognition by repair enzymes.

Adenine↗
PubMedDate not supplied

On the dependence of molecular conformation on the type of solvent environment: a molecular dynamics study of cyclosporin A.

The dependence of the conformation of cyclosporin A (CPA), a cyclic undecapeptide with potent immunosuppressive activity, on the type of solvent environment is examined using the computer simulation method of molecular dynamics (MD). Conformational and dynamic properties of CPA in aqueous solution are obtained from MD simulations of a CPA molecule dissolved in a box with water molecules. Corresponding properties of CPA in apolar solution are obtained from MD simulations of CPA in a box with carbontetrachloride. The results of these simulations in H2O and in CCl4 are compared to each other and to those of previous simulations of crystalline CPA and of an isolated CPA molecule. The conformation of the backbone of the cyclic polypeptide is basically independent of the type of solvent. In aqueous solution the beta-pleated sheet is slightly weaker and the gamma-turn is a bit less pronounced than in apolar solution. Side chains may adopt different conformations in different solvents. In apolar solution the hydrophobic side chain of the MeBmt residue is in an extended conformation with its hydroxyl group hydrogen bonded to the backbone carbonyl group. In aqueous solution this hydrophobic side chain folds over the core of the molecule and the mentioned hydrogen bond is broken in favor of hydrogen bonding to water molecules. The conformation obtained from the MD simulation in CCl4 nicely agrees with experimental atom-atom distance data as obtained from nmr experiments in chloroform. In aqueous solution the relaxation of atomic motion tends to be slower than in apolar solution.

Chemical Phenomena↗