Understanding protein hydrogen bond formation with kinetic H/D amide isotope effects

Through the development of a procedure to measure when hydrogen bonds form under two-state folding conditions, alpha-helices have been determined to form proportionally to denaturant-sensitive surface area buried in the transition state. Previous experiments assessing H/D isotope effects are applied to various model proteins, including lambda and Arc repressor variants, a coiled coil domain, cytochrome c, colicin immunity protein 7, proteins L and G, acylphosphatase, chymotrypsin inhibitor II and a Src SH3 domain. The change in free energy accompanied by backbone deuteration is highly correlated to secondary structure composition when hydrogen bonds are divided into two classes. The number of helical hydrogen bonds correlates with an average equilibrium isotope effect of 8.6 +/- 0.9 cal x mol(-1) x site(-1). However, beta-sheet and long-range hydrogen bonds have little isotope effect. The kinetic isotope effects support our hypothesis that, for helical proteins, hydrophobic association cannot be separated from helix formation in the transition state. Therefore, folding models that describe an incremental build-up of structure in which hydrophobic burial and hydrogen bond formation occur commensurately are more consistent with the data than are models that posit the extensive formation of one quantity before the other.     

Equilibrium hydrogen exchange reveals extensive hydrogen bonded secondary structure in the on-pathway intermediate of Im7

The four-helical immunity protein Im7 folds through an on-pathway intermediate that has a specific, but partially misfolded, hydrophobic core. In order to gain further insight into the structure of this species, we have identified the backbone hydrogen bonds formed in the ensemble by measuring the amide exchange rates (under EX2 conditions) of the wild-type protein and a variant, I72V. In this mutant the intermediate is significantly destabilised relative to the unfolded state (deltadeltaG(ui) = 4.4 kJ/mol) but the native state is only slightly destabilised (deltadeltaG(nu) = 1.8 kJ/mol) at 10 degrees C in 2H2O, pH* 7.0 containing 0.4 M Na2SO4, consistent with the view that this residue forms significant non-native stabilising interactions in the intermediate state. Comparison of the hydrogen exchange rates of the two proteins, therefore, enables the state from which hydrogen exchange occurs to be identified. The data show that amides in helices I, II and IV in both proteins exchange slowly with a free energy similar to that associated with global unfolding, suggesting that these helices form highly protected hydrogen-bonded helical structure in the intermediate. By contrast, amides in helix III exchange rapidly in both proteins. Importantly, the rate of exchange of amides in helix III are slowed substantially in the Im7* variant, I72V, compared with the wild-type protein, whilst other amides exchange more rapidly in the mutant protein, in accord with the kinetics of folding/unfolding measured using chevron analysis. These data demonstrate, therefore, that local fluctuations do not dominate the exchange mechanism and confirm that helix III does not form stable secondary structure in the intermediate. By combining these results with previously obtained Phi-values, we show that the on-pathway folding intermediate of Im7 contains extensive, stable hydrogen-bonded structure in helices I, II and IV, and that this structure is stabilised by both native and non-native interactions involving amino acid side-chains in these helices.     

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

The Ribbon of Hydrogen Bonds in the Three-Dimensional Structure of Globular Proteins

We report quantum mechanical computations and experimental evidence which suggest that the backbone conformation of globular proteins depends generally on the conservation of that part of the hydrogen bond network or ribbon which is joined, in general, directly to the backbone and is largely independent of the remainder of this whole network of hydrogen bonds. The familiar hydrogen bonds of the helix and the sheet form about one-half of this ribbon of hydrogen bonds. Both water molecules and hydrogen bonding side chain groups are involved in the formation of the ribbon.This view of the three-dimensional structure of globular proteins in terms of the `molecule' allows us to deal with the non-secondary structure as well as with the familiar secondary structure. It also suggests that the ribbon contains approximately the same number of hydrogen bonds within all three structures – the helix, the sheet and the coil – and that this is the reason for the ease of interconversion of these three structures.The quantum mechanical computations on hydrogen bonding suggest that delocalised water molecules which have substantial mobility are an essential part of the ribbon. This situation arises because the hydrogen bonding groups of the protein molecule are not free to move to optimise the hydrogen bonding geometries as are the oxygen atoms in the waters and ices. Such delocalised water molecules either have high B values or are invisible in the X-ray data and yet are able to form a structure which is as strong as a normal hydrogen bond.The experimental data on the point mutations of the THRI57 residue of the T4 phage lysome provides an initial test of this model. Both the local backbone conformation and the ribbon of hydrogen bonds are conserved throughout all the mutations of residue 157,providing that the delocalised water molecules are accepted as a genuine part of the structure. These mutations include the introduction of hydrocarbon side chains at position 157 when water molecules or other side chain groups take over the formation of the hydrogen bonds.We suggest that, provided steric effects are not important, many point mutations succeed because they leave the ribbon of hydrogen bonds (and so the backbone conformation) largely unchanged.     

Helix capping.

Helix-capping motifs are specific patterns of hydrogen bonding and hydrophobic interactions found at or near the ends of helices in both proteins and peptides. In an alpha-helix, the first four >N-H groups and last four >C=O groups necessarily lack intrahelical hydrogen bonds. Instead, such groups are often capped by alternative hydrogen bond partners. This review enlarges our earlier hypothesis (Presta LG, Rose GD. 1988. Helix signals in proteins. Science 240:1632-1641) to include hydrophobic capping. A hydrophobic interaction that straddles the helix terminus is always associated with hydrogen-bonded capping. From a global survey among proteins of known structure, seven distinct capping motifs are identified-three at the helix N-terminus and four at the C-terminus. The consensus sequence patterns of these seven motifs, together with results from simple molecular modeling, are used to formulate useful rules of thumb for helix termination. Finally, we examine the role of helix capping as a bridge linking the conformation of secondary structure to supersecondary structure.     

Invariant features of globin primary structure and coding of their secondary structure

Invariant features of the primary structure of 67 globins are analysed. These features may be responsible for the formation of the secondary structure of these proteins at the first stage of self-organization (in the unfolded chain). It is shown that in primary structures of globins there are 11 sites or regions of one to four residues in which at least one of the residues Asn, Asp, His, Pro, Ser or Thr is located in every globin (haem-linking His residues are excluded from these sites). An unambiguous correlation exists between the position of these regions and the secondary structure of globins: all these regions (except one) are located near the ends of helices in globins whose three-dimensional structure is known and the ends of all helices (except for the helix F) are coded by such regions. A decrease in the set of residues listed above leads to a sharp drop in the number of regions invariantly occupied by the residues, while an addition of residues such as Tyr and Gly to this set does not eventually increase the number of invariant regions. Five residues (Asn, Asp, His, Ser and Thr) of the six that code the ends of helical regions have polar side groups with a small number of degrees of freedom capable of forming hydrogen bonds with atoms of the backbone with a relatively small loss of entropy. One residue (Pro) has no NH-group and, therefore, has less chance of participating in the formation of hydrogen bonds between atoms of the backbone. This corroborates the hypothesis that competition between hydrogen bonds of short polar side groups and hydrogen bonds in the backbone is essential for the formation of the secondary structure in unfolded protein chains. Amino acid replacements in hydrophobic cores of the 67 globins are considered in the Appendix.     

Helices in peptoids of alpha- and beta-peptides

Peptoids of alpha- and beta-peptides (alpha- and beta-peptoids) can be obtained by shifting the amino acid side chains from the backbone carbon atoms of the monomer constituents to the peptide nitrogen atoms. They are, therefore, N-substituted poly-glycines and poly-beta-alanines, respectively. Due to the substituted nitrogen atoms, the ability for hydrogen bond formation between peptide bonds gets lost. It may be very interesting to see whether such non-natural oligomers could be regarded as foldamers, which fold into definite backbone conformers. In this paper, we provide a complete overview on helix formation in alpha- and beta-peptoids on the basis of systematic theoretical conformational analyses employing the methods of ab initio molecular orbital (MO) theory. It can be shown that the alpha- and beta-peptoid structures form helical structures with both trans and cis peptide bonds despite the missing hydrogen bonds. Obviously, the conformational properties of the backbone are more important for folding than the possibility of hydrogen bonding. There are close relationships between the helices of alpha-peptoids and poly-glycine and poly-proline helices of alpha-peptides, whereas the helices of beta-peptoids correspond to the well-known helical structures of beta-peptides as, for instance, the 3(1)-helix of beta-peptides with 14-membered hydrogen-bonded rings. Thus, alpha- and beta-peptoids enrich the field of foldamers and may be used as useful tools in peptide and protein design.     

Stereochemical analysis of the secondary structure of polypeptide chains with the aid of Courtauld three-dimensional models. II. Hydrogen and hydrophobic bonds.

Sterically permissible hydrogen bonds between side chains and the backbone, which fix a small number of angles of internal rotation have been identified. The hydrogen bonds on the ends and within the secondary structures, as well as the most important bonds in the alpha helix and in irregular structures, have been considered.     

Statistical and molecular dynamics studies of buried waters in globular proteins

Buried solvent molecules are common in the core of globular proteins and contribute to structural stability. Folding necessitates the burial of polar backbone atoms in the protein core, whose hydrogen-bonding capacities should be satisfied on average. Whereas the residues in alpha-helices and beta-sheets form systematic main-chain hydrogen bonds, the residues in turns, coils and loops often contain polar atoms that fail to form intramolecular hydrogen bonds. The statistical analysis of 842 high resolution protein structures shows that well-resolved, internal water molecules preferentially reside near residues without alpha-helical and beta-sheet secondary structures. These buried waters most often form primary hydrogen bonds to main-chain atoms not involved in intramolecular hydrogen bonds, providing strong evidence that hydrating main-chain atoms is a key structural role of buried water molecules. Additionally, the average B-factor of protein atoms hydrogen-bonded to waters is smaller than that of protein atoms forming intramolecular hydrogen bonds, and the average B-factor of water molecules involved in primary hydrogen bonds with main-chain atoms is smaller than the average B-factor of water molecules involved in secondary hydrogen bonds to protein atoms that form concurrent intramolecular hydrogen bonds. To study the structural coupling between internal waters and buried polar atoms in detail we simulated the dynamics of wild-type FKBP12, in which a buried water, Wat137, forms one side-chain and multiple main-chain hydrogen bonds. We mutated E60, whose side-chain hydrogen bonds with Wat137, to Q, N, S or A, to modulate the multiplicity and geometry of hydrogen bonds to the water. Mutating E60 to a residue that is unable to form a hydrogen bond with Wat137 results in reorientation of the water molecule and leads to a structural readjustment of residues that are both near and distant to the water. We predict that the E60A mutation will result in a significantly reduced affinity of FKBP12 for its ligand FK506. The propensity of internal waters to hydrogen bond to buried polar atoms suggests that ordered water molecules may constitute fundamental structural components of proteins, particularly in regions where alpha-helical or beta-sheet secondary structure is not present.     

High-resolution study of the three-dimensional structure of horse heart metmyoglobin

The three-dimensional structure of horse heart metmyoglobin has been refined to a final R-factor of 15.5% for all observed data in the 6.0 to 1.9 A resolution range. The final model consists of 1242 non-hydrogen protein atoms, 154 water molecules and one sulfate ion. This structure has nearly ideal bonding and bond angle geometry. A Luzzati plot of the variation in R-factor with resolution yields an estimated mean co-ordinate error of 0.18 A. An extensive analysis of the pattern of hydrogen bonds formed in horse heart metmyoglobin has been completed. Over 80% of the polypeptide chain is involved in eight helical segments, of which seven are composed mainly of alpha-helical (3.6(13))-type hydrogen bonds; the remaining helix is composed entirely of 3(10) hydrogen bonds. Altogether, of 102 hydrogen bonds between main-chain atoms only six are not involved in helical structures, and four of these six occur within beta-turns. The majority of water molecules in horse heart metmyoglobin are found in solvent networks that range in size from two to 35 members. The size of water molecule networks can be rationalized on the basis of three factors: the number of hydrogen bonds to the protein surface, the presence of charged side-chain atoms, and the ability to bridge to neighboring molecules in the crystal lattice. Bridging water networks form the dominant intermolecular interactions. The backbone conformation of horse heart metmyoglobin is very similar to sperm whale metmyoglobin, with significant differences in secondary structure occurring only near residues 119 and 120, where residues 120 to 123 in sperm whale form a distorted type I reverse turn and the horse heart protein has a type II turn at residues 119 to 122. Nearly all of the hydrogen bonds between main-chain atoms (occurring mainly in helical regions) are common to both proteins, and more than half of the hydrogen bonds involving side-chain atoms observed in horse heart are also found in sperm whale metmyoglobin. Unlike sperm whale metmyoglobin, the heme iron atom in horse heart metmyoglobin is not significantly displaced from the plane of the heme group.     

Accessible surfaces of beta proteins increase with increasing protein molecular mass more rapidly than those of other proteins

Here we present a systematic analysis of accessible surface areas and hydrogen bonds of 2554 globular proteins from four structural classes (all-α, all-β, α/β and α+β proteins) that is aimed to learn in which structural class the accessible surface area increases with increasing protein molecular mass more rapidly than in other classes, and what structural peculiarities are responsible for this effect. The beta structural class of proteins was found to be the leader, with the following possible explanations of this fact. First, in beta structural proteins, the fraction of residues not included in the regular secondary structure is the largest, and second, the accessible surface area of packaged elements of the beta-structure increases more rapidly with increasing molecular mass in comparison with the alpha-structure. Moreover, in the beta structure, the probability of formation of backbone hydrogen bonds is higher than that in the alpha helix for all residues of α+β proteins (the average probability is 0.73±0.01 for the beta-structure and 0.60±0.01 for the alpha-structure without proline) and α/β proteins, except for asparagine, aspartic acid, glycine, threonine, and serine (0.70±0.01 for the beta-structure and 0.60±0.01 for the alpha-structure without the proline residue). There is a linear relationship between the number of hydrogen bonds and the number of amino acid residues in the protein (Number of hydrogen bonds=0.678·number of residues-3.350).     

The influence of protonation states on the dynamics of the NhaA antiporter from Escherichia coli

The crystal structure of NhaA Na(+)/H(+) antiporter of Escherichia coli has provided a basis to explore the mechanism of Na(+) and H(+) exchange and its regulation by pH. However, the dynamics and nature of the pH-induced changes in the proteins remained unknown. Using molecular mechanics methods, we studied the dynamic behavior of the hydrogen-bonded network in NhaA on shifting the pH from 4 to 8. The helical regions preserved the general architecture of NhaA throughout the pH change. In contrast, large conformational drifts occurred at pH 8 in the loop regions, and an increased flexibility of helix IVp was observed on the pH shift. A remarkable pH-induced conformational reorganization was found: at acidic pH helix X is slightly curved, whereas at alkaline pH, it is kinked around residue Lys(300). The barrier that exists between the cytoplasmic and periplasmic funnels at low pH is removed, and the two funnels are bridged by hydrogen bonds between water molecules and residues located in the TMSs IV/XI assembly and helix X at alkaline pH. In the variant Gly(338)Ser that lost pH control, a hydrogen-bonded chain between Ser(338) and Lys(300) was found to block the pH-induced conformational reorganization of helix X.     

Cooperative effects in hydrogen-bonding of protein secondary structure elements: a systematic analysis of crystal data using Secbase

A systematic analysis of the hydrogen-bonding geometry in helices and beta sheets has been performed. The distances and angles between the backbone carbonyl O and amide N atoms were correlated considering more than 1500 protein chains in crystal structures determined to a resolution better than 1.5 A. They reveal statistically significant trends in the H-bond geometry across the different secondary structural elements. The analysis has been performed using Secbase, a modular extension of Relibase (Receptor Ligand Database) which integrates information about secondary structural elements assigned to individual protein structures with the various search facilities implemented into Relibase. A comparison of the mean hydrogen-bond distances in alpha helices and 3(10) helices of increasing length shows opposing trends. Whereas in alpha helices the mean H-bond distance shrinks with increasing helix length and turn number, the corresponding mean dimension in 3(10) helices expands in a comparable series. Comparing similarly the hydrogen-bond lengths in beta sheets there is no difference to be found between the mean H-bond length in antiparallel and parallel beta sheets along the strand direction. In contrast, an interesting systematic trend appears to be given for the hydrogen bonds perpendicular to the strands bridging across an extended sheet. With increasing number of accumulated strands, which results in a growing number of back-to-back piling hydrogen bonds across the strands, a slight decrease of the mean H-bond distance is apparent in parallel beta sheets whereas such trends are obviously not given in antiparallel beta sheets. This observation suggests that cooperative effects mutually polarizing spatially well-aligned hydrogen bonds are present either in alpha helices and parallel beta sheets whereas such influences seem to be lacking in 3(10) helices and antiparallel beta sheets.     

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.     

Investigation of the backbone dynamics of the IgG-binding domain of streptococcal protein G by heteronuclear two-dimensional 1H-15N nuclear magnetic resonance spectroscopy.

The backbone dynamics of the immunoglobulin-binding domain (B1) of streptococcal protein G, uniformly labeled with 15N, have been investigated by two-dimensional inverse detected heteronuclear 1H-15N NMR spectroscopy at 500 and 600 MHz. 15N T1, T2, and nuclear Overhauser enhancement data were obtained for all 55 backbone NH vectors of the B1 domain at both field strengths. The overall correlation time obtained from an analysis of the T1/T2 ratios was 3.3 ns at 26 degrees C. Overall, the B1 domain is a relatively rigid protein, consistent with the fact that over 95% of the residues participate in secondary structure, comprising a four-stranded sheet arranged in a -1, +3x, -1 topology, on top of which lies a single helix. Residues in the turns and loops connecting the elements of secondary structure tend to exhibit a higher degree of mobility on the picosecond time scale, as manifested by lower values of the overall order parameter. A number of residues at the ends of the secondary structure elements display two distinct internal motions that are faster than the overall rotational correlation time: one is fast (< 20 ps) and lies in the extreme narrowing limit, whereas the other is one to two orders of magnitude slower (1-3 ns) and lies outside the extreme narrowing limit. The slower motion can be explained by large-amplitude (20-40 degrees) jumps in the N-H vectors between states with well-defined orientations that are stabilized by hydrogen bonds.(ABSTRACT TRUNCATED AT 250 WORDS)     

High-Resolution Crystal Structures of Protein Helices Reconciled with Three-Centered Hydrogen Bonds and Multipole Electrostatics

Theoretical and experimental evidence for non-linear hydrogen bonds in protein helices is ubiquitous. In particular, amide three-centered hydrogen bonds are common features of helices in high-resolution crystal structures of proteins. These high-resolution structures (1.0 to 1.5 Å nominal crystallographic resolution) position backbone atoms without significant bias from modeling constraints and identify Φ = -62°, ψ = -43 as the consensus backbone torsional angles of protein helices. These torsional angles preserve the atomic positions of α-β carbons of the classic Pauling α-helix while allowing the amide carbonyls to form bifurcated hydrogen bonds as first suggested by Némethy et al. in 1967. Molecular dynamics simulations of a capped 12-residue oligoalanine in water with AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Applications), a second-generation force field that includes multipole electrostatics and polarizability, reproduces the experimentally observed high-resolution helical conformation and correctly reorients the amide-bond carbonyls into bifurcated hydrogen bonds. This simple modification of backbone torsional angles reconciles experimental and theoretical views to provide a unified view of amide three-centered hydrogen bonds as crucial components of protein helices. The reason why they have been overlooked by structural biologists depends on the small crankshaft-like changes in orientation of the amide bond that allows maintenance of the overall helical parameters (helix pitch (p) and residues per turn (n)). The Pauling 3.6₁₃ α-helix fits the high-resolution experimental data with the minor exception of the amide-carbonyl electron density, but the previously associated backbone torsional angles (Φ, Ψ) needed slight modification to be reconciled with three-atom centered H-bonds and multipole electrostatics. Thus, a new standard helix, the 3.6_13/10-, Némethy- or N-helix, is proposed. Due to the use of constraints from monopole force fields and assumed secondary structures used in low-resolution refinement of electron density of proteins, such structures in the PDB often show linear hydrogen bonding.     

Computer modeling of polyleucine-based coiled coil dimers in a realistic membrane environment: insight into helix-helix interactions in membrane proteins

Simulated annealing was performed to model parallel dimers of alpha-helical transmembrane peptides with the sequence L(11)XL(12), predicting left-handed coiled coil geometry in all cases. Insertion of peptides containing threonine, asparagine, alanine, phenylalanine, and leucine in position 12 into realistic model membranes showed these structures were stable for 20 ns of molecular dynamics simulation time. Threonine could participate in intermolecular hydrogen bonds, but predominantly formed hydrogen bonds to the backbone of the helix it resided on. These hydrogen bonds, although infrequent, appeared to promote closer association of polyleucine helices. Asparagine participated in multiple, rapidly fluctuating intermolecular and intramolecular hydrogen bonds, and may have slightly destabilized optimum van der Waals packing in favor of optimum hydrogen bonding. Coordinated rotations of transmembrane helices about their axes were observed, indicating helices may rotate around one another during the folding of membrane proteins or other processes. These rotations were inhibited by phenylalanine, suggesting a role for bulky residues in modulating membrane protein dynamics.     

Stabilization of conformationally dynamic helices by covalently attached acyl chains

Acylation of proteins is known to mediate membrane attachment and to influence subcellular sorting. Here, we report that acylation can stabilize secondary structure. Circular dichroism spectroscopy showed that N‐terminal attachment of acyl chains decreases the ability of an intrinsically flexible hydrophobic model peptide to refold from an α‐helical state to β‐sheet in response to changing solvent conditions. Acylation also stabilized the membrane‐embedded α‐helix. This increase of global helix stability did not result from decreased local conformational dynamics of the helix backbone as assessed by deuterium/hydrogen‐exchange experiments. We concluded that acylation can stabilize the structure of intrinsically dynamic helices and may thus prevent misfolding.     

Energy profile of maltooligosaccharide permeation through maltoporin as derived from the structure and from a statistical analysis of saccharide-protein interactions.

The crystal structure of the maltodextrin-specific porin from Salmonella typhimurium ligated with a maltotrioside at the pore eyelet is known at 2.4 A resolution. The three glucose units assume a conformation close to the natural amylose helix. The pore eyelet fits exactly the cross-section of a maltooligosaccharide chain and thus functions as a constraining orifice. The oligomer permeates the membrane by screwing along the amylose helix through this orifice. Because each glucose glides along the given helix, its interactions can be sampled at any point along the pathway. The interactions are mostly hydrogen bonds, but also contacts to aromatic rings at one side of the pore. We have derived the energy profile of a gliding maltooligosaccharide by following formation and breakage of hydrogen bonds and by assessing the saccharide-aromatics interactions from a statistical analysis of saccharide binding sites in proteins. The resulting profile indicates smooth permeation despite extensive hydrogen bonding at the orifice.     

Probing backbone hydrogen bonds in the hydrophobic core of GCN4

Backbone amide hydrogen bonds play a central role in protein secondary and tertiary structure. Previous studies have shown that substitution of a backbone ester (-COO-) in place of a backbone amide (-CONH-) can selectively destabilize backbone hydrogen bonds in a protein while maintaining a similar conformation to the native backbone structure. The majority of these studies have focused on backbone substitutions that were accessible to solvent. The GCN4 coiled coil domain is an example of a stable alpha-helical dimer that possesses a well-packed hydrophobic core. Amino acids in the a and d positions of the GCN4 helix, which pack the hydrophobic core, were replaced with the corresponding alpha-hydroxy acids in the context of a chemoselectively ligated heterodimer. While the overall structure and oligomerization state of the heterodimer were maintained, the overall destabilization of the ester analogues was greater (average DeltaDeltaG of 3+ kcal mol(-1)) and more variable than previous studies. Since burial of the more hydrophobic ester should stabilize the backbone and reduce the DeltaDeltaG, the increased destabilization must come from another source. However, the observed destabilization is correlated with the protection factors for individual amide hydrogens from previous hydrogen exchange experiments. Therefore, our results suggest that backbone engineering through ester substitution is a useful approach for probing the relative strength of backbone hydrogen bonds.