Understanding the frustration arising from the competition between function, misfolding, and aggregation in a globular protein

Biomolecular Condensates Can Enhance Pathological RNA Clustering

Intracellular aggregation of repeat expanded RNA has been implicated in many neurological disorders. Here, we study the role of biomolecular condensates on irreversible RNA clustering. We find that physiologically relevant and disease-associated repeat RNAs spontaneously undergo an age-dependent percolation transition inside multi-component protein-nucleic acid condensates to form nanoscale clusters. Homotypic RNA clusters drive the emergence of multiphasic condensate structures with an RNA-rich solid core surrounded by an RNA-depleted fluid shell. The timescale of the RNA clustering, which drives a liquid-to-solid transition of biomolecular condensates, is determined by the sequence features, stability of RNA secondary structure, and repeat length. Importantly, G3BP1, the core scaffold of stress granules, introduces heterotypic buffering to homotypic RNA-RNA interactions and impedes intra-condensate RNA clustering in an ATP-independent manner. Our work suggests that biomolecular condensates can act as sites for RNA aggregation. It also highlights the functional role of RNA-binding proteins in suppressing aberrant RNA phase transitions.

The Mechanism of Folding of Human Frataxin in Comparison to the Yeast Homologue - Broad Energy Barriers and the General Properties of the Transition State

The funneled energy landscape theory suggests that the folding pathway of homologous proteins should converge at the late stages of folding. In this respect, proteins displaying a broad energy landscape for folding are particularly instructive, allowing inferring both the early, intermediate and late stages of folding. In this paper we explore the folding mechanisms of human frataxin, an essential mitochondrial protein linked to the neurodegenerative disorder Friedreich's ataxia. Building upon previous studies on the yeast homologue, the folding pathway of human frataxin is thoroughly examined, revealing a mechanism implying the presence of a broad energy barrier, reminiscent of the yeast counterpart. Through an extensive site-directed mutagenesis, we employed a Φ -value analysis to map native-like contacts in the folding transition state. The presence of a broad energy barrier facilitated the exploration of such contacts in both early and late folding events. We compared results from yeast and human frataxin providing insights into the impact of native topology on the folding mechanism and elucidating the properties of the underlying free energy landscape. The findings are discussed in the context of the funneled energy landscape theory of protein folding.

In silico design of misfolding resistant proteins: the role of structural similarity of a competing conformational ensemble in the optimization of frustration

Most state-of-the-art in silico design methods fail due to misfolding of designed sequences to a conformation other than the target. Thus, a method to design misfolding resistant proteins will provide a better understanding of the misfolding phenomenon and will also increase the success rate of in silico design methods. In this work, we optimize the conformational ensemble to be selected for negative design purposes based on the similarity of the conformational ensemble to the target. Five ensembles with different degrees of similarity to the target are created and destabilized and the target is stabilized while designing sequences using mean field theory and Monte Carlo simulation methods. The results suggest that the degree of similarity of the non-native conformations to the target plays a prominent role in designing misfolding resistant protein sequences. The design procedures that destabilize the conformational ensemble with moderate similarity to the target have proven to be more promising. Incorporation of either highly similar or highly dissimilar conformations to the target conformation into the non-native ensemble to be destabilized may lead to sequences with a higher misfolding propensity. This will significantly reduce the conformational space to be considered in any protein design procedure. Interestingly, the results suggest that a sequence with higher frustration in the target structure does not necessarily lead to a misfold prone sequence. A successful design method may purposefully choose a frustrated sequence in the target conformation if that sequence is even more frustrated in the competing non-native conformations.

Computational insights into dynamics and conformational stability of N-acetylmannosamine kinase mutations

The activity of UDP-GlcNAc 2-epimerase/ManNAc kinase (GNE) is essential for the biosynthesis of sialic acid, which is involved in cellular processes in health and diseases. GNE contains an N-terminal epimerase domain and a C-terminal kinase domain (N-acetylmannosamine kinase, MNK). Mutations of the GNE protein led to hypoactivity of the enzyme and cause sialurea or autosomal recessive inclusion body myopathy/Nonaka myopathy. Here, we used all-atom molecular dynamics (MD) simulations to comprehend the folding, dynamics and conformational stability of MNK variants, including the wild type (WT) and three mutants (H677R, V696M and H677R/V696M). The deleterious and destabilizing nature of MNK mutants were predicted using different prediction tools. Results predicted that mutations modulate the stability, flexibility and function of MNK. The effect of mutations on the conformational stability and dynamics of MNK was next studied through the free-energy landscape (FEL), hydrogen-bonds and secondary structure changes. The FEL results show that the mutations interfere with various conformational transitions in both WT and mutants, exposing the structural underpinnings of protein destabilization and unfolding brought on by mutation. We discover that, when compared to the other two mutations, V696M and H677R/V696M, H677R has the most harmful effects. These findings have a strong correlation with published experimental studies that demonstrate how these mutations disrupt MNK activity. Hence, this computational study describes the structural details to unravel the mutant effects at the atomistic resolution and has implications for understanding the GNE's physiological and pathological role.Communicated by Ramaswamy H. Sarma.

Automated optimisation of solubility and conformational stability of antibodies and proteins

Biologics, such as antibodies and enzymes, are crucial in research, biotechnology, diagnostics, and therapeutics. Often, biologics with suitable functionality are discovered, but their development is impeded by developability issues. Stability and solubility are key biophysical traits underpinning developability potential, as they determine aggregation, correlate with production yield and poly-specificity, and are essential to access parenteral and oral delivery. While advances for the optimisation of individual traits have been made, the co-optimization of multiple traits remains highly problematic and time-consuming, as mutations that improve one property often negatively impact others. In this work, we introduce a fully automated computational strategy for the simultaneous optimisation of conformational stability and solubility, which we experimentally validate on six antibodies, including two approved therapeutics. Our results on 42 designs demonstrate that the computational procedure is highly effective at improving developability potential, while not affecting antigen-binding. We make the method available as a webserver at www-cohsoftware.ch.cam.ac.uk.

An Outlook on the Complexity of Protein Morphogenesis in Health and Disease

The study of the mechanisms whereby proteins achieve their native functionally competent conformation has been a key issue in molecular biosciences over the last 6 decades. Nevertheless, there are several debated issues and open problems concerning some aspects of this fundamental problem. By considering the emerging complexity of the so-called “native state,” we attempt hereby to propose a personal account on some of the key topics in the field, ranging from the relationships between misfolding and diseases to the significance of protein disorder. Finally, we briefly describe the recent and exciting advances in predicting protein structures from their amino acid sequence.

Characterization of early and late transition states of the folding pathway of a SH2 domain

Albeit SH2 domains are abundant protein-protein interaction modules with fundamental roles in the regulation of several physiological and molecular pathways in the cell, the available information about the determinants of their thermodynamic stability and folding properties are still very limited. In this work, we provide a quantitative characterization of the folding pathway of the C-terminal SH2 domain of SHP2, conducted through a combination of site-directed mutagenesis and kinetic (un)folding experiments (Φ-value analysis). The energetic profile of the folding reaction of the C-SH2 domain is described by a three-state mechanism characterized by the presence of two transition states and a high-energy intermediate. The production of 29 site-directed variants allowed us to calculate the degree of native-like interactions occurring in the early and late events of the folding reaction. Data analysis highlights the presence of a hydrophobic folding nucleus surrounded by a lower degree of structure in the early events of folding, further consolidated as the reaction proceeds towards the native state. Interestingly, residues physically located in the functional region of the domain reported unusual Φ-values, a hallmark of the presence of transient misfolding. We compared our results with previous ones obtained for the N-terminal SH2 domain of SHP2. Notably, a conserved complex folding mechanism implying the presence of a folding intermediate arise from comparison, and the relative stability of such intermediate appears to be highly sequence dependent. Data are discussed under the light of previous works on SH2 domains.

Unveiling induced folding of intrinsically disordered proteins - Protein engineering, frustration and emerging themes

Intrinsically disordered proteins (IDPs) can be generally described as a class of proteins that lack a well-defined ordered structure in isolation at physiological conditions. Upon binding to their physiological ligands, IDPs typically undergo a disorder-to-order transition, which may or may not lead to the complete folding of the IDP. In this short review, we focus on some of the key findings pertaining to the mechanisms of such induced folding. In particular, first we describe the general features of the reaction; then, we discuss some of the most remarkable findings obtained from applying protein engineering in synergy with kinetic studies to induced folding; and finally, we offer a critical view on some of the emerging themes when considering the structural heterogeneity of IDPs vis-à-vis to their inherent frustration.

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.

Dissecting Inter-domain Cooperativity in the Folding of a Multi Domain Protein

Correct protein folding underlies all cellular functions. While there are detailed descriptions and a good understanding of protein folding pathways for single globular domains there is a paucity of quantitative data regarding folding of multidomain proteins. We have here investigated the folding of a three-domain supramodule from the protein PSD-95, consisting of one PDZ domain, one SH3 domain and one guanylate kinase-like (GK) domain. This supramodule has previously been shown to work as one functional unit with regard to ligand binding. We used equilibrium and kinetic folding experiments to demonstrate that the PDZ domain folds faster and independently from the SH3-GK tandem, which folds as one cooperative unit. However, concurrent folding of the PDZ domain slows down folding of SH3-GK by non-native interactions, resulting in an off-pathway folding intermediate. Our data contribute to an emerging description of multidomain protein folding in which individual domains cannot a priori be viewed as separate folding units.

Probing the Effects of Local Frustration in the Folding of a Multidomain Protein

Our current knowledge of protein folding is primarily based on experimental data obtained from isolated domains. In fact, because of their complexity, multidomain proteins have been elusive to the experimental characterization. Thus, the folding of a domain in isolation is generally assumed to resemble what should be observed for more complex structural architectures. Here we compared the folding mechanism of a protein domain in isolation and in the context of its supramodular multidomain structure. By carrying out an extensive mutational analysis we illustrate that while the early events of folding are malleable and influenced by the absence/presence of the neighboring structures, the late events appear to be more robust. These effects may be explained by analyzing the local frustration patterns of the domain, providing critical support for the funneled energy landscape theory of protein folding, and highlighting the role of protein frustration in sculpting the early events of the reaction.

Converging experimental and computational views of the knotting mechanism of a small knotted protein

MJ0366 from Methanocaldococcus jannaschii is the smallest topologically knotted protein known to date. 92 residues in length, MJ0366 ties a trefoil (31) knot by threading its C-terminal helix through a buttonhole formed by the remainder of the secondary structure elements. By generating a library of point mutations at positions pertinent to the knot formation, we systematically evaluated the contributions of individual residues to the folding stability and kinetics of MJ0366. The experimental Φ-values were used as restraints to computationally generate an ensemble of conformations that correspond to the transition state of MJ0366, which revealed several nonnative contacts. The importance of these nonnative contacts in stabilizing the transition state of MJ0366 was confirmed by a second round of mutagenesis, which also established the pivotal role of F15 in stapling the network of hydrophobic interactions around the threading C-terminal helix. Our converging experimental and computational results show that, despite the small size, the transition state of MJ0366 is formed at a very late stage of the folding reaction coordinate, following a polarized pathway. Eventually, the formation of extensive native contacts, as well as a number of nonnative ones, leads to the threading of the C-terminal helix that defines the topological knot.

Fuzziness and Frustration in the Energy Landscape of Protein Folding, Function, and Assembly

Are all protein interactions fully optimized? Do suboptimal interactions compromise specificity? What is the functional impact of frustration? Why does evolution not optimize some contacts? Proteins and their complexes are best described as ensembles of states populating an energy landscape. These ensembles vary in breadth from narrow ensembles clustered around a single average X-ray structure to broader ensembles encompassing a few different functional "taxonomic" states on to near continua of rapidly interconverting conformations, which are called "fuzzy" or even "intrinsically disordered". Here we aim to provide a comprehensive framework for confronting the structural and dynamical continuum of protein assemblies by combining the concepts of energetic frustration and interaction fuzziness. The diversity of the protein structural ensemble arises from the frustrated conflicts between the interactions that create the energy landscape. When frustration is minimal after folding, it results in a narrow ensemble, but residual frustrated interactions result in fuzzy ensembles, and this fuzziness allows a versatile repertoire of biological interactions. Here we discuss how fuzziness and frustration play off each other as proteins fold and assemble, viewing their significance from energetic, functional, and evolutionary perspectives.We demonstrate, in particular, that the common physical origin of both concepts is related to the ruggedness of the energy landscapes, intramolecular in the case of frustration and intermolecular in the case of fuzziness. Within this framework, we show that alternative sets of suboptimal contacts may encode specificity without achieving a single structural optimum. Thus, we demonstrate that structured complexes may not be optimized, and energetic frustration is realized via different sets of contacts leading to multiplicity of specific complexes. Furthermore, we propose that these suboptimal, frustrated, or fuzzy interactions are under evolutionary selection and expand the biological repertoire by providing a multiplicity of biological activities. In accord, we show that non-native interactions in folding or interaction landscapes can cooperate to generate diverse functional states, which are essential to facilitate adaptation to different cellular conditions. Thus, we propose that not fully optimized structures may actually be beneficial for biological activities of proteins via an alternative set of suboptimal interactions. The importance of such variability has not been recognized across different areas of biology.This account provides a modern view on folding, function, and assembly across the protein universe. The physical framework presented here is applicable to the structure and dynamics continuum of proteins and opens up new perspectives for drug design involving not fully structured, highly dynamic protein assemblies.

Studying GGDEF Domain in the Act: Minimize Conformational Frustration to Prevent Artefacts

GGDEF-containing proteins respond to different environmental cues to finely modulate cyclic diguanylate (c-di-GMP) levels in time and space, making the allosteric control a distinctive trait of the corresponding proteins. The diguanylate cyclase mechanism is emblematic of this control: two GGDEF domains, each binding one GTP molecule, must dimerize to enter catalysis and yield c-di-GMP. The need for dimerization makes the GGDEF domain an ideal conformational switch in multidomain proteins. A re-evaluation of the kinetic profile of previously characterized GGDEF domains indicated that they are also able to convert GTP to GMP: this unexpected reactivity occurs when conformational issues hamper the cyclase activity. These results create new questions regarding the characterization and engineering of these proteins for in solution or structural studies.

The Effect of Proline cis-trans Isomerization on the Folding of the C-Terminal SH2 Domain from p85

SH2 domains are protein domains that modulate protein-protein interactions through a specific interaction with sequences containing phosphorylated tyrosines. In this work, we analyze the folding pathway of the C-terminal SH2 domain of the p85 regulatory subunit of the protein PI3K, which presents a proline residue in a cis configuration in the loop between the βE and βF strands. By employing single and double jump folding and unfolding experiments, we demonstrate the presence of an on-pathway intermediate that transiently accumulates during (un)folding. By comparing the kinetics of folding of the wild-type protein to that of a site-directed variant of C-SH2 in which the proline was replaced with an alanine, we demonstrate that this intermediate is dictated by the peptidyl prolyl cis-trans isomerization. The results are discussed in the light of previous work on the effect of peptidyl prolyl cis-trans isomerization on folding events.

Folding and Dynamics Are Strongly pH-Dependent in a Psychrophile Frataxin

Protein dynamics, folding, and thermodynamics represent a central aspect of biophysical chemistry. pH, temperature, and denaturant perturbations inform our understanding of diverse contributors to stability and rates. In this work, we performed a thermodynamic analysis using a combined experimental and computational approach to gain insights into the role of electrostatics in the folding reaction of a psychrophile frataxin variant from Psychromonas ingrahamii. This folding reaction is strongly modulated by pH with a single, narrow, and well-defined transition state with ∼80% compactness, ∼70% electrostatic interactions, and ∼60% hydration shell compared to the native state (αD = 0.82, αH = 0.67, and αΔCp = 0.59). Our results are best explained by a two-proton/two-state model with very different pKa values of the native and denatured states (∼5.5 and ∼8.0, respectively). As a consequence, the stability strongly increases from pH 8.0 to 6.0 (|ΔΔG°| = 5.2 kcal mol-1), mainly because of a decrease in the TΔS°. Variation of ΔH° and ΔS° at pH below 7.0 is dominated by a change in ΔHf⧧ and ΔSf⧧, while at pH above 7.0, it is governed by ΔHu⧧ and ΔSu⧧. Molecular dynamics simulations showed that these pH modulations could be explained by the fluctuations of two regions, rich in electrostatic contacts, whose dynamics are pH-dependent and motions are strongly correlated. Results presented herein contribute to the understanding of the stability and dynamics of this frataxin variant, pointing to an intrinsic feature of the family topology to support different folding mechanisms.

Diverse Structural Conversion and Aggregation Pathways of Alzheimer's Amyloid-β (1-40)

Complex amyloid aggregation of amyloid-β (1-40) (Aβ1-40) in terms of monomer structures has not been fully understood. Herein, we report the microscopic mechanism and pathways of Aβ1-40 aggregation with macroscopic viewpoints through tuning its initial structure and solubility. Partial helical structures of Aβ1-40 induced by low solvent polarity accelerated cytotoxic Aβ1-40 amyloid fibrillation, while predominantly helical folds did not aggregate. Changes in the solvent polarity caused a rapid formation of β-structure-rich protofibrils or oligomers via aggregation-prone helical structures. Modulation of the pH and salt concentration transformed oligomers to protofibrils, which proceeded to amyloid formation. We reveal diverse molecular mechanisms underlying Aβ1-40 aggregation with conceptual energy diagrams and propose that aggregation-prone partial helical structures are key to inducing amyloidogenesis. We demonstrate that context-dependent protein aggregation is comprehensively understood using the macroscopic phase diagram, which provides general insights into differentiation of amyloid formation and phase separation from unfolded and folded structures.

Localizing Frustration in Proteins Using All-Atom Energy Functions

The problems of protein folding and protein design are two sides of the same coin. Protein folding involves exploring a protein's configuration space given a fixed sequence, whereas protein design involves searching in sequence space given a particular target structure. For a protein to fold quickly and reliably, its energy landscape must be biased toward the folded ensemble throughout its configuration space and must lack deep kinetic traps that would otherwise frustrate folding. Evolution has "designed" the sequences of many naturally occurring proteins, through an eons-long process of random mutation and selection, to yield landscapes with a minimal degree of frustration. The task facing humans hoping to design protein sequences that fold into particular structures is to use the available approximate energy functions to sculpt funneled landscapes that work in the laboratory. In this work, we demonstrate how to calculate several localized frustration measures using an all-atom energy function. Specifically, we employ the Rosetta energy function, which has been used successfully to design proteins and which has a natural pairwise decomposition that is suitably solvent-averaged. We calculate these newly developed frustration measures for both a mutated WW domain, FiP35, and a three-helix bundle that was designed completely by humans, Alpha3D. The structure of FiP35 exhibits less localized frustration than that of Alpha3D. A mutation toward the consensus sequence for WW domains in FiP35, which has been shown unexpectedly in experiment to disrupt folding, induces localized frustration by disrupting the hydrophobic core. By performing a limited redesign on the sequence of Alpha3D, we show that some, but not all, mutations that lower the energy also result in decreased frustration. The results suggest that, in addition to being useful for detecting residual frustration in protein structures, optimizing the localized frustration measures presented here may be a useful and automatic means of balancing positive and negative design in protein design tasks.

Third generation antibody discovery methods: in silico rational design

Owing to their outstanding performances in molecular recognition, antibodies are extensively used in research and applications in molecular biology, biotechnology and medicine. Recent advances in experimental and computational methods are making it possible to complement well-established in vivo (first generation) and in vitro (second generation) methods of antibody discovery with novel in silico (third generation) approaches. Here we describe the principles of computational antibody design and review the state of the art in this field. We then present Modular, a method that implements the rational design of antibodies in a modular manner, and describe the opportunities offered by this approach.

A structurally heterogeneous transition state underlies coupled binding and folding of disordered proteins

Many intrinsically disordered proteins (IDPs) attain a well-defined structure in a coupled folding and binding reaction with another protein. Such reactions may involve early to late formation of different native structural regions along the reaction pathway. To obtain insights into the transition state for a coupled binding and folding reaction, we performed restrained molecular dynamics simulations using previously determined experimental binding Φ_b values of the interaction between two IDP domains: the activation domain from the p160 transcriptional co-activator for thyroid hormone and retinoid receptors (ACTR) and the nuclear co-activator binding domain (NCBD) of CREB-binding protein, each forming three well-defined α-helices upon binding. These simulations revealed that both proteins are largely disordered in the transition state for complex formation, except for two helices, one from each domain, that display a native-like structure. The overall transition state structure was extended and largely dynamic with many weakly populated contacts. To test the transition state model, we combined site-directed mutagenesis with kinetic experiments, yielding results consistent with overall diffuse interactions and formation of native intramolecular interactions in the third NCBD helix during the binding reaction. Our findings support the view that the transition state and, by inference, any encounter complex in coupled binding and folding reactions are structurally heterogeneous and largely independent of specific interactions. Furthermore, experimental Φ_b values and Brønsted plots suggested that the transition state is globally robust with respect to most mutations but can display more native-like features for some highly destabilizing mutations, possibly because of Hammond behavior or ground-state effects.

Elucidation of the structural stability and dynamics of heterogeneous intermediate ensembles in unfolding pathway of the N-terminal domain of TDP-43

The N-terminal domain of the RNA binding protein TDP-43 (NTD) is essential to both physiology and proteinopathy; however, elucidation of its folding/unfolding still remains a major quest. In this study, we have investigated the biophysical behavior of intermediate ensembles employing all-atom molecular dynamics simulations in 8 M urea accelerated with high temperatures to achieve unfolded states in a confined computation time. The cumulative results of the 2.75 μs simulations show that unfolding of the NTD at 350 K evolves through different stable and meta-stable intermediate states. The free-energy landscape reveals two meta-stable intermediates (IN and IU) stabilized by non-native interactions, which are largely hydrophilic and highly energetically frustrated. A single buried tryptophan residue, W80, undergoes solvent exposure to different extents during unfolding; this suggests a structurally heterogeneous population of intermediate ensembles. Furthermore, the structure properties of the IN state show a resemblance to the molten globule (MG) state with most of the secondary structures intact. The unfolding of the NTD is initiated by the loss of β-strands, and the unfolded (U) states exhibit a population of non-native α-helices. These non-native unfolded intermediate ensembles may mediate protein oligomerization, leading to the formation of pathological, irreversible aggregates, characteristics of disease pathogenesis.

A method for partitioning the information contained in a protein sequence between its structure and function

Proteins employ the information stored in the genetic code and translated into their sequences to carry out well-defined functions in the cellular environment. The possibility to encode for such functions is controlled by the balance between the amount of information supplied by the sequence and that left after that the protein has folded into its structure. We study the amount of information necessary to specify the protein structure, providing an estimate that keeps into account the thermodynamic properties of protein folding. We thus show that the information remaining in the protein sequence after encoding for its structure (the 'information gap') is very close to what needed to encode for its function and interactions. Then, by predicting the information gap directly from the protein sequence, we show that it may be possible to use these insights from information theory to discriminate between ordered and disordered proteins, to identify unknown functions, and to optimize artificially-designed protein sequences.

Stability of an aggregation-prone partially folded state of human profilin-1 correlates with aggregation propensity

A set of missense mutations in the gene encoding profilin-1 has been linked to the onset of familial forms of ALS (fALS), also known as Lou Gehrig's disease. The pathogenic potential of these mutations is linked to the formation of intracellular inclusions of the mutant proteins and correlates with the mutation-induced destabilization of its native, fully folded state. However, the mechanism by which these mutations promote misfolding and self-assembly is yet unclear. Here, using temperature-jump and stopped-flow kinetic measurements, we show that, during refolding, WT profilin-1 transiently populates a partially folded (PF) state endowed with hydrophobic clusters exposed to the solvent and with no detectable secondary structure. We observed that this conformational state is marginally stable at neutral pH but becomes significantly populated at mildly acidic pH. Interestingly, the fALS-associated mutations did not cause a change in the refolding mechanism of profilin-1, but induced a stabilization of the PF state. In the presence of preformed profilin-1 aggregates, the PF state, unlike the unfolded and folded states, could interact with these aggregates via nonspecific hydrophobic interactions and also increase thioflavin-T fluorescence, revealing its amyloidogenic potential. Moreover, in the variants tested, we found a correlation between conformational stability of PF and aggregation propensity, defining this conformational state as an aggregation-prone folding intermediate. In conclusion, our findings indicate that mutation-induced stabilization of a partially folded state can enhance profilin-1 aggregation and thereby contribute to the pathogenicity of the mutations.

Deciphering the dark proteome of Chikungunya virus

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus. The outbreak of CHIKV infection has been seen in many tropical and subtropical regions of the biosphere. Current reports evidenced that after outbreaks in 2005-06, the fitness of this virus propagating in Aedes albopictus enhanced due to the epistatic mutational changes in its envelope protein. In our study, we evaluated the prevalence of intrinsically disordered proteins (IDPs) and IDP regions (IDPRs) in CHIKV proteome. IDPs/IDPRs are known as members of a 'Dark Proteome' that defined as a set of polypeptide segments or whole protein without unique three-dimensional structure within the cellular milieu but with significant biological functions, such as cell cycle regulation, control of signaling pathways, and maintenance of viral proteomes. However, the intrinsically disordered aspects of CHIKV proteome and roles of IDPs/IDPRs in the pathogenic mechanism of this important virus have not been evaluated as of yet. There are no existing reports on the analysis of intrinsic disorder status of CHIKV. To fulfil this goal, we have analyzed the abundance and functionality of IDPs/IDPRs in CHIKV proteins, involved in the replication and maturation. It is likely that these IDPs/IDPRs can serve as novel targets for disorder based drug design.

Frustration, function and folding

Natural protein molecules are exceptional polymers. Encoded in apparently random strings of amino-acids, these objects perform clear physical tasks that are rare to find by simple chance. Accurate folding, specific binding, powerful catalysis, are examples of basic chemical activities that the great majority of polypeptides do not display, and are thought to be the outcome of the natural history of proteins. Function, a concept genuine to Biology, is at the core of evolution and often conflicts with the physical constraints. Locating the frustration between discrepant goals in a recurrent system leads to fundamental insights about the chances and necessities that shape the encoding of biological information.

The Folding Pathway of the KIX Domain

The KIX domain is an 89-residues globular domain with an important role in mediating protein-protein interactions. The presence of two distinct binding sites in such a small domain makes KIX a suitable candidate to investigate the effect of the potentially divergent demands between folding and function. Here, we report an extensive mutational analysis of the folding pathway of the KIX domain, based on 30 site-directed mutants, which allow us to assess the structures of both the transition and denatured states. Data reveal that, while the transition state presents mostly native-like interactions, the denatured state is somewhat misfolded. We mapped some of the non-native contacts in the denatured state using a second round of mutagenesis, based on double mutant cycles on 15 double mutants. Interestingly, such a misfolding arises from non-native interactions involving the residues critical for the function of the protein. The results described in this work appear to highlight the diverging demands between folding and function that may lead to misfolding, which may be observed in the early stages of folding.

Kinetics and polymorphs of yeast prion Sup35NM amyloidogenesis

Amyloidogenic proteins often form many types of aggregates, which are a critical determinant of cytotoxicity and tissue specificity. However, the molecular mechanisms underlying the generation of distinct amyloids and their influence on cells remain largely unknown. We herein investigated the polymorphic amyloid formation of the yeast prion protein, Sup35NM, an intrinsically disordered N-terminal fragment of Sup35, under various conditions and its potential relationship to cytotoxicity. Sup35NM aggregated to amyloid fibrils with distinct kinetics, structures, morphologies, tinctorial properties, and conformational stabilities depending on the concentration of NaCl, pH, and temperature, indicating the polymorphic amyloidogenesis of Sup35NM. Detailed kinetic analyses of Sup35NM amyloid formation revealed a strong inverse correlation between the lag time and elongation rate without a correlation between kinetic and structural parameters. These results suggest that kinetic polymorphisms due to distinct nucleation and elongation rates result in structural polymorphs of amyloid fibrils, and also that conditions that enhance or inhibit the nucleation of Sup35NM promote or delay fibril growth. The deleterious effects of polymorphic Sup35NM amyloid fibrils on membrane integrity and cell vitality were minimal. We hypothesize that the innocuous polymorphic nature of Sup35NM amyloid fibrils may be beneficial for gaining time for prion infection prior to cell death.

Smoothing a rugged protein folding landscape by sequence-based redesign

The rugged folding landscapes of functional proteins puts them at risk of misfolding and aggregation. Serine protease inhibitors, or serpins, are paradigms for this delicate balance between function and misfolding. Serpins exist in a metastable state that undergoes a major conformational change in order to inhibit proteases. However, conformational labiality of the native serpin fold renders them susceptible to misfolding, which underlies misfolding diseases such as α₁-antitrypsin deficiency. To investigate how serpins balance function and folding, we used consensus design to create conserpin, a synthetic serpin that folds reversibly, is functional, thermostable, and polymerization resistant. Characterization of its structure, folding and dynamics suggest that consensus design has remodeled the folding landscape to reconcile competing requirements for stability and function. This approach may offer general benefits for engineering functional proteins that have risky folding landscapes, including the removal of aggregation-prone intermediates, and modifying scaffolds for use as protein therapeutics.

Protein Domain-Swapping Can Be a Consequence of Functional Residues

Monomer topology has been implicated in domain-swapping, a potential first step on the route to disease-causing protein aggregation. Despite having the same topology (β1-α1-β2-β3-β4-β5), the cysteine protease inhibitor stefin-B domain swaps more readily than a single-chain variant of the heterodimeric sweet protein monellin (scMn). Here, we computationally study the folding of stefin-B and scMn in order to understand the molecular basis for the difference in their domain-swapping propensities. In agreement with experiments, our structure-based simulations show that scMn folds cooperatively without the population of an intermediate while stefin-B populates an equilibrium intermediate state. Since the simulation intermediate has only one domain structured (β3-β4-β5), it can directly lead to domain-swapping. Using computational variants of stefin-B, we show that the population of this intermediate is caused by regions of stefin-B that have been implicated in protease inhibition. We also find that the protease-binding regions are located on two structural elements and localized in space. In contrast, the residues that contribute to the sweetness of monellin are not localized to a few structural elements but are distributed over the protein fold. We conclude that the distributed functional residues of monellin do not induce large local perturbations in the protein structure, eliminating the formation of folding intermediates and in turn domain-swapping. On the other hand, the localized protease-binding regions of stefin-B promote the formation of a folding intermediate which can lead to domain-swapping. Thus, domain-swapping can be a direct consequence of the constraints that function imposes on the protein structure.

Free energy surface of an intrinsically disordered protein: comparison between temperature replica exchange molecular dynamics and bias-exchange metadynamics

Intrinsically disordered proteins (IDPs), which are expected to be largely unstructured under physiological conditions, make up a large fraction of eukaryotic proteins. Molecular dynamics simulations have been utilized to probe structural characteristics of these proteins, which are not always easily accessible to experiments. However, exploration of the conformational space by brute force molecular dynamics simulations is often limited by short time scales. Present literature provides a number of enhanced sampling methods to explore protein conformational space in molecular simulations more efficiently. In this work, we present a comparison of two enhanced sampling methods: temperature replica exchange molecular dynamics and bias exchange metadynamics. By investigating both the free energy landscape as a function of pertinent order parameters and the per-residue secondary structures of an IDP, namely, human islet amyloid polypeptide, we found that the two methods yield similar results as expected. We also highlight the practical difference between the two methods by describing the path that we followed to obtain both sets of data.

Consensus protein design

A popular and successful strategy in semi-rational design of protein stability is the use of evolutionary information encapsulated in homologous protein sequences. Consensus design is based on the hypothesis that at a given position, the respective consensus amino acid contributes more than average to the stability of the protein than non-conserved amino acids. Here, we review the consensus design approach, its theoretical underpinnings, successes, limitations and challenges, as well as providing a detailed guide to its application in protein engineering.

Rational design of mutations that change the aggregation rate of a protein while maintaining its native structure and stability

A wide range of human diseases is associated with mutations that, destabilizing proteins native state, promote their aggregation. However, the mechanisms leading from folded to aggregated states are still incompletely understood. To investigate these mechanisms, we used a combination of NMR spectroscopy and molecular dynamics simulations to compare the native state dynamics of Beta-2 microglobulin (β2m), whose aggregation is associated with dialysis-related amyloidosis, and its aggregation-resistant mutant W60G. Our results indicate that W60G low aggregation propensity can be explained, beyond its higher stability, by an increased average protection of the aggregation-prone residues at its surface. To validate these findings, we designed β2m variants that alter the aggregation-prone exposed surface of wild-type and W60G β2m modifying their aggregation propensity. These results allowed us to pinpoint the role of dynamics in β2m aggregation and to provide a new strategy to tune protein aggregation by modulating the exposure of aggregation-prone residues.

Protein Frustratometer 2: a tool to localize energetic frustration in protein molecules, now with electrostatics

The protein frustratometer is an energy landscape theory-inspired algorithm that aims at localizing and quantifying the energetic frustration present in protein molecules. Frustration is a useful concept for analyzing proteins’ biological behavior. It compares the energy distributions of the native state with respect to structural decoys. The network of minimally frustrated interactions encompasses the folding core of the molecule. Sites of high local frustration often correlate with functional regions such as binding sites and regions involved in allosteric transitions. We present here an upgraded version of a webserver that measures local frustration. The new implementation that allows the inclusion of electrostatic energy terms, important to the interactions with nucleic acids, is significantly faster than the previous version enabling the analysis of large macromolecular complexes within a user-friendly interface. The webserver is freely available at URL: http://frustratometer.qb.fcen.uba.ar.

On the reliability of peptide nonplanarity seen in ultra-high resolution crystal structures

Ultra-high resolution protein crystal structures have been considered as relatively reliable sources for defining details of protein geometry, such as the extent to which the peptide unit deviates from planarity. Chellapa and Rose (Proteins 2015; 83:1687) recently called this into question, reporting that for a dozen representative protein structures determined at ∼ 1 Å resolution, the diffraction data could be equally well fit with models restrained to have highly planar peptides, i.e. having a standard deviation of the ω torsion angles of only ∼ 1° instead of the typically observed value of ∼ 6°. Here, we document both conceptual and practical shortcomings of that study and show that the more tightly restrained models are demonstrably incorrect and do not fit the diffraction data equally well. We emphasize the importance of inspecting electron density maps when investigating the agreement between a model and its experimental data. Overall, this report reinforces that modern standard refinement protocols have been well-conceived and that ultra-high resolution protein crystal structures, when evaluated carefully and used with an awareness of their levels of coordinate uncertainty, are powerful sources of information for providing reliable information about the details of protein geometry.

High-Resolution Mapping of the Folding Transition State of a WW Domain

Fast-folding WW domains are among the best-characterized systems for comparing experiments and simulations of protein folding. Recent microsecond-resolution experiments and long duration (totaling milliseconds) single-trajectory modeling have shown that even mechanistic changes in folding kinetics due to mutation can now be analyzed. Thus, a comprehensive set of experimental data would be helpful to benchmark the predictions made by simulations. Here, we use T-jump relaxation in conjunction with protein engineering and report mutational Φ-values (Φ(M)) as indicators for folding transition-state structure of 65 side chain, 7 backbone hydrogen bond, and 6 deletion and /or insertion mutants within loop 1 of the 34-residue hPin1 WW domain. Forty-five cross-validated consensus mutants could be identified that provide structural constraints for transition-state structure within all substructures of the WW domain fold (hydrophobic core, loop 1, loop 2, β-sheet). We probe the robustness of the two hydrophobic clusters in the folding transition state, discuss how local backbone disorder in the native-state can lead to non-classical Φ(M)-values (Φ(M) > 1) in the rate-determining loop 1 substructure, and conclusively identify mutations and positions along the sequence that perturb the folding mechanism from loop 1-limited toward loop 2-limited folding.

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.

Dissection of the water cavity of yeast thioredoxin 1: the effect of a hydrophobic residue in the cavity

The water cavity of yeast thioredoxin 1 (yTrx1) is an ancestral, conserved structural element that is poorly understood. We recently demonstrated that the water cavity is involved in the complex protein dynamics that are responsible for the catalytically relevant event of coupling hydration, proton exchange, and motion at the interacting loops. Its main feature is the presence of the conserved polar residue, Asp24, which is buried in a hydrophobic cavity. Here, we evaluated the role of the solvation of Asp24 as the main element that is responsible for the formation of the water cavity in thioredoxins. We showed that the substitution of Asp24 with a hydrophobic residue (D24A) was not sufficient to completely close the cavity. The dynamics of the D24A mutant of yTrx1 at multiple time scales revealed that the D24A mutant presents motions at different time scales near the active site, interaction loops, and water cavity, revealing the existence of a smaller dissected cavity. Molecular dynamics simulation, along with experimental molecular dynamics, allowed a detailed description of the water cavity in wild-type yTrx1 and D24A. The cavity connects the interacting loops, the central β-sheet, and α-helices 2 and 4. It is formed by three contiguous lobes, which we call lobes A-C. Lobe A is hydrophilic and the most superficial. It is formed primarily by the conserved Lys54. Lobe B is the central lobe formed by the catalytically important residues Cys33 and Asp24, which are strategically positioned. Lobe C is the most hydrophobic and is formed by the conserved cis-Pro73. The central lobe B is closed upon introduction of the D24A mutation, revealing that independent forces other than solvation of Asp24 maintain lobes A and C in the open configuration. These data allow us to better understand the properties of this enzyme.

Energy landscapes of functional proteins are inherently risky

Evolutionary pressure for protein function leads to unavoidable sampling of conformational states that are at risk of misfolding and aggregation. The resulting tension between functional requirements and the risk of misfolding and/or aggregation in the evolution of proteins is becoming more and more apparent. One outcome of this tension is sensitivity to mutation, in which only subtle changes in sequence that may be functionally advantageous can tip the delicate balance toward protein aggregation. Similarly, increasing the concentration of aggregation-prone species by reducing the ability to control protein levels or compromising protein folding capacity engenders increased risk of aggregation and disease. In this Perspective, we describe examples that epitomize the tension between protein functional energy landscapes and aggregation risk. Each case illustrates how the energy landscapes for the at-risk proteins are sculpted to enable them to perform their functions and how the risks of aggregation are minimized under cellular conditions using a variety of compensatory mechanisms.