Inherent structure analysis of protein folding

An analysis in terms of the inherent structures (IS, local minima) of the multidimensional potential energy landscape is applied to proteins. Detailed calculations are performed for the 46 bead BLN model, which folds into a four-stranded beta-barrel. Enhanced sampling has allowed determination of 239 199 IS states, believed to encompass nearly all the compact, low-energy states, and of well-averaged thermodynamic quantities at low temperature. The density of states shows distinct lobes for compact and extended states, and entropic barriers for the collapse and local ordering transitions. A two-dimensional scatterplot or density of states clearly shows the multifunnel structure of the energy landscape. The anharmonic vibrational free energy is found to play a crucial role in protein folding. The problem of determining the folding transition in a multifunnel system is discussed, and novel indicators of folding are introduced. A particularly clear picture is obtained through the occupation probabilities, pi, of individual low-lying IS, which become finite below the collapse temperature; it is suggested that poor foldability corresponds to a large "misfolding interval" where the excited state pi>0 exceeds that of the native state p0.     

A strategy to find minimal energy nanocluster structures

An unbiased strategy to search for the global and local minimal energy structures of free standing nanoclusters is presented. Our objectives are twofold: to find a diverse set of low lying local minima, as well as the global minimum. To do so, we use massively the fast inertial relaxation engine algorithm as an efficient local minimizer. This procedure turns out to be quite efficient to reach the global minimum, and also most of the local minima. We test the method with the Lennard–Jones (LJ) potential, for which an abundant literature does exist, and obtain novel results, which include a new local minimum for LJ₁₃, 10 new local minima for LJ₁₄, and thousands of new local minima for . Insights on how to choose the initial configurations, analyzing the effectiveness of the method in reaching low‐energy structures, including the global minimum, are developed as a function of the number of atoms of the cluster. Also, a novel characterization of the potential energy surface, analyzing properties of the local minima basins, is provided. The procedure constitutes a promising tool to generate a diverse set of cluster conformations, both two‐ and three‐dimensional, that can be used as an input for refinement by means of ab initio methods. © 2013 Wiley Periodicals, Inc.     

Particle rearrangements during transitions between local minima of the potential energy landscape of a binary Lennard-Jones liquid

The potential energy landscape (PEL) of binary Lennard-Jones (BLJ) mixtures exhibits local minima, or inherent structures (IS), which are organized into metabasins (MBs). We study the particle rearrangements related to transitions between both successive IS and successive MB for a small 80:20 BLJ system near the mode-coupling temperature TMCT. The analysis includes the displacements of individual particles, the localization of the rearrangements, and the relevance of string-like motion. We find that the particle rearrangements during IS and MB transitions do not change significantly at TMCT. In particular, an onset of single particle hopping on the length scale of the interparticle distance is not observed. Further, it is demonstrated that IS and MB dynamics are spatially heterogeneous and facilitated by string-like motion. To investigate the mechanism of string-like motion, we follow the particle rearrangements during suitable sequences of IS transitions. We find that most strings observed after a series of transitions do not move coherently during a single transition, but subunits of different sizes are active at different times. Several findings suggest that, though string-like motion is of comparable relevance when the system explores a MB and when it moves from one MB to another, the occurrence of a successful string enables the system to exit a MB. Moreover, we show that the particle rearrangements during two consecutive MB transitions are basically uncorrelated. In particular, different groups of particles are highly mobile. We further find the positions of strings during successive MB transitions weakly but positively correlated, supporting the idea of dynamic facilitation. Finally, the relation between the features of the potential energy landscape and the relaxation processes in supercooled liquids is discussed.     

Monte Carlo method for computing density of states and quench probability of potential energy and enthalpy landscapes

The thermodynamics and kinetics of a many-body system can be described in terms of a potential energy landscape in multidimensional configuration space. The partition function of such a landscape can be written in terms of a density of states, which can be computed using a variety of Monte Carlo techniques. In this paper, a new self-consistent Monte Carlo method for computing density of states is described that uses importance sampling and a multiplicative update factor to achieve rapid convergence. The technique is then applied to compute the equilibrium quench probability of the various inherent structures (minima) in the landscape. The quench probability depends on both the potential energy of the inherent structure and the volume of its corresponding basin in configuration space. Finally, the methodology is extended to the isothermal-isobaric ensemble in order to compute inherent structure quench probabilities in an enthalpy landscape.     

Preferential attachment during the evolution of a potential energy landscape

It has previously been shown that the network of connected minima on a potential energy landscape is scale-free, and that this reflects a power-law distribution for the areas of the basins of attraction surrounding the minima. Here, the aim is to understand more about the physical origins of these puzzling properties by examining how the potential energy landscape of a 13-atom cluster evolves with the range of the potential. In particular, on decreasing the range of the potential the number of stationary points increases and thus the landscape becomes rougher and the network gets larger. Thus, the evolution of the potential energy landscape can be followed from one with just a single minimum to a complex landscape with many minima and a scale-free pattern of connections. It is found that during this growth process, new edges in the network of connected minima preferentially attach to more highly connected minima, thus leading to the scale-free character. Furthermore, minima that appear when the range of the potential is shorter and the network is larger have smaller basins of attraction. As there are many of these smaller basins because the network grows exponentially, the observed growth process thus also gives rise to a power-law distribution for the hyperareas of the basins.     

Free energy landscape of A-DNA to B-DNA conversion in aqueous solution

The interconversion between the well-characterized A- and B-forms of DNA is a structural transition for which the intermediate states and the free energy difference between the two endpoints are not known precisely. In the present study, the difference between the Root Mean Square Distance (RMSD) from canonical A-form and B-form DNA is used as an order parameter to characterize this free energy difference using umbrella sampling molecular dynamics (MD) simulations with explicit solvent. The constraint imposed along this order parameter allows relatively unrestricted evolution of the intermediate structures away from both canonical A- and B-forms. The free energy difference between the A- and B-forms for the hexamer DNA sequence CTCGAG in aqueous solution is conservatively estimated to be at least 2.8 kcal/mol. A continuum of intermediate structures with no well-defined local minima links the two forms. The absence of any major barriers in the free energy surface is consistent with spontaneous conversion of the A-form DNA to B-form DNA in unconstrained simulations. The extensive sampling in the MD simulations (>0.1 mus) also allowed quantitative energetic characterization of local backbone conformational variables such as sugar pseudorotation angles and BI/BII state equilibria and their dependence on base identity. The absolute minimum in the calculated free energy profile corresponds closely to the crystal structure of the hexamer sequence, indicating that the present method has the potential to identify the most stable state for an arbitrary DNA sequence in water.     

On the ion‐pair dissociation mechanisms in the small NaCl·(H2O)6 cluster: A perspective from reaction path search calculations

We performed reaction path search calculations for the NaCl·(H₂O)₆ cluster using the global reaction route mapping (GRRM) code to understand the atomic‐level mechanisms of the NaCl → Na⁺ + Cl⁻ ionic dissociation induced by water solvents. Low‐lying minima, transition states connecting two local minima and corresponding intrinsic reaction coordinates on the potential energy surface are explored. We found that the Na Cl distances at the transitions states for the dissociation pathways were distributed in a relatively wide range of 2.7–3.7 Å and that the Na Cl distance at the transition state did not correlate with the commonly used solvation coordinates. This suggests that the definition of the transition states with specific structures as well as good reaction coordinate is very difficult for the ionic dissociation process even in a small water cluster. © 2018 Wiley Periodicals, Inc.     

Theoretical studies of the reaction channels on the SO2/OH/NO singlet potential energy surface

Ab initio MP2/6-311++G(2d,2p) investigation of the SO2/OH/NO singlet potential energy surface (PES) has been performed with the aim to localize and describe the existing minima and transition states linking them. The systematic studies have revealed seven minima, with the trans-HONO-SO2 complex (1t) being the global minimum. Eight transition states between minima or between minima and the relevant reactant species have been described. Several available izomerization and dissociation routes have been identified and discussed. The most favorable association of HOSO2 and NO was found to be a barrierless process forming nitrososulfonic acids. Isomerizations between trans-, cis-, and gauche- nitrososulfonic acids (2t, 2c, and 2g) are possible with low-energy barriers. The HOSO2 and NO species can also react via another channels involving high-energy transition states to produce the HOSO-NO2 (3) and HNO-SO3 (4) complexes.     

Free energy landscape analysis of two-dimensional dipolar solvent model at temperatures below and above the rotational freezing point

Ionic solvation in a polar solvent is modeled by a central charge surrounded by dipolar molecules posted on two-dimensional distorted lattice sites with simple rotational dynamics. Density of states is calculated by applying the Wang-Landau algorithm to both the energy and polarization states. The free energy landscapes of solvent molecules as a function of polarization are depicted to explore the competition between the thermal fluctuation and solvation energy. Without a central charge, for temperatures higher than the energy scale of the dipole-dipole interactions, the energy landscape for the small polarization region exhibits a parabolic shape as predicted by Marcus [Rev. Mod. Phys. 65, 599 (1993)] for electron transfer reaction, while there is an additional quartic contribution to the landscape for the large polarization region. When the temperature drops, the simulated free energy landscapes are no longer smooth due to the presence of multiple local minima arising from the frustrated interaction among the dipoles. The parabolic contribution becomes negligible and the energy landscape becomes quartic in shape. For a strong central charge, the energy landscape exhibits an asymmetric profile due to the contributions of linear and cubic terms that arise from the charge-dipole interactions.     

A realistic potential model for N-H vector diffusion in proteins

A realistic model for the potential energy for the diffusion of N-H vectors in a protein is proposed, massively modifying the simplistic models currently used in the literature. In particular, a quantitative and analytical connection between the order parameter of the N-H vector diffusion in a protein and the number of potential minima is established, offering a significant insight into the longstanding question of how protein dynamics is affected by the potential-energy landscape. The largest number of potential minima in a protein is estimated to be no more than around 25. In addition, the conformational entropies derived from classical statistical mechanics and quantum statistical mechanics are proved to be identical. Based on the presented theoretical formula, the number of potential minima for each residue of five representative proteins is evaluated and shows a good correlation between local structural flexibility and the number of potential minima.     

Conformational flexibility and the mechanisms of allosteric transitions in topologically similar proteins

Conformational flexibility plays a central role in allosteric transition of proteins. In this paper, we extend the analysis of our previous study [S. Tripathi and J. J. Portman, Proc. Natl. Acad. Sci. U.S.A. 106, 2104 (2009)] to investigate how relatively minor structural changes of the meta-stable states can significantly influence the conformational flexibility and allosteric transition mechanism. We use the allosteric transitions of the domains of calmodulin as an example system to highlight the relationship between the transition mechanism and the inter-residue contacts present in the meta-stable states. In particular, we focus on the origin of transient local unfolding (cracking), a mechanism that can lower free energy barriers of allosteric transitions, in terms of the inter-residue contacts of the meta-stable states and the pattern of local strain that develops during the transition. We find that the magnitude of the local strain in the protein is not the sole factor determining whether a region will ultimately crack during the transition. These results emphasize that the residue interactions found exclusively in one of the two meta-stable states is the key in understanding the mechanism of allosteric conformational change.     

Hybridizing rapidly exploring random trees and basin hopping yields an improved exploration of energy landscapes

The number of local minima of the potential energy landscape (PEL) of molecular systems generally grows exponentially with the number of degrees of freedom, so that a crucial property of PEL exploration algorithms is their ability to identify local minima, which are low lying and diverse. In this work, we present a new exploration algorithm, retaining the ability of basin hopping (BH) to identify local minima, and that of transition based rapidly exploring random trees (T‐RRT) to foster the exploration of yet unexplored regions. This ability is obtained by interleaving calls to the extension procedures of BH and T‐RRT, and we show tuning the balance between these two types of calls allows the algorithm to focus on low lying regions. Computational efficiency is obtained using state‐of‐the art data structures, in particular for searching approximate nearest neighbors in metric spaces. We present results for the BLN69, a protein model whose conformational space has dimension 207 and whose PEL has been studied exhaustively. On this system, we show that the propensity of our algorithm to explore low lying regions of the landscape significantly outperforms those of BH and T‐RRT. © 2015 Wiley Periodicals, Inc.     

Protein structure prediction using a combination of sequence homology and global energy minimization I. Global energy minimization of surface loops

A procedure has been developed for global energy minimization of surface loops of proteins in the presence of a fixed core. The ECEPP potential function has been modified to allow more accurate representations of hydrogen bond interactions and intrinsic torsional energies. A computationally efficient representation of hydration free energy has been introduced. A local minimization procedure has been developed that uses a cutoff distance, minimization with respect to subsets of degrees of freedom, analytical second derivatives, and distance constraints between rigid segments to achieve efficiency in applications to surface loops. Efficient procedures have been developed for deforming segments of the initial backbone structure and for removing overlaps. Global energy minimization of a surface loop is accomplished by generating a sequence (or a trajectory) of local minima, the component steps of which are generated by searching collections of local minima obtained by deforming seven-residue segments of the surface loop. The search at each component step consists of the following calculations: (1) A large collection of backbone structures is generated by deforming a seven-residue segment of the initial backbone structure. (2) A collection of low-energy backbone structures is generated by applying local energy minimization to the resulting collection of backbone structures (interactions involving side chains that will be searched in this component step are not included in the energy). (3) One low-energy side-chain structure is generated for each of the resulting low-energy backbone structures. (4) A collection of low-energy local minima is generated by applying local energy minimization to the resulting collection of structures. (5) The local minimum with the lowest energy is retained as the next point of the trajectory. Applications of our global search procedure to surface segments of bovine pancreatic trypsin inhibitor (BPTI) and bovine trypsin suggest that component-step searches are reasonably complete. The computational efficiency of component-step searches is such that trajectories consisting of about 10 component steps are feasible using an FPS-5200 array processor. Our procedure for global energy minimization of surface loops is being used to identify and correct problems with the potential function and to calculate protein structure using a combination of sequence homology and global energy minimization.     

Tetrahelix conformations and transformation pathways in Pt1Pd12 clusters

The threshold method is used to explore the potential energy surface of the Pt(1)Pd(12) bimetallic cluster, defined by the Gupta semiempirical potential. A set of helical structures, which follow a Bernal tetrahelix pattern, correspond to local minima for the Pt(1)Pd(12) cluster, characterizing the region of the energy landscape where these structures are present. Both right-handed and left-handed chiral forms were discovered in our searches. Energetic and structural details of each of the tetrahelices are reported as well as the corresponding transition probabilities between these structures and with respect to the icosahedron-shaped global minimum structure via a disconnectivity graph analysis.     

Amyloid Polymorphism in the Protein Folding and Aggregation Energy Landscape

Protein folding involves a large number of steps and conformations in which the folding protein samples different thermodynamic states characterized by local minima. Kinetically trapped on‐ or off‐pathway intermediates are metastable folding intermediates towards the lowest absolute energy minima, which have been postulated to be the natively folded state where intramolecular interactions dominate, and the amyloid state where intermolecular interactions dominate. However, this view largely neglects the rich polymorphism found within amyloid species. We review the protein folding energy landscape in view of recent findings identifying specific transition routes among different amyloid polymorphs. Observed transitions such as twisted ribbon→crystal or helical ribbon→nanotube, and forbidden transitions such helical ribbon?crystal, are discussed and positioned within the protein folding and aggregation energy landscape. Finally, amyloid crystals are identified as the ground state of the protein folding and aggregation energy landscape.     

A study of density of states and ground states in hydrophobic-hydrophilic protein folding models by equi-energy sampling

We propose an equi-energy (EE) sampling approach to study protein folding in the two-dimensional hydrophobic-hydrophilic (HP) lattice model. This approach enables efficient exploration of the global energy landscape and provides accurate estimates of the density of states, which then allows us to conduct a detailed study of the thermodynamics of HP protein folding, in particular, on the temperature dependence of the transition from folding to unfolding and on how sequence composition affects this phenomenon. With no extra cost, this approach also provides estimates on global energy minima and ground states. Without using any prior structural information of the protein the EE sampler is able to find the ground states that match the best known results in most benchmark cases. The numerical results demonstrate it as a powerful method to study lattice protein folding models.     

Internal rotation in NH4(+)-Rg dimers (Rg = He, Ne, Ar): potential energy surfaces and IR spectra of the nu 3 band

The intermolecular potential energy surfaces for the electronic ground states of the ammonium ion-rare gas dimers NH4(+)-He and NH4(+)-Ne are calculated at the MP2 and CCSD(T)/aug-cc-pVXZ (X = D/T/Q) levels of theory. The global minima of both potentials correspond to proton (vertex)-bound structures, Re = 3.13 A, De = 171 cm-1 (He) and Re = 3.21 A, De = 302 cm-1 (Ne). The face- and edge-bound structures are local minima and transition states for the internal rotation dynamics, corresponding to barriers of approximately 20 (He) and 50 cm-1 (Ne). The ab initio potentials are employed in numerical solutions to the rotation-intermolecular vibration Hamiltonian to determine the term values and the rotational and distortion constants for the lowest bound levels in the intramolecular ground vibrational state of both complexes. The results are used to assess the accuracy of two-dimensional (fixed-R) representations of the potentials for determining the internal rotor levels in the ground and nu 3 vibrational states. This model is employed to produce simulations of the IR nu 3 transitions, which are compared to the experimental spectra recorded using photofragmentation spectroscopy. In the case of NH4(+)-Ne the potential parameters are least-squares fitted to the experimental spectrum. The trends within the NH4(+)-Rg series (Rg = He, Ne, Ar) revealed by both the IR spectra and theoretical calculations are discussed.     

Finding pathways between distant local minima

We report a new algorithm for constructing pathways between local minima that involve a large number of intervening transition states on the potential energy surface. A significant improvement in efficiency has been achieved by changing the strategy for choosing successive pairs of local minima that serve as endpoints for the next search. We employ Dijkstra's algorithm [E. W. Dijkstra, Numer. Math. 1, 269 (1959)] to identify the "shortest" path corresponding to missing connections within an evolving database of local minima and the transition states that connect them. The metric employed to determine the shortest missing connection is a function of the minimized Euclidean distance. We present applications to the formation of buckminsterfullerene and to the folding of various biomolecules: the B1 domain of protein G, tryptophan zippers, and the villin headpiece subdomain. The corresponding pathways contain up to 163 transition states and will be used in future discrete path sampling calculations.     

Conformational ensembles and sampled energy landscapes: Analysis and comparison

We present novel algorithms and software addressing four core problems in computational structural biology, namely analyzing a conformational ensemble, comparing two conformational ensembles, analyzing a sampled energy landscape, and comparing two sampled energy landscapes. Using recent developments in computational topology, graph theory, and combinatorial optimization, we make two notable contributions. First, we present a generic algorithm analyzing height fields. We then use this algorithm to perform density‐based clustering of conformations, and to analyze a sampled energy landscape in terms of basins and transitions between them. In both cases, topological persistence is used to manage (geometric) frustration. Second, we introduce two algorithms to compare transition graphs. The first is the classical earth mover distance metric which depends only on local minimum energy configurations along with their statistical weights, while the second incorporates topological constraints inherent to conformational transitions. Illustrations are provided on a simplified protein model (BLN69), whose frustrated potential energy landscape has been thoroughly studied. The software implementing our tools is also made available, and should prove valuable wherever conformational ensembles and energy landscapes are used. © 2015 Wiley Periodicals, Inc.