The Role of Solvent in Protein Folding and in Aggregation

We discuss features of the effect of solvent on protein folding andaggregation, highlighting the physics related to the particulate nature and the peculiar structure of the aqueous solvent, and the biological significance of interactions between solvent and proteins. To this purpose we use a generalized energy landscape of extended dimensionality. A closer look at the properties of solvent induced interactions and forces proves useful for understanding the physical grounds of `ad hoc' interactions and for devising realistic ways of accounting for solvent effects. The solvent has long been known to be a crucially important part of biological systems, and times appear mature for it to be adequately accounted for in the protein folding problem. Use of the extended dimensionality energy landscape helpseliciting the possibility of coupling among conformational changes and aggregation, such as proved by experimental data in the literature.     

Heterogeneous Packing in the Folding/Unfolding Intermediate State of Bitter Gourd Trypsin Inhibitor

The conformation and dynamics of a protein are essential in characterizing the protein folding/unfolding intermediate state. They are closely involved in the packing and site-specific interactions of peptide elements to build and stabilize the tertiary structure of the protein. In this study, it was confirmed that trypsin inhibitor obtained from seeds of bitter gourd (BGTI) adopted a peculiar but plausible conformation and dynamics in the unfolding intermediate state. The fluorescence spectrum of one of two tryptophan residues of BGTI, Trp9, shifted to the blue side in the presence of 2–3 M guanidine hydrochloride, although the other, Trp54, did not show this spectral shift. At the same time, the motional freedom of Trp9 revealed by a time-resolved fluorescence study decreased, suggesting that the segmental motion of this residue was more restricted. These results indicate that BGTI takes such a conformation state that the hydrophobic core and loop domains arranging Trp9 and Trp54 respectively are heterogeneously packed in the unfolding intermediate state.     

Folding Nuclei in Proteins

When a protein folds or unfolds, it passes through many half-folded microstates. Only a few of them can accumulate and be seen experimentally, and this happens only when the folding (or unfolding) occurs far from the point of thermodynamic equilibrium between the native and denatured states. The universal features of folding, though, are observed in the vicinity of the equilibrium point. Here the two-state transition proceeds without any accumulation of metastable intermediates, and only the transition state (folding nucleus) is outlined by its key influence on the folding/unfolding kinetics. This review covers recent experimental and theoretical studies of folding nuclei.     

Evolutionary Aspects of Protein Structure and Folding

Four basic stages of evolution of protein structure are described, basing on recent work of the authors aimed specifically to reconstruct the earliest events in the protein evolution. According to this reconstruction, the initial stage of short peptides comprising, probably, only a few amino acid residues had been followed by formation of closed loops of 25–30 residues, which corresponds to the polymer-statistically optimal ring closure size for mixed polypeptide chains. The next stage involved fusion of relatively small linear genes and formation of protein structures consisting of several closed loops of a nearly standard size, with 4–6 loops (100–200 amino acid residues) in a typical protein fold. The last, modern stage began with combinatorial fusion of the presumably circular 300–600 bp DNA units and, accordingly, formation of multidomain proteins.     

Folding cooperativity in a three-stranded beta-sheet model

The thermodynamic behavior of a previously designed three-stranded beta-sheet was studied via several microseconds of standard and replica exchange molecular dynamics simulations. The system is shown to populate at least four thermodynamic minima, including two partially folded states in which only a single hairpin is formed. Simulated melting curves show different profiles for the C and N-terminal hairpins, consistent with differences in secondary structure content in published NMR and CD/FTIR measurements, which probed different regions of the chain. Individual beta-hairpins that comprise the three-stranded beta-sheet are observed to form cooperatively. Partial folding cooperativity between the component hairpins is observed, and good agreement between calculated and experimental values quantifying this cooperativity is obtained when similar analysis techniques are used. However, the structural detail in the ensemble of conformations sampled in the simulations permits a more direct analysis of this cooperativity than has been performed on the basis of experimental data. The results indicate the actual folding cooperativity perpendicular to strand direction is significantly larger than the lower bound obtained previously.     

Polymer Uncrossing and Knotting in Protein Folding,and Their Role in Minimal Folding Pathways

We introduce a method for calculating the extent to which chain non-crossing is important in the most efficient, optimal trajectories or pathways for a protein to fold. This involves recording all unphysical crossing events of a ghost chain, and calculating the minimal uncrossing cost that would have been required to avoid such events. A depth-first tree search algorithm is applied to find minimal transformations to fold , , , and knotted proteins. In all cases, the extra uncrossing/non-crossing distance is a small fraction of the total distance travelled by a ghost chain. Different structural classes may be distinguished by the amount of extra uncrossing distance, and the effectiveness of such discrimination is compared with other order parameters. It was seen that non-crossing distance over chain length provided the best discrimination between structural and kinetic classes. The scaling of non-crossing distance with chain length implies an inevitable crossover to entanglement-dominated folding mechanisms for sufficiently long chains. We further quantify the minimal folding pathways by collecting the sequence of uncrossing moves, which generally involve leg, loop, and elbow-like uncrossing moves, and rendering the collection of these moves over the unfolded ensemble as a multiple-transformation “alignment”. The consensus minimal pathway is constructed and shown schematically for representative cases of an , , and knotted protein. An overlap parameter is defined between pathways; we find that proteins have minimal overlap indicating diverse folding pathways, knotted proteins are highly constrained to follow a dominant pathway, and proteins are somewhere in between. Thus we have shown how topological chain constraints can induce dominant pathway mechanisms in protein folding.     

Folding of viral envelope glycoproteins in the endoplasmic reticulum

Viral glycoproteins fold and oligomerize in the endoplasmic reticulum of the host cell. They employ the cellular machinery and receive assistance from cellular folding factors. During the folding process, they are retained in the compartment and their structural quality is checked by the quality control system of the endoplasmic reticulum. A special characteristic that distinguishes viral fusion proteins from most cellular proteins is the extensive conformational change they undergo during fusion of the viral and cellular membrane. Many viral proteins fold in conjunction with and dependent on a viral partner protein, sometimes even synthesized from the same mRNA. Relevant for folding is that viral glycoproteins from the same or related virus families may consist of overlapping sets of domain modules. The consequences of these features for viral protein folding are at the heart of this review.     

Loop Folds in Proteins and Evolutionary Conservation of Folding Nuclei

Abstract

We show that loops of close contacts involving hydrophobic residues are important in protein folding. Contrary to Berezovsky and Trifonov (J. Biomol. Struct. Dyn. 20, 5–6, 2002) the loops important in protein folding usually are much larger in size than 23–31 residues, being instead comparable to the size of the protein for single domain proteins. Additionally what is important are not single loop contacts, but a highly interconnected network of such loop contacts, which provides extra stability to a protein fold and which leads to their conservation in evolution.     

Glycoprotein Folding in the Endoplasmic Reticulum

Our understanding of eukaryotic protein folding in the endoplasmic reticulum has increased enormously over the last 5 years. In this review, we summarize some of the major research themes that have captivated researchers in this field during the last years of the 20th century. We follow the path of a typical protein as it emerges from the ribosome and enters the reticular environment. While many of these events are shared between different polypeptide chains, we highlight some of the numerous differences between proteins, between cell types, and between the chaperones utilized by different ER glycopro-teins. Finally, we consider the likely advances in this field as the new century unfolds and we address the prospect of a unified understanding of how protein folding, degradation, and translation are coordinated within a cell.     

Recent Advances in GFP Folding Reporter and Split-GFP Solubility Reporter Technologies. Application to Improving the Folding and Solubility of Recalcitrant Proteins from Mycobacterium tuberculosis

We have improved our green fluorescent protein (GFP) folding reporter technology [Waldo et al., (1999) Nat. Biotechnol. 17, 691–695] to evolve recalcitrant proteins from Mycobacterium tuberculosis. The target protein is inserted into the scaffolding of the GFP, eliminating false-positive artifacts caused by expression of truncated protein variants from internal cryptic ribosome binding sites in the target RNA. In parallel, we have developed a new quantitative fluorescent protein tagging and detection system based on micro-domains of GFP. This split-GFP system, which works both in vivo and in vitro, is amenable to high-throughput assays of protein expression and solubility [Cabantous et al., (2005) Nat. Biotechnol. 23, 102–107]. Together, the GFP folding reporter and split-GFP technologies offer a comprehensive system for manipulating and improving protein folding and solubility.     

Folding of the yeast prion protein Ure2: kinetic evidence for folding and unfolding intermediates.

The Saccharomyces cerevisiae non-Mendelian factor [URE3] propagates by a prion-like mechanism, involving aggregation of the chromosomally encoded protein Ure2. The N-terminal prion domain (PrD) of Ure2 is required for prion activity in vivo and amyloid formation in vitro. However, the molecular mechanism of the prion-like activity remains obscure. Here we measure the kinetics of folding of Ure2 and two N-terminal variants that lack all or part of the PrD. The kinetic folding behaviour of the three proteins is identical, indicating that the PrD does not change the stability, rates of folding or folding pathway of Ure2. Both unfolding and refolding kinetics are multiphasic. An intermediate is populated during unfolding at high denaturant concentrations resulting in the appearance of an unfolding burst phase and roll-over in the denaturant dependence of the unfolding rate constants. During refolding the appearance of a burst phase indicates formation of an intermediate during the dead-time of stopped-flow mixing. A further fast phase shows second-order kinetics, indicating formation of a dimeric intermediate. Regain of native-like fluorescence displays a distinct lag due to population of this on-pathway dimeric intermediate. Double-jump experiments indicate that isomerisation of Pro166, which is cis in the native state, occurs late in refolding after regain of native-like fluorescence. During protein refolding there is kinetic partitioning between productive folding via the dimeric intermediate and a non-productive side reaction via an aggregation prone monomeric intermediate. In the light of this and other studies, schemes for folding, aggregation and prion formation are proposed.     

Scaling of Folding Properties in Go Models of Proteins

Insights about scaling of folding properties of proteins are obtained bystudying folding in heteropolymers described by Go-like Hamiltonians. Bothlattice and continuum space models are considered. In the latter case, themonomer-monomer interactions correspond to the Lennard-Jones potential.Several statistical ensembles of the two- and three-dimensional targetnative conformations are considered. Among them are maximally compactconformations which are confined to a lattice and those which are obtainedeither through quenching or annealing of homopolymers to their compactlocal energy minima. Characteristic folding times are found to grow aspower laws with the system size. The corresponding exponents are notuniversal. The size related deterioration of foldability is found to beconsistent with the scaling behavior of the characteristic temperatures:asymptotically, the folding temperature becomes much lower than thetemperature at which glassy kinetics become important. The helicalconformations are found to have the lowest overall scaling exponent andthe best foldability among the classes of conformations studied. Thescaling properties of the Go-like models of the protein conformationsstored in the Protein Data Bank suggest that proteins are not optimizedkinetically.     

Hydrogen Bonds,Hydrophobicity Forces and the Character of the Collapse Transition

We study the thermodynamic behavior of a model protein with 54 amino acidsthat is designed to form a three-helix bundle in its native state. The model contains three types of amino acids and five to six atoms per amino acid, and has the Ramachandran torsion angles as its only degrees of freedom.The force field is based on hydrogen bonds and effective hydrophobicity forces. We study how the character of the collapse transition depends on the strengths of these forces. For a suitable choice of these two parameters, it is found that the collapse transition is first-order-like and coincides with the folding transition. Also shown is that the corresponding one- and two-helix segments make less stable secondary structure than the three-helix sequence.     

A Comparison of Experimental and Computational Methods for Mapping the Interactions Present in the Transition State for Folding of FKBP12

The folding pathway of FKBP12, a 107 residue / protein, has been characterised in detail using a combination of experimental and computational techniques. FKBP12 follows a two-state model of folding in which only the denatured and native states are significantly populated; no intermediate states are detected. The refolding rate constant in water is 4 s^-1 at 25 °C. Two different experimental strategies were employed for studying the transition state for folding. In the first case, a non-mutagenic approach was used and the unfolding and refolding of the wild-type protein measured as a function of experimental conditions such as temperature, denaturant, ligand and trifluoroethanol (TFE) concentration. These data suggest a compact transition state relative to the unfolded state with some 70% of the surface area buried. The ligand-binding site, whichis mainly formed by two long loops, is largely unstructured in the transition state. TFE experiments suggest that the -helix may be formed in the transition state. The second experimental approach involved using protein engineering techniques with -value analysis. Residue-specific information on the structure and energetics of the transition state can be obtained by this method. 34 mutations were made at sites throughout the protein to probe the extent of secondary and tertiary structure in the transition state. In contrast to some other proteins of this size, no element of structure is fully formed in the transition state, instead, the transition state is similar to that found for smaller, single-domain proteins, such as chymotrypsin inhibitor 2 and the SH3 domainfrom -spectrin. For FKBP12, the central three strands of the -sheet (2, 4 and 5), comprise the most structured region of the transition state. In particular Val 101, which is one of the most highly buried residues and located in the middle of the central -strand,makes approximately 60% of its native interactions. The outer -strands, and the ends of the central -strands are formed to a lesser degree. The short -helix is largely unstructured in the transition state as are the loops. The data are consistent with a nucleation-condensation model of folding, the nucleus of which is formed by side chains within -strands 2, 4 and 5 and the C-terminus of the -helix. These residues are distant in the primary sequence, demonstrating the importance of tertiary interactions in the transition state. High-temperature molecular dynamic simulations on the unfoldingpathway of FKBP12 are in good agreement with the experimental results.     

Folding kinetics of an alpha helical membrane protein in phospholipid bilayer vesicles

We report a detailed kinetic study of the folding of an alpha-helical membrane protein in a lipid bilayer environment. SDS denatured bacteriorhodopsin was folded directly into phosphatidylcholine lipid vesicles by stopped-flow mixing. The folding kinetics were monitored with millisecond time resolution by time-resolving changes in protein fluorescence as well as in the absorption of the retinal chromophore. The kinetics were similar to those previously reported for folding bacteriorhodopsin in detergent or lipid micelles, except for the presence of an additional apoprotein intermediate. We suggest this intermediate is a result of the greater internal two-dimensional pressure present in these lipid vesicles as compared to micelles. These results lay the groundwork for future studies aimed at understanding the mechanistic origin of the effect of lipid bilayer properties on protein folding. Furthermore, the use of biologically relevant phosphatidylcholine lipids, together with a straightforward rapid mixing process to initiate the folding reaction, means the method is generally applicable, and thus paves the way for an improved understanding of the in vitro folding of transmembrane alpha-helical proteins.     

Enumerating Designing Sequences in the HP Model

The hydrophobic/polar HP model on the square lattice has been widely used toinvestigate basics of protein folding. In the cases where all designing sequences (sequences with unique ground states) were enumerated without restrictions on the number of contacts, the upper limit on the chain length N has been 18–20 because of the rapid exponential growth of thenumbers of conformations and sequences. We show how a few optimizations push this limit by about 5 units. Based on these calculations, we study the statistical distribution of hydrophobicity along designing sequences. We find that the average number of hydrophobic and polar clumps along the chains is larger for designing sequences than for random ones, which is in agreement with earlier findings for N 18 and with results for real enzymes. We also show that this deviation from randomness disappears if the calculations are restricted to maximally compact structures.     

Thermodynamic analysis of protein unfolding in aqueous solutions as a multisite reaction of protein with water and solute molecules

Thermal unfolding of ribonuclease A, lysozyme, and chymotrypsinogen A was analyzed as a multisite reaction of a protein molecule with water and solute molecules. The protein unfolding process in various solutions of sugars and denaturants was described well by the van't Hoff equation. The reciprocal form of the Wyman-Tanford equation, which describes the unfolded-to-folded protein ratio as a function of water activity, was successfully applied to obtain a good linear relationship. From this analysis, the role of water activity on protein stability was clearly explained and the contributions of hydration and solute binding to protein molecule were separately discussed in protein unfolding. General solution for the free energy of protein stability was obtained as a simple function of solute concentration.     

The solvation interface is a determining factor in peptide conformational preferences

The 21 residue polyalanine-based F(s) peptide was studied using thousands of long, explicit solvent, atomistic molecular dynamics simulations that reached equilibrium at the ensemble level. Peptide conformational preference as a function of hydrophobicity was examined using a spectrum of explicit solvent models, and the peptide length-dependence of the hydrophilic and hydrophobic components of solvent-accessible surface area for several ideal conformational types was considered. Our results demonstrate how the character of the solvation interface induces several conformational preferences, including a decrease in mean helical content with increased hydrophilicity, which occurs predominantly through reduced nucleation tendency and, to a lesser extent, destabilization of helical propagation. Interestingly, an opposing effect occurs through increased propensity for 3(10)-helix conformations, as well as increased polyproline structure. Our observations provide a framework for understanding previous reports of conformational preferences in polyalanine-based peptides including (i) terminal 3(10)-helix prominence, (ii) low pi-helix propensity, (iii) increased polyproline conformations in short and unfolded peptides, and (iv) membrane helix stability in the presence and absence of water. These observations provide physical insight into the role of water in peptide conformational equilibria at the atomic level, and expand our view of the complexity of even the most "simple" of biopolymers. Whereas previous studies have focused predominantly on hydrophobic effects with respect to tertiary structure, this work highlights the need for consideration of such effects at the secondary structural level.