The performance of fine‐grained and coarse‐grained elastic network models and its dependence on various factors

In a recent work we developed a method for deriving accurate simplified models that capture the essentials of conventional all‐atom NMA and identified two best simplified models: ssNMA and eANM, both of which have a significantly higher correlation with NMA in mean square fluctuation calculations than existing elastic network models such as ANM and ANMr2, a variant of ANM that uses the inverse of the squared separation distances as spring constants. Here, we examine closely how the performance of these elastic network models depends on various factors, namely, the presence of hydrogen atoms in the model, the quality of input structures, and the effect of crystal packing. The study reveals the strengths and limitations of these models. Our results indicate that ssNMA and eANM are the best fine‐grained elastic network models but their performance is sensitive to the quality of input structures. When the quality of input structures is poor, ANMr2 is a good alternative for computing mean‐square fluctuations while ANM model is a good alternative for obtaining normal modes. Proteins 2015; 83:1273–1283. © 2015 Wiley Periodicals, Inc.     

Bridging between normal mode analysis and elastic network models

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 can be greatly enhanced by including an additional torsional term and a geometry term. Proteins 2014; 82:2157–2168. © 2014 Wiley Periodicals, Inc.     

Anharmonic Normal Mode Analysis of Elastic Network Model Improves the Modeling of Atomic Fluctuations in Protein Crystal Structures

Protein conformational dynamics, despite its significant anharmonicity, has been widely explored by normal mode analysis (NMA) based on atomic or coarse-grained potential functions. To account for the anharmonic aspects of protein dynamics, this study proposes, and has performed, an anharmonic NMA (ANMA) based on the C_α-only elastic network models, which assume elastic interactions between pairs of residues whose C_α atoms or heavy atoms are within a cutoff distance. The key step of ANMA is to sample an anharmonic potential function along the directions of eigenvectors of the lowest normal modes to determine the mean-squared fluctuations along these directions. ANMA was evaluated based on the modeling of anisotropic displacement parameters (ADPs) from a list of 83 high-resolution protein crystal structures. Significant improvement was found in the modeling of ADPs by ANMA compared with standard NMA. Further improvement in the modeling of ADPs is attained if the interactions between a protein and its crystalline environment are taken into account. In addition, this study has determined the optimal cutoff distances for ADP modeling based on elastic network models, and these agree well with the peaks of the statistical distributions of distances between C_α atoms or heavy atoms derived from a large set of protein crystal structures.     

Anisotropic elastic network modeling of entire microtubules

Microtubules are supramolecular structures that make up the cytoskeleton and strongly affect the mechanical properties of the cell. Within the cytoskeleton filaments, the microtubule (MT) exhibits by far the highest bending stiffness. Bending stiffness depends on the mechanical properties and intermolecular interactions of the tubulin dimers (the MT building blocks). Computational molecular modeling has the potential for obtaining quantitative insights into this area. However, to our knowledge, standard molecular modeling techniques, such as molecular dynamics (MD) and normal mode analysis (NMA), are not yet able to simulate large molecular structures like the MTs; in fact, their possibilities are normally limited to much smaller protein complexes. In this work, we developed a multiscale approach by merging the modeling contribution from MD and NMA. In particular, MD simulations were used to refine the molecular conformation and arrangement of the tubulin dimers inside the MT lattice. Subsequently, NMA was used to investigate the vibrational properties of MTs modeled as an elastic network. The coarse-grain model here developed can describe systems of hundreds of interacting tubulin monomers (corresponding to up to 1,000,000 atoms). In particular, we were able to simulate coarse-grain models of entire MTs, with lengths up to 350 nm. A quantitative mechanical investigation was performed; from the bending and stretching modes, we estimated MT macroscopic properties such as bending stiffness, Young modulus, and persistence length, thus allowing a direct comparison with experimental data.     

Analysis of conformational motions and related key residue interactions responsible for a specific function of proteins with elastic network model

Protein collective motions play a critical role in many biochemical processes. How to predict the functional motions and the related key residue interactions in proteins is important for our understanding in the mechanism of the biochemical processes. Normal mode analysis (NMA) of the elastic network model (ENM) is one of the effective approaches to investigate the structure-encoded motions in proteins. However, the motion modes revealed by the conventional NMA approach do not necessarily correspond to a specific function of protein. In the present work, a new analysis method was proposed to identify the motion modes responsible for a specific function of proteins and then predict the key residue interactions involved in the functional motions by using a perturbation approach. In our method, an internal coordinate that accounts for the specific function was introduced, and the Cartesian coordinate space was transformed into the internal/Cartesian space by using linear approximation, where the introduced internal coordinate serves as one of the axes of the coordinate space. NMA of ENM in this internal/Cartesian space was performed and the function-relevant motion modes were identified according to their contributions to the specific function of proteins. Then the key residue interactions important for the functional motions of the protein were predicted as the interactions whose perturbation largely influences the fluctuation along the internal coordinate. Using our proposed methods, the maltose transporter (MalFGK₂) from E. Coli was studied. The functional motions and the key residue interactions that are related to the channel-gating function of this protein were successfully identified.     

Predicting the order in which contacts are broken during single molecule protein stretching experiments

We combine two methods to enable the prediction of the order in which contacts are broken under external stretching forces in single molecule experiments. These two methods are Gō-like models and elastic network models. The Gō-like models have shown remarkable success in representing many aspects of protein behavior, including the reproduction of experimental data obtained from atomic force microscopy. The simple elastic network models are often used successfully to predict the fluctuations of residues around their mean positions, comparing favorably with the experimentally measured crystallographic B-factors. The behavior of biomolecules under external forces has been demonstrated to depend principally on their elastic properties and the overall shape of their structure. We have studied in detail the muscle protein titin and green fluorescent protein and tested for ten other proteins. First, we stretch the proteins computationally by performing stochastic dynamics simulations with the Gō-like model. We obtain the force-displacement curves and unfolding scenarios of possible mechanical unfolding. We then use the elastic network model to calculate temperature factors (B-factors) and compare the slowest modes of motion for the stretched proteins and compare them with the predicted order of breaking contacts between residues in the Gō-like model. Our results show that a simple Gaussian network model is able to predict contacts that break in the next time stage of stretching. Additionally, we have found that the contact disruption is strictly correlated with the highest force exerted by the backbone on these residues. Our prediction of bond-breaking agrees well with the unfolding scenario obtained with the Gō-like model. We anticipate that this method will be a useful new tool for interpreting stretching experiments.     

Validating and improving elastic network models with molecular dynamics simulations

Elastic network models (ENMs) are a class of simple models intended to represent the collective motions of proteins. In contrast to all‐atom molecular dynamics simulations, the low computational investment required to use an ENM makes them ideal for speculative hypothesis‐testing situations. Historically, ENMs have been validated via comparison to crystallographic B‐factors, but this comparison is relatively low‐resolution and only tests the predictions of relative flexibility. In this work, we systematically validate and optimize a number of ENM‐type models by quantitatively comparing their predictions to microsecond‐scale all‐atom simulations of three different G protein coupled receptors. We show that, despite their apparent simplicity, well‐optimized ENMs perform remarkably well, reproducing the protein fluctuations with an accuracy comparable to what one would expect from all‐atom simulations run for several hundred nanoseconds. Proteins 2010. © 2010 Wiley‐Liss, Inc.     

Effects of surface water on protein dynamics studied by a novel coarse-grained normal mode approach

Normal mode analysis (NMA) has received much attention as a direct approach to extract the collective motions of macromolecules. However, the stringent requirement of computational resources by classical all-atom NMA limits the size of the macromolecules to which the method is normally applied. We implemented a novel coarse-grained normal mode approach based on partitioning the all-atom Hessian matrix into relevant and nonrelevant parts. It is interesting to note that, using classical all-atom NMA results as a reference, we found that this method generates more accurate results than do other coarse-grained approaches, including elastic network model and block normal mode approaches. Moreover, this new method is effective in incorporating the energetic contributions from the nonrelevant atoms, including surface water molecules, into the coarse-grained protein motions. The importance of such improvements is demonstrated by the effect of surface water to shift vibrational modes to higher frequencies and by an increase in overlap of the coarse-grained eigenvector space (the motion directions) with that obtained from molecular dynamics simulations of solvated protein in a water box. These results not only confirm the quality of our method but also point out the importance of incorporating surface structural water in studying protein dynamics.     

Protein structural variation in computational models and crystallographic data

Normal mode analysis offers an efficient way of modeling the conformational flexibility of protein structures. We use anisotropic displacement parameters from crystallography to test the quality of prediction of both the magnitude and directionality of conformational flexibility. Normal modes from four simple elastic network model potentials and from the CHARMM force field are calculated for a data set of 83 diverse, ultrahigh-resolution crystal structures. While all five potentials provide good predictions of the magnitude of flexibility, all-atom potentials have a clear edge at prediction of directionality, and the CHARMM potential has the highest prediction quality. The low-frequency modes from different potentials are similar, but those computed from the CHARMM potential show the greatest difference from the elastic network models. The comprehensive evaluation demonstrates the costs and benefits of using normal mode potentials of varying complexity.     

Opening and closing of a toroidal group II chaperonin revealed by a symmetry constrained elastic network model

Recently, the atomic structures of both the closed and open forms of Group 2 chaperonin protein Mm‐cpn were revealed through crystallography and cryo‐electron microscopy. This toroidal‐like chaperonin is composed of two eightfold rings that face back‐to‐back. To gain a computational advantage, we used a symmetry constrained elastic network model (SCENM), which requires only a repeated subunit structure and its symmetric connectivity to neighboring subunits to simulate the entire system. In the case of chaperonin, only six subunits (i.e., three from each ring) were used out of the eight subunits comprising each ring. A smooth and symmetric pathway between the open and closed conformations was generated by elastic network interpolation (ENI). To support this result, we also performed a symmetry‐constrained normal mode analysis (NMA), which revealed the intrinsic vibration features of the given structures. The NMA and ENI results for the representative single subunit were duplicated according to the symmetry pattern to reconstruct the entire assembly. To test the feasibility of the symmetry model, its results were also compared with those obtained from the full model. This study allowed the folding mechanism of chaperonin Mm‐cpn to be elucidated by SCENM in a timely manner.     

Use of normal modes for structural modeling of proteins: the case study of rat heme oxygenase 1

We present an original approach based on full-atom normal mode analysis (NMA) aimed to expand the general framework of homology modeling. Using the rat heme-free oxygenase 1 as a case system, we show how NMA can be used to model different physiologically relevant conformations of the same protein. Starting from a unique heme-bound X-ray structure, and using two structural templates corresponding to a human and an incomplete rat heme-free structures, we generate models of the rat unbound species with open and closed conformations. Less than 100 lowest frequency modes of the target were sufficient to obtain the heme-free conformations, the closest to the templates. The rat HO-1 model built for the open form shows features similar to the open form of the human heme-free oxygenase, and the one built for the closed form was similar to the incompletely resolved X-ray structure of the same protein available in the Protein DataBank. In the latter case, the use of NMA was particularly useful since it allowed to build a complete structure and therefore to discuss on the reason of the structural differences between open and closed forms. This study shows that the amount of main chain flexibility provided by the normal modes can lead to major improvements in homology modeling approaches. Such applications will allow the characterization of alternative conformations of a target protein with respect to the templates and/or the construction of good quality 3D models based on existing templates with unresolved parts in their tertiary structure.     

A Coarse-Grained Elastic Network Atom Contact Model and Its Use in the Simulation of Protein Dynamics and the Prediction of the Effect of Mutations

Normal mode analysis (NMA) methods are widely used to study dynamic aspects of protein structures. Two critical components of NMA methods are coarse-graining in the level of simplification used to represent protein structures and the choice of potential energy functional form. There is a trade-off between speed and accuracy in different choices. In one extreme one finds accurate but slow molecular-dynamics based methods with all-atom representations and detailed atom potentials. On the other extreme, fast elastic network model (ENM) methods with C_α−only representations and simplified potentials that based on geometry alone, thus oblivious to protein sequence. Here we present ENCoM, an Elastic Network Contact Model that employs a potential energy function that includes a pairwise atom-type non-bonded interaction term and thus makes it possible to consider the effect of the specific nature of amino-acids on dynamics within the context of NMA. ENCoM is as fast as existing ENM methods and outperforms such methods in the generation of conformational ensembles. Here we introduce a new application for NMA methods with the use of ENCoM in the prediction of the effect of mutations on protein stability. While existing methods are based on machine learning or enthalpic considerations, the use of ENCoM, based on vibrational normal modes, is based on entropic considerations. This represents a novel area of application for NMA methods and a novel approach for the prediction of the effect of mutations. We compare ENCoM to a large number of methods in terms of accuracy and self-consistency. We show that the accuracy of ENCoM is comparable to that of the best existing methods. We show that existing methods are biased towards the prediction of destabilizing mutations and that ENCoM is less biased at predicting stabilizing mutations.     

WEBnm@ v2.0: Web server and services for comparing protein flexibility

Background

Normal mode analysis (NMA) using elastic network models is a reliable and cost-effective computational method to characterise protein flexibility and by extension, their dynamics. Further insight into the dynamics–function relationship can be gained by comparing protein motions between protein homologs and functional classifications. This can be achieved by comparing normal modes obtained from sets of evolutionary related proteins.

Results

We have developed an automated tool for comparative NMA of a set of pre-aligned protein structures. The user can submit a sequence alignment in the FASTA format and the corresponding coordinate files in the Protein Data Bank (PDB) format. The computed normalised squared atomic fluctuations and atomic deformation energies of the submitted structures can be easily compared on graphs provided by the web user interface. The web server provides pairwise comparison of the dynamics of all proteins included in the submitted set using two measures: the Root Mean Squared Inner Product and the Bhattacharyya Coefficient. The Comparative Analysis has been implemented on our web server for NMA, WEBnm@, which also provides recently upgraded functionality for NMA of single protein structures. This includes new visualisations of protein motion, visualisation of inter-residue correlations and the analysis of conformational change using the overlap analysis. In addition, programmatic access to WEBnm@ is now available through a SOAP-based web service. Webnm@ is available at http://apps.cbu.uib.no/webnma.

Conclusion

WEBnm@ v2.0 is an online tool offering unique capability for comparative NMA on multiple protein structures. Along with a convenient web interface, powerful computing resources, and several methods for mode analyses, WEBnm@ facilitates the assessment of protein flexibility within protein families and superfamilies. These analyses can give a good view of how the structures move and how the flexibility is conserved over the different structures.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0427-6) contains supplementary material, which is available to authorized users.     

Why are large conformational changes well described by harmonic normal modes?

Low-frequency normal modes generated by elastic network models tend to correlate strongly with large conformational changes of proteins, despite their reliance on the harmonic approximation, which is only valid in close proximity of the native structure. We consider 12 variants of the torsional network model (TNM), an elastic network model in torsion angle space, that adopt different sets of torsion angles as degrees of freedom and reproduce with similar quality the thermal fluctuations of proteins but present drastic differences in their agreement with conformational changes. We show that these differences are related to the extent of the deviations from the harmonic approximation, assessed through an anharmonic energy function whose harmonic approximation coincides with the TNM. Our results indicate that mode anharmonicity is more strongly related to its collectivity, i.e., the number of atoms displaced by the mode, than to its amplitude; low-frequency modes can remain harmonic even at large amplitudes, provided they are sufficiently collective. Finally, we assess the potential benefits of different strategies to minimize the impact of anharmonicity. The reduction of the number of degrees of freedom or their regularization by a torsional harmonic potential significantly improves the collectivity and harmonicity of normal modes and the agreement with conformational changes. In contrast, the correction of normal mode frequencies to partially account for anharmonicity does not yield substantial benefits. The TNM program is freely available at https://github.com/ugobas/tnm.     

Application of Elastic Network Models to Proteins in the Crystalline State

Normal mode analysis using elastic network models has grown popular for probing the low-frequency collective dynamics of proteins and other biomolecular assemblies. In most previous studies, these models were validated by comparing calculated atomic fluctuations for isolated proteins with experimental temperature factors determined in the crystalline state, although there were also hints that including crystal contacts in the calculations has an impact on the comparison. In this study, a set of 83 ultra-high resolution crystal structures with experimentally determined anisotropic displacement parameters is used to evaluate several C_α-based elastic network models that either ignore or treat the crystal environment in different ways; the latter include using periodic boundary conditions defined with respect to the asymmetric unit or the primitive unit cell as well as using the Born-von Kármán boundary condition that accounts for lattice vibrations. For all elastic network models, treating the crystal environment leads to better agreement with experimental anisotropic displacement parameters with the Born-von Kármán boundary condition giving the best agreement. Atomic correlations over the entire protein are clearly affected by the presence of the crystal contacts and fairly sensitive to the way that the crystal environment is treated. These observations highlight the importance of properly treating the protein system in an environment consistent with experiment when either evaluating approximate protein models or using approximate dynamic models in structural refinement application types. Finally, investigation of the scaling behaviors of the cumulative density of states and the heat capacity indicates that there are still gaps between simplified elastic models and all-atom models for proteins.     

Fast assessment of structural models of ion channels based on their predicted current–voltage characteristics

Computational prediction of protein structures is a difficult task, which involves fast and accurate evaluation of candidate model structures. We propose to enhance single‐model quality assessment with a functionality evaluation phase for proteins whose quantitative functional characteristics are known. In particular, this idea can be applied to evaluation of structural models of ion channels, whose main function ‐ conducting ions ‐ can be quantitatively measured with the patch‐clamp technique providing the current–voltage characteristics. The study was performed on a set of KcsA channel models obtained from complete and incomplete contact maps. A fast continuous electrodiffusion model was used for calculating the current–voltage characteristics of structural models. We found that the computed charge selectivity and total current were sensitive to structural and electrostatic quality of models. In practical terms, we show that evaluating predicted conductance values is an appropriate method to eliminate models with an occluded pore or with multiple erroneously created pores. Moreover, filtering models on the basis of their predicted charge selectivity results in a substantial enrichment of the candidate set in highly accurate models. Tests on three other ion channels indicate that, in addition to being a proof of the concept, our function‐oriented single‐model quality assessment method can be directly applied to evaluation of structural models of some classes of protein channels. Finally, our work raises an important question whether a computational validation of functionality should be included in the evaluation process of structural models, whenever possible. Proteins 2016; 84:217–231. © 2015 Wiley Periodicals, Inc.     

Two life prediction models for optoelectronic displays without conventional life tests

To precisely predict the life of optoelectronic displays over a short time, two life prediction models were established based on the three‐parameter Weibull right approximation method (TPWRAM). In Model I, the acceleration life under each stress was obtained using TPWRAM, the data points formed by acceleration life and acceleration stress were, respectively, fitted by three extrapolation functions and the optimal extrapolation function was determined by comparing the fitting determination coefficient and root‐mean‐square error. In Model II, after introducing an acceleration parameter, the luminance attenuation data under conventional stress were calculated directly by combining the ones obtained using TPWRAM under each acceleration stress. The luminance attenuation test data from the vacuum fluorescent display (VFD) were collected through four groups of constant‐stress accelerated degradation tests (ADT), and the two models were applied to the life prediction of VFD. The results indicated that the designed ADT scheme was feasible, Model I revealed the changing law of life with stress and simplified the process of life prediction, and Model II made it possible to obtain the luminance attenuation formula at conventional stress without conducting a conventional life test, overcoming the shortcomings of long time‐consuming traditional life tests. It was verified by comparing the life prediction values that the two models had very high precision, and both of these not only achieve the accurate estimation of optoelectronic product life without resorting to conventional life test, but also improved the method of life prediction and perfect its theoretical system.     

Normal mode analysis with molecular geometry restraints: bridging molecular mechanics and elastic models

A new method for normal mode analysis is reported for all-atom structures using molecular geometry restraints (MGR). Similar to common molecular mechanics force fields, the MGR potential contains short- and long-range terms. The short-range terms are defined by molecular geometry, i.e., bond lengths, angles and dihedrals; the long-range term is similar to that in elastic network models. Each interaction term uses a single force constant parameter, and is determined by fitting against a set of known structures. Tests on proteins/non-proteins show that MGR can produce low frequency eigenvectors closer to all-atom force-field-based methods than conventional elastic network models. Moreover, the “tip effect”, found in low frequency eigenvectors in elastic network models, is reduced in MGR to the same level of the modes produced by force-field-based methods. The results suggest that molecular geometry plays an important role, in addition to molecular shape, in determining low frequency deformational motions. MGR does not require initial energy minimization, and is applicable to almost any structure, including the one with missing atoms, bad contacts, or bad geometries, frequently observed in low-resolution structure determination and refinement. The method bridges the two major representations in normal mode analyses, i.e., the molecular mechanics models and elastic network models.     

Elastic network normal mode dynamics reveal the GPCR activation mechanism

G‐protein‐coupled receptors (GPCR) are a family of membrane‐embedded metabotropic receptors which translate extracellular ligand binding into an intracellular response. Here, we calculate the motion of several GPCR family members such as the M2 and M3 muscarinic acetylcholine receptors, the A_2A adenosine receptor, the β₂‐adrenergic receptor, and the CXCR4 chemokine receptor using elastic network normal modes. The normal modes reveal a dilation and a contraction of the GPCR vestibule associated with ligand passage, and activation, respectively. Contraction of the vestibule on the extracellular side is correlated with cavity formation of the G‐protein binding pocket on the intracellular side, which initiates intracellular signaling. Interestingly, the normal modes of rhodopsin do not correlate well with the motion of other GPCR family members. Electrostatic potential calculation of the GPCRs reveal a negatively charged field around the ligand binding site acting as a siphon to draw‐in positively charged ligands on the membrane surface. Altogether, these results expose the GPCR activation mechanism and show how conformational changes on the cell surface side of the receptor are allosterically translated into structural changes on the inside. Proteins 2014; 82:579–586. © 2013 Wiley Periodicals, Inc.     

A Comparison Between Elastic Network Interpolation and MD Simulation of 16S Ribosomal RNA

Abstract

In this paper a coarse-grained method called elastic network interpolation (ENI) is used to generate feasible transition pathways between two given conformations of the core central domain of 16S Ribosomal RNA (16S rRNA). The two given conformations are the extremes generated by a molecular dynamics (MD) simulation, which differ from each other by 10Å in root-mean-square deviation (RMSD). It takes only several hours to build an ENI pathway on a 1.5GHz Pentium with 512 MB memory, while the MD takes several weeks on high-performance multi-processor servers such as the SGI ORIGIN 2000/2100. It is shown that multiple ENI pathways capture the essential anharmonic motions of millions of timesteps in a particular MD simulation. A coarse-grained normal mode analysis (NMA) is performed on each intermediate ENI conformation, and the lowest 1% of the normal modes (representing about 40 degrees of freedom (DOF)) are used to parameterize fluctuations. This combined ENI/NMA method captures all intermediate conformations in the MD run with 1.5Å RMSD on average. In addition, if we restrict attention to the time interval of the MD run between the two extreme conformations, the RMSD between the closest ENI/NMA pathway and the MD results is about 1Å. These results may serve as a paradigm for reducedDOF dynamic simulations of large biological macromolecules as well as a method for the reduced-parameter interpretation of massive amounts of MD data.