Hydrogen bonding of nucleotide base pairs: Application of the PM3 method

The ability of the pm3 semiempirical quantum mechanical method to reproduce hydrogen bonding in nucleotide base pairs was assessed. Results of pm3 calculations on the nucleotides 2′-deoxyadenosine 5′-monophosphate (pdA), 2′-deoxyguanosine 5′-monophosphate (pdG), 2′-deoxycytidine 5′-monophosphate (pdC), and 2′-deoxythymidine 5′-monophosphate (pdT) and the base pairs pdA–pdT, pdG–pdC, and pdG(syn)–pdC are presented and discussed. The pm3 method is the first of the parameterized nddo quantum mechanical models with any ability to reproduce hydrogen bonding between nucleotide base pairs. Intermolecular hydrogen bond lengths between nucleotides displaying Watson–Crick base pairing are 0.1–0.2 Å less than experimental results. Nucleotide bond distances, bond angles, and torsion angles about the glycosyl bond (χ), the C4′C5′ bond (γ), and the C5′O5′ bond (β) agree with experimental results. There are many possible conformations of nucleotides. pm3 calculations reveal that many of the most stable conformations are stabilized by intramolecular CHO hydrogen bonds. These interactions disrupt the usual sugar puckering. The stacking interactions of a dT–pdA duplex are examined at different levels of gradient optimization. The intramolecular hydrogen bonds found in the nucleotide base pairs disappear in the duplex, as a result of the additional constraints on the phosphate group when part of a DNA backbone. Sugar puckering is reproduced by the pm3 method for the four bases in the dT–pdA duplex. pm3 underestimates the attractive stacking interactions of base pairs in a B-DNA helical conformation. The performance of the pm3 method implemented in SPARTAN is contrasted with that implemented in MOPAC. At present, accurate ab initio calculations are too timeconsuming to be of practical use, and molecular mechanics methods cannot be used to determine quantum mechanical properties such as reaction-path calculations, transition-state structures, and activation energies. The pm3 method should be used with extreme caution for examination of small DNA systems. Future parameterizations of semiempirical methods should incorporate base stacking interactions into the parameterization data set to enhance the ability of these methods.

Ability of the PM3 Quantum Mechanical Method to Model Intermolecular Hydrogen Bonding between Neutral Molecules

The PM3 semiempirical quantum-mechanical method was found to systematically describe intermolecular hydrogen bonding in small polar molecules. PM3 shows charge transfer from the donor to acceptor molecules on the order of 0.02-0.06 units of charge when strong hydrogen bonds are formed. The PM3 method is predictive; calculated hydrogen bond energies with an absolute magnitude greater than 2 kcal mol-' suggest that the global minimum is a hydrogen bonded complex; absolute energies less than 2 kcal mol-' imply that other van der Waals complexes are more stable. The geometries of the PM3 hydrogen bonded complexes agree with high-resolution spectroscopic observations, gas electron diffraction data, and high-level ab initio calculations. The main limitations in the PM3 method are the underestimation of hydrogen bond lengths by 0.1-0.2 for some systems and the underestimation of reliable experimental hydrogen bond energies by approximately 1-2 kcal mol-l. The PM3 method predicts that ammonia is a good hydrogen bond acceptor and a poor hydrogen donor when interacting with neutral molecules. Electronegativity differences between F, N, and 0 predict that donor strength follows the order F > 0 > N and acceptor strength follows the order N > 0 > F. In the calculations presented in this article, the PM3 method mirrors these electronegativity differences, predicting the F-H- - -N bond to be the strongest and the N-H- - -F bond the weakest. It appears that the PM3 Hamiltonian is able to model hydrogen bonding because of the reduction of two-center repulsive forces brought about by the parameterization of the Gaussian core-core interactions. The ability of the PM3 method to model intermolecular hydrogen bonding means reasonably accurate quantum-mechanical calculations can be applied to small biologic systems.

Direct measurement of hydrogen bonding in DNA nucleotide bases by atomic force microscopy.

We have used self-assembled purines and pyrimidines on planar gold surfaces and on gold-coated atomic force microscope (AFM) tips to directly probe intermolecular hydrogen bonds. Electron spectroscopy for chemical analysis (ESCA) and thermal programmed desorption (TPD) measurements of the molecular layers suggested monolayer coverage and a desorption energy of about 25 kcal/mol. Experiments were performed under water, with all four DNA bases immobilized on AFM tips and flat surfaces. Directional hydrogen-bonding interaction between the tip molecules and the surface molecules could be measured only when opposite base-pair coatings were used. The directional interactions were inhibited by excess nucleotide base in solution. Nondirectional van der Waals forces were present in all other cases. Forces as low as two interacting base pairs have been measured. With coated AFM tips, surface chemistry-sensitive recognition atomic force microscopy can be performed.

Toxin-induced pore formation is hindered by intermolecular hydrogen bonding in sphingomyelin bilayers

Sticholysin I and II (StnI and StnII) are pore-forming toxins that use sphingomyelin (SM) for membrane binding. We examined how hydrogen bonding among membrane SMs affected the StnI- and StnII-induced pore formation process, resulting in bilayer permeabilization. We compared toxin-induced permeabilization in bilayers containing either SM or dihydro-SM (lacking the trans 4 double bond of the long-chain base), since their hydrogen-bonding properties are known to differ greatly. We observed that whereas both StnI and StnII formed pores in unilamellar vesicles containing palmitoyl-SM or oleoyl-SM, the toxins failed to similarly form pores in vesicles prepared from dihydro-PSM or dihydro-OSM. In supported bilayers containing OSM, StnII bound efficiently, as determined by surface plasmon resonance. However, StnII binding to supported bilayers prepared from dihydro-OSM was very low under similar experimental conditions. The association of the positively charged StnII (at pH 7.0) with unilamellar vesicles prepared from OSM led to a concentration-dependent increase in vesicle charge, as determined from zeta-potential measurements. With dihydro-OSM vesicles, a similar response was not observed. Benzyl alcohol, which is a small hydrogen-bonding compound with affinity to lipid bilayer interfaces, strongly facilitated StnII-induced pore formation in dihydro-OSM bilayers, suggesting that hydrogen bonding in the interfacial region originally prevented StnII from membrane binding and pore formation. We conclude that interfacial hydrogen bonding was able to affect the membrane association of StnI- and StnII, and hence their pore forming capacity. Our results suggest that other types of protein interactions in bilayers may also be affected by hydrogen-bonding origination from SMs.

Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

Interactions of phenyldithioesters with gold nanoparticles (AuNPs): Implications for AuNP functionalization and molecular barcoding of AuNP assemblies

The interactions of phenyldithioesters with gold nanoparticles (AuNPs) have been studied by monitoring changes in the surface plasmon resonance (SPR), depolarised light scattering, and surface enhanced Raman spectroscopy (SERS). Changes in the SPR indicated that an AuNP-phenyldithioester charge transfer complex forms in equilibrium with free AuNPs and phenyldithioester. Analysis of the Langmuir binding isotherms indicated that the equilibrium adsorption constant, Kads, was 2.3 ± 0.1 × 106 M−1, which corresponded to a free energy of adsorption of 36 ± 1 kJ mol−1. These values are comparable to those reported for interactions of aryl thiols with gold and are of a similar order of magnitude to moderate hydrogen bonding interactions. This has significant implications in the application of phenyldithioesters for the functionalization of AuNPs. The SERS results indicated that the phenyldithioesters interact with AuNPs through the C═S bond, and the molecules do not disassociate upon adsorption to the AuNPs. The SERS spectra are dominated by the portions of the molecule that dominate the charge transfer complex with the AuNPs. The significance of this in relation to the use of phenyldithioesters for molecular barcoding of nanoparticle assemblies is discussed.

Hydrogen bonding in the water stabilized structure of the modified unsymmetrically substituted viologen chromophore, N-ethyl-N'-(2-phosphonomethyl)-4,4'-bipyridylium dichloride

The crystal structure of the modified unsymmetrically N, N'-substituted viologen chromophore, N-ethyl- N'-(2-phosphonoethyl)-4, 4'-bipyridinium dichloride 0.75 hydrate. (1) has been determined. Crystals are triclinic, space group P-1 with Z = 2 in a cell with a = 7.2550(1), b = 13.2038(5), c = 18.5752(7) Å, α = 86.495(3), β = 83.527(2), γ = 88.921(2)o. The two independent but pseudo-symmetrically related cations in the asymmetric unit form one-dimensional hydrogen-bonded chains through short homomeric phosphonic acid O-H...O links [2.455(4), 2.464(4)A] while two of the chloride anions are similarly strongly linked to phosphonic acid groups [O-H…Cl, 2.889(4), 2.896(4)Å]. The other two chloride anions together with the two water molecules of solvation (one with partial occupancy) form unusual cyclic hydrogen-bonded bis(Cl...water) dianion units which lie between the layers of bipyridylium rings of the cation chain structures with which they are weakly associated.

Structure of selected basic copper (II) sulphate minerals based upon spectroscopy : implications for hydrogen bonding

The application of near-infrared and infrared spectroscopy has been used for identification and distinction of basic Cu-sulphates that include devilline, chalcoalumite and caledonite. Near-infrared spectra of copper sulphate minerals confirm copper in divalent state. Jahn-Teller effect is more significant in chalcoalumite where 2B1g ® 2B2g transition band shows a larger splitting (490 cm-1) confirming more distorted octahedral coordination of Cu2+ ion. One symmetrical band at 5145 cm-1 with shoulder band 5715 cm-1 result from the absorbed molecular water in the copper complexes are the combinations of OH vibrations of H2O. One sharp band at around 3400 cm-1 in IR common to the three complexes is evidenced by Cu-OH vibrations. The strong absorptions observed at 1685 and 1620 cm-1 for water bending modes in two species confirm strong hydrogen bonding in devilline and chalcoalumite. The multiple bands in v3 and v4(SO4)2- stretching regions are attributed to the reduction of symmetry to the sulphate ion from Td to C2V. Chalcoalumite, the excellent IR absorber over the range 3800-500 cm-1 is treated as most efficient heat insulator among the Cu-sulphate complexes.

Spectroscopy of selected copper group minerals : chalcophyllite and chenevixite-implications for hydrogen bonding

NIR and IR spectroscopy has been applied for detection of chemical species and the nature of hydrogen bonding in arsenate complexes. The structure and spectral properties of copper(II) arsenate minerals chalcophyllite and chenevixite are compared with copper(II) sulphate minerals devilline, chalcoalumite and caledonite. Split NIR bands in the electronic spectrum of two ranges 11700-8500 cm-1 and 8500-7200 cm-1 confirm distortion of octahedral symmetry for Cu(II) in the arsenate complexes. The observed bands with maxima at 9860 and 7750 cm-1 are assigned to Cu(II) transitions 2B1g ® 2B2g and 2B1g ® 2A1g. Overlapping bands in the NIR region 4500-4000 cm-1 is the effect of multi anions OH-, (AsO4)3- and (SO4)2-. The observation of broad and diffuse bands in the range 3700-2900 cm-1 confirms strong hydrogen bonding in chalcophyllite relative to chenevixite. The position of the water bending vibrations indicates the water is strongly hydrogen bonded in the mineral structure. The strong absorption feature centred at 1644 cm-1 in chalcophyllite indicates water is strongly hydrogen bonded in the mineral structure. The H2O-bending vibrations shift to low wavenumbers in chenevixite and an additional band observed at 1390 cm-1 is related to carbonate impurity. The characterisation of IR spectra by ν3 antisymmetric stretching vibrations of (SO4)2- and (AsO4)3 ions near 1100 and 800 cm-1 respectively is the result of isomorphic substitution for arsenate by sulphate in both the minerals of chalcophyllite and chenevixite.

Three-dimensional hydrogen-bonded structures in the guanidinium salts of the monocyclic dicarboxylic acids rac-trans-cyclohexane-1,2-dicarboxylic acid (2:1, anhydrous), isophthalic acid (1:1, monohydrate) and terephthalic acid (2:1, trihydrate).

The structures of bis(guanidinium)rac-trans-cyclohexane-1,2-dicarboxylate, 2(CH6N3+) C8H10O4- (I), guanidinium 3-carboxybenzoate monohydrate CH6N3+ C8H5O4- . H2O (II) and bis(guanidinium) benzene-1,4-dicarboxylate trihydrate, 2(CH6N3+) C8H4O4^2- . 3H2O (III) have been determined and the hydrogen bonding in each examined. All three compounds form three-dimensional hydrogen-bonded framework structures. In anhydrous (I), both guanidinium cations give classic cyclic R2/2(8) N--H...O,O'(carboxyl) and asymmetric cyclic R1/2(6) hydrogen-bonding interactions while one cation gives an unusual enlarged cyclic interaction with O acceptors of separate ortho-related carboxyl groups [graph set R2/2(11)]. Cations and anions also associate across inversion centres giving cyclic R2/4(8) motifs. In the 1:1 guanidinium salt (II), the cation gives two separate cyclic R1/2(6) interactions, one with a carboxyl O-acceptor, the other with the water molecule of solvation. The structure is unusual in that both carboxyl groups give short inter-anion O...H...O contacts, one across a crystallographic inversion centre [2.483(2)\%A], the other about a two-fold axis of rotation [2.462(2)\%A] with a half-occupancy hydrogen delocalized on the symmetry element in each. The water molecule links the cation--anion ribbon structures into a three-dimensional framework. In (III), the repeating molecular unit comprises a benzene-1,4-dicarboxylate dianion which lies across a crystallographic inversion centre, two guanidinium cations and two water molecules of solvation (each set related by two-fold rotational symmetry), and a single water molecule which lies on a two-fold axis. Each guanidinium cation gives three types of cyclic interactions with the dianions: one R^1^~2~(6), the others R2/3(8) and R3/3(10) (both of these involving the water molecules), giving a three-dimensional structure through bridges down the b cell direction. The water molecule at the general site also forms an unusual cyclic R2/2(4) homodimeric association across an inversion centre [O--H...O, 2.875(2)\%A]. The work described here provides further examples of the common cyclic guanidinium cation...carboxylate anion hydrogen-bonding associations as well as featuring other less common cyclic motifs.

Hydrogen bonding in proton-transfer compounds of 5-sulfosalicylic acid with a series of aliphatic nitrogen Lewis bases

The crystal structures of the proton-transfer compounds of 5-sulfosalicylic acid (3-carboxy-4-hydroxybenzenesulfonic acid) with the aliphatic nitrogen Lewis bases, hydroxylamine, triethylamine, pyrrolidine, morpholine, N-methylmorpholine and piperazine, viz. hydroxyammonium 3-carboxy-4-hydroxybenzenesulfonate (1), triethylaminium 3-carboxy-4-hydroxybenzenesulfonate (2), pyrrolidinium 3-carboxy-4-hydroxybenzenesulfonate monohydrate (3), morpholinium 3-carboxy-4-hydroxybenzenesulfonate monohydrate (4), N-methylmorpholinium 3-carboxy-4-hydroxybenzenesulfonate monohydrate (5) and piperazine-1,4-diium bis(3-carboxy-4-hydroxybenzenesulfonate) hexahydrate (6) have been determined and their comparative structural features and hydrogen-bonding patterns described. Crystals of 4 are triclinic, space group P-1 while the remainder are monoclinic with space group either P21/c (1 - 3) or P21/n (5, 6). Unit cell dimensions and contents are: for 1, a = 5.0156(3), b = 10.5738(6), c = 18.4785(9) Å, β = 96.412(5)o, Z = 4; for 2, a = 8.4998(4), b = 12.3832(6), c = 15.4875(9) Å, β = 102.411(5)o, Z = 4; for 3, a = 6.8755(2), b = 15.5217(4), c = 12.8335(3) Å, β = 92.074(2)o, Z = 4; for 4, a = 6.8397(2), b = 12.9756(5), c = 15.8216(6) Å, α = 90.833(3), β = 95.949(3), γ = 92.505(3)o, Z = 4; for 5, a = 7.0529(3), b = 13.8487(7), c = 15.6448(6) Å, β = 90.190(6)o, Z = 4; for 6, a = 7.0561(2), b = 15.9311(4), c = 12.2102(3) Å, β = 100.858(3)o, Z = 2. The hydrogen bonding generates structures which are either two-dimensional (2 and 5) or three-dimensional (1, 3, 4 and 6). Compound 6 represents the third reported structure of a salt of 5-sulfosalicylic acid having a dicationic piperazine species.

Structure of selected basic zinc/copper (II) phosphate minerals based upon near-infrared spectroscopy : implications for hydrogen bonding

The NIR spectra of reichenbachite, scholzite and parascholzite have been studied at 298 K. The spectra of the minerals are different, in line with composition and crystal structural variations. Cation substitution effects are significant in their electronic spectra and three distinctly different electronic transition bands are observed in the near-infrared spectra at high wavenumbers in the 12000-7600 cm-1 spectral region. Reichenbachite electronic spectrum is characterised by Cu(II) transition bands at 9755 and 7520 cm-1. A broad spectral feature observed for ferrous ion in the 12000-9000 cm-1 region both in scholzite and parascholzite. Some what similarities in the vibrational spectra of the three phosphate minerals are observed particularly in the OH stretching region. The observation of strong band at 5090 cm-1 indicates strong hydrogen bonding in the structure of the dimorphs, scholzite and parascholzite. The three phosphates exhibit overlapping bands in the 4800-4000 cm-1 region resulting from the combinations of vibrational modes of (PO4)3- units.

The hydrogen-bonded coordination polymer tetracesium bis(5-nitroisophthalate) heptahydrate

The coordination polymer complex tetracesium bis(5-nitroisophthalate) heptahydrate [Cs4(C8H3NO6)2 (H2O)7]n has been synthesized and characterized using single-crystal X-ray diffraction. Crystals are monoclinic, space group P21/c, with Z = 4 in a cell with dimensions a = 12.3213(3), b =6.7557(2) c = 36.2020(9) Å, β = 90.548(2)o. The complex is based on a repeating unit comprising four independent and different Cs coordination centres, two 6-coordinate, and two 8-coordinate [Cs-O, range 2.959(5)-3.386(5)Å], and seven water molecules, two of which are monodentate and the other five bridging, while all other oxygen atoms in the structure, including those of the nitro groups form inter-Cs bridges. Extensive water O-H…O hydrogen-bonding interactions give a three-dimensional framework. This structure represents the first of an alkali metal compound of 5-nitroisophthalic acid that has been reported.