Effects of hydrogen-bonding and stacking interactions with amino acids on the acidity of uracil

The effects of hydrogen-bonding interactions with amino acids on the (N1) acidity of uracil are evaluated using (B3LYP) density functional theory. Many different binding arrangements of each amino acid to three uracil binding sites are considered. The effects on the uracil acidity are found to significantly depend upon the nature of the amino acid and the binding orientation, but weakly depend on the binding site. Our results reveal that in some instances small models for the amino acids can be used, while for other amino acids larger models are required to properly describe the binding to uracil. The gas-phase acidity of uracil is found to increase by up to approximately 60 kJ mol(-1) due to discrete hydrogen-bonding interactions. Although (MP2) stacking interactions with aromatic amino acids decrease the acidity of uracil, unexpected increases in the acidity are found when any of the aromatic amino acids, or the backbone, hydrogen bond to uracil. Consideration of enzymatic and aqueous environments leads to decreases in the effects of the amino acids on the acidity of uracil. However, we find that the magnitude of the decrease varies with the nature of the molecule bound, as well as the (gas-phase) binding orientations and strengths, and therefore solvation effects should be considered on a case-by-case basis in future work. Nevertheless, the effects of amino acid interactions within enzymatic environments are as much as approximately 35 kJ mol(-1). The present study has general implications for understanding the nature of active site amino acids in enzymes, such as DNA repair enzymes, that catalyze reactions involving anionic nucleobase intermediates.     

Configurational specificity of stacking interactions in DNA base pairs: A computational analysis

A theoretical study of stacking patterns of various hydrogen-bonded base pair complexes has been undertaken. Modified Rayleigh-Schrodinger perturbation theory for intermediate range interactions, has been employed to evaluate the stacking interactions using multicentered-multipole expansion method. Net atomic charge and corresponding dipole components located at each of the atomic centers have been computed by CNDO/2 method. An analysis of the intermolecular forces involved in the stable formation of the various base pair complexes, has been presented and the results have been discussed in the light of experimental as well as other theoretical observations. The possibility of relative preference of the left-handed configuration for alternating sequences has been quantitatively explored.     

A computational study on the stacking interaction in quinhydrone

The stability and electron density topology of quinhydrone complex was studied using multiple computational levels, including MPW1B95 Truhlar's density functional. The QTAIM analysis demonstrates that an electron population transfer from hydroquinone to quinone monomer accompanies the complex formation. The variations undergone by atomic populations indicate that the electron transfer through HOMO LUMO overlap is combined with a reorganization of the electron density within each monomer. Variations of two- and six-center delocalization indices show a small reduction of electron delocalization in the hydroquinone ring upon complex formation.     

A combined computational and experimental study of the hydrogen-bonded dimers of xanthine and hypoxanthine

In addition to uracil, the noncanonical nucleobases xanthine and hypoxanthine are important lesions that are formed from the canonical bases when a cell is under oxidative stress. It is known that they lead to point mutations; however, more detailed information about their ability to form hydrogen-bonded complexes is not available. In the present paper such information is obtained by a combined experimental and theoretical approach. Accurate association constants of xanthosine and inosine dimers are determined by concentration dependent 1H NMR experiments, and a structural characterization of individual complexes formed in solution is performed through measurements under slow exchange conditions at very low temperatures. An interpretation of the experimental data concerning complex geometries becomes possible through a comparison of measured and computed NMR chemical shifts. Further qualitative insights into the hydrogen bonding abilities of xanthine and hypoxanthine are obtained by a theoretical characterization of all possible pairing modes of xanthine and hypoxanthine dimers and by a comparison with simplified model systems. The influence of a polar medium on the bonding properties is also estimated and the importance of the various effects is discussed. Our analysis shows to what extent secondary electronic and electrostatic effects influence the hydrogen bonding properties of xanthine and hypoxanthine in the gas phase and in polar solvents.     

The effect of methylation on the hydrogen-bonding and stacking interaction of nucleic acid bases

Methylated nucleosides play an important role in DNA/RNA function, and may affect republication by interrupting the base-pairing and base-stacking. In order to investigate the effect of methylation on the interaction between nucleic acid bases, this work presents the hydrogen-bonding and stacking interactions between 5-methylcytosine and guanine (G), cytosine (C) and G, 1-methyladenine and thymine (T), as well as adenine and T. Geometry optimization and potential energy surface scan have been performed for the involved complexes by MP2 calculations. The interaction energies, which were corrected for the basis-set superposition error by the full Boys–Bernardi counterpoise correction scheme, were used to evaluate the interaction intensity of these nucleic acid bases. The atoms in molecules theory and natural bond orbital analysis have been performed to study the hydrogen bonds in these complexes. The result shows that the methyl substitute contributes the stability to these complexes because it enhances either the hydrogen bonding or the staking interaction between nucleic acid bases studied.     

Interplay between edge-to-face aromatic and hydrogen-bonding interactions

The interplay between two important noncovalent interactions involving aromatic rings is studied by means of MP2/6-31++G** ab initio calculations. They indicate that synergistic effects are present in complexes where edge-to-face aromatic interactions and hydrogen-bonding interactions coexist. These synergistic effects have been studied bu using the atoms in molecules theory and the molecular interaction potential with polarization partition scheme. Experimental evidence for such interactions has been obtained from the Cambridge Structural Database.     

The combination of transition metal ions and hydrogen-bonding interactions

This feature article presents an overview of the types of hydrogen bonding interactions involving metal complexes and their functional effects. It shows with recent examples why hydrogen bonds have become a crucial functional and structural element in modern inorganic chemistry. The relevance of this combination in tackling current chemistry challenges such as energy production and the development of new materials and more effective catalysts, sensors and medicines is illustrated.     

Electrostatic control of aromatic stacking interactions

A supramolecular approach has been used to investigate the free energies of intermolecular aromatic stacking interactions. Chemical double mutant cycles have been used to measure the effect of a range of substituents on face-to-face stacking interactions with phenyl and pentafluorophenyl rings. Electrostatic effects dominate the trends in interaction energy.     

Thermodynamics of stacking interactions in proteins

Using a database of 6166 experimental structures taken from the Protein Data Bank, we have studied pair interactions between planar residues (Phe, Tyr, His, Arg, Glu and Asp) in proteins, known as pi-pi interactions. On the basis of appropriate coordinates defining the mutual arrangement of two residues, we have calculated 2-D potentials of mean force aimed at determining the stability of the most probable structures for aromatic-aromatic, aromatic-cation and aromatic-anion bound pairs. Our analysis reveals the thermodynamic relevance and the ubiquity of stacked complexes in proteins.     

Efficient modulation of hydrogen-bonding interactions by remote substituents

[structure: see text] A series of tetralactam macrocycles having different substituents were prepared, and their binding affinities for an adipamide guest were investigated in CDCl3 by 1H NMR titrations. The association constants strongly depend on the substituents, varying up to DeltaDeltaG = 3.4 kcal/mol; electron-donating substituents (OMe, NMe2) decrease the binding affinity, while electron-withdrawing groups (Cl, NO2) increase it. These large substituent effects have been rationalized by secondary repulsions and partial perturbations of intramolecular hydrogen bonds.     

A computational proposal for the experimentally observed discriminatory behavior of hypoxanthine, a weak universal nucleobase

A computational model composed of six nucleobases was used to investigate why hypoxanthine does not yield duplexes of equal stability when paired opposite each of the natural DNA nucleobases. The magnitudes of all nearest-neighbor interactions in a DNA helix were calculated, including hydrogen-bonding, intra- and interstrand stacking interactions, as well as 1-3 intrastrand stacking interactions. Although the stacking interactions in DNA relevant arrangements are significant and account for at least one third of the total stabilization energy in our nucleobase complexes, the trends in the magnitude of the stacking interactions cannot explain the relative experimental melting temperatures previously reported in the literature. Furthermore, although the total hydrogen-bonding interactions explain why hypoxanthine preferentially pairs with cytosine, the experimental trend for the remaining nucleobases (A, T, G) is not explained. In fact, the calculated pairing preference of hypoxanthine matches that determined experimentally only when the sum of all types of nearest-neighbor interactions is considered. This finding highlights a strong correlation between the relative magnitude of the total nucleobase-nucleobase interactions and measured melting temperatures for DNA strands containing hypoxanthine despite the potential role of other factors (including hydration, temperature, sugar-phosphate backbone). By considering a large range of sequence combinations, we reveal that the binding preference of hypoxanthine is strongly dependent on the nucleobase sequence, which may explain the varied ability of hypoxanthine to universally bind to the natural bases. As a result, we propose that future work should closely examine the interplay between the dominant nucleobase-nucleobase interactions and the overall strand stability to fully understand how sequence context affects the universal binding properties of modified bases and to aid the design of new molecules with ambiguous pairing properties.     

Competition of hydrogen-bonding and electrostatic interactions within hybrid polymer multilayers

Using a layer-by-layer sequential adsorption technique, we report the construction of hybrid films in which layers of hydrogen-bonded polymers are embedded within electrostatically associated polyelectrolytes. The components of the hybrid film include a neutral hydrogen-bonding polymer, a weak polycarboxylic acid, and a strong polycation. Depending on the pH value used for the deposition of the electrostatic film, we found two distinctive regimes of film growth. At pHs lower than a critical value, deposition of electrostatic layers occurred on top of hydrogen-bonded stacks to produce hybrid, three-component films. At pHs higher than a critical value, neutral, hydrogen-bonded chains were displaced by the adsorbing chains of the polycation, producing two-component films. The property of the hydrogen-bonded stacks of hybrid films to be selectively dissolved by exposing them to a high pH makes these films promising candidates for producing free polyelectrolyte films.     

A heterodifunctionalised ferrocene derivative that self-assembles in solution through complementary hydrogen-bonding interactions

A heterodifunctionalised ferrocene, containing a carboxylic acid and an amidopyridine unit, self-assembles in organic solvents through complementary hydrogen bonds.     

Femtosecond pump-probe measurements of solvation by hydrogen-bonding interactions.

An additional ultrafast blue shift in the transient absorption spectra of hydrogen-bonding complexes of a strong photoacid, 8-hydroxypyrene 1,3,6-trisdimethylsulfonamide (HPTA), over the solvation response of the uncomplexed HPTA and also over that of the methoxy derivative of the photoacid (MPTA) in the presence of the hydrogen-bonding base was observed on optical excitation of the photoacid. The additional 55 +/- 10 fs solvation response was found to be about 35 % and 19% of the total C(t) of HPTA in dichloromethane (DCM) when it was hydrogen-bonded to dimethylsulfoxide (DMSO) and dioxane, respectively, and about 29% of the total C(t) of HPTA in dichloroethane (DCE) when it was hydrogen-bonded to DMSO. We have assigned this additional dynamic spectral shift to a transient change in the hydrogen bond (O-H...O) that links HPTA to the complexing base, after the electronic excitation of the photoacid.     

Local nature of substituent effects in stacking interactions

Popular explanations of substituent effects in π-stacking interactions hinge upon substituent-induced changes in the aryl π-system. This entrenched view has been used to explain substituent effects in countless stacking interactions over the past 2 decades. However, for a broad range of stacked dimers, it is shown that substituent effects are better described as arising from local, direct interactions of the substituent with the proximal vertex of the other ring. Consequently, substituent effects in stacking interactions are additive, regardless of whether the substituents are on the same or opposite rings. Substituent effects are also insensitive to the introduction of heteroatoms on distant parts of either stacked ring. This local, direct interaction viewpoint provides clear, unambiguous explanations of substituent effects for myriad stacking interactions that are in accord with robust computational data, including DFT-D and new benchmark CCSD(T) results. Many of these computational results cannot be readily explained using traditional π-polarization-based models. Analyses of stacking interactions based solely on the sign of the electrostatic potential above the face of an aromatic ring or the molecular quadrupole moment face a similar fate. The local, direct interaction model provides a simple means of analyzing substituent effects in complex aromatic systems and also offers simple explanations of the crystal packing of fluorinated benzenes and the recently published dependence of the stability of protein-RNA complexes on the regiochemistry of fluorinated base analogues [J. Am. Chem. Soc.2011, 133, 3687-3689].     

Infrared and Raman spectroscopic characterization of the hydrogen-bonding network in l-serine crystal

The IR spectra (4000–400 cm⁻¹) of neat and isotopically substituted (ND/OD ≤ 10% D and ≅30% D) polycrystalline l-serine (α-amino-β-hydroxypropionic acid; HO–CH₂–CH(NH₃)⁺–COO⁻) were recorded in the temperature range 300–10 K and assigned. The isotopic-doping/low-temperature methodology, which allows for decoupling of individual proton vibrational modes from the crystal bulk vibrations, was used for estimating the lengths and energies of the different H-bonds present in l-serine crystal. To this end, the frequency shifts observed in both the NH/OH stretching and out-of-plane bending spectral regions (relatively to reference values for these vibrations in non-hydrogen-bonded l-serine molecules) were used, together with previously developed empirical correlations between these spectral parameters and the H-bond properties. In addition, the room-temperature Raman spectrum (4000–150 cm⁻¹) of a single crystal of neat l-serine was also recorded and interpreted. A systematic comparison was made between the spectroscopic data obtained currently for l-serine and previously for dl-serine, revealing that the vibrational spectra of the two crystals reflect well the different characteristics of their hydrogen-bond networks, and also correlate accurately with the different susceptibility of the two crystals to pressure-induced strain.     

A computational study of interactions between acetic acid and water molecules

Density functional theory calculations were performed for the title reactions to elucidate the difference between the strong cyclic hydrogen bond of (Me-COOH)(2) and the electrolytic dissociation, MeCOOH <==> Me-COO(-) + H(+), as a weak acid. The association of water clusters with acetic acid dimers strengthens the cyclic hydrogen bond. A nucleophilic attack of the carboxylic carbon by a water cluster leads to a first zwitterionic intermediate, MeCOO(-) + H(3)O(+) + (HO)(3)C-Me. The intermediate is unstable and is isomerized to a neutral interacting system, MeCOOH...(HO)(3)C-Me + H(2)O. The ethanetriol, (HO)(3)-CMe is transformed to an acetic acid monomer. The monomer may be dissociated to give a second zwitterionic intermediate with reasonable proton-relay patterns and energy changes. In proton relay reaction channels, H in MeCOOH is not an acidic proton but is always a hydroxy proton.     

Insights into hydrogen bonding and stacking interactions in cellulose

In this quantum chemical study, we explore hydrogen bonding (H-bonding) and stacking interactions in different crystalline cellulose allomorphs; namely, cellulose I(β) and cellulose III(I). We consider a model system representing a cellulose crystalline core made from six cellobiose units arranged in three layers with two chains per layer. We calculate the contributions of intrasheet and intersheet interactions to the structure and stability in both cellulose I(β) and cellulose III(I) crystalline cores. Reference structures for this study were generated from molecular dynamics simulations of water-solvated cellulose I(β) and III(I) fibrils. A systematic analysis of various conformations describing different mutual orientations of cellobiose units is performed using the hybrid density functional theory with the M06-2X with 6-31+G(d,p) basis sets. We dissect the nature of the forces that stabilize the cellulose I(β) and cellulose III(I) crystalline cores and quantify the relative strength of H-bonding and stacking interactions. Our calculations demonstrate that individual H-bonding interactions are stronger in cellulose I(β) than in cellulose III(I); however, the total H-bonding contribution to stabilization is larger in cellulose III(I) because of the highly cooperative nature of the H-bonding network. In addition, we observe a significant contribution from cooperative stacking interactions to the stabilization of cellulose I(β). The theory of atoms-in-molecules (AIM) has been employed to characterize and quantify these intermolecular interactions. AIM analyses highlight the role of nonconventional CH···O H-bonding in the cellulose assemblies. Finally, we calculate molecular electrostatic potential maps for the cellulose allomorphs that capture the differences in chemical reactivity of the systems considered in our study.     

Self-assembled aggregates originated from the balance of hydrogen-bonding, electrostatic, and hydrophobic interactions

Rich phase behavior was observed in salt-free cationic and anionic (catanionic) mixtures of a double-tailed surfactant, di(2-ethylhexyl)phosphoric acid (abbreviated as DEHPA), and tetradecyldimethylamine oxide (C(14)DMAO) in water. At a fixed C(14)DMAO concentration, phase transition from L(1) phase to L(α) phase occurs with increasing amounts of DEHPA. Moreover, in the L(α) phase, with the increase in DEHPA concentration, a gradual transition process from vesicle phase (L(αv)) to stacked lamellar phase (L(αl)) was determined by cryo- and FF-TEM observations combining with (2)H NMR measurements. The rheological data show that the viscosity increases with DEHPA amounts for L(αv) phase samples because of the increase in vesicle density. At a certain molar ratio of DEHPA to C(14)DMAO, i.e., 80:250, the samples are with the highest viscoelasticity, indicating the existence of densely packed vesicles. While for L(αl) phase samples, with increasing DEHPA amount, a decrease of bilayer curvature was induced, leading to a decrease of viscosity obviously. Compared with general catanionic surfactant mxitures, in addition to the electrostatic interaction of ion pairs, the transition of the microstructures is also ascribed to the formation of the hydrogen bonding (-N(+)-O-H···O-N-) between C(14)DMAO molecules and protonated C(14)DMAOH(+), which induces the growth of aggregates and the decrease of aggregate curvatures.     

Organometallic ruthenium(II) diamine anticancer complexes: arene-nucleobase stacking and stereospecific hydrogen-bonding in guanine adducts

Organometallic ruthenium(II) arene anticancer complexes of the type [(eta(6)-arene)Ru(II)(en)Cl][PF(6)] (en = ethylenediamine) specifically target guanine bases of DNA oligomers and form monofunctional adducts (Morris, R., et al. J. Med. Chem. 2001). We have determined the structures of monofunctional adducts of the "piano-stool" complexes [(eta(6)-Bip)Ru(II)(en)Cl][PF(6)] (1, Bip = biphenyl), [(eta(6)-THA)Ru(II)(en)Cl][PF(6)] (2, THA = 5,8,9,10-tetrahydroanthracene), and [(eta(6)-DHA)Ru(II)(en)Cl][PF(6)] (3, DHA = 9,10-dihydroanthracene) with guanine derivatives, in the solid state by X-ray crystallography, and in solution using 2D [(1)H,(1)H] NOESY and [(1)H,(15)N] HSQC NMR methods. Strong pi-pi arene-nucleobase stacking is present in the crystal structures of [(eta(6)-C(14)H(14))Ru(en)(9EtG-N7)][PF(6)](2).(MeOH) (6) and [(eta(6)-C(14)H(12))Ru(en)(9EtG-N7)][PF(6)](2).2(MeOH) (7) (9EtG = 9-ethylguanine). The anthracene outer ring (C) stacks over the purine base at distances of 3.45 A for 6 and 3.31 A for 7, with dihedral angles of 3.3 degrees and 3.1 degrees, respectively. In the crystal structure of [(eta(6)-biphenyl)Ru(en)(9EtG-N7)][PF(6)](2).(MeOH) (4), there is intermolecular stacking between the pendant phenyl ring and the purine six-membered ring at a distance of 4.0 A (dihedral angle 4.5 degrees). This stacking stabilizes a cyclic tetramer structure in the unit cell. The guanosine (Guo) adduct [(eta(6)-biphenyl)Ru(en)(Guo-N7)][PF(6)](2).3.75(H(2)O) (5) exhibits intramolecular stacking of the pendant phenyl ring with the purine five-membered ring (3.8 A, 23.8 degrees) and intermolecular stacking of the purine six-membered ring with an adjacent pendant phenyl ring (4.2 A, 23.0 degrees). These occur alternately giving a columnar-type structure. A syn orientation of arene and purine is present in the crystal structures 5, 6, and 7, while the orientation is anti for 4. However, in solution, a syn orientation predominates for all the biphenyl adducts 4, 5, and the guanosine 5'-monophosphate (5'-GMP) adduct 8 [(eta(6)-biphenyl)Ru(II)(en)(5'-GMP-N7)], as revealed by NMR NOE studies. The predominance of the syn orientation both in the solid state and in solution can be attributed to hydrophobic interactions between the arene and purine rings. There are significant reorientations and conformational changes of the arene ligands in [(eta(6)-arene)Ru(II)(en)(G-N7)] complexes in the solid state, with respect to those of the parent chloro-complexes [(eta(6)-arene)Ru(II)(en)Cl](+). The arene ligands have flexibility through rotation around the arene-Ru pi-bonds, propeller twisting for Bip, and hinge-bending for THA and DHA. Thus propeller twisting of Bip decreases by ca. 10 degrees so as to maximize intra- or intermolecular stacking with the purine ring, and stacking of THA and DHA with the purine is optimized when their tricyclic ring systems are bent by ca. 30 degrees, which involves increased bending of THA and a flattening of DHA. This flexibility makes simultaneous arene-base stacking and N7-covalent binding compatible. Strong stereospecific intramolecular H-bonding between an en NH proton oriented away from the arene (en NH(d)) and the C6 carbonyl of G (G O6) is present in the crystal structures of 4, 5, 6, and 7 (average N...O distance 2.8 A, N-H...O angle 163 degrees ). NMR studies of the 5'-GMP adduct 8 provided evidence that en NH(d) protons are involved in strong H-bonding with the 5'-phosphate and O6 of 5'-GMP. The strong H-bonding from G O6 to en NH(d) protons partly accounts for the high preference for binding of [(eta(6)-arene)Ru(II)en](2+) to G versus A (adenine). These studies suggest that simultaneous covalent coordination, intercalation, and stereospecific H-bonding can be incorporated into Ru(II) arene complexes to optimize their DNA recognition behavior, and as potential drug design features.