Normal mode analysis (NMA) has been a powerful tool for studying protein dynamics. Elastic network models (ENM), through their simplicity, have made normal mode computations accessible to a much broader research community and for many more biomolecular systems. The drawback of ENMs, however, is that they are less accurate than NMA. In this work, through steps of simplification that starts with NMA and ends with ENMs we build a tight connection between NMA and ENMs. In the process of bridging between the two, we have also discovered several high‐quality simplified models. Our best simplified model has a mean correlation with the original NMA that is as high as 0.88. In addition, the model is force‐field independent and does not require energy minimization, and thus can be applied directly to experimental structures. Another benefit of drawing the connection is a clearer understanding why ENMs work well and how it can be further improved. We discovered that ANM can be greatly enhanced by including an additional torsional term and a geometry term. Proteins 2014; 82:2157–2168. © 2014 Wiley Periodicals, Inc.
Abstract. It is shown that the density of modes of the vibrational spectrum of globular proteins is universal, i.e., regardless of the protein in question, it closely follows one universal curve. The present study, including 135 proteins analyzed with a full atomic empirical potential (CHARMM22) and using the full complement of all atoms Cartesian degrees of freedom, goes far beyond previous claims of universality, confirming that universality holds even in the frequency range that is well above 100 cm −1 (300 -4000 cm −1 ), where peaks and turns in the density of states are faithfully reproduced from one protein to the next. We also characterize fluctuations of the spectral density from the average, paving the way to a meaningful discussion of rare, unusual spectra and the structural reasons for the deviations in such "outlier" proteins. Since the method used for the derivation of the vibrational modes (potential energy formulation, set of degrees of freedom employed, etc.) has a dramatic effect on the spectral density, another significant implication of our findings is that the universality can provide an exquisite tool for assessing and improving the quality of potential functions and the quality of various models used for NMA computations. Finally, we show that the arXiv:1507.07491v2 [physics.bio-ph] 28 Jan 2016 input configuration too affects the density of modes, thus emphasizing the importance of simplified potential energy formulations that are minimized at the outset. In summary, our findings call for a serious two-way dialogue between theory and experiment: Experimental spectra of proteins could now guide the fine tuning of theoretical empirical potentials, and the various features and peaks observed in theoretical studies -being universal, and hence now rising in importance -would hopefully spur experimental confirmation.
Normal mode analysis (NMA) is an important tool for studying protein dynamics. Because of the complexity of conventional NMA that uses an all-atom model and a semi-empirical force field, many simplified NMA models have been developed, some of which are known as elastic network models. The quality of these simplified NMA models was assessed mostly by evaluating their predictions against experimental B-factors, and rarely by comparing them with the original NMA. In this work, we take the effort to create a publicly accessible dataset of proteins with their minimized structures, NMA modes, and mean-square fluctuations. Then, for the first time, we evaluate the quality of individual normal modes of several widely used elastic network models by comparing them with the conventional NMA. Our results demonstrate that the conventional NMA presents a better and more complete evaluation measure of the quality of elastic network models. This realization should be very helpful in improving current or designing new, higher quality elastic network models. Moreover, using the conventional NMA as the standard of evaluation, a number of interesting and significant insights into the elastic network models are gained.
Dynamics can provide deep insights into the functional mechanisms of proteins and protein complexes. For large protein complexes such as GroEL/GroES with more than 8,000 residues, obtaining a fine-grained all-atom description of its normal mode motions can be computationally prohibitive and is often unnecessary. For this reason, coarse-grained models have been used successfully. However, most existing coarse-grained models use extremely simple potentials to represent the interactions within the coarse-grained structures and as a result, the dynamics obtained for the coarse-grained structures may not always be fully realistic. There is a gap between the quality of the dynamics of the coarse-grained structures given by all-atom models and that by coarse-grained models. In this work, we resolve an important question in protein dynamics computations—how can we efficiently construct coarse-grained models whose description of the dynamics of the coarse-grained structures remains as accurate as that given by all-atom models? Our method takes advantage of the sparseness of the Hessian matrix and achieves a high efficiency with a novel iterative matrix projection approach. The result is highly significant since it can provide descriptions of normal mode motions at an all-atom level of accuracy even for the largest biomolecular complexes. The application of our method to GroEL/GroES offers new insights into the mechanism of this biologically important chaperonin, such as that the conformational transitions of this protein complex in its functional cycle are even more strongly connected to the first few lowest frequency modes than with other coarse-grained models.
It is shown that the density of modes of the vibrational spectrum of globular proteins is universal, i.e., regardless of the protein in question it closely follows one universal curve. The present study, including 135 proteins analyzed with a full atomic empirical potential (CHARMM22) and using the full complement of all atoms Cartesian degrees of freedom, goes far beyond previous claims of universality, confirming that universality holds even in the high-frequency range (300- 4000 1/cm), where peaks and turns in the density of states are faithfully reproduced from one protein to the next. We also characterize fluctuations of the spectral density from the average, paving the way to a meaningful discussion of rare, unusual spectra and the structural reasons for the deviations in such "outlier" proteins. Since the method used for the derivation of the vibrational modes (potential energy formulation, set of degrees of freedom employed, etc.) has a dramatic effect on the spectral density, another significant implication of our findings is that the universality can provide an exquisite tool for assessing and improving the quality of various models used for NMA computations. Finally, we show that the input configuration too affects the density of modes, thus emphasizing the importance of simplified potential energy formulations that are minimized at the outset
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.