CH...O hydrogen bonds at protein-protein interfaces

For the first time, a statistical potential has been developed to quantitatively describe the CH.O hydrogen bonding interaction at the protein-protein interface. The calculated energies of the CH.O pair interaction show a favorable valley at approximately 3.3 A, exhibiting a feature typical of an H-bond and similar to the ab initio quantum calculation result (Scheiner, S., Kar, T., and Gu, Y. (2001) J. Biol. Chem. 276, 9832-9837). The potentials have been applied to a set of 469 protein-protein complexes to calculate the contribution of different types of interactions to each protein complex: the average energy contribution of a conventional H-bond is approximately 30%; that of a CH.O H-bond is 17%; and that of a hydrophobic interaction is 50%. In some protein-protein complexes, the contribution of the CH.O H-bond can reach as high as approximately 40-50%, indicating the importance of the CH.O H-bond at the protein interface. At the interfaces of these complexes, C(alpha)H.O H-bonds frequently occur between adjacent strands in both parallel and antiparallel orientations, having the obvious structural motif of bifurcated H-bonds. Our study suggests that the weak CH.O H-bond makes an important contribution to the association and stability of protein complexes and needs more attention in protein-protein interaction studies.     

Twist and shear in beta-sheets and beta-ribbons

The structures of the beta-sheets and the beta-ribbons have been analysed using high-resolution protein structure data. Systematic asymmetries measured in both parallel and antiparallel beta-structures include the sheet twist and the strand shear. In order to determine the origin of these asymmetries, numerous interactions and correlations were examined. The strongest correlations are observed for residues in antiparallel beta-sheets and beta-ribbons that form non-H-bonded pairs. For these residues, the sheet twist is correlated to the backbone phi angle but not to the psi angle. Our analysis supports the existence of an inter-strand C(alpha)H(alpha)...O weak H-bond, which, together with the CO...HN H-bond, constitutes a bifurcated H-bond that links neighbouring beta-strands. Residues of beta-sheets and beta-ribbons in high-resolution protein structures form a distinct region of the Ramachandran plot, which is determined by the formation of the bifurcated H-bond, the formation of an intra-strand O...H(alpha) non-bonded polar interaction, and an intra-strand O...C(beta) steric clash. Using beta-strands parameterised by phi-psi values from the allowed beta-sheet region of the Ramachandran plot, the shear and the right-hand twist can be reproduced in a simple model of the antiparallel and parallel beta-ribbon that models the bifurcated H-bonds specifically. The conformations of interior residues of beta-sheets are shown to be subsets of the conformations of residues of beta-ribbons.     

Pex, analytical tools for PDB files. II. H-Pex: noncanonical H-bonds in alpha-helices

We use the H-Pex (Thomas et al., this issue) to analyze the main chain interactions in 131 proteins. In antiparallel beta-sheets, the geometry of the N...O bond is: median N...O distances, 2.9 SA, C==O...N angles at 154 degrees and the C alpha--C==O...H angles are dispersed around 3 degrees. In some instances, the other side of the C==O axis is occupied by a HC alpha. As recently supported by Vargas et al. (J Am Chem Soc 2000;122:4750-4755) C alpha H...O and NH...O could cooperate to sheet stability. In alpha-helices, the main chain C==O interact with the NH of their n + 4 neighbor on one side, and with a C beta H or C gamma H on the other side. The median O...N distance (3.0 A) and C==N angle (147 degrees) suggest a canonical H-bond, but the C alpha--C==O...H dihedral angle invalidates this option, since the hydrogen attacks the oxygen at 122 degrees, i.e., between the sp(2) and pi orbitals. This supports that the H-bond is noncanonical. In many instances, the C gamma H or the C beta H of the n + 4 residue stands opposite to the NH with respect to the oxygen. Therefore, we propose that, in alpha-helices, the C gamma H or C beta H and the NH of the n + 4 residue hold the oxygen like an electrostatic pincher. Proteins 2001;43:37-44.     

Conformational analysis of biologically active thyroliberin analogs by two-dimensional NMR spectroscopy]

Preferable conformations of thyrotropin-releasing hormone (TRH, Glp-His-Pro-NH2) and its analogues Glp-Glu(R)-Pro-NH2 (R = NHCH(CH3)CH2Ar), Glp-Gln-Abu-NH2, Dho-Gln-Abu-NH2 in DMSO solution are determined using two-dimensional 1H NMR spectroscopy (delta-J-correlated, COSY and NOESY). Torsion angles psi i and chi i for every amino acid were calculated on the basis of the spin-spin coupling constants 3JNH-C alpha H and 3JC alpha H-C beta H values. The NOESY data were used for selecting the peptide conformations realized in solution. Distances between protons interacting by the dipole mechanism (d-contacts) were calculated using NOE values. These experiments allow one to estimate the torsion angles psi (between C alpha H-CO). TRH has an intramolecular H-bond between NH2-protons and His carbonyl with the torsion angles omega 3 = 180 degrees and psi 3 = 0 degrees. It is formation of this H-bond that apparently promotes the domination of the trans configuration of the His-Pro peptide bond. An intramolecular NH2-C alpha CO (Glp) H-bonding is revealed in other investigated compounds. It is known that a similar conformation of the TRH is realized in the course of its interaction with receptor.     

Demonstration of a tertiary interaction in solution between the extra arm and the D-stem in two different transfer RNA''s by NMR.

According to the X-ray structure of yeast tRNAPhe at 2.5 A resolution, a hydrogen bond is formed between m7G46 and G22. By removal of this m7G46-residue we demonstrate that this interaction is present in solution as well. Comparison of the 1H 360 MHz NMR spectra of intact yeast tRNAPhe and its m7G-excised derivative locates the position of this tertiary H-bond at 12.5 ppm downfield from DSS. Additional evidence for the presence of this interaction in solution comes from a comparison of 1H NMR spectra of E. coli tRNAf1Met and E. coli tRNAf3Met, which differ only in a single position in the extra arm. In tRNAf3Met residue 47 is a m7G-residue, whereas in tRNAf3Met it is A, resulting in the absence of the m7G47 - G23 - C13 triple interaction, characteristic of tRNAf1Met. The resonance position of this tertiary interaction in tRNAf1Met is located around -13.6 ppm, a chemical shift difference of 1.1 ppm with respect to the position observed for tRNAPhe. The origin of this chemical shift difference is discussed in relation to the structure of their respective augmented D-helices.     

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.     

Peptides containing the sulfonamide junction. 2. Structure and conformation of Z-Tau-Pro-D-Phe-NHiPr

The taurine (Tau) containing N-protected pseudotripeptide isopropylamide Z-Tau-Pro-D-Phe-NHiPr (1) has been specifically designed and synthesized as suitable model to test the ability of the sulfonamido group to participate as H-bond acceptor to a type II beta-turn and to get information on the preferred rotameric conformation around the S-N bond and the hybridization state of the nitrogen atom. The present structural investigation reveals that, although the sulfonamide junction is invariably folded in a gauche mode, the beta-turn structure, stabilized by the 4 --> 1 hydrogen bond, is not found in the crystal and the sulfonamido oxygen atoms are not involved in any intra- or intermolecular hydrogen-bond interaction. More than one conformer populates the CDCl(3) solution with only a minor contribution by the expected beta-turn. The Pro nitrogen is significantly pyramidalized and the nitrogen lone pair points in opposite direction to that of the Pro C(alpha)H bond thus adopting R chirality, in an arrangement practically identical to that found in the previously studied homochiral analogue Z-Tau-Pro-Phe-NHiPr.     

Contribution of hydrogen bonding to protein stability estimated from isotope effects

An unresolved issue in structural biology concerns the relative contribution of H bonds to protein stability. We use the small molecules 4-acetamidobenzoic acid and N-acetylanthranilic acid as model compounds to relate the energetic contribution from hydrogen bonds (H bonds) to the deuterium/hydrogen amide isotope effect. N-Acetylanthranilic acid models carbonyl-amide H bonds formed during protein folding; 4-acetamidobenzoic acid models the unfolded state in which the amide H bonds to water. NMR is used to measure shifts in the pK(a) of the ionizable carboxyl group when the amides of the compounds are either protonated or deuterated. From the pK(a) shift, we obtain a quantitative scale factor: SF = partial partial differential(DeltaG(HB))/partial partial differential(RT ln Phi), where DeltaG(HB) is the change in free energy of an H bond upon isotope substitution and Phi is the fractionation factor. Isotope effect data also are reported for a small globular protein, lambda repressor, using the "C(m) experiment". The protein's isotope effect, which reports on the shape of the energy well, is converted to H-bonding free energy by applying the scale factor. We estimate that amide-related H bonds (amide-carbonyl and amide-water) contribute favorably to protein stability by approximately 30-50 kcal/mol in lambda repressor, GCN4 coiled coil, and cytochrome c but unfavorably by approximately 6 kcal/mol in ubiquitin. The results indicate that H-bond strength varies from one protein to another and presumably at different sites within the same protein.     

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.     

An additional H-bond in the alpha 1 beta 2 interface as the structural basis for the low oxygen affinity and high cooperativity of a novel recombinant hemoglobin (beta L105W)

Site-directed mutagenesis has been used to construct three recombinant mutant hemoglobins (rHbs), rHb(beta L105W), rHb(alpha D94A/betaL105W), and rHb(alpha D94A). rHb(beta L105W) is designed to form a new hydrogen bond from beta 105Trp to alpha 94Asp in the alpha(1)beta(2) subunit interface to lower the oxygen binding affinity by stabilizing the deoxy quaternary structure. We have found that rHb(beta L105W) does indeed possess a very low oxygen affinity and maintains normal cooperativity (P(50) = 28.2 mmHg, n(max) = 2.6 in 0.1 M sodium phosphate at pH 7.4) compared to those of Hb A (P(50) = 9.9 mmHg, n(max) = 3.2 at pH 7.4). rHb(alpha D94A/beta L105W) and rHb(alpha D94A) are expressed to provide evidence that rHb(betaL 105W) does form a new H-bond from beta 105Trp to alpha 94Asp in the alpha(1)beta(2) subunit interface of the deoxy quaternary structure. Our multinuclear, multidimensional nuclear magnetic resonance (NMR) studies on (15)N-labeled rHb(beta L105W) have identified the indole nitrogen-attached (1)H resonance of beta 105Trp for rHb(beta L105W). (1)H NMR studies on Hb A and mutant rHbs have been used to investigate the structural basis for the low O(2) affinity of rHb(beta L105W). Our NMR results provide evidence that rHb(beta L105W) forms a new H-bond from beta 105Trp to alpha 94Asp in the alpha(1)beta(2) subunit interface of the deoxy quaternary structure. The NMR results also show that these three rHbs can switch from the R quaternary structure to the T quaternary structure in their ligated state upon addition of an allosteric effector, inositol hexaphosphate. We propose that the low O(2) affinity of rHb(beta L105W) is due to the formation of a new H-bond between alpha 105Trp and alpha 94Asp in the deoxy quaternary structure.     

Quantification of H/D Isotope Effects on Protein Hydrogen-bonds by h3JNC′ and 1JNC′ Couplings and Peptide Group 15N and 13C′ Chemical Shifts

The effect of hydrogen/deuterium exchange on protein hydrogen bond coupling constants (h3)J(NC') has been investigated in the small globular protein ubiquitin. The couplings across deuterated or protonated hydrogen bonds were measured by a long-range quantitative HA(CACO)NCO experiment. The analysis is combined with a determination of the H(N)/D(N) isotope effect on the amide group (1)J(NC') couplings and the (15)N and (13)C' chemical shifts. On average, H-bond deuteration exchange weakens (h3)J(NC') and strengthens (1)J(NC') couplings. A correlation is found between the size of the (15)N isotope shift, the (15)N chemical shift, and the (h3)J(NC') coupling constants. The data are consistent with a reduction of donor-acceptor overlap as expected from the classical Ubbelohde effect and the common understanding that H(N)/D(N) exchange leads to a shortening of the N-hydron bond length.     

Direct detection of N-H[...]O=C hydrogen bonds in biomolecules by NMR spectroscopy

A nuclear magnetic resonance (NMR) experiment is described for the direct detection of N-H[...]O=C hydrogen bonds (H-bonds) in 15N and 13C isotope-labeled biomolecules. This quantitative 'long-range' HNCO-COSY (correlation spectroscopy) experiment detects and quantifies electron-mediated scalar couplings across the H-bond (H-bond scalar couplings), which connect the magnetically active (15)N and (13)C nuclei on both sides of the H-bond. Detectable H-bonds comprise the canonical backbone H-bonds in proteins as well as other H-bonds in proteins and nucleic acids with N-H donors and O=C (carbonylic or carboxylic) acceptors. Unlike other NMR observables, which provide only indirect evidence of the presence of H-bonds, the H-bond scalar couplings identify all partners of the H-bond, the donor, the donor proton and the acceptor, in a single experiment. The size of the scalar couplings can be related to H-bond geometries. The time required to detect the N-H[...]O=C H-bonds in small proteins (< or = approximately 10 kDa) is typically on the order of 1 d at millimolar concentrations, whereas H-bond detection for larger proteins (< or = approximately 30 kDa) may be possible within several days depending on concentration, isotope composition, magnetic field strength and molecular weight. The proteins ubiquitin (8.6 kDa), dimeric RANTES (2 x 8.5 kDa) and MAP30 (30 kDa) are used as examples to illustrate this procedure.     

An unusual C-H...O hydrogen bond mediated reversal of polypeptide chain direction in a synthetic peptide helix

An unusual C-terminal conformation has been detected in a synthetic decapeptide designed to analyze the stereochemistry of helix termination in polypeptides. The crystal structure of the decapeptide Boc-Leu-Aib-Val-Ala-Leu-Aib-Val-(D)Ala-(D)Leu-Aib-OMe reveals a helical segment spanning residues 1-7 and helix termination by formation of a Schellman motif, generated by (D)Ala(8) adopting the left-handed helical (alpha(L)) conformation. The extended conformation at (D)Leu(9) results in a compact folded structure, stabilized by a potentially strong C-H. O hydrogen bond between Ala(4) C(alpha)H and (D)Leu(9) CO. The parameters for C-H. O interaction are Ala(4) C(alpha)H. O=C (D)Leu(9) distance 3.27 A, C(alpha)-H. O angle 176 degrees, and O. H(alpha) distance 2.29 A. This structure suggests that insertion of contiguous D-residues may provide a handle for the generation of designed structures containing more than one helical segment folded in a compact manner.     

A 1H and 13C NMR study on the role of salt-bridges in the formation of a type I beta-turn in N-acetyl-L-Asp-L-Glu-L-Lys-L-Ser-NH2

The conformation of the tetrapeptide N-Acetyl-Asp7-Glu8-Lys9-Ser10-NH2, a fragment of the type I collagen alpha-1 chain N-telopeptide, has been studied by 1H and 13C NMR and circular dichroism spectroscopy. The spectroscopic evidence, based on two-dimensional, phase-sensitive NMR techniques such as COSY, ROESY, proton-carbon shift correlation and selective COLOC, indicates a strong dependence of the conformation on the experimental conditions. In CD3OH/H2O (60/40) at ca. neutral pH the tetrapeptide forms a beta-turn, stabilized by a hydrogen bond between NH(S10) and CO(D7) and a strong salt-bridge between COO-(E8) and NH3+(K9). The beta-turn is type I and appears to coexist with a non-hydrogen-bonded structure. The coexistence of these two conformers is proven by proton NMR data such as NH-NH ROEs, reduced NH-H alpha (E8) coupling constant, NH(E8) low-field shift and the temperature coefficient of NH(S10), whereas the conclusion regarding the salt-bridge is based on 13C results. In the same solvent, at a pH below the pKa of the carboxyl groups, no evidence for a conformation other than extended can be found. In aqueous solution at approximately neutral pH, evidence for the E8-K9 charge interaction is observed, but not for a hydrogen bond anywhere in the molecule.     

Quantitative comparison of the hydrogen bond network of A-state and native ubiquitin by hydrogen bond scalar couplings

The backbone hydrogen bond (H-bond) network of the partially folded A-state of ubiquitin (60% methanol, 40% water, pH 2) has been characterized quantitatively by (h3)J(NC)(') H-bond scalar couplings between the (15)N nuclei of amino acid H-bond donors and the (13)C carbonyl nuclei of the acceptors. Results on (h3)J(NC)(') couplings and the amide proton ((1)H(N)) chemical shifts for the A-state are compared quantitatively to the native state. The (h3)J(NC)(') correlations of the A-state show intact, nativelike H-bonds of the first beta-hairpin beta1/beta2 and the alpha-helix, albeit at lower strength, whereas the H-bonds in the C-terminal part change from a pure beta-structure to an all alpha-helical H(N)(i)-->O(i-4) connectivity pattern. A residue-specific analysis reveals that the conformations within the conserved secondary structure segments are much more homogeneous in the A-state than in the native state. Thus, the strong asymmetry of (h3)J(NC)(') couplings and (1)H(N) chemical shifts between the interior and exterior sides of the native state alpha-helix vanishes in the A-state. This indicates that the bend of this helix around the native state hydrophobic core is released in the homogeneous solvent environment of the A-state. Similarly, an irregularity in the behavior of H-bond I3-->L15 in hairpin beta1/beta2, which results from strong contacts to strand beta5 in the native state, is absent in the A-state. These findings rationalize the behavior of the (1)H(N) chemical shifts in both states and indicate that the A-state is in many aspects similar to the onset of thermal denaturation of the native state.     

DFT studies on the intrinsic conformational properties of non-ionic pyrrolysine in gas phase

B3LYP/6-31G(d,p) level of theory is used to carry out a detailed gas phase conformational analysis of non-ionized (neutral) pyrrolysine molecule about its nine internal back-bone torsional angles. A total of 13 minima are detected from potential energy surface exploration corresponding to the nine internal back-bone torsional angles. These minima are then subjected to full geometry optimization and vibrational frequency calculations at B3LYP/6-31++G(d,p) level. Characteristic intramolecular hydrogen bonds present in each conformer, their relative energies, theoretically predicted vibrational spectra, rotational constants and dipole moments are systematically reported. Single point calculations are carried out at B3LYP/6-311++G(d,p) and MP2/6-31++G(d,p) levels. Six types of intramolecular H-bonds, viz. O…H–O, N…H-O, O…H–N, N…H–N, O…H–C and N…H–C, are found to exist in the pyrrolysine conformers; all of which contribute to the stability of the conformers. The vibrational frequencies are found to shift invariably toward the lower side of frequency scale corresponding to the presence of intramolecular H-bond interactions in the conformers.     

MP2, density functional theory, and molecular mechanical calculations of C-H···π and hydrogen bond interactions in a cellulose-binding module-cellulose model system

Exploring non-covalent interactions, such as C-H···π stacking and classical hydrogen bonding (H-bonding), between carbohydrates and carbohydrate-binding modules (CBMs) is an important task in glycobiology. The present study focuses on intermolecular interactions, such as C-H?π (sugar-aromatic stacking) and H-bonds, between methyl β-d-glucopyranoside and l-tyrosine—a proxy model system for a cellulose-CBM complex. This work has made use of various types of quantum mechanics (QM) and molecular mechanics (MM) methods to determine which is the most accurate and computationally efficient. The calculated interaction potential energies ranged between −24 and −38 kJ/mol. The larger interaction energy is due to H-bonding between the phenyl hydroxyl of tyrosine and the O4 of the sugar. Density functional theory (DFT) methods, such as BHandHLYP and B3LYP, exaggerate the H-bond. Although one of the MM methods (viz. MM+) considered in this study does maintain the C-H?π stacking configuration, it underestimates the interaction energy due to the loss of the H-bond. When the O-H bond vector is in the vicinity of O4 (O-H?O4 ≈ 2 Å, e.g., in the case of MP2/6-31G(d)), the torsional energy drops to a minimum. For this configuration, natural bond orbital (NBO) analysis also supports the presence of this H-bond which arises due to orbital interaction between one lone pair of the sugar O4 and the σ∗(O-H) orbital of the phenyl group of tyrosine. The stabilization energy due to orbital delocalization of the H-bonded system is ∼13 kJ/mol. This H-bond interaction plays an important role in controlling the CH/π interaction geometry. Therefore, the C-H?π dispersive interaction is the secondary force, which supports the stabilization of the complex. The meta-hybrid DFT method, M05-2X, with the 6-311++G(d,p) basis set agrees well with the MP2 results and is less computationally expensive. However, the M05-2X method is strongly basis set dependent in describing this CH/π interaction. Computed IR spectra with the MP2/6-31G(d) method show blue shifts for C1-H, C3-H, and C5-H stretching frequencies due to the C-H?π interaction. However, the M05-2X/6-311++G(d,p) method shows a small red shift for the C1-H stretching region and blue shifts for the C2-H and C3-H stretches. For the aromatic tyrosine C_δ1-C_ε1 and C_δ2-C_ε2 bonds in the complex, the calculated IR spectra show red shifts of 12 cm⁻¹ (MP2/6-31G(d)) and 5 cm⁻¹ (M05-2X/6-311++G(d,p)). This study also reports the upfield shifts of computed ¹H NMR chemical shifts due to the C-H?π interaction.     

The influence of Li+, Na+, Mg2+, Ca2+, and Zn2+ ions on the hydrogen bonds of the Watson-Crick base pairs

The interaction of mono- and divalent metal ions with the nucleic acid base pairs A:T and G:C has been studied using ab initio self-consistent field Hartree-Fock computations with minimal basis sets. Energy-optimized structures of the two base pairs with a final base-base distance of L = 10.35 A have been determined and were further used in calculations on ternary complexes Mn+ - A:B together with previously computed coordination geometries of the cations at adenine (Ade), thymine (Thy), and guanine (Gua). Besides the binding energy of the various metal ions to the base pairs, changes in the stability of the H bonds between Ade and Thy or Gua and Cyt have been determined. Polarization effects of the metal ion on the ligand turned out to increase the binding between complementary bases. Regardless of the metal species, cation binding to Gua N(3) and Thy O(2) leads to a special increase in H-bond stability, whereas binding to Ade N(3) changes the H-bond stability least. Situated in between are the stabilizing effects caused by Gua and Ade N(7) coordination. A remarkable relation between the stability of the H bond and the distance from metal binding site to H bonds was found. This relationship has been rationalized in terms of partial charges of the atoms participating in H bonding, which can reveal the trend in the electrostatic part of total H bond energy. It can be shown that a short distance between coordination site and acceptor hydrogen increases the H-bond strength substantially, while a long distance shows minor effects as supposed. On the other hand, the opposite effect is observed for the influence of the distance between binding site and donor atom. A comparison of our findings with a new model of transition metal ion facilitated rewinding of denatured DNA proposed by S. Miller, D. VanDerveer, and L. Marzilli is given [(1985) J. Am. Chem. Soc. 107, 1048-1055].     

A C-H triplebond O hydrogen bond stabilized polypeptide chain reversal motif at the C terminus of helices in proteins

The serendipitous observation of a C-H cdots, three dots, centered O hydrogen bond mediated polypeptide chain reversal in synthetic peptide helices has led to a search for the occurrence of a similar motif in protein structures. From a dataset of 634 proteins, 1304 helices terminating in a Schellman motif have been examined. The C-H triplebond O interaction between the T-4 C(alpha)H and T+1 Cz doublebond O group (C triplebond O< or =3.5A) becomes possible only when the T+1 residue adopts an extended beta conformation (T is defined as the helix terminating residue adopting an alpha(L) conformation). In all, 111 examples of this chain reversal motif have been identified and the compositional and conformational preferences at positions T-4, T, and T+1 determined. A marked preference for residues like Ser, Glu and Gln is observed at T-4 position with the motif being further stabilized by the formation of a side-chain-backbone O triplebond H-N hydrogen bond involving the side-chain of residue T-4 and the N-H group of residue T+3. In as many as 57 examples, the segment following the helix was extended with three to four successive residues in beta conformation. In a majority of these cases, the succeeding beta strand lies approximately antiparallel with the helix, suggesting that the backbone C-H triplebond O interactions may provide a means of registering helices and strands in an antiparallel orientation. Two examples were identified in which extended registry was detected with two sets of C-H cdots, three dots, centered O hydrogen bonds between (T-4) C(alpha)H triplebond O (T+1) and (T-8) C(alpha)H triplebondC doublebond O (T+3).     

Keratan sulfates from bovine tracheal cartilage structural studies of intact polymer chains using H and 13C NMR spectroscopy.

Intact keratan sulfate chains derived from bovine tracheal cartilage have been examined using both one-dimensional methods and the two-dimensional experiments COSY-45 and TOCSY for homonuclear shift correlations and a modified COLOC (correlated spectroscopy for long-range couplings) approach for 13C-1H shift correlations. Partial 1H and 13C NMR signal assignments for residues within the intact polymer chain are reported; data derived from the repeat region signals and from chain cap residues are assigned by comparison with published data derived from oligosaccharides obtained through cleavage of keratan sulfate polymer chains using keratanase and keratanase II and are discussed in detail. The one-dimensional spectra for both 1H and 13C nuclei contain highly crowded signal clusters for which data analysis is not directly possible. COSY-45 analysis allow the correlation and assignment of many proton resonances located within the 3.4-4.8 p.p.m. chemical shift region while from the C/H correlation spectrum data are assignable for some signals within the complex set of carbon resonances which fall in the region between 68 and 86 p.p.m., This work using material from tracheal cartilage has permitted the first detailed combined 1H and 13C NMR examination of the primary keratan sulfate polymer structure; this sequence forms the basis for the more complex members of the keratan sulfate family present in other tissues such as articular cartilage and cornea where further residues such as (alpha1-3)-linked fucose and (alpha2-6)-linked N-acetylneuraminic acid are also present. This nondestructive method of analysis complements the currently available degradative methods for structure determination which may then subsequently be utilized.