Understanding the mechanism of beta-hairpin folding via phi-value analysis

The folding kinetics of a 16-residue beta-hairpin (trpzip4) and five mutants were studied by a laser-induced temperature-jump infrared method. Our results indicate that mutations which affect the strength of the hydrophobic cluster lead to a decrease in the thermal stability of the beta-hairpin, as a result of increased unfolding rates. For example, the W45Y mutant has a phi-value of approximately zero, implying a folding transition state in which the native contacts involving Trp45 are not yet formed. On the other hand, mutations in the turn or loop region mostly affect the folding rate. In particular, replacing Asp46 with Ala leads to a decrease in the folding rate by roughly 9 times. Accordingly, the phi-value for D46A is determined to be approximately 0.77, suggesting that this residue plays a key role in stabilizing the folding transition state. This is most likely due to the fact that the main chain and side chain of Asp46 form a characteristic hydrogen bond network with other residues in the turn region. Taken together, these results support the folding mechanism we proposed before, which suggests that the turn formation is the rate-limiting step in beta-hairpin folding and, consequently, a stronger turn-promoting sequence increases the stability of a beta-hairpin primarily by increasing its folding rate, whereas a stronger hydrophobic cluster increases the stability of a beta-hairpin primarily by decreasing its unfolding rate. In addition, we have examined the compactness of the thermally denatured and urea-denatured states of another 16-residue beta-hairpin, using the method of fluorescence resonance energy transfer. Our results show that the thermally denatured state of this beta-hairpin is significantly more compact than the urea-denatured state, suggesting that the very first step in beta-hairpin folding, when initiated from an extended conformation, probably corresponds to a process of hydrophobic collapse.     

Conformational studies of the C-terminal 16-amino-acid-residue fragment of the B3 domain of the immunoglobulin binding protein G from Streptococcus

The structure and stability of the 16-amino-acid-residue fragment [IG(46-61)] corresponding to the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus was investigated by means of CD and NMR spectroscopy and by differential scanning calorimetry. The CD and 2D NMR experiments were carried out (i) in water at different temperatures and (ii) at one temperature (305 K), with only CD, at different TFE concentrations. Our results show that the IG(46-61) peptide possesses organized three-dimensional structure at all investigated temperatures. The three-dimensional structure of the IG(46-61) peptide resembles the general shape of a beta-hairpin that is also observed for this peptide in the experimental structure of the B3 domain in the whole G protein; the structure is stabilized by hydrophobic interactions between nonpolar side chains. Our study shows that the melting temperature of the IG(46-61) peptide is about 320 K which supports the hypothesis that the investigated peptide can serve as a folding initiation site of the B3 domain of the immunoglobulin binding protein G.     

Context-dependent effects of proline residues on the stability and folding pathway of ubiquitin.

Substitution of trans-proline at three positions in ubiquitin (residues 19, 37 and 38) produces significant context-dependent effects on protein stability (both stabilizing and destabilizing) that reflect changes to a combination of parameters including backbone flexibility, hydrophobic interactions, solvent accessibility to polar groups and intrinsic backbone conformational preferences. Kinetic analysis of the wild-type yeast protein reveals a predominant fast-folding phase which conforms to an apparent two-state folding model. Temperature-dependent studies of the refolding rate reveal thermodynamic details of the nature of the transition state for folding consistent with hydrophobic collapse providing the overall driving force. Br?nsted analysis of the refolding and unfolding rates of a family of mutants with a variety of side chain substitutions for P37 and P38 reveals that the two prolines, which are located in a surface loop adjacent to the C terminus of the main alpha-helix (residues 24-33), are not significantly structured in the transition state for folding and appear to be consolidated into the native structure only late in the folding process. We draw a similar conclusion regarding position 19 in the loop connecting the N-terminal beta-hairpin to the main alpha-helix. The proline residues of ubiquitin are passive spectators in the folding process, but influence protein stability in a variety of ways.     

The mechanism of beta-hairpin formation

Beta-hairpins constitute an important class of connecting protein secondary structures. Several groups have postulated that such structures form early in the folding process and serve to nucleate the formation of extended beta-sheet structures. Despite the importance of beta-hairpins in protein folding, little is known about the mechanism of formation of these structures. While it is well established that there is a complex interplay between the stability of a beta-hairpin and loop conformational propensity, loop length, and the formation of stabilizing cross-strand interactions (H-bonds and hydrophobic interactions), the influence of these factors on the folding rate is poorly understood. Peptide models provide a simple framework for exploring the molecular details of the formation of beta-hairpin structures. We have explored the fundamental processes of folding in two linear peptides that form beta-hairpin structures, having a stabilizing hydrophobic cluster connected by loops of differing lengths. This approach allows us to evaluate existing models of the mechanism of beta-hairpin formation. We find a substantial acceleration of the folding rate when the connecting loop is made shorter (i.e., the hydrophobic cluster is moved closer to the turn). Analysis of the folding kinetics of these two peptides reveals that this acceleration is a direct consequence of the reduced entropic cost of the smaller loop search.     

Three-dimensional structure of horse liver alcohol dehydrogenase at 2-4 A resolution.

The crystal structure analysis of horse liver alcohol dehydrogenase has been extended to 2.4 Å resolution. From the corresponding electron density map of the apoenzyme we have determined the positions of the 374 amino acids in the polypeptide chain of each subunit.The coenzyme binding domain of the subunit comprises residues 176 to 318. 45% of these residues are helical and 32% are in the central six-stranded pleated sheet structure. The positions and orientations of the helices with respect to the pleated sheet indicate a possible folding mechanism for this part of the subunit structure. The coenzyme analogue ADP-ribose binds to this domain in a position and orientation very similar to coenzyme binding to lactate dehydrogenase. The adenine part binds in a hydrophobic pocket, the adenosine ribose is hydrogen-bonded to the side chain of Asp223, the pyrophosphate is positioned by interaction with Arg47 and the nicotinamide ribose is 6Å away from the catalytic zinc atom.The catalytic domain is mainly built up from three distinct antiparallel pleated-sheet regions. Residues within this domain provide ligands to the catalytic zinc atom; Cys46, His67 and Cys174. An approximate tetrahedral coordination of this zinc is completed by a water molecule or hydroxyl ion depending on the pH. Residues 95 to 113 form a lobe that binds the second zinc atom of the subunit. This zinc is liganded in a distorted tetrahedral arrangement by four sulphur atoms from the cysteine residues 97, 100, 103 and 111. The lobe forms one side of a significant cleft in the enzyme surface suggesting that this region might constitute a second catalytic centre of unknown function.The two domains of the subunit are separated by a crevice that contains a wide and deep hydrophobic pocket. The catalytic zinc atom is at the bottom of this pocket, with the zinc-bound water molecule projecting out into the pocket. This water molecule is hydrogen-bonded to the side chain of Ser48 which in turn is hydrogen-bonded to His51. The pocket which in all probability is the binding site for the substrate and the nicotinamide moiety of the coenzyme, is lined almost exclusively with hydrophobic side chains. Both subunits contribute residues to each of the two substrate binding pockets of the molecule. The only accessible polar groups in the vicinity of the catalytic centre are Ser48 and Thr178 apart from zinc and the zinc-bound water molecule.     

Hydrogen bonds between short polar side chains and peptide backbone: prevalence in proteins and effects on helix-forming propensities

A survey of 322 proteins showed that the short polar (SP) side chains of four residues, Thr, Ser, Asp, and Asn, have a very strong tendency to form hydrogen bonds with neighboring backbone amides. Specifically, 32% of Thr, 29% of Ser, 26% of Asp, and 19% of Asn engage in such hydrogen bonds. When an SP residue caps the N terminal of a helix, the contribution to helix stability by a hydrogen bond with the amide of the N3 or N2 residue is well established. When an SP residue is in the middle of a helix, the side chain is unlikely to form hydrogen bonds with neighboring backbone amides for steric and geometric reasons. In essence the SP side chain competes with the backbone carbonyl for the same hydrogen-bonding partner (i.e., the backbone amide) and thus SP residues tend to break backbone carbonyl-amide hydrogen bonds. The proposition that this is the origin for the low propensities of SP residues in the middle of alpha helices (relative to those of nonpolar residues) was tested. The combined effects of restricting side-chain rotamer conformations (documented by Creamer and Rose, Proc Acad Sci USA, 1992;89:5937-5941; Proteins, 1994;19:85-97) and excluding side- chain to backbone hydrogen bonds by the helix were quantitatively analyzed. These were found to correlate strongly with four experimentally determined scales of helix-forming propensities. The correlation coefficients ranged from 0.72 to 0.87, which are comparable to those found for nonpolar residues (for which only the loss of side-chain conformational entropy needs to be considered).     

Molecular dynamics simulations of a beta-hairpin fragment of protein G: balance between side-chain and backbone forces

How is the native structure encoded in the amino acid sequence? For the traditional backbone centric view, the dominant forces are hydrogen bonds (backbone) and phi-psi propensity. The role of hydrophobicity is non-specific. For the side-chain centric view, the dominant force of protein folding is hydrophobicity. In order to understand the balance between backbone and side-chain forces, we have studied the contributions of three components of a beta-hairpin peptide: turn, backbone hydrogen bonding and side-chain interactions, of a 16-residue fragment of protein G. The peptide folds rapidly and cooperatively to a conformation with a defined secondary structure and a packed hydrophobic cluster of aromatic side-chains. Our strategy is to observe the structural stability of the beta-hairpin under systematic perturbations of the turn region, backbone hydrogen bonds and the hydrophobic core formed by the side-chains, respectively. In our molecular dynamics simulations, the peptides are solvated. with explicit water molecules, and an all-atom force field (CFF91) is used. Starting from the original peptide (G41EWTYDDATKTFTVTE56), we carried out the following MD simulations. (1) unfolding at 350 K; (2) forcing the distance between the C(alpha) atoms of ASP47 and LYS50 to be 8 A; (3) deleting two turn residues (Ala48 and Thr49) to form a beta-sheet complex of two short peptides, GEWTYDD and KTFTVTE; (4) four hydrophobic residues (W43, Y45, F52 and T53) are replaced by a glycine residue step-by-step; and (5) most importantly, four amide hydrogen atoms (T44, D46, T53, and T55, which are crucial for backbone hydrogen bonding), are substituted by fluorine atoms. The fluorination not only makes it impossible to form attractive hydrogen bonding between the two beta-hairpin strands, but also introduces a repulsive force between the two strands due to the negative charges on the fluorine and oxygen atoms. Throughout all simulations, we observe that backbone hydrogen bonds are very sensitive to the perturbations and are easily broken. In contrast, the hydrophobic core survives most perturbations. In the decisive test of fluorination, the fluorinated peptide remains folded under our simulation conditions (5 ns, 278 K). Hydrophobic interactions keep the peptide folded, even with a repulsive force between the beta-strands. Thus, our results strongly support a side-chain centric view for protein folding.     

A set of homology models of pore loop domain of six eukaryotic voltage-gated potassium channels Kv1.1-Kv1.6

Homology models of the pore loop domain of six eukaryotic potassium channels, Kv1.1-Kv1.6, were generated based on the crystallographic structure of KcsA. The results of amino acid sequence alignment indicate that these Kv channels are composed of two structurally and functionally independent domains: the N-terminal 'voltage sensor' domain and the C-terminal 'pore loop' domain. The homology models reveal that the pore loop domains of these Kv channels exhibit similar folds to those of KcsA. The structural features and specific packing of aromatic residues around the selectivity filter of these Kv channels are nearly identical to those of KcsA, whereas most of the structural variations occur in the turret as well as in the inner and outer helices. The distribution of polar and nonpolar side chains on the surfaces of the KcsA and Kv channels reveals that they exhibit a segregation of side chains common to most integral membrane proteins. As the hydrogen bond between Glu71 and Asp80 in KcsA plays an important role in stabilizing the channel, the substituted Val residue in the Kv family corresponding to Glu71 of KcsA stabilizes the channel by making hydrophobic contact with Tyr residue from the signature sequence of the selectivity filter. The homology models of these Kv channels provide particularly attractive subjects for further structure-based studies.     

A novel tripartite motif involved in aquaporin topogenesis, monomer folding and tetramerization

Aquaporin (AQP) folding in the endoplasmic reticulum is characterized by two distinct pathways of membrane insertion that arise from divergent residues within the second transmembrane segment. We now show that in AQP1 these residues (Asn49 and Lys51) interact with Asp185 at the C terminus of TM5 to form a polar, quaternary structural motif that influences multiple stages of folding. Asn49 and Asp185 form an intramolecular hydrogen bond needed for proper helical packing, monomer formation and function. In contrast, Lys51 interacts with Asp185 on an adjacent monomer to stabilize the AQP1 tetramer. Although these residues are unique to AQP1, they share a highly conserved architecture whose functional properties can be transferred to other family members. These findings suggest a general mechanism by which evolutionary divergence of membrane proteins can confer new functional properties via alternative folding pathways that give rise to a common final structure.     

Structural analysis of the N-terminal domain of the human T-cell leukemia virus capsid protein

The N-terminal domain of the retroviral capsid (CA) protein is one of the least conserved regions encoded in the genome. Surprisingly, the three-dimensional structures of the CA from different genera exhibit alpha-helical structural features that are highly conserved. The N-terminal residues of the human immunodeficiency virus type 1 (HIV-1) and Rous sarcoma virus (RSV) capsid proteins form a beta-hairpin. To determine if this feature is conserved in the retroviral family, we cloned, expressed, purified, and solved the structure of a N-terminal 134 amino acid fragment (CA(134)) from the human T-cell leukemia virus type 1 (HTLV-I) using high resolution nuclear magnetic resonance (NMR) spectroscopy. The CA(134) fragment contains an N-terminal beta-hairpin and a central coiled-coil-like structure composed of six alpha-helices. The N-terminal Pro1 residue contacts Asp54 in the helical cluster through a salt bridge. Thus, the beta-hairpin is conserved and the helical cluster is structurally similar to other retroviral CA domains. However, although the same Asp residue defines the orientation of the hairpin in both the HTLV-1 and HIV-1 CA proteins, the HTLV-I hairpin is oriented away, rather than towards, the helical core. Significant differences were also detected in the spatial orientation and helical content of the long centrally located loop connecting the helices in the core. It has been proposed that the salt bridge allows the formation of a CA-CA interface that is important for the assembly of the conical cores that are characteristic of HIV-1. As HTLV-I forms spherical cores, the salt-bridge feature is apparently not conserved for this function although its role in determining the orientation of the beta-hairpin may be critical, along with the central loop. Comparison of three-dimensional structures is expected to elucidate the relationships between the retroviral capsid protein structure and its function.     

Structural basis for the unusual specificity of Escherichia coli aminopeptidase N

Aminopeptidase N from Escherichia coli is a M1 class aminopeptidase with the active-site region related to that of thermolysin. The enzyme has unusual specificity, cleaving adjacent to the large, nonpolar amino acids Phe and Tyr but also cleaving next to the polar residues Lys and Arg. To try to understand the structural basis for this pattern of hydrolysis, the structure of the enzyme was determined in complex with the amino acids L-arginine, L-lysine, L-phenylalanine, L-tryptophan, and L-tyrosine. These amino acids all bind with their backbone atoms close to the active-site zinc ion and their side chain occupying the S1 subsite. This subsite is in the form of a cylinder, about 10 A in cross-section and 12 A in length. The bottom of the cylinder includes the zinc ion and a number of polar side chains that make multiple hydrogen-bonding and other interactions with the alpha-amino group and the alpha-carboxylate of the bound amino acid. The walls of the S1 cylinder are hydrophobic and accommodate the nonpolar or largely nonpolar side chains of Phe and Tyr. The top of the cylinder is polar in character and includes bound water molecules. The epsilon-amino group of the bound lysine side chain and the guanidinium group of arginine both make multiple hydrogen bonds to this part of the S1 site. At the same time, the hydrocarbon part of the lysine and arginine side chains is accommodated within the nonpolar walls of the S1 cylinder. This combination of hydrophobic and hydrophilic binding surfaces explains the ability of ePepN to cleave Lys, Arg, Phe, and Tyr. Another favored substrate has Ala at the P1 position. The short, nonpolar side chain of this residue can clearly be bound within the hydrophobic part of the S1 cylinder, but the reason for its facile hydrolysis remains uncertain.     

Role of hydrophilic and hydrophobic contacts in folding of the second beta-hairpin fragment of protein G: molecular dynamics simulation studies of an all-atom model

Predicting the folding mechanism of the second beta-hairpin fragment of the Ig-binding domain B of streptococcal protein G is unexpectedly challenging for simplified reduced models because the models developed so far indicated a different folding mechanism from what was suggested from high-temperature unfolding and equilibrium free-energy surface analysis based on established all-atom empirical force fields in explicit or implicit solvent. This happened despite the use of empirical residue-based interactions, multibody hydrophobic interactions, and inclusions of hydrogen bonding effects in the simplified models. This article employs a recently developed all-atom (except nonpolar hydrogens) model interacting with simple square-well potentials to fold the peptide fragment by molecular dynamics simulation methods. In this study, 193 out of 200 trajectories are folded at two reduced temperatures (3.5 and 3.7) close to the transition temperature T* approximately 4.0. Each simulation takes <7 h of CPU time on a Pentium 800-MHz PC. Folding of the new all-atom model is found to be initiated by collapse before the formation of main-chain hydrogen bonds. This verifies the mechanism proposed from previous all-atom unfolding and equilibrium simulations. The new model further predicts that the collapse is initiated by two nucleation contacts (a hydrophilic contact between D46 and T49 and a hydrophobic contact between Y45 and F52), in agreement with recent NMR measurements. The results suggest that atomic packing and native contact interactions play a dominant role in folding mechanism.     

Refined 1.2 A crystal structure of the complex formed between subtilisin Carlsberg and the inhibitor eglin c. Molecular structure of eglin and its detailed interaction with subtilisin.

The crystal structure of the complex formed between eglin c, an elastase inhibitor from the medical leech, and subtilisin Carlsberg has been determined at 1.2 A resolution by a combination of Patterson search methods and isomorphous replacement techniques. The structure has been refined to a crystallographic R-value of 0.18 (8-1.2 A). Eglin consists of a four-stranded beta-sheet with an alpha-helical segment and the protease-binding loop fixed on opposite sides. This loop, which contains the reactive site Leu45I--Asp46I, is mainly held in its conformation by unique electrostatic/hydrogen bond interactions of Thr44I and Asp46I with the side chains of Arg53I and Arg51I which protrude from the hydrophobic core of the molecule. The conformation around the reactive site is similar to that found in other proteinase inhibitors. The nine residues of the binding loop Gly40I--Arg48I are involved in direct contacts with subtilisin. In this interaction, eglin segment Pro42I--Thr44I forms a three-stranded anti-parallel beta-sheet with subtilisin segments Gly100--Gly102 and Ser125--Gly127. The reactive site peptide bond of eglin is intact, and Ser221 OG of the enzyme is 2.81 A apart from the carbonyl carbon.     

Conformation of the diphtheria toxin T domain in membranes: a site-directed spin-labeling study of the TH8 helix and TL5 loop.

The isolated T domain of diphtheria toxin was mutated by cysteine-scanning mutagenesis at 28 consecutive sites (residues 328-355) that comprise the TH8 helix and the TL5 interhelical loop in the native toxin. After derivatizing the mutant proteins with a sulfhydryl-selective nitroxide reagent, we examined the mobility of each nitroxide and its accessibility to polar and nonpolar paramagnetic reagents, before and after insertion into phospholipid bilayers. The data obtained with the proteins in solution at pH 8 are generally consistent with predictions from the crystal structure of the toxin. Upon membrane binding at pH 4.6, a major structural reorganization of the domain was seen, which dramatically reduced the accessibility of most residues in this region to the polar reagent nickel(II)-ethylenediaminediacetate complex (NiEDDA). Many of these residues also showed reduced accessibility to the nonpolar reagent O(2). Periodic accessibility of the nitroxide side chains along the sequence to these reagents shows that TH8 remains largely helical in the membrane-bound state, with one surface associated with protein and the other facing the hydrophobic interior of the bilayer. In addition, the TL5 loop also appears to become alpha-helical in the membrane, with one surface in contact with protein and the other in contact with the bilayer interior. These findings provide a structural framework for understanding how the T domain forms a transmembrane channel and mediates translocation of diphtheria toxin's enzymic moiety across a membrane.     

Comparison of the hemocyanin beta-barrel with other Greek key beta-barrels: possible importance of the "beta-zipper" in protein structure and folding.

The Greek key beta-barrel topology is a folding motif observed in many proteins of widespread evolutionary origin. The arthropodan hemocyanins also have such a Greek key beta-barrel, which forms the core of the third domain of this protein. The hemocyanin beta-barrel was found to be structurally very similar to the beta-barrels of the immunoglobulin domains, Cu,Zn-superoxide dismutase and the chromophore carrying antitumor proteins. The structural similarity within this group of protein families is not accompanied by an evolutionary or functional relationship. It is therefore possible to study structure-sequence relations without bias from nonstructural constraints. The present study reports a conserved pattern of features in these Greek key beta-barrels that is strongly suggestive of a folding nucleation site. This proposed nucleation site, which we call a "beta-zipper," shows a pattern of well-conserved, large hydrophobic residues on two sequential beta-strands joined by a short loop. Each beta-zipper strand is near the center of one of the beta-sheets, so that the two strands face each other from opposite sides of the barrel and interact through their hydrophobic side chains, rather than forming a hydrogen-bonded beta-hairpin. Other protein families with Greek key beta-barrels that do not as strongly resemble the immunoglobulin fold--such as the azurins, plastocyanins, crystallins, and prealbumins--also contain the beta-zipper pattern, which might therefore be a universal feature of Greek key beta-barrel proteins.     

Secondary structures without backbone: an analysis of backbone mimicry by polar side chains in protein structures.

Backbone mimicry by the formation of closed-loop C7, C10 and C13 (mimics of gamma-, beta- and alpha-turns) conformations through side chain-main chain hydrogen bonds by polar groups is a frequent observation in protein structures. A data set of 250 non-homologous and high-resolution protein crystal structures was used to analyze these conformations for their characteristic features. Seven out of the nine polar residues (Ser, Thr, Asn, Asp, Gln, Glu and His) have hydrogen bonding groups in their side chains which can participate in such mimicry and as many as 15% of all these polar residues engage in such conformations. The distributions of dihedral angles of these mimics indicate that only certain combinations of the dihedral angles involved aid the formation of these mimics. The observed examples were categorized into various classes based on these combinations, resulting in well defined motifs. Asn and Asp residues show a very high capability to perform such backbone secondary structural mimicry. The most highly mimicked backbone structure is of the C10 conformation by the Asx residues. The mimics formed by His, Ser, Thr and Glx residues are also discussed. The role of such conformations in initiating the formation of regular secondary structures during the course of protein folding seems significant.     

Sequential unfolding of the hemolysin two‐partner secretion domain from Proteus mirabilis

Protein secretion is a major contributor to Gram‐negative bacterial virulence. Type Vb or two‐partner secretion (TPS) pathways utilize a membrane bound β‐barrel B component (TpsB) to translocate large and predominantly virulent exoproteins (TpsA) through a nucleotide independent mechanism. We focused our studies on a truncated TpsA member termed hemolysin A (HpmA265), a structurally and functionally characterized TPS domain from Proteus mirabilis. Contrary to the expectation that the TPS domain of HpmA265 would denature in a single cooperative transition, we found that the unfolding follows a sequential model with three distinct transitions linking four states. The solvent inaccessible core of HpmA265 can be divided into two different regions. The C‐proximal region contains nonpolar residues and forms a prototypical hydrophobic core as found in globular proteins. The N‐proximal region of the solvent inaccessible core, however, contains polar residues. To understand the contributions of the hydrophobic and polar interiors to overall TPS domain stability, we conducted unfolding studies on HpmA265 and site‐specific mutants of HpmA265. By correlating the effect of individual site‐specific mutations with the sequential unfolding results we were able to divide the HpmA265 TPS domain into polar core, nonpolar core, and C‐terminal subdomains. Moreover, the unfolding studies provide quantitative evidence that the folding free energy for the polar core subdomain is more favorable than for the nonpolar core and C‐terminal subdomains. This study implicates the hydrogen bonds shared among these conserved internal residues as a primary means for stabilizing the N‐proximal polar core subdomain.     

Amino-acid solvation structure in transmembrane helices from molecular dynamics simulations

Understanding the solvation of amino acids in biomembranes is an important step to better explain membrane protein folding. Several experimental studies have shown that polar residues are both common and important in transmembrane segments, which means they have to be solvated in the hydrophobic membrane, at least until helices have aggregated to form integral proteins. In this work, we have used computer simulations to unravel these interactions on the atomic level, and classify intramembrane solvation properties of amino acids. Simulations have been performed for systematic mutations in poly-Leu helices, including not only each amino acid type, but also every z-position in a model helix. Interestingly, many polar or charged residues do not desolvate completely, but rather retain hydration by snorkeling or pulling in water/headgroups--even to the extent where many of them exist in a microscopic polar environment, with hydration levels corresponding well to experimental hydrophobicity scales. This suggests that even for polar/charged residues a large part of solvation cost is due to entropy, not enthalpy loss. Both hydration level and hydrogen bonding exhibit clear position-dependence. Basic side chains cause much less membrane distortion than acidic, since they are able to form hydrogen bonds with carbonyl groups instead of water or headgroups. This preference is supported by sequence statistics, where basic residues have increased relative occurrence at carbonyl z-coordinates. Snorkeling effects and N-/C-terminal orientation bias are directly observed, which significantly reduces the effective thickness of the hydrophobic core. Aromatic side chains intercalate efficiently with lipid chains (improving Trp/Tyr anchoring to the interface) and Ser/Thr residues are stabilized by hydroxyl groups sharing hydrogen bonds to backbone oxygens.     

Folding specificity induced by loop stiffness

To test the importance of loop stiffness in restricting the heterogeneity of transition state ensemble, we relaxed the distal loop of 10 unstable redesigned hydrophobic core mutants of alpha-spectrin SH3 domain. This was achieved by replacing Asp48 by Gly at the tip of the distal hairpin. Although the change was local in nature, the effect on stabilization was not uniform across the core mutants tested. There is an inverse rough correlation between the stabilization and the increase in buried hydrophobic volume, with respect to the wild type. Interestingly enough, proteins that although unstable are properly folded become molten globule-like after relaxation of the distal loop. These results highlight the importance of stiffness in restricting the conformational heterogeneity of a protein during the folding reaction. An interplay between unspecific hydrophobic interactions and constraint induced by polar interactions, or in this case local stiffness, is essential to achieve a well-ordered folded structure.     

Structural characterization of a mutant peptide derived from ubiquitin: implications for protein folding

The formation of the N-terminal beta-hairpin of ubiquitin is thought to be an early event in the folding of this small protein. Previously, we have shown that a peptide corresponding to residues 1-17 of ubiquitin folds autonomously and is likely to have a native-like hairpin register. To investigate the causes of the stability of this fold, we have made mutations in the amino acids at the apex of the turn. We find that in a peptide where Thr9 is replaced by Asp, U(1-17)T9D, the native conformation is stabilized with respect to the wild-type sequence, so much so that we are able to characterize the structure of the mutant peptide fully by NMR spectroscopy. The data indicate that U(1-17)T9D peptide does indeed form a hairpin with a native-like register and a type I turn with a G1 beta-bulge, as in the full-length protein. The reason for the greater stability of the U(1-17)T9D mutant remains uncertain, but there are nuclear Overhauser effects between the side chains of Asp9 and Lys 11, which may indicate that a charge-charge interaction between these residues is responsible.