The native-state ensemble of proteins provides clues for folding, misfolding and function

The predominant equilibrium in proteins is not between native and unfolded states, it is between the native and multiple partially unfolded forms. Some of these partially unfolded forms can be energetically close to the native state and, therefore, have the potential to become appreciably populated. This could have an important role in protein function or misfolding diseases. The recent identification and characterization of the partially unfolded forms of apoflavodoxin furthers our understanding of their formation.     

Folding of proteins with diverse folds

Using parallel tempering simulations with high statistics, we investigate the folding and thermodynamic properties of three small proteins with distinct native folds: the all-helical 1RIJ, the all-sheet beta3s, and BBA5, which has a mixed helix-sheet fold. In all three cases, simulations with our energy function find the native structures as global minima in free energy at experimentally relevant temperatures. However, the folding process strongly differs for the three molecules, indicating that the folding mechanism is correlated with the form of the native structure.     

Denatured state is critical in determining the properties of model proteins designed on different folds

The thermodynamics of proteins designed on three common folds (SH3, chymotrypsin inhibitor 2 [CI2], and protein G) is studied with a simplified C(alpha) model and compared with the thermodynamics of proteins designed on random-generated folds. The model allows to design sequences to fold within a dRMSD ranging from 1.2 to 4.2 A from the crystallographic native conformation and to study properties that are hard to be measured experimentally. It is found that the denatured state of all of them is not random but is, to different extents, partially structured. The degree of structure is more abundant for SH3 and protein G, giving rise to a weaker stability but a more efficient folding kinetics than CI2 and, even more, than the random-generated folds. Consequently, the features of the unfolded state seem to be as important in the determination of the thermodynamic properties of these proteins as the features of the native state.     

Circular permuted proteins in the universe of protein folds

Finding and identifying circular permuted protein pairs (CPP) is one of the harder tasks for structure alignment programs, because of the different location of the break in the polypeptide chain connectivity. The protein structure alignment tool GANGSTA+ was used to search for CPPs in a database of nearly 10,000 protein structures. It also allows determination of the statistical significance of the occurrence of circular permutations in the protein universe. The number of detected CPPs was found to be higher than expected, raising questions about the evolutionary processes leading to CPPs. The GANGSTA+ protein structure alignment tool is available online via the web server at http://gangsta.chemie.fu‐berlin.de . On the same webpage the complete data base of similar protein structure pairs based on the ASTRAL40 set of protein domains is provided and one can select CPPs specifically. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Role of native-state structure in rubredoxin native-state hydrogen exchange

Amide exchange rates were measured for Pyrococcus furiosus (Pf) rubredoxin substituted with either Zn(II), Ga(III), or Ge(IV). Base-catalyzed exchange rate constants increase up to 3000-fold per unit charge for the highly protected amides surrounding the active site metal, yielding apparent residue-specific conformational energy decreases of more than 8 kcal/mol in a comparison of the Zn(II)- and Ge(IV)-substituted proteins. However, the exchange kinetics for many of the other amides of the protein are insensitive to these metal substitutions. These differential rates are inversely correlated with the distance between the amide nitrogen and the metal in the X-ray structure, out to a distance of at least 12 A, consistent with an electrostatic potential-dependent shifting of the amide nitrogen pK. This strongly correlated distance dependence is consistent with a nativelike structure for the exchange-competent conformations. The electric field potential within the interior of the rubredoxin structure gives rise to a change of as much as a million-fold in the rate for the exchange-competent state of the individual amide hydrogens. Nevertheless, the strength of these electrostatic interactions in Pf rubredoxin appears to be comparable to those previously reported within other proteins. As a result, contrary to the conventional analysis of hydrogen exchange data, for exchange processes that occur via nonglobal transitions, the residual conformational structure will often modulate the observed rates. Although this necessarily complicates the estimation of the conformational equilibria of these exchange-competent states, this dependence on residual structure can provide insight into the conformation of these transient states.     

Effects of sucrose on conformational equilibria and fluctuations within the native-state ensemble of proteins

Osmolytes increase the thermodynamic conformational stability of proteins, shifting the equilibrium between native and denatured states to favor the native state. However, their effects on conformational equilibria within native-state ensembles of proteins remain controversial. We investigated the effects of sucrose, a model osmolyte, on conformational equilibria and fluctuations within the native-state ensembles of bovine pancreatic ribonuclease A and S and horse heart cytochrome c. In the presence of sucrose, the far- and near-UV circular dichroism spectra of all three native proteins were slightly altered and indicated that the sugar shifted the native-state ensemble toward species with more ordered, compact conformations, without detectable changes in secondary structural contents. Thermodynamic stability of the proteins, as measured by guanidine HCl-induced unfolding, increased in proportion to sucrose concentration. Native-state hydrogen exchange (HX) studies monitored by infrared spectroscopy showed that addition of 1 M sucrose reduced average HX rate constants at all degrees of exchange of the proteins, for which comparison could be made in the presence and absence of sucrose. Sucrose also increased the exchange-resistant core regions of the proteins. A coupling factor analysis relating the free energy of HX to the free energy of unfolding showed that sucrose had greater effects on large-scale than on small-scale fluctuations. These results indicate that the presence of sucrose shifts the conformational equilibria toward the most compact protein species within native-state ensembles, which can be explained by preferential exclusion of sucrose from the protein surface.     

Topological accessibility shows a distinct asymmetry in the folds of betaalpha proteins

A novel measure, called "topological accessibility" quantifies how easy it is to reconstruct a protein structure using only local contacts when starting at any point on the chain. Plotting this measure for all points in the chain gives a picture of how accessible the fold is. Simple folds are accessible from all positions, others are accessible only from limited positions while the most complex folds are not accessible from any position. The distribution of topological accessibility along the chain was found to be completely symmetric for the all-alpha and all-beta protein classes. However, for the betaalpha class, a distinct asymmetry was found (with probability 10(-30) of being due to chance). Examination of the proteins contributing to this signal indicated many that have an ancient origin. This suggests that the folds of these proteins may have become fixed under the influence of amino-terminal folding before the advent of chaperone assisted folding.     

Effect of the ventricular folds in a synthetic larynx model

Within the human larynx, the ventricular folds serve primarily as a protecting valve during swallowing. They are located directly above the sound-generating vocal folds. During normal phonation, the ventricular folds are passive structures that are not excited to periodical oscillations. However, the impact of the ventricular folds on the phonation process has not yet been finally clarified.An experimental synthetic human larynx model was used to investigate the effect of the ventricular folds on the phonation process. The model includes self-oscillating vocal fold models and allows the comparison of the pressure distribution at multiple locations in the larynx for configurations with and without ventricular folds.The results indicate that the ventricular folds increase the efficiency of the phonation process by reducing the phonation threshold level of the pressure below the vocal folds. Two effects caused by the ventricular folds could be identified as reasons: (1) a decrease in the mean pressure level in the region between vocal and ventricular folds (ventricles) and (2) an increase in the glottal flow resistance.The reason for the first effect is a reduction of the pressure level in the ventricles due to the jet entrainment and the low static pressure in the glottal jet. The second effect results from an increase in the glottal flow resistance that enhances the aerodynamic energy transfer into the vocal folds. This effect reduces the onset threshold of the pressure difference across the glottis.     

Explicit-chain model of native-state hydrogen exchange: implications for event ordering and cooperativity in protein folding

Native-state hydrogen exchange experiments on several proteins have revealed partially unfolded conformations with diverse stabilities. These equilibrium observations have been used to support kinetic arguments that folding proceeds via a sequential "pathway." This interpretative logic is evaluated here by analyzing the relationship between thermodynamic behavior and folding kinetics in a class of simplified lattice protein models. The chain models studied have varying degrees of cooperative interplay (coupling) between local helical conformational preference and favorable nonlocal interactions. When model cooperativity is high, as native conditions are weakened, "isotherms" of free energy of exchange for residues belonging to the same helix merge together before global unfolding. The point of merger depends on the model energetic favorability of the helix. This trend is similar to the corresponding experimental observations. Kinetically, we find that the ordering of helix formation in the very last stage of native core assembly tends to follow the stabilities of their converged isotherms. In a majority (but not all) of folding trajectories, the final assembly of helices that are thermodynamically more stable against exchange precedes that of helices that are less stable against exchange. These model features are in partial agreement with common experimental interpretations. However, the model results also underscore the ensemble nature of the folding process: the kinetics of helix formation is not a discrete, strictly "all-or-none" process as that envisioned by certain non-explicit-chain models. Helices generally undergo many cycles of partial formation and dissolution before their conformations are fixed in the final assembly stage of folding, a kinetic stage that takes up only approximately 2% of the average folding time in the present model; and the ordering of the helices' final assembly in some trajectories can be different from the dominant ordering stipulated by the exchange isotherms.     

Loop folds in proteins and evolutionary conservation of folding nuclei

We show that loops of close contacts involving hydrophobic residues are important in protein folding. Contrary to Berezovsky 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.     

Geometrical parameters and shape anisotropy in globular proteins

The anisotropy in the shape of globular proteins is derived by a comparative analysis of three types of geometrical parameters. The role of secondary structures in the design of the shape of globular proteins is worked out.     

Relationship between sequence determinants of stability for two natural homologous proteins with different folds

In the Cro protein family, an evolutionary change in secondary structure has converted an alpha-helical fold to a mixture of alpha-helix and beta-sheet. P22 Cro and lambda Cro represent the ancestral all-alpha and descendant alpha+beta folds, respectively. The major structural differences between these proteins are at the C-terminal end of the domain (residues 34-56), where two alpha-helices in P22 Cro align with two beta-strands in lambda Cro. We sought to assess the possibility that smooth evolutionary transitions could have converted the all-alpha structure to the alpha+beta structure through sequences that could adopt both folds. First, we used scanning mutagenesis to identify and compare patterns of key stabilizing residues in the C-terminal regions of both P22 Cro and lambda Cro. These patterns exhibited little similarity to each other, with structurally important residues in the two proteins most often occurring at different sequence positions. Second, "hybrid scanning" studies, involving replacement of each wild-type residue in P22 Cro with the aligned wild-type residue in lambda Cro and vice versa, revealed five or six residues in each protein that strongly destabilized the other. These results suggest that key stability determinants for each Cro fold are quite different and that the P22 Cro sequence strongly favors the all-alpha structure while the lambda Cro sequence strongly favors the alpha+beta structure. Nonetheless, we were able to design a "structurally ambivalent" sequence fragment (SASF1), which corresponded to residues 39-56 and simultaneously incorporated most key stabilizing residues for both P22 Cro and lambda Cro. NMR experiments showed SASF1 to stably fold as a beta-hairpin when incorporated into the lambda Cro sequence but as a pair of alpha-helices when incorporated into P22 Cro.     

Mycobacterial PE_PGRS proteins contain calcium-binding motifs with parallel beta-roll folds

The PE_PGRS family of proteins unique to mycobacteria is demonstrated to contain multiple calcium-binding and glycine-rich sequence motifs GGXGXD/NXUX. This sequence repeat constitutes a calcium-binding parallel beta-roll or parallel beta-helix structure and is found in RTX toxins secreted by many Gram-negative bacteria. It is predicted that the highly homologous PE PGRS proteins containing multiple copies of the nona-peptide motif could fold into similar calcium-binding structures. The implication of the predicted calcium-binding property of PE PGRS proteins in the light of macrophage-pathogen interaction and pathogenesis is presented.     

Folding mechanisms of proteins with high sequence identity but different folds

The problem of how a protein folds from a linear chain of amino acids to the three-dimensional structure necessary for function is often investigated using proteins with a low degree of sequence identity that adopt different folds. The design of pairs of proteins with a high degree of sequence identity but different folds offers the opportunity for a complementary study; in two highly similar sequences, which residues are the most important in directing folding to a particular structure? Here we use molecular dynamics simulations to characterize the folding-unfolding pathways of a pair of proteins designed by Bryan and co-workers [Alexander, P. A., et al. (2005) Biochemistry 44, 14045-14054; He, Y. N., et al. (2005) Biochemistry 44, 14055-14061]. Despite being 59% identical, the two protein sequences fold to two different structures. The first sequence folds to the alpha+beta protein G structure and the second to the all-alpha-helical protein A structure. We show that the final protein structure is determined early along the folding pathway. In folding to the protein G structure, the single alpha-helix (alpha1) and the beta3-beta4 turn fold early. Formation of the hairpin turn essentially prevents folding to helical structure in this region of the protein. This early structure is then consolidated by formation of long-range hydrophobic interactions between alpha1 and the beta3-beta4 turn. The protein A sequence differs both in the residues that form the beta3-beta4 turn and also in many of the residues that form the early hydrophobic interactions in the protein G structure. Instead, in the protein A sequence, a more hierarchical mechanism is observed, with helices folding before many of the tertiary interactions are formed. We find that small, but critical, sequence differences determine the topology of the protein early along the folding pathway, which help to explain the process by which one fold can evolve into another.     

Detecting native folds in mixtures of proteins that contain disulfide bonds

High-throughput in vitro refolding of proteins that contain disulfide bonds, for which soluble expression is particularly difficult, is severely impeded by the absence of effective methods for detecting their native forms. We demonstrate such a method, which combines mass spectrometry with mild reductions, requires no prior experimentation or knowledge of proteins' physicochemical characteristics, function or activity, and is amenable to automation. These are necessary criteria for structural genomics and proteomics applications.     

Similar folds with different stabilization mechanisms: the cases of prion and doppel proteins

Background  

Protein misfolding is the main cause of a group of fatal neurodegenerative diseases in humans and animals. In particular, in Prion-related diseases the normal cellular form of the Prion Protein PrP (PrP ^C ) is converted into the infectious PrP ^Sc through a conformational process during which it acquires a high β-sheet content. Doppel is a protein that shares a similar native fold, but lacks the scrapie isoform. Understanding the molecular determinants of these different behaviours is important both for biomedical and biophysical research.     

Distance geometry generates native-like folds for small helical proteins using the consensus distances of predicted protein structures.

For successful ab initio protein structure prediction, a method is needed to identify native-like structures from a set containing both native and non-native protein-like conformations. In this regard, the use of distance geometry has shown promise when accurate inter-residue distances are available. We describe a method by which distance geometry restraints are culled from sets of 500 protein-like conformations for four small helical proteins generated by the method of Simons et al. (1997). A consensus-based approach was applied in which every inter-Calpha distance was measured, and the most frequently occurring distances were used as input restraints for distance geometry. For each protein, a structure with lower coordinate root-mean-square (RMS) error than the mean of the original set was constructed; in three cases the topology of the fold resembled that of the native protein. When the fold sets were filtered for the best scoring conformations with respect to an all-atom knowledge-based scoring function, the remaining subset of 50 structures yielded restraints of higher accuracy. A second round of distance geometry using these restraints resulted in an average coordinate RMS error of 4.38 A.