Rigid-cluster models of conformational transitions in macromolecular machines and assemblies

Protein Fluctuations in Response to Random External Forces

Elastic network models (ENMs) have been widely used in the last decades to investigate protein motions and dynamics. There the intrinsic fluctuations based on the isolated structures are obtained from the normal modes of these elastic networks, and they generally show good agreement with the B-factors extracted from X-ray crystallographic experiments, which are commonly considered to be indicators of protein flexibility. In this paper, we propose a new approach to analyze protein fluctuations and flexibility, which has a more appropriate physical basis. It is based on the application of random forces to the protein ENM to simulate the effects of collisions of solvent on a protein structure. For this purpose, we consider both the Cα-atom coarse-grained anisotropic network model (ANM) and an elastic network augmented with points included for the crystallized waters. We apply random forces to these protein networks everywhere, as well as only on the protein surface alone. Despite the randomness of the directions of the applied perturbations, the computed average displacements of the protein network show a remarkably good agreement with the experimental B-factors. In particular, for our set of 919 protein structures, we find that the highest correlation with the B-factors is obtained when applying forces to the external surface of the water-augmented ANM (an overall gain of 3% in the Pearson’s coefficient for the entire dataset, with improvements up to 30% for individual proteins), rather than when evaluating the fluctuations obtained from the normal modes of a standard Cα-atom coarse-grained ANM. It follows that protein fluctuations should be considered not just as the intrinsic fluctuations of the internal dynamics, but also equally well as responses to external solvent forces, or as a combination of both.

Disorder-Order Transitions in Conformational Selection of a Peptide by Ebola Virus Nucleoprotein

This study presents parallel-tempering lattice Monte Carlo simulations based on the side-chain-only (SICHO) model for calculating the conformational landscape of a 28-residue intrinsically disordered peptide extracted from the Ebola virus protein VP35. The central issue is the applicability of the SICHO potential energy function and in general coarse-grained (CG) representations of intermediate resolution for modeling large-scale conformational heterogeneity that includes both folded and unstructured peptide states. Crystallographic data shows that the peptide folds in a 4₁₀-helix-turn-3₁₀-helix topology upon complex formation with the Ebola virus nucleoprotein, whereas in isolation, the peptide transitions to a disordered conformational ensemble as observed in circular dichroism experiments. The simulation reveals a potential of mean force that displays conformational diversity along the helix-forming reaction coordinate consistent with disorder-order transitions, yet unexpectedly the bound topology is poorly sampled, and a population shift to an unstructured state incurs a significant free-energy penalty. Applying an elastic network interpolation model suggests a hybrid binding mechanism through conformational selection of the 4₁₀-helix followed by an induced fit of the 3₁₀-helix. A comparison of the CG model with previously reported all-atom CHARMM-based simulations highlights a lattice-based approach that is computationally fast and with the correct parameterization yields good resolution to modeling conformational plasticity.

Conformational Selection of a Polyproline Peptide by Ebola Virus VP30

An adaptive temperature-based replica-exchange simulation of a peptide extracted from the Ebola virus nucleoprotein containing a polyproline sequence motif is reported. The simulation results of applying the CHARMM36m force field with a generalized Born solvent model is presented. Conformational heterogeneity is described by potentials of mean force (PMFs) for a set of reaction coordinates that define the topological fold space. Starting from an extended backbone conformation of the peptide observed in an X-ray crystallographic assembly with the Ebola virus protein VP30, the PMFs report a conformational landscape populated by chain excursions to collapsed states with limited transitions to either an extended fold or a canonical polyproline type II helix. Clustering of the conformations and applying an elastic network interpolation model yield a multistep pathway of conformational selection that minimizes the net transition-state cost from the population hub to the bound state. Related difference between the pathway endpoints taken from the PMFs reveal a significant free-energy penalty in reaching a population shift. To evaluate sequence fitness of the Ebola virus peptide in generating probability distributions, two human sequence variants are modeled and are found to produce profiles that show extensive deviations, thus suggesting either dissimilar binding mechanisms or the lack of recognition by VP30.

Theoretical framework for analyzing structural compliance properties of proteins

We propose methods for directly analyzing structural compliance (SC) properties of elastic network models of proteins, and we also propose methods for extracting information about motion properties from the SC properties. The analysis of SC properties involves describing the relationships between the applied forces and the deformations. When decomposing the motion according to the magnitude of SC (SC mode decomposition), we can obtain information about the motion properties under the assumption that the lower SC mode motions or the softer motions occur easily. For practical applications, the methods are formulated in a general form. The parts where forces are applied and those where deformations are evaluated are separated from each other for enabling the analyses of allosteric interactions between the specified parts. The parts are specified not only by the points but also by the groups of points (the groups are treated as flexible bodies). In addition, we propose methods for quantitatively evaluating the properties based on the screw theory and the considerations of the algebraic structures of the basic equations expressing the SC properties. These methods enable quantitative discussions about the relationships between the SC mode motions and the motions estimated from two different conformations; they also help identify the key parts that play important roles for the motions by comparing the SC properties with those of partially constrained models. As application examples, lactoferrin and ATCase are analyzed. The results show that we can understand their motion properties through their lower SC mode motions or the softer motions.

Symmetrical Parameterization of Rigid Body Transformations for Biomolecular Structures

Assessing preferred relative rigid body position and orientation is important in the description of biomolecular structures (such as proteins) and their interactions. In this article, we extend and apply the "symmetrical parameterization," which we recently introduced in the kinematics community, to address problems in structural biology. We also review parameterization methods that are widely used in structural biology to describe relative rigid body motions (in particular, orientations) as a basis for comparison. The new symmetrical parameterization is useful in describing the relative biomolecular rigid body motions, where the parameters are symmetrical in the sense that the subunits of a complex biomolecular structure are described in the same way for the corresponding motion and its inverse. The properties of this new parameterization, singularity analysis, and inverse kinematics are also investigated in more detail. Finally, parameterization is applied to real biomolecular structures and a potential application to structure modeling of symmetric macromolecules to show the efficacy of the symmetrical parameterization in the field of computational structural biology.

Normal mode-guided transition pathway generation in proteins

The biological function of proteins is closely related to its structural motion. For instance, structurally misfolded proteins do not function properly. Although we are able to experimentally obtain structural information on proteins, it is still challenging to capture their dynamics, such as transition processes. Therefore, we need a simulation method to predict the transition pathways of a protein in order to understand and study large functional deformations. Here, we present a new simulation method called normal mode-guided elastic network interpolation (NGENI) that performs normal modes analysis iteratively to predict transition pathways of proteins. To be more specific, NGENI obtains displacement vectors that determine intermediate structures by interpolating the distance between two end-point conformations, similar to a morphing method called elastic network interpolation. However, the displacement vector is regarded as a linear combination of the normal mode vectors of each intermediate structure, in order to enhance the physical sense of the proposed pathways. As a result, we can generate more reasonable transition pathways geometrically and thermodynamically. By using not only all normal modes, but also in part using only the lowest normal modes, NGENI can still generate reasonable pathways for large deformations in proteins. This study shows that global protein transitions are dominated by collective motion, which means that a few lowest normal modes play an important role in this process. NGENI has considerable merit in terms of computational cost because it is possible to generate transition pathways by partial degrees of freedom, while conventional methods are not capable of this.

Minimum action transition paths connecting minima on an energy surface

Dynamics is essential to the biological functions of many bio-molecules, yet our knowledge of dynamics remains fragmented. Experimental techniques for studying bio-molecules either provide high resolution information on static conformations of the molecule or provide low-resolution, ensemble information that does not shed light on single molecule dynamics. In parallel, bio-molecular dynamics occur at time scale that are not yet attainable through detailed simulation methods. These limitations are especially noticeable when studying transition paths. To address this issue, we report in this paper two methods that derive meaningful trajectories for proteins between two of their conformations. The first method, MinActionPath, uses approximations of the potential energy surface for the molecule to derive an analytical solution of the equations of motion related to the concept of minimum action path. The second method, RelaxPath, follows the same principle of minimum action path but implements a more sophisticated potential, including a mixed elastic potential and a collision term to alleviate steric clashes. Using this new potential, the equations of motion cannot be solved analytically. We have introduced a relaxation method for solving those equations. We describe both the theories behind the two methods and their implementations, focusing on the specific techniques we have used that make those implementations amenable to study large molecular systems. We have illustrated the performance of RelaxPath on simple 2D systems. We have also compared MinActionPath and RelaxPath to other methods for generating transition paths on a well suited test set of large proteins, for which the end points of the trajectories as well as an intermediate conformation between those end points are known. We have shown that RelaxPath outperforms those other methods, including MinActionPath, in its ability to generate trajectories that get close to the known intermediates. We have also shown that the structures along the RelaxPath trajectories remain protein-like. Open source versions of the two programs MinActionPath and RelaxPath are available by request.

A better prediction of conformational changes of proteins using minimally connected network models

Elastic network models have recently been used for studying low-frequency collective motions of proteins. These models simplify the complexity that arises from normal mode analysis by considering a simplified potential involving a few parameters. Two common parameters in most of the elastic network models are cutoff radius and force constant. Although the latter has been studied extensively and even elaborate new models were introduced, for the former usually an ad-hoc cutoff radius is considered. Moreover, the quality of the network models is usually assessed by evaluating their prediction against experimental B-factors. In this work, we consider various common elastic network models with different cutoff radii and assess them by their ability to predict conformational changes of proteins in complexes from unbound to bound state. This prediction is performed by perturbing the unbound structure using a number of low-frequency normal modes of its network model to optimally fit the bound structure. We evaluated a dataset of 30 proteins with distinct unbound and bound structures using this criterion. The results showed that, opposed to the common calibration process based on B-factors, a meaningful relationship exists between the quality of the prediction and model parameters. It was shown that the cutoff radius has a major role in this prediction and minimally connected network models, which use the shortest cutoff radius for which the network is stable, give the best results. It was also shown that by considering the first ten normal modes, the conformational changes can be predicted by about 25 percent. Hence, the evaluation process was extended to the case of considering the contribution of all normal modes in the prediction. The results indicated that minimally connected network models are superior, despite their simplicity, when any number of modes are considered in the prediction.

ClustENM: ENM-Based Sampling of Essential Conformational Space at Full Atomic Resolution

Accurate sampling of conformational space and, in particular, the transitions between functional substates has been a challenge in molecular dynamic (MD) simulations of large biomolecular systems. We developed an Elastic Network Model (ENM)-based computational method, ClustENM, for sampling large conformational changes of biomolecules with various sizes and oligomerization states. ClustENM is an iterative method that combines ENM with energy minimization and clustering steps. It is an unbiased technique, which requires only an initial structure as input, and no information about the target conformation. To test the performance of ClustENM, we applied it to six biomolecular systems: adenylate kinase (AK), calmodulin, p38 MAP kinase, HIV-1 reverse transcriptase (RT), triosephosphate isomerase (TIM), and the 70S ribosomal complex. The generated ensembles of conformers determined at atomic resolution show good agreement with experimental data (979 structures resolved by X-ray and/or NMR) and encompass the subspaces covered in independent MD simulations for TIM, p38, and RT. ClustENM emerges as a computationally efficient tool for characterizing the conformational space of large systems at atomic detail, in addition to generating a representative ensemble of conformers that can be advantageously used in simulating substrate/ligand-binding events.

Ligand Docking to Intermediate and Close-To-Bound Conformers Generated by an Elastic Network Model Based Algorithm for Highly Flexible Proteins

Incorporating receptor flexibility in small ligand-protein docking still poses a challenge for proteins undergoing large conformational changes. In the absence of bound structures, sampling conformers that are accessible by apo state may facilitate docking and drug design studies. For this aim, we developed an unbiased conformational search algorithm, by integrating global modes from elastic network model, clustering and energy minimization with implicit solvation. Our dataset consists of five diverse proteins with apo to complex RMSDs 4.7-15 Å. Applying this iterative algorithm on apo structures, conformers close to the bound-state (RMSD 1.4-3.8 Å), as well as the intermediate states were generated. Dockings to a sequence of conformers consisting of a closed structure and its "parents" up to the apo were performed to compare binding poses on different states of the receptor. For two periplasmic binding proteins and biotin carboxylase that exhibit hinge-type closure of two dynamics domains, the best pose was obtained for the conformer closest to the bound structure (ligand RMSDs 1.5-2 Å). In contrast, the best pose for adenylate kinase corresponded to an intermediate state with partially closed LID domain and open NMP domain, in line with recent studies (ligand RMSD 2.9 Å). The docking of a helical peptide to calmodulin was the most challenging case due to the complexity of its 15 Å transition, for which a two-stage procedure was necessary. The technique was first applied on the extended calmodulin to generate intermediate conformers; then peptide docking and a second generation stage on the complex were performed, which in turn yielded a final peptide RMSD of 2.9 Å. Our algorithm is effective in producing conformational states based on the apo state. This study underlines the importance of such intermediate states for ligand docking to proteins undergoing large transitions.

Ribosome Mechanics Informs about Mechanism

The essential aspects of the ribosome's mechanism can be extracted from coarse-grained simulations, including the ratchet motion, the movement together of critical bases at the decoding center, and movements of the peptide tunnel lining that assist in the expulsion of the synthesized peptide. Because of its large size, coarse graining helps to simplify and to aid in the understanding of its mechanism. Results presented here utilize coarse-grained elastic network modeling to extract the dynamics, and both RNAs and proteins are coarse grained. We review our previous results, showing the well-known ratchet motions and the motions in the peptide tunnel and in the mRNA tunnel. The motions of the lining of the peptide tunnel appear to assist in the expulsion of the growing peptide chain, and clamps at the ends of the mRNA tunnel with three proteins ensure that the mRNA is held tightly during decoding and essential for the helicase activity at the entrance. The entry clamp may also assist in base recognition to ensure proper selection of the incoming tRNA. The overall precision of the ribosome machine-like motions is remarkable.

Human Inducible Hsp70: Structures, Dynamics, and Interdomain Communication from All-Atom Molecular Dynamics Simulations

The 70 kDa human heat shock protein is a major molecular chaperone involved in de novo folding of proteins in vivo and refolding of proteins under stress conditions. Hsp70 is related to several "misfolding diseases" and other major pathologies, such as cancer, and is a target for new therapies. Hsp70 is comprised of two main domains: an N-terminal nucleotide binding domain (NBD) and a C-terminal substrate protein binding domain (SBD). The chaperone function of Hsp70 is based on an allosteric mechanism. Binding of ATP in NBD decreases the affinity of the substrate for SBD, and hydrolysis of ATP is promoted by binding of polypeptide segments in the SBD. No complete structure of human Hsp70 is known. Here, we report two models of human Hsp70, constructed by homology with Saccharomyces cerevisiae cochaperone protein Hsp110 (open model) and with Escherichia coli 70 kDa DnaK (closed model) and relaxed for several tens to hundreds of nanoseconds by using all-atom molecular dynamics simulations in explicit solvent. We obtain two stable states, Hsp70 with SBD open and SBD closed, which agree with experimental and structural information for ATP-Hsp70 and ADP-Hsp70, respectively. The dynamics of the transition from the open to closed states is investigated with a coarse-grained model and normal-mode analysis. The results show that the conformational change between the two states can be represented by a relatively small number of collective modes which involved major conformational changes in the two domains. These modes provide a mechanistic representation of the communication between NBD and SBD and allow us to identify subdomains and residues that appear to have a critical role in the conformational change mechanism that guides the chaperoning cycle of Hsp70.

Multi-scaled normal mode analysis method for dynamics simulation of protein-membrane complexes: A case study of potassium channel gating motion correlations

Membrane proteins play critically important roles in many cellular activities such as ions and small molecule transportation, signal recognition, and transduction. In order to fulfill their functions, these proteins must be placed in different membrane environments and a variety of protein-lipid interactions may affect the behavior of these proteins. One of the key effects of protein-lipid interactions is their ability to change the dynamics status of membrane proteins, thus adjusting their functions. Here, we present a multi-scaled normal mode analysis (mNMA) method to study the dynamics perturbation to the membrane proteins imposed by lipid bi-layer membrane fluctuations. In mNMA, channel proteins are simulated at all-atom level while the membrane is described with a coarse-grained model. mNMA calculations clearly show that channel gating motion can tightly couple with a variety of membrane deformations, including bending and twisting. We then examined bi-channel systems where two channels were separated with different distances. From mNMA calculations, we observed both positive and negative gating correlations between two neighboring channels, and the correlation has a maximum as the channel center-to-center distance is close to 2.5 times of their diameter. This distance is larger than recently found maximum attraction distance between two proteins embedded in membrane which is 1.5 times of the protein size, indicating that membrane fluctuation might impose collective motions among proteins within a larger area. The hybrid resolution feature in mNMA provides atomic dynamics information for key components in the system without costing much computer resource. We expect it to be a conventional simulation tool for ordinary laboratories to study the dynamics of very complicated biological assemblies. The source code is available upon request to the authors.

Computational Analysis of SAXS Data Acquisition

Small-angle x-ray scattering (SAXS) is an experimental biophysical method used for gaining insight into the structure of large biomolecular complexes. Under appropriate chemical conditions, the information obtained from a SAXS experiment can be equated to the pair distribution function, which is the distribution of distances between every pair of points in the complex. Here we develop a mathematical model to calculate the pair distribution function for a structure of known density, and analyze the computational complexity of these calculations. Efficient recursive computation of this forward model is an important step in solving the inverse problem of recovering the three-dimensional density of biomolecular structures from their pair distribution functions. In particular, we show that integrals of products of three spherical-Bessel functions arise naturally in this context. We then develop an algorithm for the efficient recursive computation of these integrals.

Features of large hinge-bending conformational transitions. Prediction of closed structure from open state

We performed a detailed analysis of conformational transition pathways for a set of 10 proteins, which undergo large hinge-bending-type motions with 4-12 Å RMSD (root mean-square distance) between open and closed crystal structures. Anisotropic network model-Monte Carlo (ANM-MC) algorithm generates a targeted pathway between two conformations, where the collective modes from the ANM are used for deformation at each iteration and the conformational energy of the deformed structure is minimized via an MC algorithm. The target structure was approached successfully with an RMSD of 0.9-4.1 Å when a relatively low cutoff radius of 10 Å was used in ANM. Even though one predominant mode (first or second) directed the open-to-closed conformational transition, changes in the dominant mode character were observed for most cases along the transition. By imposing radius of gyration constraint during mode selection, it was possible to predict the closed structure for eight out of 10 proteins (with initial 4.1-7.1 Å and final 1.7-2.9 Å RMSD to target). Deforming along a single mode leads to most successful predictions. Based on the previously reported free energy surface of adenylate kinase, deformations along the first mode produced an energetically favorable path, which was interestingly facilitated by a change in mode shape (resembling second and third modes) at key points. Pathway intermediates are provided in our database of conformational transitions (http://safir.prc.boun.edu.tr/anmmc/method/1).

Generating conformational transitions using the euclidean distance matrix

Elastic network interpolation (ENI) is an efficient method for generating intermediate conformations between two end protein conformations. Its current formulation uses interatomic distance. We show how this can be generalized to interatomic distances-squared. This generalization is part of an effort to study protein dynamics on the set of positive semidefinite (PSD) matrices, which has a rich mathematical structure. We use lattice structures to test this interpolation scheme, and discuss some limitations observed. We conclude with some suggestions for future research.

Insights into mechanism kinematics for protein motion simulation

Background

The high demanding computational requirements necessary to carry out protein motion simulations make it difficult to obtain information related to protein motion. On the one hand, molecular dynamics simulation requires huge computational resources to achieve satisfactory motion simulations. On the other hand, less accurate procedures such as interpolation methods, do not generate realistic morphs from the kinematic point of view. Analyzing a protein’s movement is very similar to serial robots; thus, it is possible to treat the protein chain as a serial mechanism composed of rotational degrees of freedom. Recently, based on this hypothesis, new methodologies have arisen, based on mechanism and robot kinematics, to simulate protein motion. Probabilistic roadmap method, which discretizes the protein configurational space against a scoring function, or the kinetostatic compliance method that minimizes the torques that appear in bonds, aim to simulate protein motion with a reduced computational cost.

Results

In this paper a new viewpoint for protein motion simulation, based on mechanism kinematics is presented. The paper describes a set of methodologies, combining different techniques such as structure normalization normalization processes, simulation algorithms and secondary structure detection procedures. The combination of all these procedures allows to obtain kinematic morphs of proteins achieving a very good computational cost-error rate, while maintaining the biological meaning of the obtained structures and the kinematic viability of the obtained motion.

Conclusions

The procedure presented in this paper, implements different modules to perform the simulation of the conformational change suffered by a protein when exerting its function. The combination of a main simulation procedure assisted by a secondary structure process, and a side chain orientation strategy, allows to obtain a fast and reliable simulations of protein motion.

Exploring the conformational transitions of biomolecular systems using a simple two-state anisotropic network model

Biomolecular conformational transitions are essential to biological functions. Most experimental methods report on the long-lived functional states of biomolecules, but information about the transition pathways between these stable states is generally scarce. Such transitions involve short-lived conformational states that are difficult to detect experimentally. For this reason, computational methods are needed to produce plausible hypothetical transition pathways that can then be probed experimentally. Here we propose a simple and computationally efficient method, called ANMPathway, for constructing a physically reasonable pathway between two endpoints of a conformational transition. We adopt a coarse-grained representation of the protein and construct a two-state potential by combining two elastic network models (ENMs) representative of the experimental structures resolved for the endpoints. The two-state potential has a cusp hypersurface in the configuration space where the energies from both the ENMs are equal. We first search for the minimum energy structure on the cusp hypersurface and then treat it as the transition state. The continuous pathway is subsequently constructed by following the steepest descent energy minimization trajectories starting from the transition state on each side of the cusp hypersurface. Application to several systems of broad biological interest such as adenylate kinase, ATP-driven calcium pump SERCA, leucine transporter and glutamate transporter shows that ANMPathway yields results in good agreement with those from other similar methods and with data obtained from all-atom molecular dynamics simulations, in support of the utility of this simple and efficient approach. Notably the method provides experimentally testable predictions, including the formation of non-native contacts during the transition which we were able to detect in two of the systems we studied. An open-access web server has been created to deliver ANMPathway results.

Many biomolecules are like tiny molecular machines that need to change their shapes and visit many states to perform their biological functions. For a complete molecular understanding of a biological process, one needs to have information on the relevant stable states of the system in question, as well as the pathways by which the system travels from one state to another. We report here an efficient computational method that uses the knowledge of experimental structures of a pair of stable states in order to construct an energetically favoravle pathway between them. We adopt a simple representation of the molecular system by replacing the atoms with beads connected by springs and constructing an energy function with two minima around the end-states. We searched for the structure with highest energy that the system is most likely to visit during the transition and created two paths starting from this structure and proceeding toward the end-states. The combined result of these two paths is the minimum energy pathway between the two stable states. We apply this method to study important structural changes in one enzyme and three large proteins that transport small molecules and ions across the cell membrane.

Multiscale Factorization Method for Simulating Mesoscopic Systems with Atomic Precision

Mesoscopic N-atom systems derive their structural and dynamical properties from processes coupled across multiple scales in space and time. A multiscale method for simulating these systems in the friction dominated regime from the underlying N-atom formulation is presented. The method integrates notions of multiscale analysis, Trotter factorization, and a hypothesis that the momenta conjugate to coarse-grained variables constitute a stationary process on the time scale of coarse-grained dynamics. The method is demonstrated for lactoferrin, nudaurelia capensis omega virus, and human papillomavirus to assess its accuracy.

Modeling protein conformational transitions by a combination of coarse-grained normal mode analysis and robotics-inspired methods

BACKGROUND

Obtaining atomic-scale information about large-amplitude conformational transitions in proteins is a challenging problem for both experimental and computational methods. Such information is, however, important for understanding the mechanisms of interaction of many proteins.

METHODS

This paper presents a computationally efficient approach, combining methods originating from robotics and computational biophysics, to model protein conformational transitions. The ability of normal mode analysis to predict directions of collective, large-amplitude motions is applied to bias the conformational exploration performed by a motion planning algorithm. To reduce the dimension of the problem, normal modes are computed for a coarse-grained elastic network model built on short fragments of three residues. Nevertheless, the validity of intermediate conformations is checked using the all-atom model, which is accurately reconstructed from the coarse-grained one using closed-form inverse kinematics.

RESULTS

Tests on a set of ten proteins demonstrate the ability of the method to model conformational transitions of proteins within a few hours of computing time on a single processor. These results also show that the computing time scales linearly with the protein size, independently of the protein topology. Further experiments on adenylate kinase show that main features of the transition between the open and closed conformations of this protein are well captured in the computed path.

CONCLUSIONS

The proposed method enables the simulation of large-amplitude conformational transitions in proteins using very few computational resources. The resulting paths are a first approximation that can directly provide important information on the molecular mechanisms involved in the conformational transition. This approximation can be subsequently refined and analyzed using state-of-the-art energy models and molecular modeling methods.

Efficient prediction of protein conformational pathways based on the hybrid elastic network model

Various computational models have gained immense attention by analyzing the dynamic characteristics of proteins. Several models have achieved recognition by fulfilling either theoretical or experimental predictions. Nonetheless, each method possesses limitations, mostly in computational outlay and physical reality. These limitations remind us that a new model or paradigm should advance theoretical principles to elucidate more precisely the biological functions of a protein and should increase computational efficiency. With these critical caveats, we have developed a new computational tool that satisfies both physical reality and computational efficiency. In the proposed hybrid elastic network model (HENM), a protein structure is represented as a mixture of rigid clusters and point masses that are connected with linear springs. Harmonic analyses based on the HENM have been performed to generate normal modes and conformational pathways. The results of the hybrid normal mode analyses give new physical insight to the 70S ribosome. The feasibility of the conformational pathways of hybrid elastic network interpolation (HENI) was quantitatively evaluated by comparing three different overlap values proposed in this paper. A remarkable observation is that the obtained mode shapes and conformational pathways are consistent with each other. Our timing results show that HENM has some advantage in computational efficiency over a coarse-grained model, especially for large proteins, even though it takes longer to construct the HENM. Consequently, the proposed HENM will be one of the best alternatives to the conventional coarse-grained ENMs and all-atom based methods (such as molecular dynamics) without loss of physical reality.

Nonicosahedral pathways for capsid expansion

For a significant number of viruses a structural transition of the protein container that encapsulates the viral genome forms an important part of the life cycle and is a prerequisite for the particle becoming infectious. Despite many recent efforts the mechanism of this process is still not fully understood, and a complete characterization of the expansion pathways is still lacking. We present here a coarse-grained model that captures the essential features of the expansion process and allows us to investigate the conditions under which a viral capsid becomes unstable. Based on this model we demonstrate that the structural transitions in icosahedral viral capsids are likely to occur through a low-symmetry cascade of local expansion events spreading in a wavelike manner over the capsid surface.

A mass weighted chemical elastic network model elucidates closed form domain motions in proteins

An elastic network model (ENM), usually Cα coarse-grained one, has been widely used to study protein dynamics as an alternative to classical molecular dynamics simulation. This simple approach dramatically saves the computational cost, but sometimes fails to describe a feasible conformational change due to unrealistically excessive spring connections. To overcome this limitation, we propose a mass-weighted chemical elastic network model (MWCENM) in which the total mass of each residue is assumed to be concentrated on the representative alpha carbon atom and various stiffness values are precisely assigned according to the types of chemical interactions. We test MWCENM on several well-known proteins of which both closed and open conformations are available as well as three α-helix rich proteins. Their normal mode analysis reveals that MWCENM not only generates more plausible conformational changes, especially for closed forms of proteins, but also preserves protein secondary structures thus distinguishing MWCENM from traditional ENMs. In addition, MWCENM also reduces computational burden by using a more sparse stiffness matrix.

Elastic network normal modes provide a basis for protein structure refinement

It is well recognized that thermal motions of atoms in the protein native state, the fluctuations about the minimum of the global free energy, are well reproduced by the simple elastic network models (ENMs) such as the anisotropic network model (ANM). Elastic network models represent protein dynamics as vibrations of a network of nodes (usually represented by positions of the heavy atoms or by the C(α) atoms only for coarse-grained representations) in which the spatially close nodes are connected by harmonic springs. These models provide a reliable representation of the fluctuational dynamics of proteins and RNA, and explain various conformational changes in protein structures including those important for ligand binding. In the present paper, we study the problem of protein structure refinement by analyzing thermal motions of proteins in non-native states. We represent the conformational space close to the native state by a set of decoys generated by the I-TASSER protein structure prediction server utilizing template-free modeling. The protein substates are selected by hierarchical structure clustering. The main finding is that thermal motions for some substates, overlap significantly with the deformations necessary to reach the native state. Additionally, more mobile residues yield higher overlaps with the required deformations than do the less mobile ones. These findings suggest that structural refinement of poorly resolved protein models can be significantly enhanced by reduction of the conformational space to the motions imposed by the dominant normal modes.

Discovering free energy basins for macromolecular systems via guided multiscale simulation

An approach for the automated discovery of low free energy states of macromolecular systems is presented. The method does not involve delineating the entire free energy landscape but proceeds in a sequential free energy minimizing state discovery; i.e., it first discovers one low free energy state and then automatically seeks a distinct neighboring one. These states and the associated ensembles of atomistic configurations are characterized by coarse-grained variables capturing the large-scale structure of the system. A key facet of our approach is the identification of such coarse-grained variables. Evolution of these variables is governed by Langevin dynamics driven by thermal-average forces and mediated by diffusivities, both of which are constructed by an ensemble of short molecular dynamics runs. In the present approach, the thermal-average forces are modified to account for the entropy changes following from our knowledge of the free energy basins already discovered. Such forces guide the system away from the known free energy minima, over free energy barriers, and to a new one. The theory is demonstrated for lactoferrin, known to have multiple energy-minimizing structures. The approach is validated using experimental structures and traditional molecular dynamics. The method can be generalized to enable the interpretation of nanocharacterization data (e.g., ion mobility-mass spectrometry, atomic force microscopy, chemical labeling, and nanopore measurements).

KOSMOS: a universal morph server for nucleic acids, proteins and their complexes

KOSMOS is the first online morph server to be able to address the structural dynamics of DNA/RNA, proteins and even their complexes, such as ribosomes. The key functions of KOSMOS are the harmonic and anharmonic analyses of macromolecules. In the harmonic analysis, normal mode analysis (NMA) based on an elastic network model (ENM) is performed, yielding vibrational modes and B-factor calculations, which provide insight into the potential biological functions of macromolecules based on their structural features. Anharmonic analysis involving elastic network interpolation (ENI) is used to generate plausible transition pathways between two given conformations by optimizing a topology-oriented cost function that guarantees a smooth transition without steric clashes. The quality of the computed pathways is evaluated based on their various facets, including topology, energy cost and compatibility with the NMA results. There are also two unique features of KOSMOS that distinguish it from other morph servers: (i) the versatility in the coarse-graining methods and (ii) the various connection rules in the ENM. The models enable us to analyze macromolecular dynamics with the maximum degrees of freedom by combining a variety of ENMs from full-atom to coarse-grained, backbone and hybrid models with one connection rule, such as distance-cutoff, number-cutoff or chemical-cutoff. KOSMOS is available at http://bioengineering.skku.ac.kr/kosmos.

Rigid-CLL: avoiding constant-distance computations in cell linked-lists algorithms

Many of the existing molecular simulation tools require the efficient identification of the set of nonbonded interacting atoms. This is necessary, for instance, to compute the energy values or the steric contacts between atoms. Cell linked-lists can be used to determine the pairs of atoms closer than a given cutoff distance in asymptotically optimal time. Despite this long-term optimality, many spurious distances are anyway computed with this method. Therefore, several improvements have been proposed, most of them aiming to refine the volume of influence for each atom. Here, we suggest a different improvement strategy based on avoiding to fill cells with those atoms that are always at a constant distance of a given atom. This technique is particularly effective when large groups of the particles in the simulation behave as rigid bodies as it is the case in simplified models considering only few of the degrees of freedom of the molecule. In these cases, the proposed technique can reduce the number of distance computations by more than one order of magnitude, as compared with the standard cell linked-list technique. The benefits of this technique are obtained without incurring in additional computation costs, because it carries out the same operations as the standard cell linked-list algorithm, although in a different order. Since the focus of the technique is the order of the operations, it might be combined with existing improvements based on bounding the volume of influence for each atom.

Generating stereochemically acceptable protein pathways

We describe a new method for rapidly generating stereochemically acceptable pathways in proteins. The method, called geometric targeting, is publicly available at the webserver http://pathways.asu.edu, and includes tools for visualization of the pathway and creating movie files for use in presentations. The user submits an initial structure and a target structure, and a pathway between the two input states is generated automatically. Besides visualization, the structural quality of the pathways makes them useful as input pathways into pathway refinement techniques and further computations. The approach in geometric targeting is to gradually change the system's RMSD relative to the target structure while enforcing a set of geometric constraints. The generated pathways are not minimum free energy pathways, but they are geometrically plausible pathways that maintain good covalent bond distances and angles, keep backbone dihedral angles in allowed Ramachandran regions, avoid eclipsed side-chain torsion angles, avoid non-bonded overlap, and maintain a set of hydrogen bonds and hydrophobic contacts. Resulting pathways for over 20 proteins featuring a wide variety of conformational changes are reported here, including the very large GroEL complex.

Modeling loop entropy

Proteins fold from a highly disordered state into a highly ordered one. Traditionally, the folding problem has been stated as one of predicting "the" tertiary structure from sequential information. However, new evidence suggests that the ensemble of unfolded forms may not be as disordered as once believed, and that the native form of many proteins may not be described by a single conformation, but rather an ensemble of its own. Quantifying the relative disorder in the folded and unfolded ensembles as an entropy difference may therefore shed light on the folding process. One issue that clouds discussions of "entropy" is that many different kinds of entropy can be defined: entropy associated with overall translational and rotational Brownian motion, configurational entropy, vibrational entropy, conformational entropy computed in internal or Cartesian coordinates (which can even be different from each other), conformational entropy computed on a lattice, each of the above with different solvation and solvent models, thermodynamic entropy measured experimentally, etc. The focus of this work is the conformational entropy of coil/loop regions in proteins. New mathematical modeling tools for the approximation of changes in conformational entropy during transition from unfolded to folded ensembles are introduced. In particular, models for computing lower and upper bounds on entropy for polymer models of polypeptide coils both with and without end constraints are presented. The methods reviewed here include kinematics (the mathematics of rigid-body motions), classical statistical mechanics, and information theory.

Group theory and biomolecular conformation: I. Mathematical and computational models

Biological macromolecules, and the complexes that they form, can be described in a variety of ways ranging from quantum mechanical and atomic chemical models, to coarser grained models of secondary structure and domains, to continuum models. At each of these levels, group theory can be used to describe both geometric symmetries and conformational motion. In this survey, a detailed account is provided of how group theory has been applied across computational structural biology to analyze the conformational shape and motion of macromolecules and complexes.

Optimal identification of semi-rigid domains in macromolecules from molecular dynamics simulation

Biological function relies on the fact that biomolecules can switch between different conformations and aggregation states. Such transitions involve a rearrangement of parts of the biomolecules involved that act as dynamic domains. The reliable identification of such domains is thus a key problem in biophysics. In this work we present a method to identify semi-rigid domains based on dynamical data that can be obtained from molecular dynamics simulations or experiments. To this end the average inter-atomic distance-deviations are computed. The resulting matrix is then clustered by a constrained quadratic optimization problem. The reliability and performance of the method are demonstrated for two artificial peptides. Furthermore we correlate the mechanical properties with biological malfunction in three variants of amyloidogenic transthyretin protein, where the method reveals that a pathological mutation destabilizes the natural dimer structure of the protein. Finally the method is used to identify functional domains of the GroEL-GroES chaperone, thus illustrating the efficiency of the method for large biomolecular machines.

Exploring the contribution of collective motions to the dynamics of forced-unfolding in tubulin

Decomposition of the intrinsic dynamics of proteins into collective motions among distant regions of the protein structure provides a physically appealing approach that couples the dynamics of the system with its functional role. The cellular functions of microtubules (an essential component of the cytoskeleton in all eukaryotic cells) depend on their dynamic instability, which is altered by various factors among which applied forces are central. To shed light on the coupling between forces and the dynamic instability of microtubules, we focus on the investigation of the response of the microtubule subunits (tubulin) to applied forces. We address this point by adapting an approach designed to survey correlations for the equilibrium dynamics of proteins to the case of correlations for proteins forced-dynamics. The resulting collective motions in tubulin have a number of functional implications, such as the identification of long-range couplings with a role in blocking the dynamic instability of microtubules. A fundamental implication of our study for the life of a cell is that, to increase the likelihood of unraveling of large cytoskeletal filaments under physiological forces, molecular motors must use a combination of pulling and torsion rather than just pulling.

Multiscale network model for large protein dynamics

Protein dynamics is essential for gaining insight into biological functions of proteins. Although protein dynamics is well delineated by molecular model, the molecular model is computationally prohibited for simulating large protein structures. In this work, we provide a multiscale network model (MNM) that allows the efficient computation on low-frequency normal modes related to structural deformation of proteins as well as dynamic behavior of functional sites. Specifically, MNM consists of two regions, one of which is described as a low-resolution structure, while the other is dictated by a high-resolution structure. The high-resolution regions using all alpha carbons of the protein are mainly binding site parts, which play a critical function in molecules, while the low-resolution parts are constructed from a further coarse-grained model (not using all alpha carbons). The feasibility of MNM to observe the cooperative motion of a protein structure was validated. It was shown that the MNM enables us to understand functional motion of proteins with computational efficiency.

The energy profiles of atomic conformational transition intermediates of adenylate kinase

The elastic network interpolation (ENI) (Kim et al., Biophys J 2002;83:1620-1630) is a computationally efficient and physically realistic method to generate conformational transition intermediates between two forms of a given protein. However it can be asked whether these calculated conformations provide good representatives for these intermediates. In this study, we use ENI to generate conformational transition intermediates between the open form and the closed form of adenylate kinase (AK). Based on C(alpha)-only intermediates, we construct atomic intermediates by grafting all the atoms of known AK structures onto the C(alpha) atoms and then perform CHARMM energy minimization to remove steric conflicts and optimize these intermediate structures. We compare the energy profiles for all intermediates from both the CHARMM force-field and from knowledge-based energy functions. We find that the CHARMM energies can successfully capture the two energy minima representing the open AK and closed AK forms, while the energies computed from the knowledge-based energy functions can detect the local energy minimum representing the closed AK form and show some general features of the transition pathway with a somewhat similar energy profile as the CHARMM energies. The combinatorial extension structural alignment (Shindyalov et al., 1998;11:739-747) and the k-means clustering algorithm are then used to show that known PDB structures closely resemble computed intermediates along the transition pathway.

A novel approach for large-scale polypeptide folding based on elastic networks using continuous optimization

We present a new computationally efficient method for large-scale polypeptide folding using coarse-grained elastic networks and gradient-based continuous optimization techniques. The folding is governed by minimization of energy based on Miyazawa-Jernigan contact potentials. Using this method we are able to substantially reduce the computation time on ordinary desktop computers for simulation of polypeptide folding starting from a fully unfolded state. We compare our results with available native state structures from Protein Data Bank (PDB) for a few de-novo proteins and two natural proteins, Ubiquitin and Lysozyme. Based on our simulations we are able to draw the energy landscape for a small de-novo protein, Chignolin. We also use two well known protein structure prediction software, MODELLER and GROMACS to compare our results. In the end, we show how a modification of normal elastic network model can lead to higher accuracy and lower time required for simulation.

Comparisons of experimental and computed protein anisotropic temperature factors

Because of its appealing simplicity, the anisotropic network model (ANM) has been widely accepted and applied to study many molecular motion problems: such as ribosome motions, the molecular mechanisms of GroEL-GroES function, allosteric changes in hemoglobin, motor-protein motions, and conformational changes in general. However, the validity of the ANM has not been closely examined. In this work, we use ANM to predict the anisotropic temperature factors of proteins obtained from X-ray and NMR data. The rich, directional anisotropic temperature factor data available for hundreds of proteins in the protein data bank are used as validation data to closely test the ANM model. The significance of this work is that it presents a timely, important evaluation of the model, shows the extent of its accuracy in reproducing experimental anisotropic temperature factors, and suggests ways to improve the model. An improved model will help us better understand the internal dynamics of proteins, which in turn can greatly expand the usefulness of the models, which has already been demonstrated in many applications.

Protein elastic network models and the ranges of cooperativity

Elastic network models (ENMs) are entropic models that have demonstrated in many previous studies their abilities to capture overall the important internal motions, with comparisons having been made against crystallographic B-factors and NMR conformational variabilities. ENMs have become an increasingly important tool and have been widely used to comprehend protein dynamics, function, and even conformational changes. However, reliance upon an arbitrary cutoff distance to delimit the range of interactions has presented a drawback for these models, because the optimal cutoff values can differ somewhat from protein to protein and can lead to quirks such as some shuffling in the order of the normal modes when applied to structures that differ only slightly. Here, we have replaced the requirement for a cutoff distance and introduced the more physical concept of inverse power dependence for the interactions, with a set of elastic network models that are parameter-free, with the distance cutoff removed. For small fluctuations about the native forms, the power dependence is the inverse square, but for larger deformations, the power dependence may become inverse 6th or 7th power. These models maintain and enhance the simplicity and generality of the original ENMs, and at the same time yield better predictions of crystallographic B-factors (both isotropic and anisotropic) and of the directions of conformational transitions. Thus, these parameter-free ENMs can be models of choice whenever elastic network models are used.

Large Protein Dynamics Described by Hierarchical-Component Mode Synthesis

Protein dynamics has played a pivotal role in understanding the biological function of protein. For investigation of such dynamics, normal-mode analysis (NMA) has been broadly employed with atomistic model and/or coarse-grained models such as elastic network model (ENM). For large protein complexes, NMA with even ENM encounters the expensive computational process such as diagonalization of Hessian (stiffness) matrix. Here, we suggest the hierarchical-component mode synthesis (hCMS), which allows the fast computation of low-frequency normal modes related to conformational change. Specifically, a large protein structure is regarded as a combination of several structural units, for which the eigen-value problem is utilized for obtaining the frequencies and their normal modes for each structural unit, and consequently, such frequencies and normal modes are assembled with geometrical constraint for interface between structural units in order to find the low-frequency normal modes of a large protein complex. It is shown that hCMS is able to provide the normal modes with accuracy, quantitatively comparable to those of original NMA. This implies that hCMS may enable the computationally efficient analysis of large protein dynamics.

Building alternate protein structures using the elastic network model

We describe a method for efficiently generating ensembles of alternate, all-atom protein structures that (a) differ significantly from the starting structure, (b) have good stereochemistry (bonded geometry), and (c) have good steric properties (absence of atomic overlap). The method uses reconstruction from a series of backbone framework structures that are obtained from a modified elastic network model (ENM) by perturbation along low-frequency normal modes. To ensure good quality backbone frameworks, the single force parameter ENM is modified by introducing two more force parameters to characterize the interaction between the consecutive carbon alphas and those within the same secondary structure domain. The relative stiffness of the three parameters is parameterized to reproduce B-factors, while maintaining good bonded geometry. After parameterization, violations of experimental Calpha-Calpha distances and Calpha-Calpha-Calpha pseudo angles along the backbone are reduced to less than 1%. Simultaneously, the average B-factor correlation coefficient improves to R = 0.77. Two applications illustrate the potential of the approach. (1) 102,051 protein backbones spanning a conformational space of 15 A root mean square deviation were generated from 148 nonredundant proteins in the PDB database, and all-atom models with minimal bonded and nonbonded violations were produced from this ensemble of backbone structures using the SCWRL side chain building program. (2) Improved backbone templates for homology modeling. Fifteen query sequences were each modeled on two targets. For each of the 30 target frameworks, dozens of improved templates could be produced In all cases, improved full atom homology models resulted, of which 50% could be identified blind using the D-Fire statistical potential.

Iterative cluster-NMA: A tool for generating conformational transitions in proteins

Computational models provide insight into the structure-function relationship in proteins. These approaches, especially those based on normal mode analysis, can identify the accessible motion space around a given equilibrium structure. The large magnitude, collective motions identified by these methods are often well aligned with the general direction of the expected conformational transitions. However, these motions cannot realistically be extrapolated beyond the local neighborhood of the starting conformation. In this article, the iterative cluster-NMA (icNMA) method is presented for traversing the energy landscape from a starting conformation to a desired goal conformation. This is accomplished by allowing the evolving geometry of the intermediate structures to define the local accessible motion space, and thus produce an appropriate displacement. Following the derivation of the icNMA method, a set of sample simulations are performed to probe the robustness of the model. A detailed analysis of beta1,4-galactosyltransferase-T1 is also given, to highlight many of the capabilities of icNMA. Remarkably, during the transition, a helix is seen to be extended by an additional turn, emphasizing a new unknown role for secondary structures to absorb slack during transitions. The transition pathway for adenylate kinase, which has been frequently studied in the literature, is also discussed.

Conformational transition pathways explored by Monte Carlo simulation integrated with collective modes

Conformational transitions between open/closed or free/bound states in proteins possess functional importance. We propose a technique in which the collective modes obtained from an anisotropic network model (ANM) are used in conjunction with a Monte Carlo (MC) simulation approach, to investigate conformational transition pathways and pathway intermediates. The ANM-MC technique is applied to adenylate kinase (AK) and hemoglobin. The iterative method, in which normal modes are continuously updated during the simulation, proves successful in accomplishing the transition between open-closed conformations of AK and tense-relaxed forms of hemoglobin (C(alpha)-root mean square deviations between two end structures of 7.13 A and 3.55 A, respectively). Target conformations are reached by root mean-square deviations of 2.27 A and 1.90 A for AK and hemoglobin, respectively. The intermediate conformations overlap with crystal structures from the AK family within a 3.0-A root mean-square deviation. In the case of hemoglobin, the transition of tense-to-relaxed passes through the relaxed state. In both cases, the lowest-frequency modes are effective during transitions. The targeted Monte Carlo approach is used without the application of collective modes. Both the ANM-MC and targeted Monte Carlo techniques can explore sequences of events in transition pathways with an efficient yet realistic conformational search.

Protein conformational transitions explored by mixed elastic network models

We develop a mixed elastic network model (MENM) to study large-scale conformational transitions of proteins between two (or more) known structures. Elastic network potentials for the beginning and end states of a transition are combined, in effect, by adding their respective partition functions. The resulting effective MENM energy function smoothly interpolates between the original surfaces, and retains the beginning and end structures as local minima. Saddle points, transition paths, potentials of mean force, and partition functions can be found efficiently by largely analytic methods. To characterize the protein motions during a conformational transition, we follow "transition paths" on the MENM surface that connect the beginning and end structures and are invariant to parameterizations of the model and the mathematical form of the mixing scheme. As illustrations of the general formalism, we study large-scale conformation changes of the motor proteins KIF1A kinesin and myosin II. We generate possible transition paths for these two proteins that reveal details of their conformational motions. The MENM formalism is computationally efficient and generally applicable even for large protein systems that undergo highly collective structural changes.

Coarse-graining of protein structures for the normal mode studies

The coarse-grained structural model such as Gaussian network has played a vital role in the normal mode studies for understanding protein dynamics related to biological functions. However, for the large proteins, the Gaussian network model is computationally unfavorable for diagonalization of Hessian (stiffness) matrix for the normal mode studies. In this article, we provide the coarse-graining method, referred to as "dynamic model condensation," which enables the further coarse-graining of protein structures consisting of small number of residues. It is shown that the coarser-grained structures reconstructed by dynamic model condensation exhibit the dynamic characteristics, such as low-frequency normal modes, qualitatively comparable to original structures. This sheds light on that dynamic model condensation and may enable one to study the large protein dynamics for gaining insight into biological functions of proteins.

UMass Morph Server: Macromolecular Dynamics Analyses Using Elastic Network Models

The geometry-based mechanical models called elastic network models (ENMs) in various resolutions have been developed for the study of macromolecular motions. In a coarse-grained ENM, a biological system is represented as a network of springs connecting representative points. They range from single atoms to functional domains depending on the level of details in modeling. In this paper presented are the various kinds of coarse-graining methods such as symmetry-constrained, rigid-cluster, and hybrid ENMs. They enable us to overcome the computational burden and memory limitation in the conventional molecular dynamics (MD) simulations and full-atom normal mode analysis (NMA) without loss of generality. For the broad impact of this work on the structural biology area we also develop the UMass Morph Server (UMMS). Based on the requests from online users, UMMS does not only serve a harmonic NMA that describes thermal behaviors (i.e., fluctuations) of a macromolecule around its equilibrium state, but also generates anharmonic transition pathways between two end conformations by using the elastic network interpolation (ENI) also developed by the author. In addition, UMMS can provide two unique features as follows: (i) interpretation of massive MD data by finding essential pathways (ii) the conformation prediction incorporated with time-resolved information such as FRET data. Many example movies and numeric data can be downloadable at http://biomechanics.ecs.umass.edu/umms.html.