Folding pathways of proteins with increasing degree of sequence identities but different structure and function

Exploring the Effect of Enhanced Sampling on Protein Stability Prediction

Changes in protein stability due to side-chain mutations are evaluated by alchemical free-energy calculations based on classical molecular dynamics (MD) simulations in explicit solvent using the GROMOS force field. Three proteins of different complexity with a total number of 93 single-point mutations are analyzed, and the relative free-energy differences are discussed with respect to configurational sampling and (dis)agreement with experimental data. For the smallest protein studied, a 34-residue WW domain, the starting structure dependence of the alchemical free-energy changes, is discussed in detail. Deviations from previous simulations for the two other proteins are shown to result from insufficient sampling in the earlier studies. Hamiltonian replica exchange in combination with multiple starting structures and sufficient sampling time of more than 100 ns per intermediate alchemical state is required in some cases to reach convergence.

Exposing the distinctive modular behavior of β-strands and α-helices in folded proteins

Although folded proteins are commonly depicted as simplistic combinations of β-strands and α-helices, the actual properties and functions of these secondary-structure elements in their native contexts are just partly understood. The principal reason is that the behavior of individual β- and α-elements is obscured by the global folding cooperativity. In this study, we have circumvented this problem by designing frustrated variants of the mixed α/β-protein S6, which allow the structural behavior of individual β-strands and α-helices to be targeted selectively by stopped-flow kinetics, X-ray crystallography, and solution-state NMR. Essentially, our approach is based on provoking intramolecular "domain swap." The results show that the α- and β-elements have quite different characteristics: The swaps of β-strands proceed via global unfolding, whereas the α-helices are free to swap locally in the native basin. Moreover, the α-helices tend to hybridize and to promote protein association by gliding over to neighboring molecules. This difference in structural behavior follows directly from hydrogen-bonding restrictions and suggests that the protein secondary structure defines not only tertiary geometry, but also maintains control in function and structural evolution. Finally, our alternative approach to protein folding and native-state dynamics presents a generally applicable strategy for in silico design of protein models that are computationally testable in the microsecond-millisecond regime.

Fast and flexible coarse-grained prediction of protein folding routes using ensemble modeling and evolutionary sequence variation

MOTIVATION

Protein folding is a dynamic process through which polypeptide chains reach their native 3D structures. Although the importance of this mechanism is widely acknowledged, very few high-throughput computational methods have been developed to study it.

RESULTS

In this paper, we report a computational platform named P3Fold that combines statistical and evolutionary information for predicting and analyzing protein folding routes. P3Fold uses coarse-grained modeling and efficient combinatorial schemes to predict residue contacts and evaluate the folding routes of a protein sequence within minutes or hours. To facilitate access to this technology, we devise graphical representations and implement an interactive web interface that allows end-users to leverage P3Fold predictions. Finally, we use P3Fold to conduct large and short scale experiments on the human proteome that reveal the broad conservation and variations of structural intermediates within protein families.

AVAILABILITY AND IMPLEMENTATION

A Web server of P3Fold is freely available at http://csb.cs.mcgill.ca/P3Fold.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Effects of Topology and Sequence in Protein Folding Linked via Conformational Fluctuations

Experiments have compared the folding of proteins with different amino acid sequences but the same basic structure, or fold. Results indicate that folding is robust to sequence variations for proteins with some nonlocal folds, such as all-β, whereas the folding of more local, all-α proteins typically exhibits a stronger sequence dependence. Here, we use a coarse-grained model to systematically study how variations in sequence perturb the folding energy landscapes of three model sequences with 3α, 4β + α, and β-barrel folds, respectively. These three proteins exhibit folding features in line with experiments, including expected rank order in the cooperativity of the folding transition and stability-dependent shifts in the location of the free-energy barrier to folding. Using a generalized-ensemble simulation approach, we determine the thermodynamics of around 2000 sequence variants representing all possible hydrophobic or polar single- and double-point mutations. From an analysis of the subset of stability-neutral mutations, we find that folding is perturbed in a topology-dependent manner, with the β-barrel protein being the most robust. Our analysis shows, in particular, that the magnitude of mutational perturbations of the transition state is controlled in part by the size or "width" of the underlying conformational ensemble. This result suggests that the mutational robustness of the folding of the β-barrel protein is underpinned by its conformationally restricted transition state ensemble, revealing a link between sequence and topological effects in protein folding.

Investigating the Molecular Basis of the Aggregation Propensity of the Pathological D76N Mutant of Beta-2 Microglobulin: Role of the Denatured State

Beta-2 microglobulin (β2m) is a protein responsible for a pathologic condition, known as dialysis-related amyloidosis (DRA), caused by its aggregation and subsequent amyloid formation. A naturally occurring mutation of β2m, D76N, presents a higher amyloidogenic propensity compared to the wild type counterpart. Since the three-dimensional structure of the protein is essentially unaffected by the mutation, the increased aggregation propensity of D76N has been generally ascribed to its lower thermodynamic stability and increased dynamics. In this study we compare the equilibrium unfolding and the aggregation propensity of wild type β2m and D76N variant at different experimental conditions. Our data revealed a surprising effect of the D76N mutation in the residual structure of the denatured state, which appears less compact than that of the wild type protein. A careful investigation of the structural malleability of the denatured state of wild type β2m and D76N pinpoint a clear role of the denatured state in triggering the amyloidogenic propensity of the protein. The experimental results are discussed in the light of the previous work on β2m and its role in disease.

Kinetic and thermodynamic studies reveal chemokine homologues CC11 and CC24 with an almost identical tertiary structure have different folding pathways

BACKGROUND

Proteins with low sequence identity but almost identical tertiary structure and function have been valuable to uncover the relationship between sequence, tertiary structure, folding mechanism and functions. Two homologous chemokines, CCL11 and CCL24, with low sequence identity but similar tertiary structure and function, provide an excellent model system for respective studies.

RESULTS

The kinetics and thermodynamics of the two homologous chemokines were systematically characterized. Despite their similar tertiary structures, CCL11 and CCL24 show different thermodynamic stability in guanidine hydrochloride titration, with D_50% = 2.20 M and 4.96 M, respectively. The kinetics curves clearly show two phases in the folding/unfolding processes of both CCL11 and CCL24, which suggests the existence of an intermediate state in their folding/unfolding processes. The folding pathway of both CCL11 and CCL24 could be well described using a folding model with an on-pathway folding intermediate. However, the folding kinetics and stability of the intermediate state of CCL11 and CCL24 are obviously different.

CONCLUSION

Our results suggest homologous proteins with low sequence identity can display almost identical tertiary structure, but very different folding mechanisms, which applies to homologues in the chemokine protein family, extending the general applicability of the above observation.

Protein folding, misfolding and aggregation: The importance of two-electron stabilizing interactions

Proteins associated with neurodegenerative diseases are highly pleiomorphic and may adopt an all-α-helical fold in one environment, assemble into all-β-sheet or collapse into a coil in another, and rapidly polymerize in yet another one via divergent aggregation pathways that yield broad diversity of aggregates’ morphology. A thorough understanding of this behaviour may be necessary to develop a treatment for Alzheimer’s and related disorders. Unfortunately, our present comprehension of folding and misfolding is limited for want of a physicochemical theory of protein secondary and tertiary structure. Here we demonstrate that electronic configuration and hyperconjugation of the peptide amide bonds ought to be taken into account to advance such a theory. To capture the effect of polarization of peptide linkages on conformational and H-bonding propensity of the polypeptide backbone, we introduce a function of shielding tensors of the C^α atoms. Carrying no information about side chain-side chain interactions, this function nonetheless identifies basic features of the secondary and tertiary structure, establishes sequence correlates of the metamorphic and pH-driven equilibria, relates binding affinities and folding rate constants to secondary structure preferences, and manifests common patterns of backbone density distribution in amyloidogenic regions of Alzheimer’s amyloid β and tau, Parkinson’s α-synuclein and prions. Based on those findings, a split-intein like mechanism of molecular recognition is proposed to underlie dimerization of Aβ, tau, αS and PrP^C, and divergent pathways for subsequent association of dimers are outlined; a related mechanism is proposed to underlie formation of PrP^Sc fibrils. The model does account for: (i) structural features of paranuclei, off-pathway oligomers, non-fibrillar aggregates and fibrils; (ii) effects of incubation conditions, point mutations, isoform lengths, small-molecule assembly modulators and chirality of solid-liquid interface on the rate and morphology of aggregation; (iii) fibril-surface catalysis of secondary nucleation; and (iv) self-propagation of infectious strains of mammalian prions.

How cooperative are protein folding and unfolding transitions?

A thermodynamically and kinetically simple picture of protein folding envisages only two states, native (N) and unfolded (U), separated by a single activation free energy barrier, and interconverting by cooperative two-state transitions. The folding/unfolding transitions of many proteins occur, however, in multiple discrete steps associated with the formation of intermediates, which is indicative of reduced cooperativity. Furthermore, much advancement in experimental and computational approaches has demonstrated entirely non-cooperative (gradual) transitions via a continuum of states and a multitude of small energetic barriers between the N and U states of some proteins. These findings have been instrumental towards providing a structural rationale for cooperative versus noncooperative transitions, based on the coupling between interaction networks in proteins. The cooperativity inherent in a folding/unfolding reaction appears to be context dependent, and can be tuned via experimental conditions which change the stabilities of N and U. The evolution of cooperativity in protein folding transitions is linked closely to the evolution of function as well as the aggregation propensity of the protein. A large activation energy barrier in a fully cooperative transition can provide the kinetic control required to prevent the accumulation of partially unfolded forms, which may promote aggregation. Nevertheless, increasing evidence for barrier-less "downhill" folding, as well as for continuous "uphill" unfolding transitions, indicate that gradual non-cooperative processes may be ubiquitous features on the free energy landscape of protein folding.

Molecular dynamics simulations on the conformational transitions from the GA 98 (GA 88) to GB 98 (GB 88) proteins

We performed conventional and targeted molecular dynamics simulations to address the dynamic transition mechanisms of the conformational transitions from the G_A 98 protein with only 1 mutation of Leu45Tyr to G_B 98 and from the G_A 88 protein with 7 mutations of Gly24Ala, Ile25Thr, Ile30Phe, Ile33Tyr, Leu45Tyr, Ile49Thr, and Leu50Lys to G_B 88. The results show that the conformational transition mechanism from the mutated 3α G_A 98 (G_A 88) state to the α+4β G_B 98 (G_B 88) state via several intermediate conformations involves the bending of loops at the N and C termini firstly, the unfolding of αA and αC, then the traversing of αB, and the formation of the 4β layer with the conversion of the hydrophobic core. The bending of loops at the N and C termini and the formation of the crucial transition conformation with the full unfolded structure are key factors in their transition processes. The communication of the interaction network, the bending directions of loops, and the traversing site of αB in the transition of G_A 98 to G_B 98 are markedly different from those in G_A 88 to G_B 88 because of the different mutated residues. The analysis of the correlations and the calculated mass center distances between some segments further supported their conformational transition mechanisms. These results could help people to better understand the Paracelsus challenge. Copyright © 2016 John Wiley & Sons, Ltd.

Theoretical Insights into the Biophysics of Protein Bi-stability and Evolutionary Switches

Deciphering the effects of nonsynonymous mutations on protein structure is central to many areas of biomedical research and is of fundamental importance to the study of molecular evolution. Much of the investigation of protein evolution has focused on mutations that leave a protein’s folded structure essentially unchanged. However, to evolve novel folds of proteins, mutations that lead to large conformational modifications have to be involved. Unraveling the basic biophysics of such mutations is a challenge to theory, especially when only one or two amino acid substitutions cause a large-scale conformational switch. Among the few such mutational switches identified experimentally, the one between the G_A all-α and G_B α+β folds is extensively characterized; but all-atom simulations using fully transferrable potentials have not been able to account for this striking switching behavior. Here we introduce an explicit-chain model that combines structure-based native biases for multiple alternative structures with a general physical atomic force field, and apply this construct to twelve mutants spanning the sequence variation between G_A and G_B. In agreement with experiment, we observe conformational switching from G_A to G_B upon a single L45Y substitution in the GA98 mutant. In line with the latent evolutionary potential concept, our model shows a gradual sequence-dependent change in fold preference in the mutants before this switch. Our analysis also indicates that a sharp G_A/G_B switch may arise from the orientation dependence of aromatic π-interactions. These findings provide physical insights toward rationalizing, predicting and designing evolutionary conformational switches.

The biological functions of globular proteins are intimately related to their folded structures and their associated conformational fluctuations. Evolution of new structures is an important avenue to new functions. Although many mutations do not change the folded state, experiments indicate that a single amino acid substitution can lead to a drastic change in the folded structure. The physics of this switch-like behavior remains to be elucidated. Here we develop a computational model for the relevant physical forces, showing that mutations can lead to new folds by passing through intermediate sequences where the old and new folds occur with varying probabilities. Our approach helps provide a general physical account of conformational switching in evolution and mutational effects on conformational dynamics.

Elucidating the Key Determinants of Structure, Folding, and Stability for the ( 4β+ α ) Conformation of the B1 Domain of Protein G Using Bioinformatics Approaches

The B1 domain of protein G (GB1) is a small, 56 amino acid bacterial immunoglobulin-binding protein with a 4β+ α fold. Architecturally, it is composed of a two-layer sandwich consisting of a four-stranded β -sheet that packs against an α -helix. Using several bioinformatics approaches, we investigated which residues may be key determinants of this fold. We identified nine structurally conserved amino acids using a conservation analysis and propose they are critical to forming and stabilizing the fold. The nine conserved residues form a predominantly hydrophobic nucleus within the core of GB1. A network analysis of all the long-range interactions in the structure of GB1 in concert with a betweenness centrality analysis revealed the relative significance of each conserved amino acid residue based on the number and location of the interactions. This bioinformatics analysis provides an important foundation for the design and interpretation of both computational and experimental work which may be helpful in solving the protein folding problem.

Where the complex things are: single molecule and ensemble spectroscopic investigations of protein folding dynamics

Progress in our understanding of the simple folding dynamics of small proteins and the complex dynamics of large proteins is reviewed. Recent characterizations of the folding transition path of small proteins revealed a simple dynamics explainable by the native centric model. In contrast, the accumulated data showed the substates containing residual structures in the unfolded state and partially populated intermediates, causing complexity in the early folding dynamics of small proteins. The size of the unfolded proteins in the absence of denaturants is likely expanded but still controversial. The steady progress in the observation of folding of large proteins has clarified the rapid formation of long-range contacts that seem inconsistent with the native centric model, suggesting that the folding strategy of large proteins is distinct from that of small proteins.

Protein folding of the SAP domain, a naturally occurring two-helix bundle

•

Thol SAP domain is one of the smallest model protein folding domains.

•

SAP domain folds through a diffuse transition state in which helix 1 is most formed.

•

Native state stability is dominated by contacts formed after the transition state.

The SAP domain from the Saccharomyces cerevisiae Tho1 protein is comprised of just two helices and a hydrophobic core and is one of the smallest proteins whose folding has been characterised. Φ-value analysis revealed that Tho1 SAP folds through a transition state where helix 1 is the most extensively formed element of secondary structure and flickering native-like core contacts from Leu35 are also present. The contacts that contribute most to native state stability of Tho1 SAP are not formed in the transition state.

Foldase and inhibitor functionalities of the pepsinogen prosegment are encoded within discrete segments of the 44 residue domain

Pepsin is initially produced as the zymogen pepsinogen, containing a 44 residue prosegment (PS) domain. When folded without the PS, pepsin forms a thermodynamically stable denatured state (refolded pepsin, Rp). To guide native folding, the PS binds to Rp, stabilizes the folding transition state, and binds tightly to native pepsin (Np), thereby driving the folding equilibrium to favor the native state. It is unknown whether these functionalities of the PS are encoded within the entire sequence or within discrete segments. PS residues 1p-29p correspond to a highly conserved region in pepsin-like aspartic proteases and we hypothesized that this segment is critical to PS-catalyzed folding. This notion was tested in the present study by characterizing the ability of various truncated PS peptides to bind Rp, catalyze folding from Rp to Np, and to inhibit Np. Four PS truncations were examined, corresponding to PS residues 1p-16p (PS1-16), 1p-29p (PS1-29), 17p-44p (PS17-44) and 30p-44p (PS30-44). The three PS functionalities could be ascribed primarily to discrete regions within the highly conserved motif: 1p-16p dictated Rp binding, 17p-29p dictated Np binding/inhibition, while the entire 1p-29p dictated transition state binding/catalyzing folding. Conversely, PS30-44 played no obvious role in PS-catalyzed folding; it is hypothesized that this more variable region may serve as a linker between PS1-29 and the mature domain. The high sequence conservation of PS1-29 and its role in catalyzing pepsin folding strongly suggest that there is a conserved PS-catalyzed folding mechanism shared by pepsin-like aspartic proteases with this motif.

Better theoretical models and protein design experiments can help to understand protein folding

In our study, we have concluded that two proteins with 88% homology choose different energetically favorable pathways in the very early stage of the folding process to attain their native folds. Subsequent reports from other investigators by performing folding and unfolding kinetics experiments concur with our findings. We herewith discuss the key papers revealing computational and experimental analysis of two designed proteins with similar sequence distant folds. Further we suggest that the theoretical/computational analysis of protein sequences and structures along with the relevant experiments provide a better understanding of the relationship between protein sequence, folding, and structure.

Take home lessons from studies of related proteins

The 'Fold Approach' involves a detailed analysis of the folding of several topologically, structurally and/or evolutionarily related proteins. Such studies can reveal determinants of the folding mechanism beyond the gross topology, and can dissect the residues required for folding from those required for stability or function. While this approach has not yet matured to the point where we can predict the native conformation of any polypeptide chain in silico, it has been able to highlight, amongst others, the specific residues that are responsible for nucleation, pathway malleability, kinetic intermediates, chain knotting, internal friction and Paracelsus switches. Some of the most interesting discoveries have resulted from the attempt to explain differences between homologues.