The relationship between sequence and structure in elementary folding units

Computational scrutiny of the effect of N-terminal proline and residue stereochemistry in the nucleation of α-helix fold

N-Terminal l- to d-residue mutation nucleate helical fold in Ac–^DAla–^LAla₃–NHMe (Ib, m2), Ac–^DPro–^LAla₃–NHMe (IIb, m1), and Ac–^DPro–^LPro–^LAla₂–NHMe (IIIb, m2) peptides.

Design of monomeric water-soluble β-hairpin and β-sheet peptides

Since the first report in 1993 (JACS 115, 5887-5888) of a peptide able to form a monomeric β-hairpin structure in aqueous solution, the design of peptides forming either β-hairpins (two-stranded antiparallel β-sheets) or three-stranded antiparallel β-sheets has become a field of growing interest and activity. These studies have yielded great insights into the principles governing the stability and folding of β-hairpins and antiparallel β-sheets. This chapter provides an overview of the reported β-hairpin/β-sheet peptides focussed on the applied design criteria, reviews briefly the factors contributing to β-hairpin/β-sheet stability, and describes a protocol for the de novo design of β-sheet-forming peptides based on them. Guidelines to select appropriate turn and strand residues and to avoid self-association are provided. The methods employed to check the success of new designed peptides are also summarized. Since NMR is the best technique to that end, NOEs and chemical shifts characteristic of β-hairpins and three-stranded antiparallel β-sheets are given.

A transient α-helical molecular recognition element in the disordered N-terminus of the Sgs1 helicase is critical for chromosome stability and binding of Top3/Rmi1

The RecQ-like DNA helicase family is essential for the maintenance of genome stability in all organisms. Sgs1, a member of this family in Saccharomyces cerevisiae, regulates early and late steps of double-strand break repair by homologous recombination. Using nuclear magnetic resonance spectroscopy, we show that the N-terminal 125 residues of Sgs1 are disordered and contain a transient α-helix that extends from residue 25 to 38. Based on the residue-specific knowledge of transient secondary structure, we designed proline mutations to disrupt this α-helix and observed hypersensitivity to DNA damaging agents and increased frequency of genome rearrangements. In vitro binding assays show that the defects of the proline mutants are the result of impaired binding of Top3 and Rmi1 to Sgs1. Extending mutagenesis N-terminally revealed a second functionally critical region that spans residues 9–17. Depending on the position of the proline substitution in the helix functional impairment of Sgs1 function varied, gradually increasing from the C- to the N-terminus. The multiscale approach we used to interrogate structure/function relationships in the long disordered N-terminal segment of Sgs1 allowed us to precisely define a functionally critical region and should be generally applicable to other disordered proteins.

Study protein folding and aggregation using nonnatural amino acid p-cyanophenylalanine as a sensitive optical probe

Incorporation of nonnatural amino acids with a variety of special side groups into protein sequences has substantially expanded the experimental means of exploring protein structures and functions. Recently, p-cyanophenylalanine (PheCN), the nitrile analogue of phenylalanine, has been used as a novel optical probe for protein binding and folding studies. The fluorescence emission of PheCN is sensitive to solvent and local environment of the residue, making it a useful fluorescent probe of protein structural change at residue-specific resolution. Moreover, the utility of PheCN is increased by its ability to excite tryptophan fluorescence via the mechanism of fluorescence resonance energy transfer. PheCN could be applied to study a variety of biological problems, e.g., protein folding/unfolding and protein aggregation.

Dependence of α-helical and β-sheet amino acid propensities on the overall protein fold type

Background

A large number of studies have been carried out to obtain amino acid propensities for α-helices and β-sheets. The obtained propensities for α-helices are consistent with each other, and the pair-wise correlation coefficient is frequently high. On the other hand, the β-sheet propensities obtained by several studies differed significantly, indicating that the context significantly affects β-sheet propensity.

Results

We calculated amino acid propensities for α-helices and β-sheets for 39 and 24 protein folds, respectively, and addressed whether they correlate with the fold. The propensities were also calculated for exposed and buried sites, respectively. Results showed that α-helix propensities do not differ significantly by fold, but β-sheet propensities are diverse and depend on the fold. The propensities calculated for exposed sites and buried sites are similar for α-helix, but such is not the case for the β-sheet propensities. We also found some fold dependence on amino acid frequency in β-strands. Folds with a high Ser, Thr and Asn content at exposed sites in β-strands tend to have a low Leu, Ile, Glu, Lys and Arg content (correlation coefficient = −0.90) and to have flat β-sheets. At buried sites in β-strands, the content of Tyr, Trp, Gln and Ser correlates negatively with the content of Val, Ile and Leu (correlation coefficient = −0.93). "All-β" proteins tend to have a higher content of Tyr, Trp, Gln and Ser, whereas "α/β" proteins tend to have a higher content of Val, Ile and Leu.

Conclusions

The α-helix propensities are similar for all folds and for exposed and buried residues. However, β-sheet propensities calculated for exposed residues differ from those for buried residues, indicating that the exposed-residue fraction is one of the major factors governing amino acid composition in β-strands. Furthermore, the correlations we detected suggest that amino acid composition is related to folding properties such as the twist of a β-strand or association between two β sheets.

Exploring the folding free energy landscape of a β-hairpin miniprotein, chignolin, using multiscale free energy landscape calculation method

The folding process for a β-hairpin miniprotein, chignolin, was investigated by free energy landscape (FEL) calculations using the recently proposed multiscale free energy landscape calculation method (MSFEL). First, coarse-grained molecular dynamics simulations searched a broad conformational space, then multiple independent, all-atom molecular dynamics simulations with explicit solvent determined the detailed local FEL using massively distributed computing. The combination of the two models enabled efficient calculation of the free energy landscapes. The MSFEL analysis showed that chignolin has an intermediate state as well as a misfolded state. The folding process is initiated by the formation of a β-hairpin turn, followed by the formation of contacts in the hydrophobic core between Tyr2 and Trp9. Furthermore, mutation of Tyr2 shifts the population to the misfolded conformation. The results indicate that the hydrophobic core plays an important role in stabilizing the native state of chignolin.

Beta-strand flipping and slipping triggered by turn replacement reveal the opportunistic nature of beta-strand pairing

We investigated how the register between adjacent beta-strands is specified using a series of mutants of the single-layer beta-sheet (SLB) in Borrelia OspA. The single-layer architecture of this system eliminates structural restraints imposed by a hydrophobic core, enabling us to address this question. A critical turn (turn 9/10) in the SLB was replaced with a segment with an intentional structural mismatch. Its crystal structure revealed a one-residue insertion into the central beta-strand (strand 9) of the SLB. This insertion triggered a surprisingly large-scale structural rearrangement: (i) the central strand (strand 9) was shifted by one residue, causing the strand to flip with respect to the adjacent beta-strands and thus completely disrupting the native side-chain contacts; (ii) the three-residue turn located on the opposite end of the beta-strand (turn 8/9) was pushed into its preceding beta-strand (strand 8); (iii) the register between strands 8 and 9 was shifted by three residues. Replacing the original sequence for turn 8/9 with a stronger turn motif restored the original strand register but still with a flipped beta-strand 9. The stability differences of these distinct structures were surprisingly small, consistent with an energy landscape where multiple low-energy states with different beta-sheet configurations exist. The observed conformations can be rationalized in terms of maximizing the number of backbone H-bonds. These results suggest that adjacent beta-strands "stick" through the use of factors that are not highly sequence specific and that beta-strands could slide back and forth relatively easily in the absence of external elements such as turns and tertiary packing.

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.

Folding amphipathic helices into membranes: amphiphilicity trumps hydrophobicity

High amphiphilicity is a hallmark of interfacial helices in membrane proteins and membrane-active peptides, such as toxins and antimicrobial peptides. Although there is general agreement that amphiphilicity is important for membrane-interface binding, an unanswered question is its importance relative to simple hydrophobicity-driven partitioning. We have examined this fundamental question using measurements of the interfacial partitioning of a family of 17-residue amidated-acetylated peptides into both neutral and anionic lipid vesicles. Composed only of Ala, Leu, and Gln residues, the amino acid sequences of the peptides were varied to change peptide amphiphilicity without changing total hydrophobicity. We found that peptide helicity in water and interface increased linearly with hydrophobic moment, as did the favorable peptide partitioning free energy. This observation provides simple tools for designing amphipathic helical peptides. Finally, our results show that helical amphiphilicity is far more important for interfacial binding than simple hydrophobicity.

Hydrophobic surface burial is the major stability determinant of a flat, single-layer beta-sheet

Formation of a flat beta-sheet is a fundamental event in beta-sheet-mediated protein self-assembly. To investigate the contributions of various factors to the stability of flat beta-sheets, we performed extensive alanine-scanning mutagenesis experiments on the single-layer beta-sheet segment of Borrelia outer surface protein A (OspA). This beta-sheet segment consists of beta-strands with highly regular geometries that can serve as a building block for self-assembly. Our Ala-scanning approach is distinct from the conventional host-guest method, in that it introduces only conservative, truncation mutations that should minimize structural perturbation. Our results showed very weak correlation with experimental beta-sheet propensity scales, statistical beta-sheet propensity scales, or cross-strand pairwise correlations. In contrast, our data showed strong positive correlation with the change in buried non-polar surface area. Polar interactions including prominent Glu-Lys cross-strand pairs contribute marginally to the beta-sheet stability. These results were corroborated by results from additional non-Ala mutations. Taken together, these results demonstrate the dominant contribution of non-polar surface burial to flat beta-sheet stability even at solvent-exposed positions. The OspA single-layer beta-sheet achieves efficient hydrophobic surface burial without forming a hydrophobic core by a strategic placement of a variety of side-chains. These findings further suggest the importance of hydrophobic interactions within a beta-sheet layer in peptide self-assembly.

Characterization of β-turn and Asx-turns mimicry in a model peptide: Stabilization via CH · · · O interaction

The chemical synthesis and single-crystal X-ray diffraction analysis of a model peptide, Boc-Thr-Thr-NH2 (1) comprised of proteinogenic residues bearing an amphiphilic Cbeta -stereogenic center, has been described. Interestingly, the analysis of its molecular structure revealed the existence of a distinct conformation that mimics a typical beta-turn and Asx-turns, i.e., the two Thr residues occupy the left- and right-corner positions. The main-chain torsion angles of the N- and C-terminal residues i.e., semiextended: phi = -68.9 degrees , psi = 128.6 degrees ; semifolded: phi = -138.1 degrees , psi = 2.5 degrees conformations, respectively, in conjunction with a gauche- disposition of the obligatory C-terminus Thr CgammaH3 group, characterize the occurrence of the newly described beta-turn- and Asx-turns-like topology. The preferred molecular structure is suggested to be stabilized by an effective nonconventional main-chain to side-chain Ci=O . . . H--Cgamma(i+2)-type intraturn hydrogen bond. Noteworthy, the observed topology of the resulting 10-membered hydrogen-bonded ring is essentially similar to the one perceived for a classical beta-turn and the Asx-turns, stabilized by a conventional intraturn hydrogen bond. Considering the signs as well as magnitudes of the backbone torsion angles and the orientation of the central peptide bond, the overall mimicked topology resembles the type II beta-turn or type II Asx-turns. An analysis of Xaa-Thr sequences in high-resolution X-ray elucidated protein structures revealed the novel topology prevalence in functional proteins (unpublished). In view of indubitable structural as well as functional importance of nonconventional interactions in bioorganic and biomacromolecules, we intend to highlight the participation of Thr CgammaH in the creation of a short-range C=O . . . H--Cgamma -type interaction in peptides and proteins.

On the Influence of Charged Side Chains on the Folding–Unfolding Equilibrium of β-Peptides: A Molecular Dynamics Simulation Study

The influence of charged side chains on the folding-unfolding equilibrium of beta-peptides was investigated by means of molecular dynamics simulations. Four different peptides containing only negatively charged side chains, positively charged side chains, both types of charged side chains (with the ability to form stabilizing salt bridges) or no charged side chains were studied under various conditions (different simulation temperatures, starting structures and solvent environment). The NMR solution structure in methanol of one of the peptides (A) has already been published; the synthesis and NMR analysis of another peptide (B) is described here. The other peptides (C and D) studied herein have hitherto not been synthesized. All four peptides A-D are expected to adopt a left-handed 3(14)-helix in solution as well as in the simulations. The resulting ensembles of structures were analyzed in terms of conformational space sampled by the peptides, folding behavior, structural properties such as hydrogen bonding, side chain-side chain and side chain-backbone interactions and in terms of the level of agreement with the NMR data available for two of the peptides. It was found that the presence of charged side chains significantly slows down the folding process in methanol solution due to the stabilization of intermediate conformers with side chain-backbone interactions. In water, where the solvent competes with the solute-solute polar interactions, the folding process to the 3(14)-helix is faster in the simulations.

Peptide fragment studies on the folding elements of dihydrofolate reductase from Escherichia coli

One of the necessary conditions for a protein to be foldable is the presence of a complete set of "folding elements" (FEs) that are short, contiguous peptide segments distributed over an amino acid sequence. The FE-assembly model of protein folding has been proposed, in which the FEs play a role in guiding structure formation through FE-FE interactions early in folding. However, two major issues remain to be clarified regarding the roles of the FEs in determining protein foldability. Are the FEs AFUs that can form nativelike structures in isolation? Is the presence of only the FEs without mutual connections a sufficient condition for a protein to be foldable? Here, we address these questions using peptide fragments corresponding to the FEs of DHFR from Escherichia coli. We show by CD measurement that the FE peptides are unfolded under the native conditions, and some of them have the propensities toward non-native helices. MD simulations also show the non-native helical propensities of the peptides, and the helix contents estimated from the simulations are well correlated with those estimated from the CD in TFE. Thus, the FEs of DHFR are not AFUs, suggesting the importance of the FEs in nonlocal interactions. We also show that equimolar mixtures of the FE peptides do not induce any structural formation. Therefore, mutual connections between the FEs, which should strengthen the nonlocal FE-FE interactions, are also one of the necessary conditions for a protein to be foldable.

β-Hairpin formation in aqueous solution and in the presence of trifluoroethanol: A1H and13C nuclear magnetic resonance conformational study of designed peptides

In order to check our current knowledge on the principles involved in beta-hairpin formation, we have modified the sequence of a 3:5 beta-hairpin forming peptide with two different purposes, first to increase the stability of the formed 3:5 beta-hairpin, and second to convert the 3:5 beta-hairpin into a 2:2 beta-hairpin. The conformational behavior of the designed peptides was investigated in aqueous solution and in 30% trifluoroethanol (TFE) by analysis of the following nuclear magnetic resonance (NMR) parameters: nuclear Overhauser effect (NOE) data, and C(alpha)H, (13)C(alpha), and (13)C(beta) conformational shifts. From the differences in the ability to adopt beta-hairpin structures in these peptides, we have arrived to the following conclusions: (i) beta-Hairpin population increases with the statistical propensity of residues to occupy each turn position. (ii) The loop length, and in turn, the beta-hairpin type, can be modified as a function of the type of turn favored by the loop sequence. These two conclusions reinforce previous results about the importance of beta-turn sequence in beta-hairpin folding. (iii) Side-chain packing on each face of the beta-sheet may play a major role in beta-hairpin stability; hence simplified analysis in terms of isolated pair interactions and intrinsic beta-sheet propensities is insufficient. (iv) Contributions to beta-hairpin stability of turn and strand sequences are not completely independent. (v) The burial of hydrophobic surface upon beta-hairpin formation that, in turn, depends on side-chain packing also contributes to beta-hairpin stability. (vi) As previously observed, TFE stabilizes beta-hairpin structures, but the extent of the contribution of different factors to beta-hairpin formation is sometimes different in aqueous solution and in 30% TFE.

A Comparative Study of Amyloid Fibril Formation by Residues 15–19 of the Human Calcitonin Hormone: A Single β-Sheet Model with a Small Hydrophobic Core

Experimentally, the human calcitonin hormone (hCT) can form highly stable amyloid protofibrils. Further, a peptide consisting of hCT residues 15-19, DFNKF, was shown to create highly ordered fibrils, similar to those formed by the entire hormone sequence. However, there are limited experimental data regarding the detailed 3D arrangement of either of these fibrils. We have modeled the DFNKF protofibril, using molecular dynamics simulations. We tested the stabilities of single sheet and of various multi sheet models. Remarkably, our most ordered and stable model consists of a parallel-stranded, single beta-sheet with a relatively insignificant hydrophobic core. We investigate the chemical and physical interactions responsible for the high structural organization of this single beta-sheet amyloid fibril. We observe that the most important chemical interactions contributing to the stability of the DFNKF organization are electrostatic, specifically between the Lys and the C terminus, between the Asp and N terminus, and a hydrogen bond network between the Asn side-chains of adjacent strands. Additionally, we observe hydrophobic and aromatic pi stacking interactions. We further simulated truncated filaments, FNKF and DFNK. Our tetra-peptide mutant simulations assume models similar to the penta-peptide. Experimentally, the FNKF does not create fibrils while DFNK does, albeit short and less ordered than DFNKF. In the simulations, the FNKF system was less stable than the DFNK and DFNKF. DFNK also lost many of its original interactions becoming less organized, however, many contacts were maintained. Thus, our results emphasize the role played by specific amino acid interactions. To further study specific interactions, we have mutated the penta-peptide, simulating DANKF, DFNKA and EFNKF. Here we describe the model, its relationship to experiment and its implications to amyloid organization.

An autonomously folding ?-hairpin derived from the human YAP65 WW domain: Attempts to define a minimum ligand-binding motif

WW domains are broadly distributed among natural proteins; these modules play a role in bringing specific proteins together. The ligands recognized by WW domains are short segments rich in proline residues. We have tried to identify the minimum substructure within a WW domain that is required for ligand binding. WW domains typically comprise ca. 40 residues and fold to a three-stranded beta-sheet. Structural data for several WW domain/ligand complexes suggest that most or all of the intermolecular contacts involve beta-strands 2 and 3. We have developed a 16-residue peptide that folds to a beta-hairpin conformation that appears to mimic beta-strands 2 and 3 of the human YAP65 WW domain, but this peptide does not bind to known ligands. Thus, the minimum binding domain is larger than the latter two strands of the WW domain beta-sheet.

Effects of alternative side chain pairings and reverse turn sequences on antiparallel sheet structure in beta-peptide hairpins

[structure: see text] We describe a series of beta-peptide hexamers that allow us to explore relationships between sequence and hairpin folding. Different reverse turn segments are compared at the central two positions, and the outer residues allow a variety of interstrand side chain-side chain pairings. NMR analysis in methanol demonstrates that several reverse turn and side chain pairing arrangements are compatible with antiparallel beta-peptide sheet structure; however, none of the beta-peptides folds in water.

The protein folding 'speed limit'

How fast can a protein possibly fold? This question has stimulated experimentalists to seek fast folding proteins and to engineer them to fold even faster. Proteins folding at or near the speed limit are prime candidates for all-atom molecular dynamics simulations. They may also have no free energy barrier, allowing the direct observation of intermediate structures on the pathways from the unfolded to the folded state. Both experimental and theoretical approaches predict a speed limit of approximately N/100micros for a generic N-residue single-domain protein, with alpha proteins folding faster than beta or alphabeta. The predicted limits suggest that most known ultrafast folding proteins can be engineered to fold more than ten times faster.

10 Residue Folded Peptide Designed by Segment Statistics

We have designed a peptide termed chignolin, consisting of only 10 amino acid residues (GYDPETGTWG), on the basis of statistics derived from more than 10,000 protein segments. The peptide folds into a unique structure in water and shows a cooperative thermal transition, both of which may be hallmarks of a protein. Also, the experimentally determined beta-hairpin structure was very close to what we had targeted. The performance of the short peptide not only implies that the methodology employed here can contribute toward development of novel techniques for protein design, but it also yields insights into the raison d'etre of an autonomous element involved in a natural protein. This is of interest for the pursuit of folding mechanisms and evolutionary processes of proteins.

Side chain interactions determine the amyloid organization: a single layer -sheet molecular structure of the calcitonin peptide segment 15–19

In this paper we present a detailed atomic model for a protofilament, the most basic organization level, of the amyloid fibre formed by the peptide DFNKF. This pentapeptide is a segment derived from the human calcitonin, a natural amyloidogenic protein. Our model, which represents the outcome of extensive explicit solvent molecular dynamics (MD) simulations of different strand/sheet organizations, is a single beta-sheet filament largely without a hydrophobic core. Nevertheless, this structure is capable of reproducing the main features of the characteristic amyloid fibril organization and provides clues to the molecular basis of its experimental aggregation behaviour. Our results show that the side chains' chemical diversity induces the formation of a complex network of interactions that finally determine the microscopic arrangement of the strands at the protofilament level. This network of interactions, consisting of both side chain-side chain and backbone-side chain interactions, confers on the final single beta-sheet arrangement an unexpected stability, both by enhancing the association of related chemical groups and, at the same time, by shielding the hydrophobic segments from the polar solvent. The chemical physical characterization of this protofilament provides hints to the possible thermodynamical basis of the supra molecular organization that allows the formation of the filaments by lateral association of the preformed protofibrils. Its regular, highly polarized structure shows how other protofilaments can assemble. In terms of structural biology, our results clearly indicate that an amyloid organization implies a degree of complexity far beyond a simple nonspecific association of peptide strands via amide hydrogen bonds.

General synthesis of unsymmetrical norbornane scaffolds as inducers for hydrogen bond interactions in peptides

Starting from readily accessible endo-cis-(2S,3R)-norbornene dicarboxylic acid benzyl monoester, a general and efficient synthetic approach toward unsymmetrical two-stranded peptidic structures was developed. In these structures the peptide strands are oriented in a parallel geometry. Their synthesis is easily applicable to a variety of amino acids and peptides. Specifically, a norbornane template as molecular scaffold induces hydrogen bonding between the adjacent peptide strands. The specific hydrogen bonding patterns between these strands were revealed by detailed NMR analysis including TOCSY/NOE experiments.

Factors involved in the stability of isolated beta-sheets: Turn sequence, beta-sheet twisting, and hydrophobic surface burial

We have recently reported on the design of a 20-residue peptide able to form a significant population of a three-stranded up-and-down antiparallel beta-sheet in aqueous solution. To improve our beta-sheet model in terms of the folded population, we have modified the sequences of the two 2-residue turns by introducing the segment DPro-Gly, a sequence shown to lead to more rigid type II' beta-turns. The analysis of several NMR parameters, NOE data, as well as Deltadelta(CalphaH), DeltadeltaC(beta), and Deltadelta(Cbeta) values, demonstrates that the new peptide forms a beta-sheet structure in aqueous solution more stable than the original one, whereas the substitution of the DPro residues by LPro leads to a random coil peptide. This agrees with previous results on beta-hairpin-forming peptides showing the essential role of the turn sequence for beta-hairpin folding. The well-defined beta-sheet motif calculated for the new designed peptide (pair-wise RMSD for backbone atoms is 0.5 +/- 0.1 A) displays a high degree of twist. This twist likely contributes to stability, as a more hydrophobic surface is buried in the twisted beta-sheet than in a flatter one. The twist observed in the up-and-down antiparallel beta-sheet motifs of most proteins is less pronounced than in our designed peptide, except for the WW domains. The additional hydrophobic surface burial provided by beta-sheet twisting relative to a "flat" beta-sheet is probably more important for structure stability in peptides and small proteins like the WW domains than in larger proteins for which there exists a significant contribution to stability arising from their extensive hydrophobic cores.

Local sequence information in cellular retinoic acid-binding protein I: Specific residue roles in β-turns*

We have recently shown that two of the beta-turns (III and IV) in the ten-stranded, beta-clam protein, cellular retinoic acid-binding protein I (CRABP I), are favored in short peptide fragments, arguing that they are encoded by local interactions (K. S. Rotondi and L. M. Gierasch, Biochemistry, 2003, Vol. 42, pp. 7976-7985). In this paper we examine these turns in greater detail to dissect the specific local interactions responsible for their observed native conformational biases. Conformations of peptides corresponding to the turn III and IV fragments were examined under conditions designed to selectively disrupt stabilizing interactions, using pH variation, chaotrope addition, or mutagenesis to probe specific side-chain influences. We find that steric constraints imposed by excluded volume effects between near neighbor residues (i,i+2), favorable polar (i,i+2) interactions, and steric permissiveness of glycines are the principal factors accounting for the observed native bias in these turns. Longer-range stabilizing interactions across the beta-turns do not appear to play a significant role in turn stability in these short peptides, in contrast to their importance in hairpins. Additionally, our data add to a growing number of examples of the 3:5 type I turn with a beta-bulge as a class of turns with high propensity to form locally defined structure. Current work is directed at the interplay between the local sequence information in the turns and more long-range influences in the mechanism of folding of this predominantly beta-sheet protein.

Stereospecific interactions of proline residues in protein structures and complexes

The constrained backbone torsion angle of a proline (Pro) residue has usually been invoked to explain its three-dimensional context in proteins. Here we show that specific interactions involving the pyrrolidine ring atoms also contribute to its location in a given secondary structure and its binding to another molecule. It is adept at participating in two rather non-conventional interactions, C-H...pi and C-H...O. The geometry of interaction between the pyrrolidine and aromatic rings, vis-à-vis the occurrence of the C-H...pi interactions has been elucidated. Some of the secondary structural elements stabilized by Pro-aromatic interactions are beta-turns, where a Pro can interact with an adjacent aromatic residue, and in antiparallel beta-sheet, where a Pro in an edge strand can interact with an aromatic residue in the adjacent strand at a non-hydrogen-bonded site. The C-H groups at the Calpha and Cdelta positions can form strong C-H...O interactions (as seen from the clustering of points) and such interactions involving a Pro residue at C' position relative to an alpha-helix can cap the hydrogen bond forming potentials of the free carbonyl groups at the helix C terminus. Functionally important Pro residues occurring at the binding site of a protein almost invariably engage aromatic residues (with one of them being held by C-H...pi interaction) from the partner molecule in the complex, and such aromatic residues are highly conserved during evolution.

Insights into the determinants of beta-sheet stability: 1H and 13C NMR conformational investigation of three-stranded antiparallel beta-sheet-forming peptides

In a previous study we designed a 20-residue peptide able to adopt a significant population of a three-stranded antiparallel beta-sheet in aqueous solution (de Alba et al. [1999]Protein Sci.8, 854-865). In order to better understand the factors contributing to beta-sheet folding and stability we designed and prepared nine variants of the parent peptide by substituting residues at selected positions in its strands. The ability of these peptides to form the target motif was assessed on the basis of NMR parameters, in particular NOE data and 13Calpha conformational shifts. The populations of the target beta-sheet motif were lower in the variants than in the parent peptide. Comparative analysis of the conformational behavior of the peptides showed that, as expected, strand residues with low intrinsic beta-sheet propensities greatly disfavor beta-sheet folding and that, as already found in other beta-sheet models, specific cross-strand side chain-side chain interactions contribute to beta-sheet stability. More interestingly, the performed analysis indicated that the destabilization effect of the unfavorable strand residues depends on their location at inner or edge strands, being larger at the latter. Moreover, in all the cases examined, favorable cross-strand side chain-side chain interactions were not strong enough to counterbalance the disfavoring effect of a poor beta-sheet-forming residue, such as Gly.

Influence of strand number on antiparallel beta-sheet stability in designed three- and four-stranded beta-sheets

We describe experiments that probe whether antiparallel beta-sheet secondary structure becomes more stable as the number of strands increases. Several groups, including ours, have explored this issue with peptides designed to adopt three-stranded beta-sheet conformations, but the conclusions have not been consistent. In this study, we examine the effect on conformational stability of beta-sheet lengthening perpendicular to the strand direction via analysis of designed peptides that adopt three-stranded or four-stranded antiparallel beta-sheet conformations in aqueous solution. The findings reported here, along with the context provided by earlier studies, suggest that antiparallel beta-sheet does, in general, become more stable when the number of strands is increased from two to three. We show that this conclusion is not influenced by the rigidity of the loop segment used to link adjacent beta-strands (D-Pro-Gly versus Asn-Gly). We show that further extension, from three strands to four, leads to a further increase in antiparallel beta-sheet stability.

Identification of conserved residue patterns in small beta-barrel proteins

Our abilities to predict three-dimensional conformation of a polypeptide, given its amino acid sequence, remain limited despite advances in structure analysis. Analysis of structures and sequences of protein families with similar secondary structural elements, but varying topologies, might help in addressing this problem. We have studied the small beta-barrel class of proteins characterized by four strands (n = 4) and a shear number of 8 (S = 8) to understand the principles of barrel formation. Multiple alignments of the various protein sequences were generated for the analysis. Positional entropy, as a measure of residue conservation, indicated conservation of non-polar residues at the core positions. The presence of a type II beta-turn among the various barrel proteins considered was another strikingly invariant feature. A conserved glycyl-aspartyl dipeptide at the beta-turn appeared to be important in guiding the protein sequence into the barrel fold. Molecular dynamics simulations of the type II beta-turn peptide suggested that aspartate is a key residue in the folding of the protein sequence into the barrel. Our study suggests that the conserved type II beta-turn and the non-polar residues in the barrel core are crucial for the folding of the protein's primary sequence into the beta-barrel conformation.

Mechanism of fast protein folding

An explosion of in vitro experimental data on the folding of proteins has revealed many examples of folding in the millisecond or faster timescale, often occurring in the absence of stable intermediate states. We review experimental methods for measuring fast protein folding kinetics, and then discuss various analytical models used to interpret these data. Finally, we classify general mechanisms that have been proposed to explain fast protein folding into two catagories, heterogeneous and homogeneous, reflecting the nature of the transition state. One heterogeneous mechanism, the diffusion-collision mechanism, can be used to interpret experimental data for a number of proteins.

Analysis of the factors that stabilize a designed two-stranded antiparallel beta-sheet

Autonomously folding beta-hairpins (two-strand antiparallel beta-sheets) have become increasingly valuable tools for probing the forces that control peptide and protein conformational preferences. We examine the effects of variations in sequence and solvent on the stability of a previously designed 12-residue peptide (1). This peptide adopts a beta-hairpin conformation containing a two-residue loop (D-Pro-Gly) and a four-residue interstrand sidechain cluster that is observed in the natural protein GB1. We show that the conformational propensity of the loop segment plays an important role in beta-hairpin stability by comparing 1 with (D)P--> N mutant 2. In addition, we show that the sidechain cluster contributes both to conformational stability and to folding cooperativity by comparing 1 with mutant 3, in which two of the four cluster residues have been changed to serine. Thermodynamic analysis suggests that the high loop-forming propensity of the (D)PG segment decreases the entropic cost of beta-hairpin formation relative to the more flexible NG segment, but that the conformational rigidity of (D)PG may prevent optimal contacts between the sidechains of the GB1-derived cluster. The enthalpic favorability of folding in these designed beta-hairpins suggests that they are excellent scaffolds for studying the fundamental mechanisms by which amino acid sidechains interact with one another in folded proteins.

Aromatic-aromatic interactions in and around alpha-helices

To understand the role of aromatic-aromatic interactions in imparting specificity to the folding process, the geometries of four aromatic residues with different sequence spacing, located in alpha-helices or five residues from helical ends, interacting with each other have been elucidated. The geometry is found to depend on the sequence difference. Specific interactions (C-H...pi and N-H...pi) which result from this geometry may cause a given pair of residues (such as Phe-His) with a particular sequence difference to occur more than expected. The most conspicuous residue in an aromatic pair in the context of helix stability is His, which is found at the last (C1) position or the two positions (Ncap and Ccap) immediately flanking the helix. An alpha-helix and a contiguous 3(10)-helix or two helices separated by a non-helical residue can have interacting aromatic pairs, the geometry of interaction and the relative orientation between the helices being rather fixed. Short helices can also have interacting residues from either side.

Domain architecture of a high mobility group A-type bacterial transcriptional factor

Myxococcus xanthus transcriptional factor CarD participates in carotenogenesis and fruiting body formation. It is the only reported prokaryotic protein having adjacent "AT-hook" DNA-binding and acidic regions characteristic of eukaryotic high mobility group A (HMGA) proteins. The latter are small, unstructured, nonhistone nuclear proteins that function as architectural factors to remodel DNA and chromatin structure and modulate various DNA binding activities. We find CarD to be predominantly dimeric with two stable domains: (a) an N-terminal domain of defined secondary and tertiary structure which is absent in eukaryotic HMGA proteins; (b) a C-terminal domain formed by the acidic and AT-hook segments and lacking defined structure. CarD, like HMGA proteins, binds specifically to the minor-groove of AT-rich DNA present in two appropriately spaced tracts. As in HMGA proteins, casein kinase II can phosphorylate the CarD acidic region, and this dramatically decreases the DNA binding affinity of CarD. The acidic region, in addition to modulating DNA binding, confers structural stability to CarD. We discuss how the structural and functional plasticity arising from domain organization in CarD could be linked to its role as a general transcriptional factor in M. xanthus.

Length-dependent stability and strand length limits in antiparallel beta -sheet secondary structure

Designed peptides that fold autonomously to specific conformations in aqueous solution are useful for elucidating protein secondary structural preferences. For example, autonomously folding model systems have been essential for establishing the relationship between alpha-helix length and alpha-helix stability, which would be impossible to probe with alpha-helices embedded in folded proteins. Here, we use designed peptides to examine the effect of strand length on antiparallel beta-sheet stability. alpha-Helices become more stable as they grow longer. Our data show that a two-stranded beta-sheet ("beta-hairpin") becomes more stable when the strands are lengthened from five to seven residues, but that further strand lengthening to nine residues does not lead to further beta-hairpin stabilization for several extension sequences examined. (In one case, all-threonine extension, there may be an additional stabilization on strand lengthening from seven to nine residues.) These results suggest that there may be an intrinsic limit to strand length for most sequences in antiparallel beta-sheet secondary structure.

Novel inactivating mutations of transforming growth factor-beta type I receptor gene in head-and-neck cancer metastases

Carcinoma cell lines are frequently refractory to transforming growth factor-beta (TGF beta)-mediated cell cycle arrest. Whether and how TGF beta signaling is disrupted in the majority of human tumors, however, remains unclear. To investigate whether TGF beta signaling might be disrupted by inactivation of the key signaling molecule, the TGF beta type I (T beta R-I) receptor, and whether or not T beta R-I inactivation is associated with late stage disease, we conducted a comprehensive structural analysis of the T beta R-I gene in fine-needle aspirates of 23 head-&-neck cancer metastases. We encountered 4 different mutations of T beta R-I, 3 of which have not been previously identified. In 1 case, we found a somatic intragenic 4-bp deletion predicting for a truncation of the receptor protein. This is the first example of a true loss-of-function mutation of T beta R-I in a human epithelial neoplasm. In 2 other cases, we identified missense mutations located between the juxtamembrane- and serine-threonine kinase domains. One of these resulted in an alanine-to-threonine substitution (A230T), which disrupts receptor signaling activity by causing rapid protein degradation within the endoplasmatic reticulum. This represents a novel mechanism of inactivation of a TGF beta signaling intermediate. Finally, we identified a serine-to-tyrosine substitution at codon 387 (S387Y) in a metastasis but not in the corresponding primary tumor. We had previously shown this S387Y mutant to be predominantly associated with breast cancer metastases and to have a diminished ability to mediate TGF beta-dependent signaling. In aggregate, these findings provide further support for the hypothesis that inactivation of the TGF beta signaling pathway occurs in a significant subset of human cancers.

The interrelationships of side-chain and main-chain conformations in proteins

The accurate determination of a large number of protein structures by X-ray crystallography makes it possible to conduct a reliable statistical analysis of the distribution of the main-chain and side-chain conformational angles, how these are dependent on residue type, adjacent residue in the sequence, secondary structure, residue-residue interactions and location at the polypeptide chain termini. The interrelationship between the main-chain (phi, psi) and side-chain (chi 1) torsion angles leads to a classification of amino acid residues that simplify the folding alphabet considerably and can be a guide to the design of new proteins or mutational studies. Analyses of residues occurring with disallowed main-chain conformation or with multiple conformations shed some light on why some residues are less favoured in thermophiles.