Engineering folding dynamics from two-state to downhill: application to λ-repressor

From covalent transition states in chemistry to noncovalent in biology: from β- to Φ-value analysis of protein folding

Solving the mechanism of a chemical reaction requires determining the structures of all the ground states on the pathway and the elusive transition states linking them. 2024 is the centenary of Brønsted's landmark paper that introduced the β-value and structure-activity studies as the only experimental means to infer the structures of transition states. It involves making systematic small changes in the covalent structure of the reactants and analysing changes in activation and equilibrium-free energies. Protein engineering was introduced for an analogous procedure, Φ-value analysis, to analyse the noncovalent interactions in proteins central to biological chemistry. The methodology was developed first by analysing noncovalent interactions in transition states in enzyme catalysis. The mature procedure was then applied to study transition states in the pathway of protein folding - 'part (b) of the protein folding problem'. This review describes the development of Φ-value analysis of transition states and compares and contrasts the interpretation of β- and Φ-values and their limitations. Φ-analysis afforded the first description of transition states in protein folding at the level of individual residues. It revealed the nucleation-condensation folding mechanism of protein domains with the transition state as an expanded, distorted native structure, containing little fully formed secondary structure but many weak tertiary interactions. A spectrum of transition states with various degrees of structural polarisation was then uncovered that spanned from nucleation-condensation to the framework mechanism of fully formed secondary structure. Φ-analysis revealed how movement of the expanded transition state on an energy landscape accommodates the transition from framework to nucleation-condensation mechanisms with a malleability of structure as a unifying feature of folding mechanisms. Such movement follows the rubric of analysis of classical covalent chemical mechanisms that began with Brønsted. Φ-values are used to benchmark computer simulation, and Φ and simulation combine to describe folding pathways at atomic resolution.

The influence of intrinsic folding mechanism of an unfolded protein on the coupled folding-binding process during target recognition

Intrinsically disordered proteins (IDPs) are extensively involved in dynamic signaling processes which require a high association rate and a high dissociation rate for rapid binding/unbinding events and at the same time a sufficient high affinity for specific recognition. Although the coupled folding-binding processes of IDPs have been extensively studied, it is still impossible to predict whether an unfolded protein is suitable for molecular signaling via coupled folding-binding. In this work, we studied the interplay between intrinsic folding mechanisms and coupled folding-binding process for unfolded proteins through molecular dynamics simulations. We first studied the folding process of three representative IDPs with different folded structures, that is, c-Myb, AF9, and E3 rRNase. We found the folding free energy landscapes of IDPs are downhill or show low barriers. To further study the influence of intrinsic folding mechanism on the binding process, we modulated the folding mechanism of barnase via circular permutation and simulated the coupled folding-binding process between unfolded barnase permutant and folded barstar. Although folding of barnase was coupled to target binding, the binding kinetics was significantly affected by the intrinsic folding free energy barrier, where reducing the folding free energy barrier enhances binding rate up to two orders of magnitude. This accelerating effect is different from previous results which reflect the effect of structure flexibility on binding kinetics. Our results suggest that coupling the folding of an unfolded protein with no/low folding free energy barrier with its target binding may provide a way to achieve high specificity and rapid binding/unbinding kinetics simultaneously.

pH-Dependent cooperativity and existence of a dry molten globule in the folding of a miniprotein BBL

Solution pH plays an important role in protein dynamics, stability, and folding; however, detailed mechanisms remain poorly understood. Here we use continuous constant pH molecular dynamics in explicit solvent with pH replica exchange to describe the pH profile of the folding cooperativity of a miniprotein BBL, which has drawn intense debate in the past. Our data reconciled the two opposing hypotheses (downhill vs. two-state) and uncovered a sparsely populated unfolding intermediate. As pH is lowered from 7 to 5, the folding barrier vanishes. As pH continues to decrease, the unfolding barrier lowers and denaturation is triggered by the protonation of Asp162, consistent with experimental evidence. Interestingly, unfolding proceeded via an intermediate, with intact secondary structure and a compact, unlocked hydrophobic core shielded from solvent, lending support to the recent hypothesis of a universal dry molten globule in protein folding. Our work demonstrates that constant pH molecular dynamics is a unique tool for testing this and other hypotheses to advance the knowledge in protein dynamics, stability, and folding.

Molecular dynamics approach to understand the denaturing effect of a millimolar concentration of dodine on a λ-repressor and counteraction by trehalose

Here, we use a molecular dynamics approach to calculate the spatial distribution function of the ternary water–dodine–trehalose (1.0 M) system.

Markov state models of protein misfolding

Markov state models (MSMs) are an extremely useful tool for understanding the conformational dynamics of macromolecules and for analyzing MD simulations in a quantitative fashion. They have been extensively used for peptide and protein folding, for small molecule binding, and for the study of native ensemble dynamics. Here, we adapt the MSM methodology to gain insight into the dynamics of misfolded states. To overcome possible flaws in root-mean-square deviation (RMSD)-based metrics, we introduce a novel discretization approach, based on coarse-grained contact maps. In addition, we extend the MSM methodology to include "sink" states in order to account for the irreversibility (on simulation time scales) of processes like protein misfolding. We apply this method to analyze the mechanism of misfolding of tandem repeats of titin domains, and how it is influenced by confinement in a chaperonin-like cavity.

Reconciling Intermediates in Mechanical Unfolding Experiments with Two-State Protein Folding in Bulk

Most experimentally well-characterized single domain proteins of less than 100 residues have been found to be two-state folders. That is, only two distinct populations can explain both equilibrium and kinetic measurements. Results from single molecule force spectroscopy, where a protein is unfolded by applying a mechanical pulling force to its ends, have largely confirmed this description for proteins found to be two-state in ensemble experiments. Recently, however, stable intermediates have been reported in mechanical unfolding experiments on a cold-shock protein previously found to be a prototypical two-state folder. Here, we tackle this discrepancy using free energy landscapes and Markov state models derived from coarse-grained molecular simulations. We show that protein folding intermediates can be selectively stabilized by the pulling force and that the populations of these intermediates vary in a force-dependent manner. Our model qualitatively captures the experimental results and suggests a possible origin of the apparent discrepancy.

Dependence of internal friction on folding mechanism

An outstanding challenge in protein folding is understanding the origin of "internal friction" in folding dynamics, experimentally identified from the dependence of folding rates on solvent viscosity. A possible origin suggested by simulation is the crossing of local torsion barriers. However, it was unclear why internal friction varied from protein to protein or for different folding barriers of the same protein. Using all-atom simulations with variable solvent viscosity, in conjunction with transition-path sampling to obtain reaction rates and analysis via Markov state models, we are able to determine the internal friction in the folding of several peptides and miniproteins. In agreement with experiment, we find that the folding events with greatest internal friction are those that mainly involve helix formation, while hairpin formation exhibits little or no evidence of friction. Via a careful analysis of folding transition paths, we show that internal friction arises when torsion angle changes are an important part of the folding mechanism near the folding free energy barrier. These results suggest an explanation for the variation of internal friction effects from protein to protein and across the energy landscape of the same protein.