Dynamics of internal water in fatty acid binding protein: computer simulations and comparison with experiments

Multiple molecular dynamics (MD) simulations of fully solvated rat intestinal fatty acid binding protein (I-FABP) were conducted to investigate the dynamics of internal water molecules. Although the long time average of the number of internal water molecules in I-FABP is 22 as shown by the X-ray crystal structure, MD simulations predict large variations in the instantaneous number of internal water molecules on the nanosecond time scale. The computational model employed predicts that w135 (internal) and w217 (located on the protein surface) may be the water molecules with long residence times observed in previously reported magnetic relaxation dispersion studies. The average residence time of approximately 20 internal water molecules occupying the fatty acid binding cavity is estimated to be between 0.6 and 2.0 nanoseconds. Exchange of internal water in I-FABP appears to occur almost exclusively through the interface of beta-strands EF with the rest of the protein, which has significant implications for the pathways of the fatty acid entry and exit from the binding cavity. Proteins 2001;43:65-72.     

Structure and dynamics of the fatty acid binding cavity in apo rat intestinal fatty acid binding protein.

The structure and dynamics of the fatty acid binding cavity in I-FABP (rat intestinal fatty acid binding protein) were analyzed. In the crystal structure of apo I-FABP, the probe occupied cavity volume and surface are 539+/-8 A3 and 428 A2, respectively (1.4 A probe). A total of 31 residues contact the cavity with their side chains. The side-chain cavity surface is partitioned according to the residue type as follows: 36-39% hydrophobic, 21-25% hydrophilic, and 37-43% neutral or ambivalent. Thus, the cavity surface is neither like a typical protein interior core, nor is like a typical protein external surface. All hydrophilic residues that contact the cavity-with the exception of Asp74-are clustered on the one side of the cavity. The cavity appears to expand its hydrophobic surface upon fatty acid binding on the side opposite to this hydrophilic patch. In holo I-FABP the fatty acid chain interactions with the hydrophilic side chains are mediated by water molecules. Molecular dynamics (MD) simulation of fully solvated apo I-FABP showed global conformational changes of I-FABP, which resulted in a large, but seemingly transient, exposure of the cavity to the external solvent. The packing density of the side chains lining the cavity, studied by Voronoi volumes, showed the presence of two distinctive small hydrophobic cores. The MD simulation predicts significant structural perturbations of the cavity on the subnanosecond time scale, which are capable of facilitating exchange of I-FABP internal water.     

Dramatic differences in the motions of the mouth of open and closed cytochrome P450BM-3 by molecular dynamics simulations

Molecular dynamics trajectories were calculated separately for each of the two molecules in the asymmetric unit of the crystal structure of the hemoprotein domain of cytochrome P450BM-3. Each simulation was 200 ps in length and included a 10 Å layer of explicit solvent. The simulated time-average structure of each P450BM-3 molecule is closer to its crystal structure than the two molecular dynamics time-averaged structures are to each other. In the crystal structure, molecule 2 has a more accessible substrate binding pocket than molecule 1, and this difference is maintained throughout the simulations presented here. In particular, the substrate docking regions of molecule 1 and molecule 2 diverge in the solution state simulations. The mouth of the substrate binding pocket is significantly more mobile in the simulation of molecule 2 than in the simulation of molecule 1. For molecule 1, the width of the mouth is only slightly larger than its X-ray value of 8.7 Å and undergoes fluctuations of about 1 Å. However, in molecule 2, the mouth of the substrate binding pocket is dramatically more open in the time-average molecular dynamics structure (14.7 Å) than in the X-ray structure (10.9 Å). Furthermore, this region of the protein undergoes large amplitude motions during the trajectory that are not seen in the trajectory of molecule 1, repeatedly opening and closing up to 7 Å. Presumably, the binding of different substrates will induce the mouth region to adopt different conformations from within the wide range of structures that are accessible. © 1995 Wiley-Liss, Inc.     

Binding of Net-Fla,a Netropsin-Flavin Hybrid Molecule,to DNA: Molecular Mechanics and Dynamics Studies in vacuo and in Water Solution

Abstract

We have studied the binding of the hybrid netropsin-flavin (Net-Fla) molecule onto four sequences containing four A.T base pairs. Molecular mechanics minimizations in vacuo show numerous minimal conformations separated by one base pair. 400 ps molecular dynamics simulations in vacuo have been performed using the lowest minima as the starting conformations. During these simulations, the flavin moiety of the drug makes two hydrogen bonds with an amino group of a neighboring guanine. A 200 ps molecular dynamics simulation in explicit water solution suggests that the binding of Net-Fla upon the DNA substrate is enhanced by water bridges. A water molecule bridging the amidinium of Net-Fla to the N3 atom of an adenine seems to be stuck in the dmg-DNA complex during the whole simulation. The fluctuations of the DNA helical parameters and of the torsion angles of the sugar-phosphate backbone are very similar in the simulations in vacuo and in water. The time auto-correlation functions for the DNA helical parameters decrease rapidly in the picosecond range in vacuo. The same functions computed from the water solution molecular dynamics simulations seem to have two modes: the rapid mode is similar to the behavior in vacuo, and is followed by a slower mode in the 10 ps range.     

Water penetration and escape in proteins

The kinetics of water penetration and escape in cytochrome c (cyt c) is studied by molecular dynamics (MD) simulations at various temperatures. Water molecules that penetrate the protein interior during the course of an MD simulation are identified by monitoring the number of water molecules in the first coordination shell (within 3.5 A) of each water molecule in the system. Water molecules in the interior of cyt c have 0-3 water molecules in their first hydration shell and this coordination number persists for extended periods of time. At T = 300 K we identify over 200 events in which water molecules penetrate the protein and reside inside for at least 5 picoseconds (ps) within a 1.5 nanoseconds (ns) time period. Twenty-seven (27) water molecules reside for at least 300 ps, 17 water molecules reside in the protein interior for times longer than 500 ps, and two interior water molecules do not escape; at T = 360 K one water molecule does not escape; at 430 K all water molecules exchange. Some of the internal water molecules show mean square displacements (MSD) of 1 A2 characteristic of structural waters. Others show MSD as large as 12 A2, suggesting that some of these water molecules occupy transient cavities and diffuse extensively within the protein. Motions of protein-bound water molecules are rotationally hindred, but show large librations. Analysis of the kinetics of water escape in terms of a survival time correlation function shows a power law behavior in time that can be interpreted in terms of a broad distribution of energy barriers, relative to kappa BT, for water exchange. At T = 300 K estimates of the roughness of the activation energy distribution is 4-10 kJ/mol (2-4 kappa BT). Activation enthalpies for water escape are 6-23 kJ/mol. The difference in activation entropies between fast exchanging (0.01 ns) and slow exchanging (0.1-1 ns) water molecules is -27 J/K/mol. Dunitz (Science 1997;264:670.) has estimated the maximum entropy loss of a water molecule due to binding to be 28 J/K/mol. Therefore, our results suggest that the entropy of interior water molecules is similar to entropy of bulk water.     

'Wave-type' structure of a synthetic hexaglycosylated decapeptide: A part of the extracellular domain of human glycophorin A

The three-dimensional structure of a glycopeptide, His-Thr*-Ser*-Thr*-Ser*-Ser*-Ser*-Val-Thr-Lys, with 2-acetamido-2-deoxy--D-galactose (GalNAc) residues linked to six adjacent amino acids from Thr-10 to Ser-15, was studied by NMR spectroscopy and molecular dynamics (MD) simulations. The hexaglycosylated decapeptide is part of the extracellular domain of human glycophorin A and shows an extended structure of the peptide backbone due to O-glycosylation. Furthermore, each GalNAc residue exhibits one and only one NOE contact from the NHAc proton to the backbone amide proton of the amino acid that the sugar is directly bound to. This indicates a strong preference for the orientation of all GalNAc residues towards the N-terminus. NOE build-up curves were used to determine 42 inter-proton distances that, in connection with angles of the peptide backbone obtained from 3J-coupling constants, resulted in constraints for a MD simulation in water. The NMR data and the MD simulations show a preference for an extended backbone structure. The GalNAc residues are located alternatingly on opposite sides of the backbone and reduce the flexibility of the peptide backbone. The conformation of the molecule is relatively rigid and shows a 'wave-type' 3D structure of the peptide backbone within the glycosylation cluster. This new structural element is also supported by the unusual CD spectrum of the glycopeptide.     

Solution structure of human intestinal fatty acid binding protein: Implications for ligand entry and exit

The human intestinal fatty acid binding protein (I-FABP) is a small (131 amino acids) proteinwhich binds dietary long-chain fatty acids in the cytosol of enterocytes. Recently, an alanineto threonine substitution at position 54 in I-FABP has been identified which affects fatty acidbinding and transport, and is associated with the development of insulin resistance in severalpopulations including Mexican-Americans and Pima Indians. To investigate the molecularbasis of the binding properties of I-FABP, the 3D solution structure of the more commonform of human I-FABP (Ala54) was studied by multidimensional NMR spectroscopy.Recombinant I-FABP was expressed from E. coli in the presence and absence of 15N-enriched media. The sequential assignments for non-delipidated I-FABP were completed byusing 2D homonuclear spectra (COSY, TOCSY and NOESY) and 3D heteronuclear spectra(NOESY-HMQC and TOCSY-HMQC). The tertiary structure of human I-FABP wascalculated by using the distance geometry program DIANA based on 2519 distance constraintsobtained from the NMR data. Subsequent energy minimization was carried out by using theprogram SYBYL in the presence of distance constraints. The conformation of human I-FABPconsists of 10 antiparallel -strands which form two nearly orthogonal -sheets offive strands each, and two short -helices that connect the -strands A and B. Theinterior of the protein consists of a water-filled cavity between the two -sheets. TheNMR solution structure of human I-FABP is similar to the crystal structure of rat I-FABP.The NMR results show significant conformational variability of certain backbone segmentsaround the postulated portal region for the entry and exit of fatty acid ligand.     

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.     

Insights into the function of the zinc hydroxide-Thr199-Glu106 hydrogen bonding network in carbonic anhydrases

The exact functional role of the zinc hydroxide (water)-Thr199-Glu106 hydrogen bond network in the carbonic anhydrases is unknown. However, from the results of molecular dynamics simulations (MD) we are able to better define its function. From computer graphics analysis and MD simulations on the zinc hydroxide form of human carbonic anhydrase II we find that this interaction forces the hydroxide hydrogen atom to be in a "down" position relative to the deep water-binding pocket. From previous work we have found that this pocket is a high-affinity binding site for CO2. We also note that during the timescale of our simulation (126 ps) the hydrogen bonds between the hydroxide hydrogen atom and Thr199 and the one between Thr199 and Glu106 are not fluxional. We propose that the role of the zinc hydroxide (water)-Thr199-Glu106 hydrogen bond network is to lock the hydrogen atom in the down position in order to expose the CO2 molecule bound in the deep water pocket to a lone pair of the hydroxide oxygen atom. This would allow for the rapid reaction of the CO2 molecule around the zinc ion. Furthermore, if the hydroxide hydrogen atom were not locked in the down position the binding of CO2 to the deep water pocket could be interfered with by the unrestrained hydroxide hydrogen atom (e.g. the N-Zn-O-H torsion could undergo rotational transitions that would partially block the deep water pocket). In summary, the roles we ascribe to this hydrogen bonding network are (1) to allow for facile access of CO2 to the deep water pocket and (2) to allow for maximal exposure of a hydroxide oxygen lone pair to the CO2 carbon atom.     

Exploration of the conformational space of myosin recovery stroke via molecular dynamics

Muscle contractions are driven by cyclic conformational changes of myosin, whose molecular mechanisms of operation are being elucidated by recent advances in crystallographic studies and single molecule experiments. To complement such structural studies and consider the energetics of the conformational changes of myosin head, umbrella sampling molecular dynamics (MD) simulations were performed with the all-atom model of the scallop myosin sub-fragment 1 (S1) with a bound ATP in solution in explicit water using the crystallographic near-rigor and transition state conformations as two references. The constraints on RMSD reaction coordinates used for the umbrella sampling were found to steer the conformational changes efficiently, and relatively close correlations have been observed between the set of characteristic structural changes including the lever arm rotation and the closing of the nucleotide binding pocket. The lever arm angle and key residue interaction distances in the nucleotide binding pocket and the relay helix show gradual changes along the recovery stroke reaction coordinate, consistent with previous crystallographic and computational minimum energy studies. Thermal fluctuations, however, appear to make the switch-2 coordination of ATP more flexible than suggested by crystal structures. The local solvation environment of the fluorescence probe, Trp 507 (scallop numbering), also appears highly mobile in the presence of thermal fluctuations.     

Receptor rigidity and ligand mobility in trypsin-ligand complexes

The trypsin-like serine proteases comprise a structurally similar family of proteins with a wide diversity of biological functions. Members of this family play roles in digestion, hemostasis, immune responses, and cancer metastasis. Bovine trypsin is an archetypical member of this family that has been extensively characterized both functionally and structurally, and that preferentially hydrolyzes Arg/Lys-Xaa peptide bonds. We have used molecular dynamics (MD) simulations to study bovine trypsin complexed with the two noncovalent small-molecule ligands, benzamidine and tranylcypromine, that have the same hydrogen-bond donating moieties as Arg and Lys side-chains, respectively. Multiple (10) simulations ranging from 1 ns to 2.2 ns, with explicit water molecules and periodic boundary conditions, were performed. The simulations reveal that the trypsin binding pocket residues are relatively rigid regardless of whether there is no ligand, a high-affinity ligand (benzamidine), or a low-affinity ligand (tranylcypromine). The thermal average of the conformations sampled by benzamidine bound to trypsin is planar and consistent with the planar internal geometry of the benzamidine crystallographic model coordinates. However, the most probable bound benzamidine conformations are +/-25 degrees out of plane, implying that the observed X-ray electron density represents an average of densities from two mirror symmetric, nonplanar conformations. Solvated benzamidine has free energy minima at +/-45 degrees , and the induction of a more planar geometry upon binding is associated with approximately 1 kcal/mol of intramolecular strain. Tranylcypromine's hydrogen-bonding pattern in the MD differs substantially from that inferred from the X-ray electron density. Early in simulations of this system, tranylcypromine adopts an alternative binding conformation, changing from the crystallographic conformation, with a direct hydrogen bond between its amino moiety and the backbone oxygen of Gly219, to one having a bridging water molecule. This result is consistently seen with the CHARMM22, Amber, or OPLS-AA force fields. The trypsin-tranylcypromine hydrogen-bonding pattern observed in the simulations also occurs as the crystallographic binding mode of the Lys15 side-chain of bovine pancreatic trypsin inhibitor bound to trypsin. In this latter cocrystal, a bridging crystallographic water does reside between the side-chain's amino group and the trypsin Gly219 backbone oxygen. Furthermore, the trypsin-tranylcypromine simulations sample two different stable noncrystallographic binding poses. These data suggest that some of the electron density ascribed to tranylcypromine in the X-ray model is rather due to a bound water molecule, and that multiple tranylcypromine binding conformations (crystallographic disorder) may be the cause of ambiguous electron density. The combined trypsin-benzamidine and trypsin- tranylcypromine results highlight the ability of simulations to augment protein-ligand complex structural data by deconvoluting the effects of thermal and structural averaging, and by finding energetically optimal ligand and bound water positions for weakly bound ligands.     

The interaction of lipophilic drugs with intestinal fatty acid-binding protein

Intestinal fatty acid-binding protein (I-FABP) is a small protein that binds long-chain dietary fatty acids in the cytosol of the columnar absorptive epithelial cells (enterocytes) of the intestine. The binding cavity of I-FABP is much larger than is necessary to bind a fatty acid molecule, which suggests that the protein may be able to bind other hydrophobic and amphipathic ligands such as lipophilic drugs. Herein we describe the binding of three structurally diverse lipophilic drugs, bezafibrate, ibuprofen (both R- and S-isomers) and nitrazepam to I-FABP. The rank order of affinity for I-FABP determined for these compounds was found to be R-ibuprofen approximately bezafibrate > S-ibuprofen > nitrazepam. The binding affinities were not directly related to aqueous solubility or partition coefficient of the compounds; however, the freely water-soluble drug diltiazem showed no affinity for I-FABP. Drug-I-FABP interaction interfaces were defined by analysis of chemical shift perturbations in NMR spectra, which revealed that the drugs bound within the central fatty acid binding cavity. Each drug participated in a different set of interactions within the cavity; however, a number of common contacts were observed with residues also involved in fatty acid binding. These data suggest that the binding of non-fatty acid lipophilic drugs to I-FABP may increase the cytosolic solubility of these compounds and thereby facilitate drug transport from the intestinal lumen across the enterocyte to sites of distribution and metabolism.     

Determination of the Backbone Mobility of Ribonuclease T1 and its 2′GMP Complex Using Molecular Dynamics Simulations and NMR Relaxation Data

Abstract

The results of 1-nanosecond molecular dynamics simulations of the enzyme ribonuclease T1 and its 2′GMP complex in water are presented. A classification of the angular reorientations of the backbone amide groups is achieved via a transformation of NH-vector trajectories into several coordinate frames, thus unravelling contributions of NH-bond librations and backbone dihedral angle fluctuations.

The former turned out to be similar for all amides, as characterized by correlation times of librational motions in a subpicosecond scale, angular amplitudes of about 10–12° for out-of-peptide-plane displacements of the NH-bond and 3–5° for the in-plane displacements, whereas the contributions of much slower backbone dihedral angle fluctuations strongly depend on the secondary structure. Correlation functions relevant for NMR were obtained and analyzed utilizing the ‘model-free’ approach (Lipari, G. and Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546–4559,4559-4570; Clore et al., (1990) J. Am. Chem. Soc. 112, 4989–4991). The dependence of the amplitude of local motion on the residue location in the backbone is in good agreement with the results of NMR relaxation measurements and X-ray data. The protein dynamics is characterized by a highly restricted local motion of those parts of the backbone with defined secondary structure as well as by a high flexibility in loop regions. The comparison of results derived from different periods of the trajectory (of 50 ps and 1 ns duration, 1000 points sampled) reveals a dependence of the observed dynamic picture on the characteristic time scale of the experimental method used. Comparison of the MD data for the free and liganded enzyme clearly indicates a restriction of the mobility within certain regions of the backbone upon inhibitor binding.     

Polyene fatty acid interactions with recombinant intestinal and liver fatty acid-binding proteins. Spectroscopic studies

Binding and proximity relationships of fatty acids with recombinant rat liver fatty acid-binding protein (L-FABP) and intestinal fatty acid-binding protein (I-FABP) were studied with absorption and fluorescence spectroscopy. Protein aromatic amino acids were examined in the absence and presence of bound fatty acid. Second derivative absorbance spectroscopy of the apo- and holoproteins suggested that fatty acid binding altered the conformation of L-FABP, but not of I-FABP. Fatty acid binding also blocked the accessibility of L-FABP tyrosine and I-FABP tryptophan to Stern-Volmer quenching by acrylamide, indicating that these amino acids were present in the fatty acid-binding pocket. Forster energy transfer from I-FABP tryptophan to bound cis-parinaric acid resulted in quenching of tryptophan lifetime and appearance of sensitized lifetime of bound cis-parinaric acid. The calculated donor-acceptor distances were 16.9 +/- 0.6 and 19.2 +/- 0.3 A for I-FABP and L-FABP, respectively. Absorbance spectral shifts and ratios of fluorescence excitation maxima indicated that the parinaric acid microenvironment in the fatty acid-binding site of I-FABP was much less polar than that of L-FABP. Parinaric acids displayed similar rotational correlation time and limiting anisotropy when bound to I-FABP and to L-FABP. These results are consistent with a close proximity of bound fatty acids to the tyrosine and tryptophan residues and with immobilization of the polyene fatty acids in the fatty acid-binding site(s) of L-FABP and I-FABP. The two proteins differ in that only L-FABP has two fatty acid-binding sites and appears to undergo significant conformational change upon fatty acid binding.     

Starting‐structure dependence of nanosecond timescale intersubstate transitions and reproducibility of MD‐derived order parameters

Factors affecting the accuracy of molecular dynamics (MD) simulations are investigated by comparing generalized order parameters for backbone NH vectors of the B3 immunoglobulin‐binding domain of streptococcal protein G (GB3) derived from simulations with values obtained from NMR spin relaxation (Yao L, Grishaev A, Cornilescu G, Bax A, J Am Chem Soc 2010;132:4295‐4309.). Choices for many parameters of the simulations, such as buffer volume, water model, or salt concentration, have only minor influences on the resulting order parameters. In contrast, seemingly minor conformational differences in starting structures, such as orientations of sidechain hydroxyl groups, resulting from applying different protonation algorithms to the same structure, have major effects on backbone dynamics. Some, but not all, of these effects are mitigated by increased sampling in simulations. Most discrepancies between simulated and experimental results occur for residues located at the ends of secondary structures and involve large amplitude nanosecond timescale transitions between distinct conformational substates. These transitions result in autocorrelation functions for bond vector reorientation that do not converge when calculated over individual simulation blocks, typically of length similar to the overall rotational diffusion time. A test for convergence before averaging the order parameters from different blocks results in better agreement between order parameters calculated from different sets of simulations and with NMR‐derived order parameters. Thus, MD‐derived order parameters are more strongly affected by transitions between conformational substates than by fluctuations within individual substates themselves, while conformational differences in the starting structures affect the frequency and scale of such transitions. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Structure, dynamics and hydration of the nogalamycin-d(ATGCAT)2Complex determined by NMR and molecular dynamics simulations in solution.

The structure of the 1:1 nogalamycin:d(ATGCAT)2 complex has been determined in solution from high-resolution NMR data and restrained molecular dynamics (rMD) simulations using an explicit solvation model. The antibiotic intercalates at the 5'-TpG step with the nogalose lying along the minor groove towards the centre of the duplex. Many drug-DNA nuclear Overhauser enhancements (NOEs) in the minor groove are indicative of hydrophobic interactions over the TGCA sequence. Steric occlusion prevents a second nogalamycin molecule from binding at the symmetry-related 5'-CpA site, leading to the conclusion that the observed binding orientation in this complex is the preferred orientation free of the complication of end-effects (drug molecules occupy terminal intercalation sites in all X-ray structures) or steric interactions between drug molecules (other NMR structures have two drug molecules bound in close proximity), as previously suggested. Fluctuations in key structural parameters such as rise, helical twist, slide, shift, buckle and sugar pucker have been examined from an analysis of the final 500 ps of a 1 ns rMD simulation, and reveal that many sequence-dependent structural features previously identified by comparison of different X-ray structures lie within the range of dynamic fluctuations observed in the MD simulations. Water density calculations on MD simulation data reveal a time-averaged pattern of hydration in both the major and minor groove, in good agreement with the extensive hydration observed in two related X-ray structures in which nogalamycin is bound at terminal 5'-TpG sites. However, the pattern of hydration determined from the sign and magnitude of NOE and ROE cross-peaks to water identified in 2D NOESY and ROESY experiments identifies only a few "bound" water molecules with long residence times. These solvate the charged bicycloaminoglucose sugar ring, suggesting an important role for water molecules in mediating drug-DNA electrostatic interactions within the major groove. The high density of water molecules found in the minor groove in X-ray structures and MD simulations is found to be associated with only weakly bound solvent in solution.     

Dynamics of Ca2+-saturated calmodulin D129N mutant studied by multiple molecular dynamics simulations

Fifteen independent 1-nsec MD simulations of fully solvated Ca(2+) saturated calmodulin (CaM) mutant D129N were performed from different initial conditions to provide a sufficient statistical basis to gauge the significance of observed dynamical properties. In all MD simulations the four Ca(2+) ions remained in their binding sites, and retained a single water ligand as observed in the crystal structure. The coordination of Ca(2+) ions in EF-hands I, II, and III was sevenfold. In EF-hand IV, which was perturbed by the mutation of a highly conserved Asp129, an anomalous eightfold Ca(2+) coordination was observed. The Ca(2+) binding loop in EF-hand II was observed to dynamically sample conformations related to the Ca(2+)-free form. Repeated MD simulations implicate two well-defined conformations of Ca(2+) binding loop II, whereas similar effect was not observed for loops I, III, and IV. In 8 out of 15 MD simulations Ca(2+) binding loop II adopted an alternative conformation in which the Thr62 >C=O group was displaced from the Ca(2+) coordination by a water molecule, resulting in the Ca(2+) ion ligated by two water molecules. The alternative conformation of the Ca(2+) binding loop II appears related to the "closed" state involved in conformational exchange previously detected by NMR in the N-terminal domain fragment of CaM and the C-terminal domain fragment of the mutant E140Q. MD simulations suggest that conformations involved in microsecond exchange exist partially preformed on the nanosecond time scale.     

Locally accessible conformations of proteins: multiple molecular dynamics simulations of crambin.

Multiple molecular dynamics (MD) simulations of crambin with different initial atomic velocities are used to sample conformations in the vicinity of the native structure. Individual trajectories of length up to 5 ns sample only a fraction of the conformational distribution generated by ten independent 120 ps trajectories at 300 K. The backbone atom conformational space distribution is analyzed using principal components analysis (PCA). Four different major conformational regions are found. In general, a trajectory samples only one region and few transitions between the regions are observed. Consequently, the averages of structural and dynamic properties over the ten trajectories differ significantly from those obtained from individual trajectories. The nature of the conformational sampling has important consequences for the utilization of MD simulations for a wide range of problems, such as comparisons with X-ray or NMR data. The overall average structure is significantly closer to the X-ray structure than any of the individual trajectory average structures. The high frequency (less than 10 ps) atomic fluctuations from the ten trajectories tend to be similar, but the lower frequency (100 ps) motions are different. To improve conformational sampling in molecular dynamics simulations of proteins, as in nucleic acids, multiple trajectories with different initial conditions should be used rather than a single long trajectory.     

Accuracy and precision of NMR relaxation experiments and MD simulations for characterizing protein dynamics

Model-free parameters obtained from nuclear magnetic resonance (NMR) relaxation experiments and molecular dynamics (MD) simulations commonly are used to describe the intramolecular dynamical properties of proteins. To assess the relative accuracy and precision of experimental and simulated model-free parameters, three independent data sets derived from backbone ¹⁵N NMR relaxation experiments and two independent data sets derived from MD simulations of Escherichia coli ribonuclease HI are compared. The widths of the distributions of the differences between the order parameters for pairs of NMR data sets are congruent with the uncertainties derived from statistical analyses of individual data sets; thus, current protocols for analyzing NMR data encapsulate random uncertainties appropriately. Large differences in order parameters for certain residues are attributed to systematic differences between samples for intralaboratory comparisons and unknown, possibly magnetic field-dependent, experimental effects for interlaboratory comparisons. The widths of distributions of the differences between the order parameters for two NMR sets are similar to widths of distributions for an NMR and an MD set or for two MD sets. The linear correlations between the order parameters for an MD set and an NMR set are within the range of correlations observed between pairs of NMR sets. These comparisons suggest that the NMR and MD generalized order parameters for the backbone amide N—H bond vectors are of comparable accuracy for residues exhibiting motions on a fast time scale (<100 ps). Large discrepancies between NMR and MD order parameters for certain residues are attributed to the occurrence of “rare” motional events over the simulation trajectories, the disruption of an element of secondary structure in one of the simulations, and lack of consensus among the experimental data sets. Consequently, (easily detectable) severe distortions of local protein structure and infrequent motional events in MD simulations appear to be the most serious artifacts affecting the accuracy and precision, respectively, of MD order parameters relative to NMR values. In addition, MD order parameters for motions on a fast (<100 ps) timescale are more precisely determined than their NMR counterparts, thereby permitting more detailed dynamic characterization of biologically important residues by MD simulation than is sometimes possible by experimental methods. Proteins 28:481–493, 1997. © 1997 Wiley-Liss, Inc.