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.     

Qualitative aspects of hydrogen-deuterium exchange in the 1H, 13C, and 15N nuclear magnetic resonance spectra of viomycin in aqueous solution.

The 1H, 13C, and 15N high field nuclear magnetic resonance spectra of the cyclic peptide viomycin have been fully assigned using homo- and heteronuclear double resonance experiments and pH effects. In addition it is shown how the two- and three-bond H-D isotope effects upon carbonyl resonances may assist in their assignment. The resistance to exchange with solvent water of the amide proton involved in the transannular hydrogen bond is observed directly in the 1H spectra, via the isotope effect on a carbonyl resonance in the 13C spectra, and via the one-bond 1H couppling in the 15N spectra.     

Hydrogen exchange of primary amide protons in basic pancreatic trypsin inhibitor: evidence for NH2 group rotation in buried asparagine side chains

Hydrogen-deuterium exchange is measured for the buried primary amide groups of Asn-43 and Asn-44 in bovine pancreatic trypsin inhibitor. Amide protons trans and cis to the amide carbonyl oxygen (HE and HZ, respectively) exchange at indistinguishable rates. Uncorrelated exchange of HE and HZ is established for both residues by following the nuclear Overhauser enhancement from HE to HZ during the deuterium exchange. The exchange of Asn-43 and Asn-44 side-chain protons differs qualitatively from exchange of primary amide groups in fully solvated model compounds, for which HE generally exchanges faster than HZ. The equal rates for the buried primary amide HE and HZ in BPTI are not a consequence of coupled exchange. The data indicate rapid rotation around the CO-NH2 bond for both Asn-43 and Asn-44 and suggest considerable lability of intramolecular hydrogen bonds. The side chain of Asn-43 has all of its polar atoms integrated into the very stable hydrogen-bonded structure of the protein. Asn-44 is hydrogen-bonded to side chains and to a buried water molecule. Solvent isotope exchange is several orders of magnitude more restricted by protein secondary and tertiary structure than the CO-NH2 rotation, indicating that N delta H2 groups flip many times before hydrogen isotope exchange occurs.     

Intramolecular interactions in chemically modified Escherichia coli thioredoxin monitored by hydrogen/deuterium exchange and electrospray ionization mass spectrometry.

Site specific amide hydrogen/deuterium content of oxidized and reduced Escherichia colithioredoxin, and alkylated derivatives, Cys-32-ethylglutathionylated and Cys-32-ethylcysteinylated thioredoxins are measured, after exposure for 20 s to D(2)O/phosphate buffer (pH 5.7), by electrospray mass spectrometry. The degree of deuteration of Oxi-TRX and Red-TRX correlated with the rates of H/D exchange measured previously by NMR. The ethylcysteinyl modification was shown to minimally perturb the active site of the reduced protein, but showed more global effects on structures of alpha-helices and beta-strands distant from the site of modification. In contrast, the larger ethylglutathionyl group had little effect on the protein's overall conformation, but significantly affected the structure of loops close to the active site. A molecular model of GS-ethyl-TRX derived from molecular simulation allowed the H/D exchange results to be interpreted in terms of specific interactions between the alkyl chain and the protein surface. The specific conformation of the ethylglutathione modification was predicted to be fixed by salt bridges between the carboxylates of the gamma-Glu and Gly of glutathione and the guanidinium of Arg-73 and epsilon-amino group of Lys-90 of the protein. Specific hydrogen bonding interactions between the glutathione carbonyl oxygens and the amide protons of thioredoxin residues Ile-75 and Ala-93 were predicted. The H/D exchange studies showed low levels of deuterium incorporation at backbone nitrogens of these residues. The data also provided evidence for an unusual amide proton-amide nitrogen hydrogen bond within the ethylglutathionylated chain. These same sets of electrostatic and hydrogen bonding interactions were not predicted or observed for the smaller alkyl modification in Cys-ethyl-TRX.     

N.m.r. study on conformation of Ac-Thr(alpha-GalNAc)-Ala-Ala-OMe as a model for mucin type glycoprotein

This study report on the results of high resolution 1H n.m.r. investigations on Ac-Thr(alpha-GalNAc)-Ala-Ala-OMe 1 as a mucin type model glycopeptide of antifreeze glycoprotein (AFGP) in both dimethyl sulfoxide (DMSO) and H2O. The temperature dependence of amide proton chemical shifts strongly suggested the presence of the intramolecular hydrogen bond between the amide proton of GalNAc and the carbonyl oxygen of the Thr residues. Due to this bond, the orientation of the sugar residue of 1 appears to be fairly restricted relative to its peptide backbone. Despite the lack of the clear evidence for such intramolecular hydrogen bond in H2O, 1H coupling constant data suggested the structural similarity of 1 in DMSO and H2O, indicating the presence of the intramolecular hydrogen bond even in H2O, which may play an important role in determining the orientation of the sugar moiety with respect to the peptide backbone in glycoprotein.     

Conformation-Determining role for the N-acetyl group in the O-glycosidic linkage, α-GalNAc-Thr

An ¹H-nmr study of 2-acetamido-2-deoxy-3,4,6-tri-O-acetyl-D-galactopyranose (AcGalNAc) glycosylated Thr-containing tripeptides in Me₂SO-d₆ solution reveals two mutually exclusive intramolecular hydrogen bonds. In Z-Thr(AcGalNAc)-Ala-Ala-OMe, there is an intramolecular hydrogen bond between the Thr amide proton and the sugar N-acetyl carbonyl oxygen. The strength of this hydrogen bond will be dependent on the amino acid residues on the Thr C terminal side to some undetermined distance. In Ac-Thr(AcGalNAc)-Ala-Ala-OMe, a different intramolecular hydrogen bond between the sugar N-acetyl amide proton and the Thr carbonyl oxygen exists. The choice of hydrogen bonds seems dependent on the bulkiness of the residues on the Thr N terminal side. The consequence of such strong hydrogen bonds is a clearly defined orientation of the sugar moiety with respect to the peptide backbone. In the former, the plane of the sugar pyranose ring is roughly oriented perpendicularly to the peptide backbone. The latter orientation is where the plane of the sugar ring is roughly in line with the peptide backbone. In both orientations, the sugar moiety can increase the shielding of the neighboring amino acid residues from the solvent. The idea that the amino acid residues near the glycosylated Thr influence orientation of the sugar moiety with respect to the peptide backbone and in turn possibly hinder peptide backbone flexibility has interesting implications in the conformational as well as the biological role of O-glycoproteins.     

Kinetic isotope effects reveal the presence of significant secondary structure in the transition state for the folding of the N-terminal domain of L9

Our present understanding of the nature of the transition state for protein folding depends predominantly on studies where individual side-chain contributions are mapped out by mutational analysis (phi value analysis). This approach, although extremely powerful, does not in general provide direct information about the formation of backbone hydrogen bonds. Here, we report the results of amide H/D isotope effect studies that probe the development of hydrogen bonded interactions in the transition state for the folding of a small alpha-beta protein, the N-terminal domain of L9. Replacement of amide protons by deuterons in a solvent of constant isotopic composition destabilized the domain, decreasing both its T(m) and Delta G(0) of unfolding. The folding rate also decreased. The parameter Phi(H/D), defined as the ratio of the effect of isotopic substitution upon the activation free energy to the equilibrium free energy was determined to be 0.6 in a D(2)O background and 0.75 in a H(2)O background, indicating that significant intraprotein hydrogen bond interactions are developed in the transition state for the folding of NTL9. The value is in remarkably good agreement with more traditional measures of the position of the transition state, which report on the relative burial of surface area. The results provide a picture of a compact folding transition state containing significant secondary structure. Indirect analysis argues that the bulk of the kinetic isotope effect arises from the beta-sheet-rich region of the protein, and suggests that the development of intraprotein hydrogen bonds in this region plays a critical role in the folding of NTL9.     

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

Carbon-13 NMR method for the detection of correlated hydrogen exchange at adjacent backbone peptide amides and its application to hydrogen exchange in five antiparallel beta strands within the hydrophobic core of Streptomyces subtilisin inhibitor (SSI)

A novel method for monitoring proton-deuteron (H/D) exchange at backbone amides is based on the observation of H/D isotope effects on the (13)C NMR signals from peptide carbonyls. The line shape of the carbonyl (13)C(i) signal is influenced by differential H/D occupancy at the two adjacent amides: the H(N)(i)(+1) (beta site) and the H(N)(i) (gamma site). At a carbon frequency of 75.4 MHz, the H --> D isotope shifts on the (13)C signal are about 5-7 Hz for exchange at the beta site and 2 Hz or less for exchange at the gamma site. Because the effects at the two sites are additive, the time dependence of the line shape of a particular carbonyl resonance can report not only the exchange rates at the individual sites but also the level of dual exchange. Therefore, the data can be analyzed to determine the rate (k(c)) and degree of correlated exchange (X(betagamma)) at the two sites. We have applied this approach to the investigation of the pH dependence of hydrogen exchange at several adjacent residues in Streptomyces subtilisin inhibitor (SSI). Two selectively labeled SSI proteins were produced: one with selective (13)C' labeling at all valyl residues and one with selective (13)C' labeling at all leucyl residues. This permitted the direct observation by one-dimensional (13)C NMR of selected carbonyl signals from residues with slowly exchanging amides at the i and i + 1 positions. The residues investigated were located in an alpha helix and in a five-stranded antiparallel beta sheet. Samples of the two labeled proteins were prepared at various pH values, and (13)C NMR spectra were collected at 50 degrees C prior to and at various times after transferring the sample from H(2)O to (2)H(2)O. Most of the slowly exchanging amides studied were intramolecular hydrogen-bond donors. In agreement with prior studies, the results indicated that the exchange rates of the amide hydrogens in proteins are governed not only by hydrogen bonding but also by other factors. For example, the amide hydrogen of Thr34 exchanges rapidly even though it is an intramolecular hydrogen-bond donor. Over nearly the whole pH range studied, the apparent rates of uncorrelated exchange (k(beta) and k(gamma)) were proportional to [OH(-)] and the apparent rates of correlated exchange at two adjacent sites (k(c)) were roughly proportional to [OH(-)](2). This enabled us to extract the pH-independent exchange rates (k(beta) degrees , k(gamma) degrees , and k(c) degrees ). In all cases in which correlated exchange could be measured, the observed sigmoidal pH dependence of X(betagamma) could be replicated roughly from the derived pH-independent rates.     

Deuterium isotope effects on <Superscript>15</Superscript>N backbone chemical shifts in proteins

Quantum mechanical calculations are presented that predict that one-bond deuterium isotope effects on the ¹⁵N chemical shift of backbone amides of proteins, ¹Δ¹⁵N(D), are sensitive to backbone conformation and hydrogen bonding. A quantitative empirical model for ¹Δ¹⁵N(D) including the backbone dihedral angles, Φ and Ψ, and the hydrogen bonding geometry is presented for glycine and amino acid residues with aliphatic side chains. The effect of hydrogen bonding is rationalized in part as an electric-field effect on the first derivative of the nuclear shielding with respect to N–H bond length. Another contributing factor is the effect of increased anharmonicity of the N–H stretching vibrational state upon hydrogen bonding, which results in an altered N–H/N–D equilibrium bond length ratio. The N–H stretching anharmonicity contribution falls off with the cosine of the N–H···O bond angle. For residues with uncharged side chains a very good prediction of isotope effects can be made. Thus, for proteins with known secondary structures, ¹Δ¹⁵N(D) can provide insights into hydrogen bonding geometries. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

C‐H…O hydrogen bonds in FK506‐binding protein–ligand interactions

Hydrogen bonds are important interaction forces observed in protein structures. They can be classified as stronger or weaker depending on their energy, thereby reflecting on the type of donor. The contribution of weak hydrogen bonds is deemed as an important factor toward structure stability along with the stronger bonds. One such bond, the C‐H…O type hydrogen bond, is shown to make a contribution in maintaining three dimensional structures of proteins. Apart from their presence within protein structures, the role of these bonds in protein–ligand interactions is also noteworthy. In this study, we present a statistical analysis on the presence of C‐H…O hydrogen bonds observed between FKBPs and their cognate ligands. The FK506‐binding proteins (FKBPs) carry peptidyl cis–trans isomerase activity apart from the immunosuppressive property by binding to the immunosuppressive drugs FK506 or rapamycin. Because the active site of FKBPs is lined up by many hydrophobic residues, we speculated that the prevalence of C‐H…O hydrogen bonds will be considerable. In a total of 25 structures analyzed, a higher frequency of C‐H…O hydrogen bonds is observed in comparison with the stronger hydrogen bonds. These C‐H…O hydrogen bonds are dominated by a highly conserved donor, the C^α/β of Val55 and an acceptor, the backbone oxygen of Glu54. Both these residues are positioned in the β4‐α1 loop, whereas the other residues Tyr26, Phe36 and Phe99 with higher frequencies are lined up at the opposite face of the active site. These preferences could be implicated in FKBP pharmacophore models toward enhancing the ligand affinity. This study could be a prelude to studying other proteins with hydrophobic pockets to gain better insights into ligand recognition. Copyright © 2013 John Wiley & Sons, Ltd.     

On the pH dependence of amide proton exchange rates in proteins.

We have analyzed the pH dependencies of published amide proton exchange rates (kex) in three proteins: bovine pancreatic trypsin inhibitor (BPTI), bull seminal plasma proteinase inhibitor IIA (BUSI IIA), and calbindin D9K. The base-catalyzed exchange rate constants (kOH) of solvent exposed amides in BPTI are lower for residues with low peptide carbonyl exposure, showing that the environment around the carbonyl oxygen influences kOH. We also examined the possible importance of an exchange mechanism that involves formations of imidic acid intermediates along chains of hydrogen-bonded peptides in the three proteins. By invoking this "relayed imidic acid exchange mechanism," which should be essentially acid-catalyzed, we can explain the surprisingly high pHmin (the pH value at which kex reaches a minimum) found for the non-hydrogen-bonded amide protons in the beta-sheet in BPTI. The successive increase of pHmin along a chain of hydrogen-bonded peptides from the free amide to the free carbonyl, observed in BPTI, can be explained as an increasing contribution of the proposed mechanism in this direction of the chain. For BUSI IIA (pH 4-5) and calbindin D9K (pH 6-7) the majority of amide protons with negative pH dependence of kex are located in chains of hydrogen-bonded peptides; this situation is shown to be consistent with the proposed mechanism.     

Pressure response of protein backbone structure. Pressure-induced amide 15N chemical shifts in BPTI.

The effect of pressure on amide 15N chemical shifts was studied in uniformly 15N-labeled basic pancreatic trypsin inhibitor (BPTI) in 90%1H2O/10%2H2O, pH 4.6, by 1H-15N heteronuclear correlation spectroscopy between 1 and 2,000 bar. Most 15N signals were low field shifted linearly and reversibly with pressure (0.468 +/- 0.285 ppm/2 kbar), indicating that the entire polypeptide backbone structure is sensitive to pressure. A significant variation of shifts among different amide groups (0-1.5 ppm/2 kbar) indicates a heterogeneous response throughout within the three-dimensional structure of the protein. A tendency toward low field shifts is correlated with a decrease in hydrogen bond distance on the order of 0.03 A/2 kbar for the bond between the amide nitrogen atom and the oxygen atom of either carbonyl or water. The variation of 15N shifts is considered to reflect site-specific changes in phi, psi angles. For beta-sheet residues, a decrease in psi angles by 1-2 degrees/2 kbar is estimated. On average, shifts are larger for helical and loop regions (0.553 +/- 0.343 and 0.519 +/- 0.261 ppm/2 kbar, respectively) than for beta-sheet (0.295 +/- 0.195 ppm/2 kbar), suggesting that the pressure-induced structural changes (local compressibilities) are larger in helical and loop regions than in beta-sheet. Because compressibility is correlated with volume fluctuation, the result is taken to indicate that the volume fluctuation is larger in helical and loop regions than in beta-sheet. An important aspect of the volume fluctuation inferred from pressure shifts is that they include motions in slower time ranges (less than milliseconds) in which many biological processes may take place.     

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.     

Trans-hydrogen bond deuterium isotope effects of A:T base pairs in DNA

The chemical shifts of (13)C2 of adenosine residues of DNA were observed to experience a through-space or trans-hydrogen bond isotope effect as a result of deuterium substitution at the imino hydrogen site of base-paired thymidine residues. NMR measurements of several self-complementary DNA duplexes at natural abundance (13)C in 50% H(2)O, 50% D(2)O solvent mixtures yielded an average trans-hydrogen bond isotope effect, (2h)Delta(13)C2, of -47 ppb. The data suggest that stronger hydrogen bonds have more negative (2h)Delta(13)C2 values, which means that A:T N1.H3 hydrogen bonds increase the anharmonicity of the effective vibrational potential of H3. However, (2h)Delta(13)C2 values do not correlate with intra-residue (2)Delta(13)C4 values of thymidine observed here and earlier (Vakonakis et al., 2003), which suggests that (2h)Delta(13)C2 is not determined entirely by hydrogen bond strength. Instead, the variations observed in (2h)Delta(13)C2 values suggest that they may also be sensitive to base pair geometry.     

The hydration of amides in helices; a comprehensive picture from molecular dynamics,IR, and NMR

We examined the hydration of amides of alpha(3)D, a simple, designed three-helix bundle protein. Molecular dynamics calculations show that the amide carbonyls on the surface of the protein tilt away from the helical axis to interact with solvent water, resulting in a lengthening of the hydrogen bonds on this face of the helix. Water molecules are bonded to these carbonyl groups with partial occupancy ( approximately 50%-70%), and their interaction geometries show a large variation in their hydrogen bond lengths and angles on the nsec time scale. This heterogeneity is reflected in the carbonyl stretching vibration (amide I' band) of a group of surface Ala residues. The surface-exposed amides are broad, and shift to lower frequency (reflecting strengthening of the hydrogen bonds) as the temperature is decreased. By contrast, the amide I' bands of the buried (13)C-labeled Leu residues are significantly sharper and their frequencies are consistent with the formation of strong hydrogen bonds, independent of temperature. The rates of hydrogen-deuterium exchange and the proton NMR chemical shifts of the helical amide groups also depend on environment. The partial occupancy of the hydration sites on the surface of helices suggests that the interaction is relatively weak, on the order of thermal energy at room temperature. One unexpected feature that emerged from the dynamics calculations was that a Thr side chain subtly disrupted the helical geometry 4-7 residues N-terminal in sequence, which was reflected in the proton chemical shifts and the rates of amide proton exchange for several amides that engage in a mixed 3(10)/alpha/pi-helical conformation.     

Synthesis of 17O isotope labeled Leu-enkephalin and 17O n.m.r

Oxygen-17 isotope was introduced into the alpha-carboxyl group of glycine, 1-phenylalanine, 1-leucine and 1-tyrosine by acid catalyzed exchange of 17O from H2O(17) or by acid hydrolysis of respective amino acid methyl esters in H2O(17). Quantitative enrichment of glycine was achieved by acid hydrolysis of amino acetonitrile in H2O(17). For alpha-amino protection in amino acids t-butoxycarbonyl (Boc) group was employed for 17O labeled enkephalin synthesis. Five analogues of Leu-enkephalins (I-V) labeled with 17O at different amino acid residues were synthesized by solid phase method. 17O n.m.r. spectra were measured at 24.4 and 67.8 MHz for Leu-enkephalins 17O labeled at Gly2 and Phe4 positions. A downfield shift was observed for 17O labeled Gly2 Leu-enkephalin upon heating. This shift is indicative of the rupture of intramolecular hydrogen bonds. The preliminary results confirm the hypothesis that an intramolecular hydrogen bond exists between the carbonyl group of Gly2 and NH group of Leu5.     

Conformational preferences and prolyl cis–trans isomerization of phosphorylated Ser/Thr‐Pro motifs

The conformational study on Ac‐pSer‐Pro‐NHMe and Ac‐pThr‐Pro‐NHMe peptides has been carried out using hybrid density functional methods with the implicit solvation reaction field theory at the B3LYP/ 6‐311++G(d,p)//B3LYP/6‐31+G(d) level of theory in the gas phase and in solution (chloroform and water). For both pSer‐Pro and pThr‐Pro peptides in the gas phase and in chloroform, the most preferred conformation has the α‐helical structure for the pSer/pThr residue, the down‐puckered polyproline I structure for the Pro residue, and the cis prolyl peptide bond between the two residues, in which two hydrogen bonds between the phosphate oxygens with the backbone N? H groups seem to play a role. However, the trans conformations that have a single hydrogen bond of the phosphate oxygen with either of two backbone N? H groups become most preferred for both peptides in water. This is because the hydration free energy of the anionic oxygen of the phosphate group is expected to dramatically decrease for the cis conformation upon formation of the hydrogen bond with the backbone N? H groups. These calculated results are consistent with the observations by NMR and IR experiments, suggesting the existence of hydrogen bonds between the charged phosphoryl group and the backbone amide protons in solution. The calculated cis populations of 14.7 and 14.2% and rotational barriers of 19.87 and 20.57 kcal/mol to the cis‐to‐trans isomerization for pSer‐Pro and pThr‐Pro peptides in water, respectively, are consistent with the observed values for pSer‐Pro and pThr‐Pro containing peptides from NMR experiments. However, the hydrogen bond between the prolyl nitrogen and the following amide N? H group, which was suggested to be capable of catalyzing the prolyl isomerization, does not play a role in stabilizing the preferred transition state for the pSer/pThr‐Pro peptides in water. Instead, the amide hydrogen of the NHMe group is involved in a bifurcated hydrogen bond with the anionic oxygen and phosphoester oxygen of the phosphate group. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 330–339, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.     

Solvated helical backbones: x-ray diffraction study of Boc-Ala-Leu-Aib-Ala-Leu-Aib-OMe.H2O

A second example of insertion of a water molecule into the helical backbone of an apolar peptide is presented here and compared to a similar occurrence in a longer peptide with the same type of sequence of residues, i.e., Boc-Aib-(Ala-Leu-Aib)3-OMe. The backbone of the title compound assumes an approximate 3(10)-helical form with three 4----1 hydrogen bonds. In the place of a fourth 4----1 hydrogen bond, a water molecule is inserted between O(1) and N(4), and acts as a bridge by forming hydrogen bonds N(4) ... W(1) (2.95 A) and W(1) ... O(1) (2.81 A). The water molecule participates in a third hydrogen bond with a neighboring peptide molecule, W(1) ... O(4) (2.91 A). The insertion of the water molecule causes the apolar peptide to mimic an amphiphilic helix. Crystals grown from ethyl acetate/petroleum ether (reported here) or from methanol/water solution are in space group P2(1)2(1)2(1) with a = 12.024(4) A, b = 15.714(6) A, c = 21.411(7) A, Z = 4 and dcalc = 1.124 g/cm3 for C32H58N6O9.H2O. The overall agreement factor R is 6.3% for 2707 reflections observed with intensities greater than 3 sigma(F) and the resolution is 0.90 A.