The protein-folding problem: the native fold determines packing, but does packing determine the native fold?

Dry Molten Globule-Like Intermediates in Protein Folding, Function, and Disease

The performance of a protein depends on its correct folding to the final functional native form. Hence, understanding the process of protein folding has remained an important field of research for the scientific community for the past five decades. Two important intermediate states, namely, wet molten globule (WMG) and dry molten globule (DMG), have emerged as critical milestones during protein folding-unfolding reactions. While much has been discussed about WMGs as a common unfolding intermediate, the evidence for DMGs has remained elusive owing to their near-native features, which makes them difficult to probe using global structural probes. This Review puts together the available literature and new evidence on DMGs to give a broader perspective on the universality of DMGs and discuss their significance in protein folding, function, and disease.

A driving force for polypeptide and protein collapse

Polypeptide collapse is driven by the solvent-excluded volume decrease, the presence of nonpolar side chains is not so important.

Evidence for Dry Molten Globule-Like Domains in the pH-Induced Equilibrium Folding Intermediate of a Multidomain Protein

The role of van der Waals (vdW) packing interactions compared to the hydrophobic effect in stabilizing the functional structure of proteins is poorly understood. Here we show, using fluorescence resonance energy transfer, dynamic fluorescence quenching, red-edge excitation shift, and near- and far-UV circular dichroism, that the pH-induced structural perturbation of a multidomain protein leads to the formation of a state in which two out of the three domains have characteristics of dry molten globules, that is, the domains are expanded compared to the native protein with disrupted packing interactions but have dry cores. We quantitatively estimate the energetic contribution of vdW interactions and show that they play an important role in the stability of the native state and cooperativity of its structural transition, in addition to the hydrophobic effect. Our results also indicate that during the pH-induced unfolding, side-chain unlocking and hydrophobic solvation occur in two distinct steps and not in a concerted manner, as commonly believed.

The Unfolding MD Simulations of Cyclophilin: Analyzed by Surface Contact Networks and Their Associated Metrics

Currently, considerable interest exists with regard to the dissociation of close packed aminoacids within proteins, in the course of unfolding, which could result in either wet or dry moltenglobules. The progressive disjuncture of residues constituting the hydrophobic core ofcyclophilin from L. donovani (LdCyp) has been studied during the thermal unfolding of the molecule, by molecular dynamics simulations. LdCyp has been represented as a surface contactnetwork (SCN) based on the surface complementarity (Sm) of interacting residues within themolecular interior. The application of Sm to side chain packing within proteins make it a very sensitive indicator of subtle perturbations in packing, in the thermal unfolding of the protein. Network based metrics have been defined to track the sequential changes in the disintegration ofthe SCN spanning the hydrophobic core of LdCyp and these metrics prove to be highly sensitive compared to traditional metrics in indicating the increased conformational (and dynamical) flexibility in the network. These metrics have been applied to suggest criteria distinguishing DMG, WMG and transition state ensembles and to identify key residues involved in crucial conformational/topological events during the unfolding process.

The food and environmental safety of Bt crops

Bacillus thuringiensis (Bt) microbial pesticides have a 50-year history of safety in agriculture. Cry proteins are among the active insecticidal ingredients in these pesticides, and genes coding for Cry proteins have been introduced into agricultural crops using modern biotechnology. The Cry gene sequences are often modified to enable effective expression in planta and several Cry proteins have been modified to increase biological activity against the target pest(s). Additionally, the domains of different but structurally conserved Cry proteins can be combined to produce chimeric proteins with enhanced insecticidal properties. Environmental studies are performed and include invertebrates, mammals, and avian species. Mammalian studies used to support the food and feed safety assessment are also used to support the wild mammal assessment. In addition to the NTO assessment, the environmental assessment includes a comparative assessment between the Bt crop and the appropriate conventional control that is genetically similar but lacks the introduced trait to address unintended effects. Specific phenotypic, agronomic, and ecological characteristics are measured in the Bt crop and the conventional control to evaluate whether the introduction of the insect resistance has resulted in any changes that might cause ecological harm in terms of altered weed characteristics, susceptibility to pests, or adverse environmental impact. Additionally, environmental interaction data are collected in field experiments for Bt crop to evaluate potential adverse effects. Further to the agronomic and phenotypic evaluation, potential movement of transgenes from a genetically modified crop plants into wild relatives is assessed for a new pest resistance gene in a new crop. This review summarizes the evidence for safety of crops containing Cry proteins for humans, livestock, and other non-target organisms.

The food and environmental safety of Bt crops

Bacillus thuringiensis (Bt) microbial pesticides have a 50-year history of safety in agriculture. Cry proteins are among the active insecticidal ingredients in these pesticides, and genes coding for Cry proteins have been introduced into agricultural crops using modern biotechnology. The Cry gene sequences are often modified to enable effective expression in planta and several Cry proteins have been modified to increase biological activity against the target pest(s). Additionally, the domains of different but structurally conserved Cry proteins can be combined to produce chimeric proteins with enhanced insecticidal properties. Environmental studies are performed and include invertebrates, mammals, and avian species. Mammalian studies used to support the food and feed safety assessment are also used to support the wild mammal assessment. In addition to the NTO assessment, the environmental assessment includes a comparative assessment between the Bt crop and the appropriate conventional control that is genetically similar but lacks the introduced trait to address unintended effects. Specific phenotypic, agronomic, and ecological characteristics are measured in the Bt crop and the conventional control to evaluate whether the introduction of the insect resistance has resulted in any changes that might cause ecological harm in terms of altered weed characteristics, susceptibility to pests, or adverse environmental impact. Additionally, environmental interaction data are collected in field experiments for Bt crop to evaluate potential adverse effects. Further to the agronomic and phenotypic evaluation, potential movement of transgenes from a genetically modified crop plants into wild relatives is assessed for a new pest resistance gene in a new crop. This review summarizes the evidence for safety of crops containing Cry proteins for humans, livestock, and other non-target organisms.

Kinetic evidence for a two-stage mechanism of protein denaturation by guanidinium chloride

Dry molten globular (DMG) intermediates, an expanded form of the native protein with a dry core, have been observed during denaturant-induced unfolding of many proteins. These observations are counterintuitive because traditional models of chemical denaturation rely on changes in solvent-accessible surface area, and there is no notable change in solvent-accessible surface area during the formation of the DMG. Here we show, using multisite fluorescence resonance energy transfer, far-UV CD, and kinetic thiol-labeling experiments, that the guanidinium chloride (GdmCl)-induced unfolding of RNase H also begins with the formation of the DMG. Population of the DMG occurs within the 5-ms dead time of our measurements. We observe that the size and/or population of the DMG is linearly dependent on [GdmCl], although not as strongly as the second and major step of unfolding, which is accompanied by core solvation and global unfolding. This rapid GdmCl-dependent population of the DMG indicates that GdmCl can interact with the protein before disrupting the hydrophobic core. These results imply that the effect of chemical denaturants cannot be interpreted solely as a disruption of the hydrophobic effect and strongly support recent computational studies, which hypothesize that chemical denaturants first interact directly with the protein surface before completely unfolding the protein in the second step (direct interaction mechanism).

Fast closure of long loops at the initiation of the folding transition of globular proteins studied by time-resolved FRET-based methods

The protein folding problem would be considered “solved” when it will be possible to “read genes”, i.e., to predict the native fold of proteins, their dynamics, and the mechanism of fast folding based solely on sequence data. The long-term goal should be the creation of an algorithm that would simulate the stepwise mechanism of folding, which constrains the conformational space and in which random search for stable interactions is possible. Here, we focus attention on the initial phases of the folding transition starting with the compact disordered collapsed ensemble, in search of the initial sub-domain structural biases that direct the otherwise stochastic dynamics of the backbone. Our studies are designed to test the “loop hypothesis”, which suggests that fast closure of long loop structures by non-local interactions between clusters of mainly non-polar residues is an essential conformational step at the initiation of the folding transition of globular proteins. We developed and applied experimental methods based on time-resolved resonance excitation energy transfer (trFRET) measurements combined with fast mixing methods and studied the initial phases of the folding of

Toxicological evaluation of proteins introduced into food crops

This manuscript focuses on the toxicological evaluation of proteins introduced into GM crops to impart desired traits. In many cases, introduced proteins can be shown to have a history of safe use. Where modifications have been made to proteins, experience has shown that it is highly unlikely that modification of amino acid sequences can make a non-toxic protein toxic. Moreover, if the modified protein still retains its biological function, and this function is found in related proteins that have a history of safe use (HOSU) in food, and the exposure level is similar to functionally related proteins, then the modified protein could also be considered to be "as-safe-as" those that have a HOSU. Within nature, there can be considerable evolutionary changes in the amino acid sequence of proteins within the same family, yet these proteins share the same biological function. In general, food crops such as maize, soy, rice, canola etc. are subjected to a variety of processing conditions to generate different food products. Processing conditions such as cooking, modification of pH conditions, and mechanical shearing can often denature proteins in these crops resulting in a loss of functional activity. These same processing conditions can also markedly lower human dietary exposure to (functionally active) proteins. Safety testing of an introduced protein could be indicated if its biological function was not adequately characterized and/or it was shown to be structurally/functionally related to proteins that are known to be toxic to mammals.

Packing in molten globules and native states

Close packing of hydrophobic residues in the protein interior is an important determinant of protein stability. Cavities introduced by large to small substitutions are known to destabilize proteins. Conversely, native states of proteins and protein fragments can be stabilized by filling in existing cavities. Molten globules (MGs) were initially used to describe a state of protein which has well-defined secondary structure but little or no tertiary packing. Subsequent studies have shown that MGs do have some degree of native-like topology and specific packing. Wet molten globules (WMGs) with hydrated cores and considerably decreased packing relative to the native state have been studied extensively. Recently there has been renewed interest in identification and characterization of dry molten globules (DMGs). These are slightly expanded forms of the native state which show increased conformational flexibility, native-like main-chain hydrogen bonding and dry interiors. The generality of occurrence of DMGs during protein unfolding and the extent and nature of packing in DMGs remain to be elucidated. Packing interactions in native proteins and MGs can be probed through mutations. Next generation sequencing technologies make it possible to determine relative populations of mutants in a large pool. When this is coupled to phenotypic screens or cell-surface display, it becomes possible to rapidly examine large panels of single-site or multi-site mutants. From such studies, residue specific contributions to protein stability and function can be estimated in a highly parallelized fashion. This complements conventional biophysical methods for characterization of packing in native states and molten globules.

Dry molten globule intermediates and the mechanism of protein unfolding

New experimental results show that either gain or loss of close packing can be observed as a discrete step in protein folding or unfolding reactions. This finding poses a significant challenge to the conventional two-state model of protein folding. Results of interest involve dry molten globule (DMG) intermediates, an expanded form of the protein that lacks appreciable solvent. When an unfolding protein expands to the DMG state, side chains unlock and gain conformational entropy, while liquid-like van der Waals interactions persist. Four unrelated proteins are now known to form DMGs as the first step of unfolding, suggesting that such an intermediate may well be commonplace in both folding and unfolding. Data from the literature show that peptide amide protons are protected in the DMG, indicating that backbone structure is intact despite loss of side-chain close packing. Other complementary evidence shows that secondary structure formation provides a major source of compaction during folding. In our model, the major free-energy barrier separating unfolded from native states usually occurs during the transition between the unfolded state and the DMG. The absence of close packing at this barrier provides an explanation for why phi-values, derived from a Brønsted-Leffler plot, depend primarily on structure at the mutational site and not on specific side-chain interactions. The conventional two-state folding model breaks down when there are DMG intermediates, a realization that has major implications for future experimental work on the mechanism of protein folding.

Minimal energy packings and collapse of sticky tangent hard-sphere polymers

We enumerate all minimal energy packings (MEPs) for small single linear and ring polymers composed of spherical monomers with contact attractions and hard-core repulsions and compare them to corresponding results for monomer packings. We define and identify "dividing surfaces" in polymer packings, which reduce the number of arrangements that satisfy hard-sphere and covalent-bond constraints. Compared to monomer MEPs, polymer MEPs favor intermediate structural symmetry. We also examine the packing-preparation dependence for longer single chains using molecular dynamics simulations. For slow temperature quenches, chains form crystallites with close-packed cores. As the quench rate increases, the core size decreases and the exterior becomes more disordered. By examining the contact number, we connect the suppression of crystallization to the onset of isostaticity in disordered packings.

Revisiting the myths of protein interior: studying proteins with mass-fractal hydrophobicity-fractal and polarizability-fractal dimensions

A robust marker to describe mass, hydrophobicity and polarizability distribution holds the key to deciphering structural and folding constraints within proteins. Since each of these distributions is inhomogeneous in nature, the construct should be sensitive in describing the patterns therein. We show, for the first time, that the hydrophobicity and polarizability distributions in protein interior follow fractal scaling. It is found that (barring ‘all-α’) all the major structural classes of proteins have an amount of unused hydrophobicity left in them. This amount of untapped hydrophobicity is observed to be greater in thermophilic proteins, than that in their (structurally aligned) mesophilic counterparts. ‘All-β’(thermophilic, mesophilic alike) proteins are found to have maximum amount of unused hydrophobicity, while ‘all-α’ proteins have been found to have minimum polarizability. A non-trivial dependency is observed between dielectric constant and hydrophobicity distributions within (α+β) and ‘all-α’ proteins, whereas absolutely no dependency is found between them in the ‘all-β’ class. This study proves that proteins are not as optimally packed as they are supposed to be. It is also proved that origin of α-helices are possibly not hydrophobic but electrostatic; whereas β-sheets are predominantly hydrophobic in nature. Significance of this study lies in protein engineering studies; because it quantifies the extent of packing that ensures protein functionality. It shows that myths regarding protein interior organization might obfuscate our knowledge of actual reality. However, if the later is studied with a robust marker of strong mathematical basis, unknown correlations can still be unearthed; which help us to understand the nature of hydrophobicity, causality behind protein folding, and the importance of anisotropic electrostatics in stabilizing a highly complex structure named ‘proteins’.

Satisfaction of hydrogen-bonding potential influences the conservation of polar sidechains

Although polar amino acids tend to be found on the surface of proteins due to their hydrophilic nature, their important roles within the core of proteins are now becoming better recognized. It has long been understood that a significant number of mainchain functions will not achieve hydrogen bond satisfaction through the formation of secondary structures; in these circumstances, it is generally buried polar residues that provide hydrogen bond satisfaction. Here, we describe an analysis of the hydrogen-bonding of polar amino acids in a set of structurally aligned protein families. This allows us not only to calculate the conservation of each polar residue but also to assess whether conservation is correlated with the hydrogen-bonding potential of polar sidechains. We show that those polar sidechains whose hydrogen-bonding potential is satisfied tend to be more conserved than their unsatisfied or nonhydrogen-bonded counterparts, particularly when buried. Interestingly, these buried and satisfied polar residues are significantly more conserved than buried hydrophobic residues. Forming hydrogen bonds to mainchain amide atoms also influences conservation, with those satisfied buried polar residues that form two hydrogen bonds to mainchain amides being significantly more conserved than those that form only one or none. These results indicate that buried polar residues whose hydrogen-bonding potential is satisfied are likely to have important roles in maintaining protein structure.

Orientational distributions of contact clusters in proteins closely resemble those of an icosahedron

The orientational geometry of residue packing in proteins was studied in the past by superimposing clusters of neighboring residues with several simple lattices (Bagci et al., Proteins 2003;53:56-67; Raghunathan et al., Protein Sci 1997;6:2072-2083). In this work, instead of a lattice we use the regular polyhedron, the icosahedron, as the model to describe the orientational distribution of contacts in clusters derived from a high-resolution protein dataset (522 protein structures with high resolution < 1.5 A). We find that the order parameter (orientation function) measuring the angular overlap of directions in coordination clusters with directions of the icosahedron is 0.91, which is a significant improvement in comparison with the value 0.82 for the order parameter with the face-centered cubic (fcc) lattice. Close packing tendencies and patterns of residue packing in proteins are considered in detail and a theoretical description of these packing regularities is proposed.

A heterospecific leucine zipper tetramer

Protein-protein interactions dictate the assembly of the macromolecular complexes essential for functional networks and cellular behavior. Elucidating principles of molecular recognition governing important interfaces such as coiled coils is a challenging goal for structural and systems biology. We report here that two valine-containing mutants of the GCN4 leucine zipper that fold individually as four-stranded coiled coils associate preferentially in mixtures to form an antiparallel, heterotetrameric structure. X-ray crystallographic analysis reveals that the coinciding hydrophobic interfaces of the hetero- and homotetramers differ in detail, explaining their partnering and structural specificity. Equilibrium disulfide exchange and thermal denaturation experiments show that the 50-fold preference for heterospecificity results from a combination of preferential packing and hydrophobicity. The extent of preference is sensitive to the side chains comprising the interface. Thus, heterotypic versus homotypic interaction specificity in coiled coils reflects a delicate balance in complementarity of shape and chemistry of the participating side chains.

Thermodynamic properties of BPTI variants with highly simplified amino acid sequences

We report the first detailed thermodynamic analysis of simplified proteins by differential scanning calorimetry (DSC). The experiments were carried out with five simplified BPTI variants, whose structures and activities have been reported, in which several residues not essential for specifying the tertiary structure were replaced by alanine. In most aspects, the thermodynamics of simplified proteins were very similar to, if not essentially identical with, those of natural proteins. In particular, they undergo a highly cooperative two-state thermal unfolding process with a large enthalpy change, which is a thermodynamic hallmark of the native state of natural globular proteins. Furthermore, the specific enthalpy and entropy changes upon unfolding at 110 degrees C were close to values invariably observed for small natural globular proteins (55 J g(-1) and ~16 J K(-1) g(-1), respectively). On the other hand, two simplified BPTI variants, BPTI-21 and BPTI-22 (containing 21 and 22 alanine residues), were enthalpically stabilized while entropically destabilized with respect to the reference BPTI-[5,55] molecule. This peculiar type of entropy-enthalpy compensation is in sharp contrast to the usual enthalpy destabilization/entropy stabilization observed in mutational studies of natural proteins. Overall, we conclude that a thermodynamic native state can be achieved by proteins encoded with extensively simplified sequences.

Geometry of nonbonded interactions involving planar groups in proteins

Although hydrophobic interaction is the main contributing factor to the stability of the protein fold, the specificity of the folding process depends on many directional interactions. An analysis has been carried out on the geometry of interaction between planar moieties of ten side chains (Phe, Tyr, Trp, His, Arg, Pro, Asp, Glu, Asn and Gln), the aromatic residues and the sulfide planes (of Met and cystine), and the aromatic residues and the peptide planes within the protein tertiary structures available in the Protein Data Bank. The occurrence of hydrogen bonds and other nonconventional interactions such as C-H...pi, C-H...O, electrophile-nucleophile interactions involving the planar moieties has been elucidated. The specific nature of the interactions constraints many of the residue pairs to occur with a fixed sequence difference, maintaining a sequential order, when located in secondary structural elements, such as alpha-helices and beta-turns. The importance of many of these interactions (for example, aromatic residues interacting with Pro or cystine sulfur atom) is revealed by the higher degree of conservation observed for them in protein structures and binding regions. The planar residues are well represented in the active sites, and the geometry of their interactions does not deviate from the general distribution. The geometrical relationship between interacting residues provides valuable insights into the process of protein folding and would be useful for the design of protein molecules and modulation of their binding properties.

All-atom ab initio folding of a diverse set of proteins

Natural proteins fold to a unique, thermodynamically dominant state. Modeling of the folding process and prediction of the native fold of proteins are two major unsolved problems in biophysics. Here, we show successful all-atom ab initio folding of a representative diverse set of proteins by using a minimalist transferable-energy model that consists of two-body atom-atom interactions, hydrogen bonding, and a local sequence-energy term that models sequence-specific chain stiffness. Starting from a random coil, the native-like structure was observed during replica exchange Monte Carlo (REMC) simulation for most proteins regardless of their structural classes; the lowest energy structure was close to native-in the range of 2-6 A root-mean-square deviation (rmsd). Our results demonstrate that the successful folding of a protein chain to its native state is governed by only a few crucial energetic terms.

The jigsaw puzzle model: search for conformational specificity in protein interiors

The jigsaw puzzle model postulates that the predominant factor relating primary sequence to three-dimensional fold lies in the stereospecific packing of interdigitating side-chains within densely packed protein interiors. An attempt has been made to check the validity of the model by means of a surface complementarity function. Out of a database of 100 highly resolved protein structures the contacts between buried hydrophobic residues (Leu, Ile, Val, Phe) and their neighbours have been categorized in terms of the extent of side-chain surface area involved in a contact (overlap) and their steric fit (Sm). The results show that the majority of contacts between a buried residue and its immediate neighbours (side-chains) are of high steric fit and in the case of extended overlap at least one of the angular parameters characterizing interresidue geometry to have pronounced deviation from a random distribution, estimated by chi(2). The calculations thus tend to support the "jigsaw puzzle" model in that 75-85% of the contacts involving hydrophobic residues are of high surface complementarity, which, coupled to high overlap, exercise fairly stringent constraints over the possible geometrical orientations between interacting residues. These constraints manifest in simple patterns in the distributions of orientational angles. Approximately 60-80% of the buried side-chain surface packs against neighbouring side-chains, the rest interacting with main-chain atoms. The latter partition of the surface maintains an equally high steric fit (relative to side-chain contacts) emphasizing a non-trivial though secondary role played by main-chain atoms in interior packing. The majority of this class of contacts, though of high complementarity, is of reduced overlap. All residues whether hydrophobic or polar/charged show similar surface complementarity measures upon burial, indicating comparable competence of all amino acids in packing effectively with their atomic environments. The specificity thus appears to be distributed over the entire network of contacts within proteins. The study concludes with a proposal to classify contacts as specific and non-specific (based on overlap and fit), with the former perhaps contributing more to the specificity between sequence and fold than the latter.

The origin and extent of coarse-grained regularities in protein internal packing

Despite the suitability of various lattice geometries for coarse-grained modeling of proteins, the actual packing geometry of residues in folded structures has remained largely unexplored. A strong tendency to assume a regular packing geometry is shown here by optimally reorienting and superimposing clusters of neighboring residues from databank structures examined on a coarse-grained (single-site-per-residue) scale. The orientation function (or order parameter) of the examined coordination clusters with respect to fcc lattice directions is found to be 0.82. The observed geometry, which may be termed an incomplete distorted face-centered cubic (fcc) packing, is apparently favored by the drive to maximize packing density, in a fashion analogous to the way identical spheres pack densely and follow fcc geometry. About 2/3 of all residues obey this packing geometry, while the remainder occupy other context-dependent positions. The preferred coordination directions show relatively small variations over the various amino acid types, consistent with uniform residue viewpoint. Both the extremes of solvent-exposed and completely buried residue neighborhoods approximate the same generic packing, the only difference being in the numbers (and not the orientations) of coordination sites that are occupied (or left void for solvent occupancy). We observe the prevalence of a rather uniform (tight) residue packing density throughout the structure, including even the residues packed near solvent-exposed regions. The observed orientation distribution reveals an underlying, intrinsic orientation lattice for proteins.

Application of a novel graph-theoretic folding degree index to the study of steroid-DB3 antibody binding affinity

A novel folding degree index, together with other macromolecular descriptors, is used to study steroid-DB3 antibody interactions. This index is based on graph spectral moments of a matrix representing the dihedral angles of a protein backbone. The causes influencing the different order of binding affinity of steroids to DB3 antibody are identified. It is shown that the changes in the chain compactness of the DB3 antibody with respect to its center of mass (radius of gyration) is compensated by a change in the folding degree index in the contrary sense. In fact, the increment in compactness of chain L and the lower increment in the folding degree index of chain H are able to explain the variations in affinity for DB3 of the steroids studied. Consequently, the highest binding affinities are reached by increasing the compactness of chain L in DB3 at the same time that producing the smallest increment in the folding degree of chain H. This study shows the possibilities of application for the graph-theoretic folding degree index in studying drug-protein interactions.

Stably folded de novo proteins from a designed combinatorial library

Binary patterning of polar and nonpolar amino acids has been used as the key design feature for constructing large combinatorial libraries of de novo proteins. Each position in a binary patterned sequence is designed explicitly to be either polar or nonpolar; however, the precise identities of these amino acids are varied extensively. The combinatorial underpinnings of the "binary code" strategy preclude explicit design of particular side chains at specified positions. Therefore, packing interactions cannot be specified a priori. To assess whether the binary code strategy can nonetheless produce well-folded de novo proteins, we constructed a second-generation library based upon a new structural scaffold designed to fold into 102-residue four-helix bundles. Characterization of five proteins chosen arbitrarily from this new library revealed that (1) all are alpha-helical and quite stable; (2) four of the five contain an abundance of tertiary interactions indicative of well-ordered structures; and (3) one protein forms a well-folded structure with native-like features. The proteins from this new 102-residue library are substantially more stable and dramatically more native-like than those from an earlier binary patterned library of 74-residue sequences. These findings demonstrate that chain length is a crucial determinant of structural order in libraries of de novo four-helix bundles. Moreover, these results show that the binary code strategy--if applied to an appropriately designed structural scaffold--can generate large collections of stably folded and/or native-like proteins.

A procedure for detection and quantitation of cavity volumes proteins. Application to measure the strength of the hydrophobic driving force in protein folding

Accurate identification of cavities is important in the study of protein structure, stability, design, and ligand binding. Identification and quantitation of cavities is a nontrivial problem because most cavities are connected to the protein exterior. We describe a computational procedure for quantitating cavity volumes and apply this to derive an estimate of the hydrophobic driving force in protein folding. A grid-based Monte Carlo procedure is used to position water molecules on the surface of a protein. A Voronoi procedure is used to identify and quantitate empty space within the solvated protein. Additional cavities not detected by other existing procedures can be identified. Most of these are close to surface concavities. Residue volumes for both the interior and the surface residues as well as cavity volumes are in good agreement with volumes calculated from fully hydrated protein structures obtained from molecular dynamic simulations. We show that the loss of stability because of cavity-creating mutations correlates better with cavity volumes determined by this procedure than with cavity volumes determined by other methods. Available structural and thermodynamic data for a number of cavity-containing mutants were analyzed to obtain estimates of 26.1 cal x mol(-1) x A(-3) and 18.5 cal x mol(-1) x A(-2) for the relative contributions of cavity formation and the hydrophobic effect to the observed stability changes. The present estimate for the hydrophobic driving force is at the lower end of estimates derived from model compound studies and considerably lower than previous estimates of approximately 50 cal x mol(-1) x A(-2) derived from protein mutational data. In the absence of structural rearrangement, on average, deletion of a single methylene group is expected to result in losses in stability of 0.41 and 0.70 kcal x mol(-1) resulting from decrease in hydrophobicity and packing, respectively.

Four-body potentials reveal protein-specific correlations to stability changes caused by hydrophobic core mutations

Mutational experiments show how changes in the hydrophobic cores of proteins affect their stabilities. Here, we estimate these effects computationally, using four-body likelihood potentials obtained by simplicial neighborhood analysis of protein packing (SNAPP). In this procedure, the volume of a known protein structure is tiled with tetrahedra having the center of mass of one amino acid side-chain at each vertex. Log-likelihoods are computed for the 8855 possible tetrahedra with equivalent compositions from structural databases and amino acid frequencies. The sum of these four-body potentials for tetrahedra present in a given protein yields the SNAPP score. Mutations change this sum by changing the compositions of tetrahedra containing the mutated residue and their related potentials. Linear correlation coefficients between experimental mutational stability changes, Delta(DeltaG(unfold)), and those based on SNAPP scoring range from 0.70 to 0.94 for hydrophobic core mutations in five different proteins. Accurate predictions for the effects of hydrophobic core mutations can therefore be obtained by virtual mutagenesis, based on changes to the total SNAPP likelihood potential. Significantly, slopes of the relation between Delta(DeltaG(unfold)) and DeltaSNAPP for different proteins are statistically distinct, and we show that these protein-specific effects can be estimated using the average SNAPP score per residue, which is readily derived from the analysis itself. This result enhances the predictive value of statistical potentials and supports previous suggestions that "comparable" mutations in different proteins may lead to different Delta(DeltaG(unfold)) values because of differences in their flexibility and/or conformational entropy.

The slow folding reaction of barstar: the core tryptophan region attains tight packing before substantial secondary and tertiary structure formation and final compaction of the polypeptide chain

The slow folding of a single tryptophan-containing mutant of barstar has been studied in the presence of 2 M urea at 10 degrees C, using steady state and time-resolved fluorescence methods and far and near-UV CD measurements. The protein folds in two major phases: a fast phase, which is lost in the dead time of measurement during which the polypeptide collapses to a compact form, is followed by a slow observable phase. During the fast phase, the rotational correlation time of Trp53 increases from 2.2 ns to 7.2 ns, and its mean fluorescence lifetime increases from 2.3 ns to 3.4 ns. The fractional changes in steady-state fluorescence, far-UV CD, and near-UV CD signals, which are associated with the fast phase are, respectively, 36 %, 46 %, and 16 %. The product of the fast phase can bind the hydrophobic dye ANS. These observations together suggest that the folding intermediate accumulated at the end of the fast phase has: (a) about 20 % of the native-state secondary structure, (b) marginally formed or disordered tertiary structure, (c) a water-intruded and mobile protein interior; and (d) solvent-accessible patches of hydrophobic groups. Measurements of the anisotropy decay of Trp53 suggest that it undergoes two types of rotational motion in the intermediate: (i) fast (tau(r) approximately 1 ns) local motion of its indole side-chain, and (ii) a slower (tau(r) approximately 7.2 ns) motion corresponding to global tumbling of the entire protein molecule. The ability of the Trp53 side-chain to undergo fast local motion in the intermediate, but not in the fully folded protein where it is completely buried in the hydrophobic core, suggests that the core of the intermediate is still poorly packed. The global tumbling time of the fully folded protein is faster at 5.6 ns, suggesting that the volume of the intermediate is 25 % more than that of the fully folded protein. The rate of folding of this intermediate to the native state, measured by steady-state fluorescence, far-UV CD, and near-UV CD, is 0.07(+/-0.01) min(-1) This rate compares to a rate of folding of 0.03(+/-0.005) min(-1), determined by double-jump experiments which monitor directly formation of native protein; and to a rate of folding of 0.05 min(-1), when determined from time-resolved anisotropy measurements of the long rotational correlation time, which relaxes from an initial value of 7.2 ns to a final value of 5. 6 ns as the protein folds. On the other hand, the amplitude of the short correlation time decreases rapidly with a rate of 0.24(+/-0.06) min(-1). These results suggest that tight packing of residues in the hydrophobic core occurs relatively early during the observable slow folding reaction, before substantial secondary and tertiary structure formation and before final compaction of the protein.

Effects of side-chain characteristics on stability and oligomerization state of a de novo-designed model coiled-coil: 20 amino acid substitutions in position "d"

We describe the de novo design and biophysical characterization of a model coiled-coil protein in which we have systematically substituted 20 different amino acid residues in the central "d" position. The model protein consists of two identical 38 residue polypeptide chains covalently linked at their N termini via a disulfide bridge. The hydrophobic core contained Val and Ile residues at positions "a" and Leu residues at positions "d". This core allowed for the formation of both two-stranded and three-stranded coiled-coils in benign buffer, depending on the substitution at position "d". The structure of each analog was analyzed by CD spectroscopy and their relative stability determined by chemical denaturation using GdnHCI (all analogs denatured from the two-stranded state). The oligomeric state(s) was determined by high-performance size-exclusion chromatography and sedimentation equilibrium analysis in benign medium. Our results showed a thermodynamic stability order (in order of decreasing stability) of: Leu, Met, Ile, Tyr, Phe, Val, Gln, Ala, Trp, Asn, His, Thr, Lys, Ser, Asp, Glu, Arg, Orn, and Gly. The Pro analog prevented coiled-coil formation. The overall stability range was 7.4 kcal/mol from the lowest to the highest analog, indicating the importance of the hydrophobic core and the dramatic effect a single substitution in the core can have upon the stability of the protein fold. In general, the side-chain contribution to the level of stability correlated with side-chain hydrophobicity. Molecular modelling studies, however, showed that packing effects could explain deviations from a direct correlation. In regards to oligomerization state, eight analogs demonstrated the ability to populate exclusively one oligomerization state in benign buffer (0.1 M KCl, 0.05 M K(2)PO(4)(pH 7)). Ile and Val (the beta-branched residues) induced the three-stranded oligomerization state, whereas Tyr, Lys, Arg, Orn, Glu and Asp induced the two-stranded state. Asn, Gln, Ser, Ala, Gly, Phe, Leu, Met and Trp analogs were indiscriminate and populated two-stranded and three-stranded states. Comparison of these results with similar substitutions in position "a" highlights the positional effects of individual residues in defining the stability and numbers of polypeptide chains occurring in a coiled-coil structure. Overall, these results in conjunction with other work now generate a relative thermodynamic stability scale for 19 naturally occurring amino acid residues in either an "a" or "d" position of a two-stranded coiled-coil. Thus, these results will aid in the de novo design of new coiled-coil structures, a better understanding of their structure/function relationships and the design of algorithms to predict the presence of coiled-coils within native protein sequences.

Structural cassette mutagenesis in a de novo designed protein: proof of a novel concept for examining protein folding and stability

The solution to the protein folding problem lies in defining the relative energetic contributions of short-range and long-range interactions. In other words, the tendency of a stretch of amino acids to adopt a final secondary structural fold is context dependent. Our approach to this problem is to address whether an amino acid sequence, a "cassette," with a defined secondary structure in the three-dimensional structure of a native protein, can adopt a different conformation when placed into a different protein environment. Thus, we designed de novo a disulfide-bridged two-stranded alpha-helical parallel coiled coil, where each polypeptide chain consisted of 39 residues, as a "cassette holder." The 11-residue cassette would be inserted into the center of each polypeptide chain between the two nucleating alpha-helices to replace the control sequence. This Structural Cassette Mutagenesis model permits the analysis of short-range interactions within the inserted cassette as well as long-range interactions between the nucleating helices and the cassette region. The cassette holder, with a control sequence as the cassette, had a GdnHCl transition midpoint during denaturation of 5.6M. To demonstrate the feasibility of our model, an 11-residue beta-strand cassette from an immunoglobulin fold was inserted. The cassette was fully induced into the alpha-helical conformation with a [GdnHCl]1/2 value of 3.2M. To demonstrate the importance of short-range interactions (beta-sheet/alpha-helical propensities of amino acid side chains) in modulating structure and stability, a series of 1-5 threonine residues (highest beta-sheet propensity) were substituted into the solvent-exposed portions of the cassette in the alpha-helical conformation. Each successive substitution systematically decreased the stability of the coiled coil with peptide T4b (4 Thr residues) having a [GdnHCl]1/2 value of 2.2M. The single substitution of Ile in the hydrophobic core of the cassette with Ala or Thr had the most dramatic effect on protein stability (peptide 120T, [GdnHCl]1/2 value of 1.4M). Though these substitutions were able to modulate stability, they were not able to disrupt the alpha-helical conformation of the cassette, showing the importance of the nucleating alpha-helices on either side of the cassette in controlling conformation of the cassette. We have demonstrated the feasibility of our model protein to accept a beta-strand cassette. The effect of cassettes containing other beta-strands, beta-turns, loops, regions of undefined structure, and helical segments on conformation and stability of our model protein will also be determined.

Optimal region of average side-chain entropy for fast protein folding

Search and study the general principles that govern kinetics and thermodynamics of protein folding generates new insight into the factors that control this process. Here, we demonstrate based on the known experimental data and using theoretical modeling of protein folding that side-chain entropy is one of the general determinants of protein folding. We show for proteins belonging to the same structural family that there exists an optimal relationship between the average side-chain entropy and the average number of contacts per residue for fast folding kinetics. Analysis of side-chain entropy for proteins that fold without additional agents demonstrates that there exists an optimal region of average side-chain entropy for fast folding. Deviation of the average side-chain entropy from the optimal region results in an anomalous protein folding process (prions, alpha-lytic protease, subtilisin, some DNA-binding proteins). Proteins with high or low side-chain entropy would have extended unfolded regions and would require some additional agents for complete folding. Such proteins are common in nature, and their structure properties have biological importance.

De novo protein design. I. In search of stability and specificity

We have developed a fully automated protein design strategy that works on the entire sequence of the protein and uses a full atom representation. At each step of the procedure, an all-atom model of the protein is built using the template protein structure and the current designed sequence. The energy of the model is used to drive a Monte Carlo optimization in sequence space: random moves are either accepted or rejected based on the Metropolis criterion. We rely on the physical forces that stabilize native protein structures to choose the optimum sequence. Our energy function includes van der Waals interactions, electrostatics and an environment free energy. Successful protein design should be specific and generate a sequence compatible with the template fold and incompatible with competing folds. We impose specificity by maintaining the amino acid composition constant, based on the random energy model. The specificity of the optimized sequence is tested by fold recognition techniques. Successful sequence designs for the B1 domain of protein G, for the lambda repressor and for sperm whale myoglobin are presented. We show that each additional term of the energy function improves the performance of our design procedure: the van der Waals term ensures correct packing, the electrostatics term increases the specificity for the correct native fold, and the environment solvation term ensures a correct pattern of buried hydrophobic and exposed hydrophilic residues. For the globin family, we show that we can design a protein sequence that is stable in the myoglobin fold, yet incompatible with the very similar hemoglobin fold.

Effects of core-packing on the structure, function, and mechanics of a four-helix-bundle protein ROP

The effects of core-packing on the structure, function and mechanics of the RNA-binding 4-helix-bundle Rop have been studied by molecular dynamics simulations. The structural, dynamical and geometrical properties of the Rop homodimer, (formed by the antiparallel juxtaposition of two helix-turn-helix motifs), have been compared with those of three protein variants described by Munson et al. (Protein Sci, 5:1584-1593, 1996), where the core of the native protein has been systematically repacked using a two-amino acid alphabet: Ala(2)Leu(2)-8, Ala(2)Leu(2)-8-rev, and Leu(2)Ala(2)-8. The results showed that it was possible to readily distinguish the inactive protein Leu(2)Ala(2)-8 from the other functionally active systems based on tertiary and quaternary structure criteria. Structural properties such as native secondary structure content did not correlate with biological activity. Biological activity was related in part to the relative arrangement of the residues within the binding site. But, more global aspects, related to the overall topology of the helical bundle, accounted for the small functional differences between Ala(2)Leu(2)-8 and Ala(2)Leu(2)-8-rev. Mechanically, the 4-helix-bundle absorbed core mutations by altering the local structure at the sequence termini and in the turns that join the two helices of each monomer, and by changing the overall orientation and separation of the extremely rigid helices. Proteins 1999;36:436-446.

Solution structure and dynamics of a de novo designed three-helix bundle protein

Although de novo protein design is an important endeavor with implications for understanding protein folding, until now, structures have been determined for only a few 25- to 30-residue designed miniproteins. Here, the NMR solution structure of a complex 73-residue three-helix bundle protein, alpha3D, is reported. The structure of alpha3D was not based on any natural protein, and yet it shows thermodynamic and spectroscopic properties typical of native proteins. A variety of features contribute to its unique structure, including electrostatics, the packing of a diverse set of hydrophobic side chains, and a loop that incorporates common capping motifs. Thus, it is now possible to design a complex protein with a well defined and predictable three-dimensional structure.

Visualizing and quantifying molecular goodness-of-fit: small-probe contact dots with explicit hydrogen atoms

The technique of small-probe contact dot surfaces is described as a method for calculating and displaying the detailed atomic contacts inside or between molecules. It allows one both to measure and to visualize directly the goodness-of-fit of packing interactions. It requires both highly accurate structures and also the explicit inclusion of all hydrogen atoms and their van der Waals interactions. A reference dataset of 100 protein structures was chosen on the basis of resolution (1.7 A or better), crystallographic R-value, non-homology, and the absence of any unusual problems. Hydrogen atoms were added in standard geometry and, where needed, with rotational optimization of OH, SH, and NH+3 positions. Side-chain amide orientations were corrected where required by NH van der Waals clashes, as described in the accompanying paper. It was determined that, in general, methyl groups pack well in the default staggered conformation, except for the terminal methyl groups of methionine residues, which required rotational optimization. The distribution of serious clashes (i.e. non-H-bond overlap of >/=0.4 A) was studied as a function of resolution, alternate conformations, and temperature factor (B), leading to the decision that packing and other structural features would not be analyzed for residues in 'b' alternate conformations or with B-factors of 40 or above. At the level of the fine details analyzed here, structural accuracy improves quite significantly over the range from 1.7 to 1.0 A resolution. These high-resolution structures show impressively well-fitted packing interactions, with some regions thoroughly interdigitated and other regions somewhat sparser. Lower-resolution structures or model structures could undoubtedly be improved in accuracy by the incorporation of this additional information: for example, nucleic acid structures in non-canonical conformations are often very accurate for the bases and much less reliable for the backbone, whose conformation could be specified better by including explicit H atom geometry and contacts. The contact dots are an extremely sensitive method of finding problem areas, and often they can suggest how to make improvements. They can also provide explanations for structural features that have been described only as empirical regularities, which is illustrated by showing that the commonest rotamer of methionine (a left-handed spiral, with all chi values near -60 degrees) is preferred because it provides up to five good H atom van der Waals contacts. This methodology is thus applicable in two different ways: (1) for finding and correcting errors in structure models (either experimental or theoretical); and (2) for analyzing interaction patterns in the molecules themselves.

Packing of sidechains in low-resolution models for proteins

BACKGROUND

Atomic level rotamer libraries for sidechains in proteins have been proposed by several groups. Conformations of side groups in coarse-grained models, on the other hand, have not yet been analyzed, although low resolution approaches are the only efficient way to explore global structural features.

RESULTS

A residue-specific backbone-dependent library for sidechain isomers, compatible with a coarse-grained model, is proposed. The isomeric states are utilized in packing sidechains of known backbone structures. Sidechain positions are predicted with a root-mean-square deviation (r.m.s.d.) of 2.40 A with respect to crystal structure for 50 test proteins. The rmsd for core residues is 1.60 A and decreases to 1.35 A when conformational correlations and directional effects in inter-residue couplings are considered.

CONCLUSIONS

An automated method for assigning sidechain positions in coarse-grained model proteins is proposed and made available on the internet; the method accounts satisfactorily for sidechain packing, particularly in the core.

Hydrophobic core packing and protein design

Over the past few years, we have witnessed exciting advances in protein design. Several groups have reported success in the design of hydrophobic cores, and the principles developed in these studies have been recently applied to the full sequence design of a small protein motif and the design of a catalytically active metal center. These successes suggest that designing large, functional proteins in computero is more feasible than ever before.

A test of the relationship between sequence and structure in proteins: excision of the heme binding site in apocytochrome b5

The water-soluble domain of rat hepatic holocytochrome b5 is an alphabeta protein containing elements of secondary structure in the sequence beta1-alpha1-beta4-beta3-alpha2-alpha3-beta5- alpha4-alpha5-beta2-alpha6. The heme group is enclosed by four helices, a2, a3, a4, and a5. To test the hypothesis that a small b hemoprotein can be constructed in two parts, one forming the heme site, the other an organizing scaffold, a protein fragment corresponding to beta1-alpha1-beta4-beta3-lambda-beta2-alpha6 was prepared, where lambda is a seven-residue linker bypassing the heme binding site. The fragment ("abridged b5") was found to contain alpha and beta secondary structure by circular dichroism spectroscopy and tertiary structure by Trp fluorescence emission spectroscopy. NMR data revealed a species with spectral properties similar to those of the full-length apoprotein. This folded form is in slow equilibrium on the chemical shift time scale with other less folded species. Thermal denaturation, as monitored by circular dichroism, absorption, and fluorescence spectroscopy, as well as size-exclusion chromatography-fast protein liquid chromatography (SEC-FPLC), confirmed the coexistence of at least two distinct conformational ensembles. It was concluded that the protein fragment is capable of adopting a specific fold likely related to that of cytochrome b5, but does not achieve high thermodynamic stability and cooperativity. Abridged b5 demonstrates that the spliced sequence contains the information necessary to fold the protein. It suggests that the dominating influence to restrict the conformational space searched by the chain is structural propensities at a local level rather than internal packing. The sequence also holds the properties necessary to generate a barrier to unfolding.

Simplified proteins: minimalist solutions to the 'protein folding problem'

Recent research has suggested that stable, native proteins may be encoded by simple sequences of fewer than the full set of 20 proteogenic amino acids. Studies of the ability of simple amino acid sequences to encode stable, topologically complex, native conformations and to fold to these conformations in a biologically relevant time frame have provided insights into the sequence determinants of protein structure and folding kinetics. They may also have important implications for protein design and for theories of the origins of protein synthesis itself.

A folding pathway for betapep-4 peptide 33mer: from unfolded monomers and beta-sheet sandwich dimers to well-structured tetramers

It was recently reported that a de novo designed peptide 33mer, betapep-4, can form well-structured beta-sheet sandwich tetramers (Ilyina E, Roongta V, Mayo KH, 1997b, Biochemistry 36:5245-5250). For insight into the pathway of betapep-4 folding, the present study investigates the concentration dependence of betapep-4 self-association by using 1H-NMR pulsed-field gradient (PFG)-NMR diffusion measurements, and circular dichroism. Downfield chemically shifted alphaH resonances, found to arise only from the well-structured betapep-4 tetramer state, yield the fraction of tetramer within the oligomer equilibrium distribution. PFG-NMR-derived diffusion coefficients, D, provide a means for deriving the contribution of monomer and other oligomer states to this distribution. These data indicate that tetramer is the highest oligomer state formed, and that inclusion of monomer and dimer states in the oligomer distribution is sufficient to explain the concentration dependence of D values for betapep-4. Equilibrium constants calculated from these distributions [2.5 x 10(5) M(-1) for M-D and 1.2 x 10(4) M(-1) for D-T at 313 K] decrease only slightly, if at all, with decreasing temperature indicating a hydrophobically mediated, entropy-driven association/folding process. Conformational analyses using NMR and CD provide a picture where "random coil" monomers associate to form molten globule-like beta-sheet sandwich dimers that further associate and fold as well-structured tetramers. Betapep-4 folding is thermodynamically linked to self-association. As with folding of single-chain polypeptides, the final folding step to well-structured tetramer betapep-4 is rate limiting.

Hydrophobic sequence minimization of the alpha-lactalbumin molten globule

The molten globule, a widespread protein-folding intermediate, can attain a native-like backbone topology, even in the apparent absence of rigid side-chain packing. Nonetheless, mutagenesis studies suggest that molten globules are stabilized by some degree of side-chain packing among specific hydrophobic residues. Here we investigate the importance of hydrophobic side-chain diversity in determining the overall fold of the alpha-lactalbumin molten globule. We have replaced all of the hydrophobic amino acids in the sequence of the helical domain with a representative amino acid, leucine. Remarkably, the minimized molecule forms a molten globule that retains many structural features characteristic of a native alpha-lactalbumin fold. Thus, nonspecific hydrophobic interactions may be sufficient to determine the global fold of a protein.

Non-randomness in side-chain packing: the distribution of interplanar angles

We analyze the distributions of interplanar angles between interacting side chains with well-defined planar regions, to see whether these distributions correspond to random packing or alternatively show orientational preferences. We use a non-homologous set of 79 high-resolution protein chain structures to show that the observed distributions are significantly different from the sinusoidal one expected for random packing. Overall, we see a relative excess of small angles and a paucity of large interplanar angles; the difference between the expected and observed distributions can be described as a shift of 5% of the interplanar angles from large (> or = 60 degrees) to small (< 30 degrees) values. By grouping the residue pairs into categories based on chemical similarity, we find that some categories have very non-sinusoidal interplanar angle distributions, whereas other categories have distributions that are close to sinusoidal. For a few categories, observed deviations from a sinusoidal distribution can be explained by the electrostatic anisotropy of the isolated pair potential energy. In other cases, the observed distributions reflect the longer range effects of different possible interaction geometries. In particular, geometries that disrupt external hydrogen bonding are disfavored.

Ideal architecture of residue packing and its observation in protein structures

A simple model of sphere packing has been investigated as an ideal model for long-range interactions for the packing of non-bonded residues in protein structures. By superposing all residues, the geometry of packing around a central residue is investigated. It is found that all residues conform almost perfectly to this lattice model for sphere packing when a radius of 6.5 A is used to define non-bonded (virtual) interacting residues. Side-chain positions with respect to sequential backbone segments are relatively regular as well. This lattice can readily be used in conformation simulations to reduce the conformational space.