Folding Of Single-domain Proteins Research Articles

Molecular Information Theory Meets Protein Folding.

We propose an application of molecular information theory to analyze the folding of single domain proteins. We analyze results from various areas of protein science, such as sequence-based potentials, reduced amino acid alphabets, backbone configurational entropy, secondary structure content, residue burial layers, and mutational studies of protein stability changes. We found that the average information contained in the sequences of evolved proteins is very close to the average information needed to specify a fold ∼2.2 ± 0.3 bits/(site·operation). The effective alphabet size in evolved proteins equals the effective number of conformations of a residue in the compact unfolded state at around 5. We calculated an energy-to-information conversion efficiency upon folding of around 50%, lower than the theoretical limit of 70%, but much higher than human-built macroscopic machines. We propose a simple mapping between molecular information theory and energy landscape theory and explore the connections between sequence evolution, configurational entropy, and the energetics of protein folding.

Piecewise All-Atom SMD Simulations Reveal Key Secondary Structures in Luciferase Unfolding Pathway

Although the folding of single-domain proteins is well characterized theoretically and experimentally, the folding of large multidomain proteins is less well known. Firefly luciferase, a 550 residue three-domain protein, has been commonly used as a substrate to study chaperone reactions and as a model system for the study of folding of long polypeptide chains, including related phenomena such as cotranslational folding. Despite being characterized by various experimental techniques, the atomic-level contributions of various secondary structures of luciferase to its fold’s mechanical stability remain unknown. Here, we developed a piecewise approach for all-atom steered molecular dynamics simulations to examine specific secondary structures that resist mechanical unfolding while minimizing the amount of computational resources required by the large water box of standard all-atom steered molecular dynamics simulations. We validated the robustness of this approach with a small NI3C protein and used our approach to elucidate the specific secondary structures that provide the largest contributions to luciferase mechanostability. In doing so, we show that piecewise all-atom steered molecular dynamics simulations can provide novel atomic resolution details regarding mechanostability and can serve as a platform for novel mutagenesis studies as well as a point for comparison with high-resolution force spectroscopy experiments.

Confinement and Crowding Effects on Folding of a Multidomain Y-Family DNA Polymerase.

Proteins in vivo endure highly various interactions from the luxuriant surrounding macromolecular cosolutes. Confinement and macromolecular crowding are the two major effects that should be considered while comparing the results of protein dynamics from in vitro to in vivo. However, efforts have been largely focused on single domain protein folding up to now, and the quantifications of the in vivo effects in terms of confinements and crowders on modulating the structure and dynamics as well as the physical understanding of the underlying mechanisms on multidomain protein folding are still challenging. Here we developed a topology-based model to investigate folding of a multidomain Y-family DNA polymerase (DPO4) within spherical confined space and in the presence of repulsive and attractive crowders. We uncovered that the entropic component of the thermodynamic driving force led by confinements and repulsive crowders increases the stability of folded states relative to the folding intermediates and unfolded states, while the enthalpic component of the thermodynamic driving force led by attractive crowders gives rise to the opposite effects with less stability. We found that the shapes of DPO4 conformations influenced by the confinements and the crowders are quite different even when only the entropic component of the thermodynamic driving force is considered. We uncovered that under all in vivo conditions, the folding cooperativity of DPO4 decreases compared to that in bulk. We showed that the loss of folding cooperativity can promote the sequential domain-wise folding, which was widely found in cotranslational multidomain protein folding, and effectively prohibit the backtracking led by topological frustrations during multidomain protein folding processes.

Measuring the (Good) Solvent Quality of Water for Disordered Proteins from a Single SAXS Measurement

Significance: Proteins adopt non-native geometries during their synthesis, folding, transport, and turnover and a significant fraction of the proteome is intrinsically disordered. The characterization of disordered proteins is challenging in part due to difficulties in the interpretation of techniques including SAXS and FRET. We have developed a method to extract the dimensions and solvent quality from a single SAXS measurement. We apply it to intrinsically disordered proteins (IDPs), finding them to be expanded in water, which behaves as a good solvent.Abstract: SAXS has been used to measure the properties of denatured proteins and IDPs. Many ensembles are expanded in water, with little dependence on denaturant. These findings appear at odds with FRET studies, suggesting issues with data interpretation. We conduct SAXS measurements on IDPs up to 334 residues using in-line size-exclusion chromatography. The IDPs are expanded and the Flory exponent, describing the solvent quality and a chain's collapse propensity, is ν∼0.53. This value is between the ideal self-avoiding random walk (ν=0.59) and the theta solvent limit (ν=½) where intra-chain attraction exactly counterbalances excluded volume terms. Addition of denaturants results in a mild expansion, most noticeable in longer sequences, and an increase in ν to 0.59. Simulations generate a molecular form factor applicable to IDPs that can be used to determine the Rg and ν from a single measurement. The methodology is robust (hetero- or homopolymers, energy function, Ramachandran maps). This new analysis procedure and the dependence on chain length, along with recent FRET studies, largely resolves the controversy concerning the dimensions disordered proteins in water. These studies reaffirm that an early stable collapse often is not the first step in folding of single domain proteins upon dilution of denaturant.

A comparison of the folding kinetics of a small, artificially selected DNA aptamer with those of equivalently simple naturally occurring proteins

The folding of larger proteins generally differs from the folding of similarly large nucleic acids in the number and stability of the intermediates involved. To date, however, no similar comparison has been made between the folding of smaller proteins, which typically fold without well-populated intermediates, and the folding of small, simple nucleic acids. In response, in this study, we compare the folding of a 38-base DNA aptamer with the folding of a set of equivalently simple proteins. We find that, as is true for the large majority of simple, single domain proteins, the aptamer folds through a concerted, millisecond-scale process lacking well-populated intermediates. Perhaps surprisingly, the observed folding rate falls within error of a previously described relationship between the folding kinetics of single-domain proteins and their native state topology. Likewise, similarly to single-domain proteins, the aptamer exhibits a relatively low urea-derived Tanford β, suggesting that its folding transition state is modestly ordered. In contrast to this, however, and in contrast to the behavior of proteins, ϕ-value analysis suggests that the aptamer's folding transition state is highly ordered, a discrepancy that presumably reflects the markedly more important role that secondary structure formation plays in the folding of nucleic acids. This difference notwithstanding, the similarities that we observe between the two-state folding of single-domain proteins and the two-state folding of this similarly simple DNA presumably reflect properties that are universal in the folding of all sufficiently cooperative heteropolymers irrespective of their chemical details.

Folding Upon Binding - Not a Simple Protein Folding Problem

Although there have been many thermodynamic studies of the association of intrinsically disordered proteins (IDPs) with their binding partners, there have been relatively few kinetic studies.We have analysed the association of two different IDP systems, which form helical structure upon binding. These systems have markedly different behaviour. One associates very slowly and weakly, and the other very rapidly and strongly. The kinetic data therefore need to be analysed in a different manner to each other, and to the folding of single domain proteins. We discuss why standard approaches developed for folded protein-protein interactions may not be appropriate for reactions where one or both proteins are disordered. In addition, we show how protein engineering can be used to provide insights into the pathways/mechanisms of folding.Finally we have analysed the kinetics of the folding of an intrinsically disordered region (IDR) within a protein - a domain that only folds when attached to a neighbouring folded domain. This multidomain protein folds in a very different manner compared to those formed from adjacent domains which can fold in isolation.

Local and non-local native topologies reveal the underlying folding landscape of proteins

Due to Plaxco, Simons, Baker and others, it is now well known that the two-state single domain protein folding rate is fairly well predicted from knowledge of the topology of the native structure. Plaxco et al found that the folding rates of two-state proteins correlate with the average degree to which native contacts are ‘local’ within the chain sequence: fast-folders usually have mostly local structures. Here, we dissected the native topology further by focusing on non-local and local contacts using lower and upper bounds of allowable sequence separation in computing the average contact order. We analyzed non-local and local contacts of 82 two-state proteins whose experimental folding rates span over six orders of magnitude. We observed that both the number of non-local contacts and the average sequence separation of non-local contacts (non-local CO) are both negatively correlated with the folding rate, showing that the non-local contacts dominate the barrier-crossing process. Surprisingly, the local contact orders of the proteins also correlate with the folding rates. However, this correlation shows a strong positive trend indicating the role of a diffusive search in the denatured basin.

The Folding of Single Domain Proteins - have we Reached a Consensus?

I will discuss fundamental topics in protein folding where disagreement remains and where adequate experimental data is lacking. These topics include properties of the denatured state, the nature of the fundamental steps and barriers and the extent of pathway heterogeneity.View Large Image | View Hi-Res Image | Download PowerPoint Slide

Direct Observation of Tug-of-War during the Folding of a Mutually Exclusive Protein

Although most protein folding studies are carried out on single-domain proteins, over two-thirds of proteins in proteomes are built up from multiple individually folded domains. A significant fraction of these multidomain proteins are domain-insertion proteins, in which one guest domain is inserted into a surface loop of a host protein. Intricate thermodynamic and kinetic coupling between the two domains can have a profound impact on their folding dynamics. Here we use an engineered mutually exclusive protein as a model system to directly illustrate one such complex dynamic process: the "tug-of-war" process during protein folding. By inserting a guest protein I27w34f into a host protein GB1-L5 (GL5), we engineered a novel, mutually exclusive protein, GL5/I27w34f, in which only one domain can remain folded at any given time due to topological constraints imposed by the folded structures. Using stopped-flow techniques, we obtained the first kinetic evidence that the guest and host domains engage in a folding tug-of-war as they attempt to fold, in which the host domain folds rapidly into its three-dimensional structure and is then automatically unfolded, driven by the folding of the guest domain. Our results provided direct evidence that protein folding can generate sufficient mechanical strain to unravel a host protein. Using single-molecule atomic force microscopy, we provide direct evidence for the existence of a conformational equilibrium between the two mutually exclusive conformations. Our results highlight important roles played by the intricate coupling between folding kinetics, thermodynamic stability, and mechanical strain in the folding of complex multidomain proteins, which cannot be addressed in traditional single-domain protein folding studies.

Lattice simulations of cotranslational folding of single domain proteins

We use lattice protein models and Monte Carlo simulations to study cotranslational folding of small single domain proteins. We show that the assembly of native structure begins during late extrusion stages, but final formation of native state occurs during de novo folding, when all residues are extruded. There are three main results in our study. First, for the sequences displaying two-state refolding mechanism de novo cotranslational folding pathway differs from that sampled in in vitro refolding. The change in folding pathways is due to partial assembly of native interactions during extrusion that results in different starting conditions for in vitro refolding and for de novo cotranslational folding. For small single domain proteins cotranslational folding is slower than in vitro refolding, but is generally fast enough to be completed before the release from a ribosome. Second, we found that until final stages of biosynthesis cotranslational folding is essentially equilibrium. This observation is explained by low stability of structured states for partially extruded chains. Finally, our data suggest that the proteins, which refold in vitro slowly via intermediates, complete their de novo folding after the release from a ribosome. Comparison of our lattice cotranslational simulations with recent experimental and computational studies is discussed.

A folding inhibitor of the HIV‐1 protease

Because the human immunodeficiency virus type 1 protease (HIV-1-PR) is an essential enzyme in the viral life cycle, its inhibition can control AIDS. The folding of single-domain proteins, like each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES, folding units stabilized by strongly interacting, highly conserved, as a rule hydrophobic, amino acids). These LES have evolved over myriad generations to recognize and strongly attract each other, so as to make the protein fold fast and be stable in its native conformation. Consequently, peptides displaying a sequence identical to those segments of the monomers associated with LES are expected to act as competitive inhibitors and thus destabilize the native structure of the enzyme. These inhibitors are unlikely to lead to escape mutants as they bind to the protease monomers through highly conserved amino acids, which play an essential role in the folding process. The properties of one of the most promising inhibitors of the folding of the HIV-1-PR monomers found among these peptides are demonstrated with the help of spectrophotometric assays and circular dichroism spectroscopy.

Design of HIV‐1‐PR inhibitors that do not create resistance: Blocking the folding of single monomers

The main problems found in designing drugs are those of optimizing the drug-target interaction and of avoiding the insurgence of resistance. We suggest a scheme for the design of inhibitors that can be used as leads for the development of a drug and that do not face either of these problems, and then apply it to the case of HIV-1-PR. It is based on the knowledge that the folding of single-domain proteins, such as each of the monomers forming the HIV-1-PR homodimer, is controlled by local elementary structures (LES), stabilized by local contacts among hydrophobic, strongly interacting, and highly conserved amino acids that play a central role in the folding process. Because LES have evolved over many generations to recognize and strongly interact with each other so as to make the protein fold fast and avoid aggregation with other proteins, highly specific (and thus little toxic) as well as effective folding-inhibitor molecules suggest themselves: short peptides (or eventually their mimetic molecules) displaying the same amino acid sequence of that of LES (p-LES). Aside from being specific and efficient, these inhibitors are expected not to induce resistance; in fact, mutations in HIV-1-PR that successfully avoid the action of p-LES imply the destabilization of one or more LES and thus should lead to protein denaturation. Making use of Monte Carlo simulations, we first identify the LES of the HIV-1-PR and then show that the corresponding p-LES peptides act as effective inhibitors of the folding of the protease.

Apparent Debye−Huckel Electrostatic Effects in the Folding of a Simple, Single Domain Protein

We have monitored the effects of salts and denaturants on the folding of the simple, two-state protein FynSH3. As predicted by Debye-Huckel limiting law, both the stability and (log) folding rate of FynSH3 increase nearly perfectly linearly (r(2)> 0.99) with the square root of ionic strength upon increasing concentrations of the relatively nonchaotropic salt sodium chloride. The stability of FynSH3 is also linear in square root ionic strength when the relatively nonchaotropic salts sodium bromide, potassium bromide, and potassium chloride are employed. Comparison of the kinetic and equilibrium effects of sodium chloride suggests that the electrostatic interactions formed in the folding transition state are approximately 50% as destabilizing as those formed in the native state, presumably reflecting the more compact nature of the latter. In contrast, the relationship between concentration and folding kinetics is more complex when the highly chaotropic salt guanidine hydrochloride (GuHCl) is employed. At moderate to high GuHCl concentrations the net effect of the linear, presumably chaotrope-induced deceleration and the presumed, square root-dependent ionic strength-induced acceleration is well approximated as linear, thus accounting for the observation of "chevron behavior" (log folding rate linear in denaturant concentration) typically reported for the folding of single domain proteins. At very low GuHCl concentrations, however, significant kinetic rollover is observed. This rollover is reasonably well fitted as a sum of a linear, presumably chaotropic effect and a square root-dependent, presumably electrostatic effect. These results thus not only provide insight into the nature of the folding transition state but also suggest that caution is in order when extrapolating GuHCl-based chevrons to estimate folding rates in the absence of denaturant and in interpreting deviations from chevron linearity as evidence for non-two-state kinetics.

Chain length dependence of apomyoglobin folding: structural evolution from misfolded sheets to native helices.

Very little is known about how protein structure evolves during the polypeptide chain elongation that accompanies cotranslational protein folding. This in vitro model study is aimed at probing how conformational space evolves for purified N-terminal polypeptides of increasing length. These peptides are derived from the sequence of an all-alpha-helical single domain protein, Sperm whale apomyoglobin (apoMb). Even at short chain lengths, ordered structure is found. The nature of this structure is strongly chain length dependent. At relatively short lengths, a predominantly non-native beta-sheet conformation is present, and self-associated amyloid-like species are generated. As chain length increases, alpha-helix progressively takes over, and it replaces the beta-strand. The observed trends correlate with the specific fraction of solvent-accessible nonpolar surface area present at different chain lengths. The C-terminal portion of the chain plays an important role by promoting a large and cooperative overall increase in helical content and by consolidating the monomeric association state of the full-length protein. Thus, a native-like energy landscape develops late during apoMb chain elongation. This effect may provide an important driving force for chain expulsion from the ribosome and promote nearly-posttranslational folding of single domain proteins in the cell. Nature has been able to overcome the above intrinsic misfolding trends by modulating the composition of the intracellular environment. An imbalance or improper functioning by the above modulating factors during translation may play a role in misfolding-driven intracellular disorders.

Design and folding of dimeric proteins

In a similar way in which the folding of single-domain proteins provides an important test in the study of self-organization, the folding of homodimers constitutes a basic challenge in the quest for the mechanisms that are the basis of biological recognition. Dimerization is studied by following the evolution of two identical 20-letter amino acid chains within the framework of a lattice model and using Monte Carlo simulations. It is found that when design (evolution pressure) selects few, strongly interacting (conserved) amino acids to control the process, a three-state folding scenario follows, where the monomers first fold forming the halves of the eventual dimeric interface independently of each other, and then dimerize ("lock and key" kind of association). On the other hand, if design distributes the control of the folding process on a large number of (conserved) amino acids, a two-state folding scenario ensues, where dimerization takes place at the beginning of the process, resulting in an "induced type" of association. Making use of conservation patterns of families of analogous dimers, it is possible to compare the model predictions with the behavior of real proteins. It is found that theory provides an overall account of the experimental findings.

Consistency in structural energetics of protein folding and peptide recognition.

We report a new free energy decomposition that includes structure-derived atomic contact energies for the desolvation component, and show that it applies equally well to the analysis of single-domain protein folding and to the binding of flexible peptides to proteins. Specifically, we selected the 17 single-domain proteins for which the three-dimensional structures and thermodynamic unfolding free energies are available. By calculating all terms except the backbone conformational entropy change and comparing the result to the experimentally measured free energy, we estimated that the mean entropy gain by the backbone chain upon unfolding (delta Sbb) is 5.3 cal/K per mole of residue, and that the average backbone entropy for glycine is 6.7 cal/K. Both numbers are in close agreement with recent estimates made by entirely different methods, suggesting a promising degree of consistency between data obtained from disparate sources. In addition, a quantitative analysis of the folding free energy indicates that the unfavorable backbone entropy for each of the proteins is balanced predominantly by favorable backbone interactions. Finally, because the binding of flexible peptides to receptors is physically similar to folding, the free energy function should, in principle, be equally applicable to flexible docking. By combining atomic contact energies, electrostatics, and sequence-dependent backbone entropy, we calculated a priori the free energy changes associated with the binding of four different peptides to HLA-A2, 1 MHC molecule and found agreement with experiment to within 10% without parameter adjustment.