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

Cyclosporin A: Conformational Complexity and Chameleonicity

The chameleonic behavior of cyclosporin A (CsA) was investigated through conformational ensembles employing multicanonical molecular dynamics simulations that could sample the cis and trans isomers of N-methylated amino acids; these assessments were conducted in explicit water, dimethyl sulfoxide, acetonitrile, methanol, chloroform, cyclohexane (CHX), and n-hexane (HEX) using AMBER ff03, AMBER10:EHT, AMBER12:EHT, and AMBER14:EHT force fields. The conformational details were discussed employing the free-energy landscapes (FELs) at T = 300 K; it was observed that the experimentally determined structures of CsA were only a part of the conformational space. Comparing the ROESY measurements in CHX-d12 and HEX-d14, the major conformations in those apolar solvents were essentially the same as that in CDCl₃ except for the observation of some sidechain rotamers. The effects of the metal ions on the conformations, including the cis/trans isomerization, were also investigated. Based on the analysis of FELs, it was concluded that the AMBER ff03 force field best described the experimentally derived conformations, indicating that CsA intrinsically formed membrane-permeable conformations and that the metal ions might be the key to the cis/trans isomerization of N-methylated amino acids before binding a partner protein.

Toward managing chronic rejection after lung transplant: the fate and effects of inhaled cyclosporine in a complex environment

The fate and effects of inhaled cyclosporine A (CsA) are considered after deposition on the lung surface. Special emphasis is given to a post-lung transplant environment and to the potential effects of the drug on the various cell types it is expected to encounter. The known stability, metabolism, pharmacokinetics and pharmacodynamics of the drug have been reviewed and discussed in the context of the lung microenvironment. Arguments support the contention that the immuno-inhibitory and anti-inflammatory effects of CsA are not restricted to T-cells. It is likely that pharmacologically effective concentrations of CsA can be sustained in the lungs but due to the complexity of uptake and action, the elucidation of effective posology must ultimately rely on clinical evidence.

Probing the bioactive conformation of an archetypal natural product HDAC inhibitor with conformationally homogeneous triazole-modified cyclic tetrapeptides

Fooling enzymes with mock amides: Analogues of apicidin, a cyclic-tetrapeptide inhibitor of histone deacetylase (HDAC), were designed with a 1,4- or 1,5-disubstituted 1,2,3-triazole in place of a backbone amide bond to fix the bond in question in either a trans-like or a cis-like configuration. Thus, the binding affinity of distinct peptide conformations (see picture) could be probed. One analogue proved in some cases to be superior to apicidin as an HDAC inhibitor.

Structure Refinement of Cyclosporin A in Chloroform by Using RDCs Measured in a Stretched PDMS-Gel

New developments concerning alignment media for apolar solvents like chloroform make it possible to measure anisotropic parameters such as residual dipolar couplings (RDCs) at relatively low concentrations and natural isotopic abundance. As RDCs provide structural restraints with respect to an external coordinate system, long-range structural arrangements of the time-averaged structure can be determined with high precision. The method is demonstrated on the well-studied cyclo-undecapeptide Cyclosporin A (CsA), for which crystal and conventionally derived NMR structures are available. Neither crystal nor NMR structure are consistent with heteronuclear D(CH) RDCs measured in a stretched poly(dimethylsiloxane) gel, and refinement by using the anisotropic parameter results in a highly defined structure with a slightly changed backbone conformation. The applied methods and interpretation of the structural model are discussed.

A study of the binding requirements of calyculin A and dephosphonocalyculin A with PP1, development of a molecular recognition model for the binding interactions of the okadaic acid class of compounds with PP1

The interactions of the okadaic acid class of compounds, with special emphasis on the solution structures of calyculin A and dephosphonocalyculin A with PP1 are reported. After examination of the interactions of all docked structures, a receptor based pharmacophore model for the interactions of the protein phosphatase inhibitors has been developed. Calyculin A or dephosphonocalyculin A can interact with the enzyme in either a manner similar to the reported crystal structure, or in an extended form. The inhibitors require two essential regions interacting with the hydrophobic region and the central metal binding regions of the enzyme. This simplified model is consistent with previously published models of the okadaic acid class of compounds with PP1.

Structure investigation of amphiphilic cyclopeptides in isotropic and anisotropic environments-A model study simulating peptide-membrane interactions

It has been proposed that the membrane allows a much more efficient binding of certain small or medium-sized amphiphilic messenger molecules to their receptor, not only by accumulation of the drug, but also by induction of orientations and conformations that are much more favorable for receptor docking than structures adopted in isotropic phases. A series of eight amphiphilic cyclic peptides containing lipophilic (L-alpha-aminodecanoic acid = Ada, L-alpha-aminohexadecanoic acid = Ahd, Nhdg = N-hexadecylglycine) and hydrophilic (Lys, Asp) amino acids were synthesized and examined by means of NMR spectroscopy and molecular dynamics (MD) simulations in isotropic (CDCl3) and membrane-mimicking anisotropic (SDS/H2O) solvents to study the influence of the environment on their individual conformations. NMR data of cyclo(-Gly1-D-Asp2-Ahd3-Ahd4-Asp5-Gly6+ ++-) (C4), cyclo(-Lys1-D-Pro2-Lys3-Ada4-Pro5-Ada6-) (C5) and cyclo(-Lys1-Pro2-Lys3-Ada4-D-Pro5-Ada6-) (C6) clearly indicate that those compounds are too rigid to perform a conformational change upon transition from an isotropic to an anisotropic environment. On the other hand, the experimental data of cyclo (-Gly1-Asp2-Ahd3-Ahd4-Asp5-Gly6-) (C1), cyclo(-Asp1-Ala2-Nhdg3-Ala4-D-Asp5-) (C7), and cyclo(-D-Asp1-Ala2-Nhdg3-Ala4-Asp5-) (C8) suggest highly flexible unstructured molecules in both environments. However, for cyclo(-Asp1-Asp2-Gly3-Ahd4-Ahd5-Gly6-) (C2) we observed a structure inducing effect of a membrane-like environment. The compound populates three different conformations in SDS/H2O, whereas in CDCI3 no preferred conformation can be detected. cyclo(-D-Asp1-Asp2-Gly3-Ahd4-Ahd5-Gly6-) (C3) clearly exhibits two different conformations with a shifted beta,beta-turn motif in CDCI3 and SDS/H2O solutions. The conformational change could be reproduced in a restraint-free MD simulation using the biphasic membrane mimetic CCl4/H2O. Our results give clear evidence that membrane interactions may not only lead to structure inductions, but can also induce major conformational changes in compounds already exhibiting a defined structure in isotropic solution.

Analysis of the dynamics of rhizomucor miehei lipase at different temperatures

The dynamics of Rhizomucor miehei lipase has been studied by molecular dynamics simulations at temperatures ranging from 200-500K. Simulations carried out in periodic boundary conditions and using explicit water molecules were performed for 400 ps at each temperature. Our results indicate that conformational changes and internal motions in the protein are significantly influenced by the temperature increase. With increasing temperature, the number of internal hydrogen bonds decreases, while surface accessibility, radius of gyration and the number of residues in random coil conformation increase. In the temperature range studied, the motions can be described in a low dimensional subspace, whose dimensionality decreases with increasing temperature. Approximately 80% of the total motion is described by the first (i) 80 eigenvectors at T=200K, (ii) 30 eigenvectors at T=300K and (iii) 10 eigenvectors at T=400K. At high temperature, the alpha-helix covering the active site in the native Rhizomucor miehei lipase, the helix at which end the active site is located, and in particular, the loop (Gly35-Lys50) show extensive flexibility.

The solution structures of calyculin A and dephosphonocalyculin A by NMR

The NMR solution structure of calyculin A (1) in chloroform exhibits intramolecular interactions, resembling the original crystal structure. In methanol, calyculin A has the hydrogen bonding moieties solvent exposed. Dephosphonocalyculin A in chloroform resembles calyculin A in chloroform and the crystal structure of calyculin A. Dephosphonocalyculin A in methanol resembles calyculin A in methanol.

Dynamics of proteins in different solvent systems: analysis of essential motion in lipases

We have investigated the effect of different solvents on the dynamics of Rhizomucor miehei lipase. Molecular dynamics simulations were performed in water, methyl hexanoate, and cyclohexane. Analysis of the 400-ps trajectories showed that the solvent has a pronounced effect on the geometrical properties of the protein. The radius of gyration and total accessibility surface decrease in organic solvents, whereas the number of hydrogen bonds increases. The essential motions of the protein in different solvents can be described in a low-dimensional "essential subspace," and the dynamic behavior in this subspace correlates with the polarity of the solvent. Methyl hexanoate, which is a substrate for R. miehei lipase, significantly increases the fluctuations in the active-site loop. During the simulation, a methyl hexanoate entered the active-site groove. This observation provides insight into the possible docking mechanism of the substrate.

New approaches in molecular structure prediction

In the past years, much effort has been put on the development of new methodologies and algorithms for the prediction of protein secondary and tertiary structures from (sequence) data; this is reviewed in detail. New approaches for these predictions such as neural network methods, genetic algorithms, machine learning, and graph theoretical methods are discussed. Secondary structure prediction algorithms were improved mostly by considering families of related proteins; however, for the reliable tertiary structure modeling of proteins, knowledge-based techniques are still preferred. Methods and examples with more or less successful results are described. Also, programs and parameterizations for energy minimisations, molecular dynamics, and electrostatic interactions have been improved, especially with respect to their former limits of applicability. Other topics discussed in this review include the use of traditional and on-line databases, the docking problem and surface properties of biomolecules, packing of protein cores, de novo design and protein engineering, prediction of membrane protein structures, the verification and reliability of model structures, and progress made with currently available software and computer hardware. In summary, the prediction of the structure, function, and other properties of a protein is still possible only within limits, but these limits continue to be moved.

Modeling conformational changes in cyclosporin A

NMR and X-ray structures for the immunosuppressant cyclosporin A (CsA) reveal a remarkable difference between the unbound (free) conformation in organic solvents and the conformation bound to cyclophilin. We have performed computer simulations of the molecular dynamics of CsA under a variety of conditions and confirmed the stability of these two conformations at room temperature in water and in vacuum. However, when the free conformation was modeled in vacuum at 600 K, a transition pathway leading to the bound conformation was observed. This involved a change in the cis MeLeu-9 peptide bond to a trans conformation and the movement of the side chains forming the dominant hydrophobic cluster (residues MeBmt-1, MeLeu-4, MeLeu-6, and MeLeu-10) to the opposite side of the plane formed by the backbone atoms in the molecular ring. The final conformation had a backbone RMS deviation from the bound conformation of 0.53 A and was as stable in dynamics simulations as the bound conformation. Our calculations allowed us to make a detailed analysis of a transition pathway between the free and the bound conformations of CsA and to identify two distinct regions of coordinated movement in CsA, both of which underwent transitions independently.

The effect of environment on the stability of an integral membrane helix: molecular dynamics simulations of surfactant protein C in chloroform, methanol and water

A series of three molecular dynamics simulations at 300 K in explicit solvent environments of chloroform, methanol and water has been performed on the pulmonary surfactant lipoprotein, SP-C, comprising several consecutive valine residues in order to investigate the stability of the alpha-helical conformation. Two additional simulations were performed on truncated SP-C with a five-residue N-terminal deletion at 300 K and 500 K in water, the high temperature run in order to increase the rate of peptide denaturation. Indications of destabilization appear in chloroform during 1 ns while the SP-C alpha-helix is remarkably stable during 1 ns in methanol and water. In particular the polyvalyl part comprising residues Val15 to Val21 remains intact even at elevated temperature, and the valines do not disrupt the alpha-helical conformation. The valyl-rotamer sampling is partly restricted. Unfolding takes place successively along the primary sequence starting from the C-terminal end. Factors affecting polypeptide stability in molecular dynamics simulations are addressed. The intrinsic helix-forming tendency of valine residues and its dependence on the sequence context, and the role of the solvent environment in stabilizing or destabilizing an alpha-helical fold, are discussed.

Conformation and molecular dynamics calculations on uteroglobin fragment 18-47

The conformational properties of fragment 18-47 of rabbit uteroglobin in aqueous solution containing SDS micelles were investigated by two-dimensional nmr spectroscopy and molecular dynamics calculations. The fragment comprises helices II and III and the beta-turn connecting the two helices. The nmr results and nmr-restrained molecular dynamics calculations showed that in the isolated fragment the elements of secondary structure present in the intact protein are preserved only in part. Specifically, a well-defined alpha-helix was found in the sequence 33-44, corresponding to helix III of uteroglobin, while the regions of helix II and beta-turn are characterized by high flexibility in the fragment.

Conformational changes in the cyclic undecapeptide cyclosporin induced by interaction with metal ions. An FTIR study

Infra-red spectra have been measured for the cyclic undecapeptide cyclosporin A (CsA) and three analogues CsC, CsD and CsH in acetonitrile and in the presence of 10:1 molar excess of Mg2+, Ca2+, Na+ or Li+ in the same solvent. Interaction with each of these ions is suggested by marked changes in band positions over the amide I region (1600-1700 cm-1). The formation of complexes of cyclosporin with calcium and magnesium ions is indicated by the presence of C = O stretching bands well outside the range normally expected for the amide I absorptions of free peptides. Although they share this characteristic, the spectra indicate that the mode and/or strength of Ca2+ binding is quite different from that of Mg2+ binding. In contrast, the two monovalent ions interact with CsA, CsC and CsD to yield spectra that are very similar to one another. The spectra are consistent with binding of the monovalent ions simultaneously to several carbonyl groups of the loop structure.

Improved molecular dynamics simulations for the determination of peptide structures

In this article a few methods or modifications proven to be useful in the conformational examination of peptides and related molecules by molecular dynamics are illustrated. The first is the explicit use of organic solvents in the simulations. For many cases such solvents are appropriate since the nmr measurements (or other experimental observations) were carried out in the same solvent. Here, the use of dimethylsulfoxide and chloroform in molecular dynamics is described, with some advantages of the use of these solvents high-lighted. A constant allowing for the scaling of the nonbonded interactions of the force field, an idea previously employed in distance geometry and simulated annealing, has been implemented. The usefulness of this method is that when the nonbonded term is turned to zero, atoms can pass through each other, while the connectivity of the molecule is maintained. It will be shown that such simulations, if a sufficient driving force is present (i.e., nuclear Overhauser effects restraints), can produce the correct stereoconfiguration (i.e., chiral center) as well as configurational isomer (i.e., cis/trans isomers). Lastly, a penalty term for coupling constants directly related to the Karplus curve has been implemented into the potential energy force field. The advantages of this method over the commonly used dihedral angle restraining are discussed. In particular, it is shown that with more than one coupling constant about a dihedral angle a great reduction of the allowed conformational space is obtained.

Complexation and medium effects on the conformation of cyclosporin A studied by NMR spectroscopy and molecular dynamics calculations

The conformation of cyclosporin A was investigated in different environments. Whereas there was no distinct difference between the structure of CSA in CDCl3 and THF, the addition of LiCl caused a conformational change. Also, using polar solvents (DMSO, water) induced conformational changes in the backbone as well as the side chains.