Exploring the Denatured State Ensemble by Single-Molecule Chemo-Mechanical Unfolding: The Effect of Force, Temperature, and Urea

Intermolecular interactions underlie protein/peptide phase separation irrespective of sequence and structure at crowded milieu

Liquid-liquid phase separation (LLPS) has emerged as a crucial biological phenomenon underlying the sequestration of macromolecules (such as proteins and nucleic acids) into membraneless organelles in cells. Unstructured and intrinsically disordered domains are known to facilitate multivalent interactions driving protein LLPS. We hypothesized that LLPS could be an intrinsic property of proteins/polypeptides but with distinct phase regimes irrespective of their sequence and structure. To examine this, we studied many (a total of 23) proteins/polypeptides with different structures and sequences for LLPS study in the presence and absence of molecular crowder, polyethylene glycol (PEG-8000). We showed that all proteins and even highly charged polypeptides (under study) can undergo liquid condensate formation, however with different phase regimes and intermolecular interactions. We further demonstrated that electrostatic, hydrophobic, and H-bonding or a combination of such intermolecular interactions plays a crucial role in individual protein/peptide LLPS.

Engineering shape memory and morphing protein hydrogels based on protein unfolding and folding

Engineering shape memory/morphing materials have achieved considerable progress in polymer-based systems with broad potential applications. However, engineering protein-based shape memory/morphing materials remains challenging and under-explored. Here we report the design of a bilayer protein-based shape memory/morphing hydrogel based on protein folding-unfolding mechanism. We fabricate the protein-bilayer structure using two tandem modular elastomeric proteins (GB1)₈ and (FL)₈. Both protein layers display distinct denaturant-dependent swelling profiles and Young’s moduli. Due to such protein unfolding-folding induced changes in swelling, the bilayer hydrogels display highly tunable and reversible bidirectional bending deformation depending upon the denaturant concentration and layer geometry. Based on these programmable and reversible bending behaviors, we further utilize the protein-bilayer structure as hinge to realize one-dimensional to two-dimensional and two-dimensional to three-dimensional folding transformations of patterned hydrogels. The present work will offer new inspirations for the design and fabrication of novel shape morphing materials.

Engineering shape memory and morphing materials achieved considerable progress in polymer-based systems, but protein-based shape memory and morphing materials remain less investigated. Here, the authors report the engineering of protein-based shape memory and morphing hydrogels using protein folding-unfolding as a general mechanism to trigger shape morphing in protein-bilayer structures.

Using Optical Tweezers to Monitor Allosteric Signals Through Changes in Folding Energy Landscapes

Signaling proteins are composed of conserved protein interaction domains that serve as allosteric regulatory elements of enzymatic or binding activities. The ubiquitous, structurally conserved cyclic nucleotide binding (CNB) domain is found covalently linked to proteins with diverse folds that perform multiple biological functions. Given that the structures of cAMP-bound CNB domains in different proteins are very similar, it remains a challenge to determine how this domain allosterically regulates such diverse protein functions and folds. Instead of a structural perspective, we focus our attention on the energy landscapes underlying the CNB domains and their responses to cAMP binding. We show that optical tweezers is an ideal tool to investigate how cAMP binding coupled to interdomain interactions remodel the energy landscape of the regulatory subunit of protein kinase A (PKA), which harbors two CNB domains. We mechanically manipulate and probe the unfolding and refolding behavior of the CNB domains as isolated structures or selectively as part of the PKA regulatory subunit, and in the presence and absence of cAMP. Optical tweezers allows us to dissect the changes in the energy landscape associated with cAMP binding, and to examine the allosteric interdomain interactions triggered by the cyclic nucleotide. This single molecule approach can be used to study other modular, multidomain signaling proteins found in nature.

From folding to function: complex macromolecular reactions unraveled one-by-one with optical tweezers

Single-molecule manipulation with optical tweezers has uncovered macromolecular behaviour hidden to other experimental techniques. Recent instrumental improvements have made it possible to expand the range of systems accessible to optical tweezers. Beyond focusing on the folding and structural changes of isolated single molecules, optical tweezers studies have evolved into unraveling the basic principles of complex molecular processes such as co-translational folding on the ribosome, kinase activation dynamics, ligand-receptor binding, chaperone-assisted protein folding, and even dynamics of intrinsically disordered proteins (IDPs). In this mini-review, we illustrate the methodological principles of optical tweezers before highlighting recent advances in studying complex protein conformational dynamics - from protein synthesis to physiological function - as well as emerging future issues that are beginning to be addressed with novel approaches.

Lessons about Protein Folding and Binding from Archetypal Folds

The function of proteins as biological nanomachines relies on their ability to fold into complex 3D structures, bind selectively to partners, and undergo conformational changes on cue. The native functional structures, and the rates of interconversion between conformational states (folded-unfolded, bound-free), are all encoded in the physical chemistry of their amino acid sequence. However, despite extensive research over decades, this code has proven difficult to fully crack, in terms of both prediction and understanding the molecular mechanisms at play.Earlier work on single-domain proteins reported a commonality of slow rates (10^-2-10² s^-1) and simple behavior in both kinetic and thermodynamic unfolding experiments, which suggested the process was all-or-none and thereby analogous to a chemical reaction (e.g., A ⇄ B). In the absence of a first-principles pre-exponential factor for protein (un)folding dynamics, the rates could only be interpreted in relative terms, e.g., the changes induced by mutation, and hence, neither the height of nor the entropic contribution to the free energy barriers was known. The rates were also many orders of magnitude too slow for direct atomistic simulations, and the computational focus was on predicting rate changes induced by mutation via coarse grained simulations. However, even the effects of mutation proved to be strikingly homogeneous with all experimental data clustering at ∼1/3 of the free energy perturbation recovered on folding and ∼2/3 on unfolding.The implementation of ultrafast kinetic methods turned the field upside down because they allowed researchers to measure the time scales of elementary (un)folding motions, which set the pre-exponential factor for protein conformational transitions at ∼1 μs. In parallel, we and others set out to investigate the simplest possible protein structures capable of autonomous folding, which we defined as archetypal folds. The rationale was to recapitulate the hierarchical organization of protein structure, starting from the bottom up. The study of fold archetypes ended up opening new research avenues in protein (un)folding, but also making unexpected connections with the folding upon binding of intrinsically disordered proteins and suggesting their functioning as conformational rheostats.This Account describes our work on the kinetic, thermodynamic, mechanistic, and functional analysis of fold archetypes. We first discuss the kinetic studies, emphasizing their impact on our understanding of (un)folding rates, of barrierless (downhill) folding, and as benchmarks for atomistic simulations. We continue with the thermodynamic analysis, introducing the differential scanning calorimetry, multiprobe, and NMR approaches that we developed to dissect their gradual, minimally cooperative (un)folding transitions and to probe the underlying mechanisms with unprecedented detail. The last two sections cover single-molecule analyses and some recent, mostly computational, results on the exploration of possible biological and technological roles for the gradual conformational transitions of fold archetypes.

Recent advances in understanding catalysis of protein folding by molecular chaperones

Molecular chaperones are highly conserved proteins that promote proper folding of other proteins in vivo. Diverse chaperone systems assist de novo protein folding and trafficking, the assembly of oligomeric complexes, and recovery from stress-induced unfolding. A fundamental function of molecular chaperones is to inhibit unproductive protein interactions by recognizing and protecting hydrophobic surfaces that are exposed during folding or following proteotoxic stress. Beyond this basic principle, it is now clear that chaperones can also actively and specifically accelerate folding reactions in an ATP-dependent manner. We focus on the bacterial Hsp70 and chaperonin systems as paradigms, and review recent work that has advanced our understanding of how these chaperones act as catalysts of protein folding.

On the nature of the optimal form of the holdase-type chaperone stress response

The holdase paradigm of chaperone action involves preferential binding by the chaperone to the unfolded protein state, thereby preventing it from either, associating with other unstable proteins (to form large dysfunctional aggregates), or being degraded by the proteolytic machinery of the cell/organism. In this paper, we examine the necessary physical constraints imposed upon the holdase chaperone response in a cell-like environment and use these limitations to comment on the likely nature of the optimal form of chaperone response in vivo.

Computational Estimation of Microsecond to Second Atomistic Folding Times

Despite the development of massively parallel computing hardware including inexpensive graphics processing units (GPUs), it has remained infeasible to simulate the folding of atomistic proteins at room temperature using conventional molecular dynamics (MD) beyond the microsecond scale. Here, we report the folding of atomistic, implicitly solvated protein systems with folding times τ ranging from ∼10 μs to ∼100 ms using the weighted ensemble (WE) strategy in combination with GPU computing. Starting from an initial structure or set of structures, WE organizes an ensemble of GPU-accelerated MD trajectory segments via intermittent pruning and replication events to generate statistically unbiased estimates of rate constants for rare events such as folding; no biasing forces are used. Although the variance among atomistic WE folding runs is significant, multiple independent runs are used to reduce and quantify statistical uncertainty. Folding times are estimated directly from WE probability flux and from history-augmented Markov analysis of the WE data. Three systems were examined: NTL9 at low solvent viscosity (yielding τf = 0.8-9 μs), NTL9 at water-like viscosity (τf = 0.2-2 ms), and Protein G at low viscosity (τf = 3-200 ms). In all cases, the folding time, uncertainty, and ensemble properties could be estimated from WE simulation; for Protein G, this characterization required significantly less overall computing than would be required to observe a single folding event with conventional MD simulations. Our results suggest that the use and calibration of force fields and solvent models for precise estimation of kinetic quantities is becoming feasible.

Using Single-Molecule Chemo-Mechanical Unfolding to Simultaneously Probe Multiple Structural Parameters in Protein Folding

While single-molecule force spectroscopy has greatly advanced the study of protein folding, there are limitations to what can be learned from studying the effect of force alone. We developed a novel technique, chemo-mechanical unfolding, that combines multiple perturbants—force and chemical denaturant—to more fully characterize the folding process by simultaneously probing multiple structural parameters—the change in end-to-end distance, and solvent accessible surface area. Here, we describe the theoretical background, experimental design, and data analysis for chemo-mechanical unfolding experiments probing protein folding thermodynamics and kinetics. This technique has been applied to characterize parallel protein folding pathways, the protein denatured state, protein folding on the ribosome, and protein folding intermediates.

A small single-domain protein folds through the same pathway on and off the ribosome

In vivo, proteins fold and function in a complex environment subject to many stresses that can modulate a protein's energy landscape. One aspect of the environment pertinent to protein folding is the ribosome, since proteins have the opportunity to fold while still bound to the ribosome during translation. We use a combination of force and chemical denaturant (chemomechanical unfolding), as well as point mutations, to characterize the folding mechanism of the src SH3 domain both as a stalled ribosome nascent chain and free in solution. Our results indicate that src SH3 folds through the same pathway on and off the ribosome. Molecular simulations also indicate that the ribosome does not affect the folding pathway for this small protein. Taken together, we conclude that the ribosome does not alter the folding mechanism of this small protein. These results, if general, suggest the ribosome may exert a bigger influence on the folding of multidomain proteins or protein domains that can partially fold before the entire domain sequence is outside the ribosome exit tunnel.

Local and non-local topological information in the denatured state ensemble of a β-barrel protein

The folding of predominantly β-sheet proteins is complicated by the presence of a large number of non-local interactions in their native states, which increase the ruggedness of their folding energy landscapes. However, forming non-local contacts early in folding or even in the unfolded state can smooth the energy landscape and facilitate productive folding. We report that several sequence regions of a β-barrel protein, cellular retinoic acid-binding protein 1 (CRABP1), populate native-like secondary structure to a significant extent in the denatured state in 8 M urea. In addition, we provide evidence for both local and non-local interactions in the denatured state of CRABP1. NMR chemical shift perturbations (CSPs) under denaturing conditions upon substitution of single residues by mutation support the presence of several non-local interactions in topologically key sites, arguing that the denatured state is conformationally restricted and contains topological information for the native fold. Among the most striking non-local interactions are those between the N- and C-terminal regions, which are involved in closure of the native β-barrel. In addition, CSPs support the presence of two features in the denatured state: a major hydrophobic cluster involving residues from various parts of the sequence and a native-like interaction similar to one identified in previous studies as forming early in folding (Budyak et al., Structure 21, 476 [2013]). Taken together, our data support a model in which transient structures involving nonlocal interactions prime early folding interactions in CRABP1, determine its barrel topology, and may protect this predominantly β-sheet protein against aggregation.

Modulating charge patterning and ionic strength as a strategy to induce conformational changes in intrinsically disordered proteins

We present an analytical theory to describe conformational changes as a function of salt for polymers with a given sequence of charges. We apply this model to describe Intrinsically Disordered Proteins (IDPs) by explicitly accounting for charged residues and their exact placement in the primary sequence while approximating the effect of non-electrostatic interactions at a mean-field level by effective short-range (two body and three-body) interaction parameters. The effect of ions is introduced by treating electrostatic interactions within Debye-Huckle approximation. Using typical values of the short-range mean-field parameters derived from all-atom Monte Carlo simulations (at zero salt), we predict the conformational changes as a function of salt concentration. We notice that conformational transitions in response to changes in ionic strength strongly depend on sequence specific charge patterning. For example, globule to coil transition can be observed upon increasing salt concentration, in stark contrast to uniformly charged polyelectrolyte theories based on net charge only. In addition, it is possible to observe non-monotonic behavior with salt as well. Drastic differences in salt-induced conformational transitions is also evident between two doubly phosphorylated sequences-derived from the same wild type sequence-that only differ in the site of phosphorylation. Similar effects are also predicted between two sequences derived from the same parent sequence differing by a single site mutation where a negative charge is replaced by a positive charge. These effects are purely a result of charge decoration and can only be understood in terms of metrics based on specific placement of charges, and cannot be explained by models based on charge composition alone. Identifying sequences and hot spots within sequences-for post translational modification or charge mutation-using our high-throughput theory will yield fundamental insights into design and biological regulation mediated by phosphorylation and/or local changes in salt concentration.