Simulations of a protein fold switch reveal crowding-induced population shifts driven by disordered regions.

Macromolecular crowding effects on globular proteins, which usually adopt a single stable fold, have been widely studied. However, little is known about crowding effects on fold-switching proteins, which reversibly switch between distinct folds. Here we study the mutationally driven switch between the folds of GA and GB, the two 56-amino acid binding domains of protein G, using a structure-based dual-basin model. We show that, in the absence of crowders, the fold populations PA and PB can be controlled by the strengths of contacts in the two folds, κA and κB. A population balance, PA ≈ PB, is obtained for κB/κA = 0.92. The resulting model protein is subject to crowding at different packing fractions, ϕc. We find that crowding increases the GB population and reduces the GA population, reaching PB/PA ≈ 4 at ϕc = 0.44. We analyze the ϕc-dependence of the crowding-induced GA-to-GB switch using scaled particle theory, which provides a qualitative, but not quantitative, fit of our data, suggesting effects beyond a spherical description of the folds. We show that the terminal regions of the protein chain, which are intrinsically disordered only in GA, play a dominant role in the response of the fold switch to crowding effects.

Crowding-induced protein destabilization in the absence of soft attractions.

It is generally assumed that volume exclusion by macromolecular crowders universally stabilizes the native states of proteins and destabilization suggests soft attractions between crowders and protein. Here we show that proteins can be destabilized even by crowders that are purely repulsive. With a coarse-grained sequence-based model, we study the folding thermodynamics of two sequences with different native folds, a helical hairpin and a ß-barrel, in a range of crowder volume fractions, φc. We find that the native state, N, remains structurally unchanged under crowded conditions, while the size of the unfolded state, U, decreases monotonically with φc. Hence, for all φc>0, U is entropically disfavored relative to N. This entropy-centric view holds for the helical hairpin protein, which is stabilized under all crowded conditions as quantified by changes in either the folding midpoint temperature, Tm, or the free energy of folding. We find, however, that the ß-barrel protein is destabilized under low-T, low-φc conditions. This destabilization can be understood from two characteristics of its folding: 1) a relatively compact U at T

The C-terminal domain of transcription factor RfaH: Folding, fold switching and energy landscape.

We simulate the folding and fold switching of the C-terminal domain (CTD) of the transcription factor RfaH using an all-atom physics-based model augmented with a dual-basin structure-based potential energy term. We show that this hybrid model captures the essential thermodynamic behavior of this metamorphic domain, that is, a change in the global free energy minimum from an α-helical hairpin to a 5-stranded ß-barrel upon the dissociation of the CTD from the rest of the protein. Using Monte Carlo sampling techniques, we then analyze the energy landscape of the CTD in terms of progress variables for folding toward the two folds. We find that, below the folding transition, the energy landscape is characterized by a single, dominant funnel to the native ß-barrel structure. The absence of a deep funnel to the α-helical hairpin state reflects a negligible population of this fold for the isolated CTD. We observe, however, a higher α-helix structure content in the unfolded state compared to results from a similar but fold switch-incompetent version of our model. Moreover, in folding simulations started from an extended chain conformation we find transiently formed α-helical structure, occurring early in the process and disappearing as the chain progresses toward the thermally stable ß-barrel state.

Structural fluctuations and mechanical stabilities of the metamorphic protein RfaH.

RfaH is a compact two-domain bacterial transcription factor that functions both as a regulator of transcription and an enhancer of translation. Underpinning the dual functional roles of RfaH is a partial but dramatic fold switch, which completely transforms the ~50-amino acid C-terminal domain (CTD) from an all-α state to an all-ß state. The fold switch of the CTD occurs when RfaH binds to RNA polymerase (RNAP), however, the details of how this structural transformation is triggered is not well understood. Here we use all-atom Monte Carlo simulations to characterize structural fluctuations and mechanical stability properties of the full-length RfaH and the CTD as an isolated fragment. In agreement with experiments, we find that interdomain contacts are crucial for maintaining a stable, all-α CTD in free RfaH. To probe mechanical properties, we use pulling simulations to measure the work required to inflict local deformations at different positions along the chain. The resulting mechanical stability profile reveals that free RfaH can be divided into a "rigid" part and a "soft" part, with a boundary that nearly coincides with the boundary between the two domains. We discuss the potential role of this feature for how fold switching may be triggered by interaction with RNAP.

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.

Phase behavior of blocky charge lattice polymers: Crystals, liquids, sheets, filaments, and clusters.

Motivated by the idea that intrinsically disordered proteins (IDPs) condense into liquidlike droplets within cells, we carry out Monte Carlo simulations of a polymer lattice model to study the relationship between charge patterning and phase separation. Polymer chains containing neutral, positively charged, and negatively charged monomers are placed on a cubic lattice. Only nearest-neighbor interactions between charges are considered. We determine the phase diagram for a systematically varied set of sequences. We observe homogeneous fluids, liquid condensation, cluster phases, filaments, and crystal states. Of the six sequences we study, three form crystals at low temperatures. The other three sequences, which have lower charge densities, instead collapse into gel-like networks or unconnected finite clusters. Longer neutral patches along the sequence sterically limit the size and shape of low-energy structures, which is analogous to the effect of charge or limited valence in attractive colloids. Only one sequence clearly exhibits liquid behavior; this sequence has a reduced tendency to individually fold and crystallize compared to others of similar charge density and draws parallels to real IDP behavior.

Characterization of a peptide containing the major heparin binding domain of human hepatic lipase.

Human hepatic lipase (hHL) is a cell surface associated enzyme that hydrolyzes triacylglycerols and phospholipids within circulating lipoproteins. We hypothesized that an amino acid sequence mimicking the major heparin binding domain (HBD) of hHL will displace hHL from cell surfaces. To test this hypothesis, we generated a recombinant protein of thioredoxin linked with a cleavable, tagged sequence containing amino acids 442 to 476 of the mature hHL sequence, which contains the major HBD of hHL. The recombinant protein associated with heparin-sepharose, and its peak elution from heparin-sepharose occurred in the presence of 0.5 M NaCl. We cleaved and purified the tagged sequence containing the HBD from the recombinant protein and tested the ability of the peptide to displace full-length hHL from HEK-293 cells. The peptide indeed displaced hHL from cell surfaces, while no significant displacement was observed in the presence of a peptide with a scrambled sequence. Finally, we obtained structural information for the peptide containing the HBD. 1 H- and 15 N-NMR spectra of the peptide indicate the peptide is largely unstructured, although not completely random coil. The addition of heparin to the peptide induced some changes in chemical shift, suggesting changes in peptide structure and/or specific interactions with heparin. Molecular simulations confirm the largely unstructured nature of the isolated peptide, but they also indicate weak tendencies for both α- and ß-structure formation in different parts of the chain. Overall, these data provide a proof-of-principle for the use of mimetic peptides for the displacement of cell surface associated lipases.

Structures of apolipoprotein A-I in high density lipoprotein generated by electron microscopy and biased simulations.

BACKGROUND: Apolipoprotein A-I (apoA-I) in high-density lipoprotein (HDL) is a key protein for the transport of cholesterol from the vascular wall to the liver. The formation and structure of nascent HDL, composed of apoA-I and phospholipids, is critical to this process. METHODS: The HDL was assembled in vitro from apoA-I, cholesterol and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) at a 1:4:50 molar ratio. The structure of HDL was investigated in vitreous samples, frozen at cryogenic temperatures, as well as in negatively stained samples by transmission electron microscopy. Low resolution electron density maps were next used as restraints in biased Monte Carlo simulations of apolipoprotein A-I dimers, with an initial structure derived from atomic resolution X-ray structures. RESULTS: Two final apoA-I structure models for the full-length structure of apoA-I dimer in the lipid bound conformation were generated, showing a nearly circular, flat particle with an uneven particle thickness. CONCLUSIONS: The generated structures provide evidence for the discoidal, antiparallel arrangement of apoA-I in nascent HDL, and propose two preferred conformations of the flexible N-termini. GENERAL SIGNIFICANCE: The novel full-length structures of apoA-I dimers deepens the understanding to the structure-function relationship of nascent HDL with significance for the prevention of lipoprotein-related disease. The biased simulation method used in this study provides a powerful and convenient modelling tool with applicability for structural studies and modelling of other proteins and protein complexes.

Binding Specificity Profiles from Computational Peptide Screening.

The computational peptide screening method is a Monte Carlo-based procedure to systematically characterize the specificity of a peptide-binding site. The method is based on a generalized-ensemble algorithm in which the peptide sequence has become a dynamic variable, i.e., molecular simulations with ordinary conformational moves are enhanced with a type of "mutational" move such that proper statistics are achieved for multiple sequences in a single run. The peptide screening method has two main steps. In the first, reference simulations of the unbound state are performed and used to parametrize a linear model of the unbound state free energy, determined by requiring that the marginal distribution of peptide sequences is approximately flat. In the second step, simulations of the bound state are performed. By using the linear model as a free energy reference point, the marginal distribution of peptide sequences becomes skewed towards sequences with higher binding free energies. From analyses of the sequences generated in the second step and their conformational ensembles, information on peptide binding specificity, relative binding affinities, and the molecular basis of specificity can be achieved. Here we demonstrate how the algorithm can be implemented and applied to determine the peptide binding specificity of a PDZ domain from the protein GRIP1.

Local versus global fold switching in protein evolution: insight from a three-letter continuous model.

Recent design experiments have demonstrated that some proteins can switch their folds in response to a small number of point mutations either directly, in a single mutational step, or via intermediate bistable sequences, which populate two different folds simultaneously. Here we explore the hypothesis that bistable intermediates are more common in switches between structurally similar folds while direct switches are more common between dissimilar folds. To this end, we use a reduced model with seven atoms per amino acid and three amino acid types as a biophysical basis for protein folding and stability. We compare a set of mutational pathways, selected for optimal stability properties, that lead to switches between ß-hairpin and α-helix folds with 16 amino acids and between [Formula: see text] and [Formula: see text] folds with 35 amino acids, respectively. Fold switching in each case is sharp, taking only a few mutations to be completed. While the sharpness of mutationally driven protein fold switching can be traced to a shift in the energy balance of the two native states, conformational entropy contributes to determining the point at which fold switching occurs along a pathway.

Smooth functional transition along a mutational pathway with an abrupt protein fold switch.

Recent protein design experiments have demonstrated that proteins can migrate between folds through the accumulation of substitution mutations without visiting disordered or nonfunctional points in sequence space. To explore the biophysical mechanism underlying such transitions we use a three-letter continuous protein model with seven atoms per amino acid to provide realistic sequence-structure and sequence-function mappings through explicit simulation of the folding and interaction of model sequences. We start from two 16-amino-acid sequences folding into an α-helix and a ß-hairpin, respectively, each of which has a preferred binding partner with 35 amino acids. We identify a mutational pathway between the two folds, which features a sharp fold switch. By contrast, we find that the transition in function is smooth. Moreover, the switch in preferred binding partner does not coincide with the fold switch. Discovery of new folds in evolution might therefore be facilitated by following fitness slopes in sequence space underpinned by binding-induced conformational switching.

Conformational and aggregation properties of the 1-93 fragment of apolipoprotein A-I.

Several disease-linked mutations of apolipoprotein A-I, the major protein in high-density lipoprotein (HDL), are known to be amyloidogenic, and the fibrils often contain N-terminal fragments of the protein. Here, we present a combined computational and experimental study of the fibril-associated disordered 1-93 fragment of this protein, in wild-type and mutated (G26R, S36A, K40L, W50R) forms. In atomic-level Monte Carlo simulations of the free monomer, validated by circular dichroism spectroscopy, we observe changes in the position-dependent ß-strand probability induced by mutations. We find that these conformational shifts match well with the effects of these mutations in thioflavin T fluorescence and transmission electron microscopy experiments. Together, our results point to molecular mechanisms that may have a key role in disease-linked aggregation of apolipoprotein A-I.

Exploring Protein-Peptide Binding Specificity through Computational Peptide Screening.

The binding of short disordered peptide stretches to globular protein domains is important for a wide range of cellular processes, including signal transduction, protein transport, and immune response. The often promiscuous nature of these interactions and the conformational flexibility of the peptide chain, sometimes even when bound, make the binding specificity of this type of protein interaction a challenge to understand. Here we develop and test a Monte Carlo-based procedure for calculating protein-peptide binding thermodynamics for many sequences in a single run. The method explores both peptide sequence and conformational space simultaneously by simulating a joint probability distribution which, in particular, makes searching through peptide sequence space computationally efficient. To test our method, we apply it to 3 different peptide-binding protein domains and test its ability to capture the experimentally determined specificity profiles. Insight into the molecular underpinnings of the observed specificities is obtained by analyzing the peptide conformational ensembles of a large number of binding-competent sequences. We also explore the possibility of using our method to discover new peptide-binding pockets on protein structures.

Aggregate geometry in amyloid fibril nucleation.

We present and study a minimal structure-based model for the self-assembly of peptides into ordered ß-sheet-rich fibrils. The peptides are represented by unit-length sticks on a cubic lattice and interact by hydrogen bonding and hydrophobicity forces. Using Monte Carlo simulations with >10(5) peptides, we show that fibril formation occurs with sigmoidal kinetics in the model. To determine the mechanism of fibril nucleation, we compute the joint distribution in length and width of the aggregates at equilibrium, using an efficient cluster move and flat-histogram techniques. This analysis, based on simulations with 256 peptides in which aggregates form and dissolve reversibly, shows that the main free-energy barriers that a nascent fibril has to overcome are associated with changes in width.

In silico peptide-binding predictions of passerine MHC class I reveal similarities across distantly related species, suggesting convergence on the level of protein function.

The major histocompatibility complex (MHC) genes are the most polymorphic genes found in the vertebrate genome, and they encode proteins that play an essential role in the adaptive immune response. Many songbirds (passerines) have been shown to have a large number of transcribed MHC class I genes compared to most mammals. To elucidate the reason for this large number of genes, we compared 14 MHC class I alleles (α1-α3 domains), from great reed warbler, house sparrow and tree sparrow, via phylogenetic analysis, homology modelling and in silico peptide-binding predictions to investigate their functional and genetic relationships. We found more pronounced clustering of the MHC class I allomorphs (allele specific proteins) in regards to their function (peptide-binding specificities) compared to their genetic relationships (amino acid sequences), indicating that the high number of alleles is of functional significance. The MHC class I allomorphs from house sparrow and tree sparrow, species that diverged 10 million years ago (MYA), had overlapping peptide-binding specificities, and these similarities across species were also confirmed in phylogenetic analyses based on amino acid sequences. Notably, there were also overlapping peptide-binding specificities in the allomorphs from house sparrow and great reed warbler, although these species diverged 30 MYA. This overlap was not found in a tree based on amino acid sequences. Our interpretation is that convergent evolution on the level of the protein function, possibly driven by selection from shared pathogens, has resulted in allomorphs with similar peptide-binding repertoires, although trans-species evolution in combination with gene conversion cannot be ruled out.

Epitope-specificity of recombinant antibodies reveals promiscuous peptide-binding properties.

Protein-peptide interactions are a common occurrence and essential for numerous cellular processes, and frequently explored in broad applications within biology, medicine, and proteomics. Therefore, understanding the molecular mechanism(s) of protein-peptide recognition, specificity, and binding interactions will be essential. In this study, we report the first detailed analysis of antibody-peptide interaction characteristics, by combining large-scale experimental peptide binding data with the structural analysis of eight human recombinant antibodies and numerous peptides, targeting tryptic mammalian and eukaryote proteomes. The results consistently revealed that promiscuous peptide-binding interactions, that is, both specific and degenerate binding, were exhibited by all antibodies, and the discovery was corroborated by orthogonal data, indicating that this might be a general phenomenon for low-affinity antibody-peptide interactions. The molecular mechanism for the degenerate peptide-binding specificity appeared to be executed through the use of 2-3 semi-conserved anchor residues in the C-terminal part of the peptides, in analogue to the mechanism utilized by the major histocompatibility complex-peptide complexes. In the long-term, this knowledge will be instrumental for advancing our fundamental understanding of protein-peptide interactions, as well as for designing, generating, and applying peptide specific antibodies, or peptide-binding proteins in general, in various biotechnical and medical applications.

Binding of two intrinsically disordered peptides to a multi-specific protein: a combined Monte Carlo and molecular dynamics study.

The unique ability of intrinsically disordered proteins (IDPs) to fold upon binding to partner molecules makes them functionally well-suited for cellular communication networks. For example, the folding-binding of different IDP sequences onto the same surface of an ordered protein provides a mechanism for signaling in a many-to-one manner. Here, we study the molecular details of this signaling mechanism by applying both Molecular Dynamics and Monte Carlo methods to S100B, a calcium-modulated homodimeric protein, and two of its IDP targets, p53 and TRTK-12. Despite adopting somewhat different conformations in complex with S100B and showing no apparent sequence similarity, the two IDP targets associate in virtually the same manner. As free chains, both target sequences remain flexible and sample their respective bound, natively [Formula: see text]-helical states to a small extent. Association occurs through an intermediate state in the periphery of the S100B binding pocket, stabilized by nonnative interactions which are either hydrophobic or electrostatic in nature. Our results highlight the importance of overall physical properties of IDP segments, such as net charge or presence of strongly hydrophobic amino acids, for molecular recognition via coupled folding-binding.

Coupled folding-binding in a hydrophobic/polar protein model: impact of synergistic folding and disordered flanks.

Coupled folding-binding is central to the function of many intrinsically disordered proteins, yet not fully understood. With a continuous three-letter protein model, we explore the free-energy landscape of pairs of interacting sequences and how it is impacted by 1), variations in the binding mechanism; and 2), the addition of disordered flanks to the binding region. In particular, we focus on two sequences, one with 16 and one with 35 amino acids, which make a stable dimeric three-helix bundle at low temperatures. Three distinct binding mechanisms are realized by altering the stabilities of the individual monomers: docking, coupled folding-binding of a single α-helix, and synergistic folding and binding. Compared to docking, the free-energy barrier for binding is reduced when the single α-helix is allowed to fold upon binding, but only marginally. A greater reduction is found for synergistic folding, which in addition results in a binding transition state characterized by very few interchain contacts. Disordered flanking chain segments attached to the α-helix sequence can, despite a negligible impact on the dimer stability, lead to a downhill free-energy surface in which the barrier for binding is eliminated.

Binding free energy landscape of domain-peptide interactions.

Peptide recognition domains (PRDs) are ubiquitous protein domains which mediate large numbers of protein interactions in the cell. How these PRDs are able to recognize peptide sequences in a rapid and specific manner is incompletely understood. We explore the peptide binding process of PDZ domains, a large PRD family, from an equilibrium perspective using an all-atom Monte Carlo (MC) approach. Our focus is two different PDZ domains representing two major PDZ classes, I and II. For both domains, a binding free energy surface with a strong bias toward the native bound state is found. Moreover, both domains exhibit a binding process in which the peptides are mostly either bound at the PDZ binding pocket or else interact little with the domain surface. Consistent with this, various binding observables show a temperature dependence well described by a simple two-state model. We also find important differences in the details between the two domains. While both domains exhibit well-defined binding free energy barriers, the class I barrier is significantly weaker than the one for class II. To probe this issue further, we apply our method to a PDZ domain with dual specificity for class I and II peptides, and find an analogous difference in their binding free energy barriers. Lastly, we perform a large number of fixed-temperature MC kinetics trajectories under binding conditions. These trajectories reveal significantly slower binding dynamics for the class II domain relative to class I. Our combined results are consistent with a binding mechanism in which the peptide C terminal residue binds in an initial, rate-limiting step.

Cooperativity, local-nonlocal coupling, and nonnative interactions: principles of protein folding from coarse-grained models.

Coarse-grained, self-contained polymer models are powerful tools in the study of protein folding. They are also essential to assess predictions from less rigorous theoretical approaches that lack an explicit-chain representation. Here we review advances in coarse-grained modeling of cooperative protein folding, noting in particular that the Levinthal paradox was raised in response to the experimental discovery of two-state-like folding in the late 1960s, rather than to the problem of conformational search per se. Comparisons between theory and experiment indicate a prominent role of desolvation barriers in cooperative folding, which likely emerges generally from a coupling between local conformational preferences and nonlocal packing interactions. Many of these principles have been elucidated by native-centric models, wherein nonnative interactions may be treated perturbatively. We discuss these developments as well as recent applications of coarse-grained chain modeling to knotted proteins and to intrinsically disordered proteins.