The gaussian network model: precise prediction of residue fluctuations and application to binding problems

The single-parameter Gamma matrix of force constants proposed by the Gaussian Network Model (GNM) is iteratively modified to yield native state fluctuations that agree exactly with experimentally observed values. The resulting optimized Gamma matrix contains residue-specific force constants that may be used for an accurate analysis of ligand binding to single or multiple sites on proteins. Bovine Pancreatic Trypsin Inhibitor (BPTI) is used as an example. The calculated off-diagonal elements of the Gamma matrix, i.e., the optimized spring constants, obey a Lorentzian distribution. The mean value of the spring constants is approximately -0.1, a value much weaker than -1 of the GNM. Few of the spring constants are positive, indicating repulsion between residues. Residue pairs with large number of neighbors have spring constants around the mean, -0.1. Large negative spring constants are between highly correlated pairs of residues. The fluctuations of the distance between anticorrelated pairs of residues are subject to smaller spring constants. The importance of the number of neighbors of residue pairs in determining the elements of the Gamma matrix is pointed out. Allosteric effects of binding on a single or multiple residues of BPTI are illustrated and discussed. Comparison of the predictions of the present model with those of the standard GNM shows that the two models agree at lower modes, i.e., those relating to global motions, but they disagree at higher modes. In the higher modes, the present model points to the important contributions from specific residues whereas the standard GNM fails to do so.     

Tensorial elastic network model for protein dynamics: Integration of the anisotropic network model with bond‐bending and twist elasticities

We present a tensorial elastic network model (TNM) to describe the equilibrium fluctuations of proteins near their native fold structure. The model combines the anisotropic network model (ANM), bond bending elasticity, and backbone twist elasticity, and can predict both the isotropic fluctuations, similar to the Gaussian network model (GNM), and anisotropic fluctuations, similar to the ANM. TNM performs equally well for B‐factor predictions as GNM and predicts the anisotropy of B‐factors better than ANM. The model also outperforms the ANM in its predictability of the complete anisotropic displacement parameters. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

vGNM: a better model for understanding the dynamics of proteins in crystals

The dynamics of proteins are important for understanding their functions. In recent years, the simple coarse-grained Gaussian Network Model (GNM) has been fairly successful in interpreting crystallographic B-factors. However, the model clearly ignores the contribution of the rigid body motions and the effect of crystal packing. The model cannot explain the fact that the same protein may have significantly different B-factors under different crystal packing conditions. In this work, we propose a new GNM, called vGNM, which takes into account both the contribution of the rigid body motions and the effect of crystal packing, by allowing the amplitude of the internal modes to be variables. It hypothesizes that the effect of crystal packing should cause some modes to be amplified and others to become less important. In doing so, vGNM is able to resolve the apparent discrepancy in experimental B-factors among structures of the same protein but with different crystal packing conditions, which GNM cannot explain. With a small number of parameters, vGNM is able to reproduce experimental B-factors for a large set of proteins with significantly better correlations (having a mean value of 0.81 as compared to 0.59 by GNM). The results of applying vGNM also show that the rigid body motions account for nearly 60% of the total fluctuations, in good agreement with previous findings.     

Dynamics of proteins predicted by molecular dynamics simulations and analytical approaches: application to alpha-amylase inhibitor

The dynamics of alpha-amylase inhibitors has been investigated using molecular dynamics (MD) simulations and two analytical approaches, the Gaussian network model (GNM) and anisotropic network model (ANM). MD simulations use a full atomic approach with empirical force fields, while the analytical approaches are based on a coarse-grained single-site-per-residue model with a single-parameter harmonic potential between sufficiently close (r 相似文献   

Insights into equilibrium dynamics of proteins from comparison of NMR and X-ray data with computational predictions

For a representative set of 64 nonhomologous proteins, each containing a structure solved by NMR and X-ray crystallography, we analyzed the variations in atomic coordinates between NMR models, the temperature (B) factors measured by X-ray crystallography, and the fluctuation dynamics predicted by the Gaussian network model (GNM). The NMR and X-ray data exhibited a correlation of 0.49. The GNM results, on the other hand, yielded a correlation of 0.59 with X-ray data and a distinctively better correlation (0.75) with NMR data. The higher correlation between GNM and NMR data, compared to that between GNM and X-ray B factors, is shown to arise from the differences in the spectrum of modes accessible in solution and in the crystal environment. Mainly, large-amplitude motions sampled in solution are restricted, if not inaccessible, in the crystalline environment of X-rays. Combined GNM and NMR analysis emerges as a useful tool for assessing protein dynamics.     

Structural features that predict real-value fluctuations of globular proteins

It is crucial to consider dynamics for understanding the biological function of proteins. We used a large number of molecular dynamics (MD) trajectories of nonhomologous proteins as references and examined static structural features of proteins that are most relevant to fluctuations. We examined correlation of individual structural features with fluctuations and further investigated effective combinations of features for predicting the real value of residue fluctuations using the support vector regression (SVR). It was found that some structural features have higher correlation than crystallographic B‐factors with fluctuations observed in MD trajectories. Moreover, SVR that uses combinations of static structural features showed accurate prediction of fluctuations with an average Pearson's correlation coefficient of 0.669 and a root mean square error of 1.04 Å. This correlation coefficient is higher than the one observed in predictions by the Gaussian network model (GNM). An advantage of the developed method over the GNMs is that the former predicts the real value of fluctuation. The results help improve our understanding of relationships between protein structure and fluctuation. Furthermore, the developed method provides a convienient practial way to predict fluctuations of proteins using easily computed static structural features of proteins. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

A comparative analysis of the equilibrium dynamics of a designed protein inferred from NMR,X‐ray,and computations

A detailed analysis of high‐resolution structural data and computationally predicted dynamics was carried out for a designed sugar‐binding protein. The mean‐square deviations in the positions of residues derived from nuclear magnetic resonance (NMR) models and those inferred from X‐ray crystallographic B‐factors for two different crystal forms were compared with the predictions based on the Gaussian Network Model (GNM) and the results from molecular dynamics (MD) simulations. GNM systematically yielded a higher correlation than MD, with experimental data, suggesting that the lack of atomistic details in the coarse‐grained GNM is more than compensated for by the mathematically exact evaluation of fluctuations using the native contacts topology. Evidence is provided that particular loop motions are curtailed by intermolecular contacts in the crystal environment causing a discrepancy between theory and experiments. Interestingly, the information conveyed by X‐ray crystallography becomes more consistent with NMR models and computational predictions when ensembles of X‐ray models are considered. Less precise (broadly distributed) ensembles indeed appear to describe the accessible conformational space under native state conditions better than B‐factors. Our results highlight the importance of using multiple conformations obtained by alternative experimental methods, and analyzing results from both coarse‐grained models and atomic simulations, for accurate assessment of motions accessible to proteins under native state conditions. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

A large-scale comparison of computational models on the residue flexibility for NMR-derived proteins

As an alternative to X-ray crystallography, nuclear magnetic resonance (NMR) has also emerged as the method of choice for studying both protein structure and dynamics in solution. However, little work using computational models such as Gaussian network model (GNM) and machine learning approaches has focused on NMR-derived proteins to predict the residue flexibility, which is represented by the root mean square deviation (RMSD) with respect to the average structure. We provide a large-scale comparison of computational models, including GNM, parameter-free GNM and several linear regression models using local solvent exposures as inputs, based on a dataset of 1609 protein chains whose structures were resolved by NMR. The result again confirmed that the correlation of GNM outputs with raw RMSD values was better than that using B-factors of X-ray data. Nevertheless, it was also concluded that the parameter-free GNM and the solvent exposure based linear regression models performed worse than GNM when predicting RMSD, contrary to results using X-ray data. The discrepancy of residue flexibility prediction between NMR and X-ray data is likely attributable to a combination of their physical and methodological differences.     

Identification of kinetically hot residues in proteins.

A number of recent studies called attention to the presence of kinetically important residues underlying the formation and stabilization of folding nuclei in proteins, and to the possible existence of a correlation between conserved residues and those participating in the folding nuclei. Here, we use the Gaussian network model (GNM), which recently proved useful in describing the dynamic characteristics of proteins for identifying the kinetically hot residues in folded structures. These are the residues involved in the highest frequency fluctuations near the native state coordinates. Their high frequency is a manifestation of the steepness of the energy landscape near their native state positions. The theory is applied to a series of proteins whose kinetically important residues have been extensively explored: chymotrypsin inhibitor 2, cytochrome c, and related C2 proteins. Most of the residues previously pointed out to underlie the folding process of these proteins, and to be critically important for the stabilization of the tertiary fold, are correctly identified, indicating a correlation between the kinetic hot spots and the early forming structural elements in proteins. Additionally, a strong correlation between kinetically hot residues and loci of conserved residues is observed. Finally, residues that may be important for the stability of the tertiary structure of CheY are proposed.     

Coupling of structural fluctuations to deamidation reaction in triosephosphate isomerase by Gaussian network model

We study the structural fluctuations of triosephosphate isomerase (TIM) by an elastic model, namely, the Gaussian network model (GNM), to identify a network of coupled motions in the allosteric communication between its deamidation and catalytic sites, and the promoting motions for the deamidation activity. For this, three TIM structures have been studied: one crystal structure and two model structures designed to describe different putative models for the deamidation reaction taking place at the subunit interface. The structural fluctuations have been mapped on the functional properties; then the differences in the fluctuations between the two models in relation to the deamidation reaction have been considered. The results demonstrate that the qualitative picture of the mean-square fluctuations and the correlations between the fluctuations are similar in both, but the differences may affect the observed barrier height of the deamidation reaction. The higher packing density at regions close to deamidation sites, reflected by the high-frequency fluctuating residues in the respective regions, the stronger positive correlation between the fluctuations of the deamidation sites, and enhanced positive correlation of the primary deamidation site with the extended vicinity of the catalytic region on the juxtaposed unit promote the probability of the deamidation reaction. The results in general emphasize the importance of structural fluctuations in enzyme reactions, as well as proposing the present methodology as a plausible approach for studies on the network of coupled promoting motions in protein functions.     

Functional motions of influenza virus hemagglutinin: a structure-based analytical approach.

Influenza virus hemagglutinin (HA), a homotrimeric integral membrane glycoprotein essential for viral infection, is engaged in two biological functions: recognition of target cells' receptor proteins and fusion of viral and endosomal membranes, both requiring substantial conformational flexibility from the part of the glycoprotein. The different modes of collective motions underlying the functional mobility/adaptability of the protein are determined in the present study using an extension of the Gaussian network model (GNM) to treat concerted anisotropic motions. We determine the molecular mechanisms that may underlie HA function, along with the structural regions or residues whose mutations are expected to impede function. Good agreement between theoretically predicted fluctuations of individual residues and corresponding x-ray crystallographic temperature factors is found, which lends support to the GNM elucidation of the conformational dynamics of HA by focusing upon a subset of dominant modes. The lowest frequency mode indicates a global torsion of the HA trimer about its longitudinal axis, accompanied by a substantial mobility at the viral membrane connection. This mode is proposed to constitute the dominant molecular mechanism for the translocation and aggregation of HAs, and for the opening and dilation of the fusion pore. The second and third collective modes indicate a global bending, allowing for a large lateral surface exposure, which is likely to facilitate the close association of the viral and endosomal membranes before pore opening. The analysis of kinetically hot residues, in contrast, reveals a localization of energy centered around the HA2 residue Asp112, which apparently triggers the solvent exposure of the fusion peptide.     

Generalized spring tensor models for protein fluctuation dynamics and conformation changes

Background

In the last decade, various coarse-grained elastic network models have been developed to study the large-scale motions of proteins and protein complexes where computer simulations using detailed all-atom models are not feasible. Among these models, the Gaussian Network Model (GNM) and Anisotropic Network Model (ANM) have been widely used. Both models have strengths and limitations. GNM can predict the relative magnitudes of protein fluctuations well, but due to its isotropy assumption, it can not be applied to predict the directions of the fluctuations. In contrast, ANM adds the ability to do the latter, but loses a significant amount of precision in the prediction of the magnitudes.

Results

In this article, we develop a single model, called generalized spring tensor model (STeM), that is able to predict well both the magnitudes and the directions of the fluctuations. Specifically, STeM performs equally well in B-factor predictions as GNM and has the ability to predict the directions of fluctuations as ANM. This is achieved by employing a physically more realistic potential, the Gō-like potential. The potential, which is more sophisticated than that of either GNM or ANM, though adds complexity to the derivation process of the Hessian matrix (which fortunately has been done once for all and the MATLAB code is freely available electronically at http://www.cs.iastate.edu/~gsong/STeM), causes virtually no performance slowdown.

Conclusions

Derived from a physically more realistic potential, STeM proves to be a natural solution in which advantages that used to exist in two separate models, namely GNM and ANM, are achieved in one single model. It thus lightens the burden to work with two separate models and to relate the modes of GNM with those of ANM at times. By examining the contributions of different interaction terms in the Gō potential to the fluctuation dynamics, STeM reveals, (i) a physical explanation for why the distance-dependent, inverse distance square (i.e., Open image in new window ) spring constants perform better than the uniform ones, and (ii), the importance of three-body and four-body interactions to properly modeling protein dynamics.

     

On the relationship between the protein structure and protein dynamics

Recently, we have developed a method (Shih et al., Proteins: Structure, Function, and Bioinformatics 2007;68: 34-38) to compute correlation of fluctuations of proteins. This method, referred to as the protein fixed-point (PFP) model, is based on the positional vectors of atoms issuing from the fixed point, which is the point of the least fluctuations in proteins. One corollary from this model is that atoms lying on the same shell centered at the fixed point will have the same thermal fluctuations. In practice, this model provides a convenient way to compute the average dynamical properties of proteins directly from the geometrical shapes of proteins without the need of any mechanical models, and hence no trajectory integration or sophisticated matrix operations are needed. As a result, it is more efficient than molecular dynamics simulation or normal mode analysis. Though in the previous study the PFP model has been successfully applied to a number of proteins of various folds, it is not clear to what extent this model will be applied. In this article, we have carried out the comprehensive analysis of the PFP model for a dataset comprising 972 high-resolution X-ray structures with pairwise sequence identity or=0.5. Our result shows that the fixed-point model is indeed quite general and will be a useful tool for high throughput analysis of dynamical properties of proteins.     

Network models reveal stability and structural rearrangement of signal recognition particle

The signal recognition particle (SRP) and its receptors (SR) mediate the cotranslational targeting of the membrane and secretory proteins in all cells. In Escherichia coli, SRP is composed of the Ffh protein and the 4.5S SRP RNA. Ffh is a multidomain protein comprising a methionine-rich (M) domain, a helical N domain, and a Ras-like guanine triphosphatase (GTPase) (G) domain. The N and G domains are commonly referred to as one structural unit, the NG domain. In this article, the complex structure of SRP and SR is investigated with the Gaussian network model (GNM) and anisotropic network model (ANM). GNM provides the information of structure stability. It is found that the intermolecular interactions between SRP and SR can obviously decrease the fluctuation of NG domains. Nevertheless, the large structural rearrangement will take place during the cotranslational protein targeting cycle. Hence, the moving directions of fluctuation regions are further ascertained by using cross-correlation analysis and the ANM. The NG domain of Ffh undergoes a clockwise rotation around the GM linker and the M domain of Ffh shows an opposite direction to the NG domain. These functional movements will facilitate the SRP structure to transform into the free form and the sequence-bound form. These simple coarse-grained analyses can be used as a general and quick method for the mechanism studies of protein assembly and supramolecular systems.     

Anisotropy and anharmonicity of atomic fluctuations in proteins: analysis of a molecular dynamics simulation

Positional probability density functions (pdf) for the atomic fluctuations are determined from a molecular dynamics simulation for hen egg-white lysozyme. Most atoms are found to have motions that are highly anisotropic but only slightly anharmonic. The largest deviations from harmonic motion are in the direction of the largest rms fluctuations in the local principal axis frame. Backbone atoms tend to be more nearly harmonic than sidechain atoms. The atoms with the largest anharmonicities tend to have pdfs with multiple peaks, each of which is close to harmonic. Several model pdfs are evaluated on the basis of how well they fit probability densities from the dynamics simulations when parameterized in terms of the moments of the distribution. Gram-Charlier and Edgeworth perturbation expansions, which have been successful in describing the motions of small molecules in crystals, are shown to be inadequate for the distributions found in the dynamics of proteins. Multipeaked distribution functions are found to be more appropriate.     

Escherichia coli adenylate kinase dynamics: comparison of elastic network model modes with mode-coupling (15)N-NMR relaxation data

The dynamics of adenylate kinase of Escherichia coli (AKeco) and its complex with the inhibitor AP(5)A, are characterized by correlating the theoretical results obtained with the Gaussian Network Model (GNM) and the anisotropic network model (ANM) with the order parameters and correlation times obtained with Slowly Relaxing Local Structure (SRLS) analysis of (15)N-NMR relaxation data. The AMPbd and LID domains of AKeco execute in solution large amplitude motions associated with the catalytic reaction Mg(+2)*ATP + AMP --> Mg(+2)*ADP + ADP. Two sets of correlation times and order parameters were determined by NMR/SRLS for AKeco, attributed to slow (nanoseconds) motions with correlation time tau( perpendicular) and low order parameters, and fast (picoseconds) motions with correlation time tau( parallel) and high order parameters. The structural connotation of these patterns is examined herein by subjecting AKeco and AKeco*AP(5)A to GNM analysis, which yields the dynamic spectrum in terms of slow and fast modes. The low/high NMR order parameters correlate with the slow/fast modes of the backbone elucidated with GNM. Likewise, tau( parallel) and tau( perpendicular) are associated with fast and slow GNM modes, respectively. Catalysis-related domain motion of AMPbd and LID in AKeco, occurring per NMR with correlation time tau( perpendicular), is associated with the first and second collective slow (global) GNM modes. The ANM-predicted deformations of the unliganded enzyme conform to the functional reconfiguration induced by ligand-binding, indicating the structural disposition (or potential) of the enzyme to bind its substrates. It is shown that NMR/SRLS and GNM/ANM analyses can be advantageously synthesized to provide insights into the molecular mechanisms that control biological function.     

Network models reveal stability and structural rearrangement of signal recognition particle

The signal recognition particle (SRP) and its receptors (SR) mediate the cotranslational targeting of the membrane and secretory proteins in all cells. In Escherichia coli, SRP is composed of the Ffh protein and the 4.5S SRP RNA. Ffh is a multidomain protein comprising a methionine-rich (M) domain, a helical N domain, and a Ras-like guanine triphosphatase (GTPase) (G) domain. The N and G domains are commonly referred to as one structural unit, the NG domain. In this article, the complex structure of SRP and SR is investigated with the Gaussian network model (GNM) and anisotropic network model (ANM). GNM provides the information of structure stability. It is found that the intermolecular interactions between SRP and SR can obviously decrease the fluctuation of NG domains. Nevertheless, the large structural rearrangement will take place during the cotranslational protein targeting cycle. Hence, the moving directions of fluctuation regions are further ascertained by using cross-correlation analysis and the ANM. The NG domain of Ffh undergoes a clockwise rotation around the GM linker and the M domain of Ffh shows an opposite direction to the NG domain. These functional movements will facilitate the SRP structure to transform into the free form and the sequence-bound form. These simple coarse-grained analyses can be used as a general and quick method for the mechanism studies of protein assembly and supramolecular systems.