Folding Pathway Research Articles

Elucidation of the folding pathway of a circular permutant of topologically knotted YbeA by tryptophan substitutions

CP74 is an engineered circular permutant of a deep trefoil knotted SpoU-TrmD (SPOUT) RNA methyl transferase protein YbeA from E. coli. We have previously established that the circular permutation unties the knotted topology of YbeA and CP74 forms a domain-swapped dimer with a large dimeric interface of ca. 4600 Å2. To understand the impact of domain-swapping and the newly formed hinge region joining the two folded domains on the folding and stability of CP74, the five equally spaced tryptophan residues were individually substituted into phenylalanine to monitor their conformational and stability changes by a battery of biophysical tools. Far-UV circular dichroism, intrinsic fluorescence, and small-angle X-ray scattering dictated minimal global conformational perturbations to the native structures in the tryptophan variants. The structures of the tryptophan variants also showed the conservation of the domain-swapped ternary structure with the exception that the W72F exhibited significant asymmetry in the α-helix 5. Comparative global thermal and chemical stability analyses indicated the pivotal role of W100 in the folding of CP74 followed by W19 and W72. Solution-state NMR spectroscopy and hydrogen-deuterium exchange mass spectrometry further revealed the accumulation of a native-like intermediate state in which the hinge region made important contributions to maintain the domain-swapped ternary structure of CP74.

Wave space sonification of the folding pathways of protein molecules modeled as hyper-redundant robotic mechanisms

Wave space sonification of the folding pathways of protein molecules modeled as hyper-redundant robotic mechanisms

Robust membrane protein tweezers reveal the folding speed limit of helical membrane proteins.

Single-molecule tweezers, such as magnetic tweezers, are powerful tools for probing nm-scale structural changes in single membrane proteins under force. However, the weak molecular tethers used for the membrane protein studies have limited the observation of long-time, repetitive molecular transitions due to force-induced bond breakage. The prolonged observation of numerous transitions is critical in reliable characterizations of structural states, kinetics, and energy barrier properties. Here, we present a robust single-molecule tweezer method that uses dibenzocyclooctyne cycloaddition and traptavidin binding, enabling the estimation of the folding 'speed limit' of helical membrane proteins. This method is >100 times more stable than a conventional linkage system regarding the lifetime, allowing for the survival for ~12 hr at 50 pN and ~1000 pulling cycle experiments. By using this method, we were able to observe numerous structural transitions of a designer single-chained transmembrane homodimer for 9 hr at 12 pN and reveal its folding pathway including the hidden dynamics of helix-coil transitions. We characterized the energy barrier heights and folding times for the transitions using a model-independent deconvolution method and the hidden Markov modeling analysis, respectively. The Kramers rate framework yields a considerably low-speed limit of 21 ms for a helical hairpin formation in lipid bilayers, compared to μs scale for soluble protein folding. This large discrepancy is likely due to the highly viscous nature of lipid membranes, retarding the helix-helix interactions. Our results offer a more valid guideline for relating the kinetics and free energies of membrane protein folding.

Thermodynamic and kinetic analysis of the LAO binding protein and its isolated domains reveal non-additivity in stability, folding and function.

Substrate-binding proteins (SBPs) are used by organisms from the three domains of life for transport and signalling. SBPs are composed of two domains that collectively trap ligands with high affinity and selectivity. To explore the role of the domains and the integrity of the hinge region between them in the function and conformation of SBPs, here, we describe the ligand binding, conformational stability and folding kinetics of the Lysine Arginine Ornithine (LAO) binding protein from Salmonella thiphimurium and constructs corresponding to its two independent domains. LAO is a class II SBP formed by a continuous and a discontinuous domain. Contrary to the expected behaviour based on their connectivity, the discontinuous domain shows a stable native-like structure that binds l-arginine with moderate affinity, whereas the continuous domain is barely stable and shows no detectable ligand binding. Regarding folding kinetics, studies of the entire protein revealed the presence of at least two intermediates. While the unfolding and refolding of the continuous domain exhibited only a single intermediate and simpler and faster kinetics than LAO, the folding mechanism of the discontinuous domain was complex and involved multiple intermediates. These findings suggest that in the complete protein the continuous domain nucleates folding and that its presence funnels the folding of the discontinuous domain avoiding nonproductive interactions. The strong dependence of the function, stability and folding pathway of the lobes on their covalent association is most likely the result of the coevolution of both domains as a single unit.

Hsp90 mutants with distinct defects provide novel insights into cochaperone regulation of the folding cycle.

Molecular chaperones play a key role in maintaining proteostasis and cellular health. The abundant, essential, cytosolic Hsp90 (Heat shock protein, 90 kDa) facilitates the folding and activation of hundreds of newly synthesized or misfolded client proteins in an ATP-dependent folding pathway. In a simplified model, Hsp70 first helps load client onto Hsp90, ATP binding results in conformational changes in Hsp90 that result in the closed complex, and then less defined events result in nucleotide hydrolysis, client release and return to the open state. Cochaperones bind and assist Hsp90 during this process. We previously identified a series of yeast Hsp90 mutants that appear to disrupt either the 'loading', 'closing' or 'reopening' events, and showed that the mutants had differing effects on activity of some clients. Here we used those mutants to dissect Hsp90 and cochaperone interactions. Overexpression or deletion of HCH1 had dramatically opposing effects on the growth of cells expressing different mutants, with a phenotypic shift coinciding with formation of the closed conformation. Hch1 appears to destabilize Hsp90-nucleotide interaction, hindering formation of the closed conformation, whereas Cpr6 counters the effects of Hch1 by stabilizing the closed conformation. Hch1 and the homologous Aha1 share some functions, but the role of Hch1 in inhibiting progression through the early stages of the folding cycle is unique. Sensitivity to the Hsp90 inhibitor NVP-AUY922 also correlates with the conformational cycle, with mutants defective in the loading phase being most sensitive and those defective in the reopening phase being most resistant to the drug. Overall, our results indicate that the timing of transition into and out of the closed conformation is tightly regulated by cochaperones. Further analysis will help elucidate additional steps required for progression through the Hsp90 folding cycle and may lead to new strategies for modulating Hsp90 function.

Deep Boosted Molecular Dynamics: Accelerating Molecular Simulations with Gaussian Boost Potentials Generated Using Probabilistic Bayesian Deep Neural Network.

We have developed a new deep boosted molecular dynamics (DBMD) method. Probabilistic Bayesian neural network models were implemented to construct boost potentials that exhibit Gaussian distribution with minimized anharmonicity, thereby allowing for accurate energetic reweighting and enhanced sampling of molecular simulations. DBMD was demonstrated on model systems of alanine dipeptide and the fast-folding protein and RNA structures. For alanine dipeptide, 30 ns DBMD simulations captured up to 83-125 times more backbone dihedral transitions than 1 μs conventional molecular dynamics (cMD) simulations and were able to accurately reproduce the original free energy profiles. Moreover, DBMD sampled multiple folding and unfolding events within 300 ns simulations of the chignolin model protein and identified low-energy conformational states comparable to previous simulation findings. Finally, DBMD captured a general folding pathway of three hairpin RNAs with the GCAA, GAAA, and UUCG tetraloops. Based on a deep learning neural network, DBMD provides a powerful and generally applicable approach to boosting biomolecular simulations. DBMD is available with open source in OpenMM at https://github.com/MiaoLab20/DBMD/.

Is Posttranslational Folding More Efficient Than Refolding from a Denatured State: A Computational Study.

The folding of proteins into their native conformation is a complex process that has been extensively studied over the past half-century. The ribosome, the molecular machine responsible for protein synthesis, is known to interact with nascent proteins, adding further complexity to the protein folding landscape. Consequently, it is unclear whether the folding pathways of proteins are conserved on and off the ribosome. The main question remains: to what extent does the ribosome help proteins fold? To address this question, we used coarse-grained molecular dynamics simulations to compare the mechanisms by which the proteins dihydrofolate reductase, type III chloramphenicol acetyltransferase, and d-alanine-d-alanine ligase B fold during and after vectorial synthesis on the ribosome to folding from the full-length unfolded state in bulk solution. Our results reveal that the influence of the ribosome on protein folding mechanisms varies depending on the size and complexity of the protein. Specifically, for a small protein with a simple fold, the ribosome facilitates efficient folding by helping the nascent protein avoid misfolded conformations. However, for larger and more complex proteins, the ribosome does not promote folding and may contribute to the formation of intermediate misfolded states cotranslationally. These misfolded states persist posttranslationally and do not convert to the native state during the 6 μs runtime of our coarse-grain simulations. Overall, our study highlights the complex interplay between the ribosome and protein folding and provides insight into the mechanisms of protein folding on and off the ribosome.

Pathway and mechanism of tubulin folding mediated by TRiC/CCT along its ATPase cycle revealed using cryo-EM

The eukaryotic chaperonin TRiC/CCT assists the folding of about 10% of cytosolic proteins through an ATP-driven conformational cycle, and the essential cytoskeleton protein tubulin is the obligate substrate of TRiC. Here, we present an ensemble of cryo-EM structures of endogenous human TRiC throughout its ATPase cycle, with three of them revealing endogenously engaged tubulin in different folding stages. The open-state TRiC-tubulin-S1 and -S2 maps show extra density corresponding to tubulin in the cis-ring chamber of TRiC. Our structural and XL-MS analyses suggest a gradual upward translocation and stabilization of tubulin within the TRiC chamber accompanying TRiC ring closure. In the closed TRiC-tubulin-S3 map, we capture a near-natively folded tubulin—with the tubulin engaging through its N and C domains mainly with the A and I domains of the CCT3/6/8 subunits through electrostatic and hydrophilic interactions. Moreover, we also show the potential role of TRiC C-terminal tails in substrate stabilization and folding. Our study delineates the pathway and molecular mechanism of TRiC-mediated folding of tubulin along the ATPase cycle of TRiC, and may also inform the design of therapeutic agents targeting TRiC-tubulin interactions.

Non-Native Structures of Apomyoglobin and Apoleghemoglobin in Folding Intermediates Related to the Protein Misfolding.

Protein folding is essential for a polypeptide chain to acquire its proper structure and function. Globins are a superfamily of ubiquitous heme-binding α-helical proteins whose function is principally to regulate oxygen homoeostasis. In this review, we explore the hierarchical helical formation in the globin proteins apomyoglobin and leghemoglobin, and we discuss the existence of non-native and misfolded structures occurring during the course of folding to its native state. This review summarizes the research aimed at characterizing and comparing the equilibrium and kinetic intermediates, as well as delineating the complete folding pathway at a molecular level, in order to answer the following questions: "What is the mechanism of misfolding via a folding intermediate? Does the non-native structure stabilize the contemporary intermediate structure? Does the non-native structure induce slower folding?" The role of the non-native structures in the folding intermediate related to misfolding is also discussed.

A Novel Peptide Reagent for Investigating Disulfide-Coupled Folding Intermediates of Mid-Size Proteins.

Investigations of protein folding have largely involved the use of disulfide-containing proteins, since the disulfide-coupled folding of proteins allows folding intermediates to be trapped and their conformations determined. However, studies of the folding mechanisms of mid-size proteins face several problems, one of which is that detecting folding intermediates is difficult. Therefore, to solve this issue, a novel peptide reagent, maleimidohexanoyl-Arg5-Tyr-NH2, was designed and applied to the detection of folding intermediates of model proteins. BPTI was chosen as a model small protein to estimate the ability of the novel reagent to detect folding intermediates. In addition, a precursor protein (prococoonase) of Bombyx mori cocoonase was used as a model mid-size protein. Cocoonase is classified as a serine protease and has a high homology with trypsin. We recently found that the propeptide sequence of prococoonase (proCCN) is important for the folding of cocoonase. However, it was difficult to study the folding pathway of proCCN since the folding intermediates could not be separated on a reversed-phase HPLC (RP-HPLC). Therefore, to separate the folding intermediates by RP-HPLC, the novel labeling reagent was used to accomplish this for proCCN. The results indicated that the peptide reagent allowed the intermediates to be captured, separated on SDS-PAGE, and analyzed by RP-HPLC without the occurrence of undesirable disulfide-exchange reactions during the labeling reactions. The peptide reagent reported herein is a practical tool for investigating the mechanisms of disulfide-coupled folding of mid-size proteins.

N-glycosylation mediated folding and quality control in serine proteases of the hepsin family.

N-linked glycans are specifically attached to asparagine residues in a N-X-S/T motif of secretory pathway glycoproteins. N-glycosylation of newly synthesized glycoproteins directs their folding via the lectin chaperones calnexin and calreticulin that are associated with protein-folding enzymes and glycosidases of the endoplasmic reticulum (ER). Misfolded glycoproteins are retained in the ER by the same lectin chaperones. The work by Sun etal. (FEBS J 2023, 10.1111/febs.16757) in this issue focusses on hepsin, a serine protease on the surface of liver and other organs. The authors deduce that spatial positioning of N-glycans on one side of a conserved domain of hepsin, known as the scavenger receptor-rich cysteine domain, regulates calnexin selection for hepsin maturation and transport through the secretory pathway. If N-glycosylation is elsewhere on hepsin, then it is misfolded and has a prolonged accumulation with calnexin and BiP. This association coincides with the engagement of stress response pathways thatsense glycoprotein misfolding. The topological considerations of N-glycosylation dissected by Sun etal. may help unravel how key sites of N-glycosylation sites required for protein folding and transport have evolved to select the lectin chaperone calnexin pathway for folding and quality control.

Generating Protein Folding Trajectories Using Contact-Map-Driven Directed Walks.

Recent advances in machine learning methods have had a significant impact on protein structure prediction, but accurate generation and characterization of protein-folding pathways remains intractable. Here, we demonstrate how protein folding trajectories can be generated using a directed walk strategy operating in the space defined by the residue-level contact-map. This double-ended strategy views protein folding as a series of discrete transitions between connected minima on the potential energy surface. Subsequent reaction-path analysis for each transition enables thermodynamic and kinetic characterization of each protein-folding path. We validate the protein-folding paths generated by our discretized-walk strategy against direct molecular dynamics simulations for a series of model coarse-grained proteins constructed from hydrophobic and polar residues. This comparison demonstrates that ranking discretized paths based on the intermediate energy barriers provides a convenient route to identifying physically sensible folding ensembles. Importantly, by using directed walks in the protein contact-map space, we circumvent several of the traditional challenges associated with protein-folding studies, namely, long time scales required and the choice of a specific order parameter to drive the folding process. As such, our approach offers a useful new route for studying the protein-folding problem.

Impact of the yeast S0/uS2-cluster ribosomal protein rpS21/eS21 on rRNA folding and the architecture of small ribosomal subunit precursors.

RpS0/uS2, rpS2/uS5, and rpS21/eS21 form a cluster of ribosomal proteins (S0-cluster) at the head-body junction near the central pseudoknot of eukaryotic small ribosomal subunits (SSU). Previous work in yeast indicated that S0-cluster assembly is required for the stabilisation and maturation of SSU precursors at specific post-nucleolar stages. Here, we analysed the role of S0-cluster formation for rRNA folding. Structures of SSU precursors isolated from yeast S0-cluster expression mutants or control strains were analysed by cryogenic electron microscopy. The obtained resolution was sufficient to detect individual 2'-O-methyl RNA modifications using an unbiased scoring approach. The data show how S0-cluster formation enables the initial recruitment of the pre-rRNA processing factor Nob1 in yeast. Furthermore, they reveal hierarchical effects on the pre-rRNA folding pathway, including the final maturation of the central pseudoknot. Based on these structural insights we discuss how formation of the S0-cluster determines at this early cytoplasmic assembly checkpoint if SSU precursors further mature or are degraded.

Protein structure and folding pathway prediction based on remote homologs recognition using PAthreader

Recognition of remote homologous structures is a necessary module in AlphaFold2 and is also essential for the exploration of protein folding pathways. Here, we propose a method, PAthreader, to recognize remote templates and explore folding pathways. Firstly, we design a three-track alignment between predicted distance profiles and structure profiles extracted from PDB and AlphaFold DB, to improve the recognition accuracy of remote templates. Secondly, we improve the performance of AlphaFold2 using the templates identified by PAthreader. Thirdly, we explore protein folding pathways based on our conjecture that dynamic folding information of protein is implicitly contained in its remote homologs. The results show that the average accuracy of PAthreader templates is 11.6% higher than that of HHsearch. In terms of structure modelling, PAthreader outperform AlphaFold2 and ranks first on the CAMEO blind test for the latest three months. Furthermore, we predict protein folding pathways for 37 proteins, in which the results of 7 proteins are almost consistent with those of biological experiments, and the other 30 human proteins have yet to be verified by biological experiments, revealing that folding information can be exploited from remote homologous structures.

Functional, pathogenic, and pharmacological roles of protein folding intermediates.

Protein expression and function in eukaryotic cells are tightly harmonized processes modulated by the combination of different layers of regulation, including transcription, processing, stability, and translation of messenger RNA, as well as assembly, maturation, sorting, recycling, and degradation of polypeptides. Integrating all these pathways and the protein quality control machinery, deputed to avoid the production and accumulation of aberrantly folded proteins, determines protein homeostasis. Over the last decade, the combined development of accurate time-resolved experimental techniques and efficient computer simulations has opened the possibility of investigating biological mechanisms at atomic resolution with physics-based models. A meaningful example is the reconstruction of protein folding pathways at atomic resolution, which has enabled the characterization of the folding kinetics of biologically relevant globular proteins consisting of a few hundred amino acids. Combining these innovative computational technologies with rigorous experimental approaches reveals the existence of non-native metastable states transiently appearing along the folding process of such proteins. Here, we review the primary evidence indicating that these protein folding intermediates could play roles in disparate biological processes, from the posttranslational regulation of protein expression to disease-relevant protein misfolding mechanisms. Finally, we discuss how the information encoded into protein folding pathways could be exploited to design an entirely new generation of pharmacological agents capable of promoting the selective degradation of protein targets.

Investigating the role of computer simulations in modeling complex biological systems

The use of computer simulations in modeling complex biological systems has become an important tool for biotechnology research. This study investigates the role of computer simulations in modeling complex biological systems and their potential applications in biotechnology. The focus of the study was on protein folding, a complex biological process that plays a critical role in the structure and function of proteins. The simulation was carried out using molecular dynamics techniques that allowed the simulation of the movement of the protein molecules over time. The simulation results were compared to experimental data and showed a good agreement with the actual protein structure. Furthermore, the simulation provided insights into the dynamics of the protein folding process, such as the folding pathway and the stability of the intermediate states. These insights can be used to design new drugs and therapies that target specific protein structures and functions. The study concludes that computer simulations are a valuable tool in biotechnology research and have the potential to revolutionize the field by providing accurate and detailed information about complex biological systems. Future research could investigate the use of advanced computer simulations and machine learning algorithms in biotechnology applications to improve their accuracy and efficiency.

Inferring Pathways of Oxidative Folding from Prefolding Free Energy Landscapes of Disulfide-Rich Toxins.

Short, cysteine-rich peptides can exist in stable or metastable structural ensembles due to the number of possible patterns of formation of their disulfide bonds. One interesting subset of this peptide group is the conotoxins, which are produced by aquatic snails in the family Conidae. The μ conotoxins, which are antagonists and blockers of the voltage-gated sodium channel, exist in a folding spectrum: on one end of the spectrum are more hirudin-like folders, which form disulfide bonds and then reshuffle them, leading to an ensemble of kinetically trapped isomers, and on the other end are more BPTI-like folders, which form the native disulfide bonds one by one in a particular order, leading to a preponderance of conformations existing in a single stable state. In this Article, we employ the composite diffusion map approach to study the unified free energy surface of prefolding μ-conotoxin equilibrium. We identify the two most important nonlinear collective modes of the unified folding landscape and demonstrate that in the absence of their disulfides, the conotoxins can be thought of as largely disordered polymers. A small increase in the number of hydrophobic residues in the protein shifts the free energy landscape toward hydrophobically collapsed coil conformations responsible for cysteine proximity in hirudin-like folders, compared to semiextended coil conformations with more distal cysteines in BPTI-like folders. Overall, this work sheds important light on the folding processes and free energy landscapes of cysteine-rich peptides and demonstrates the extent to which sequence and length contribute to these landscapes.

Biophysical Correlates of Enhanced Immunogenicity of a Stabilized Variant of the Receptor Binding Domain of SARS-CoV-2.

The receptor binding domain (RBD) of SARS-CoV-2 is the primary target of neutralizing antibodies. We have previously reported the design and characterization of a mammalian cell expressed RBD derivative, mRBD1-3.2, that has higher thermal stability and greatly enhanced immunogenicity relative to the wild type mRBD. The protein is highly thermotolerant and immunogenic and is being explored for use in room temperature stable Covid-19 vaccine formulations. In the current study, we have investigated the folding pathway of both WT and stabilized RBD. It was found that chemical denaturation of RBD proceeds through a stable equilibrium intermediate. Thermal and chemical denaturation is reversible, as assayed by binding to the receptor ACE2. Unusually, in its native state, RBD binds to the hydrophobic probe ANS, and enhanced ANS binding is observed for the equilibrium intermediate state. Further characterization of the folding of mRBD1-3.2, both in solution and after reconstitution of lyophilized protein stored for a month at 37 °C, revealed a higher stability represented by higher Cm, faster refolding, slower unfolding, and enhanced resistance to proteolytic cleavage relative to WT. In contrast to WT RBD, the mutant showed decreased interaction with the hydrophobic moiety linoleic acid. Collectively, these data suggest that the enhanced immunogenicity results from reduced conformational fluctuations that likely enhance in vivo half-life as well as reduce the exposure of irrelevant non-neutralizing epitopes to the immune system.

Learning to self-fold at a bifurcation.

Disordered mechanical systems can deform along a network of pathways that branch and recombine at special configurations called bifurcation points. Multiple pathways are accessible from these bifurcation points; consequently, computer-aided design algorithms have been sought to achieve a specific structure of pathways at bifurcations by rationally designing the geometry and material properties of these systems. Here, we explore an alternative physical training framework in which the topology of folding pathways in a disordered sheet is changed in a desired manner due to changes in crease stiffnesses induced by prior folding. We study the quality and robustness of such training for different "learning rules," that is, different quantitative ways in which local strain changes the local folding stiffness. We experimentally demonstrate these ideas using sheets with epoxy-filled creases whose stiffnesses change due to folding before the epoxy sets. Our work shows how specific forms of plasticity in materials enable them to learn nonlinear behaviors through their prior deformation history in a robust manner.

Protein folding pathway and intermediates on the dimensionality reduced landscape with UMAP and native contact likelihood.

To understand protein folding mechanisms from molecular dynamics (MD) simulations, it is important to explore not only folded/unfolded states but also representative intermediate structures on the conformational landscape. Generally, such a conformational landscape is obtained by linear dimensionality reduction methods, principal component analysis (PCA) or time-lagged independent component analysis (tICA), for instance. However, it is often difficult to apply these linear methods to large-scale, nonlinear conformational changes such as protein folding. Here, we propose a novel approach to construct the landscape using the uniform manifold approximation and projection (UMAP) method, which nonlinearly reduces the dimensionality without losing data-point proximity. In the approach, native contact likelihood, a measure for native contact formation, is used as feature variables rather than the conventional Cartesian coordinates or dihedral angles of protein structures. We tested the performance of UMAP for coarse-grained MD simulation trajectories of B1 domain in protein G and observed on-pathway transient structures and other metastable states on the UMAP conformational landscape. In contrast, these structures were not clearly distinguished on the dimensionality reduced landscape using PCA or tICA. This approach is also useful to obtain dynamical information through Markov State Modeling and would be applicable to large-scale conformational changes in many other biomacromolecules.