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.

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.

Calculating three-dimensional molecular structure from atom-atom distance information: cyclosporin A.

In recent years methods for deriving spatial molecular structure from atom-atom distance information have gained in importance due to the emergence of two-dimensional nuclear magnetic resonance (n.m.r) techniques, which make it possible to obtain such distance information for polypeptides, small proteins, sugars, and DNA fragments in solution. Distance geometry (DG) and restrained molecular dynamics (MD) refinement are applied to a cyclic polypeptide, the immunosuppressive drug cyclosporin A, and the results are compared. Two different procedures, DG followed by restrained MD, and straightforward restrained MD starting from the X-ray structure, both lead to a unique conformation that satisfies the 58 experimentally determined distance constraints. The results nicely show the relative merits of DG and restrained MD techniques for determining spatial molecular structure from distance information.

Molecular dynamics simulations of cyclosporin A: the crystal structure and dynamic modelling of a structure in apolar solution based on NMR data.

The conformation of the immunosuppressive drug cyclosporin A (CPA), both in apolar solution and in crystalline state, has been studied by computer simulation techniques. Three molecular dynamics (MD) simulations have been performed: one modelling the crystal structure and two modelling the structure in apolar solution, using a restrained MD approach in which data from nuclear magnetic resonance (NMR) and infrared (IR) spectroscopy are taken into account. The simulation of the crystalline state (MDC) concerns a system of 4 unit cells containing 16 cyclosporin A molecules and 22 water molecules, which is simulated using crystalline periodic boundary conditions. The simulations modelling the apolar solvent conformation (MDS) concern one isolated cyclosporin A molecule. In these simulations an extra term in the interatomic potential function is used, which forces the molecule to satisfy a set of 57 atom-atom distance constraints originating from nuclear Overhauser effects (NOEs) obtained from NMR spectroscopy and one distance constraint deduced from IR spectroscopy. From a comparison of the results of the crystal simulation to those of the X-ray experiment in terms of structure, atomic fluctuations, hydrogen bond pattern, etc., it is concluded that the force field that is used yields an adequate representation of crystalline cyclosporin A. Secondly, it is shown that the dynamic modelling technique that is used to obtain a structure in a polar solution from NMR distance information works well. Starting from initial conformations which have a root mean square difference of 0.14 nm both distance restrained MD simulations converge to the same final solution structure. A comparison of the crystal structure of cyclosporin A and the one in apolar solution shows that there are significant differences. The overall difference in atomic positions is 0.09 nm for the C alpha atoms and 0.17 nm for all atoms. In apolar solution, the molecule is slightly more bent and the side chains of 1 MeBmt and 10 MeLeu adopt a different conformation.

Conformational analysis of a cyclic thymopoietin-analogue by 1H n.m.r. spectroscopy and restrained molecular dynamics simulations.

The internuclear distances of the cyclic thymopoietin derivative c[D-Val-Tyr-Arg-Lys-Glu] have been determined using two-dimensional nuclear Overhauser n.m.r. spectroscopy. These distances are used as constraints for a restrained Molecular Dynamics (MD) simulation. The two starting structures used for the calculations consist of a beta and gamma turn for model 1 and two gamma turns for model 2. The rms difference in atomic positions of the two conformations is 0.242 nm. They converge during the restrained MD simulation to the same final structure. The positional rms difference of the time averaged (5-14 ps) conformations is 0.011 nm. The hydrogen bond pattern is similar to that of model 1, but in addition we find three more gamma turns. The vicinal NH-C alpha H couplings agree well with those calculated from the time averaged structures.