 Research
 Open Access
 Published:
Generalized spring tensor models for protein fluctuation dynamics and conformation changes
BMC Structural Biology volume 10, Article number: S3 (2010)
Abstract
Background
In the last decade, various coarsegrained elastic network models have been developed to study the largescale motions of proteins and protein complexes where computer simulations using detailed allatom 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 Bfactor 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 distancedependent, inverse distance square (i.e., ) spring constants perform better than the uniform ones, and (ii), the importance of threebody and fourbody interactions to properly modeling protein dynamics.
Introduction
It is now well accepted that the functions of a protein are closely related to not only its structure but also its dynamics. With the advancement of the computational power and increasing availability of computational resources, functionrelated protein dynamics, such as largescale conformation transitions, has been probed by various computational methods at multiple scales. Among these computational methods, coarsegrained models play an important role since many functional processes take place over time scales that are well beyond the capacity of allatom simulations [1]. One type of coarsegrained models, the elastic network models (ENMs), have been particularly successful and widely used in studying protein dynamics and in relating the intrinsic motions of a protein with its functionalrelated conformation changes over the last decade [2–5].
The reason why ENMs have been well received as compared to the conventional normal mode analysis (NMA) lies at its simplicity to use. ENMs do not require energy minimization and therefore can be applied directly to crystal structures to compute the modes of motions. In contrast, minimization is required for carrying out the conventional normal mode analysis (NMA). The problematic aspect of energy minimization is that it usually shifts the protein molecule away from its crystal conformation by about 2 Å. In addition, in ENMs analytical solutions to residue fluctuations and motion correlations can be easily derived. On the other hand, the simplicity of ENMs leaves much room for improvement and many new models have been proposed [6–12].
The two most widely used ENM models are Gaussian Network Model (GNM) and Anisotropic Network Model (ANM). They have been used to predict the magnitudes or directions of the residue fluctuations from a single structure and have been applied in many research areas [4, 5], such as domain decomposition [13] and allosteric communication [14–17]. Both models have their own advantages and disadvantages. GNM can predict the relative magnitudes of the 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 it loses a significant amount of precision in the prediction of the magnitudes.
Gaussian network model. Gaussian Network Model (GNM) was first introduced in [2] under the assumption that the separation between a pair of residues in the folded protein is Gaussianly distributed. Given its simplicity, the model performs extremely well in predicting the experimental Bfactors. The model represents a protein structure using its C_{ α } atoms. The connectivity among the C_{ α } 's is expressed in Kirchhoff matrix Γ (see Eq. (1). Two C_{ α } 's are considered to be in contact if their distance falls within a certain cutoff distance. The cutoff distance between a pair of residues is the only parameter in the model and is normally set to be 7 Å to 8 Å. Let Δr_{ i } and Δr_{ j } represent the instantaneous fluctuations from equilibrium positions of residues i and j and r_{ ij } and r_{0}_{,ij} be the respective instantaneous and equilibrium distances between residues i and j. The Kirchhoff matrix Γ is:
(1)
where i and j are the indices of the residues and r_{ c } is the cutoff distance.
The simplicity of the Kirchhoff matrix formulation results from the assumption that the fluctuations of each residue are isotropic and Gaussianly distributed along the X, Y and Z directions. The expected value of residue fluctuations, <Δr_{ i }^{2} >, and correlations, <Δr_{ i } · Δr_{ j } >, can be easily obtained from the inverse of the Kirchhoff matrix:
(2)
(3)
where k_{ B } is the Boltzmann constant and T is the temperature. γ is the spring constant. The <Δr_{ i }^{2} > term is directly proportional to the crystallographic Bfactors.
Anisotropic network model. GNM provides only the magnitudes of residue fluctuations. To study the motions of a protein in more details, especially to determine the directions of the fluctuations, normal mode analysis (NMA) is needed. Traditional NMA is allatom based and requires a structure to be first energyminimized before the Hessian matrix and normal modes can be computed, which was rather cumbersome. Even after the energy minimization, the derivation of the Hessian matrix is not easy due to the complicated allatom potential. In Tirion's pioneering work [18], the energy minimization step was removed and a much simpler Hookean potential was used, and yet it was shown that the low frequency normal modes remained mostly accurate. Since then, the Hookean spring potentials have been favored in most coarsegrained C_{α} models [3, 19, 20]. One of such models is best known as Anisotropic Network Model (ANM) [3] since it has anisotropic, directional information of the fluctuations. The potential in ANM has the simplest harmonic form. Assuming that a given structure is at equilibrium, the Hessian matrix (3N×3N) can be derived analytically from such a potential [3]. The 3N×3N Hessian matrix H_{ ANM } can be repartitioned into N×N super elements and each super element is a 3×3 tensor.
(4)
where H_{ i j } is the interaction tensor between residues i and j and can be expressed as:
(5)
Let H^{+} be the pseudo inverse of Hessian matrix H_{ ANM }. The mean square fluctuation <Δr_{ i }^{2} > and correlation <Δr_{ i } · Δr_{ j } > can be calculated by summing over the X, Y and Z components:
(6)
(7)
Strengths and limitations of GNM and ANM. The advantages of ANM or GNM over the conventional NMA lie in several aspects: (i) it is a coarsegrained model and uses the C_{ α }'s to represent the residues in a structure; (ii) it does not require energy minimization and thus can be applied directly to crystal structures to compute the modes of motions; (iii) it provides analytical solutions to the mean square fluctuations and motion correlations.
The limitations of the GNM model. GNM provides only information on the magnitudes of residue fluctuations but no directional information. Therefore, the modes of GNM should not be interpreted as protein motions or components of the motions, since the potential in GNM is not rotationally invariant [21].
The limitations of the ANM model. In contrast to that in GNM, the potential in ANM is based on simple, harmonic Hookean springs and is rotationally invariant. And thus, the modes of ANM do represent the possible modes of protein motions. In doing this, however, ANM loses a significant amount of precision in predicting the magnitudes of the fluctuations. The reason is that, in GNM, the fluctuations in the separation between a pair of residues are assumed to be Gaussianly distributed and isotropic, while in ANM, because only a Hookean spring is attached between a pair of residues i and j, the fluctuation of residue j is constrained only longitudinally along the axis from i to j. The fluctuation is unconstrained transversely. The interaction spring tensor between residues i and j in Eq. (5) becomes the following in the local frame (where the Z axis is along the direction from residues i to j):
(8)
Because the fluctuation of residue j is unconstrained transversely relative to residue i, the fluctuations given by ANM are less realistic than those by GNM, which are assumed to be isotropic. The isotropy in GNM is equivalent to an interaction spring tensor between residues i and j of the following form:
(9)
From the two tensors and given in Eqs. 8 and 9, the causes for the limitations in GNM and ANM are clearly displayed. The unrealisticness in ANM is an artifact resulting from its oversimplified potential. The isotropy assumption of GNM, on the other hand, does a better job than ANM in modeling the effect of residue interactions on the magnitudes of the fluctuations, but gives up completely on representing the anisotropic nature that is intrinsic to all physical forces and interactions, since only the magnitudes of the meansquare fluctuations and crosscorrelations were of concern when GNM was first proposed. Therefore, to overcome the limitations of GNM and ANM, what is needed is a generalized interaction spring tensor that both is anisotropic and can exert more proper constraints on the fluctuations than the ANM tensor does. This calls for a model that has a physically more realistic potential than that of ANM. Since potentials with only twobody interactions can provide only longitudinal constraints, it is necessary to include multibody interactions in the potential in order to have transversal constraints as well. The multibody interactions provide additional diagonal and offdiagonal terms to the interaction spring tensor between residues i and j. For example, by properly including threebody interactions, the interaction spring tensor may look like:
(10)
where k represent the indices of the residues that interact with both residues i and j through threebody interaction S. The first tensor on the right side of the equation represents the twobody interaction, which is similar to , except that the interaction strength T(i, j) depends on residues i and j, and thus may be distancedependent as well.
Our contributions. To overcome the limitations of ANM and GNM, we have developed a generalized spring tensor model for studying protein fluctuation dynamics and conformation changes. It is called generalized spring tensor model, or STeM, for the reason that the interaction between a pair of residues i and j is no longer a Hookean spring that has the tensor form of Eq. (8), but takes a generalized tensor form (similar to that in Eq. 10) that can provide both longitudinal and transversal constraints on a residue's fluctuations relative to its neighbours. We obtain the generalized tensor form by deriving the Hessian matrix from a physically more realistic Gōlike potential (Eq. 11), which has been successfully used in many MD simulations to study protein folding processes and conformation changes [22–24]. In additional to the Hookean spring interactions, the potential includes bond bending and torsional interactions, both of which had been found to be helpful in removing the "tip effect" of the ANM model [9]. The inclusion of the bond bending and torsional interactions is reflected in the generalized tensor spring interaction between residues i and j, in such a way that the tensor now includes not only the twobody interaction between residues i and j, but also threebody and fourbody interactions that involve residues i and j (see Eq. 10). In doing this, the STeM model is able to integrate all the aforementioned attractive features of ANM and GNM and overcome their limitations. Specifically, STeM performs equally well in Bfactors predictions as GNM and has the ability to predict the directions of the fluctuations as ANM. This is accomplished with virtually no performance slowdown. The only potential drawback of this model is the significantly increased complexity in deriving the Hessian matrix. Fortunately, this has been done once for all and the derivation results are available electronically at http://www.cs.iastate.edu/~gsong/STeM. STeM is physically more accurate by explicitly including the bond bending and torsional interactions since they capture the chain behavior of protein molecules, which are neglected in most elastic network models where a protein is treated as an elastic rubber. Therefore, we have reasons to expect this model will further distinguish itself in studying protein dynamics where a correct modeling of bond bending and/or torsional rotations is critical.
Results and discussion
Crystallographic Bfactor prediction
Table 1 shows the correlation coefficients between the experimental and calculated Bfactors of the 111 proteins in the first dataset. The mean values of the correlation coefficients of ANM, GNM, and STeM are 0.53, 0.59, and 0.60 respectively. STeM provides the directional information of the residue fluctuations as ANM and has an accuracy even slightly better than GNM in Bfactor predictions. Figure 1 shows the distributions of the correlation coefficients between the calculated Bfactors and the experimental Bfactors. STeM is the only model in which there are instances where the correlation coefficient is above 0.85 and no instances where the correlation coefficient is below 0.25. This implies that the performance of STeM is more steady than either ANM or GNM. The scatter plot of the correlation coefficients between ANM and STeM in Figure 2 shows that STeM performs better than ANM for 80% of the proteins in the dataset. Protein structures of higher resolution have more accurate data on atom coordinates and Bfactors. We investigate whether our model's performance can be further improved when the dataset used is limited to structures with higher resolution. We select the 12 structures with resolution better than 1.3 A from the first dataset. The mean values of the correlation coefficients of these 12 structures are 0.56, 0.62, and 0.63 for ANM, GNM, and STeM, respectively, which gives an improvement of about 5%6% for all of the three models. Since the improvement is based on a relatively small set of 12 structures, a larger dataset is needed to further examine this potential dependence of Bfactor prediction accuracy on structure quality.
The contributions of different interaction terms to the fluctuations
The Gōlike potential in eq. (11) has four different interaction terms, namely, bond stretching, bond bending, torsional interactions, and the nonbonded interactions. It is of great interest to investigate the relative contributions of these different terms to the agreement with experimental Bfactors. Since only the nonbonded interaction term (V_{ 4 }) is able to provide by itself enough constraints to ensure the Hessian matrix to have no more than six zero eigenvalues, V_{ 4 } is used as the base term for the evaluation of different terms' contributions to the meansquare fluctuations. The Hessian matrix of ANM, denoted by H_{ ANM }, is used as another baseline for comparison purposes. Table 2 lists the contributions of these different terms to the improvement of Bfactor predictions as they are added to the potential.
First, it is seen that the nonbonded interactions, as are present in and H_{ ANM } , play a dominant role in contributing to the Bfactors. This is not surprising since the meansquare fluctuations of a residue are mostly constrained by its interactions with its spatial neighbours, most of which are through nonbonded interactions. What is more interesting is that term alone performs better than H_{ ANM }. This is in agreement with recent results that the performance of Bfactor predictions can be improved by using distancedependent force constants [25, 26]. Particularly, the spring constants that take the form of inverse distance square have been shown to be superior in a recent exhaustive study that experimented with different distancedependent spring constants on a large dataset [10]. The Taylor expansion of the nonbonded interaction term (V_{4}) shows that it has an equivalent spring constant of the form (see Eq. 36), which is exactly proportional to the inverse of the pairwise distance square. Thus, STeM provides a physicsbased explanation for the choice of using inverse square distance spring constants.
The contribution to the improvement in Bfactor predictions from each of the bonded interactions, such as that of bond stretching, is small, as had been pointed out by Bahar et al when GNM was first proposed over a decade ago [2]. However, when the contributions of all of these four terms are added up, they together enable the STeM model to gain a significant improvement over ANM to reach the level of accuracy on a par with GNM.
Conformational change evaluation
It is known that the modes derived from the open form of a structure have better overlaps and correlations with the direction of a protein's conformation change than the ones derived from the closed form [20]. Here we apply the STeM model to study the conformation changes between the open and closed forms of 20 proteins and the open forms are used to calculate the normal modes. Table 3 lists the overlaps and correlations of the observed conformation changes and the indices of the modes that are most involved in the conformation changes. GNM is not considered since it cannot provide directional information. The mean values of the overlaps and correlation coefficients of ANM are 0.49 and 0.61 respectively, and 0.52 and 0.64 respectively for STeM. These amount to an improvement of about 5% for STeM over ANM on both overlap and correlation. Since the results are obtained based on a relatively small set of 20 protein pairs, the significance of the improvement seen here needs to be further tested by conducting a more exhaustive analysis that uses a larger set of proteins and varying parameters, and preferably taking into account the effect of crystal packing as well. We will leave this for future work. It is also worth noting that, in both the overlap and correlation calculations, the modes that are most involved in the conformation change tend to have lower indices in STeM than in ANM (see Table 3), which may imply the modes of STeM be of higher quality than those of ANM.
Conclusions
Protein meansquare fluctuations and conformation changes are two closely related aspects of protein dynamics. However, in the past, two separate groups of models were needed to best explain protein meansquare fluctuations or conformation changes. Specifically, the best models for predicting meansquare fluctuations cannot predict conformation changes, and the models that can predict conformation changes do not have the best performance in predicting meansquare fluctuations. There is thus an obvious gap between the models that work well in predicting one aspect of the dynamics and those in another. Since protein meansquare fluctuations and conformation changes are two closely related dynamic phenomena and share a similar physical origin, we reasoned that models based on a physically more accurate potential should be able to bridge the gap and predict both aspects of the protein dynamics well. Indeed, by using a Gōlike potential, we have successfully developed a spring tensor model (STeM) that is able to singly predict well both meansquare fluctuations and conformation changes. Specifically, STeM performs equally well in Bfactor predictions as GNM and has the ability to predict the directions of fluctuations as ANM.
The new STeM model does come with a cost. As is seen, the derivation process of the Hessian matrix in STeM is much more complex than models using only twobody Hookean potentials, such as those used in ANM. However, the introduced complexity in the potential is necessary in resolving the aforementioned gap that is mainly due to oversimplified potentials and in providing a single, unified model for protein dynamics. Moreover, the derivation process, though more complex, needs to be done only once. Examining the different interaction terms in the Gō potential and their contributions to the agreement with experimental Bfactors provides further benefits. Along the way, we have discovered a physical explanation for why the distancedependent, inverse distance square (i.e., ) spring constants perform better than the uniform ones. The van der Waals interaction term in the potential naturally renders inverse distance square spring constants! By including the bond bending and torsional interactions and their contributions to the improvement in Bfactor predictions, the STeM model confirms the importance of 3body and 4body potentials. The importance of multibody potentials is made even more evident when their contribution to the interaction spring tensor is examined  the multibody potentials are shown to be necessary in providing proper constraints on residue fluctuations, even transversely. It is worth noting that the 3body and 4body potentials introduced through bond bending and torsional interactions only scratch the surface of the full extensity of the multibody potentials since bond bending and torsional interactions are restricted to only consecutive residues along the protein chain. The improvement seen here calls for other generalized spring tensor models that have a thorough treatment of the multibody potentials. Chain breaking, such as that due to missing residues, has a more felt impact on STeM than on ANM or GNM, since the first, the second, and the third terms of the potential used to derive the model are all related to the continuity of the chain. We have not evaluated such impact in the current work but this could be a future research direction and our STeM model would be a proper tool for evaluating the impact of chain breaking on protein motions. STeM does not always outperform ANM in Bfactor predictions  it does better than ANM for 80% of the proteins studied. it would be interesting to find out why this is so. Crystal packing has been known to impact significantly the meansquare fluctuations. Therefore, a proper inclusion of the crystal packing effect may further enhance STeM's performance. Since STeM takes into account bond bending and torsional interactions, it is expected that it should further distinguish itself in studying protein dynamics where a correct modeling of bond bending or torsional rotations is critical, such as in predicting the S^{2} order parameters of NMR structures.
Methods
In this section we will show the derivations of the Hessian matrix from a Gōlike potential proposed by Clementi et al [22].
The Gōlike potential
The Gōlike potential in [22] takes the nonnative and native (equilibrium) conformations as input and it can be divided into four terms. The first term of this Gōlike potential (defined as V_{ 1 } for later use) preserves the chain connectivity. The second (V_{2}) and third terms (V_{3}) define the bond angle and torsional interactions respectively and the last term (V_{4}) is the nonlocal interactions. The Gōlike potential has the following expression:
(11)
In Eq. (11), r and r_{0} represent respectively the instantaneous and equilibrium lengths of the virtual bonds between the C_{ α } atoms of consecutive residues. Similarly, the θ (θ_{o}) and Φ (Φ_{0}) are respectively the instantaneous (equilibrium) virtual bond angles formed by three consecutive residues and the instantaneous (equilibrium) virtual dihedral angles formed by four consecutive residues. The r_{ ij } and r_{0,}_{ ij } represent respectively the instantaneous and equilibrium distances between two nonconsecutive residues i and j. The Gōlike potential in Eq. (11) includes several force parameters (K_{ r } , K_{ θ } , , and ε) and the values of these parameters are taken directly from [22] without any tuning. The values of these parameters are: K_{ r } = 100ε, K_{ θ } = 20ε, = ε, = 0.5ε and ε = 0.36.
Anisotropic fluctuations from the second derivative of the Gōlike potential
Similar to ANM, STeM has a 3N×3N Hessian matrix that can be decomposed into N×N superelements. Each superelement in STeM, H_{ i,j }, is a summation of four 3×3 matrices. The first 3×3 matrix is the contribution from bond stretching. The second and third 3×3 matrices are the contributions from bond bending and torsional rotations respectively. The fourth 3×3 matrix is the contribution from nonlocal contacts.
(12)
The Hessian matrix is the second derivative of the overall potential (equation 11). Let us first consider the first term of the Gōlike potential and let (X_{ i } ,Y_{ i } , Z_{ i }) and (X_{ j } , Y_{ j } , Z_{ j }) be the Cartesian coordinates of two consecutive residues i and j.
(13)
The first and second partial derivatives of V_{1} with respect to the Xdirection of residue i are
(14)
(15)
We will get similar results for the Y and Zdirections of residue i. Since we focus only on the equilibrium fluctuations, we can have r ≅ r^{0} at equilibrium and the first and second partial derivatives of V_{1} can be further simplified to the following expressions.
(16)
(17)
In a similar way, the second crossderivatives have the following form:
(18)
Equations 17 and 18 give the elements of the first 3x3 matrix of the super element H_{ ij } in equation 12. For the diagonal super elements H_{ ii }, equations 17 and 18 are substituted by the following:
(19)
(20)
Now let us consider the second term of the potential in Eq. (11) and let (X_{ i } , Y_{ i } , Z_{ i } ), (X_{ j } , Y_{ j } , Z_{ j } ) and (X_{ k } , Y_{ k } , Z_{ k }) be the Cartesian coordinates of three consecutive residues i, j and k. Suppose θ is the virtual bond angle formed by these three consecutive residues. Since the second term of the potential is V_{ 2 } = K_{ θ } (θ – θ_{0})^{2}, the first and second partial derivatives of V_{2} are
(21)
(22)
Since θ equals θ_{0} at equilibrium, can be further simplified as
(23)
Likewise, becomes
(24)
Let p = (X_{ i } — X_{ j } ,Y_{ i } — Y_{ j } , Z_{ i } — Z_{ j }) and q = (X_{ k } — X_{ j } , Y_{ k } — Y_{ j } , Z_{ k } — Z_{ j } ) and define G as the following.
(25)
The θ can be expressed as
(26)
The partial derivatives of θ are
(27)
(28)
(29)
The derivative of G is
(30)
We can also get and .
(31)
(32)
Combined eq (23),(27) and (30), we can get the following formula.
(33)
Similarly, Combined eq (24),(27), (28), (30) and (31), the second crossderivative becomes
(34)
Following a similar approach, we can get and and these second crossderivatives form the elements of the second 3×3 matrix of the super element H_{ ij } in equation 12.
Due to the complexity of the derivation process of the Hessian matrix for the third (dihedral angle) term of the potential, we omit the derivation process here.The complete derivation is given in Additional File 1.
Finally, let's consider the final (nonlocal contact) term.
(35)
A talor expansion will give us the following form.
(36)
Equation 36 has the same harmonic form as the first term but with a different force constant, so the derivation process is the same as the first term. Therefore, we give only the derivation result here.
(37)
After combining the Hessian matrices from all four terms, we can calculate the pseudo inverse of the final Hessian matrix H. The mean square displacement <Δr_{ i }^{2} > and inter residue correlation <Δr_{ i } · Δr_{ j } > can be calculated by summing the elements over the X, Y and Z directions.
(38)
(39)
The protein sets studied
To evaluate the STeM model, we apply it to compute Bfactors and to study protein conformation changes and compare the results with those computed from ANM and GNM. For Bfactors computations, the protein dataset is from [27] and contains 111 proteins. Two proteins, 1CYO and 5PTP, are removed from the dataset because they no longer exist in the current Protein Data Bank [28]. The proteins in the first dataset all have a resolution that is better than or equal to 2.0 Å. For conformation change studies, the dataset is from [20], which contains 20 pairs of protein structures. Each pair of protein structures have significantly large structure difference from each other.
Evaluation techniques
We used the same evaluation techniques as have been applied before [20, 27]. Specifically, the following three numerical measures are used.
The correlation between the experimental and calculated Bfactors
The linear correlation coefficient between the experimental and calculated Bfactors is calculated using the following formula.
(40)
where x_{ i } and y_{ i } are respectively the experimental and calculated Bfactors of the C_{ α } atom of residue i and and are the mean values. N is the number of residues.
The overlap between the experimental observed conformation changes and the calculated modes
The overlap measures the directional similarity between a conformation change and a calculated mode. The formula for calculating the overlap is
(41)
where e_{ i } is the relative displacement of residue i in a selected mode e and r_{ i } is the conformation displacement of residue i.
The correlation between the experimental observed conformation changes and the calculated modes
The correlation measures the magnitude similarity between a conformation change and a calculated mode. The formula used for calculating the correlation is the same as equation (40), with different meaning for x_{ i } and y_{ i }.
(42)
where x_{ i } is the magnitude of the displacement of residue i in the conformation change and y_{ i } is the magnitude of the displacement of residue i in the selected mode. and are the corresponding mean values.
Abbreviations
 ENM:

Elastic Network Model
 GNM:

Gaussian Network Model
 ANM:

Anisotropic Network Model
 STeM:

Spring Tensor Model
References
 1.
Voth GA: Coarsegraining of condensed phase and biomolecular systems. CRC Press; 2009.
 2.
Bahar I, Atilgan AR, Erman B: Direct evaluation of thermal fluctuations in proteins using a singleparameter harmonic potential. Fold Des 1997, 2: 173–181. 10.1016/S13590278(97)000242
 3.
Atilgan AR, Durell SR, Jernigan RL, Demirel MC, Keskin O, Bahar I: Anisotropy of fluctuation dynamics of proteins with an elastic network model. Biophys. J. 2001, 80: 505–515. 10.1016/S00063495(01)76033X
 4.
Bahar I, Rader A: Coarsegrained normal mode analysis in structural biology. Curr. Opin. Struct. Biol. 2005, 15: 586–592. 10.1016/j.sbi.2005.08.007
 5.
Ma J: Usefulness and limitations of normal mode analysis in modeling dynamics of biomolecular complexes. Structure 2005, 13: 373–380. 10.1016/j.str.2005.02.002
 6.
Ming D, Brüschweiler R: Reorientational contactweighted elastic network model for the prediction of protein dynamics: Comparison with NMR relaxation. Biophys. J. 2006, 90: 3382–3388. 10.1529/biophysj.105.071902
 7.
Song G, Jernigan RL: An enhanced elastic network model to represent the motions of domainswapped proteins. Proteins. 2006, 63: 197–209. 10.1002/prot.20836
 8.
Song G, Jernigan RL: vGNM: a better model for understanding the dynamics of proteins in crystals. J. Mol. Biol. 2007, 369: 880–893. 10.1016/j.jmb.2007.03.059
 9.
Lu M, Poon B, Ma J: A new method for coarsegrained elastic normalmode analysis. J Chem Theory Comput 2006, 2: 464–471. 10.1021/ct050307u
 10.
Yang L, Song G, Jernigan RL: Protein elastic network models and the ranges of cooperativity. Proc Natl Acad Sci U S A 2009, 106: 12347–12352. 10.1073/pnas.0902159106
 11.
Zheng W: A unification of the elastic network model and the Gaussian network model for optimal description of protein conformational motions and fluctuations. Biophys. J. 2008, 94: 3853–3857. 10.1529/biophysj.107.125831
 12.
Stember JN, Wriggers W: Bendtwiststretch model for coarse elastic network simulation of biomolecular motion. J. Chem. Phys. 2009, 131: 074112. 10.1063/1.3167410
 13.
Kundu S, Sorensen DC, Phillips GN Jr: Automatic domain decomposition of proteins by a Gaussian network model. Proteins 2007, 57: 725–733. 10.1002/prot.20268
 14.
Lin T, Song G: Predicting allosteric communication pathways using motion correlation network. In Proceedings of the 7th Asia Pacific Bioinformatics Conference (APBC). Tsinghua University, China; 2009:588–598.
 15.
Bahar I, Chennubhotla I, Tobi D: Intrinsic dynamics of enzymes in the unbound state and relation to allosteric regulation. Curr. Opin. Struct. Biol. 2007, 17: 633–640. 10.1016/j.sbi.2007.09.011
 16.
Zheng W, Brooks B: Idnetification of dynamical correlations within the myosin motor domain by the normal mode analysis of an elastic network model. J. Mol. Biol. 2005, 346: 745–759. 10.1016/j.jmb.2004.12.020
 17.
Yang Z, Májek P, Bahar I: Allosteric transitions of supramolecular systems explored by network models: Application to chaperonin GroEL. PLoS Comput. Biol. 2009, 5: e1000360. 10.1371/journal.pcbi.1000360
 18.
Tirion MM: Large amplitude elastic motions in proteins from a singleparameter, atomic analysis. Phys. Rev. Lett. 1996, 77: 1905–1908. 10.1103/PhysRevLett.77.1905
 19.
Hinsen K: Analysis of domain motions by approximate normal mode calculations. Proteins. 1998, 33: 417–429. 10.1002/(SICI)10970134(19981115)33:3<417::AIDPROT10>3.0.CO;28
 20.
Tama F, Sanejouand Y: Conformational change of proteins arising from normal mode calculations. Protein Eng. 2001, 14: 1–6. 10.1093/protein/14.1.1
 21.
Thorpe MF: Comment on elastic network models and proteins. Phys. Biol. 2007, 4: 60–63. 10.1088/14783975/4/1/N01
 22.
Clementi C, Nymeyer H, Onuchic JN: Topological and energetic factors: What determines the structural details of the transition state ensemble and "enroute" intermediates for protein folding? An investigation for small globular proteins. J. Mol. Biol. 2000, 298: 937–953. 10.1006/jmbi.2000.3693
 23.
Koga N, Takada S: Roles of native topology and chainlength scaling in protein folding: a simulation study with a G olike model. J. Mol. Biol. 2001, 313: 171–180. 10.1006/jmbi.2001.5037
 24.
Koga N, Takada S: Foldingbased molecular simulations reveal mechanisms of the rotary motor FlATPase. Proc. Natl. Acad. Sci. USA 2006, 103: 5367–5372. 10.1073/pnas.0509642103
 25.
Riccardi D, Cui Q, Phillips GN: Application of elastic network models to proteins in the crystalline state. Biophys. J. 2009, 96: 464–475. 10.1016/j.bpj.2008.10.010
 26.
Hinsen K, Petrescu AJ, Dellerue S, BellissentFunel MC, Kneller GR: Harmonicity in slow protein dynamics. Chem. Phys. 2000, 261: 25–37. 10.1016/S03010104(00)002226
 27.
Kundu S, Melton JS, Sorensen DC, Phillips GN Jr: Dynamics of proteins in crystals: Comparison of experiment with simple models. Biophys. J. 2002, 83: 723–732. 10.1016/S00063495(02)75203X
 28.
Berman HM, Westbrook J, Feng Z, Gilliland G, Bhat TN, Weissig H, Shindyalov IN, Bourne PE: The Protein Data Bank. Nucleic Acids Res. 2000, 28: 235–242. 10.1093/nar/28.1.235
Acknowledgements
This article has been published as part of BMC Structural Biology Volume 10 Supplement 1, 2010: Selected articles from the Computational Structural Bioinformatics Workshop 2009. The full contents of the supplement are available online at http://www.biomedcentral.com/14726807/10?issue=S1.
Author information
Additional information
Competing interests
The authors declare that they have no competing interests.
Authors' contributions
TuLiang Lin collaborated with Guang Song on the ideas of the Spring Tensor Model and applied the model to the predictions of the Bfactors and conformation changes. Gunag Song conceived the idea of using a more sophisticated potential and suggested the Gōlike potential. Most of the implementation was done by TuLiang Lin under the supervision of Guang Song. Both authors read and approved the manuscript.
Electronic supplementary material
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Lin, T., Song, G. Generalized spring tensor models for protein fluctuation dynamics and conformation changes. BMC Struct Biol 10, S3 (2010) doi:10.1186/1472680710S1S3
Published
DOI
Keywords
 Hessian Matrix
 Conformation Change
 Normal Mode Analysis
 Elastic Network Model
 Torsional Rotation