Serine and threonine residues bend alpha-helices in the chi(1) = g(-) conformation

The relationship between the Ser, Thr, and Cys side-chain conformation (chi(1) = g(-), t, g(+)) and the main-chain conformation (phi and psi angles) has been studied in a selection of protein structures that contain alpha-helices. The statistical results show that the g(-) conformation of both Ser and Thr residues decreases their phi angles and increases their psi angles relative to Ala, used as a control. The additional hydrogen bond formed between the O(gamma) atom of Ser and Thr and the i-3 or i-4 peptide carbonyl oxygen induces or stabilizes a bending angle in the helix 3-4 degrees larger than for Ala. This is of particular significance for membrane proteins. Incorporation of this small bending angle in the transmembrane alpha-helix at one side of the cell membrane results in a significant displacement of the residues located at the other side of the membrane. We hypothesize that local alterations of the rotamer configurations of these Ser and Thr residues may result in significant conformational changes across transmembrane helices, and thus participate in the molecular mechanisms underlying transmembrane signaling. This finding has provided the structural basis to understand the experimentally observed influence of Ser residues on the conformational equilibrium between inactive and active states of the receptor, in the neurotransmitter subfamily of G protein-coupled receptors.     

Sequence and structure patterns in proteins from an analysis of the shortest helices: implications for helix nucleation

The shortest helices (three-length 3(10) and four-length alpha), most abundant among helices of different lengths, have been analyzed from a database of protein structures. A characteristic feature of three-length 3(10)-helices is the shifted backbone conformation for the C-terminal residue (phi,psi angles: -95 degrees,0 degrees ), compared to the rest of the helix (-62 degrees,-24 degrees ). The deviation can be attributed to the release of electrostatic repulsion between the carbonyl oxygen atoms at the two C-terminal residues and further stabilization (due to a more linear geometry) of an intrahelical hydrogen bond. A consequence of this non-canonical C-terminal backbone conformation can be a potential origin of helix kinks when a 3(10)-helix is sequence-contiguous at the alpha-helix N-terminal. An analysis of hydrogen bonding, as well as hydrophobic interactions in the shortest helices shows that capping interactions, some of them not observed for longer helices, dominate at the N termini. Further, consideration of the distribution of amino acid residues indicates that the shortest helices resemble the N-terminal end of alpha-helices rather than the C terminus, implying that the folding of helices may be initiated at the N-terminal end, which does not get propagated in the case of the shortest helices. Finally, pairwise comparison of beta-turns and the shortest helices, based on correlation matrices of site-specific amino acid composition, and the relative abundance of these short secondary structural elements, leads to a helix nucleation scheme that considers the formation of an isolated beta-turn (and not an alpha-turn) as the helix nucleation step, with shortest 3(10)-helices as intermediates between the shortest alpha-helix and the beta-turn. Our results ascribe an important role played by shortest 3(10)-helices in proteins with important structural and folding implications.     

Proline kinks in transmembrane alpha-helices

Integral membrane proteins often contain proline residues in their presumably alpha-helical transmembrane segments. This is in marked contrast to globular proteins, where proline is rarely found inside alpha-helices. Proline residues cause kinks in helices, and, in addition to leaving the i-4 backbone carbonyl without its normal hydrogen bond donor, also sterically prevent the (i-3)-carbonyl-(i + l)-amide backbone hydrogen bond from forming. Here, some structural aspects of proline kinks in transmembrane helices are discussed on the basis of an analysis of Pro-kinked helices in the photosynthetic reaction center and bacteriorhodopsin, as well as results from an analysis of Pro-containing transmembrane segments identified in the NBRF Protein Sequence Databank.     

Structural and functional characterization of pi bulges and other short intrahelical deformations

We data-mined the Protein Data Bank for short intrahelical deformations, including pi bulges. These are defined as a contiguous stretch of intrahelical residues deviating from the standard alpha-helical i-->i-4 hydrogen bonding pattern, bilaterally flanked by at least one alpha-helical turn resulting in a helix kink of less than 40 degrees. We find that such motifs exist in 4.7% of a PDB subset filtered by quality metrics (resolution <2.5 A, R-factor <0.25, sequence identity <35%). These are typically characterized by at least one i-->i-5 main chain hydrogen bond, with energetically favorable main chain dihedral angles, followed by a variable number of main chain carbonyl groups that do not accept intrahelical main chain hydrogen bonds. Their stabilization commonly occurs via hydrogen bonding to water molecules or polar groups. Numerous deformations are implicated in basic yet vital functional roles, commonly as ligand binding site contributors.     

A survey of left-handed polyproline II helices

Left-handed polyproline II helices (PPII) are contiguous elements of protein secondary structure in which the phi and psi angles of constituent residues are restricted to around -75 degrees and 145 degrees, respectively. They are important in structural proteins, in unfolded states and as ligands for signaling proteins. Here, we present a survey of 274 nonhomologous polypeptide chains from proteins of known structure for regions that form these structures. Such regions are rare, but the majority of proteins contain at least one PPII helix. Most PPII helices are shorter than five residues, although the longest found contained 12 amino acids. Proline predominates in PPII, but Gln and positively charged residues are also favored. The basis of Gln's prevalence is its ability to form an i, i + 1 side-chain to main-chain hydrogen bond with the backbone carbonyl oxygen of the proceeding residue; this helps to fix the psi angle of the Gln and the phi and psi of the proceeding residue in PPII conformations and explains why Gln is favored at the first position in a PPII helix. PPII helices are highly solvent exposed, which explains why apolar amino acids are disfavored despite preferring this region of phi/psi space when in isolation. PPII helices have perfect threefold rotational symmetry and within these structures we find significant correlation between the hydrophobicity of residues at i and i + 3; thus, PPII helices in globular proteins can be considered to be amphipathic.     

Structures, basins, and energies: a deconstruction of the Protein Coil Library

Globular proteins adopt complex folds, composed of organized assemblies of alpha-helix and beta-sheet together with irregular regions that interconnect these scaffold elements. Here, we seek to parse the irregular regions into their structural constituents and to rationalize their formative energetics. Toward this end, we dissected the Protein Coil Library, a structural database of protein segments that are neither alpha-helix nor beta-strand, extracted from high-resolution protein structures. The backbone dihedral angles of residues from coil library segments are distributed indiscriminately across the phi,psi map, but when contoured, seven distinct basins emerge clearly. The structures and energetics associated with the two least-studied basins are the primary focus of this article. Specifically, the structural motifs associated with these basins were characterized in detail and then assessed in simple simulations designed to capture their energetic determinants. It is found that conformational constraints imposed by excluded volume and hydrogen bonding are sufficient to reproduce the observed ,psi distributions of these motifs; no additional energy terms are required. These three motifs in conjunction with alpha-helices, strands of beta-sheet, canonical beta-turns, and polyproline II conformers comprise approximately 90% of all protein structure.     

Localization of bound water in the solution structure of the immunoglobulin binding domain of streptococcal protein G. Evidence for solvent-induced helical distortion in solution.

The presence of bound water in the solution structure of the IgG binding domain of streptococcal protein G has been investigated by nuclear magnetic resonance using three-dimensional 1H rotating frame Overhauser 1H-15N multiple quantum coherence spectroscopy. The backbone amide protons of three residues, Ala20, Gln32 and Tyr33, are found to be in close proximity to bound water. Examination of the three-dimensional structure of the IgG binding domain indicates that in the vicinity of these three residues there are no backbone groups that do not already participate in hydrogen bonding and there are no suitably placed side-chain groups available for hydrogen bonding with water. As the lifetime of the bound water detected in this nuclear magnetic resonance experiment is greater than about one nanosecond, it is likely that the two bound water molecules participate in a bifurcating hydrogen bonding network comprising a CO-NH hydrogen bonded pair, such that the water molecule accepts a hydrogen bond from the NH proton and donates one to the carbonyl oxygen with the result that the amide proton is involved in a three center hydrogen bond. On the basis of the structure, one water molecule participates in such an interaction with the Ala20(NH)-Met1(CO) hydrogen bonded pair at the beginning of an anti-parallel beta-sheet, and the other with the Tyr33(NH)-Val29(CO) hydrogen bonded pair in the single alpha-helix. The latter, which is external and solvent accessible, is associated with a distortion in the alpha-helix centered around Tyr33 which consists of a significant increase in the CO(i-4)-N(i) and CO(i-4)-NH(i) distances relative to those in the rest of the helix, as well as a significant departure in the phi, psi angles of Tyr33 relative to regular helical geometry. Such solvent induced distortions in alpha-helices have been previously noticed in crystal structures and were postulated as possible folding intermediates for helical structures. The present observation of this phenomenon in solution indicates, however, that these water molecules are tightly bound and represent an integral part of the protein framework.     

Mimicking the action of folding chaperones in molecular dynamics simulations: Application to the refinement of homology-based protein structures

A novel method for the refinement of misfolded protein structures is proposed in which the properties of the solvent environment are oscillated in order to mimic some aspects of the role of molecular chaperones play in protein folding in vivo. Specifically, the hydrophobicity of the solvent is cycled by repetitively altering the partial charges on solvent molecules (water) during a molecular dynamics simulation. During periods when the hydrophobicity of the solvent is increased, intramolecular hydrogen bonding and secondary structure formation are promoted. During periods of increased solvent polarity, poorly packed regions of secondary structures are destabilized, promoting structural rearrangement. By cycling between these two extremes, the aim is to minimize the formation of long-lived intermediates. The approach has been applied to the refinement of structural models of three proteins generated by using the ROSETTA procedure for ab initio structure prediction. A significant improvement in the deviation of the model structures from the corresponding experimental structures was observed. Although preliminary, the results indicate computationally mimicking some functions of molecular chaperones in molecular dynamics simulations can promote the correct formation of secondary structure and thus be of general use in protein folding simulations and in the refinement of structural models of small- to medium-size proteins.     

A Ser residue influences the structure and stability of a Pro-kinked transmembrane helix dimer

When localized adjacent to a Pro-kink, Thr and Ser residues can form hydrogen bonds between their polar hydroxyl group and a backbone carbonyl oxygen and thereby modulate the actual bending angle of a distorted transmembrane α-helix. We have used the homo-dimeric transmembrane cytochrome b(559)' to analyze the potential role of a highly conserved Ser residue for assembly and stabilization of transmembrane proteins. Mutation of the conserved Ser residue to Ala resulted in altered heme binding properties and in increased stability of the holo-protein, most likely by tolerating subtle structural rearrangements upon heme binding. The results suggest a crucial impact of an intrahelical Ser hydrogen bond in defining the structure of a Pro-kinked transmembrane helix dimer.     

Alpha-helix stability in proteins. I. Empirical correlations concerning substitution of side-chains at the N and C-caps and the replacement of alanine by glycine or serine at solvent-exposed surfaces.

The importance of amino acid side-chains in helix stability has been investigated by making a series of mutations at the N-caps, C-caps and internal positions of the solvent-exposed faces of the two alpha-helices of barnase. There is a strong positional and context dependence of the effect of a particular amino acid on stability. Correlations have been found that provide insight into the physical basis of helix stabilization. The relative effects of Ala and Gly (or Ser) may be rationalized on the basis of solvent-accessible surface areas: burial of hydrophobic surface stabilizes the protein as does exposure to solvent of unpaired hydrogen bond donors or acceptors in the protein. There is a good correlation between the relative stabilizing effects of Ala and Gly at internal positions with the total change in solvent-accessible hydrophobic surface area of the folded protein on mutation of Ala----Gly. The relationship may be extended to the N and C-caps by including an extra term in hydrophilic surface area for the solvent exposure of the non-intramolecularly hydrogen-bonded main-chain CO, NH or protein side-chain hydrogen bonding groups. The requirement for solvent exposure of the C-cap main-chain CO groups may account for the strong preference for residues having positive phi and psi angles at this position, since this alpha L-conformation results in the largest solvent exposure of the C-terminal CO groups. Glycine in an alpha L-conformation results in the greatest exposure of these CO groups. Further, the side-chains of His, Asn, Arg and Lys may, with positive phi and psi-angles, form a hydrogen bond with the backbone CO of residue in position C -3 (residues are numbered relative to the C-cap). The preferences at the C-cap are Gly much greater than His greater than Asn greater than Arg greater than Lys greater than Ala approximately Ser approximately greater than Asp. The preferences at the N-cap are determined by hydrogen bonding of side-chains or solvent to the exposed backbone NH groups and are: Thr approximately Asp approximately Ser greater than Gly approximately Asn greater than Gln approximately Glu approximately His greater than Ala greater than Val much greater than Pro. These general trends may be obscured when mutation allows another side-chain to become a surrogate cap.(ABSTRACT TRUNCATED AT 400 WORDS)     

Analysis of protein main-chain solvation as a function of secondary structure

We have analysed the hydration of main-chain carbonyl and amide groups in 24 high-resolution well-refined protein structures as a function of the secondary structure in which these polar groups occur. We find that main-chain atoms in beta-sheets are as hydrated as those in alpha-helices, with most interactions involving "free" amide and carbonyl groups that do not participate in secondary structure hydrogen bonds. The distributions of water molecules around these non-bonded carbonyl groups reflect specific steric interactions due to the local secondary structure. Approximately 20% and 4%, respectively of bonded carbonyl and amide groups interact with solvent. These include interactions with carbonyl groups on the exposed faces of alpha-helices that have been correlated previously with bending of the helix. Water molecules interacting with alpha-helices occur mainly at the amino and carbonyl termini of the helices, in which case the solvent sites maintain the hydrogen bonding by bridging between residues i and i-3 or i-4 at the amino terminus and between i and i+3 or i+4 at the carbonyl terminus. We also see a number of solvent-mediated Ncap and Ccap interactions. The water molecules interacting with beta-sheets occur mainly at the edges, in which case they extend the sheet structure, or at the ends of strands, in which case they extend the beta-ladder. In summary, the solvent networks appear to extend the hydrogen-bonding structure of the secondary structures. In beta-turns, which usually occur at the surface of a protein, exposed amide and carbonyl groups are often hydrated, especially close to glycine residues. Occasionally water molecules form a bridge between residues i and i+3 in the turn and this may provide extra stabilization.     

Contemporary strategies for the stabilization of peptides in the alpha-helical conformation

Herein we review contemporary synthetic and protein design strategies to stabilize the alpha-helical motif in short peptides and miniature proteins. Advances in organometallic catalyst design, specifically for the olefin metathesis reaction, enable the use of hydrocarbon bridges to either crosslink side chains of specific residues or mimic intramolecular hydrogen bonds with carbon-carbon bonds. The resulting hydrocarbon-stapled and hydrogen bond surrogate alpha-helices provide unique synthetic ligands for targeting biomolecules. In the protein design realm, several classes of miniature proteins that display stable helical domains have been engineered and manipulated with powerful in vitro selection technologies to yield libraries of sequences that retain their helical folds. Rational re-design of these scaffolds provide distinctive reagents for the modulation of protein-protein interactions.     

Influence of solvent molecules on the stereochemical code of glycyl residues in proteins

The Ramachandran steric map and energy diagrams of the glycyl residue are symmetric. A plot of (phi,psi) angles of glycyl residues in 250 nonhomologous and high-resolution protein structures is also largely symmetric. However, there is a clear aberration in the symmetry. Although there is a cluster of points corresponding to the right-handed alpha-helical region, the "equivalent" cluster is clearly shifted to in and around the (phi,psi) values of (90 degrees, 0 degrees ) instead of being centered at the left-handed alpha-helical region of (60 degrees, 40 degrees ). This lack of symmetry exists even in the (phi,psi) distribution of residues from non-alpha-helical regions in proteins. Here we provide an explanation for this observation. An analysis of glycyl conformations in small peptide structures and in "coil" proteins, which are largely devoid of helical and sheet regions, shows that glycyl residues prefer to adopt conformations around (+/-90 degrees, 0 degrees ) instead of right- and left-handed alpha-helical regions. By using theoretical calculations, such conformations are shown to have highest solvent accessibility in a system of two-linked peptide units with glycyl residue at the central C(alpha) atom. This finding is consistent with the observations from 250 nonhomologous protein structures where glycyl residues with conformations close to (+/-90 degrees, 0 degrees ) are seen to have high solvent accessibility. Analysis of a subset of nonhomologous structures with very high resolution (1.5 A or better) shows that water molecules are indeed present at distances suitable for hydrogen bond interaction with glycyl residues possessing conformations close to (+/-90 degrees, 0 degrees ). It is suggested that water molecules play a key role in determining and stabilizing these conformations of glycyl residues and explain the aberration in the symmetry of glycyl conformations in proteins.     

Theoretical conformational analysis of methylamide of N-acetyl-L-lysine]

The spatial structure of the methylamide of N-acetyl-L-lysine has been analysed taking into account non-bonded and electrostatic interactions, torsional energy, bond angles distortion and hydrogen bonding. Conformational capacities of the backbone and mutual dependence of spatial structures of the backbone and the side chain was described by conformational maps obtained by energy minimisation, the dihedral angles and the bond angles of the side chain being varied for every phi, psi point. Every possible combination for phi, psi, x1-x5-angles was used corresponding to the stable form of the backbone and to torsion potential minima of the initial approximations in the calculation of preferred conformations of the molecule. Comparisons are made between stable forms of the methylamide of N-acetyl-L-lysine and Lys residues in proteins with known structure.     

Intrahelical hydrogen bonding of serine, threonine and cysteine residues within alpha-helices and its relevance to membrane-bound proteins

A survey of known protein structures reveals that approximately 70% of serine residues and at least 85% (potentially 100%) of threonine residues in helices make hydrogen bonds to carbonyl oxygen atoms in the preceding turn of the helix. The high frequency of intrahelical hydrogen bonding is of particular significance for intrinsic membrane-bound proteins that form transmembrane helices. Hydrogen bonding within a helix provides a way for serine, threonine and cysteine residues to satisfy their hydrogen-bonding potential permitting such residues to occur in helices buried within a hydrophobic milieu.     

Hydrogen bonding and equilibrium protium-deuterium fractionation factors in the immunoglobulin G binding domain of protein G

Protium-deuterium fractionation factors (phi) were determined for more than 85% of the backbone amide protons in the IgG binding domains of protein G, GB1 and GB2, from NMR spectra recorded over a range of H2O/D2O solvent ratios. Previous studies suggest a correlation between phi and hydrogen bond strength; amide and hydroxyl groups in strong hydrogen bonds accumulate protium (phi < 1), while weak hydrogen bonds accumulate deuterium (phi > 1). Our results show that the alpha-helical residues have slightly lower phi values (1.03 +/- 0.05) than beta-sheet residues (1.12 +/- 0.07), on average. The lowest phi value obtained (0.65) does not involve a backbone amide but rather is for the interaction between two side chains, Y45 and D47. Fractionation factors for solvent-exposed residues are between the alpha-helix and beta-sheet values, on average, and are close to those for random coil peptides. Further, the difference in phiav between alpha-helix and solvent-exposed residues is small, suggesting that differences in hydrogen bond strength for intrachain hydrogen bonds and amide...water hydrogen bonds are also small. Overall, the enrichment for deuterium suggests that most backbone...backbone hydrogen bonds are weak.     

Ser/Thr Motifs in Transmembrane Proteins: Conservation Patterns and Effects on Local Protein Structure and Dynamics

We combined systematic bioinformatics analyses and molecular dynamics simulations to assess the conservation patterns of Ser and Thr motifs in membrane proteins, and the effect of such motifs on the structure and dynamics of α-helical transmembrane (TM) segments. We find that Ser/Thr motifs are often present in β-barrel TM proteins. At least one Ser/Thr motif is present in almost half of the sequences of α-helical proteins analyzed here. The extensive bioinformatics analyses and inspection of protein structures led to the identification of molecular transporters with noticeable numbers of Ser/Thr motifs within the TM region. Given the energetic penalty for burying multiple Ser/Thr groups in the membrane hydrophobic core, the observation of transporters with multiple membrane-embedded Ser/Thr is intriguing and raises the question of how the presence of multiple Ser/Thr affects protein local structure and dynamics. Molecular dynamics simulations of four different Ser-containing model TM peptides indicate that backbone hydrogen bonding of membrane-buried Ser/Thr hydroxyl groups can significantly change the local structure and dynamics of the helix. Ser groups located close to the membrane interface can hydrogen bond to solvent water instead of protein backbone, leading to an enhanced local solvation of the peptide.     

Refinement of the structure of recombinant rat intestinal fatty acid-binding apoprotein at 1.2-A resolution.

The three-dimensional structure of the 131-residue rat intestinal fatty acid-binding protein, without bound ligand (apoI-FABP), has been refined with x-ray diffraction data to a nominal resolution of 1.19 A. The final model has a conventional crystallographic R-factor of 16.9% for 34,290 unique reflections [a root mean square (r.m.s.) deviation for bond length of 0.012 A and a r.m.s. deviation of 2.368 degrees for bond angles]. Ninety-two residues are present as components of the protein's 10 anti-parallel beta-strands while 14 residues are part of its two short alpha-helices. The beta-strands and alpha-helices are organized into two nearly orthogonal beta-sheets. Particular attention has been placed in defining solvent structure and the structures of discretely disordered groups in this protein. Two hundred thirty-seven solvent molecules have been identified; 24 are located within apoI-FABP. The refined model includes alternate conformers for 228 protein atoms (109 main-chain, 119 side-chain) and 63 solvent molecules. We have found several aromatic side-chains with multiple conformations located near, or in, the protein's ligand binding site. This observation, along with the fact that these side-chains have a temperature factor that is relatively higher than that of other aromatic residues, suggests that they may be involved in the process of noncovalent binding of fatty acid. The absence of a true hydrophobic core in I-FABP suggests that its structural integrity may be maintained primarily by a hydrogen bonding network involving protein and solvent atoms.     

Exact Solutions for Internuclear Vectors and Backbone Dihedral Angles from NH Residual Dipolar Couplings in Two Media, and their Application in a Systematic Search Algorithm for Determining Protein Backbone Structure

We have derived a quartic equation for computing the direction of an internuclear vector from residual dipolar couplings (RDCs) measured in two aligning media, and two simple trigonometric equations for computing the backbone (phi,psi) angles from two backbone vectors in consecutive peptide planes. These equations make it possible to compute, exactly and in constant time, the backbone (phi,psi) angles for a residue from RDCs in two media on any single backbone vector type. Building upon these exact solutions we have designed a novel algorithm for determining a protein backbone substructure consisting of alpha-helices and beta-sheets. Our algorithm employs a systematic search technique to refine the conformation of both alpha-helices and beta-sheets and to determine their orientations using exclusively the angular restraints from RDCs. The algorithm computes the backbone substructure employing very sparse distance restraints between pairs of alpha-helices and beta-sheets refined by the systematic search. The algorithm has been demonstrated on the protein human ubiquitin using only backbone NH RDCs, plus twelve hydrogen bonds and four NOE distance restraints. Further, our results show that both the global orientations and the conformations of alpha-helices and beta-strands can be determined with high accuracy using only two RDCs per residue. The algorithm requires, as its input, backbone resonance assignments, the identification of alpha-helices and beta-sheets as well as sparse NOE distance and hydrogen bond restraints.     

A solvent model for simulations of peptides in bilayers. II. Membrane-spanning alpha-helices

We describe application of the implicit solvation model (see the first paper of this series), to Monte Carlo simulations of several peptides in bilayer- and water-mimetic environments, and in vacuum. The membrane-bound peptides chosen were transmembrane segments A and B of bacteriorhodopsin, the hydrophobic segment of surfactant lipoprotein, and magainin2. Their conformations in membrane-like media are known from the experiments. Also, molecular dynamics study of surfactant lipoprotein with different explicit solvents has been reported (Kovacs, H., A. E. Mark, J. Johansson, and W. F. van Gunsteren. 1995. J. Mol. Biol. 247:808-822). The principal goal of this work is to compare the results obtained in the framework of our solvation model with available experimental and computational data. The findings could be summarized as follows: 1) structural and energetic properties of studied molecules strongly depend on the solvent; membrane-mimetic media significantly promote formation of alpha-helices capable of traversing the bilayer, whereas a polar environment destabilizes alpha-helical conformation via reduction of solvent-exposed surface area and packing; 2) the structures calculated in a membrane-like environment agree with the experimental ones; 3) noticeable differences in conformation of surfactant lipoprotein assessed via Monte Carlo simulation with implicit solvent (this work) and molecular dynamics in explicit solvent were observed; 4) in vacuo simulations do not correctly reproduce protein-membrane interactions, and hence should be avoided in modeling membrane proteins.