How cooperative are protein folding and unfolding transitions?

Non-Glycosylated SARS-CoV-2 Omicron BA.5 Receptor Binding Domain (RBD) with a Native-like Conformation Induces a Robust Immune Response with Potent Neutralization in a Mouse Model

The Omicron BA.5 variant of SARS-CoV-2 is known for its high transmissibility and its capacity to evade immunity provided by vaccine protection against the (original) Wuhan strain. In our prior research, we successfully produced the receptor-binding domain (RBD) of the SARS-CoV-2 spike protein in an E. coli expression system. Extensive biophysical characterization indicated that, even without glycosylation, the RBD maintained native-like conformational and biophysical properties. The current study explores the immunogenicity and neutralization capacity of the E. coli-expressed Omicron BA.5 RBD using a mouse model. Administration of three doses of the RBD without any adjuvant elicited high titer antisera of up to 7.3 × 105 and up to 1.6 × 106 after a booster shot. Immunization with RBD notably enhanced the population of CD44+CD62L+ T cells, indicating the generation of T cell memory. The in vitro assays demonstrated the antisera's protective efficacy through significant inhibition of the interaction between SARS-CoV-2 and its human receptor, ACE2, and through potent neutralization of a pseudovirus. These findings underscore the potential of our E. coli-expressed RBD as a viable vaccine candidate against the Omicron variant of SARS-CoV-2.

Spectroscopic analysis of the bacterially expressed head domain of rotavirus VP6

The rotavirus capsid protein VP6 forms the middle of three protein layers and is responsible for many critical steps in the viral life cycle. VP6 as a structural protein can be used in various applications including as a subunit vaccine component. The head domain of VP6 (VP6H) contains key sequences that allow the protein to trimerize and that represent epitopes that are recognized by human antibodies in the viral particle. The domain is rich in β-sheet secondary structures. Here, VP6H was solubilised from bacterial inclusion bodies and purified using a single affinity chromatography step. Spectral (far-UV circular dichroism and intrinsic tryptophan fluorescence) analysis revealed that the purified domain had native-like secondary and tertiary structures. The domain could maintain structure up to 44°C during thermal denaturation following which structural changes result in an intermediate forming and finally irreversible aggregation and denaturation. The chemical denaturation with urea and guanidinium hydrochloride produces intermediates that represent a loss in the cooperativity. The VP6H domain is stable and can fold to produce its native structure in the absence of the VP6 base domain but cannot be defined as an independent folding unit.

Protein G-quadruplex interactions and their effects on phase transitions and protein aggregation

The SERF family of proteins were originally discovered for their ability to accelerate amyloid formation. Znf706 is an uncharacterized protein whose N-terminus is homologous to SERF proteins. We show here that human Znf706 can promote protein aggregation and amyloid formation. Unexpectedly, Znf706 specifically interacts with stable, non-canonical nucleic acid structures known as G-quadruplexes. G-quadruplexes can affect gene regulation and suppress protein aggregation; however, it is unknown if and how these two activities are linked. We find Znf706 binds preferentially to parallel G-quadruplexes with low micromolar affinity, primarily using its N-terminus, and upon interaction, its dynamics are constrained. G-quadruplex binding suppresses Znf706's ability to promote protein aggregation. Znf706 in conjunction with G-quadruplexes therefore may play a role in regulating protein folding. RNAseq analysis shows that Znf706 depletion specifically impacts the mRNA abundance of genes that are predicted to contain high G-quadruplex density. Our studies give insight into how proteins and G-quadruplexes interact, and how these interactions affect both partners and lead to the modulation of protein aggregation and cellular mRNA levels. These observations suggest that the SERF family of proteins, in conjunction with G-quadruplexes, may have a broader role in regulating protein folding and gene expression than previously appreciated.

A Model Study to Assess Fibrillation and Product Stability to Support Peptide Drug Design

The fibrillation of therapeutic peptides can present significant quality concerns and poses challenges for manufacturing and storage. A fundamental understanding of the mechanisms of fibrillation is critical for the rational design of fibrillation-resistant peptide drugs and can accelerate product development by guiding the selection of solution-stable candidates and formulations. The studies reported here investigated the effects of structural modifications on the fibrillation of a 29-residue peptide (PepA) and two sequence modified variants (PepB, PepC). The C-terminus of PepA was amidated, whereas both PepB and PepC retained the carboxylate, and Ser16 in PepA and PepB was substituted with a helix-stabilizing residue, α-aminoisobutyric acid (Aib), in PepC. In thermal denaturation studies by far-UV CD spectroscopy and fibrillation kinetic studies by fluorescence and turbidity measurements, PepA and PepB showed heat-induced conformational changes and were found to form fibrils, whereas PepC did not fibrillate and showed only minor changes in the CD signal. Pulsed hydrogen-deuterium exchange mass spectrometry (HDX-MS) showed a high degree of protection from HD exchange in mature PepA fibrils and its proteolytic fragments, indicating that most of the sequence had been incorporated into the fibril structure and occurred nearly simultaneously throughout the sequence. The effects of the net peptide charge and formulation pH on fibrillation kinetics were investigated. In real-time stability studies of two formulations of PepA at pH's 7.4 and 8.0, analytical methods detected significant changes in the stability of the formulations at different time points during the study, which were not observed during accelerated studies. Additionally, PepA samples were withdrawn from real-time stability and subjected to additional stress (40 °C, continuous shaking) to induce fibrillation; an approach that successfully amplified oligomers or prefibrillar species previously undetected in a thioflavin T assay. Taken together, these studies present an approach to differentiate and characterize fibrillation risk in structurally related peptides under accelerated and real-time conditions, providing a model for rapid, iterative structural design to optimize the stability of therapeutic peptides.

Protein G-quadruplex interactions and their effects on phase transitions and protein aggregation

The SERF family of proteins were originally discovered for their ability to accelerate amyloid formation. Znf706 is an uncharacterized protein whose N-terminus is homologous to SERF proteins. We show here that human Znf706 can promote protein aggregation and amyloid formation. Unexpectedly, Znf706 specifically interacts with stable, non-canonical nucleic acid structures known as G-quadruplexes. G-quadruplexes can affect gene regulation and suppress protein aggregation; however, it is unknown if and how these two activities are linked. We find Znf706 binds preferentially to parallel G-quadruplexes with low micromolar affinity, primarily using its N-terminus, and upon interaction, its dynamics are constrained. G-quadruplex binding suppresses Znf706's ability to promote protein aggregation. Znf706 in conjunction with G-quadruplexes therefore may play a role in regulating protein folding. RNAseq analysis shows that Znf706 depletion specifically impacts the mRNA abundance of genes that are predicted to contain high G-quadruplex density. Our studies give insight into how proteins and G-quadruplexes interact, and how these interactions affect both partners and lead to the modulation of protein aggregation and cellular mRNA levels. These observations suggest that the SERF family of proteins, in conjunction with G-quadruplexes, may have a broader role in regulating protein folding and gene expression than previously appreciated.

Lithium Chloride Effects Field-Induced Protein Unfolding and the Transport Energetics Inside a Nanopipette

The tapered geometry of nanopipettes offers a unique perspective on protein transport through nanopores since both a gradual and fast confinement are possible depending on the translocation direction. The protein capture rate, unfolding, speed of translocation, and clogging probability are studied by toggling the LiCl concentration between 2 and 4 M. Interestingly, the proteins in this study could be transported with or against electrophoresis and offer vastly different attributes of sensing. Herein, a ruleset for studying proteins is developed that prevents irreversible pore clogging and yields upward of >100,000 events/nanopore. The extended duration of experiments further revealed that the capture rate takes ∼2 h to reach a steady state, emphasizing the importance of reaching equilibrated transport for studying the energetics and kinetics of protein transport (i.e., diffusion vs barrier-limited). Even in the equilibrated transport state, improper lowpass filtering was shown to distort the classification of diffusion-limited vs barrier-limited transport. Finally, electric-field-induced protein unfolding was found to be most prominent in electroosmotic-dominant transport, whereas electrophoretic-dominant events show no evidence of unfolding. Thus, our findings showcase the optimal conditions for protein translocations and the impact on studying protein unfolding, transporting energetics, and acquiring high bandwidth data.

Understanding the heterogeneity intrinsic to protein folding

Relating the native fold of a protein to its amino acid sequence remains a fundamental problem in biology. While computer algorithms have demonstrated recently their prowess in predicting what structure a particular amino acid sequence will fold to, an understanding of how and why a specific protein fold is achieved remains elusive. A major challenge is to define the role of conformational heterogeneity during protein folding. Recent experimental studies, utilizing time-resolved FRET, hydrogen-exchange coupled to mass spectrometry, and single-molecule force spectroscopy, often in conjunction with simulation, have begun to reveal how conformational heterogeneity evolves during folding, and whether an intermediate ensemble of defined free energy consists of different sub-populations of molecules that may differ significantly in conformation, energy and entropy.

Line Shape Analysis of 19F NMR-Monitored Chemical Denaturation of a Fold-Switching Protein RfaH Reveals Its Slow Folding Dynamics

The recent discovery of metamorphic proteins, which can switch between multiple conformations under native conditions, has challenged the well-established one sequence-one structure paradigm of protein folding. This is exemplified in the C-terminal domain of the multidomain transcription factor RfaH, which converts from an α-helical coiled-coil conformation in its autoinhibited state to a β-barrel conformation upon activation. Here, we use multisite line shape analysis of 19F NMR-monitored equilibrium chemical denaturation measurements of two 19F-labeled variants of full-length RfaH, to show that it folds/unfolds slowly on the NMR time scale, in an apparent all-or-none fashion at physiological pH and room temperature in the 3.3-4.8 M urea concentration range. The significant thermodynamic stability and slow unfolding rate (kinetic stability) are likely responsible for maintaining the closed autoinhibited state of RfaH, preventing spurious interactions with RNA polymerase (RNAP) when not functional. Our results provide a quantitative understanding of the folding-function relationship in the model fold-switching protein RfaH.

Variation of CI-2 Conformers upon Addition of Methanol to Water: An IMS-MS-Based Thermodynamic Analysis

Chymotrypsin inhibitor 2 (CI-2) is a well-studied, textbook example of a cooperative, two-state, native ↔ denatured folding transition. A recent hybrid ion mobility spectrometry (IMS)/mass spectrometry (MS) thermal denaturation study of CI-2 (the well-studied truncated 64-residue model) in water reported evidence that this two-state transition involves numerous (∼41) unique native and non-native (denatured) solution conformations. The characterization of so many, often low-abundance, states is possible because of the very high dynamic range of IMS-MS measurements of ionic species that are produced upon electrospraying CI-2 solutions from a variable temperature electrospray ionization source. A thermodynamic analysis of these states revealed large changes in enthalpy (ΔH) and entropy (ΔS) at different temperatures, and it was suggested that such variation might arise because of temperature-dependent conformational changes of the protein in response to changes in the conformational entropy and the dielectric permeability of water, which drops from a value of ε ∼ 79 at 24 °C to ∼ 60 at 82 °C. Herein, we examine how adding methanol to water influences the distributions of CI-2 conformers and their ensuing stabilities. The dielectric constant of a 60:40 water:methanol (MeOH) drops from ε ∼ 60 at 24 °C to ∼ 51 at 64 °C. Although the same set of conformers observed in water appears to be present in 60:40 water:MeOH, the abundance of each is substantially altered by the presence of methanol. Relative free energy values (ΔG) and thermodynamic values [ΔH and ΔS and heat capacities (ΔCp)] are derived from a Gibbs-Helmholtz analysis. A comparison of these data from water and water:MeOH systems allows rare insight into how variations in solvation and temperature affect many-state protein equilibria. While these studies confirm that variations in solvent dielectric constant with temperature affect the distributions of conformers that are observed, our findings suggest that other solvent differences may also affect abundances.

The temperature-dependence of host-guest binding thermodynamics: experimental and simulation studies

The thermodynamic parameters of host-guest binding can be used to describe, understand, and predict molecular recognition events in aqueous systems. However, interpreting binding thermodynamics remains challenging, even for these relatively simple molecules, as they are determined by both direct and solvent-mediated host-guest interactions. In this contribution, we focus on the contributions of water to binding by studying binding thermodynamics, both experimentally and computationally, for a series of nearly rigid, electrically neutral host-guest systems and report the temperature-dependent thermodynamic binding contributions ΔGb(T), ΔHb(T), ΔSb(T), and ΔCp,b. Combining isothermal titration calorimetry (ITC) measurements with molecular dynamics (MD) simulations, we provide insight into the binding forces at play for the macrocyclic hosts cucurbit[n]uril (CBn, n = 7-8) and β-cyclodextrin (β-CD) with a range of guest molecules. We find consistently negative changes in heat capacity on binding (ΔCp,b) for all systems studied herein - as well as for literature host-guest systems - indicating increased enthalpic driving forces for binding at higher temperatures. We ascribe these trends to solvation effects, as the solvent properties of water deteriorate as temperature rises. Unlike the entropic and enthalpic contributions to binding, with their differing signs and magnitudes for the classical and non-classical hydrophobic effect, heat capacity changes appear to be a unifying and more general feature of host-guest complex formation in water. This work has implications for understanding protein-ligand interactions and other complex systems in aqueous environments.

E. coli production of a multi-disulfide bonded SARS-CoV-2 Omicron BA.5 RBD exhibiting native-like biochemical and biophysical properties

Low-cost bacterial production of the receptor binding domain (RBD) of the SARS-CoV-2 Omicron spike protein holds significant potential in expediting the development of therapeutics against COVID-19. However, RBD contains eight cysteines forming four disulfide bonds, and expression in E. coli using standard protocols produces insoluble RBD forming non-native disulfide bonds. Here, we expressed RBD in E. coli T7 SHuffle with high aeration, which enhanced disulfide formation in the cytoplasm and reshuffling of non-native disulfide bonds, and at a low temperature of 16°C, which stabilized the native conformation and thus the formation of the native disulfide bonds. The yield of RBD was as high as 3 mg per 200 mL culture. We analyzed the conformational and biophysical properties of our E. coli-expressed RBD. First, the RP-HPLC elution profile indicated a single peak, suggesting that RBD was folded with a single disulfide bond pairing pattern. Next, circular dichroism analysis indicated a secondary structure content very close to that computed from the crystal structure. RBD's thermal denaturation monitored by CD was cooperative, strongly indicating a well-folded protein structure. Moreover, limited proteolysis showed that RBD was nearly as stable as RNase A, and the formation of native disulfide bonds was confirmed by LC-MS analysis. Furthermore, BLI analysis indicated a strong binding of RBD with the hACE2 with a dissociation constant of 0.83 nM, confirming the folded nature of RBD. Altogether, these results demonstrate that our E. coli-expression system can provide a large amount of highly purified RBD with correct disulfide bonds and native-like biochemical and biophysical properties.

Dominance vs epistasis: the biophysical origins and plasticity of genetic interactions within and between alleles

An important challenge in genetics, evolution and biotechnology is to understand and predict how mutations combine to alter phenotypes, including molecular activities, fitness and disease. In diploids, mutations in a gene can combine on the same chromosome or on different chromosomes as a "heteroallelic combination". However, a direct comparison of the extent, sign, and stability of the genetic interactions between variants within and between alleles is lacking. Here we use thermodynamic models of protein folding and ligand-binding to show that interactions between mutations within and between alleles are expected in even very simple biophysical systems. Protein folding alone generates within-allele interactions and a single molecular interaction is sufficient to cause between-allele interactions and dominance. These interactions change differently, quantitatively and qualitatively as a system becomes more complex. Altering the concentration of a ligand can, for example, switch alleles from dominant to recessive. Our results show that intra-molecular epistasis and dominance should be widely expected in even the simplest biological systems but also reinforce the view that they are plastic system properties and so a formidable challenge to predict. Accurate prediction of both intra-molecular epistasis and dominance will require either detailed mechanistic understanding and experimental parameterization or brute-force measurement and learning.

Effect of temperature and nucleotide on the binding of BiP chaperone to a protein substrate

BiP (immunoglobulin heavy-chain binding protein) is a Hsp70 monomeric ATPase motor that plays broad and crucial roles in maintaining proteostasis inside the cell. Structurally, BiP is formed by two domains, a nucleotide-binding domain (NBD) with ATPase activity connected by a flexible hydrophobic linker to the substrate-binding domain. While the ATPase and substrate binding activities of BiP are allosterically coupled, the latter is also dependent on nucleotide binding. Recent structural studies have provided new insights into BiP's allostery; however, the influence of temperature on the coupling between substrate and nucleotide binding to BiP remains unexplored. Here, we study BiP's binding to its substrate at the single molecule level using thermo-regulated optical tweezers which allows us to mechanically unfold the client protein and explore the effect of temperature and different nucleotides on BiP binding. Our results confirm that the affinity of BiP for its protein substrate relies on nucleotide binding, by mainly regulating the binding kinetics between BiP and its substrate. Interestingly, our findings also showed that the apparent affinity of BiP for its protein substrate in the presence of nucleotides remains invariable over a wide range of temperatures, suggesting that BiP may interact with its client proteins with similar affinities even when the temperature is not optimal. Thus, BiP could play a role as a "thermal buffer" in proteostasis.

Unfolding and Aggregation Pathways of Variable Domains from Immunoglobulin Light Chains

Light chain amyloidosis is the most common form of systemic amyloidosis. This disease is caused by the formation and deposition of amyloid fibers made from immunoglobulin light chains. Environmental conditions such as pH and temperature can affect protein structure and induce the development of these fibers. Several studies have shed light on the native state, stability, dynamics, and final amyloid state of these proteins; however, the initiation process and the fibril formation pathway remain poorly understood structurally and kinetically. To study this, we analyzed the unfolding and aggregation process of the 6aJL2 protein under acidic conditions, with temperature changes, and upon mutation, using biophysical and computational techniques. Our results suggest that the differences in amyloidogenicity displayed by 6aJL2 under these conditions are caused by traversing different aggregation pathways, including unfolded intermediates and the formation of oligomers.

Urea induced unfolding of rai seed cystatin: Influence of glycerol as a chemical chaperone

Cystatin superfamily members, by virtue of their thiol protease regulatory properties, show involvement in myriad physiological processes important for survival and well-being. The current study involves urea-induced denaturation of a novel variant of the cystatin superfamily, rai seed cystatin (RSC), employing a variety of biophysical assays in order to characterize different folding intermediates generated on unfolding. Urea as a denaturant presented the passage of RSC through a series of events resulting in the loss of RSC functional capability, accompanied by changes in the archetype at secondary and tertiary structural levels, as evident from protease inhibitory, UV absorption, and intrinsic fluorescence assays, respectively. ANS fluorescence also revealed routing of RSC through discrete multiple sub-states thus presenting the generation of intermediate states somewhat close to the pre-molten globule and/or molten globule forms of RSC. Furthermore, far-UV circular dichroism analysis revealed a concentration-dependent gradual loss in typical -helical RSC peaks, indicating a nearly 50 % loss in secondary structural elements around 5 M urea treatment. The study also reports the possible role of glycerol in the refolding and/or reactivation of the urea unfolded RSC form. Glycerol presented itself as a potent structural stabilizer as it assisted in the refolding and reactivation of the unfolded RSC in a dosage-dependent manner, concomitantly paving the way for unravelling the mechanistic approach involved in the phenomenon, which can facilitate future studies.

1700029I15Rik orchestrates the biosynthesis of acrosomal membrane proteins required for sperm-egg interaction

Sperm acrosomal membrane proteins, such as Izumo sperm-egg fusion 1 (IZUMO1) and sperm acrosome-associated 6 (SPACA6), play essential roles in mammalian gamete binding or fusion. How their biosynthesis is regulated during spermiogenesis has largely remained elusive. Here, we show that 1700029I15Rik knockout male mice are severely subfertile and their spermatozoa do not fuse with eggs. 1700029I15Rik is a type-II transmembrane protein expressed in early round spermatids but not in mature spermatozoa. It interacts with proteins involved in N-linked glycosylation, disulfide isomerization, and endoplasmic reticulum (ER)-Golgi trafficking, suggesting a potential role in nascent protein processing. The ablation of 1700029I15Rik destabilizes non-catalytic subunits of the oligosaccharyltransferase (OST) complex that are pivotal for N-glycosylation. The knockout testes exhibit normal expression of sperm plasma membrane proteins, but decreased abundance of multiple acrosomal membrane proteins involved in fertilization. The knockout sperm show upregulated chaperones related to ER-associated degradation (ERAD) and elevated protein ubiquitination; strikingly, SPACA6 becomes undetectable. Our results support for a specific, 1700029I15Rik-mediated pathway underpinning the biosynthesis of acrosomal membrane proteins during spermiogenesis.

Cooperative Protein Allosteric Transition Mediated by a Fluctuating Transmission Network

Allosteric communication between distant protein sites represents a key mechanism of biomolecular regulation and signal transduction. Compared to other processes such as protein folding, however, the dynamical evolution of allosteric transitions is still not well understood. As an example of allosteric coupling between distant protein regions, we consider the global open-closed motion of the two domains of T4 lysozyme, which is triggered by local motion in the hinge region. Combining extensive molecular dynamics simulations with a correlation analysis of interresidue contacts, we identify a network of interresidue distances that move in a concerted manner. The cooperative process originates from a cogwheel-like motion of the hydrophobic core in the hinge region, which constitutes an evolutionary conserved and flexible transmission network. Through rigid contacts and the protein backbone, the small local changes of the hydrophobic core are passed on to the distant terminal domains and lead to the emergence of a rare global conformational transition. As in an Ising-type model, the cooperativity of the allosteric transition can be explained via the interaction of local fluctuations.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

From Protein Design to the Energy Landscape of a Cold Unfolding Protein

Understanding protein folding is crucial for protein sciences. The conformational spaces and energy landscapes of cold (unfolded) protein states, as well as the associated transitions, are hardly explored. Furthermore, it is not known how structure relates to the cooperativity of cold transitions, if cold and heat unfolded states are thermodynamically similar, and if cold states play important roles for protein function. We created the cold unfolding 4-helix bundle DCUB1 with a de novo designed bipartite hydrophilic/hydrophobic core featuring a hydrogen bond network which extends across the bundle in order to study the relative importance of hydrophobic versus hydrophilic protein-water interactions for cold unfolding. Structural and thermodynamic characterization resulted in the discovery of a complex energy landscape for cold transitions, while the heat unfolded state is a random coil. Below ∼0 °C, the core of DCUB1 disintegrates in a largely cooperative manner, while a near-native helical content is retained. The resulting cold core-unfolded state is compact and features extensive internal dynamics. Below -5 °C, two additional cold transitions are seen, that is, (i) the formation of a water-mediated, compact, and highly dynamic dimer, and (ii) the onset of cold helix unfolding decoupled from cold core unfolding. Our results suggest that cold unfolding is initiated by the intrusion of water into the hydrophilic core network and that cooperativity can be tuned by varying the number of core hydrogen bond networks. Protein design has proven to be invaluable to explore the energy landscapes of cold states and to robustly test related theories.

The Use of the Statistical Entropy in Some New Approaches for the Description of Biosystems

Some problems of describing biological systems with the use of entropy as a measure of the complexity of these systems are considered. Entropy is studied both for the organism as a whole and for its parts down to the molecular level. Correlation of actions of various parts of the whole organism, intercellular interactions and control, as well as cooperativity on the microlevel lead to a more complex structure and lower statistical entropy. For a multicellular organism, entropy is much lower than entropy for the same mass of a colony of unicellular organisms. Cooperativity always reduces the entropy of the system; a simple example of ligand binding to a macromolecule carrying two reaction centers shows how entropy is consistent with the ambiguity of the result in the Bernoulli test scheme. Particular attention is paid to the qualitative and quantitative relationship between the entropy of the system and the cooperativity of ligand binding to macromolecules. A kinetic model of metabolism. corresponding to Schrödinger’s concept of the maintenance biosystems by “negentropy feeding”, is proposed. This model allows calculating the nonequilibrium local entropy and comparing it with the local equilibrium entropy inherent in non-living matter.

Native State Hydrogen Exchange-Mass Spectrometry Methods to Probe Protein Folding and Unfolding

Native state hydrogen exchange (HX) methods provide high-resolution structural data on the rare and transient opening motions in proteins under native conditions. Mass spectrometry-based HX methods (HX-MS) have gained popularity because of their ability to delineate population distributions, which allow a direct determination of the mechanism of inter conversion of the partially folded states under native conditions. Various technological advancements have provided further impetus to the development of HX-MS-based experiments to study protein folding. Classical HX-MS studies use proteolytic digestion to produce fragments of the protein subsequent to HX in solution, in order to obtain structural data. New chemical fragmentation methods, which achieve the same result as proteolysis and cause minimal change to the HX pattern in the protein, provide an attractive alternative to proteolysis. Moreover, when used in conjunction with proteolysis, chemical fragmentation methods have significantly increased the structural resolution afforded by HX-MS studies, even bringing them at par with the single amino acid resolution observed in NMR-based measurements. Experiments based on one such chemical fragmentation method, electron transfer dissociation (ETD), are described in this chapter. The ETD HX-MS method is introduced using data from a protein which is inherently resistant to proteolytic digestion as example of how such an experiment can provide high-resolution structural data on the folding-unfolding transitions of the protein under native conditions.

Slowest-first protein translation scheme: Structural asymmetry and co-translational folding

Proteins are translated from the N to the C terminus, raising the basic question of how this innate directionality affects their evolution. To explore this question, we analyze 16,200 structures from the Protein Data Bank (PDB). We find remarkable enrichment of α helices at the C terminus and β strands at the N terminus. Furthermore, this α-β asymmetry correlates with sequence length and contact order, both determinants of folding rate, hinting at possible links to co-translational folding (CTF). Hence, we propose the "slowest-first" scheme, whereby protein sequences evolved structural asymmetry to accelerate CTF: the slowest of the cooperatively folding segments are positioned near the N terminus so they have more time to fold during translation. A phenomenological model predicts that CTF can be accelerated by asymmetry in folding rate, up to double the rate, when folding time is commensurate with translation time; analysis of the PDB predicts that structural asymmetry is indeed maximal in this regime. This correspondence is greater in prokaryotes, which generally require faster protein production. Altogether, this indicates that accelerating CTF is a substantial evolutionary force whose interplay with stability and functionality is encoded in secondary structure asymmetry.

Folding Intermediates, Heterogeneous Native Ensembles and Protein Function

Single domain proteins fold via diverse mechanisms emphasizing the intricate relationship between energetics and structure, which is a direct consequence of functional constraints and demands imposed at the level of sequence. On the other hand, elucidating the interplay between folding mechanisms and function is challenging in large proteins, given the inherent shortcomings in identifying metastable states experimentally and the sampling limitations associated with computational methods. Here, we show that free energy profiles and surfaces of large systems (>150 residues), as predicted by a statistical mechanical model, display a wide array of folding mechanisms with ubiquitous folding intermediates and heterogeneous native ensembles. Importantly, residues around the ligand binding or enzyme active site display a larger tendency to partially unfold and this manifests as intermediates or excited states along the folding coordinate in ligand binding domains, transcription repressors, and representative enzymes from all the six classes, including the SARS-CoV-2 receptor binding domain (RBD) of the spike protein and the protease M^pro. It thus appears that it is relatively easier to distill the imprints of function on the folding landscape of larger proteins as opposed to smaller systems. We discuss how an understanding of energetic-entropic features in ordered proteins can pinpoint specific avenues through which folding mechanisms, populations of partially structured states and function can be engineered.

Correlation in Domain Fluctuations Navigates Target Search of a Viral Peptide along RNA

Biological macromolecules often exhibit correlations in fluctuations involving distinct domains. This study decodes their functional implications in RNA-protein recognition and target-specific binding. The target search of a peptide along RNA in a viral TAR-Tat complex is closely monitored using atomistic simulations, steered molecular dynamics simulations, free energy calculations, and a machine-learning-based clustering technique. An anticorrelated domain fluctuation is identified between the tetraloop and the bulge region in the apo form of TAR RNA that sets a hierarchy in the domain-specific fluctuations at each binding event and that directs the succeeding binding footsteps. Thus, at each binding footstep, the dynamic partner selects an RNA location for binding where it senses a higher fluctuation, which is conventionally reduced upon binding. This event stimulates an alternate domain fluctuation, which then dictates sequential binding footstep/s and thus the search progresses. Our cross-correlation maps show that the fluctuations relay from one domain to another specific domain until the anticorrelation between those interdomain fluctuations sustains. Artificial attenuation of that hierarchical domain fluctuation inhibits specific RNA binding. The binding is completed with the arrival of a few long-lived water molecules that mediate slightly distant RNA-protein sites and finally stabilize the overall complex. The study underscores the functional importance of naturally designed fluctuating RNA motifs (bulge, tetraloop) and their interplay in dictating the directionality of the search in a highly dynamic environment.

A hierarchy of coupling free energies underlie the thermodynamic and functional architecture of protein structures

Protein sequences and structures evolve by satisfying varied physical and biochemical constraints. This multi-level selection is enabled not just by the patterning of amino acids on the sequence, but also via coupling between residues in the native structure. Here, we employ an energetically detailed statistical mechanical model with millions of microstates to extract such long-range structural correlations, i.e. thermodynamic coupling free energies, from a diverse family of protein structures. We find that despite the intricate and anisotropic distribution of coupling patterns, the majority of residues (>70%) are only marginally coupled contributing to functional motions and catalysis. Physical origins of ‘sectors’, determinants of native ensemble heterogeneity in extant, ancient and designed proteins, and the basis for allostery emerge naturally from coupling free energies. The statistical framework highlights how evolutionary selection and optimization occur at the level of global interaction network for a given protein fold impacting folding, function, and allosteric outputs.

•

Evolution of protein structures occurs at the level of global interaction network.

•

More than 70% of the protein residues are weakly or marginally coupled.

•

Functional ‘sector’ regions are a manifestation of marginal coupling.

•

Coupling indices vary across the entire proteins in extant-ancient and natural-designed pairs.

•

The proposed methodology can be used to understand allostery and epistasis.

Computational and Experimental Assessment of Backbone Templates for Computational Redesign of the Thioredoxin Fold

Computational protein design has taken big strides in recent years; however, the tools available are still not at a state where a sequence can be designed to fold into a given protein structure at will and with high probability. We have applied here a recent release of Rosetta Design to redesign a set of structurally very similar proteins belonging to the thioredoxin fold. We used a genetic screening tool to estimate solubility/folding of the designed proteins in E. coli and to select the best hits from this for further biochemical characterization. We have previously used this set of template proteins for redesign and found that success was highly dependent on template structure, a trait which was also found in this study. Nevertheless, state-of-the-art design software is now able to predict the best template, most likely due to the introduction of an energy term that reports on stress in covalent bond lengths and angles. The template that led to the greatest fraction of successful designs was the same (a thioredoxin from spinach) as that identified in our previous study. Our previously described redesign of thioredoxin, which also used the spinach protein as a template, however also performed well as a template. In the present study, both of these templates yielded proteins with compact folded structures and enforced the conclusion that any design project must carefully consider different design templates. Fortunately, selecting designs based on energies appears to correctly identify such templates.

Biophysical and structural studies reveal marginal stability of a crucial hydrocarbon biosynthetic enzyme acyl ACP reductase

Acyl-ACP reductase (AAR) is one of the two key cyanobacterial enzymes along with aldehyde deformylating oxygenase (ADO) involved in the synthesis of long-chain alkanes, a drop-in biofuel. The enzyme is prone to aggregation when expressed in Escherichia coli, leading to varying alkane levels. The present work attempts to investigate the crucial structural aspects of AAR protein associated with its stability and folding. Characterization by dynamic light scattering experiment and intact mass spectrometry revealed that recombinantly expressed AAR in E. coli existed in multiple-sized protein particles due to diverse lipidation. Interestingly, while thermal- and urea-based denaturation of AAR showed 2-state unfolding transition in circular dichroism and intrinsic fluorescent spectroscopy, the unfolding process of AAR was a 3-state pathway in GdnHCl solution suggesting that the protein milieu plays a significant role in dictating its folding. Apparent standard free energy \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( {\Delta {\text{G}}_{{{\text{NU}}}}^{{{\text{H}}_{2} {\text{O}}}} } \right)$$\end{document}ΔGNUH2O of ~ 4.5 kcal/mol for the steady-state unfolding of AAR indicated borderline stability of the protein. Based on these evidences, we propose that the marginal stability of AAR are plausible contributing reasons for aggregation propensity and hence the low catalytic activity of the enzyme when expressed in E. coli for biofuel production. Our results show a path for building superior biocatalyst for higher biofuel production.

Pressure and Chemical Unfolding of an α-Helical Bundle Protein: The GH2 Domain of the Protein Adaptor GIPC1

When combined with NMR spectroscopy, high hydrostatic pressure is an alternative perturbation method used to destabilize globular proteins that has proven to be particularly well suited for exploring the unfolding energy landscape of small single-domain proteins. To date, investigations of the unfolding landscape of all-β or mixed-α/β protein scaffolds are well documented, whereas such data are lacking for all-α protein domains. Here we report the NMR study of the unfolding pathways of GIPC1-GH2, a small α-helical bundle domain made of four antiparallel α-helices. High-pressure perturbation was combined with NMR spectroscopy to unravel the unfolding landscape at three different temperatures. The results were compared to those obtained from classical chemical denaturation. Whatever the perturbation used, the loss of secondary and tertiary contacts within the protein scaffold is almost simultaneous. The unfolding transition appeared very cooperative when using high pressure at high temperature, as was the case for chemical denaturation, whereas it was found more progressive at low temperature, suggesting the existence of a complex folding pathway.

Selecting Subpopulations of High-Quality Protein Conformers among Conformational Mixtures of Recombinant Bovine MMP-9 Solubilized from Inclusion Bodies

A detailed workflow to analyze the physicochemical characteristics of mammalian matrix metalloproteinase (MMP-9) protein species obtained from protein aggregates (inclusion bodies-IBs) was followed. MMP-9 was recombinantly produced in the prokaryotic microbial cell factories Clearcoli (an engineered form of Escherichia coli) and Lactococcus lactis, mainly forming part of IBs and partially recovered under non-denaturing conditions. After the purification by affinity chromatography of solubilized MMP-9, four protein peaks were obtained. However, so far, the different conformational protein species forming part of IBs have not been isolated and characterized. Therefore, with the aim to link the physicochemical characteristics of the isolated peaks with their biological activity, we set up a methodological approach that included dynamic light scattering (DLS), circular dichroism (CD), and spectrofluorometric analysis confirming the separation of subpopulations of conformers with specific characteristics. In protein purification procedures, the detailed analysis of the individual physicochemical properties and the biological activity of protein peaks separated by chromatographic techniques is a reliable source of information to select the best-fitted protein populations.

Designed folding pathway of modular coiled-coil-based proteins

Natural proteins are characterised by a complex folding pathway defined uniquely for each fold. Designed coiled-coil protein origami (CCPO) cages are distinct from natural compact proteins, since their fold is prescribed by discrete long-range interactions between orthogonal pairwise-interacting coiled-coil (CC) modules within a single polypeptide chain. Here, we demonstrate that CCPO proteins fold in a stepwise sequential pathway. Molecular dynamics simulations and stopped-flow Förster resonance energy transfer (FRET) measurements reveal that CCPO folding is dominated by the effective intra-chain distance between CC modules in the primary sequence and subsequent folding intermediates, allowing identical CC modules to be employed for multiple cage edges and thus relaxing CCPO cage design requirements. The number of orthogonal modules required for constructing a CCPO tetrahedron can be reduced from six to as little as three different CC modules. The stepwise modular nature of the folding pathway offers insights into the folding of tandem repeat proteins and can be exploited for the design of modular protein structures based on a given set of orthogonal modules.

Mapping Distinct Sequences of Structure Formation Differentiating Multiple Folding Pathways of a Small Protein

To determine experimentally how the multiple folding pathways of a protein differ, in the order in which the structural parts are assembled, has been a long-standing challenge. To resolve whether structure formation during folding can progress in multiple ways, the complex folding landscape of monellin has been characterized, structurally and temporally, using the multisite time-resolved FRET methodology. After an initial heterogeneous polypeptide chain collapse, structure formation proceeds on parallel pathways. Kinetic analysis of the population evolution data across various protein segments provides a clear structural distinction between the parallel pathways. The analysis leads to a phenomenological model that describes how and when discrete segments acquire structure independently of each other in different subensembles of protein molecules. When averaged over all molecules, structure formation is seen to progress as α-helix formation, followed by core consolidation, then β-sheet formation, and last end-to-end distance compaction. Parts of the protein that are closer in the primary sequence acquire structure before parts separated by longer sequence.

Insight to Functional Conformation and Noncovalent Interactions of Protein-Protein Assembly Using MALDI Mass Spectrometry

Noncovalent interactions are the keys to the structural organization of biomolecule e.g., proteins, glycans, lipids in the process of molecular recognition processes e.g., enzyme-substrate, antigen-antibody. Protein interactions lead to conformational changes, which dictate the functionality of that protein-protein complex. Besides biophysics techniques, noncovalent interaction and conformational dynamics, can be studied via mass spectrometry (MS), which represents a powerful tool, due to its low sample consumption, high sensitivity, and label-free sample. In this review, the focus will be placed on Matrix-Assisted Laser Desorption Ionization Mass Spectrometry (MALDI-MS) and its role in the analysis of protein-protein noncovalent assemblies exploring the relationship within noncovalent interaction, conformation, and biological function.

Synthesis runs counter to directional folding of a nascent protein domain

Folding of individual domains in large proteins during translation helps to avoid otherwise prevalent inter-domain misfolding. How folding intermediates observed in vitro for the majority of proteins relate to co-translational folding remains unclear. Combining in vivo and single-molecule experiments, we followed the co-translational folding of the G-domain, encompassing the first 293 amino acids of elongation factor G. Surprisingly, the domain remains unfolded until it is fully synthesized, without collapsing into molten globule-like states or forming stable intermediates. Upon fully emerging from the ribosome, the G-domain transitions to its stable native structure via folding intermediates. Our results suggest a strictly sequential folding pathway initiating from the C-terminus. Folding and synthesis thus proceed in opposite directions. The folding mechanism is likely imposed by the final structure and might have evolved to ensure efficient, timely folding of a highly abundant and essential protein.

In vivo experiments and optical tweezers force-spectroscopy measurements assessing the co-translational folding of the G-domain from bacterial elongation factor G reveal a sequential folding pathway initiating from the C-terminus. These results suggest that protein folding and synthesis proceed in opposite directions.

A Change in Domain Cooperativity Drives the Function of Calnuc

With the increasing incidence of neurodegenerative disorders, there is an urgent need to understand the protein folding process. Examining the folding process of multidomain proteins remains a prime challenge, as their complex conformational dynamics make them highly susceptible to misfolding and/or aggregation. The presence of multiple domains in a protein can lead to interaction between the partially folded domains, thereby driving misfolding and/or aggregation. Calnuc is one such multidomain protein for which Ca²⁺ binding plays a pivotal role in governing its structural dynamics and stability and, presumably, in directing its interactions with other proteins. We demonstrate differential structural dynamics between the Ca²⁺-free and Ca²⁺-bound forms of calnuc. In the absence of Ca²⁺, full-length calnuc displays equilibrium structural transitions with four intermediate states, reporting a sum of the behavioral properties of its individual domains. Fragment-based studies illustrate the sequential events of structure adoption proceeding in the following order: EF domain followed by the NT and LZ domains in the apo state. On the other hand, Ca²⁺ binding increases domain cooperativity and enables the protein to fold as a single unit. Single-tryptophan mutant proteins, designed in a domain-dependent manner, confirm an increase in the number of interdomain interactions in the Ca²⁺-bound form as compared to the Ca²⁺-free state of the protein, thereby providing insight into its folding process. The attenuated domain crosstalk in apo-calnuc is likely to influence and regulate its physiologically important intermolecular interactions.

Mosaic Cooperativity in Slow Polypeptide Topological Isomerization Revealed by Residue-Specific NMR Thermodynamic Analysis

Slow polypeptide conformational changes on time scales of >1 s are generally assumed to be highly cooperative two-state transitions, reflecting the high energy barrier. However, few experimental characterizations have tested the validity of this assumption. We performed residue-specific NMR thermodynamic analysis of the 27-residue lantibiotic peptide, nukacin ISK-1, to characterize the isomerization between two topological states on the second time scale. Unexpectedly, the thermal transition behaviors were distinct among peptide regions, indicating that the topological isomerization process is a mosaic of different degrees of cooperativity. The conformational change path between the two NMR structures was deduced by a targeted molecular dynamics simulation. The unique side-chain threading motions through the monosulfide rings are the structural basis of the high energy barrier, and the nonlocal interactions in the hydrophobic core are the structural basis of the cooperativity. Taken together, we provide an energetic description of the topological isomerization of nukacin ISK-1.

Water as a Good Solvent for Unfolded Proteins: Folding and Collapse are Fundamentally Different

The argument that the hydrophobic effect is the primary effect driving the folding of globular proteins is nearly universally accepted (including by the authors). But does this view also imply that water is a "poor" solvent for the unfolded states of these same proteins? Here we argue that the answer is "no," that is, folding to a well-packed, extensively hydrogen-bonded native structure differs fundamentally from the nonspecific chain collapse that defines a poor solvent. Thus, the observation that a protein folds in water does not necessitate that water is a poor solvent for its unfolded state. Indeed, chain-solvent interactions that are marginally more favorable than nonspecific intrachain interactions are beneficial to protein function because they destabilize deleterious misfolded conformations and inter-chain interactions.

Stability and Conformational Resilience of Protein Disulfide Isomerase

Protein disulfide isomerase (PDI) is a redox-dependent protein with oxidoreductase and chaperone activities. It is a U-shaped protein with an abb'xa' structural organization in which the a and a' domains have CGHC active sites, the b and b' domains are involved with substrate binding, and x is a flexible linker. PDI exhibits substantial flexibility and undergoes cycles of unfolding and refolding in its interaction with cholera toxin, suggesting PDI can regain a folded, functional conformation after exposure to stress conditions. To determine whether this unfolding-refolding cycle is a substrate-induced process or an intrinsic physical property of PDI, we used circular dichroism to examine the structural properties of PDI subjected to thermal denaturation. PDI exhibited remarkable conformational resilience that is linked to its redox status. In the reduced state, PDI exhibited a 54 °C unfolding transition temperature (Tm) and regained 85% of its native structure after nearly complete thermal denaturation. Oxidized PDI had a lower Tm of 48-50 °C and regained 70% of its native conformation after 75% denaturation. Both reduced PDI and oxidized PDI were functional after refolding from these denatured states. Additional studies documented increased stability of a PDI construct lacking the a' domain and decreased thermal stability of a construct lacking the a domain. Furthermore, oxidation of the a domain limited the ability of PDI to refold. The stability and conformational resilience of PDI are thus linked to both redox-dependent and domain-specific effects. These findings document previously unrecognized properties of PDI and provide insight into the physical foundation of its biological function.

Thiamin and protein folding

A huge number of proteins that occur in the body have to be folded into a specific shape in order to become functional. Proteins are made up of chains of amino acids and the folding process is exquisitely complex. When this folding process is inhibited, the respective protein is referred to as being misfolded and nonfunctional. So the hypothesis that follows is in regard to the diseases that are caused by the misfolding of vital proteins and their reported relationship with thiamin metabolism. These diseases are termed proteopathies and there are at least 50 different conditions in which the mechanism is importantly related to a misfolded protein. In the brain, thiamin deficiency causes a cascade of events involving mild impairment of oxidative metabolism, neuroinflammation and neurodegeneration, including the pathology of Alzheimer's disease, Parkinson's and Huntington's diseases, all of which are examples of proteopathies. Prion diseases are fatal neurodegenerative disorders related to the conformational alteration of the prion protein (PrP C) into a pathogenic and protease-resistant isoform (PrPSc). The physiological form (PrP C) is a cell surface glycoprotein expressed mainly in the central nervous system. Despite numerous efforts to elucidate its role, the exact biological function remains unknown. Prion-induced diseases, due to the conformational change in the protein, are a global health problem, with lack of effective therapy and 100% mortality. Thiamin and its derivatives bind the prion protein and intermolecular actions have been noted between thiamin and other thiamin-binding proteins, although the exact importance of this is conjectural.

An arsenal of methods for the experimental characterization of intrinsically disordered proteins - How to choose and combine them?

In this review, we detail the most common experimental approaches to assess and characterize protein intrinsic structural disorder, with the notable exception of NMR and EPR spectroscopy, two ideally suited approaches that will be described in depth in two other reviews within this special issue. We discuss the advantages, the limitations, as well as the caveats of the various methods. We also describe less common and more demanding approaches that enable achieving further insights into the conformational properties of IDPs. Finally, we present recent developments that have enabled assessment of structural disorder in living cells, and discuss the currently available methods to model IDPs as conformational ensembles.

Engineering Order and Cooperativity in a Disordered Protein

Structural disorder in proteins arises from a complex interplay between weak hydrophobicity and unfavorable electrostatic interactions. The extent to which the hydrophobic effect contributes to the unique and compact native state of proteins is, however, confounded by large compensation between multiple entropic and energetic terms. Here we show that protein structural order and cooperativity arise as emergent properties upon hydrophobic substitutions in a disordered system with non-intuitive effects on folding and function. Aided by sequence-structure analysis, equilibrium, and kinetic spectroscopic studies, we engineer two hydrophobic mutations in the disordered DNA-binding domain of CytR that act synergistically, but not in isolation, to promote structure, compactness, and stability. The double mutant, with properties of a fully ordered domain, exhibits weak cooperativity with a complex and rugged conformational landscape. The mutant, however, binds cognate DNA with an affinity only marginally higher than that of the wild type, though nontrivial differences are observed in the binding to noncognate DNA. Our work provides direct experimental evidence of the dominant role of non-additive hydrophobic effects in shaping the molecular evolution of order in disordered proteins and vice versa, which could be generalized to even folded proteins with implications for protein design and functional manipulation.

Pressure-Temperature Analysis of the Stability of the CTL9 Domain Reveals Hidden Intermediates

The observation of two-state unfolding for many small single-domain proteins by denaturants has led to speculation that protein sequences may have evolved to limit the population of partially folded states that could be detrimental to fitness. How such strong cooperativity arises from a multitude of individual interactions is not well understood. Here, we investigate the stability and folding cooperativity of the C-terminal domain of the ribosomal protein L9 in the pressure-temperature plane using site-specific NMR. In contrast to apparent cooperative unfolding detected with denaturant-induced and thermal-induced unfolding experiments and stopped-flow refolding studies at ambient pressure, NMR-detected pressure unfolding revealed significant deviation from two-state behavior, with a core region that was selectively destabilized by increasing temperature. Comparison of pressure-dependent NMR signals from both the folded and unfolded states revealed the population of at least one invisible excited state at atmospheric pressure. The core destabilizing cavity-creating I98A mutation apparently increased the cooperativity of the loss of folded-state peak intensity while also increasing the population of this invisible excited state present at atmospheric pressure. These observations highlight how local stability is subtly modulated by sequence to tune protein conformational landscapes and illustrate the ability of pressure- and temperature-dependent studies to reveal otherwise hidden states.