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.     

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.     

Effect of hydrogen coverage on the buckling of penta-graphene by molecular dynamics simulation

By employing a series of MD simulations, buckling behaviour of penta-graphene and functionalised penta-graphene having different hydrogen (H) coverage is presented in this study. To this end, the buckling onset strain is determined for different systems. The results reveal that the new allotrope is slightly stiffer than graphene. Moreover, the effect of H adatoms in the range 0–56% on buckling behaviour is investigated. Finally it is shown that the H coverage has an influence on stability, and ripple-type destortion of the penta-graphene under compression.     

Essential spaces defined by NMR structure ensembles and molecular dynamics simulation show significant overlap

Large concerted motions of proteins which span its “essential space,” are an important component of protein dynamics. We investigate to what extent structure ensembles generated with standard structure calculation techniques such as simulated annealing can capture these motions by comparing them to long-time molecular dynamics (MD) trajectories. The motions are analyzed by principal component analysis and compared using inner products of eigenvectors of the respective covariance matrices. Two very different systems are studied, the β-spectrin PH domain and the single-stranded DNA binding protein (ssDBP) from the filamentous phage Pf3. A comparison of the ensembles from NMR and MD shows significant overlap of the essential spaces, which in the case of ssDBP is extraordinarily high. The influence of variations in the specifications of distance restraints is investigated. We also study the influence of the selection criterion for the final structure ensemble on the definition of mobility. The results suggest a modified criterion that improves conformational sampling in terms of amplitudes of correlated motion. Proteins 31:370–382, 1998. © 1998 Wiley-Liss, Inc.     

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.     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。     

A molecular dynamics simulation study of segment B1 of protein G

The immunoglobulin binding protein, segment B1 of protein G, has been studied experimentally as a paradigm for protein folding. This protein consists of 56 residues, includes both β sheet and α helix and contains neither disulfide bonds nor proline residues. We report an all-atom molecular dynamics study of the native manifold of the protein in explicit solvent. A 2-ns simulation starting from the nuclear magnetic resonance (NMR) structure and a 1-ns control simulation starting from the x-ray structure were performed. The difference between average structures calculated over the equilibrium portion of trajectories is smaller than the difference between their starting conformations. These simulation averages are structurally similar to the x-ray structure and differ in systematic ways from the NMR-determined structure. Partitioning of the fluctuations into fast (<20 ps) and slow (<20 ps) components indicates that the β sheet displays greater long-time mobility than does the α helix. Clore and Gronenborn [J. Mol. Biol. 223:853–856, 1992] detected two long-residence water molecules by NMR in a solution structure of segment B1 of protein G. Both molecules were found in the fully exposed regions and were proposed to be stabilized by bifurcated hydrogen bonds to the protein backbone. One of these long-residence water molecules, found near an exposed loop region, is identified in both of our simulations, and is seen to be involved in the formation of a stable water-mediated hydrogen bond bridge. The second water molecule, located near the middle of the α helix, is not seen with an exceptional residence time in either as a result of the conformation being closer to the x-ray structure in this region of the protein. Proteins 29:193–202, 1997. © 1997 Wiley-Liss, Inc.     

Accessibility of tobacco lipid transfer protein cavity revealed by 15N NMR relaxation studies and molecular dynamics simulations

Plant LTP1 are small helical proteins stabilized by four disulfide bridges and are characterized by the presence of an internal cavity, in which various hydrophobic ligands can be inserted. Recently, we have determined the solution structure of the recombinant tobacco LTP1_1. Unexpectedly, despite a global fold very similar to the structures already known for cereal seed LTP1, its binding properties are different: Tobacco LTP1_1 is able to bind only one monoacylated lipid, whereas cereal LTP1 can bind either one or two. The 3D structure of tobacco LTP1_1 revealed the presence of a hydrophobic cluster, not observed on cereal LTP1 structures, which may hinder one of the two entrances of the cavity defined for wheat LTP1. To better understand the mechanism of lipid entrance for tobacco LTP1_1 and to define the regions of the protein monitoring the accessibility of the cavity, we have complemented our structural data by the study of the internal dynamics of tobacco LTP1_1, using (15)N magnetic relaxation rate data and MD simulations at room and high temperatures. This work allowed us to define two regions of the protein experiencing the largest motions. These two regions delineate a portal that opens up during the simulation constituting a unique entrance of the hydrophobic cavity, in contrast with wheat LTP1 where two routes were detected. The hydrophobic interactions resulting from a few point mutations are strong enough to completely block the second portal so that the accessibility of the cavity is restricted to one entrance, explaining why this particular LTP1 binds only one lipid molecule.     

Refinement of homology-based protein structures by molecular dynamics simulation techniques

The use of classical molecular dynamics simulations, performed in explicit water, for the refinement of structural models of proteins generated ab initio or based on homology has been investigated. The study involved a test set of 15 proteins that were previously used by Baker and coworkers to assess the efficiency of the ROSETTA method for ab initio protein structure prediction. For each protein, four models generated using the ROSETTA procedure were simulated for periods of between 5 and 400 nsec in explicit solvent, under identical conditions. In addition, the experimentally determined structure and the experimentally derived structure in which the side chains of all residues had been deleted and then regenerated using the WHATIF program were simulated and used as controls. A significant improvement in the deviation of the model structures from the experimentally determined structures was observed in several cases. In addition, it was found that in certain cases in which the experimental structure deviated rapidly from the initial structure in the simulations, indicating internal strain, the structures were more stable after regenerating the side-chain positions. Overall, the results indicate that molecular dynamics simulations on a tens to hundreds of nanoseconds time scale are useful for the refinement of homology or ab initio models of small to medium-size proteins.     

Improving NMR protein structure quality by Rosetta refinement: a molecular replacement study

The structure of human protein HSPC034 has been determined by both solution nuclear magnetic resonance (NMR) spectroscopy and X-ray crystallography. Refinement of the NMR structure ensemble, using a Rosetta protocol in the absence of NMR restraints, resulted in significant improvements not only in structure quality, but also in molecular replacement (MR) performance with the raw X-ray diffraction data using MOLREP and Phaser. This method has recently been shown to be generally applicable with improved MR performance demonstrated for eight NMR structures refined using Rosetta (Qian et al., Nature 2007;450:259-264). Additionally, NMR structures of HSPC034 calculated by standard methods that include NMR restraints have improvements in the RMSD to the crystal structure and MR performance in the order DYANA, CYANA, XPLOR-NIH, and CNS with explicit water refinement (CNSw). Further Rosetta refinement of the CNSw structures, perhaps due to more thorough conformational sampling and/or a superior force field, was capable of finding alternative low energy protein conformations that were equally consistent with the NMR data according to the Recall, Precision, and F-measure (RPF) scores. On further examination, the additional MR-performance shortfall for NMR refined structures as compared with the X-ray structure were attributed, in part, to crystal-packing effects, real structural differences, and inferior hydrogen bonding in the NMR structures. A good correlation between a decrease in the number of buried unsatisfied hydrogen-bond donors and improved MR performance demonstrates the importance of hydrogen-bond terms in the force field for improving NMR structures. The superior hydrogen-bond network in Rosetta-refined structures demonstrates that correct identification of hydrogen bonds should be a critical goal of NMR structure refinement. Inclusion of nonbivalent hydrogen bonds identified from Rosetta structures as additional restraints in the structure calculation results in NMR structures with improved MR performance.     

Solution structure of a pentasaccharide representing the repeating unit of the O-antigen polysaccharide from Escherichia coli O142: NMR spectroscopy and molecular simulation studies

Conformational studies have been performed of a pentasaccharide derived from the O-polysaccharide from Escherichia coli O142. The polymer was selectively degraded by anhydrous hydrogen fluoride and reduced to yield an oligosaccharide model of its repeating unit, which in the branching region consists of four aminosugars. A comparison of (1)H and (13)C chemical shifts between the pentasaccharide and the polymer showed only minor differences, except where the cleavage had taken place, indicating that the oligomer is a good model of the repeating unit. Langevin dynamics and molecular dynamics simulations with explicit water molecules were carried out to sample the conformational space of the pentasaccharide. For the glycosidic linkages between the hexopyranoside residues, small but significant changes were observed between the simulation techniques. One-dimensional (1D) (1)H,(1)H double pulsed field gradient spin echo (DPFGSE) transverse rotating-frame Overhauser effect spectroscopy (T-ROESY) experiments were performed, and homonuclear cross-relaxation rates were obtained. Subsequently, a comparison of interproton distances from NMR experiment and the two simulation approaches showed that in all cases the use of explicit water in the simulations resulted in better agreement. Hydrogen-bond analysis of the trajectories from the molecular dynamics simulation revealed interresidue interactions to be important as a cluster of different hydrogen bonds and as a distinct highly populated hydrogen bond. NMR data are consistent with the presence of hydrogen bonding within the model of the repeating unit.     

In silico approaches towards understanding CALB using molecular dynamics simulation and docking

Candida antarctica lipase B (CALB), a serine protease, is involved in the hydrolysis of substrates at the aqueous lipid interface. There is a significant role played by the helices in serine proteases including acting as a flap covering the active site region. The α5 and α10 helices in the path to the active site of CALB appear to play an important role in the region. This study investigates these helices by mutational studies, docking and molecular dynamics simulations. The mutations were selected based on their proximity to the active site and their presence at the α10-helix in the path of the active site. Molecular dynamics studies reveal the flexibility, stability and hydrogen bonding ability of the α5 helix. The radius of gyration (R _g) clearly showed the compactness of the structure. Docking studies show the changes occurring at the protein's binding site before and after 15 ns of simulation. Results from the study demonstrate the importance of the two helices α5 and α10 in the stability of CALB.     

Ligand orientation of human neuroglobin obtained from solution NMR and molecular dynamics simulation as compared with X-ray crystallography

Neuroglobin, a new member of hemoprotein family, can reversibly bind oxygen and take part in many biological processes such as enzymatic reaction, signal transduction and the mitochondria function. Different from myoglobin and hemoglobin, it has a hexacoordinated heme environment, with histidyl imidazole of proximal His⁹⁶(F8) and distal His⁶⁴(E7) directly bound to the metal ion. In the present work, solution ¹H NMR spectroscopy was employed to investigate the electronic structure of heme center of wild-type met-human neuroglobin. The resonances of heme protons and key residues in the heme pocket were assigned. Two heme orientations resulting from a 180° rotation about the α-γ-meso axis with a population ratio about 2:1 were observed. Then the ¹H NMR chemical shifts of the ferriheme methyl groups were used to predict orientations of the axial ligand. The obtained axial ligand plane angle φ is consistent with that from the molecular dynamics simulation but not with those from the crystal data. Compared with mouse neuroglobin, the obtained average ligand orientation of human neuroglobin reflects the changeability of heme environment for the Ngb family.     

Release of ADP from the catalytic subunit of protein kinase A: a molecular dynamics simulation study

Substrate phosphorylation by cAMP-dependent-protein kinase A (protein kinase A, PKA) has been studied extensively. Phosphoryl transfer was found to be fast, whereas ADP release was found to be the slow, rate-limiting step. There is also evidence that ADP release may be preceded by a partially rate-limiting conformational change. However, the atomic details of the conformational change and the mode of ADP release are difficult to obtain experimentally. In this work, we studied ADP release from PKA by carrying out molecular dynamics simulations with different pulling forces applied to the ligand. The detailed ADP release pathway and the associated conformational changes were analyzed. The ADP release process was found to involve a swinging motion with the phosphate of ADP anchored to the Gly-rich loop, so that the more buried adenine base and ribose ring came out before the phosphate. In contrast to the common belief that a hinge-bending motion was responsible for the opening of the ligand-binding cleft, our simulations showed that the small lobe exhibited a large amplitude "rocking" motion when the ligand came out. The largest conformational change of the protein was observed at about the first quarter time point along the release pathway. Two prominent intermediate states were observed in the release process.     

Conformational sampling by NMR solution structures calculated with the program DIANA evaluated by comparison with long-time molecular dynamics calculations in explicit water

The NMR solution structure of bovine pancreatic trypsin inhibitor (BPTI) obtained by distance geometry calculations with the program DIANA is compared with groups of conformers generated by molecular dynamics (MD) simulations in explicit water at ambient temperature and pressure. The MD simulations started from a single conformer and were free or restrained either by the experimental NOE distance restraints or by time-averaged restraints; the groups of conformers were collected either in 10 ps intervals during 200 ps periods of simulation, or in 50 ps intervals during a 1 ns period of simulation. Overall, these comparisons show that the standard protein structure determination protocol with the program DIANA provides a picture of the protein structure that is in agreement with MD simulations using “realistic” potential functions over a nanosecond timescale. For well-constrained molecular regions there is a trend in the free MD simulation of duration 1 ns that the sampling of the conformation space is slightly increased relative to the DIANA calculations. In contrast, for surface-exposed side-chains that are less extensively constrained by the NMR data, the DIANA conformers tend to sample larger regions of conformational space than conformers selected from any of the MD trajectories. Additional insights into the behavior of surface side-chains come from comparison of the MD runs of 200 ps or 1 ns duration. In this time range the sampling of conformation space by the protein surface depends strongly on the length of the simulation, which indicates that significant side-chain transitions occur on the nanosecond timescale and that much longer simulations will be needed to obtain statistically significant data on side-chain dynamics.     

Role of the electrostatic energy between the photo-products and ionic groups on the photoinduced electron transfer rates from aromatic amino acids to the excited flavin in five single-point substitution isoforms of the charged amino acid residue-13 in the FMN-binding protein

Effect of the charge (negative, positive or neutral) of amino acid residue-13 on the photoinduced electron transfer (ET) from Trp32, Tyr35 and Trp106 to the excited isoalloxazine was evaluated in the flavin mononucleotide-binding protein from Desulfovibrio vulgaris isolate Miyazaki F (DvFBP). The protein structures of the wild type and the four isoforms where glutamic acid-13 is replaced with lysine (E13K), arginine (E13R), threonine (E13T) and glutamine (E13Q) in aqueous solution were obtained by molecular dynamics simulation. The distances between the amino acid residue-13 and isoalloxazine (Iso), and between the amino acid residue-13 and the ET donors were longer than 1 nm. The ET rates were evaluated with the Kakitani and Mataga model (KM theory) from their ultrafast fluorescence dynamics by means of a non-linear least squares method. Electrostatic (ES) energies between the photo-products and other ionic groups in the proteins markedly varied among ET donors and among the DvFBP isoforms, while the other physical quantities related to the ET rates, the solvent reorganisation and ES energies between the Iso anion and the donor cations did not vary much between the proteins and donors. A plot of the logarithmic ET rates versus either the total free energy gaps or the net ES energies between the photo-products and the other ionic groups both displayed a parabolic function and so the net ES energies are an important influential factor upon the ET rate, in addition to the donor–acceptor distance.     

A comparison of the electromechanical properties of structurally diverse proteins by molecular dynamics simulation

Proteins are subjected to electric fields both within the cell and during routine biochemical analysis. We have used atomistic molecular dynamics simulations to study conformational changes within three structurally diverse proteins subjected to high electric fields. At electric fields in excess of .5?V/nm, major structural changes were observed in all three proteins due to charge redistribution within the biomolecule. However, the electromechanical resilience was found to be highly dependent on the protein secondary structure, with α-helices showing a particularly high susceptibility to deformation by the applied electric field.     

On the ability of molecular dynamics force fields to recapitulate NMR derived protein side chain order parameters

Molecular dynamics (MD) simulations have become a central tool for investigating various biophysical questions with atomistic detail. While many different proxies are used to qualify MD force fields, most are based on largely structural parameters such as the root mean square deviation from experimental coordinates or nuclear magnetic resonance (NMR) chemical shifts and residual dipolar couplings. NMR derived Lipari–Szabo squared generalized order parameter (O²) values of amide N? H bond vectors of the polypeptide chain were also often employed for refinement and validation. However, with a few exceptions, side chain methyl symmetry axis order parameters have not been incorporated into experimental reference sets. Using a test set of five diverse proteins, the performance of several force fields implemented in the NAMDD simulation package was examined. It was found that simulations employing explicit water implemented using the TIP3 model generally performed significantly better than those using implicit water in reproducing experimental methyl symmetry axis O² values. Overall the CHARMM27 force field performs nominally better than two implementations of the Amber force field. It appeared that recent quantum mechanics modifications to side chain torsional angles of leucine and isoleucine in the Amber force field have significantly hindered proper motional modeling for these residues. There remained significant room for improvement as even the best correlations of experimental and simulated methyl group Lipari–Szabo generalized order parameters fall below an R² of 0.8.     

Exploring single-domain antibody thermostability by molecular dynamics simulation

Abstract

Single-domain antibodies also known as nanobodies are recombinant antigen-binding domains that correspond to the heavy-chain variable region of camelid antibodies. Previous experimental studies showed that the nanobodies have stable and active structures at high temperatures. In this study, the thermal stability and dynamics of nanobodies have been studied by employing molecular dynamics simulation at different temperatures. Variations in root mean square deviation, native contacts, and solvent-accessible surface area of the nanobodies during the simulation were calculated to analyze the effect of different temperatures on the overall conformation of the nanobody. Then, the thermostability mechanism of this protein was studied through calculation of dynamic cross-correlation matrix, principal component analyses, native contact analyses, and root mean square fluctuation. Our results manifest that the side chain conformation of some residues in the complementarity-determining region 3 (CDR3) and also the interaction between α-helix region of CDR3 and framework2 play a critical role to stabilize the protein at a high temperature.

Communicated by Ramaswamy H. Sarma     

Implicit solvent models for flexible protein-protein docking by molecular dynamics simulation

The suitability of three implicit solvent models for flexible protein-protein docking by procedures using molecular dynamics simulation is investigated. The three models are (i) the generalized Born (GB) model implemented in the program AMBER6.0; (ii) a distance-dependent dielectric (DDD) model; and (iii) a surface area-dependent model that we have parameterized and call the NPSA model. This is a distance-dependent dielectric model modified by neutralizing the ionizable side-chains and adding a surface area-dependent solvation term. These solvent models were first tested in molecular dynamics simulations at 300 K of the native structures of barnase, barstar, segment B1 of protein G, and three WW domains. These protein structures display a range of secondary structure contents and stabilities. Then, to investigate the performance of the implicit solvent models in protein docking, molecular dynamics simulations of barnase/barstar complexation, as well as PIN1 WW domain/peptide complexation, were conducted, starting from separated unbound structures. The simulations show that the NPSA model has significant advantages over the DDD and GB models in maintaining the native structures of the proteins and providing more accurate docked complexes.