Frustration and hydrophobicity interplay in protein folding and protein evolution

A lattice model is used to study mutations and compacting effects on protein folding rates and folding temperature. In the context of protein evolution, we address the question regarding the best scenario for a polypeptide chain to fold: either a fast nonspecific collapse followed by a slow rearrangement to form the native structure or a specific collapse from the unfolded state with the simultaneous formation of the native state. This question is investigated for optimized sequences, whose native state has no frustrated contacts between monomers, and also for mutated sequences, whose native state has some degree of frustration. It is found that the best scenario for folding may depend on the amount of frustration of the native structure. The implication of this result on protein evolution is discussed. (c) 2006 American Institute of Physics.

Analyzing the effect of homogeneous frustration in protein folding

The energy landscape theory has been an invaluable theoretical framework in the understanding of biological processes such as protein folding, oligomerization, and functional transitions. According to the theory, the energy landscape of protein folding is funneled toward the native state, a conformational state that is consistent with the principle of minimal frustration. It has been accepted that real proteins are selected through natural evolution, satisfying the minimum frustration criterion. However, there is evidence that a low degree of frustration accelerates folding. We examined the interplay between topological and energetic protein frustration. We employed a Cα structure-based model for simulations with a controlled nonspecific energetic frustration added to the potential energy function. Thermodynamics and kinetics of a group of 19 proteins are completely characterized as a function of increasing level of energetic frustration. We observed two well-separated groups of proteins: one group where a little frustration enhances folding rates to an optimal value and another where any energetic frustration slows down folding. Protein energetic frustration regimes and their mechanisms are explained by the role of non-native contact interactions in different folding scenarios. These findings strongly correlate with the protein free-energy folding barrier and the absolute contact order parameters. These computational results are corroborated by principal component analysis and partial least square techniques. One simple theoretical model is proposed as a useful tool for experimentalists to predict the limits of improvements in real proteins. © 2013 Wiley Periodicals, Inc.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Minimal and optimal mechanisms for GroE-mediated protein folding

We have analyzed the effects of different components of the GroE chaperonin system on protein folding by using a nonpermissive substrate (i.e., one that has very low spontaneous refolding yield) for which rate data can be acquired. In the absence of GroES and nucleotides, the rate of GroEL-mediated refolding of heat- and DTT-denatured mitochondrial malate dehydrogenase was extremely low, but some three times higher than the spontaneous rate. This GroEL-mediated rate was increased 17-fold by saturating concentrations of ATP, 11-fold by ADP and GroES, and 465-fold by ATP and GroES. Optimal refolding activity was observed when the dissociation of GroES from the chaperonin complex was dramatically reduced. Although GroEL minichaperones were able to bind denatured mitochondrial malate dehydrogenase, they were ineffective in enhancing the refolding rate. The spectrum of mechanisms for GroE-mediated protein folding depends on the nature of the substrate. The minimal mechanism for permissive substrates (i.e., having significant yields of spontaneous refolding), requires only binding to the apical domain of GroEL. Slow folding rates of nonpermissive substrates are limited by the transitions between high- and low-affinity states of GroEL alone. The optimal mechanism, which requires holoGroEL, physiological amounts of GroES, and ATP hydrolysis, is necessary for the chaperonin-mediated folding of nonpermissive substrates at physiologically relevant rates under conditions in which retention of bound GroES prevents the premature release of aggregation-prone folding intermediates from the chaperonin complex. The different mechanisms are described in terms of the structural features of mini- and holo-chaperones.

Configurational diffusion down a folding funnel describes the dynamics of DNA hairpins

Elucidating the mechanism of folding of polynucleotides depends on accurate estimates of free energy surfaces and a quantitative description of the kinetics of structure formation. Here, the kinetics of hairpin formation in single-stranded DNA are measured after a laser temperature jump. The kinetics are modeled as configurational diffusion on a free energy surface obtained from a statistical mechanical description of equilibrium melting profiles. The effective diffusion coefficient is found to be strongly temperature-dependent in the nucleation step as a result of formation of misfolded loops that do not lead to subsequent zipping. This simple system exhibits many of the features predicted from theoretical studies of protein folding, including a funnel-like energy surface with many folding pathways, trapping in misfolded conformations, and non-Arrhenius folding rates.

Synthesis of aromatic monothiols and aromatic dithiols to increase the folding rate and yield of disulfide containing proteins

Most pharmaceutically relevant proteins and many extracellular proteins contain disulfide bonds. Formation of the correct disulfide bonds is essential for stability in almost all cases. Disulfide containing proteins can be rapidly and inexpensively overexpressed in bacteria. However, the overexpressed proteins usually form aggregates inside the bacteria, called inclusion bodies, which contains inactive and non-native protein. To obtain native protein, inclusion bodies need to be isolated and resolubilized, and then the resulting protein refolded in vitro. In vitro protein folding is aided by the addition of a redox buffer, which is composed of a small molecule disulfide and/or a small molecule thiol. The most commonly used redox buffer contains reduced and oxidized glutathione. Recently, aliphatic dithiols and aromatic monothiols have been employed as redox buffers. Aliphatic dithiols improved the yield of native protein as compared to the aliphatic thiol, glutathione. Dithiols mimic the in vivo protein folding catalyst, protein disulfide isomerase, which has two thiols per active site. Furthermore, aromatic monothiols increased the folding rate and yield of lysozyme and RNase A relative to glutathione. By combining the beneficial properties of aliphatic dithiols and aromatic monothiols, aromatic dithiols were designed and were expected to increase in vitro protein folding rates and yields. Aromatic monothiols (1-4) and their corresponding disulfides (5-8), two series of ortho- and para-substituted ethylene glycol dithiols (9-15), and a series of aromatic quaternary ammonium salt dithiols (16-17) were synthesized on a multigram scale. Monothiols and disulfides (1-8) were utilized to fold lysozyme and bovine pancreatic trypsin inhibitor. Dithiols (11-17) were tested for their ability to fold lysozyme. At pH 7.0 and pH 8.0, and high protein concentration (1 mg/mL), aromatic dithiols (16, 17) and a monothiol (3) significantly enhanced the in vitro folding rate and yield of lysozyme relative to the aliphatic thiol, glutathione. Additionally, aromatic dithiols (16, 17) significantly enhance the folding yield as compared to the corresponding aromatic monothiol (3). Thus, the folding rate and yield enhancements achieved in in vitro protein folding at high protein concentration will decrease the volume of renaturation solution required for large scale processes and consequently reduce processing time and cost.

Synthesis of aromatic monothiols to increase the folding rate and yields of disulfide-containing proteins

Almost all pharmaceutically relevant proteins and many extracellular proteins contain disulfide bonds, which are essential for protein functions. In many cases, disulfidecontaining proteins are produced via in vitro protein folding that involves the oxidation of reduced protein to native protein, a complex process. The in vitro folding of reduced lysozyme has been extensively studied as a model system because native lysozyme is small, inexpensive, and has only four disulfide bonds. The folding of reduced lysozyme is conducted with the aid of a redox buffer consisting of a small molecule disulfide and a small molecule thiol, such as oxidized and reduced glutathione. Herein, in vitro folding rates and yields of lysozyme obtained in the presence of a series of aromatic thiols and oxidized glutathione are compared to those obtained with reduced and oxidized glutathione. Results showed that aromatic thiols significantly increase the folding rate of lysozyme compared to glutathione.

STRUCTURE AND FUNCTION OF THE GROUP III CHAPERONINS, A UNIQUE CLADE OF PROTEIN FOLDING NANOMACHINES

The survival and descent of cells is universally dependent on maintaining their proteins in a properly folded condition. It is widely accepted that the information for the folding of the nascent polypeptide chain into a native protein is encrypted in the amino acid sequence, and the Nobel Laureate Christian Anfinsen was the first to demonstrate that a protein could spontaneously refold after complete unfolding. However, it became clear that the observed folding rates for many proteins were much slower than rates estimated in vivo. This led to the recognition of required protein-protein interactions that promote proper folding. A unique group of proteins, the molecular chaperones, are responsible for maintaining protein homeostasis during normal growth as well as stress conditions. Chaperonins (CPNs) are ubiquitous and essential chaperones. They form ATP-dependent, hollow complexes that encapsulate polypeptides in two back-to-back stacked multisubunit rings, facilitating protein folding through highly cooperative allosteric articulation. CPNs are usually classified into Group I and Group II. Here, I report the characterization of a novel CPN belonging to a third Group, recently discovered in bacteria. Group III CPNs have close phylogenetic association to the Group II CPNs found in Archaea and Eukarya, and may be a relic of the Last Common Ancestor of the CPN family. The gene encoding the Group III CPN from Carboxydothermus hydrogenoformans and Candidatus Desulforudis audaxviator was cloned in E. coli and overexpressed in order to both characterize the protein and to demonstrate its ability to function as an ATPase chaperone. The opening and closing cycle of the Chy chaperonin was examined via site-directed mutations affecting the ATP binding site at R155. To relate the mutational analysis to the structure of the CPN, the crystal structure of both the AMP-PNP (an ATP analogue) and ADP bound forms were obtained in collaboration with Sun-Shin Cha in Seoul, South Korea. The ADP and ATP binding site substitutions resulted in frozen forms of the structures in open and closed conformations. From this, mutants were designed to validate hypotheses regarding key ATP interacting sites as well as important stabilizing interactions, and to observe the physical properties of the resulting complexes by calorimetry.

Stabilization of native protein fold by intein-mediated covalent cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10degreesC by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG = 2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding I measured by stopped-flow circular dichroism. rate of less than 4 min(-1) The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E. coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding. (C) 2005 Elsevier Ltd. All rights reserved.

Predicting residue-wise contact orders in proteins by support vector regression

Background: The residue-wise contact order (RWCO) describes the sequence separations between the residues of interest and its contacting residues in a protein sequence. It is a new kind of one-dimensional protein structure that represents the extent of long-range contacts and is considered as a generalization of contact order. Together with secondary structure, accessible surface area, the B factor, and contact number, RWCO provides comprehensive and indispensable important information to reconstructing the protein three-dimensional structure from a set of one-dimensional structural properties. Accurately predicting RWCO values could have many important applications in protein three-dimensional structure prediction and protein folding rate prediction, and give deep insights into protein sequence-structure relationships. Results: We developed a novel approach to predict residue-wise contact order values in proteins based on support vector regression (SVR), starting from primary amino acid sequences. We explored seven different sequence encoding schemes to examine their effects on the prediction performance, including local sequence in the form of PSI-BLAST profiles, local sequence plus amino acid composition, local sequence plus molecular weight, local sequence plus secondary structure predicted by PSIPRED, local sequence plus molecular weight and amino acid composition, local sequence plus molecular weight and predicted secondary structure, and local sequence plus molecular weight, amino acid composition and predicted secondary structure. When using local sequences with multiple sequence alignments in the form of PSI-BLAST profiles, we could predict the RWCO distribution with a Pearson correlation coefficient (CC) between the predicted and observed RWCO values of 0.55, and root mean square error (RMSE) of 0.82, based on a well-defined dataset with 680 protein sequences. Moreover, by incorporating global features such as molecular weight and amino acid composition we could further improve the prediction performance with the CC to 0.57 and an RMSE of 0.79. In addition, combining the predicted secondary structure by PSIPRED was found to significantly improve the prediction performance and could yield the best prediction accuracy with a CC of 0.60 and RMSE of 0.78, which provided at least comparable performance compared with the other existing methods. Conclusion: The SVR method shows a prediction performance competitive with or at least comparable to the previously developed linear regression-based methods for predicting RWCO values. In contrast to support vector classification (SVC), SVR is very good at estimating the raw value profiles of the samples. The successful application of the SVR approach in this study reinforces the fact that support vector regression is a powerful tool in extracting the protein sequence-structure relationship and in estimating the protein structural profiles from amino acid sequences.

Using Helix-coil Models to Study Protein Unfolded States

An abstract of a thesis devoted to using helix-coil models to study unfolded states.\\

Research on polypeptide unfolded states has received much more attention in the last decade or so than it has in the past. Unfolded states are thought to be implicated in various

misfolding diseases and likely play crucial roles in protein folding equilibria and folding rates. Structural characterization of unfolded states has proven to be

much more difficult than the now well established practice of determining the structures of folded proteins. This is largely because many core assumptions underlying

folded structure determination methods are invalid for unfolded states. This has led to a dearth of knowledge concerning the nature of unfolded state conformational

distributions. While many aspects of unfolded state structure are not well known, there does exist a significant body of work stretching back half a century that

has been focused on structural characterization of marginally stable polypeptide systems. This body of work represents an extensive collection of experimental

data and biophysical models associated with describing helix-coil equilibria in polypeptide systems. Much of the work on unfolded states in the last decade has not been devoted

specifically to the improvement of our understanding of helix-coil equilibria, which arguably is the most well characterized of the various conformational equilibria

that likely contribute to unfolded state conformational distributions. This thesis seeks to provide a deeper investigation of helix-coil equilibria using modern

statistical data analysis and biophysical modeling techniques. The studies contained within seek to provide deeper insights and new perspectives on what we presumably

know very well about protein unfolded states. \\

Chapter 1 gives an overview of recent and historical work on studying protein unfolded states. The study of helix-coil equilibria is placed in the context

of the general field of unfolded state research and the basics of helix-coil models are introduced.\\

Chapter 2 introduces the newest incarnation of a sophisticated helix-coil model. State of the art modern statistical techniques are employed to estimate the energies

of various physical interactions that serve to influence helix-coil equilibria. A new Bayesian model selection approach is utilized to test many long-standing

hypotheses concerning the physical nature of the helix-coil transition. Some assumptions made in previous models are shown to be invalid and the new model

exhibits greatly improved predictive performance relative to its predecessor. \\

Chapter 3 introduces a new statistical model that can be used to interpret amide exchange measurements. As amide exchange can serve as a probe for residue-specific

properties of helix-coil ensembles, the new model provides a novel and robust method to use these types of measurements to characterize helix-coil ensembles experimentally

and test the position-specific predictions of helix-coil models. The statistical model is shown to perform exceedingly better than the most commonly used

method for interpreting amide exchange data. The estimates of the model obtained from amide exchange measurements on an example helical peptide

also show a remarkable consistency with the predictions of the helix-coil model. \\

Chapter 4 involves a study of helix-coil ensembles through the enumeration of helix-coil configurations. Aside from providing new insights into helix-coil ensembles,

this chapter also introduces a new method by which helix-coil models can be extended to calculate new types of observables. Future work on this approach could potentially

allow helix-coil models to move into use domains that were previously inaccessible and reserved for other types of unfolded state models that were introduced in chapter 1.

Optimization of rates of protein folding: the nucleation-condensation mechanism and its implications.

Small, single-module proteins that fold in a single cooperative step may be paradigms for understanding early events in protein-folding pathways generally. Recent experimental studies of the 64-residue chymotrypsin inhibitor 2 (CI2) support a nucleation mechanism for folding, as do some computer stimulations. CI2 has a nucleation site that develops only in the transition state for folding. The nucleus is composed of a set of adjacent residues (an alpha-helix), stabilized by long-range interactions that are formed as the rest of the protein collapses around it. A simple analysis of the optimization of the rate of protein folding predicts that rates are highest when the denatured state has little residual structure under physiological conditions and no intermediates accumulate. This implies that any potential nucleation site that is composed mainly of adjacent residues should be just weakly populated in the denatured state and become structured only in a high-energy intermediate or transition state when it is stabilized by interactions elsewhere in the protein. Hierarchical mechanisms of folding in which stable elements of structure accrete are unfavorable. The nucleation-condensation mechanism of CI2 fulfills the criteria for fast folding. On the other hand, stable intermediates do form in the folding of more complex proteins, and this may be an unavoidable consequence of increasing size and nucleation at more than one site.

Numerical modelling of single-layer folding: clarification of an issue regarding the possible effect of computer codes and the influence of initial irregularities

The influence of initial perturbation geometry and material propel-ties on final fold geometry has been investigated using finite-difference (FLAC) and finite-element (MARC) numerical models. Previous studies using these two different codes reported very different folding behaviour although the material properties, boundary conditions and initial perturbation geometries were similar. The current results establish that the discrepancy was not due to the different computer codes but due to the different strain rates employed in the two previous studies (i.e. 10(-6) s(-1) in the FLAC models and 10(-14) s(-1) in the MARC models). As a result, different parts of the elasto-viscous rheological field were bring investigated. For the same material properties, strain rate and boundary conditions, the present results using the two different codes are consistent. A transition in Folding behaviour, from a situation where the geometry of initial perturbation determines final fold shape to a situation where material properties control the final geometry, is produced using both models. This transition takes place with increasing strain rate, decreasing elastic moduli or increasing viscosity (reflecting in each case the increasing influence of the elastic component in the Maxwell elastoviscous rheology). The transition described here is mechanically feasible but is associated with very high stresses in the competent layer (on the order of GPa), which is improbable under natural conditions. (C) 2000 Elsevier Science Ltd. All rights reserved.

Translation termination and protein folding pathway genes are not correlated in gastric cancer

Background: The eukaryotic release factor 3 (eRF3) has been shown to affect both tubulin and actin cytoskeleton, suggesting a role in cytoskeleton assembly, mitotic spindle formation and chromosome segregation. Also, direct interactions between eRF3 and subunits of the cytosolic chaperonin CCT have been described. Moreover, both eRF3a and CCT subunits have been described to be up-regulated in cancer tissues. Our aim was to evaluate the hypothesis that eRF3 expression levels are correlated with the expression of genes encoding proteins involved in the tubulin folding pathways. Methods: Relative expression levels of eRF1, eRF3a/GSPT1, PFDN4, CCT2, CCT4, and TBCA genes in tumour samples relative to their adjacent normal tissues were investigated using real time-polymerase chain reaction in 20 gastric cancer patients. Results: The expression levels of eRF3a/GSPT1 were not correlated with the expression levels of the other genes studied. However, significant correlations were detected between the other genes, both within intestinal and diffuse type tumours. Conclusions: eRF3a/GSPT1 expression at the mRNA level is independent from both cell translation rates and from the expression of the genes involved in tubulin-folding pathways. The differences in the patterns of expression of the genes studied support the hypothesis of genetically independent pathways in the origin of intestinal and diffuse type gastric tumours.

Coordinate and time-dependent diffusion dynamics in protein folding

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)