Are Protein Folding Intermediates the Evolutionary Consequence of Functional Constraints?

Conformational Tuning Shapes the Balance between Functional Promiscuity and Specialization in Paralogous Plasmodium Acyl-CoA Binding Proteins

Paralogous proteins confer enhanced fitness to organisms via complex sequence-conformation codes that shape functional divergence, specialization, or promiscuity. Here, we dissect the underlying mechanism of promiscuous binding versus partial subfunctionalization in paralogues by studying structurally identical acyl-CoA binding proteins (ACBPs) from Plasmodium falciparum that serve as promising drug targets due to their high expression during the protozoan proliferative phase. Combining spectroscopic measurements, solution NMR, SPR, and simulations on two of the paralogues, A16 and A749, we show that minor sequence differences shape nearly every local and global conformational feature. A749 displays a broader and heterogeneous native ensemble, weaker thermodynamic coupling and cooperativity, enhanced fluctuations, and a larger binding pocket volume compared to A16. Site-specific tryptophan probes signal a graded reduction in the sampling of substates in the holo form, which is particularly apparent in A749. The paralogues exhibit a spectrum of binding affinities to different acyl-CoAs with A749, the more promiscuous and hence the likely ancestor, binding 1000-fold stronger to lauroyl-CoA under physiological conditions. We thus demonstrate how minor sequence changes modulate the extent of long-range interactions and dynamics, effectively contributing to the molecular evolution of contrasting functional repertoires in paralogues.

Phosphorylation of Thr9 Affects the Folding Landscape of the N-Terminal Segment of Human AGT Enhancing Protein Aggregation of Disease-Causing Mutants

The mutations G170R and I244T are the most common disease cause in primary hyperoxaluria type I (PH1). These mutations cause the misfolding of the AGT protein in the minor allele AGT-LM that contains the P11L polymorphism, which may affect the folding of the N-terminal segment (NTT-AGT). The NTT-AGT is phosphorylated at T9, although the role of this event in PH1 is unknown. In this work, phosphorylation of T9 was mimicked by introducing the T9E mutation in the NTT-AGT peptide and the full-length protein. The NTT-AGT conformational landscape was studied by circular dichroism, NMR, and statistical mechanical methods. Functional and stability effects on the full-length AGT protein were characterized by spectroscopic methods. The T9E and P11L mutations together reshaped the conformational landscape of the isolated NTT-AGT peptide by stabilizing ordered conformations. In the context of the full-length AGT protein, the T9E mutation had no effect on the overall AGT function or conformation, but enhanced aggregation of the minor allele (LM) protein and synergized with the mutations G170R and I244T. Our findings indicate that phosphorylation of T9 may affect the conformation of the NTT-AGT and synergize with PH1-causing mutations to promote aggregation in a genotype-specific manner. Phosphorylation should be considered a novel regulatory mechanism in PH1 pathogenesis.

The Wako-Saitô-Muñoz-Eaton Model for Predicting Protein Folding and Dynamics

Despite the recent advances in the prediction of protein structures by deep neutral networks, the elucidation of protein-folding mechanisms remains challenging. A promising theory for describing protein folding is a coarse-grained statistical mechanical model called the Wako-Saitô-Muñoz-Eaton (WSME) model. The model can calculate the free-energy landscapes of proteins based on a three-dimensional structure with low computational complexity, thereby providing a comprehensive understanding of the folding pathways and the structure and stability of the intermediates and transition states involved in the folding reaction. In this review, we summarize previous and recent studies on protein folding and dynamics performed using the WSME model and discuss future challenges and prospects. The WSME model successfully predicted the folding mechanisms of small single-domain proteins and the effects of amino-acid substitutions on protein stability and folding in a manner that was consistent with experimental results. Furthermore, extended versions of the WSME model were applied to predict the folding mechanisms of multi-domain proteins and the conformational changes associated with protein function. Thus, the WSME model may contribute significantly to solving the protein-folding problem and is expected to be useful for predicting protein folding, stability, and dynamics in basic research and in industrial and medical applications.

Predicting and Simulating Mutational Effects on Protein Folding Kinetics

Mutational perturbations of protein structures, i.e., phi-value analysis, are commonly employed to probe the extent of involvement of a particular residue in the rate-determining step(s) of folding. This generally involves the measurement of folding thermodynamic parameters and kinetic rate constants for the wild-type and mutant proteins. While computational approaches have been reasonably successful in understanding and predicting the effect of mutations on folding thermodynamics, it has been challenging to explore the same on kinetics due to confounding structural, energetic, and dynamic factors. Accordingly, the frequent observation of fractional phi-values (mean of ~0.3) has resisted a precise and consistent interpretation. Here, we describe how to construct, parameterize, and employ a simple one-dimensional free energy surface model that is grounded in the basic tenets of the energy landscape theory to predict and simulate the effect of mutations on folding kinetics. As a proof of principle, we simulate one-dimensional free energy profiles of 806 mutations from 24 different proteins employing just the experimental destabilization as input, reproduce the relative unfolding activation free energies with a correlation of 0.91, and show that the mean phi-value of 0.3 essentially corresponds to the extent of stabilization energy gained at the barrier top while folding.

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.

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.

Thermodynamics and folding landscapes of large proteins from a statistical mechanical model

Statistical mechanical models that afford an intermediate resolution between macroscopic chemical models and all-atom simulations have been successful in capturing folding behaviors of many small single-domain proteins. However, the applicability of one such successful approach, the Wako-Saitô-Muñoz-Eaton (WSME) model, is limited by the size of the protein as the number of conformations grows exponentially with protein length. In this work, we surmount this size limitation by introducing a novel approximation that treats stretches of 3 or 4 residues as blocks, thus reducing the phase space by nearly three orders of magnitude. The performance of the 'bWSME' model is validated by comparing the predictions for a globular enzyme (RNase H) and a repeat protein (IκBα), against experimental observables and the model without block approximation. Finally, as a proof of concept, we predict the free-energy surface of the 370-residue, multi-domain maltose binding protein and identify an intermediate in good agreement with single-molecule force-spectroscopy measurements. The bWSME model can thus be employed as a quantitative predictive tool to explore the conformational landscapes of large proteins, extract the structural features of putative intermediates, identify parallel folding paths, and thus aid in the interpretation of both ensemble and single-molecule experiments.

Thermally versus Chemically Denatured Protein States

Protein unfolding thermodynamic parameters are conventionally extracted from equilibrium thermal and chemical denaturation experiments. Despite decades of work, the degree of structure and the compactness of denatured states populated in these experiments are still open questions. Here, building on previous works, we show that thermally and chemically denatured protein states are distinct from the viewpoint of far-ultraviolet circular dichroism experiments that report on the local conformational status of peptide bonds. The differences identified are independent of protein length, structural class, or experimental conditions, highlighting the presence of two sub-ensembles within the denatured states. The sub-ensembles, U_T and U_D for thermally induced and denaturant-induced unfolded states, respectively, can exclusively exchange populations as a function of temperature at high chemical denaturant concentrations. Empirical analysis suggests that chemically denatured states are ∼50% more expanded than the thermally denatured chains of the same protein. Our observations hint that the strength of protein backbone-backbone interactions and identity versus backbone-solvent/co-solvent interactions determine the conformational distributions. We discuss the implications for protein folding mechanisms, the heterogeneity in relaxation rates, and folding diffusion coefficients.

Modulation of allosteric coupling by mutations: from protein dynamics and packing to altered native ensembles and function

A large body of work has gone into understanding the effect of mutations on protein structure and function. Conventional treatments have involved quantifying the change in stability, activity and relaxation rates of the mutants with respect to the wild-type protein. However, it is now becoming increasingly apparent that mutational perturbations consistently modulate the packing and dynamics of a significant fraction of protein residues, even those that are located >10–15 Å from the mutated site. Such long-range modulation of protein features can distinctly tune protein stability and the native conformational ensemble contributing to allosteric modulation of function. In this review, I summarize a series of experimental and computational observations that highlight the incredibly pliable nature of proteins and their response to mutational perturbations manifested via the intra-protein interaction network. I highlight how an intimate understanding of mutational effects could pave the way for integrating stability, folding, cooperativity and even allostery within a single physical framework.

Tuning the Continuum of Structural States in the Native Ensemble of a Regulatory Protein

The mesoscale nature of proteins allows for an efficient coupling between environmental cues and conformational changes, enabling their function as molecular transducers. Delineating the precise structural origins of such a connection and the expected spectroscopic response has, however, been challenging. In this work, we perform a combination of urea-temperature double perturbation experiments and theoretical modeling to probe the conformational landscape of Cnu, a natural thermosensor protein. We observe unique ensemble signatures that point to a continuum of conformational substates in the native ensemble and that respond intricately to perturbations upon monitoring secondary and tertiary structures, distances between an intrinsic FRET pair, and hydrodynamic volumes. Binding assays further reveal a weakening of the Cnu functional complex with temperature, highlighting the molecular origins of signal transduction critical for pathogenic response in enterobacteriaceae.

Toward a quantitative description of microscopic pathway heterogeneity in protein folding

Experimentally consistent statistical modeling of protein folding thermodynamics reveals unprecedented complexity with numerous parallel folding routes in five different proteins.

Energetic and topological determinants of a phosphorylation-induced disorder-to-order protein conformational switch

We show that the phosphorylation of 4E-BP2 acts as a triggering event to shape its folding-function landscape that is delicately balanced between conflicting favorable energetics and intrinsically unfavorable topological connectivity. We further provide first evidence that the fitness landscapes of proteins at the threshold of disorder can differ considerably from ordered domains.

Imprints of function on the folding landscape: functional role for an intermediate in a conserved eukaryotic binding protein

In the computational characterization of single domain protein folding, the effective free energies of numerous microstates are projected onto few collective degrees of freedom that in turn serve as well-defined reaction coordinates. In this regard, one-dimensional (1D) free energy profiles are widely used mainly for their simplicity. Since folding and functional landscapes are interlinked, how well can these reduced representations capture the structural and dynamic features of functional states while being simultaneously consistent with experimental observables? We investigate this issue by characterizing the folding of the four-helix bundle bovine acyl-CoA binding protein (bACBP), which exhibits complex equilibrium and kinetic behaviours, employing an Ising-like statistical mechanical model and molecular simulations. We show that the features of the 1D free energy profile are sufficient to quantitatively reproduce multiple experimental observations including millisecond chevron-like kinetics and temperature dependence, a microsecond fast phase, barrier heights, unfolded state movements, the intermediate structure and average ϕ-values. Importantly, we find that the structural features of the native-like intermediate (partial disorder in helix 1) are intricately linked to a unique interplay between packing and electrostatics in this domain. By comparison with available experimental data, we propose that this intermediate determines the promiscuous functional behaviour of bACBP that exhibits broad substrate specificity. Our results present evidence to the possibility of employing the statistical mechanical model and the resulting 1D free energy profile to not just understand folding mechanisms but to even extract features of functionally relevant states and their energetic origins.

Human Frataxin Folds Via an Intermediate State. Role of the C-Terminal Region

The aim of this study is to investigate the folding reaction of human frataxin, whose deficiency causes the neurodegenerative disease Friedreich's Ataxia (FRDA). The characterization of different conformational states would provide knowledge about how frataxin can be stabilized without altering its functionality. Wild-type human frataxin and a set of mutants, including two highly destabilized FRDA-associated variants were studied by urea-induced folding/unfolding in a rapid mixing device and followed by circular dichroism. The analysis clearly indicates the existence of an intermediate state (I) in the folding route with significant secondary structure content but relatively low compactness, compared with the native ensemble. However, at high NaCl concentrations I-state gains substantial compaction, and the unfolding barrier is strongly affected, revealing the importance of electrostatics in the folding mechanism. The role of the C-terminal region (CTR), the key determinant of frataxin stability, was also studied. Simulations consistently with experiments revealed that this stretch is essentially unstructured, in the most compact transition state ensemble (TSE2). The complete truncation of the CTR drastically destabilizes the native state without altering TSE2. Results presented here shed light on the folding mechanism of frataxin, opening the possibility of mutating it to generate hyperstable variants without altering their folding kinetics.

Using the folding landscapes of proteins to understand protein function

Proteins fold on a biologically-relevant timescale because of a funnel-shaped energy landscape. This landscape is sculpted through evolution by selecting amino-acid sequences that stabilize native interactions while suppressing stable non-native interactions that occur during folding. However, there is strong evolutionary selection for functional residues and these cannot be chosen to optimize folding. Their presence impacts the folding energy landscape in a variety of ways. Here, we survey the effects of functional residues on folding by providing several examples. We then review how such effects can be detected computationally and be used as assays for protein function. Overall, an understanding of how functional residues modulate folding should provide insights into the design of natural proteins and their homeostasis.