Protein refolding versus aggregation: computer simulations on an intermediate-resolution protein model.

Computer simulations are performed on a system of eight model peptide chains to study how the competition between protein refolding and aggregation affects the optimal conditions for refolding of four-helix bundles. The discontinuous molecular dynamics algorithm is utilized along with an intermediate-resolution protein model that we developed for this work. Physically, the model is much more detailed than any model used to date for simulations of protein aggregation. Each model residue consists of a detailed, three-bead backbone and a simplified, single-bead side-chain. Excluded volume, hydrogen bond, and hydrophobic interactions are modeled with discontinuous (i.e. hard-sphere and square-well) potentials. Simulations efficiently sample conformational space, and complete folding trajectories from random initial configurations to two four-helix bundles are possible within two days on a single processor workstation. Folding of the bundles follows two main pathways, one through a trimeric intermediate and the other through an intermediate with two dimers. The proportion of trajectories that follow each route is significantly different for the eight-peptide system in this work than in a previously studied four-peptide system, which yields one four-helix bundle, suggesting, as our previous simulations have, that protein folding properties are strongly influenced by the presence of other proteins. Folding of the bundles is optimal within a fixed temperature range, with the high-temperature boundary a function of the complexity of the protein (or oligomer) to be folded and the low-temperature boundary a function of the complexity of the protein's environment. Above the optimal temperature range for folding, the model chains tend to unfold; below the optimal range, the model chains tend to aggregate. As has been seen previously, aggregates have substantial levels of native secondary structure, suggesting that aggregates are composed largely of partially folded intermediates, not denatured chains.     

Construction of an intermediate-resolution lattice model and re-examination of the helix-coil transition: a dynamic Monte Carlo simulation

In protein modeling, spatial resolution and computational efficiency are always incompatible. As a compromise, an intermediate-resolution lattice model has been constructed in the present work. Each residue is decomposed into four basic units, i.e. the α-carbon group, the carboxyl group, the imino group, and the side-chain group, and each basic coarse-grained unit is represented by a minimum cubic box with eight lattice sites. The spacing of the lattice is about 0.56?Å, holding the highest spatial resolution for the present lattice protein models. As the first report of this new model, the helix-coil transition of a polyalanine chain was examined via dynamic Monte Carlo simulation. The period of formed α-helix was about 3.68 residues, close to that of a natural α-helix. The resultant backbone motion was found to be in the realistic regions of the conformational space in the Ramachandran plot. Helix propagation constant and nucleation constant were further determined through the dynamic hydrogen bonding process and torsional angle variation, and the results were used to make comparison between classical Zimm-Bragg theory and Lifson-Roig theory based on the Qian-Schellman relationship. The simulation results confirmed that our lattice model can reproduce the helix-coil transition of polypeptide and construct a moderately fine α-helix conformation without significantly weakening the priority in efficiency for a lattice model.     

Competing interactions contributing to alpha-helical stability in aqueous solution.

The stability of a 15-residue peptide has been investigated using CD spectroscopy and molecular simulation techniques. The sequence of the peptide was designed to include key features that are known to stabilize alpha-helices, including ion pairs, helix dipole capping, peptide bond capping, and aromatic interactions. The degree of helicity has been determined experimentally by CD in three solvents (aqueous buffer, methanol, and trifluoroethanol) and at two temperatures. Simulations of the peptide in the aqueous system have been performed over 500 ps at the same two temperatures using a fully explicit solvent model. Consistent with the CD data, the degree of helicity is decreased at the higher temperature. Our analysis of the simulation results has focused on competition between different side-chain/side-chain and side-chain/main-chain interactions, which can, in principle, stabilize the helix. The unfolding in aqueous solution occurs at the amino terminus because the side-chain interactions are insufficient to stabilize both the helix dipole and the peptide hydrogen bonds. Loss of capping of the peptide backbone leads to water insertion within the first peptide hydrogen bond and hence unfolding. In contrast, the carboxy terminus of the alpha-helix is stable in both simulations because the C-terminal lysine residue stabilizes the helix dipole, but at the expense of an ion pair.     

Backbone-driven collapse in unfolded protein chains

Collapse of unfolded protein chains is an early event in folding. It affects structural properties of intrinsically disordered proteins, which take a considerable fraction of the human proteome. Collapse is generally believed to be driven by hydrophobic forces imposed by the presence of nonpolar amino acid side chains. Contributions from backbone hydrogen bonds to protein folding and stability, however, are controversial. To date, the experimental dissection of side-chain and backbone contributions has not yet been achieved because both types of interactions are integral parts of protein structure. Here, we realized this goal by applying mutagenesis and chemical modification on a set of disordered peptides and proteins. We measured the protein dimensions and kinetics of intra-chain diffusion of modified polypeptides at the level of individual molecules using fluorescence correlation spectroscopy, thereby avoiding artifacts commonly caused by aggregation of unfolded protein material in bulk. We found no contributions from side chains to collapse but, instead, identified backbone interactions as a source sufficient to form globules of native-like dimensions. The presence of backbone hydrogen bonds decreased polypeptide water solubility dramatically and accelerated the nanosecond kinetics of loop closure, in agreement with recent predictions from computer simulation. The presence of side chains, instead, slowed loop closure and modulated the dimensions of intrinsically disordered domains. It appeared that the transient formation of backbone interactions facilitates the diffusive search for productive conformations at the early stage of folding and within intrinsically disordered proteins.     

Long-range intra-protein communication can be transmitted by correlated side-chain fluctuations alone

Allosteric regulation is a key component of cellular communication, but the way in which information is passed from one site to another within a folded protein is not often clear. While backbone motions have long been considered essential for long-range information conveyance, side-chain motions have rarely been considered. In this work, we demonstrate their potential utility using Monte Carlo sampling of side-chain torsional angles on a fixed backbone to quantify correlations amongst side-chain inter-rotameric motions. Results indicate that long-range correlations of side-chain fluctuations can arise independently from several different types of interactions: steric repulsions, implicit solvent interactions, or hydrogen bonding and salt-bridge interactions. These robust correlations persist across the entire protein (up to 60 Å in the case of calmodulin) and can propagate long-range changes in side-chain variability in response to single residue perturbations.     

Mechanism of selection of side-chain rotamers in α-helices

In this study, a possible mechanism of selection of side-chain rotamers based on the rotamer distributions in known coiled-coil proteins is suggested. According to this mechanism, interhelical hydrophobic, polar, and packing interactions bring alpha-helices closer to each other and this effect squeezes side chains out of the helix-helix interface. As a result, in dimeric coiled coils and long alpha-alpha-hairpins where alpha-helices are packed in a face-to-face manner, most side chains occupying the a-positions have t-rotamers and those in the d-positions g(-)-rotamers. In tetramers, where alpha-helices are packed side-by-side, most side chains in the a-positions adopt g(-)-rotamers and those in the d-positions t-rotamers.     

Solvent effects on the conformational transition of a model polyalanine peptide

We have investigated the folding of polyalanine by combining discontinuous molecular dynamics simulation with our newly developed off-lattice intermediate-resolution protein model. The thermodynamics of a system containing a single Ac-KA(14)K-NH(2) molecule has been explored by using the replica exchange simulation method to map out the conformational transitions as a function of temperature. We have also explored the influence of solvent type on the folding process by varying the relative strength of the side-chain's hydrophobic interactions and backbone hydrogen bonding interactions. The peptide in our simulations tends to mimic real polyalanine in that it can exist in three distinct structural states: alpha-helix, beta-structures (including beta-hairpin and beta-sheet-like structures), and random coil, depending upon the solvent conditions. At low values of the hydrophobic interaction strength between nonpolar side-chains, the polyalanine peptide undergoes a relatively sharp transition between an alpha-helical conformation at low temperatures and a random-coil conformation at high temperatures. As the hydrophobic interaction strength increases, this transition shifts to higher temperatures. Increasing the hydrophobic interaction strength even further induces a second transition to a beta-hairpin, resulting in an alpha-helical conformation at low temperatures, a beta-hairpin at intermediate temperatures, and a random coil at high temperatures. At very high values of the hydrophobic interaction strength, polyalanines become beta-hairpins and beta-sheet-like structures at low temperatures and random coils at high temperatures. This study of the folding of a single polyalanine-based peptide sets the stage for a study of polyalanine aggregation in a forthcoming paper.     

Side-chain rotamers in protein alpha-alpha-hairpins and a mechanism of their selection

The observed features of side-chain rotamer distributions in protein alpha-alpha-hairpins are described. It was found that in left-turned alpha-alpha-hairpins most side chains occupying d-positions have t-rotamers and those in g-positions g- -rotamers. In right-turned alpha-alpha-hairpins, most side chains in a-positions adopt g- -rotamers and those in e-positions t-rotamers. Analysis of these features enables us to conclude that selection of side-chain rotamers in alpha-alpha-hairpins depends on both the type of the alpha-helix packing and the residue position. The observed features can be explained taking into account the squeezing mechanism according to which interhelical interactions bring alpha-helices closer to each other and this effect squeezes side chains out of the helix-helix interface and as a result they adopt unique conformations.     

Three-dimensional structure of proteolytic fragment 163-231 of bacterioopsin determined from nuclear magnetic resonance data in solution.

546 NOESY cross-peak volumes were measured in the two-dimensional NOESY spectrum of proteolytic fragment 163-231 of bacterioopsin in organic solution. These data and 42 detected hydrogen bonds were applied for determining the peptide spatial structure. The fold of the polypeptide chain was determined by local structure analysis, a distance geometry approach and systematic search for energetically allowed side-chain rotamers which are consistent with experimental NOESY cross-peak volumes. The effective rotational correlation time of 6 ns for the molecule was evaluated from optimization of the local structure to meet NOE data and from the dependence on mixing time of the NiH/Ci alpha H cross-peak volumes of the residues in alpha-helical conformation. The resulting structure has two well defined alpha-helical regions, 168-191 and 198-227, with root-mean-square deviation 44 pm and 69 pm, respectively, between the backbone atoms in 14 final energy refined conformations. The alpha-helices correspond to transmembrane segments F and G of bacteriorhodopsin. The segment F contains proline 186, which introduces a kink of about 25 degrees with a disruption of the hydrogen bond with the NH group of the following residue. The segments are connected by a flexible loop region 192-197. Torsion angles chi 1 are unequivocally defined for 62% of side chains in the alpha-helices but half of them differ from electron cryo-microscopy (ECM) model of bacteriorhodopsin, apparently because of the low resolution of ECM. Nevertheless, the F and G segments can be packed as in the ECM model and with side-chain conformations consistent with all NMR data in solution.     

Molecular dissection,tissue localization and Ca2+ binding of the ryanodine receptor of Caenorhabditis elegans

We investigated the possible role of residues at the Ccap position in an alpha-helix on protein stability. A set of 431 protein alpha-helices containing a C'-Gly from the Protein Data Bank (PDB) was analyzed, and the normalized frequencies for finding particular residues at the Ccap position, the average fraction of buried surface area, and the hydrogen bonding patterns of the Ccap residue side-chain were calculated. We found that on average the Ccap position is 70% buried and noted a significant correlation (R=0.8) between the relative burial of this residue and its hydrophobicity as defined by the Gibbs energy of transfer from octanol or cyclohexane to water. Ccap residues with polar side-chains are commonly involved in hydrogen bonding. The hydrogen bonding pattern is such that, the longer side-chains of Glu, Gln, Arg, Lys, His form hydrogen bonds with residues distal (>+/-4) in sequence, while the shorter side-chains of Asp, Asn, Ser, Thr exhibit hydrogen bonds with residues close in sequence (<+/-4), mainly involving backbone atoms. Experimentally we determined the thermodynamic propensities of residues at the Ccap position using the protein ubiquitin as a model system. We observed a large variation in the stability of the ubiquitin variants depending on the nature of the Ccap residue. Furthermore, the measured changes in stability of the ubiquitin variants correlate with the hydrophobicity of the Ccap residue. The experimental results, together with the statistical analysis of protein structures from the PDB, indicate that the key hydrophobic capping interactions between a helical residue (C3 or C4) and a residue outside the helix (C", C3' or C4') are frequently enhanced by the hydrophobic interactions with Ccap residues.     

Exhaustive Metropolis Monte Carlo sampling and analysis of polyalanine conformations adopted under the influence of hydrogen bonds

We propose a novel Metropolis Monte Carlo procedure for protein modeling and analyze the influence of hydrogen bonding on the distribution of polyalanine conformations. We use an atomistic model of the polyalanine chain with rigid and planar polypeptide bonds, and elastic alpha carbon valence geometry. We adopt a simplified energy function in which only hard-sphere repulsion and hydrogen bonding interactions between the atoms are considered. Our Metropolis Monte Carlo procedure utilizes local crankshaft moves and is combined with parallel tempering to exhaustively sample the conformations of 16-mer polyalanine. We confirm that Flory's isolated-pair hypothesis (the steric independence between the dihedral angles of individual amino acids) does not hold true in long polypeptide chains. In addition to 3(10)- and alpha-helices, we identify a kink stabilized by 2 hydrogen bonds with a shared acceptor as a common structural motif. Varying the strength of hydrogen bonds, we induce the helix-coil transition in the model polypeptide chain. We compare the propensities for various hydrogen bonding patterns and determine the degree of cooperativity of hydrogen bond formation in terms of the Hill coefficient. The observed helix-coil transition is also quantified according to Zimm-Bragg theory.     

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).     

Structures of N-termini of helices in proteins.

We have surveyed 393 N-termini of alpha-helices and 156 N-termini of 3(10)-helices in 85 high resolution, non-homologous protein crystal structures for N-cap side-chain rotamer preferences, hydrogen bonding patterns, and solvent accessibilities. We find very strong rotamer preferences that are unique to N-cap sites. The following rules are generally observed for N-capping in alpha-helices: Thr and Ser N-cap side chains adopt the gauche - rotamer, hydrogen bond to the N3 NH and have psi restricted to 164 +/- 8 degrees. Asp and Asn N-cap side chains either adopt the gauche - rotamer and hydrogen bond to the N3 NH with psi = 172 +/- 10 degrees, or adopt the trans rotamer and hydrogen bond to both the N2 and N3 NH groups with psi = 1-7 +/- 19 degrees. With all other N-caps, the side chain is found in the gauche + rotamer so that the side chain does not interact unfavorably with the N-terminus by blocking solvation and psi is unrestricted. An i, i + 3 hydrogen bond from N3 NH to the N-cap backbone C = O in more likely to form at the N-terminus when an unfavorable N-cap is present. In the 3(10)-helix Asn and Asp remain favorable N-caps as they can hydrogen bond to the N2 NH while in the trans rotamer; in contrast, Ser and Thr are disfavored as their preferred hydrogen bonding partner (N3 NH) is inaccessible. This suggests that Ser is the optimum choice of N-cap when alpha-helix formation is to be encouraged while 3(10)-helix formation discouraged. The strong energetic and structural preferences found for N-caps, which differ greatly from positions within helix interiors, suggest that N-caps should be treated explicitly in any consideration of helical structure in peptides or proteins.     

Deterministic features of side-chain main-chain hydrogen bonds in globular protein structures

A total of 19 835 polar residues from a data set of 250 non-homologous and highly resolved protein crystal structures were used to identify side-chain main-chain (SC-MC) hydrogen bonds. The ratio of the number of SC-MC hydrogen bonds to the total number of polar residues is close to 1:2, indicating the ubiquitous nature of such hydrogen bonds. Close to 56% of the SC-MC hydrogen bonds are local involving side-chain acceptor/donor ('i') and a main-chain donor/acceptor within the window i-5 to i+5. These short-range hydrogen bonds form well defined conformational motifs characterized by specific combinations of backbone and side-chain torsion angles. (a) The Ser/Thr residues show the greatest preference in forming intra-helical hydrogen bonds between the atoms O(gamma)(i) and O(i-4). More than half the examples of such hydrogen bonds are found at the middle of alpha-helices rather than at their ends. The most favoured motif of these examples is alpha(R)alpha(R)alpha(R)alpha(R)(g(-)). (b) These residues also show great preference to form hydrogen bonds between O(gamma)(i) and O(i-3), which are closely related to the previous type and though intra-helical, these hydrogen bonds are more often found at the C-termini of helices than at the middle. The motif represented by alpha(R)alpha(R)alpha(R)alpha(R)(g(+)) is most preferred in these cases. (c) The Ser, Thr and Glu are the most frequently found residues participating in intra-residue hydrogen bonds (between the side-chain and main-chain of the same residue) which are characterized by specific motifs of the form beta(g(+)) for Ser/Thr residues and alpha(R)(g(-)g(+)t) for Glu/Gln. (d) The side-chain acceptor atoms of Asn/Asp and Ser/Thr residues show high preference to form hydrogen bonds with acceptors two residues ahead in the chain, which are characterized by the motifs beta (tt')alphaR and beta(t)alpha(R), respectively. These hydrogen bonded segments, referred to as Asx turns, are known to provide stability to type I and type I' beta-turns. (e) Ser/Thr residues often form a combination of SC-MC hydrogen bonds, with the side-chain donor hydrogen bonded to the carbonyl oxygen of its own peptide backbone and the side-chain acceptor hydrogen bonded to an amide hydrogen three residues ahead in the sequence. Such motifs are quite often seen at the beginning of alpha-helices, which are characterized by the beta(g(+))alpha(R)alpha(R) motif. A remarkable majority of all these hydrogen bonds are buried from the protein surface, away from the surrounding solvent. This strongly indicates the possibility of side-chains playing the role of the backbone, in the protein interiors, to satisfy the potential hydrogen bonding sites and maintaining the network of hydrogen bonds which is crucial to the structure of the protein.     

Energetics and cooperativity of the hydrogen bonding and anchor interactions that bind peptides to MHC class II protein

The complexity of the interaction between major histocompatibility complex class II (MHC II) proteins and peptide ligands has been revealed through structural studies and crystallographic characterization. Peptides bind through side-chain "anchor" interactions with MHC II pockets and an extensive array of genetically conserved hydrogen bonds to the peptide backbone. Here we quantitatively investigate the kinetic hierarchy of these interactions. We present results detailing the impact of single side-chain mutations of peptide anchor residues on dissociation rates, utilizing two I-A(d)-restricted peptides, one of which has a known crystal structure, and 24 natural and non-natural amino acid mutant variants of these peptides. We find that the N-terminal P1, P4 and P6 anchor-pocket interactions can make significant contributions to binding stability. We also investigate the interactions of these peptides with four I-A(d) MHC II proteins, each mutated to disrupt conserved hydrogen bonds to the peptide backbone. These complexes exhibit kinetic behavior suggesting that binding energy is disproportionately invested near the peptide N terminus for backbone hydrogen bonds. We then evaluate the effects of simultaneously modifying both anchor and hydrogen bonding interactions. A quantitative analysis of 71 double mutant cycles reveals that there is little apparent cooperativity between anchor residue interactions and hydrogen bonds, even when they are directly adjacent (<5A).     

Side-chain entropy effects on protein secondary structure formation

Loss of conformational entropy is one of the primary factors opposing protein folding. Both the backbone and side-chain of each residue in a protein will have their freedom of motion restricted in the final folded structure. The type of secondary structure of which a residue is part will have a significant impact on how much side-chain entropy is lost. Side-chain conformational entropies have previously been determined for folded proteins, simple models of unfolded proteins, alpha-helices, and a dipeptide model for beta-strands, but not for polyproline II (PII) helices. In this work, we present side-chain conformational estimates for the three regular secondary structure types: alpha-helices, beta-strands, and PII helices. Entropies are estimated from Monte Carlo computer simulations. Beta-strands are modeled as two structures, parallel and antiparallel beta-strands. Our data indicate that restraining a residue to the PII helix or antiparallel beta-strand conformations results in side-chain entropies equal to or higher than those obtained by restraining residues to the parallel beta-strand conformation. Side-chains in the alpha-helix conformation have the lowest side-chain entropies. The observation that extended structures retain the most side-chain entropy suggests that such structures would be entropically favored in unfolded proteins under folding conditions. Our data indicate that the PII helix conformation would be somewhat favored over beta-strand conformations, with antiparallel beta-strand favored over parallel. Notably, our data imply that, under some circumstances, residues may gain side-chain entropy upon folding. Implications of our findings for protein folding and unfolded states are discussed.     

The effect of conformation on the CD of interacting helices: a theoretical study of tropomyosin

A recent report [M. E. Holtzer, et al. (1988) Biophysics Journal, 53, 96a] of the anomalous CD spectrum of the tropomyosin (TM) fragment 11TM127 motivated us to model the system as two 21-residue alpha-helices distorted to a coiled-coil conformation. We used strong-coupling exciton theory to model the optical properties of the system. Two backbone amide excited states (n pi* and pi pi*) were considered, as well as four excited states (Lb, La, Bb, Ba) for the phenolic side chain. We calculated the effect of superhelix formation on the backbone CD spectrum. The decrease in molar ellipticity of the alpha-helix parallel-polarized transition at 208 nm was found to be a simple function of superhelix tilt angle. We then modeled a coiled coil (radius = 5.5 A, pitch = -140 A) with one aromatic ring per superhelix. Steric interactions between aromatic side chains in a coiled coil were calculated as a function of side-chain conformation and heptet position. Steric interactions between phenolic rings will be significant for heptet positions a and d, but not for positions b, c, e, f, or g. We calculated the phenolic Lb transition rotational strength as a function of position within the heptet repeats, and of all possible side-chain dihedral angles, chi 1 and chi 2. When tyrosines were placed at heptet positions b, c, e, f, or g, the rotational-strength surface was nearly identical to that of a single tyrosine in an undistorted helix. In contrast, the rotational-strength surface for tyrosines in heptet positions a or d showed substantial intertyrosine coupling components. The rotational-strength surfaces for the three types of heptet positions (position a, position d, and the others) allowed an interpretation of the aromatic CD spectra of TM and its fragments. It was predicted that the three types of heptet positions will be spectroscopically distinguishable.     

Side chain burial and hydrophobic core packing in protein folding transition states

A critical step in the folding pathway of globular proteins is the formation of a tightly packed hydrophobic core. Several mutational studies have addressed the question of whether tight packing interactions are present during the rate-limiting step of folding. In some of these investigations, substituted side chains have been assumed to form native-like interactions in the transition state when the folding rates of mutant proteins correlate with their native-state stabilities. Alternatively, it has been argued that side chains participate in nonspecific hydrophobic collapse when the folding rates of mutant proteins correlate with side-chain hydrophobicity. In a reanalysis of published data, we have found that folding rates often correlate similarly well, or poorly, with both native-state stability and side-chain hydrophobicity, and it is therefore not possible to select an appropriate transition state model based on these one-parameter correlations. We show that this ambiguity can be resolved using a two-parameter model in which side chain burial and the formation of all other native-like interactions can occur asynchronously. Notably, the model agrees well with experimental data, even for positions where the one-parameter correlations are poor. We find that many side chains experience a previously unrecognized type of transition state environment in which specific, native-like interactions are formed, but hydrophobic burial dominates. Implications of these results to the design and analysis of protein folding studies are discussed.     

The occurrence of C--H...O hydrogen bonds in alpha-helices and helix termini in globular proteins

A comprehensive database analysis of C--H...O hydrogen bonds in 3124 alpha-helices and their corresponding helix termini has been carried out from a nonredundant data set of high-resolution globular protein structures resolved at better than 2.0 A in order to investigate their role in the helix, the important protein secondary structural element. The possible occurrence of 5 --> 1 C--H...O hydrogen bond between the ith residue CH group and (i - 4)th residue C==O with C...O < or = 3.8 A is studied, considering as potential donors the main-chain Calpha and the side-chain carbon atoms Cbeta, Cgamma, Cdelta and Cepsilon. Similar analysis has been carried out for 4 --> 1 C--H...O hydrogen bonds, since the C--H...O hydrogen bonds found in helices are predominantly of type 5 --> 1 or 4 --> 1. A total of 17,367 (9310 of type 5 --> 1 and 8057 of type 4 --> 1) C--H...O hydrogen bonds are found to satisfy the selected criteria. The average stereochemical parameters for the data set suggest that the observed C--H...O hydrogen bonds are attractive interactions. Our analysis reveals that the Cgamma and Cbeta hydrogen atom(s) are frequently involved in such hydrogen bonds. A marked preference is noticed for aliphatic beta-branched residue Ile to participate in 5 --> 1 C--H...O hydrogen bonds involving methylene Cgamma 1 atom as donor in alpha-helices. This may be an enthalpic compensation for the greater loss of side-chain conformational entropy for beta-branched amino acids due to the constraint on side-chain torsion angle, namely, chi1, when they occur in helices. The preference of amino acids for 4 --> 1 C--H...O hydrogen bonds is found to be more for Asp, Cys, and for aromatic residues Trp, Phe, and His. Interestingly, overall propensity for C--H...O hydrogen bonds shows that a majority of the helix favoring residues such as Met, Glu, Arg, Lys, Leu, and Gln, which also have large side-chains, prefer to be involved in such types of weak attractive interactions in helices. The amino acid side-chains that participate in C--H...O interactions are found to shield the acceptor carbonyl oxygen atom from the solvent. In addition, C--H...O hydrogen bonds are present along with helix stabilizing salt bridges. A novel helix terminating interaction motif, X-Gly with Gly at C(cap) position having 5 --> 1 Calpha--H...O, and a chain reversal structural motif having 1 --> 5 Calpha-H...O have been identified and discussed. Our analysis highlights that a multitude of local C--H...O hydrogen bonds formed by a variety of amino acid side-chains and Calpha hydrogen atoms occur in helices and more so at the helix termini. It may be surmised that the main-chain Calpha and the side-chain CH that participate in C--H...O hydrogen bonds collectively augment the cohesive energy and thereby contribute together with the classical N--H...O hydrogen bonds and other interactions to the overall stability of helix and therefore of proteins.     

A calorimetric study of the folding-unfolding of an alpha-helix with covalently closed N and C-terminal loops.

The thermal melting of a dicyclic 29-residue peptide, having helix-stabilizing side-chain to side-chain covalent links at each terminal, has been studied by circular dichroism spectropolarimetry (CD) and differential scanning calorimetry (DSC). The CD spectra for this dicyclic peptide indicate that it is monomeric, almost fully alpha-helical at -10 degrees C, and undergoes a reversible transition from the folded to the disordered state with increasing temperature. The temperature dependencies of the ellipticity at 222 nm and the excess heat capacity measured calorimetrically are well fit by a two-state model, which indicates a cooperative melting transition that is complete within the temperature ranges of these experiments (from -10 degrees C to 100 degrees C). This allows a complete analysis of the thermodynamics of helix formation. The helix unfolding is found to proceed with a small positive heat-capacity increment, consistent with the solvation of some non-polar groups upon helix unfolding. It follows that the hydrogen bonds are not the only factors responsible for the formation of the alpha-helix, and that hydrophobic interactions are also playing a role in its stabilization. At 30 degrees C, the calorimetric enthalpy and entropy values are estimated to be 650(+/-50) cal mol(-1)and 2.0(+/-0.2) cal K(-1)mole(-1), respectively, per residue of this peptide. Comparison with the thermodynamic characteristics obtained for the unfolding of double-stranded alpha-helical coiled-coils shows that at that temperature the enthalpic contribution of non-polar groups to the stabilization of the alpha-helix is insignificant and the estimated transition enthalpy can be assigned to the hydrogen bonds. With increasing temperature, the increasing magnitude of the negative enthalpy of hydration of the exposed polar groups should decrease the helix-stabilizing enthalpy of the backbone hydrogen bonds. However, the helix-stabilizing negative entropy of hydration of these groups should also increase in magnitude with increasing temperature, offsetting this effect.