Effect of hydrogen bonds on the energetics of electron transfer

For a model system consisting of a special pair of bacteriochlorophyll molecules (P) and a primary quinone with the nearest environment (QA) (which are acceptor and donor in the recombination reaction in Rhodobacter sphaeroides reaction center, respectively), energies of P+QA(-) and PQA states were calculated. Calculations were performed using several stable QA conformations differing by the positions of hydrogen bond protons. Essential influence of proton positions on the energy of vertical transition P+QA(-) --> PQA was shown.     

General features of the energetics of complex formation between ligand and nucleic acids

The analysis of the energy contributions of various physical factors to the complex formation between biologically active compounds and nucleic acids in aqueous solution was performed. A comparison of the energy parameters was made for ligand-ligand, intercalator-DNA, MGB-DNA and ligand-RNA groups. It was shown that the energetics of these reactions is of compensatory nature. Physical factors exerting the most pronounced influence on the energy parameters were identified. Correlation of the energy contributions to MGB-DNA complex formation and its biological effect was found.     

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.     

Effect of hydrogen bonds on the energetics of macromolecules in the course of electron transfer

For a model system consisting of a bacteriochlorophyll dimer (P) and a primary quinone with the nearest environment (Q_A), which are the electron donor and acceptor in the recombination reaction in the Rhodobacter spheroides reaction center, the energies of states P⁺Q _A ^? and PQ_A have been calculated at several stable conformations of Q_A that differ in the positions of the proton involved in the hydrogen bond. It is shown that the position of the proton has a considerable influence on the energy of vertical transition P⁺Q _A ^? → PQ_A.     

Conformational studies of nucleic acids. V. Sequence specificities in the conformational energetics of oligonucleotides: the homo-tetramers

The conformational energetics of the tetranucleoside triphosphates d(ApApApA), d(GpGpGpG), d(CpCpCpC), d(TpTpTpT), ApApApA, GpGpGpG, CpCpCpC, and UpUpUpU are thoroughly examined using a classical potential energy function. The sugar modeling method and multiple correlation functions derived in previous papers of this series are utilized in these examinations. The data are analyzed and compared in terms of the energy profiles for rotation about the conformation-determining torsion angles in the tetramers. Overall, the predictions are in reasonable qualitative agreement with the existing experimental data. It is found that the base type does not greatly affect the locations of the important minima in these profiles, but rather exerts a large influence on the relative depths of the minima and the barriers to conversion between them. Conformational sequence dependence is manifest to a greater extent by the DNA tetramers than the RNA tetramers. Of the DNA tetramers, d(CpCpCpC) appears, from the results presented herein, to have the greatest potential for polymorphism. This and other findings are analyzed in terms of the preferences of particular DNA sequences for either the A-, B-, or Z-conformation.     

Ligation of oligonucleotides to nucleic acids or proteins via disulfide bonds.

We have developed general methods for joining together, via cleavable disulfide bonds, either two unprotected polynucleotides or a polynucleotide and a peptide or protein. To join two oligonucleotides, each is first converted to an adduct in which cystamine is joined to the 5'-terminal phosphate of the oligonucleotide by a phosphoramidate bond. The adducts are mixed and reduced with dithiothreitol. The dithiothreitol is then removed by dialysis. Oxidation by atmospheric oxygen occurs to yield the required dimer. To join an oligonucleotide to a cysteine-containing peptide or protein, the 5'-cystamine oligomer is first converted to a 2'-pyridyldisulfide adduct and then reacted with an excess of the peptide or protein. If the peptide does not contain a free cysteine residue, it is first treated with iminothiolane to introduce one or more sulfhydryl groups. We have used these procedures to join a 16 mer deoxynucleotide probe and MDV-1 RNA, a substrate of Q beta RNA polymerase. This adduct hybridizes with a complementary target DNA. We have also joined a 16mer probe to peroxidase and MDV-1 RNA to human IgG. The probe-peroxidase adduct maintains enzymatic activity and the MDV-1 RNA-IgG adduct binds to a complementary anti-IgG.     

Model studies of interactions between nucleic acids and proteins: hydrogen bonding of amides with nucleic acid bases.

The formation of hydrogen bonded complexes between nucleic acid bases and acetamide has been studied by nuclear magnetic resonance in CDC13 at different temperatures. Pairs of hydrogen bonds are formed when acetamide binds to nucleic acid bases. Thermodynamic parameters have been computed and compared to those obtained for the association of carboxylic acids with nucleic acid bases. The role of hydrogen bonded complexes in the association of proteins with nucleic acids is discussed.     

Conformational studies of nucleic acids: IV. The conformational energetics of oligonucleotides: d(ApApApA) and ApApApA

Utilizing a new method for modeling furanose pseudorotation (D. A. Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A. Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA. The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: delta (C5'-C4'-C3'-O3')-chi (O4'-C1'-N9-C4), where the energetic basis for an approximately linear delta-chi correlation can be seen: zeta (C3'-O3'-P-O5')-delta, in which the experimentally observed linear correlation between zeta and delta in DNA(220 degrees less than zeta less than 280 degrees) is clearly predicted; zeta-epsilon (C4'-C3'-O3'-P), which shows that epsilon increases with decreasing zeta less than 260 degrees; alpha (O3'-P-O5'-C5')-gamma (O5'-C5'-C4'-C3') where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion delta is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in alpha and gamma. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.     

Chemical etiology of nucleic acids: aminopropyl nucleic acids (APNAs)

Aminopropyl nucleic acids (APNAs) are constitutionally simple nucleic acid alternatives with one stereogenic center per nucleotide, and with the potential to hybridize with RNA and to exert catalytic functions. We have developed a protecting group strategy to synthesize APNAs, although in a not very efficient way. Isolation and purification of APNAs proved to be difficult. Their structures might be more suited to function as potential catalytic polymers than as information systems that may evolve into RNA.     

Hardware and software in biology: Simultaneous origin of proteins and nucleic acids via hydrogen cyanide polymers

Hydrogen cyanide polymers—heterogeneous solids varying in color from yellow to orange to red to black—may be among the organic macromolecules most readily formed within the solar system. Current studies of these ubiquitous compounds point to the presence of polyamidine structures readily converted by water to polypeptides. Implications for prebiotic chemistry are profound. Primitive Earth may have been covered by HCN polymers as well as other organic compounds, either through bolide bombardment or by photochemical reactions in a reducing atmosphere. Membrane material—carboxylic acids, carbohydrates, polypeptides—accumulated in lakes and oceans, while in the absence of water, on land, polyamidines could have been the original dehydrating agents directing the synthesis of nucleosides and nucleotides from available sugars, phosphates and nitrogen bases. Most significant would have been the parallel synthesis of polypeptides and polynucleotides ariaing from the dehydrating action of polyamidines on nucleotides. Metabolic material—hardware—thus arose separately from genetic components—software—with subsequent interfacting producing the first replicating protocells.On our dynamic planet this polypeptide-polynucleotide symbiosis mediated by polyamidines may have set the pattern for the evolution of protein-nucleic acid systems controlled by enzymes, the mode characteristic of life today.     

Stability and fluctuations of amide hydrogen bonds in a bacterial cytochrome c: a molecular dynamics study

Molecular dynamics (MD) simulations on a bacterial cytochrome c were performed to investigate the lifetime and fluctuations of backbone hydrogen bonds and to correlate these data with protection factors for hydrogen exchange measured by NMR spectroscopy (Bartalesi et al. in Biochemistry, 42:10923–10930, 2003). The MD simulations provide a consistent pattern in that long lifetimes of hydrogen bonds go along with small amplitude fluctuations. In agreement with experiments, differences in stability were found with a rather flexible N-terminal segment as compared with a more rigid C-terminal part. Protection factors of backbone hydrogen exchange correlate strongly with the number of contacts but also with hydrogen-bond occupancy, hydrogen-bond survival times, as well as the inverse of fluctuations of backbone atoms and hydrogen-bond lengths derived from MD simulation data. We observed a conformational transition in the C-terminal loop, and significant motion in the N-terminal loop, which can be interpreted as being the structural units involved in the onset of the protein unfolding process in agreement with experimental evidence on mitochondrial cytochrome c. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users. Gernot Kieseritzky and Giulia Morra both contributed equally to this work.     

Validation of nuclear magnetic resonance structures of proteins and nucleic acids: hydrogen geometry and nomenclature

A statistical analysis is reported of 1,200 of the 1,404 nuclear magnetic resonance (NMR)-derived protein and nucleic acid structures deposited in the Protein Data Bank (PDB) before 1999. Excluded from this analysis were the entries not yet fully validated by the PDB and the more than 100 entries that contained < 95% of the expected hydrogens. The aim was to assess the geometry of the hydrogens in the remaining structures and to provide a check on their nomenclature. Deviations in bond lengths, bond angles, improper dihedral angles, and planarity with respect to estimated values were checked. More than 100 entries showed anomalous protonation states for some of their amino acids. Approximately 250,000 (1.7%) atom names differed from the consensus PDB nomenclature. Most of the inconsistencies are due to swapped prochiral labeling. Large deviations from the expected geometry exist for a considerable number of entries, many of which are average structures. The most common causes for these deviations seem to be poor minimization of average structures and an improper balance between force-field constraints for experimental and holonomic data. Some specific geometric outliers are related to the refinement programs used. A number of recommendations for biomolecular databases, modeling programs, and authors submitting biomolecular structures are given.     

Uniformity,ideality, and hydrogen bonds in transmembrane alpha-helices

Protein environments substantially influence the balance of molecular interactions that generate structural stability. Transmembrane helices exist in the relatively uniform low dielectric interstices of the lipid bilayer, largely devoid of water and with a very hydrophobic distribution of amino acid residues. Here, through an analysis of bacteriorhodopsin crystal structures and the transmembrane helix structure from M2 protein of influenza A, some helices are shown to be exceptionally uniform in hydrogen bond geometry, peptide plane tilt angle, and backbone torsion angles. Evidence from both the x-ray crystal structures and solid-state NMR structure suggests that the intramolecular backbone hydrogen bonds are shorter than their counterparts in water-soluble proteins. Moreover, the geometry is consistent with a dominance of electrostatic versus covalent contributions to these bonds. A comparison of structure as a function of resolution shows that as the structures become better characterized the helices become much more uniform, suggesting that there is a possibility that many more uniform helices will be observed, even among the moderate resolution membrane protein structures that are currently in the Protein Data Bank that do not show such features.