Low-mass molecular dynamics simulation: A simple and generic technique to enhance configurational sampling

CLN025 is one of the smallest fast-folding proteins. Until now it has not been reported that CLN025 can autonomously fold to its native conformation in a classical, all-atom, and isothermal–isobaric molecular dynamics (MD) simulation. This article reports the autonomous and repeated folding of CLN025 from a fully extended backbone conformation to its native conformation in explicit solvent in multiple 500-ns MD simulations at 277 K and 1 atm with the first folding event occurring as early as 66.1 ns. These simulations were accomplished by using AMBER forcefield derivatives with atomic masses reduced by 10-fold on Apple Mac Pros. By contrast, no folding event was observed when the simulations were repeated using the original AMBER forcefields of FF12SB and FF14SB. The results demonstrate that low-mass MD simulation is a simple and generic technique to enhance configurational sampling. This technique may propel autonomous folding of a wide range of miniature proteins in classical, all-atom, and isothermal–isobaric MD simulations performed on commodity computers—an important step forward in quantitative biology.     

Self‐assembling study of sarcolipin and its mutants in multiple molecular dynamic simulations

The Sarcolipin (SLN) is a single trans‐membrane protein that can self‐assembly to dimer and oligomer for playing importantphysiological function. In this work, we addressed the dimerization of wild type SLN (wSLN) and its mutants (mSLNs) – I17A and I20A, using both coarse‐grained (CG) and atomistic (AT) molecular dynamics (MD) simulations. Our results demonstrated that wSLN homodimer assembled as a left‐handed helical complex, while mSLNs heterodimers assembled as right‐handed complexes. Analysis of residue‐residue contacts map indicated that isoleucine (Ile)‐leucione (Leu) zipper domain played an important role in dimerization. The potential of mean force (PMF) demonstrated that wSLN homodimer was more stable than mSLNs heterodimers. Meanwhile, the mSLNs heterodimers preferred right‐handed rather than left‐handed helix. AT‐MD simulations for wSLN and mSLNs were also in line with CG‐MD simulations. These results provided the insights for understanding the mechanisms of SLNs self‐assembling. Proteins 2017; 85:1065–1077. © 2017 Wiley Periodicals, Inc.     

Solvated and generalised Born calculations differences using GPU CUDA and multi-CPU simulations of an antifreeze protein with AMBER

While there has been an increase in the number of biomolecular computational studies employing graphics processing units (GPU), results describing their use with the molecular dynamics package AMBER with the CUDA implementation are scarce. No information is available comparing MD methodologies pmemd.cuda, pmemd.mpi or sander.mpi, available in AMBER, for generalised Born (GB) simulations or with solvated systems. As part of our current studies with antifreeze proteins (AFP), and for the previous reasons, we present details of our experience comparing performance of MD simulations at varied temperatures between multi-CPU runs using sander.mpi, pmemd.mpi and pmemd.cuda with the AFP from the fish ocean pout (1KDF). We found extremely small differences in total energies between multi-CPU and GPU CUDA implementations of AMBER12 in 1ns production simulations of the solvated system using the TIP3P water model. Additionally, GPU computations achieved typical one order of magnitude speedups when using mixed precision but were similar to CPU speeds when computing with double precision. However, we found that GB calculations were highly sensitive to the choice of initial GB parametrisation regardless of the type of methodology, with substantial differences in total energies.     

Temperature and pressure denaturation of chignolin: Folding and unfolding simulation by multibaric-multithermal molecular dynamics method

A multibaric‐multithermal molecular dynamics (MD) simulation of a 10‐residue protein, chignolin, was performed. All‐atom model with the Amber parm99SB force field was used for the protein and the TIP3P model was used for the explicit water molecules. This MD simulation covered wide ranges of temperature between 260 and 560 K and pressure between 0.1 and 600 MPa and sampled many conformations without getting trapped in local‐minimum free‐energy states. Folding events to the native β‐hairpin structure occurred five times and unfolding events were observed four times. As the temperature and/or pressure increases, fraction of folded chignolin decreases. The partial molar enthalpy change ΔH and partial molar volume change ΔV of unfolding were calculated as ΔH = 24.1 ± 4.9 kJ/mol and ΔV = ?5.6 ± 1.5 cm³/mol, respectively. These values agree well with recent experimental results. Illustrating typical local‐minimum free‐energy conformations, folding and unfolding pathways were revealed. When chignolin unfolds from the β‐hairpin structure, only the C terminus or both C and N termini open first. It may undergo an α‐helix or 3₁₀‐helix structure and finally unfolds to the extended structure. Difference of the mechanism between temperature denaturation and pressure denaturation is also discussed. Temperature denaturation is caused by making the protein transferred to a higher entropy state and making it move around more with larger space. The reason for pressure denaturation is that water molecules approach the hydrophobic residues, which are not well hydrated at the folded state, and some hydrophobic contacts are broken. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Global fold of human cannabinoid type 2 receptor probed by solid‐state 13C‐, 15N‐MAS NMR and molecular dynamics simulations

The global fold of human cannabinoid type 2 (CB₂) receptor in the agonist‐bound active state in lipid bilayers was investigated by solid‐state ¹³C‐ and ¹⁵N magic‐angle spinning (MAS) NMR, in combination with chemical‐shift prediction from a structural model of the receptor obtained by microsecond‐long molecular dynamics (MD) simulations. Uniformly ¹³C‐ and ¹⁵N‐labeled CB₂ receptor was expressed in milligram quantities by bacterial fermentation, purified, and functionally reconstituted into liposomes. ¹³C MAS NMR spectra were recorded without sensitivity enhancement for direct comparison of C_α, C_β, and C?O bands of superimposed resonances with predictions from protein structures generated by MD. The experimental NMR spectra matched the calculated spectra reasonably well indicating agreement of the global fold of the protein between experiment and simulations. In particular, the ¹³C chemical shift distribution of C_α resonances was shown to be very sensitive to both the primary amino acid sequence and the secondary structure of CB₂. Thus the shape of the C_α band can be used as an indicator of CB₂ global fold. The prediction from MD simulations indicated that upon receptor activation a rather limited number of amino acid residues, mainly located in the extracellular Loop 2 and the second half of intracellular Loop 3, change their chemical shifts significantly (≥1.5 ppm for carbons and ≥5.0 ppm for nitrogens). Simulated two‐dimensional ¹³C_α(i)? ¹³C?O(i) and ¹³C?O(i)? ¹⁵NH(i + 1) dipolar‐interaction correlation spectra provide guidance for selective amino acid labeling and signal assignment schemes to study the molecular mechanism of activation of CB₂ by solid‐state MAS NMR. Proteins 2014; 82:452–465. © 2013 Wiley Periodicals, Inc.     

Microsecond simulations of mdm2 and its complex with p53 yield insight into force field accuracy and conformational dynamics

The p53‐MDM2 complex is both a major target for cancer drug development and a valuable model system for computational predictions of protein‐ligand binding. To investigate the accuracy of molecular simulations of MDM2 and its complex with p53, we performed a number of long (200 ns to 1 µs) explicit‐solvent simulations using a range of force fields. We systematically compared nine popular force fields (AMBER ff03, ff12sb, ff14sb, ff99sb, ff99sb‐ildn, ff99sb‐ildn‐nmr, ff99sb‐ildn‐phi, CHARMM22*, and CHARMM36) against experimental chemical shift data, and found similarly accurate results, with microsecond simulations achieving better agreement compared to 200‐ns trajectories. Although the experimentally determined apo structure has a closed binding cleft, simulations in all force fields suggest the apo state of MDM2 is highly flexible, and able to sample holo‐like conformations, consistent with a conformational selection model. Initial structuring of the MDM2 lid region, known to competitively bind the binding cleft, is also observed in long simulations. Taken together, these results show molecular simulations can accurately sample conformations relevant for ligand binding. We expect this study to inform future computational work on folding and binding of MDM2 ligands. Proteins 2015; 83:1665–1676. © 2015 Wiley Periodicals, Inc.     

The new program OPAL for molecular dynamics simulations and energy refinements of biological macromolecules

Summary A new program for molecular dynamics (MD) simulation and energy refinement of biological macromolecules, OPAL, is introduced. Combined with the supporting program TRAJEC for the analysis of MD trajectories, OPAL affords high efficiency and flexibility for work with diferent force fields, and offers a user-friendly interface and extensive trajectory analysis capabilities. Salient features are computational speeds of up to 1.5 GFlops on vector supercomputers such as the NEC SX-3, ellipsoidal boundaries to reduce the system size for studies in explicit solvents, and natural treatment of the hydrostatic pressure. Practical applications of OPAL are illustrated with MD simulations of pure water, energy minimization of the NMR structure of the mixed disulfide of a mutant E. coli glutaredoxin with glutathione in different solvent models, and MD simulations of a small protein, pheromone Er-2, using either instantaneous or time-averaged NMR restraints, or no restraints.Abbreviations D diffusion constant in cm²/s - Er-2 pheromone 2 from Euplotes raikovi - GFlop one billion floating point operations per second - Grx(C14S)-SG mixed disulfide between a mutant E. coli glutaredoxin, with Cys¹⁴ replaced by Ser, and glutathione - MD molecular dynamics - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - density in g/cm³     

Motional timescale predictions by molecular dynamics simulations: Case study using proline and hydroxyproline sidechain dynamics

We propose a new approach for force field optimizations which aims at reproducing dynamics characteristics using biomolecular MD simulations, in addition to improved prediction of motionally averaged structural properties available from experiment. As the source of experimental data for dynamics fittings, we use ¹³C NMR spin‐lattice relaxation times T₁ of backbone and sidechain carbons, which allow to determine correlation times of both overall molecular and intramolecular motions. For structural fittings, we use motionally averaged experimental values of NMR J couplings. The proline residue and its derivative 4‐hydroxyproline with relatively simple cyclic structure and sidechain dynamics were chosen for the assessment of the new approach in this work. Initially, grid search and simplexed MD simulations identified large number of parameter sets which fit equally well experimental J couplings. Using the Arrhenius‐type relationship between the force constant and the correlation time, the available MD data for a series of parameter sets were analyzed to predict the value of the force constant that best reproduces experimental timescale of the sidechain dynamics. Verification of the new force‐field (termed as AMBER99SB‐ILDNP) against NMR J couplings and correlation times showed consistent and significant improvements compared to the original force field in reproducing both structural and dynamics properties. The results suggest that matching experimental timescales of motions together with motionally averaged characteristics is the valid approach for force field parameter optimization. Such a comprehensive approach is not restricted to cyclic residues and can be extended to other amino acid residues, as well as to the backbone. Proteins 2014; 82:195–215. © 2013 Wiley Periodicals, Inc.     

Rapid communication capturing the destabilizing effect of dihydrouridine through molecular simulations

The structural effects of the commonly occurring modified nucleoside dihydrouridine (D) observed experimentally in model oligonucleotides include a strong destabilization of the C3′‐endo sugar conformation of D, the disruption of stacking interactions of neighboring residues with D and a possible destabilization of the C3′‐endo sugar pucker of the 5′‐neighboring nucleoside. Our simulations with a combination of a set of parameters for modified RNA residues with the recently developed AMBER FF99χ force field having reoptimized glycosidic torsion angle parameters for standard nucleosides was found to reproduce the destabilizing effect of dihydrouridine better than with the AMBER FF99 force field for nucleic acids for which the parameters for the modified residues were originally developed. © 2014 Wiley Periodicals, Inc. Biopolymers 101: 985–991, 2014.     

Relating nucleotide‐dependent conformational changes in free tubulin dimers to tubulin assembly

The complex dynamic behavior of microtubules (MTs) is believed to be primarily due to the αβ‐tubulin dimer architecture and its intrinsic GTPase activity. Hence, a detailed knowledge of the conformational variations of isolated α‐GTP‐β‐GTP‐ and α‐GTP‐β‐GDP‐tubulin dimers in solution and their implications to interdimer interactions and stability is directly relevant to understand the MT dynamics. An attempt has been made here by combining molecular dynamics (MD) simulations and protein–protein docking studies that unravels key structural features of tubulin dimer in different nucleotide states and correlates their association to tubulin assembly. Results from simulations suggest that tubulin dimers and oligomers attain curved conformations in both GTP and GDP states. Results also indicate that the tubulin C‐terminal domain and the nucleotide state are closely linked. Protein–protein docking in combination with MD simulations suggest that the GTP‐tubulin dimers engage in relatively stronger interdimer interactions even though the interdimer interfaces are bent in both GTP and GDP tubulin complexes, providing valuable insights on in vitro finding that GTP‐tubulin is a better assembly candidate than GDP‐tubulin during the MT nucleation and elongation processes. © 2012 Wiley Periodicals, Inc. Biopolymers 99: 282–291, 2013.     

Dynamic simulations of stimuli-responsive switching of azobenzene derivatives in self-assembled monolayers: reactive rotation potential and switching functions

Although the force field (FF)-based molecular dynamics (MD) simulation has been widely applied to rationalise the experimental observations and measurements in chemistry, physics, materials and life science for years, traditional FF suffers from the incapability for describing chemical reactions, which are crucial in many important transformation processes. In order to simulate the collective switching process in azobenzene-based self-assemble monolayers on Au(111) surface, reactive MD simulations with alternative FF were implemented. The classic torsion function has been modified to depict the diabatic potential energy curves for cis and trans isomers, respectively. A switching function is further introduced to connect two N = N rotation functions, and the surface hopping between the cis and trans curves is allowed. By using the reactive rotation potential and switching function, the collective effect of numerous reaction centres and the influence of environment on the quantum yield in the complex system were explored at mesoscopic dimension and timescales. The reactive FF may be also applicable for other complicated systems containing stilbene derivatives. Limitation and perspective for further developments for the other complicated reactions are also addressed.     

Improved model of hydrated calcium ion for molecular dynamics simulations using classical biomolecular force fields

Calcium ions (Ca²⁺) play key roles in various fundamental biological processes such as cell signaling and brain function. Molecular dynamics (MD) simulations have been used to study such interactions, however, the accuracy of the Ca²⁺ models provided by the standard MD force fields has not been rigorously tested. Here, we assess the performance of the Ca²⁺ models from the most popular classical force fields AMBER and CHARMM by computing the osmotic pressure of model compounds and the free energy of DNA–DNA interactions. In the simulations performed using the two standard models, Ca²⁺ ions are seen to form artificial clusters with chloride, acetate, and phosphate species; the osmotic pressure of CaAc₂ and CaCl₂ solutions is a small fraction of the experimental values for both force fields. Using the standard parameterization of Ca²⁺ ions in the simulations of Ca²⁺‐mediated DNA–DNA interactions leads to qualitatively wrong outcomes: both AMBER and CHARMM simulations suggest strong inter‐DNA attraction whereas, in experiment, DNA molecules repel one another. The artificial attraction of Ca²⁺ to DNA phosphate is strong enough to affect the direction of the electric field‐driven translocation of DNA through a solid‐state nanopore. To address these shortcomings of the standard Ca²⁺ model, we introduce a custom model of a hydrated Ca²⁺ ion and show that using our model brings the results of the above MD simulations in quantitative agreement with experiment. Our improved model of Ca²⁺ can be readily applied to MD simulations of various biomolecular systems, including nucleic acids, proteins and lipid bilayer membranes. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 752–763, 2016.     

Right‐ and left‐handed three‐helix proteins. II. Similarity and differences in mechanical unfolding of proteins

Here, we study mechanical properties of eight 3‐helix proteins (four right‐handed and four left‐handed ones), which are similar in size under stretching at a constant speed and at a constant force on the atomic level using molecular dynamics simulations. The analysis of 256 trajectories from molecular dynamics simulations with explicit water showed that the right‐handed three‐helix domains are more mechanically resistant than the left‐handed domains. Such results are observed at different extension velocities studied (192 trajectories obtained at the following conditions: v = 0.1, 0.05, and 0.01 Å ps^?1, T = 300 K) and under constant stretching force (64 trajectories, F = 800 pN, T = 300 K). We can explain this by the fact, at least in part, that the right‐handed domains have a larger number of contacts per residue and the radius of cross section than the left‐handed domains. Proteins 2014; 82:90–102. © 2013 Wiley Periodicals, Inc.     

Comparison of force fields for Alzheimer's A : A case study for intrinsically disordered proteins

Intrinsically disordered proteins are essential for biological processes such as cell signalling, but are also associated to devastating diseases including Alzheimer's disease, Parkinson's disease or type II diabetes. Because of their lack of a stable three‐dimensional structure, molecular dynamics simulations are often used to obtain atomistic details that cannot be observed experimentally. The applicability of molecular dynamics simulations depends on the accuracy of the force field chosen to represent the underlying free energy surface of the system. Here, we use replica exchange molecular dynamics simulations to test five modern force fields, OPLS, AMBER99SB, AMBER99SB*ILDN, AMBER99SBILDN‐NMR and CHARMM22*, in their ability to model Aβ₄₂, an intrinsically disordered peptide associated with Alzheimer's disease, and compare our results to nuclear magnetic resonance (NMR) experimental data. We observe that all force fields except AMBER99SBILDN‐NMR successfully reproduce local NMR observables, with CHARMM22* being slightly better than the other force fields.     

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.     

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.     

Temperature‐accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain

The catalytic domain of the adenyl cyclase (AC) toxin from Bordetella pertussis is activated by interaction with calmodulin (CaM), resulting in cAMP overproduction in the infected cell. In the X‐ray crystallographic structure of the complex between AC and the C terminal lobe of CaM, the toxin displays a markedly elongated shape. As for the structure of the isolated protein, experimental results support the hypothesis that more globular conformations are sampled, but information at atomic resolution is still lacking. Here, we use temperature‐accelerated molecular dynamics (TAMD) simulations to generate putative all‐atom models of globular conformations sampled by CaM‐free AC. As collective variables, we use centers of mass coordinates of groups of residues selected from the analysis of standard molecular dynamics (MD) simulations. Results show that TAMD allows extended conformational sampling and generates AC conformations that are more globular than in the complexed state. These structures are then refined via energy minimization and further unrestrained MD simulations to optimize inter‐domain packing interactions, thus resulting in the identification of a set of hydrogen bonds present in the globular conformations. Proteins 2014; 82:2483–2496. © 2014 Wiley Periodicals, Inc.     

Force field impact and spin-probe modeling in molecular dynamics simulations of spin-labeled T4 lysozyme

Several attempts have been made to compute electron paramagnetic resonance (EPR) spectra of biomolecules, using motional models or simulated trajectories to describe dynamics. Ideally, the simulated trajectories should capture fast (picosecond) snapshots of spin-probe rotations accurately, while being lengthy enough to ensure a proper Fourier integration of the time-domain signal. It is the interplay of the two criteria that poses computational challenges to the method. In this context, an analysis of the spin-probe and protein conformational sampling and equilibration, with different force fields and with explicit solvent, may be a useful attempt. The present work reports a comparative study of the effect of the molecular dynamics (MD) force field on conformational sampling and equilibration in two spin-labeled T4 lysozyme (T4L) variants, N40C and K48C. Ensembles of 10× 3 ns-trajectories per variant and per force field (OPLS/AMBER and AMBER99) are analyzed for a reliable assessment of convergence and sampling. It is found that subtle site-dependent differences in spin-probe rotations and torsions are more readily captured in the AMBER99 trajectories than in the OPLS/AMBER simulations. On the other hand, sampling and equilibration are found to be better with the OPLS/AMBER force field at equal trajectory lengths.Figure: Left panel: The spin-probe R1 ring and the spin-probe Euler angles , and . Middle panel: Illustration of the diffusion in a cone model for the spin-probe motion: snapshots of helix B and of the R1 ring in N40C, taken at 0.3 ps intervals from AMBER trajectory 1. Right panel: The N40C mutant with the spin label (solid mode), solvated in a cubic box.     

Dynamics in an unusual acyl carrier protein from a ladderane lipid-synthesizing organism

Anaerobic ammonium-oxidizing (anammox) bacteria express a distinct acyl carrier protein implicated in the biosynthesis of the highly unusual “ladderane” lipids these organisms produce. This “anammox-specific” ACP, or amxACP, shows several unique features such as a conserved FF motif and an unusual sequence in the functionally important helix III. Investigation of the protein's structure and dynamics, both in the crystal by ensemble refinement and by MD simulations, reveals that helix III adopts a rare six-residue-long 3₁₀-helical conformation that confers a large degree of conformational and positional variability on this part of the protein. This way of introducing structural flexibility by using the inherent properties of 3₁₀-helices appears unique among ACPs. Moreover, the structure suggests a role for the FF motif in shielding the thioester linkage between the protein's prosthetic group and its acyl cargo from hydrolysis.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.