Molecular dynamics as a tool to help design mutants

A 200 ps MD trajectory of wild type PCI and a 120 ps one for the Pro36Gly putative mutant are studied and compared with the structure of PCI in its complex with carboxypeptidase A (CPA). It is first established that the structures of PCI from X-ray and from MD simulation are essentially equal. Thereafter, data from the PCI-CPA and average MD structures together with available biochemical information are used to identify possible structural factors that may determine the inhibitory power of PCI. These structural determinants are used to analyze the mutant structure. The fold of the mutant protein shows a large degree of stability. The N-terminal tail in PCIm differs from the X-ray structure as it does in PCIw, while the mutant's C-terminal tail (which is the primary binding site with CPA) and residues 13–17 present deviations. Differences in fluctuation patterns exist between PCIm and PCIw in residues 2–4 (the N-terminal tail), 13–17, 22–23, 28–81 (the secondary contact site with CPA), and 37–38 (the C-terminal tail); the latter region is rigidified in PCIm. Results show that the MD method is able to sense long-range as well as local perturbative effects produced by amino-acid substitutions in flexible regions of this protein. The simulations suggest that the conformation of the C-terminal tail is less favorable for interaction with the target protein in the mutant than it is in the wild type protein. The Pro-36-Gly mutant is predicted to be a less potent inhibitor.Abbreviations CPA carboxypeptidase A - MD molecular dynamics - NIS non-inertial solvent - PCI potato carboxypeptidase A inhibitor - PCIm mutated inhibitor - PCIw wild inhibitor     

A molecular dynamics study of a model built Pro-36-Gly mutant derived from the potato carboxypeptidase A inhibitor protein

A 120ps molecular dynamics (MD) trajectory was calculated and analyzed for a putative Pro-36-Gly mutant of the potato carboxypeptidase A (CPA) protein inhibitor (PCIm). The mutant protein's fold shows a large degree of stability, judged from its low alpha-carbon r.m.s. deviation from the X-ray structure of the wild type PCI (PCIw). The N-terminal tail of PCIm differs slightly less from the X-ray structure than it does in PCIw, while the mutant's C-terminal tail (the primary contact site with CPA) and residues 13-17 present deviations as they approach each other. Differences in fluctuation pattern exist between PCIm and PCIw in residues 2-4 (the N-terminal tail), 13-17, 22-23, 28-31 (the secondary contact site with CPA) and 37-38 (the C-terminal tail); the latter region is rigidified in PCIm. Results show that the MD method is able to sense local perturbative effects produced by amino acid substitutions in flexible regions of protein molecules. The simulation suggests that the conformation of the C-terminal tail is less favorable for interaction with the target protein in the mutant than it is in the wild type protein. The Pro-36-Gly mutant is predicted to be a less potent inhibitor.     

Stability and fluctuations of the potato carboxypeptidase A protein inhibitor fold: a molecular dynamics study

A 120ps non-inertial solvent (NIS) molecular dynamics (MD) trajectory of the potato carboxypeptidase A protein inhibitor (PCI) was calculated and analyzed. It is shown that, in spite of a very low content of regular secondary structure, the PCI fold has a large degree of stability, judged from the fairly good agreement between the average MD and X-ray structures. The N-terminal and C-terminal regions behave differently, both in their isoatomic positional shifts with respect to the X-ray structure, and in atomic fluctuation pattern. Positional shifts up to 9A are detected in the exposed N-terminal region as it folds back on the inhibitor's core. This large deviation is most likely caused by the absence of the receptor protein or by the lack of supporting solvent molecules. In contrast, the C-terminal region, which is the primary contact site with the enzyme, has an average structure similar to the X-ray conformation; this feature is probably due to a hydrogen bond network to the central core of PCI. The C-terminal tail shows larger fluctuations than the core. The secondary contact site retains its structure in this simulation. The results evidence an intrinsically stable PCI fold which favors a spatially well defined, fairly flexible, structuration of the primary and secondary contact sites that optimizes PCI's interaction with its target enzyme.     

NMR solution structure of a novel hirudin variant HM2, N-terminal 1-47 and N64 → V + G mutant

The 64 amino acid hirudin-like peptide HM2 (Hirudinaria manillensis) is one of the agents known to specifically block the blood-clotting enzyme thrombin, and therefore is used as a potential pharmacological tool for the treatment of arterial and venous thrombosis. This peptide and its derivatives provide a new set of probes for studies aimed at elucidating the structural basis of the inhibition of α-thrombin. We used 581, 699, and 492 nmr-derived constraints respectively in a protocol employing simulated annealing, followed by restrained molecular dynamics and restrained energy minimization to derive the three-dimensional structures of HM2 and its mutants the HM2(V + G) and the HM2(1–47). HM2 consists of a well-defined core region of two double-stranded β-sheet and a disordered C-terminus. These features are shared by other members of the hirudin family. The same type of folding has also been observed for recombinant hirudins whose structure has been determined in solution by nmr spectroscopy and in the structure of the complex hirudin-thrombin determined by x-ray diffraction. Molecular dynamics (MD) simulation methods were applied in the study of the structural and dynamic fluctuation properties of the hirudin derivatives solvated by 1625 and 1276 water molecules with periodic boundary conditions for HM2 and HM2(1–47), respectively. Trajectories of 100 and 50 ps for the two unconstrained systems were generated at constant temperature and pressure. Analysis of the MD simulation shows that the structure of the peptide core is fairly rigid and stable in itself while the conformation of the C-terminal tail, which is involved in the inhibitory mechanism of thrombin, fluctuates and appears as a disordered region. © 1997 John Wiley & Sons, Inc. Biopoly 41: 731–749, 1997     

On the Entropic and Hydrophobic Properties Involved in the Inhibitory Mechanism of Carboxypeptidase A by its Natural Inhibitor from Potato

The inhibition of carboxypeptidase A (CPA) by its natural inhibitor from potato (PCI) has been widely analysed with theoretical and experimental methods. Several mutants of PCI have been obtained in order to study the physico-chemical properties related to the inhibition. Point mutations were performed in the C-tail of PCI given its fundamental role in the inhibition. The inhibition constant and the dissociation free energy of the complexes PCI-CPA was experimentally obtained for each mutant. The mutants were divided in two sets, those where the mutation was intrinsically affecting the conformation of the PCI C-tail, and those where the mutation affected the interaction between PCI and CPA. The crystallographic structure of PCI, as found in its complex with bovine carboxypeptidase A, was used to model the structure of these mutants. Two theoretical approaches were performed to explain both sets of experimental results: 1) study of the structural features of wt PCI and mutant forms by molecular dynamics (MD) simulation, and 2) modelling of the interaction of the C-tail of PCI with CPA. The first approach provides an explanation of the observed behaviour of the mutants of PCI, if the hypothesis is made of a direct relationship between the entropy of inhibition and the mobility of the C-tail of PCI. For the second set of mutants, the experimentally measured dissociation energies for the complexes PCI- CPA can be related to the theoretically estimated exposure to the solvent of the side chain of the mutated residue in the complex. In the case of the double mutation G35P+P36G, the importance of the main chain hydrogen bond between Gly 35 and Ala26, anchoring the C-tail to the core of PCI, as predicted by the MD simulations, was also supported by the experimental result. The agreement between the theoretical approaches and the experimental results shows the appropriateness of our hypotheses and also the relevance of such a combined effort of experimental and computational molecular biology in protein engineering.     

Structure and fluctuations of bacteriophage T4 glutaredoxin modelled by molecular dynamics

A 120 ps molecular dynamics (MD) trajectory for bacteriophage T4 glutaredoxin was calculated including non-inertial solvent effects. The potential energy attains an equilibrated regime after the first 20 ps. The r.m.s. difference of all non-hydrogen atoms between X-ray and average MD structures for the regular secondary structure is 0.99A which shows that the MD simulation reproduces the essentials of the structure with high accuracy. Loop displacements are detected, shown by the larger full structure all non-hydrogen atom r.m.s. difference of 1.2A. The fluctuation pattern derived from MD agrees fairly well with that derived from X-ray isotropic temperature factors. The active site is a stable structural region in this MD modellization. Structural changes are put in context with the protein's function.     

Structural analysis of substance P using molecular dynamics and NMR spectroscopy.

The present work is a combined structural study, using Nuclear Magnetic Resonance (NMR) and Molecular Dynamics(MD), of the amidated and the free acid forms of substance P in water and methanol. The results obtained using both approaches were compared in order to characterize the structural features of both peptides in solution. From the NMR experiments it was derived that the free acid form adopts an extended conformation at the N-terminus and a helical conformation at the C-terminal segment of the peptide in both water and methanol; these structural features are in qualitative agreement with the results of the MD simulations. No significant differences in behavior were observed between the amidated and the free acid forms of the peptide in the simulations and in the experiments carried out in water, suggesting that the different activities of these analogs are due to their different mode of interaction with the receptor rather than to their structural preferences. Finally, we propose that the structure of substance P can be partially inferred from its sequence due to the presence of a Pro-X-Pro motif on the N-terminus and a Gly-Leu sequence on the C-terminus.     

Structural dynamics in the C terminal domain homolog of orange carotenoid Protein reveals residues critical for carotenoid uptake

The structural features enabling carotenoid translocation between molecular entities in nature is poorly understood. Here, we present the three-dimensional X-ray structure of an expanded oligomeric state of the C-terminal domain homolog (CTDH) of the orange carotenoid protein, a key water-soluble protein in cyanobacterial photosynthetic photo-protection, at 2.9 Å resolution. This protein binds a canthaxanthin carotenoid ligand and undergoes structural reorganization at the dimeric level, which facilitates cargo uptake and delivery. The structure displays heterogeneity revealing the dynamic nature of its C-terminal tail (CTT). Molecular dynamics (MD) simulations based on the CTDH structures identified specific residues that govern the dimeric transition mechanism. Mutagenesis based on the crystal structure and these MD simulations then confirmed that these specific residues within the CTT are critical for carotenoid uptake, encapsulation and delivery processes. We present a mechanism that can be applied to other systems that require cargo uptake.     

A Molecular Dynamics Investigation of Vinculin Activation

Vinculin activation plays a critical role in focal adhesion initiation and formation. In its native state, vinculin is in an autoinhibitory conformation in which domain 1 prevents interaction of the vinculin tail domain with actin by steric hindrance. Once activated, vinculin is able to interact with both actin and talin. Several hypotheses have been put forth addressing the mechanisms of vinculin activation. One set of studies suggests that vinculin interaction with talin is sufficient to cause activation, whereas another set of studies suggests that a simultaneous interaction with several binding partners is necessary to achieve vinculin activation. Using molecular-dynamics (MD) simulations, we investigate the mechanisms of vinculin activation and suggest both a trajectory of conformational changes leading to vinculin activation, and key structural features that are likely involved in stabilizing the autoinhibited conformation. Assuming that the simultaneous interaction of vinculin with both actin and talin causes a stretching force on vinculin, and that vinculin activation results from a removal of steric hindrance blocking the actin-binding sites, we simulate with MD the stretching and activation of vinculin. The MD simulations are further confirmed by normal-mode analysis and simulation after residue modification. Taken together, the results of these simulations suggest that bending of the vinculin-binding-site region in vinculin away from the vinculin tail is the likely trajectory of vinculin activation.     

Molecular dynamics and free energy studies on the carboxypeptidases complexed with peptide/small molecular inhibitor: mechanism for drug resistance

As one potent plant protease inhibitor, potato carboxypeptidase inhibitor (PCI) can competitively inhibit insect digestive metallocarboxypeptidases (MCPs) through interfering with its digestive system that causes amino acid deficiencies and leading to serious developmental delay and mortality. However, this effective biological pest control is significantly impaired by the PCI-resistant insect MCPs. Therefore, deep understanding of the resistant mechanism of insect MCPs is particularly necessary for designing new durable pest control regimen and developing effective pesticides. In this study, the binding of PCI and small molecular inhibitor THI to insect PCI-sensitive/-resistant MCPs and human MCP was investigated by docking, molecular dynamics (MD) simulations and thermodynamic analysis. The structural analysis from MD simulations indicates that the PCI-resistant mechanism of CPBHz is mainly dominated by the Trp277A, which changes the conformation of β8-α9 loop and therefore narrow the access to the active site of CPBHz, prohibiting the entrance of the C termini tail of PCI. Additionally, the insertion of Gly247A weakens the stabilization of CPBHz and PCI through disrupting the hydrogen bond formation with its surrounding residues. Furthermore, the predicted binding free energies gives explanation of structure affinity relationship of PCI and THI with MCPs and suggest that the electrostatic energy is the main contribution term affecting the difference in binding affinities. Finally, the decomposition analysis of binding free energies infers that the key residues Glu72, Arg127, Ile247/Leu247 and Glu270 are critical for the binding of PCI/THI to MCPs.     

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.     

Atomic crystal and molecular dynamics simulation structures of human carbonic anhydrase II: insights into the proton transfer mechanism

Human carbonic anhydrase II (HCA II) is a zinc-metalloenzyme that catalyzes the reversible interconversion of CO2 and HCO3-. The rate-limiting step of this catalysis is the transfer of a proton between the Zn-bound solvent molecule and residue His64. In order to fully characterize the active site structural features implicated in the proton transfer mechanism, the refined X-ray crystal structure of uncomplexed wild type HCA II to 1.05 A resolution with an Rcryst value of 12.0% and an Rfree value of 15.1% has been elucidated. This structure provides strong clues as to the pathway of the intramolecular proton transfer between the Zn-bound solvent and His64. The structure emphasizes the role of the solvent network, the unique positioning of solvent molecule W2, and the significance of the dual conformation of His64 in the active site. The structure is compared with molecular dynamics (MD) simulation calculations of the Zn-bound hydroxyl/His64+ (charged) and the Zn-bound water/His64 (uncharged) HCA II states. A comparison of the crystallographic anisotropic atomic thermal parameters and MD simulation root-mean-square fluctuation values show excellent agreement in the atomic motion observed between the two methods. It is also interesting that the observed active site solvent positions in the crystal structure are also the most probable positions of the solvent during the MD simulations. On the basis of the comparative study of the MD simulation results, the HCA II crystal structure observed is most likely in the Zn-bound water/His64 state. This conclusion is based on the following observations: His64 is mainly (80%) orientated in an inward conformation; electron density omit maps infer that His64 is not charged in an either inward or outward conformation; and the Zn-bound solvent is most likely a water molecule.     

Mapping multiple potential ATP binding sites on the matrix side of the bovine ADP/ATP carrier by the combined use of MD simulation and docking

The mitochondrial adenosine diphosphate/adenosine triphosphate (ADP/ATP) carrier-AAC-was crystallized in complex with its specific inhibitor carboxyatractyloside (CATR). The protein consists of a six-transmembrane helix bundle that defines the nucleotide translocation pathway, which is closed towards the matrix side due to sharp kinks in the odd-numbered helices. In this paper, we describe the interaction between the matrix side of the AAC transporter and the ATP(4-) molecule using carrier structures obtained through classical molecular dynamics simulation (MD) and a protein-ligand docking procedure. Fifteen structures were extracted from a previously published MD trajectory through clustering analysis, and 50 docking runs were carried out for each carrier conformation, for a total of 750 runs ("MD docking"). The results were compared to those from 750 docking runs performed on the X-ray structure ("X docking"). The docking procedure indicated the presence of a single interaction site in the X-ray structure that was conserved in the structures extracted from the MD trajectory. MD docking showed the presence of a second binding site that was not found in the X docking. The interaction strategy between the AAC transporter and the ATP(4-) molecule was analyzed by investigating the composition and 3D arrangement of the interaction pockets, together with the orientations of the substrate inside them. A relationship between sequence repeats and the ATP(4-) binding sites in the AAC carrier structure is proposed.     

Effects of different force fields on the structural character of α synuclein β-hairpin peptide (35–56) in aqueous environment

The hallmark of Parkinson’s disease (PD) is the intracellular protein aggregation forming Lewy Bodies (LB) and Lewy neuritis which comprise mostly of a protein, alpha synuclein (α-syn). Molecular dynamics (MD) simulation methods can augment experimental techniques to understand misfolding and aggregation pathways with atomistic resolution. The quality of MD simulations for proteins and peptides depends greatly on the accuracy of empirical force fields. The aim of this work is to investigate the effects of different force fields on the structural character of β hairpin fragment of α-syn (residues 35–56) peptide in aqueous solution. Six independent MD simulations are done in explicit solvent using, AMBER03, AMBER99SB, GROMOS96 43A1, GROMOS96 53A6, OPLS-AA, and CHARMM27 force fields with CMAP corrections. The performance of each force field is assessed from several structural parameters such as root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent accessible surface area (SASA), formation of β-turn, the stability of folded β-hairpin structure, and the favourable conformations obtained for different force fields. In this study, CMAP correction of CHARMM27 force field is found to overestimate the helical conformation, while GROMOS96 53A6 is found to most successfully capture the conformational dynamics of α-syn β-hairpin fragment as elicited from NMR.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

Conformational interconversion in compstatin probed with molecular dynamics simulations

Compstatin is a 13-residue cyclic peptide that has the potential to become a therapeutic agent against unregulated complement activation. In our effort to understand the structural and dynamic characteristics of compstatin that form the basis for rational and combinatorial optimization of structure and activity, we performed 1-ns molecular dynamics (MD) simulations. We used as input in the MD simulations the ensemble of 21 lowest energy NMR structures, the average minimized structure, and a global optimization structure. At the end of the MD simulations we identified five conformations, with populations ranging between 9% and 44%. These conformations are as follows: 1) coil with alphaR-alphaR beta-turn, as was the conformation of the initial ensemble of NMR structures; 2) beta-hairpin with epsilon-alphaR beta-turn; 3) beta-hairpin with alphaR-alphaR beta-turn; 4) beta-hairpin with alphaR-beta beta-turn; and 5) alpha-helical. Conformational switch was possible with small amplitude backbone motions of the order of 0.1-0.4 A and free energy barrier crossing of 2-11 kcal/mol. All of the 21 MD structures corresponding to the NMR ensemble possessed a beta-turn, with 14 structures retaining the alphaR-alphaR beta-turn type, but the average minimized structure and the global optimization structures were converted to alpha-helical conformations. Overall, the MD simulations have aided to gain insight into the conformational space sampled by compstatin and have provided a measure of conformational interconversion. The calculated conformers will be useful as structural and possibly dynamic templates for optimization in the design of compstatin using structure-activity relations (SAR) or dynamics-activity relations (DAR).     

Structure of a novel leech carboxypeptidase inhibitor determined free in solution and in complex with human carboxypeptidase A2

Leech carboxypeptidase inhibitor (LCI) is a novel protein inhibitor present in the medicinal leech Hirudo medicinalis. The structures of LCI free and bound to carboxypeptidase A2 (CPA2)have been determined by NMR and X-ray crystallography, respectively. The LCI structure defines a new protein motif that comprises a five-stranded antiparallel beta-sheet and one short alpha-helix. This structure is preserved in the complex with human CPA2 in the X-ray structure, where the contact regions between the inhibitor and the protease are defined. The C-terminal tail of LCI becomes rigid upon binding the protease as shown in the NMR relaxation studies, and it interacts with the carboxypeptidase in a substrate-like manner. The homology between the C-terminal tails of LCI and the potato carboxypeptidase inhibitor represents a striking example of convergent evolution dictated by the target protease. These new structures are of biotechnological interest since they could elucidate the control mechanism of metallo-carboxypeptidases and could be used as lead compounds for the search of fibrinolytic drugs.     

Assessment of the molecular dynamics structure of DNA in solution based on calculated and observed NMR NOESY volumes and dihedral angles from scalar coupling constants

To assess the accuracy of the molecular dynamics (MD) models of nucleic acids, a detailed comparison between MD-calculated and NMR-observed indices of the dynamical structure of DNA in solution has been carried out. The specific focus of our comparison is the oligonucleotide duplex, d(CGCGAATTCGCG)(2), for which considerable structural data have been obtained from crystallography and NMR spectroscopy. An MD model for the structure of d(CGCGAATTCGCG)(2) in solution, based on the AMBER force field, has been extended with a 14 ns trajectory. New NMR data for this sequence have been obtained in order to allow a detailed and critical comparison between the calculated and observed parameters. Observable two-dimensional (2D) nuclear Overhauser effect spectroscopy (NOESY) volumes and scalar coupling constants were back-calculated from the MD trajectory and compared with the corresponding NMR data. The comparison of these results indicate that the MD model is in generally good agreement with the NMR data, and shows closer accord with experiment than back-calculations based on the crystal structure of d(CGCGAATTCGCG)(2) or the canonical A or B forms of the sequence. The NMR parameters are not particularly sensitive to the known deficiency in the AMBER MD model, which is a tendency toward undertwisting of the double helix when the parm.94 force field is used. The MD results are also compared with a new determination of the solution structure of d(CGCGAATTCGCG)(2) using NMR dipolar coupling data.     

EPR spectra and molecular dynamics agree that the nucleotide pocket of myosin V is closed and that it opens on binding actin

We have used EPR spectroscopy and computational modeling of nucleotide-analog spin probes to investigate conformational changes at the nucleotide site of myosin V. We find that, in the absence of actin, the mobility of a spin-labeled diphosphate analog [spin-labeled ADP (SLADP)] bound at the active site is strongly hindered, suggesting a closed nucleotide pocket. The mobility of the analog increases when the MV·SLADP complex (MV = myosin V) binds to actin, implying an opening of the active site in the A·MV·SLADP complex (A = actin). The probe mobilities are similar to those seen with myosin II, despite the fact that myosin V has dramatically altered kinetics. Molecular dynamics (MD) simulation was used to understand the EPR spectra in terms of the X-ray database. The X-ray structure of MV·ADP·BeFx shows a closed nucleotide site and has been proposed to be the detached state. The MV·ADP structure shows an open nucleotide site and has been proposed to be the A·MV·ADP state at the end of the working powerstroke. MD simulation of SLADP docked in the closed conformation gave a probe mobility comparable to that seen in the EPR spectrum of the MV·SLADP complex. The simulation of the open conformation gave a probe mobility that was 35-40° greater than that observed experimentally for the A·MV·SLADP state. Thus, EPR, X-ray diffraction, and computational analysis support the closed conformation as a myosin V state that is detached from actin. The MD results indicate that the MV·ADP crystal structure, which may correspond to the strained actin-bound post-powerstroke conformation resulting from head-head interaction in the dimeric processive motor, is superopened.     

Conformational transitions using molecular dynamics with minimum biasing

The molecular dynamics algorithm (MD), which simulates intramolecular motions on the subnanosecond timescale, has been modified to allow the investigation of slow conformational transitions that do not necessarily occur spontaneously in MD simulations. The method is designated CONTRA MD (CONformational TRAnsitions by Molecular Dynamics with minimum biasing). The method requires the prior definition of a single conformational variable that is required to vary monotonically from an initial conformation to a final target conformation. The simulation is broken up into a series of short free MD segments, and we determine, after each segment of MD, whether or not the system has evolved toward the final conformation. Those segments that do not move the system in that direction are deleted. Those that do move it toward the final conformation are patched together sequentially to generate a single representative trajectory along the transition pathway. The CONTRA MD method is demonstrated first by application to the simultaneous C2′-endo to C3′-endo repucker and anti to syn N-glycosidic torsion transitions in 2′-deoxyadenosine and then to the large-scale bending in phenylalanine transfer RNA. © 1993 John Wiley & Sons, Inc.