Occurrence of bifurcated three-center hydrogen bonds in proteins

Analysis of 13 high-resolution protein X-ray crystal structures shows that 1204 (24%) of all the 4974 hydrogen bonds are of the bifurcated three-center type with the donor X-H opposing two acceptors A1, A2. They occur systematically in alpha-helices where 90% of the hydrogen bonds are of this type; the major component is (n + 4)N-H ... O = C(n) as expected for a 3.6(13) alpha-helix, and the minor component is (n + 4)N-H ... O = C(n + 1), as observed in 3(10) helices; distortions at the C-termini of alpha-helices are stabilized by three-center bonds. In beta-sheets 40% of the hydrogen bonds are three-centered. The frequent occurrence of three-center hydrogen bonds suggests that they should not be neglected in protein structural studies.     

Double and bifurcated hydrogen bonds in alpha-helices of globular proteins]

The statistical analysis of hydrogen bonds distribution in space structures of globular proteins has been done. The parameters of H-bonds in the different secondary structures of globular proteins were collected. In alpha-helices besides the canonical 1-5 H-bonds (the mean length 3 A), 1-4 H-bonds were observed (the mean length 3.2 A). The histograms of length and angular distributions of the bonds are presented. It was found on the basis of quantum chemistry calculations that most H-bonds in alpha-helices are double or bifurcated.     

DNA curvature does not require bifurcated hydrogen bonds or pyrimidine methyl groups.

Short tracts of the homopolymer dA.dT confer intrinsic curvature on the axis of the DNA double helix. This phenomenon is assumed to be a consequence of such tracts adopting a stable B'-DNA conformation that is distinct from B-form structure normally assumed by other DNA sequences. The more stable B' structure of dA.dT tracts has been attributed to several possible stabilizing factors: (1) optimal base stacking interactions consequent upon the high propeller twist, (2) bifurcated hydrogen bonds between adjacent dA.dT base-pairs, (3) stacking interactions involving the dT methyl groups, and finally (4) a putative spine of ordered water molecules in the minor groove. DNA oligodeoxynucleotides have been synthesized that enable these hypotheses to be tested; of particular interest is the combination of effects due to bifurcation (2) and methylation of the pyrimidines nucleotides (3). The data indicate that neither bifurcated hydrogen bonds nor pyrimidine methyl groups nor both are essential for DNA curvature. The data further suggest that the influence of the minor groove spine of hydration on the B'-formation is small. The experiments favor the hypothesis that base stacking interactions are the dominant force in stabilizing the B'-form structure.     

Observation of a unique pattern of bifurcated hydrogen bonds in the crystal structures of the N-glycoprotein linkage region models

Elucidation of the intra- and intermolecular carbohydrate-protein interactions would greatly contribute toward obtaining a better understanding of the structure-function correlations of the protein-linked glycans. The weak interactions involving C-H...O have recently been attracting immense attention in the domain of biomolecular recognition. However, there has been no report so far on the occurrence of C-H...O hydrogen bonds in the crystal structures of models and analogs of N-glycoproteins. We present herein an analysis of C-H...O interactions in the crystal structures of all N-glycoprotein linkage region models and analogs. The study reveals a cooperative network of bifurcated hydrogen bonds consisting of N-H...O and C-H...O interactions seen uniquely for the models. The cooperative network consists of two antiparallel chains of bifurcated hydrogen bonds, one involving N1-H, C2'-H and O1' of the aglycon moiety and the other involving N2-H, C1-H and O1' of the sugar. Such bifurcated hydrogen bonds between the core glycan and protein are likely to play an important role in the folding and stabilization of proteins.     

Desolvation shell of hydrogen bonds in folded proteins,protein complexes and folding pathways

A few backbone hydrogen bonds (HBS) in native protein folds are poorly protected from water attack: their desolvation shell contains an inordinately low number of hydrophobic residues. Thus, an approach by solvent-structuring moieties of a binding partner should contribute significantly to enhance their stability. This effect represents an important factor in the site specificity inherent to protein binding, as inferred from a strong correlation between poorly desolvated HBs and binding sites. The desolvation shells were also examined in a dynamic context: except for a few singular under-protected bonds, the size of desolvation shells is preserved along the folding trajectory.     

Close pairs of carboxylates: a possibility of multicenter hydrogen bonds in proteins

Covalent attachment of hydrogen to the donor atom may be not an essential characteristic of stable hydrogen bonds. A positively charged particle (such as a proton), located between the two negatively charged residues, may lead to a stable interaction of the two negative residues. This paper analyzes close Asp-Glu pairs of residues in a large set of protein chains; 840 such pairs of residues were identified, of which 28% were stabilized by a metal ion, 12% by a positive residue nearby and 60% are likely to be stabilized by a proton. The absence of apparent structural constraints, secondary structure preferences, somewhat lower B-factors and a distinct correlation between pH and the minimal O-O distance in carboxylate pairs suggest that most of the abnormally close pairs could indeed be stabilized by a shared proton. Implications for protein stability and modeling are discussed.     

The criteria of identification of "critical" populations in aquatic radiochemoecology

Data on chromosome mutagenesis levels in populations of aquatic organisms in the Black and the Aegean Seas, the Danube and the Dnieper Rivers, the 30-km zone of ChNPP are presented. The highest level of mutagenesis was observed in hydrobionts populations in the 10-km zone of the ChNPP. The obvious damaged effects of ionizing radiation were noted only in these populations. The comparison of the adaptation rate of aquatic crustaceans and worms populations with different reproduction modes was made. It is found that the studied species with sexual reproduction have higher rate of adaptation to the pollution in comparison with species with prevalent asexual reproduction. Hypothetic mechanisms of population adaptation are discussed. On the basis of species and populations characteristics, the criteria for the identification of "critical" populations (species) and an algoritm of ecological risk assessment for them are proposed.     

Modelling hydrogen bonds in NN-dimethylformamide

N, N-dimethylformamide (DMF) is a ‘universal’ solvent with the simplest amide structure. DMF has different interactions with many polymers and biomolecules. It is therefore necessary to study systematically the interactions in DMF itself first. In this study, both FT-IR and two molecular theoretical methods (MP2 and DFT/B3LYP) were used to study various hydrogen bonding interactions in DMF molecules based on its weak H-bonding donors CH/CH₃ and strong H-bonding acceptor C = O. The possible H-bonding donors and acceptors in DMF molecules were first analysed followed by modelling the effect of different structural environments on vC = O bands in infrared spectra. Finally, H-bonding properties including distance, angles and the energy as well as the probability of H-bonding patterns were obtained. The results showed that there exist five possible different weak types of H-bonding dimers; among them, three dimers consist of a pair of weak H-bonds, whereas two other dimers have two pairs of H-bonds, leading to 14 (including eight different) H-bonds. Two types of dimers were dominant, whereas three others can be omitted.     

Modeling hydrogen bonds in three dimensions

The effective potential between two hydrogen bonded atoms is calculated on the basis of the Lippencott-Schroeder bent bond model, taken to be a typical model interaction. We differ from other calculations in that the minimum energy configuration for the proton is treated adiabatically, its position being recomputed at each value of the larger atoms separation. We find the typical hard core to have been a consequence of an artificial restriction of the proton to a fixed angle with the larger atom axis, basically a one-dimensional assumption. Free to move in three dimensions, the proton is squeezed off the axis as the separation narrows, and the hard core feature is gone. Depending on the degree of bond bending, the anharmonicity of the bond may be diminished, eliminated, or even reversed.     

Competition between hydrogen bonds and halogen bonds in complexes of formamidine and hypohalous acids

Quantum chemical calculations have been per-formed for the complexes of formamidine (FA) and hypohalous acid (HOX, X = F, Cl, Br, I) to study their structures, properties, and competition of hydrogen bonds with halogen bonds. Two types of complexes are formed mainly through a hydrogen bond and a halogen bond, respectively, and the cyclic structure is more stable. For the F, Cl, and Br complexes, the hydrogen-bonded one is more stable than the halogen-bonded one, while the halogen-bonded structure is favorable for the I complexes. The associated H-O and X-O bonds are elongated and exhibit a red shift, whereas the distant ones are contracted and display a blue shift. The strength of hydrogen and halogen bonds is affected by F and Li substitutents and it was found that the latter tends to smooth differences in the strength of both types of interactions. The structures, properties, and interaction nature in these complexes have been understood with natural bond orbital (NBO) and atoms in molecules (AIM) theories.     

Evaluation of C-H cdots,three dots,centered O hydrogen bonds in native and misfolded proteins

Non-traditional C-H cdots, three dots, centered Y hydrogen bonds, in which a carbon atom acts as the hydrogen donor and an electronegative atom Y (Y=N, O or S) acts as the acceptor, have been reported in proteins, but their importance in protein structures is not well established. Here, we present the results of three computational tests that examine the significance of C-H cdots, three dots, centered Y bonds involving the C(alpha) in proteins. First, we compared the number of C(alpha)-H cdots, three dots, centered Y bonds in native structures with two sets of compact, energy-minimized decoy structures. The decoy structures contain about as many C(alpha)-H cdots, three dots, centered Y bonds as the native structures, indicating that the constraints of chain connectivity and compactness can lead to incidental formation of C(alpha)-H cdots, three dots, centered Y bonds. Secondly, we examined whether short C(alpha)-H cdots, three dots, centered Y bonds have a tendency to be linear, as is expected for a cohesive hydrogen-bonding interaction. The native structures do show this trend, but so does one of the decoy sets, suggesting that this criterion is also not sufficient to indicate a stabilizing interaction. Finally, we examined the preference for C(alpha)-H cdots, three dots, centered Y bond donors to be near to strong hydrogen bond acceptors. In the native proteins, the alpha protons attract strong acceptors like oxygen atoms more than weak acceptors. In contrast, hydrogen bond donors in the decoy structures do not distinguish between strong and weak acceptors. Thus, any individual C(alpha)-H cdots, three dots, centered Y bond may be fortuitous and occur due to the polypeptide connectivity and compactness. Taken collectively, however, C(alpha)-H cdots, three dots, centered Y bonds provide a weakly cohesive force that stabilizes proteins.     

Negatively charged residues and hydrogen bonds tune the ligand histidine pKa values of Rieske iron-sulfur proteins

Rieske proteins carry a redox-active iron-sulfur cluster, which is bound by two histidine and two cysteine side chains. The reduction potential of Rieske proteins depends on pH. This pH dependence can be described by two pK(a) values, which have been assigned to the two iron-coordinating histidines. Rieske proteins are commonly grouped into two major classes: Rieske proteins from quinol-oxidizing cytochrome bc complexes, in which the ligand histidines titrate in the physiological pH range, and bacterial ferredoxin Rieske proteins, in which the ligand histidines are protonated at physiological pH. In the study presented here, we have calculated pK(a) values of the cluster ligand histidines using a combined density functional theory/continuum electrostatics approach. Experimental pK(a) values for a bc-type and a ferredoxin Rieske protein could be reproduced. We could identify functionally important differences between the two proteins: hydrogen bonds toward the cluster, which are present in bc-type Rieske proteins, and negatively charged residues, which are present in ferredoxin Rieske proteins. We removed these differences by mutating the proteins in our calculations. The Rieske centers in the mutated proteins have very similar pK(a) values. We thus conclude that the studied structural differences are the main reason for the different pH-titration behavior of the proteins. Interestingly, the shift caused by neutralizing the negative charges in ferredoxin Rieske proteins is larger than the shift caused by removing the hydrogen bonds toward the cluster in bc-type Rieske proteins.     

Binding mechanisms of TATA box-binding proteins: DNA kinking is stabilized by specific hydrogen bonds

One of the common mechanisms of DNA bending by minor groove-binding proteins is the insertion of protein side chains between basepair steps, exemplified in TBP (TATA box-binding protein)/DNA complexes. At the central basepair step of the TATA box TBP produces a noticeable decrease in twist and an increase in roll, while engaging in hydrogen bonds with the bases and sugars. This suggests a mechanism for the stabilization of DNA kinks that was explored here with ab initio quantum mechanical calculations and molecular dynamics/potential of mean force calculations. The hydrogen bonds are found to contribute the energy necessary to drive the conformational transition at the central basepair step. The Asn, Thr, and Gly residues involved in hydrogen bonding to the DNA bases and sugar oxygens form a relatively rigid motif in TBP. The interaction of this motif with DNA is found to be responsible for inducing the untwisting and rolling of the central basepair step. Notably, direct readout is shown not to be capable of discriminating between AA and AT steps, as the strength of the hydrogen bonds between TBP and the DNA are the same for both sequences. Rather, the calculated free energy cost for an equivalent conformational transition is found to be sequence-dependent, and is calculated to be higher for AA steps than for AT steps.     

Formation of disulfide bonds in proteins and peptides

For many proteins and peptides, disulfide bridges are prerequisite for their proper biological function. Many commercialized proteins are crosslinked by disulfide bridges that increase their resistance to destructive effects of extreme environment used in industrial processes or protect protein-based therapeutics from rapid proteolytic degradation. Manufacturing of these products must take into account oxidative refolding--a formation of native disulfide bonds by specific pairs of cysteines located throughout a sequence of linear protein. This review describes basic and practical aspects of oxidative folding that should be considered while designing and optimizing manufacturing of proteins using chemical synthesis, semi-synthesis and a recombinant expression.     

Direct identification of NH...N hydrogen bonds in non-canonical base pairs of RNA by NMR spectroscopy.

It is shown that the recently developed quantitative J(NN)HNN-COSY experiment can be used for the direct identification of hydrogen bonds in non-canonical base pairs in RNA. Scalar(2h)J(NN)couplings across NH.N hydrogen bonds are observed in imino hydrogen bonded GA base pairs of the hpGA RNA molecule, which contains a tandem GA mismatch, and in the reverse Hoogsteen AU base pairs of the E-loop of Escherichia coli 5S rRNA. These scalar couplings correlate the imino donor(15)N nucleus of guanine or uridine with the acceptor N1 or N7 nucleus of adenine. The values of the corresponding(2h)J(NN)coupling constants are similar in size to those observed in Watson-Crick base pairs. The reverse Hoogsteen base pairs could be directly detected for the E-loop of E.coli 5S rRNA both in the free form and in a complex with the ribosomal protein L25. This supports the notion that the E-loop is a pre-folded RNA recognition site that is not subject to significant induced conformational changes. Since Watson-Crick GC and AU base pairs are also readily detected the HNN-COSY experiment provides a useful and sensitive tool for the rapid identification of RNA secondary structure elements.