Cooperative Folding Research Articles

The folding and misfolding of multidomain proteins.

Protein folding represents a vital process for any living organism. While significant insights have been gained from studying single-domain proteins, our current knowledge on the folding mechanisms of multidomain proteins remains relatively limited, primarily due to their inherent complexity. The principal aim of this review lies in summarizing the emerging view pertaining multi-domain folding, emphasizing their modular nature, which minimizes misfolding and facilitates evolutionary innovation. We discuss the energetic interplay between domains, highlighting particularly the cases where domain interactions lead to transient misfolded intermediates. These interactions can result in diverse effects, including cooperative folding and domain-specific perturbations, which are particularly relevant to the pathogenesis of neurodegenerative diseases like polyglutamine disorders. The review underscores the critical need to understand multidomain folding, to better comprehend and potentially mitigate the molecular underpinnings of protein misfolding diseases.

Dissecting the Roles of Electrostatic Interactions in Modulating the Folding Stability and Cooperativity of Engrailed Homeodomain.

Engrailed homeodomain (EngHD), a highly charged transcription factor regulating over 200 genes, is a fast-folding protein. Recent studies have shown that the abundant charged residues in EngHD not only facilitate protein-DNA interactions but also influence the conformational disorder of its native structure. However, the mechanisms by which electrostatic interactions modulate the folding of EngHD remain unclear. Here, we employ a coarse-grained structure-based model that incorporates the salt-dependent Debye-Hückel model to investigate the (un)folding behavior of EngHD under various salt concentrations. Our findings demonstrate that increasing salt concentrations enhance both folding stability and cooperativity, while the folding barrier height remains relatively constant due to the distinct electrostatic effects on individual residues. By modulating the energetic balance between local and nonlocal interactions, we shift the folding of EngHD from a downhill process to a two-state process. Notably, we observe a nonmonotonic relationship between the strength of local interactions and residue-level coupling degree during (un)folding, likely attributed to the repulsive electrostatic interactions present in the native structure of EngHD. Additionally, we identify a critical turning point in the dependence of folding cooperativity on salt concentration, classified by the energetic balance of local and nonlocal interactions. Our results provide valuable insights into how electrostatic interactions influence the folding of EngHD, contributing to the theoretical framework for engineering highly charged proteins.

MFIB 2.0: a major update of the database of protein complexes formed by mutual folding of the constituting protein chains.

While the majority of proteins with available structures are able to fold independently and mediate interactions only after acquiring their folded state, a subset of the known protein complexes contains protein chains that are intrinsically disordered in isolation. The Mutual Folding Induced by Binding (MFIB) database collects and classifies protein complexes, wherein all constituent protein chains would be unstable/disordered in isolation but fold into a well-defined 3D complex structure upon binding. This phenomenon is often termed as cooperative folding and binding or mutual synergistic folding (MSF). Here we present a major update to the database: we collected and annotated hundreds of new protein complexes fulfilling the criteria of MSF, leading to an almost six-fold increase in the size of the database. Many novel features have also been introduced, such as clustering of the complexes based on structural similarity and domain types, assigning different evidence levels to each entry and adding the evidence coverage label that allowed us to include complexes of multi(sub)domain monomers with partial MSF. The MFIB 2.0 database is available at https://mfib.pbrg.hu.

Scaling Laws for Protein Folding under Confinement.

Spatial confinement significantly affects protein folding. Without the confinement provided by chaperones, many proteins cannot fold correctly. However, the quantitative effect of confinement on protein folding remains elusive. In this study, we observed scaling laws between the variation in folding transition temperature and the size of confinement, (Tf - Tfbulk)/Tfbulk ∼ L-ν. The scaling exponent v is significantly influenced by both the protein's topology and folding cooperativity. Specifically, for a given protein, v can decrease as the folding cooperativity of the model increases, primarily due to the heightened sensitivity of the unfolded state energy to changes in cage size. For proteins with diverse topologies, variations in topological complexity influence scaling exponents in multiple ways. Notably, v exhibits a clear positive correlation with contact order and the proportion of nonlocal contacts, as this complexity significantly enhances the sensitivity of entropy loss in the unfolded state. Furthermore, we developed a novel scaling argument yielding 5/3 ≤ ν ≤ 10/3, consistent with the simulation results.

Cooperative folding as a molecular switch in an evolved antibody binder

Designing proteins with tunable activities from easily accessible external cues remains a biotechnological challenge. Here, we set out to create a small antibody-binding domain equipped with a molecular switch inspired by the allosteric response to calcium seen in naturally derived proteins like calmodulin. We have focused on one of the three domains of Protein G that show inherent affinity to antibodies. By combining a semi-rational protein design with directed evolution, we engineered novel variants containing a calcium-binding loop rendering the inherent antibody affinity calcium-dependent. The evolved variants resulted from a designed selection strategy subjecting them to negative and positive selection pressures focused on conditional antibody-binding. Hence, these variants contained molecular “on/off” switches, controlling the target affinity towards antibody fragments simply by the presence or absence of calcium. From NMR spectroscopy we found that the molecular mechanism underlying the evolved switching behavior was a coupled calcium-binding and folding event where the target binding surface was intact and functional only in the presence of bound calcium. Notably, it was observed that the response to the employed selection pressures gave rise to the evolution of a cooperative folding mechanism. This observation illustrates why the cooperative folding reaction is an effective solution seen repeatedly in the natural evolution of fine-tuned macromolecular recognition. Engineering binding moieties to confer conditional target interaction has great potential due to the exquisite interaction control that is tunable to application requirements. Improved understanding of the molecular mechanisms behind regulated interactions is crucial to unlock how to engineer switchable proteins useful in a variety of biotechnological applications.

Dominant Chemical Interactions Governing the Folding Mechanism of Oligopeptides.

The hydrophobic effect is the main factor that drives the folding of polypeptide chains. In this study, we have examined the influence of the hydrophobic effect in the context of the main mechanical forces approach, mainly in relation to the establishment of specific interplays, such as hydrophobic and CH-π cloud interactions. By adopting three oligopeptides as model systems to assess folding features, we demonstrate herein that these finely tuned interactions dominate over electrostatic interactions, including H-bonds and electrostatic attractions/repulsions. The folding mechanism analysed here demonstrates cooperation at the single-residue level, for which we propose the terminology of "single residues cooperative folding". Overall, hydrophobic and CH-π cloud interactions produce the main output of the hydrophobic effect and govern the folding mechanism, as demonstrated in this study with small polypeptide chains, which in turn represent the main secondary structures in proteins.

The Temperature Dependence of Hydrogen Bonds Is More Uniform in Stable Proteins: An Analysis of NMR h3JNC' Couplings in Four Different Protein Structures.

Long-range HNCO NMR spectra for proteins show crosspeaks due to 1JNC', 2JNC', 3JNCγ, and h3JNC' couplings. The h3JNC' couplings are transmitted through hydrogen bonds and their sizes are correlated to hydrogen bond lengths. We collected long-range HNCO data at a series of temperatures for four protein structures. P22i and CUS-3i are six-stranded beta-barrel I-domains from phages P22 and CUS-3 that share less than 40% sequence identity. The cis and trans states of the C-terminal domain from pore-forming toxin hemolysin ΙΙ (HlyIIC) arise from the isomerization of a single G404-P405 peptide bond. For P22i and CUS-3i, hydrogen bonds detected by NMR agree with those observed in the corresponding domains from cryoEM structures of the two phages. Hydrogen bond lengths derived from the h3JNC' couplings, however, are poorly conserved between the distantly related CUS-3i and P22i domains and show differences even between the closely related cis and trans state structures of HlyIIC. This is consistent with hydrogen bond lengths being determined by local differences in structure rather than the overall folding topology. With increasing temperature, hydrogen bonds typically show an apparent increase in length that has been attributed to protein thermal expansion. Some hydrogen bonds are invariant with temperature, however, while others show apparent decreases in length, suggesting they become stabilized with increasing temperature. Considering the data for the three proteins in this study and previously published data for ubiquitin and GB3, lowered protein folding stability and cooperativity corresponds with a larger range of temperature responses for hydrogen bonds. This suggests a partial uncoupling of hydrogen bond energetics from global unfolding cooperativity as protein stability decreases.

On the role of native contact cooperativity in protein folding

The consistency of energy landscape theory predictions with available experimental data, as well as direct evidence from molecular simulations, have shown that protein folding mechanisms are largely determined by the contacts present in the native structure. As expected, native contacts are generally energetically favorable. However, there are usually at least as many energetically favorable nonnative pairs owing to the greater number of possible nonnative interactions. This apparent frustration must therefore be reduced by the greater cooperativity of native interactions. In this work, we analyze the statistics of contacts in the unbiased all-atom folding trajectories obtained by Shaw and coworkers, focusing on the unfolded state. By computing mutual cooperativities between contacts formed in the unfolded state, we show that native contacts form the most cooperative pairs, while cooperativities among nonnative or between native and nonnative contacts are typically much less favorable or even anticooperative. Furthermore, we show that the largest network of cooperative interactions observed in the unfolded state consists mainly of native contacts, suggesting that this set of mutually reinforcing interactions has evolved to stabilize the native state.

Tracing genetic diversity captures the molecular basis of misfolding disease

Genetic variation in human populations can result in the misfolding and aggregation of proteins, giving rise to systemic and neurodegenerative diseases that require management by proteostasis. Here, we define the role of GRP94, the endoplasmic reticulum Hsp90 chaperone paralog, in managing alpha-1-antitrypsin deficiency on a residue-by-residue basis using Gaussian process regression-based machine learning to profile the spatial covariance relationships that dictate protein folding arising from sequence variants in the population. Covariance analysis suggests a role for the ATPase activity of GRP94 in controlling the N- to C-terminal cooperative folding of alpha-1-antitrypsin responsible for the correction of liver aggregation and lung-disease phenotypes of alpha-1-antitrypsin deficiency. Gaussian process-based spatial covariance profiling provides a standard model built on covariant principles to evaluate the role of proteostasis components in guiding information flow from genome to proteome in response to genetic variation, potentially allowing us to intervene in the onset and progression of complex multi-system human diseases.

Assembly landscape for the bacterial large ribosomal subunit

Assembly of ribosomes in bacteria is highly efficient, taking ~2-3 min, but this makes the abundance of assembly intermediates very low, which is a challenge for mechanistic understanding. Genetic perturbations of the assembly process create bottlenecks where intermediates accumulate, facilitating structural characterization. We use cryo-electron microscopy, with iterative subclassification to identify intermediates in the assembly of the 50S ribosomal subunit from E. coli. The analysis of the ensemble of intermediates that spans the entire biogenesis pathway for the 50 S subunit was facilitated by a dimensionality reduction and cluster picking approach using PCA-UMAP-HDBSCAN. The identity of the cooperative folding units in the RNA with associated proteins is revealed, and the hierarchy of these units reveals a complete assembly map for all RNA and protein components. The assembly generally proceeds co-transcriptionally, with some flexibility in the landscape to ensure efficiency for this central cellular process under a variety of growth conditions.

Stability Islands and the Folding Cooperativity of a Seven-Repeat Array from Topoisomerase V.

Cooperativity is a central feature of protein folding, but the thermodynamic and structural origins of cooperativity remain poorly understood. To quantify cooperativity, we measured guanidine-induced unfolding transitions of single helix-hairpin-helix (HhH)2 repeats and tandem pairs from a seven-repeat segment of Methanopyrus kandleri Topoisomerase V (Topo V) to determine intrinsic repeat stability and interfacial free energies between repeats. Most single-repeat constructs are folded and stable; moreover, several pairs have unfolding midpoints that exceed midpoints of the single repeats they comprise, demonstrating favorable coupling between repeats. Analyzing unfolding transitions with a modified Ising model, we find a broad range of intrinsic and interfacial free energies. Surprisingly, the G repeat, which lacks density in the crystal structure of Topo V without DNA, is the most stable repeat in the array. Using nuclear magnetic resonance spectroscopy, we demonstrate that the isolated G repeat adopts a canonical (HhH)2 fold and forms an ordered interface with the F-repeat but not with the H repeat. Using parameters from our paired Ising fit, we built a partition function for the seven-repeat array. The multistate unfolding transition predicted from this partition function is in excellent agreement with the experimental unfolding transition, providing strong justification for the nearest-neighbor model. The seven-repeat partition function predicts a native state in which three independent segments ("stability islands") of interacting repeats are separated by two unstable interfaces. We confirm this segmented architecture by measuring the unfolding transition of an equimolar mixture of these three separate polypeptides. This segmented structural organization may facilitate wrapping around DNA.

Understanding the molecular basis of folding cooperativity through a comparative analysis of a multidomain protein and its isolated domains

Although cooperativity is a well-established and general property of folding, our current understanding of this feature in multidomain folding is still relatively limited. In fact, there are contrasting results indicating that the constituent domains of a multidomain protein may either fold independently on each other or exhibit interdependent supradomain phenomena. To address this issue, here we present the comparative analysis of the folding of a tandem repeat protein, comprising two contiguous PDZ domains, in comparison to that of its isolated constituent domains. By analyzing in detail the equilibrium and kinetics of folding at different experimental conditions, we demonstrate that despite each of the PDZ domains in isolation being capable of independent folding, at variance with previously characterized PDZ tandem repeats, the full-length construct folds and unfolds as a single cooperative unit. By exploiting quantitatively, the comparison of the folding of the tandem repeat to those observed for its constituent domains, as well as by characterizing a truncated variant lacking a short autoinhibitory segment, we successfully rationalize the molecular basis of the observed cooperativity and attempt to infer some general conclusions for multidomain systems.

GB1 hairpin kinetics: capturing the folding pathway with molecular dynamics, replica exchange and optimal dimensionality reduction

We have performed molecular dynamics (MD) and replica-exchange (REMD) simulations of folding of the GB1 hairpin peptide in aqueous solution. REMD results were consistent with a cooperative zipper folding model. 120 MD trajectories at 320 K yielded relaxation times of 1.8 and 100 ns, with the slower assigned to global folding. The MD folding/unfolding transitions also followed the cooperative zipper model, specifying nucleation at the central turn followed by consecutive hydrogen bond formation. Formation of hydrogen bonds and hydrophobic contacts were highly correlated. Coarse-grained kinetic models constructed with the Optimal Dimensionality Reduction (ODR) approach found a folding time of 3.3 and unfolding time of 4.0 Additionally, relaxation times in the 130–170 ns range could be assigned to formation of the transition state and off-path intermediates. The unfolded state was the most highly populated and, significantly, most heterogenous, assembling the largest number of microstates, primarily composed of extended and turn structures. The folded state was also heterogenous, but a to a lesser degree, involving the fully folded and partially folded in-register hairpins at early stages of the zipper pathway. The transition state corresponded to the nucleated hairpin, with central turn and first beta-sheet hydrogen bond, while the off-path intermediates were off-register partial hairpins. Our simulation results were in excellent agreement with experimental data on folded fraction, relaxation time and folding mechanism. The new findings from this work suggest a highly cooperative zipper folding mechanism, nascent hairpin transition state and ∼100 ns relaxation related to intermediate formation. Communicated by Ramaswamy H. Sarma

The cooperative folding of annexin A2 relies on a transient nonnative intermediate

Annexins (Anxs) are a family of highly homologous proteins that bind and aggregate lipid vesicles in the presence of calcium. All members of the family contain a variable N-terminus determining specific functions, followed by a conserved core region responsible for the general calcium-dependent lipid-binding property. The core structure consists of four homologous domains (DI–DIV), each consisting of a right-handed super-helix of five α-helices. We present data from a combination of site-directed mutagenesis, NMR, and circular dichroism showing that the G25–D34 region of the N-terminus as well as the contacts between residues D38A, R63A, and Q67A of AnxA2-DI are crucial for the autonomous folding and stability of DI of AnxA2. However, we also show that the folding of the full-length protein is very robust in that mutations and truncations that disrupted the folding of AnxA2-DI did not abolish the folding of full-length AnxA2, only lowering its thermal stability. This robustness of the folding of full-length AnxA2 is likely to be mediated by the existence of at least one transient nonnative intermediate as suggested by our kinetic data using stopped-flow fluorescence experiments. We also show that hydrophobic amino acids in AnxA2-DI involved in interfacial contacts with AnxA2-DIV are important for the cooperative folding and stability of the full-length protein. Mutating all of the V57E-V98R-G101Y residues in AnxA2-DI did not affect the folding of the domain, only its stability, but prevented the cooperative folding of the full-length protein. Our collective results favor a highly cooperative and robust folding process mediated by alternative intermediate steps. Since AnxA2 is a multifunctional protein involved in several steps of the progression of cell transformation, these data on structure and folding pathways are therefore crucial to designing anticancer drugs targeting AnxA2.

Loss of stability and unfolding cooperativity in hPGK1 upon gradual structural perturbation of its N-terminal domain hydrophobic core

Phosphoglycerate kinase has been a model for the stability, folding cooperativity and catalysis of a two-domain protein. The human isoform 1 (hPGK1) is associated with cancer development and rare genetic diseases that affect several of its features. To investigate how mutations affect hPGK1 folding landscape and interaction networks, we have introduced mutations at a buried site in the N-terminal domain (F25 mutants) that either created cavities (F25L, F25V, F25A), enhanced conformational entropy (F25G) or introduced structural strain (F25W) and evaluated their effects using biophysical experimental and theoretical methods. All F25 mutants folded well, but showed reduced unfolding cooperativity, kinetic stability and altered activation energetics according to the results from thermal and chemical denaturation analyses. These alterations correlated well with the structural perturbation caused by mutations in the N-terminal domain and the destabilization caused in the interdomain interface as revealed by H/D exchange under native conditions. Importantly, experimental and theoretical analyses showed that these effects are significant even when the perturbation is mild and local. Our approach will be useful to establish the molecular basis of hPGK1 genotype–phenotype correlations due to phosphorylation events and single amino acid substitutions associated with disease.

How Cooperatively Folding Are Homopolymer Molecular Knots?

Detailed thermodynamic analysis of complex systems with multiple stable configurational states allows for insight into the cooperativity of each individual transition. In this work we derive a heat capacity decomposition comprising contributions from each individual configurational state, which together sum to a baseline heat capacity, and contributions from each state-to-state transition. We apply this analysis framework to a series of replica exchange molecular dynamics simulations of linear and 1-1 coarse-grained homo-oligomer models which fold into stable, configurationally well-defined molecular knots, in order to better understand the parameters leading to stable and cooperative folding of these knots. We find that a stiff harmonic backbone bending angle potential is key to achieving knots with specific 3D structures. Tuning the backbone equilibrium angle in small increments yields a variety of knot topologies, including $3_1$, $5_1$, $7_1$, and $8_{19}$ types. Populations of different knotted states as functions of temperature can also be manipulated by tuning backbone torsion stiffness or by adding side chain beads. We find that sharp total heat capacity peaks for the homo-oligomer knots are largely due to a coil-to-globule transition, rather than a cooperative knotting step. However, in some cases the cooperativity of globule-to-knot and coil-to-globule transitions are comparable, suggesting that highly cooperative folding to knotted structures can be achieved by refining the model parameters or adding sequence specificity.

Biotinylation Eliminates the Intermediate State of Top7 Designed with an HIV-1 Epitope

Broadly neutralizing antibodies against HIV-1 are rare with the 2F5 antibody being one of the most protective. Insertion of an antibody epitope into a stable and small protein scaffold overcomes many of the obstacles found to produce antibodies. However, the design leads to grafting of epitopes that may cause protein aggregation. Here, I investigated the 2F5 epitope grafted into the Top7 as the scaffold in which the resulting immunoreactive protein precipitates along the storage time, as opposed to its completely soluble biotinylated version. Molecular dynamics showed that biotinylation eliminates the intermediate state of the scaffold-epitope Top7-2F5 by switching a noncooperative to a cooperative folding. The aggregation propensity of the Top7-designed proteins is examined in light of thermodynamic cooperativity and kinetic traps along the decreasing depth of the intermediate ensemble in the free energy landscape. This protocol may predict stable and soluble scaffold-epitopes with the purpose of composing novel therapeutic and diagnostic platforms.

Extension of the CHARMM Classical Drude Polarizable Force Field to N- and O-Linked Glycopeptides and Glycoproteins.

Molecular dynamic simulations are an effective tool to study complex molecular systems and are contingent upon the availability of an accurate and reliable molecular mechanics force field. The Drude polarizable force field, which allows for the explicit treatment of electronic polarization in a computationally efficient fashion, has been shown to reproduce experimental properties that were difficult or impossible to reproduce with the CHARMM additive force field, including peptide folding cooperativity, RNA hairpin structures, and DNA base flipping. Glycoproteins are essential components of glycoconjugate vaccines, antibodies, and many pharmaceutically important molecules, and an accurate polarizable force field that includes compatibility between the protein and carbohydrate aspect of the force field is essential to study these types of systems. In this work, we present an extension of the Drude polarizable force field to glycoproteins, including both N- and O-linked species. Parameter optimization focused on the dihedral terms using a reweighting protocol targeting NMR solution J-coupling data for model glycopeptides. Validation of the model include eight model glycopeptides and four glycoproteins with multiple N- and O-linked glycosylations. The new glycoprotein carbohydrate force field can be used in conjunction with the remainder of Drude polarizable force field through a variety of MD simulation programs including GROMACS, OPENMM, NAMD, and CHARMM and may be accessed through the Drude Prepper module in the CHARMM-GUI.

Conversion of the Native N-Terminal Domain of TDP-43 into a Monomeric Alternative Fold with Lower Aggregation Propensity.

TAR DNA-binding protein 43 (TDP-43) forms intraneuronal cytoplasmic inclusions associated with amyotrophic lateral sclerosis and ubiquitin-positive frontotemporal lobar degeneration. Its N-terminal domain (NTD) can dimerise/oligomerise with the head-to-tail arrangement, which is essential for function but also favours liquid-liquid phase separation and inclusion formation of full-length TDP-43. Using various biophysical approaches, we identified an alternative conformational state of NTD in the presence of Sulfobetaine 3-10 (SB3-10), with higher content of α-helical structure and tryptophan solvent exposure. NMR shows a highly mobile structure, with partially folded regions and β-sheet content decrease, with a concomitant increase of α-helical structure. It is monomeric and reverts to native oligomeric NTD upon SB3-10 dilution. The equilibrium GdnHCl-induced denaturation shows a cooperative folding and a somewhat lower conformational stability. When the aggregation processes were compared with and without pre-incubation with SB3-10, but at the identical final SB3-10 concentration, a slower aggregation was found in the former case, despite the reversible attainment of the native conformation in both cases. This was attributed to protein monomerization and oligomeric seeds disruption by the conditions promoting the alternative conformation. Overall, the results show a high plasticity of TDP-43 NTD and identify strategies to monomerise TDP-43 NTD for methodological and biomedical applications.

Studying protein folding landscape using NMR relaxation measurements with pressure perturbation

Repeat proteins greatly simplify the identification of the sequence determinants for protein folding cooperativity. We use the leucine rich repeat (LRR) domain of the tumor suppressor pp32 as a model. It is composed of five LRR with a capping domain on each of its termini. Previous pressure-dependent NMR experiments of cavity-containing mutants of pp32 revealed that the intermediates populated during pressure-induced unfolding depend upon the position of the cavity. In order to further characterize the nature of these intermediates, we carried out NMR relaxation experiments, relaxation dispersion (CPMG) and chemical exchange saturation transfer (CEST) experiments.