Protein folding: optimized sequences obtained by simulated breeding in a minimalist model

For a minimalist model of protein folding, which we introduced recently, we investigate various methods to obtain folding sequences. A detailed study of random sequences shows that, for this model, such sequences usually do not fold to their ground states during simulations. Straight-forward techniques for the construction of folding sequences, based solely on the target structure, fail. We describe in detail an optimization algorithm, based on genetic algorithms, for the "simulated breeding" of folding sequences in this model. We find that, for any target structure studied, there is not only a single folding sequence but a patch of sequences in sequence space that fold to this structure. In addition, we show that, much as in real proteins, nonhomologous sequences may fold to the same target structure.     

How evolution makes proteins fold quickly

Sequences of fast-folding model proteins (48 residues long on a cubic lattice) were generated by an evolution-like selection toward fast folding. We find that fast-folding proteins exhibit a specific folding mechanism in which all transition state conformations share a smaller subset of common contacts (folding nucleus). Acceleration of folding was accompanied by dramatic strengthening of interactions in the folding nucleus whereas average energy of nonnucleus interactions remained largely unchanged. Furthermore, the residues involved in the nucleus are the most conserved ones within families of evolved sequences. Our results imply that for each protein structure there is a small number of conserved positions that are key determinants of fast folding into that structure. This conjecture was tested on two protein superfamilies: the first having the classical monophosphate binding fold (CMBF; 98 families) and the second having type-III repeat fold (47 families). For each superfamily, we discovered a few positions that exhibit very strong and statistically significant "conservatism of conservatism"-amino acids in those positions are conserved within every family whereas the actual types of amino acids varied from family to family. Those amino acids are in spatial contact with each other. The experimental data of Serrano and coworkers [Lopez-Hernandez, E. & Serrano, L. (1996) Fold. Des. (London) 1, 43-55]. for one of the proteins of the CMBF superfamily (CheY) show that residues identified this way indeed belong to the folding nucleus. Further analysis revealed deep connections between nucleation in CMBF proteins and their function.     

The early stage of folding of villin headpiece subdomain observed in a 200-nanosecond fully solvated molecular dynamics simulation

A new approach in implementing classical molecular dynamics simulation for parallel computers has enabled a simulation to be carried out on a protein with explicit representation of water an order of magnitude longer than previously reported and will soon enable such simulations to be carried into the microsecond time range. We have used this approach to study the folding of the villin headpiece subdomain, a 36-residue small protein consisting of three helices, from an unfolded structure to a molten globule state, which has a number of features of the native structure. The time development of the solvation free energy, the radius of gyration, and the mainchain rms difference from the native NMR structure showed that the process can be seen as a 60-nsec "burst" phase followed by a slow "conformational readjustment" phase. We found that the burial of the hydrophobic surface dominated the early phase of the folding process and appeared to be the primary driving force of the reduction in the radius of gyration in that phase.     

An amino acid as a cofactor for a catalytic polynucleotide

Natural ribozymes require metal ion cofactors that aid both in structural folding and in chemical catalysis. In contrast, many protein enzymes produce dramatic rate enhancements using only the chemical groups that are supplied by their constituent amino acids. This fact is widely viewed as the most important feature that makes protein a superior polymer for the construction of biological catalysts. Herein we report the in vitro selection of a catalytic DNA that uses histidine as an active component for an RNA cleavage reaction. An optimized deoxyribozyme from this selection requires L-histidine or a closely related analog to catalyze RNA phosphoester cleavage, producing a rate enhancement of approximately 1-million-fold over the rate of substrate cleavage in the absence of enzyme. Kinetic analysis indicates that a DNA-histidine complex may perform a reaction that is analogous to the first step of the proposed catalytic mechanism of RNase A, in which the imidazole group of histidine serves as a general base catalyst. Similarly, ribozymes of the "RNA world" may have used amino acids and other small organic cofactors to expand their otherwise limited catalytic potential.     

Molecular mechanisms for cooperative folding of proteins

The folding of single-domain globular proteins exhibits the character of first-order or two-state thermodynamics. The origin of such high cooperativity in relatively small polymer systems is still not well understood. Recently, the statistical mechanics of protein folding has been studied extensively with simple protein models such as short cubic-lattice chains with contact-based interactions. While many valuable insights about protein folding were gained with such models, some concerns have also arisen, viz. that they lack the character of protein backbones whose interactions would limit the folding patterns of proteins. Here, a comparative study of the conventional cubic-lattice chain model and a fine-grained more realistic lattice protein model with both backbone and side-chain interactions is carried out. It is found that, even though both types of models exhibit a cooperative two-state folding transition to the native structure with optimized force fields, the character and origin of cooperativity of the two models are different. In the simple contact-based model, the free-energy barrier occurs at the low end of the energy scale, and the cooperativity arises from a concerted formation of native contacts among many residues in a compact state. In the other more complicated model, the free-energy barrier occurs in the intermediate energy region, and the folding cooperativity arises from collective orientational arrangements of locally structured units in semi-open conformational states. On the basis of these results, two limiting molecular mechanisms for protein folding emerge, which can be used for analyzing the folding process of real proteins.     

Salt effects on hydrophobic interaction and charge screening in the folding of a negatively charged peptide to a coiled coil (leucine zipper)

The stability of a coiled coil or leucine zipper is controlled by hydrophobic interactions and electrostatic forces between the constituent helices. We have designed a 30-residue peptide with the repeating seven-residue pattern of a coiled coil, (abcdefg)n, and with Glu in positions e and g of each heptad. The glutamate side chains prevented folding at pH values above 6 because of electrostatic repulsion across the helix dimer interface as well as within the individual helices. Protonation of the carboxylates changed the conformation from a random coil monomer to a coiled coil dimer. Folding at alkaline pH where the peptide had a net charge of -7e was promoted by the addition of salts. The nature of the charge screening cation was less important than that of the anion. The high salt concentrations (>1 M) necessary to induce folding indicated that the salt-induced folding resulted from alterations in the protein-water interaction. Folding was promoted by the kosmotropic anions sulfate and fluoride and to a lesser extent by the weak kosmotrope formate, whereas chloride and the strong chaotrope perchlorate were ineffective. Kosmotropes are excluded from the protein surface, which is preferentially hydrated, and this promotes folding by strengthening hydrophobic interactions at the coiled coil interface. Although charge neutralization also contributed to folding, it was effective only when the screening cation was partnered by a good kosmotropic anion. Folding conformed to a two-state transition from random coil monomer to coiled coil dimer and was enthalpy driven and characterized by a change in the heat capacity of unfolding of 3.9 +/- 1.2 kJ mol-1 K-1. The rate of folding was analyzed by fluorescence stopped-flow measurements. Folding occurred in a biphasic reaction in which the rapid formation of an initial dimer (kf = 2 x 10(7) M-1 s-1) was followed by an equally rapid concentration-independent rearrangement to the folded dimer (k > 100 s-1).     

Artificial cell containing superoxide dismutase--selection of folding aids for stabilisation of SOD

Superoxide dismutase (abbreviated as SOD) has been vigorously studied in the fields of radical chemistry and related life science. One of practical problems is how to keep its activity in certain adverse conditions causing denaturation. Artificial cell containing SOD can be prepared by polymer encapsulation or nanocapsulation which has been found to be effective to improve the stability of SOD. For construction of an ideal artificial cell system, some folding aids or aggregation inhibitors were utilised to enhance SOD stability. In this study, three groups of biopolymers are selected as folding aids or aggregation inhibitors for stabilisation of SOD, i.e. albumin, carbohydrates and glycoproteins. Results indicate that the thermostability of SOD is affected by different sort of albumin while some carbohydrates such as cyclodextrins are found to be able to enhance SOD stability. In addition, it is firstly found that selected glycoproteins such as alpha-macroglobulin and ovalbumin are several types of effective folding aids for stabilisation of SOD. They can protect SOD against denaturation even at very high temperature(over 100 degrees C). The stability was tested by the measurement of SOD activity loss using autooxidation method in different adverse conditions such as high temperature, extreme pH medium, proteolytic hydrolysis and long shelf life storage. The possible stabilisation mechanism of using cyclodextrins and glycoproteins as folding aids were discussed.     

Slow alpha helix formation during folding of a membrane protein

Very little is known about the folding of proteins within biological membranes. A "two-stage" model has been proposed on thermodynamic grounds for the folding of alpha helical, integral membrane proteins, the first stage of which involves formation of transmembrane alpha helices that are proposed to behave as autonomous folding domains. Here, we investigate alpha helix formation in bacteriorhodopsin and present a time-resolved circular dichroism study of the slow in vitro folding of this protein. We show that, although some of the protein's alpha helices form early, a significant part of the protein's secondary structure appears to form late in the folding process. Over 30 amino acids, equivalent to at least one of bacteriorhodopsin's seven transmembrane segments, slowly fold from disordered structures to alpha helices with an apparent rate constant of about 0.012 s-1 at pH 6 or 0.0077 s-1 at pH 8. This is a rate-limiting step in protein folding, which is dependent on the pH and the composition of the lipid bilayer.     

GroEL-GroES-mediated protein folding requires an intact central cavity

The chaperonin GroEL is an oligomeric double ring structure that, together with the cochaperonin GroES, assists protein folding. Biochemical analyses indicate that folding occurs in a cis ternary complex in which substrate is sequestered within the GroEL central cavity underneath GroES. Recently, however, studies of GroEL "minichaperones" containing only the apical substrate binding subdomain have questioned the functional importance of substrate encapsulation within GroEL-GroES complexes. Minichaperones were reported to assist folding despite the fact that they are monomeric and therefore cannot form a central cavity. Here we compare directly the folding activity of minichaperones with that of the full GroEL-GroES system. In agreement with earlier studies, minichaperones assist folding of some proteins. However, this effect is observed only under conditions where substantial spontaneous folding is also observed and is indistinguishable from that resulting from addition of the nonchaperone protein alpha-casein. By contrast, the full GroE system efficiently promotes folding of several substrates under conditions where essentially no spontaneous folding is observed. These data argue that the full GroEL folding activity requires the intact GroEL-GroES complex, and in light of previous studies, underscore the importance of substrate encapsulation for providing a folding environment distinct from the bulk solution.     

Recombination of protein domains facilitated by co-translational folding in eukaryotes

The evolution of complex genomes requires that new combinations of pre-existing protein domains successfully fold into modular polypeptides. During eukaryotic translation model two-domain polypeptides fold efficiently by sequential and co-translational folding of their domains. In contrast, folding of the same proteins in Escherichia coli is posttranslational, and leads to intramolecular misfolding of concurrently folding domains. Sequential domain folding in eukaryotes may have been critical in the evolution of modular polypeptides, by increasing the probability that random gene-fusion events resulted in immediately foldable protein structures.     

Refolding rate of stability-enhanced cytochrome c is independent of thermodynamic driving force

N52I iso-2 cytochrome c is a variant of yeast iso-2 cytochrome c in which asparagine substitutes for isoleucine 52 in an alpha helical segment composed of residues 49-56. The N52I substitution results in a significant increase in both stability and cooperativity of equilibrium unfolding, and acts as a "global suppressor" of destabilizing mutations. The equilibrium m-value for denaturant-induced unfolding of N52I iso-2 increases by 30%, a surprisingly large amount for a single residue substitution. The folding/unfolding kinetics for N52I iso-2 have been measured by stopped-flow mixing and by manual mixing, and are compared to the kinetics of folding/unfolding of wild-type protein, iso-2 cytochrome c. The results show that the observable folding rate and the guanidine hydrochloride dependence of the folding rate are the same for iso-2 and N52I iso-2, despite the greater thermodynamic stability of N52I iso-2. Thus, there is no linear free-energy relationship between mutation-induced changes in stability and observable refolding rates. However, for N52I iso-2 the unfolding rate is slower and the guanidine hydrochloride dependence of the unfolding rate is smaller than for iso-2. The differences in the denaturant dependence of the unfolding rates suggest that the N52I substitution decreases the change in the solvent accessible hydrophobic surface between the native state and the transition state. Two aspects of the results are inconsistent with a two-state folding/unfolding mechanism and imply the presence of folding intermediates: (1) observable refolding rate constants calculated from the two-state mechanism by combining equilibrium data and unfolding rate measurements deviate from the observed refolding rate constants; (2) kinetically unresolved signal changes ("burst phase") are observed for both N52I iso-2 and iso-2 refolding. The "burst phase" amplitude is larger for N52I iso-2 than for iso-2, suggesting that the intermediates formed during the "burst phase" are stabilized by the N52I substitution.     

Thermodynamics of a reverse turn motif. Solvent effects and side-chain packing

The linear pentapeptide, Ala-Tyr-cis-Pro-Tyr-Asp-NMA (AYPYD) is known to have a significant population of type VI turn conformers in aqueous solvent. We have carried out theoretical studies of the conformational energetics of this peptide using a potential of mean force (PMF) consisting of the AMBER/OPLS empirical potential energy function, a macroscopic electrostatic model of polar solvation, and a surface area-based model of non-polar solvation. Conformers were taken from molecular dynamics simulations reported elsewhere, or generated by a random search method reported here. The chain entropy of folding was calculated by a systematic search of accessible dihedral angle space. The intra-peptide component was found to strongly favor folding and was nearly cancelled by the polar solvation term which disfavored folding. The non-polar solvation term had little effect. Fluctuations about the average value of the PMF were small and in accord with estimates from a simple harmonic model. When applied to conformers generated by a random search, the PMF selected a conformer close to the NMR-determined structure as the lowest energy conformer. The conformer with the second-lowest energy was extended, but was found to fold rapidly to the turn state in a subsequent molecular dynamics study, and may be an important state on the folding-unfolding pathway. Averages of the PMF were combined with the entropy estimates to provide an estimate of the free energy of folding that is in reasonable agreement with experimental results. In terms of the interplay between backbone electrostatic interactions and the packing of apolar side-chains, this peptide provides a model for the energetics of protein folding, and therefore makes a useful test case for calculations.     

Brownian dynamics simulations of protein folding: access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion-collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion-collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.     

Formation of electrostatic interactions on the protein-folding pathway

We describe a novel method of obtaining information about the structures of transient conformations on the folding pathway from their ionization equilibria: the H+ -titration behavior of a protein residue is determined in detail by its environment. We follow the consolidation of electrostatic interactions in the folding process by comparing the acid-titration behavior of four conformations on the folding pathway of barnase: the denatured state (D); the folding intermediate (I); the major transition state(+); and the native state (N) in the scheme D <==>I<==>(+)<==)N. The results show that strong electrostatic interactions are present in the major transition state: some of its carboxylate groups display the highly anomalous pKA values of <2 that are found in N. However, the network of ionic surface interactions is not formed in (+), and the overall protection of titrating residues is weakened. The results are consistent with the transition state being an expanded form of the native state, with a weakened but poorly hydrated core and a loosened periphery. The surface residues in such an expanded conformation are, on average, farther apart than are those in the center of the molecule. The results concerning the folding intermediate are less clear cut. We show that the interpretation of kinetic data relating to folding intermediates depends critically on assumptions about their equilibrium with other denatured states. We have, however, characterized the pH and ionic strength dependence of an apparent stability of I, using the deviation from two-state folding behavior, which can be used to investigate electrostatic properties of folding intermediates from a variety of mechanisms. In general, the data imply that I is somewhat similar to (+). Apparently odd titration properties of I are investigated further in the accompanying paper [Oliveberg, M., & Fersht, A. (1996) Biochemistry 35, 2738-2749]. The approach in this study may be of particular use in testing theoretical results since the relationship between H+ -titration properties and protein structure can be treated by classical electrostatics.     

Kinetics of peptide folding: computer simulations of SYPFDV and peptide variants in water

The folding of Ser-Tyr-Pro-Phe-Asp-Val (SYPFDV), and sequence variants of this peptide (SYPYD and SYPFD) are studied computationally in an explicit water environment. An atomically detailed model of the peptide is embedded in a sphere of TIP3P water molecules and its optimal structure is computed by simulated annealing. At distances from the peptide that are beyond a few solvation shells, a continuum solvent model is employed. The simulations are performed using a mean field approach that enhances the efficiency of sampling peptide conformations. The computations predict a small number of conformations as plausible folded structures. All have a type VI turn conformation for the peptide backbone, similar to that found using NMR. However, some of the structures differ from the experimentally proposed ones in the packing of the proline ring with the aromatic residues. The second most populated structure has, in addition to a correctly folded backbone, the same hydrophobic packing as the conformation measured by NMR. Our simulations suggest a kinetic mechanism that consists of three separate stages. The time-scales associated with these stages are distinct and depend differently on temperature. Electrostatic interactions play an initial role in guiding the peptide chain to a roughly correct structure as measured by the end-to-end distance. At the same time or later the backbone torsions rearrange due to local tendency of the proline ring to form a turn: this step depends on solvation forces and is helped by loose hydrophobic interactions. In the final step, hydrophobic residues pack against each other. We also show the existence of an off the pathway intermediate, suggesting that even in the folding of a small peptide "misfolded" structures can form. The simulations clearly show that parallel folding paths are involved. Our findings suggest that the process of peptide folding shares many of the features expected for the significantly larger protein molecules.     

Nuclear magnetic resonance shows asymmetric loss of triple helix in peptides modeling a collagen mutation in brittle bone disease

To investigate a human folding disease, NMR studies were carried out on collagen-like peptides to define the structural consequences of a single amino acid change found in patients with osteogenesis imperfecta (OI), a disease characterized by fragile bones. One peptide included a normal collagen sequence, while a second peptide included a Gly --> Ser substitution as found in a nonlethal case of OI. Residue specific internal dynamics and conformational studies indicate that the normal collagen-like sequence forms a triple helix which is rigid along its entire length. The introduction of a Gly --> Ser substitution induces an asymmetric disruption of the uniform triple helix. While the C-terminal end of the peptide retains the triple helix, the Ser substitution site and residues N-terminal to it exhibit the mobility of a random chain. This equilibrium state indicates that a Gly substitution can terminate the C to N propagation of the triple helix and suggests that renucleation is required for folding to continue. Defective folding has been implicated in brittle bone disease, and these results begin to characterize the folding process in OI collagens. OI collagen studies may also provide insights about defective protein folding, assembly, and aggregation in other human diseases.     

Pluronic copolymer liquid crystals: unique, replaceable media for capillary gel electrophoresis

Liquid crystalline solutions of Pluronic copolymers are versatile alternatives to solutions of entangled, random coil polymers as replaceable media for capillary gel electrophoresis (CGE). Pluronic copolymers are tri-block polymers of poly(ethylene oxide) [(EO)x] and poly(propylene oxide) [(PO)y] with the general formula (EO)x(PO)y(EO)x. Large micelles form in aqueous solutions in which central, hydrophobic cores of (PO)y segments are surrounded by "brushes" of hydrated (EO)x tails. Solutions of Pluronic F127 (BASF Performance Chemicals) in a concentration range of about 18-30% are liquids at refrigerator temperatures (< or = 5 degrees C) and are easily introduced into capillaries. A self-supporting, gel-like liquid crystalline phase is formed as the temperature is raised to > or = 20 degrees C. This liquid crystalline phase consists of spherical micelles with diameters of 17-18 nm which pack with local cubic symmetry. CGE in Pluronic F127 liquid crystals separates species within several chemical classes as varied as nucleoside monophosphates and organic dyes, oligonucleotides of 4-60 nucleotides, DNA fragments of 50-3000 base pairs (bp), and supercoiled plasmid DNAs of 2000-10,000 bp. Mechanisms of molecular sieving in polymer liquid crystals must differ in fundamental ways from separations in random polymer gels because molecules move around uncrosslinked obstacles that are larger than the smallest dimensions of typical analytes. Molecular sieving in Pluronic liquid crystals is envisioned to occur as molecules squeeze between hydrated (EO)x strands of micelle brushes, or through brushtips and interstitial spaces between micelles. Small molecules such as nucleotides appear to separate by a different mechanism involving partitioning between hydrophilic and hydrophobic environments. This process is termed "hydrophobic interaction electrophoresis". The unique structures of Pluronic copolymers and their liquid crystalline phases provide new challenges and opportunities in separations science.