Theoretical description of hydrogen bonding in oxalic acid dimer and trimer based on the combined extended-transition-state energy decomposition analysis and natural orbitals for chemical valence (ETS-NOCV)

In the present study we have analyzed hydrogen bonding in dimer and trimer of oxalic acid, based on a recently proposed charge and energy decomposition scheme (ETS-NOCV). In the case of a dimer, two conformations, α and β, were considered. The deformation density contributions originating from NOCV’s revealed that the formation of hydrogen bonding is associated with the electronic charge deformation in both the σ—(Δρ_σ) and π-networks (Δρ_π). It was demonstrated that σ-donation is realized by electron transfer from the lone pair of oxygen on one monomer into the empty r_{H - O}^* \rho_{H - O}^* orbital of the second oxalic acid fragment. In addition, a covalent contribution is observed by the density transfer from hydrogen of H-O group in one oxalic acid monomer to the oxygen atom of the second fragment. The resonance assisted component (Δρ_π), is based on the transfer of electron density from the π—orbital localized on the oxygen of OH on one oxalic acid monomer to the oxygen atom of the other fragment. ETS-NOCV allowed to conclude that the σ(O---HO) component is roughly eight times as important as π (RAHB) contribution in terms of energetic estimation. The electrostatic factor (ΔE_elstat) is equally as important as orbital interaction term (ΔE_orb). Finally, comparing β-dimer of oxalic acid with trimer we found practically no difference concerning each of the O---HO bonds, neither qualitative nor quantitative.     

Probing the nature of hydrogen bonds in DNA base pairs

Energy decomposition analyses based on the block-localized wave-function (BLW-ED) method are conducted to explore the nature of the hydrogen bonds in DNA base pairs in terms of deformation, Heitler–London, polarization, electron-transfer and dispersion-energy terms, where the Heitler–London energy term is composed of electrostatic and Pauli-exchange interactions. A modest electron-transfer effect is found in the Watson–Crick adenine–thymine (AT), guanine–cytosine (GC) and Hoogsteen adenine-thymine (H-AT) pairs, confirming the weak covalence in the hydrogen bonds. The electrostatic attraction and polarization effects account for most of the binding energies, particularly in the GC pair. Both theoretical and experimental data show that the GC pair has a binding energy (−25.4 kcal mol⁻¹ at the MP2/6-31G** level) twice that of the AT (−12.4 kcal mol⁻¹) and H-AT (−12.8 kcal mol⁻¹) pairs, compared with three conventional N-H···O(N) hydrogen bonds in the GC pair and two in the AT or H-AT pair. Although the remarkably strong binding between the guanine and cytosine bases benefits from the opposite orientations of the dipole moments in these two bases assisted by the π-electron delocalization from the amine groups to the carbonyl groups, model calculations demonstrate that π-resonance has very limited influence on the covalence of the hydrogen bonds. Thus, the often adopted terminology “resonance-assisted hydrogen bonding (RHAB)” may be replaced with “resonance-assisted binding” which highlights the electrostatic rather than electron-transfer nature of the enhanced stabilization, as hydrogen bonds are usually regarded as weak covalent bonds. Figure Electron density difference (EDD) maps for the GC pair: a shows the polarization effect (isodensity 1.2×10⁻³ a.u.); b shows the charge transfer effect (isodensity 2×10⁻⁴ a.u.) Dedicated to Professor Paul von Ragué Schleyer on the occasion of his 75th birthday     

Noncovalent and covalent functionalization of a (5, 0) single-walled carbon nanotube with alanine and alanine radicals

We have systematically investigated the noncovalent and covalent adsorption of alanine and alanine radicals, respectively, onto a (5, 0) single-walled carbon nanotube using first-principles calculation. It was found that XH···π (X = N, O, C) interactions play a crucial role in the non-ovalent adsorption and that the functional group close to the carbon nanotube exhibits a significant influence on the binding strength. Noncovalent functionalization of the carbon nanotube with alanine enhances the conductivity of the metallic (5, 0) nanotube. In the covalent adsorption of each alanine radical onto a carbon nanotube, the binding energy depends on the adsorption site on CNT and the electronegative atom that binds with the CNT. The strongest complex is formed when the alanine radical interacts with a (5, 0) carbon nanotube through the amine group. In some cases, the covalent interaction of the alanine radical introduces a half-filled band at the Fermi level due to the local sp ³ hybridization, which modifies the conductivity of the tube.     

Thermodynamic and Kinetic Characterization of Duplex Formation between 2′-O, 4′-C-Methylene-modified Oligoribonucleotides,DNA and RNA

2′-O,4′-C-methylene-linked ribonucleotide derivatives, named LNA (locked nucleic acid) and BNA (bridged nucleic acid) are nucleic acid analogoues that have shown high-affinity recognition of DNA and RNA, and the employment of LNA oligomers for antisense activity, gene regulation and nucleic acid diagnostics seems promising. Here we show kinetic and thermodynamic results on the interaction of a series of 10 bases long LNA–DNA mixmers, gabmers as well as full length LNA’s with the complementary DNA, RNA and LNA oligonucleotides in the presence and absence of 10 mM Mg²⁺- ions. Our results show no significant differences in the reaction thermodynamics and kinetics between the LNA species, only a tendency to stronger duplex formation with the gabmer and mixmer. Introduction of a few LNA’s thus may be a better strategy, than using full length LNA’s to obtain an oligonucleotide that markedly increases the strength of duplexes formed with the complementary DNA and RNA.     

Molecular orbital study of the hydroxylation of benzene and monofluorobenzene catalysed by iron-oxo porphyrin π cation radical complexes

 The reaction mechanism for the hydroxylation of benzene and monofluorobenzene, catalysed by a ferryl-oxo porphyrin cation radical complex (compound) is described by electronic structure calculations in local spin density approximation. The active site of the enzyme is modelled as a six-coordinated (Por⁺)Fe(IV)O a_2u complex with imidazole or H₃CS^– as the axial ligand. The substrates under study are benzene and fluorobenzene, with the site of attack in para, meta and ortho position with respect to F. Two reaction pathways are investigated, with direct oxygen attack leading to a tetrahedral intermediate and arene oxide formation as a primary reaction step. The calculations show that the arene oxide pathway is distinctly less probable, that hydroxylation by an H₃CS^––coordinated complex is energetically favoured compared with imidazole, and that the para position with respect to F is the preferred site for hydroxylation. A partial electron transfer from the substrate to the porphyrin during the reaction is obtained in all cases. The resulting charge distribution and spin density of the substrates reveal the transition state as a combination of a cation and a radical σ-adduct intermediate with slightly more radical character in the case of H₃CS^– as axial ligand. A detailed analysis of the orbital interactions along the reaction pathway yields basically different mechanisms for the modes of substrate–porphyrin electron transfer and rupture of the Fe–O bond. In the imidazole-coordinated complex an antibonding π*(Fe–O) orbital is populated, whereas in the H₃CS^––coordinated system a shift of electron density occurs from the Fe–O bond region into the Fe–S bond. Received: 1 July 1995 / Accepted: 18 December 1995     

A hydrophobic similarity analysis of solvation effects on nucleic acid bases

We investigate the changes in the solvation properties of the natural nucleic acid bases due to the formation of the canonical Watson–Crick hydrogen-bonded complexes. To this end, the changes in the free energy of solvation of the bases induced upon hydrogen-bonded dimerization are analyzed by means of the hydrophobic similarity index, which relies on the atomic contributions to the free energy of solvation determined by the partitioning method implemented in the framework of the MST continuum model. Such an index is also used to examine the hydrophobic similarity between the canonical nucleic acid bases and a series of highly apolar analogues, which have been designed as potential candidates to expand the genetic alphabet. The ability of these analogues to be incorporated into modified DNA duplexes can be related to the large reduction in the hydrophilicity of the natural bases upon formation of the canonical hydrogen-bonded dimers. The results illustrate the suitability of the hydrophobic similarity index to rationalize the role played by solvation in molecular recognition. Proceedings of “Modeling Interactions in Biomolecules II”, Prague, September 5th–9th, 2005.     

Structural characterization of a carbon monoxide adduct of a heme–DNA complex

Abstract  

The structure of a carbon monoxide (CO) adduct of a complex between heme and a parallel G-quadruplex DNA formed from a single repeat sequence of the human telomere, d(TTAGGG), has been characterized using ¹H and ¹³C NMR spectroscopy and density function theory calculations. The study revealed that the heme binds to the 3′-terminal G-quartet of the DNA though a π–π stacking interaction between the porphyrin moiety of the heme and the G-quartet. The π–π stacking interaction between the pseudo-C ₂-symmetric heme and the C ₄-symmetric G-quartet in the complex resulted in the formation of two isomers possessing heme orientations differing by 180° rotation about the pseudo-C ₂ axis with respect to the DNA. These two slowly interconverting heme orientational isomers were formed in a ratio of approximately 1:1, reflecting that their thermodynamic stabilities are identical. Exogenous CO is coordinated to heme Fe on the side of the heme opposite the G-quartet in the complex, and the nature of the Fe–CO bond in the complex is similar to that of the Fe–CO bonds in hemoproteins. These findings provide novel insights for the design of novel DNA enzymes possessing metalloporphyrins as prosthetic groups.     

Topological analysis of aromatic halogen/hydrogen bonds by electron charge density and electrostatic potentials

In this work, the intermolecular distribution of the electronic charge density in the aromatic hydrogen/halogen bonds is studied within the framework of the atoms in molecules (AIM) theory and the molecular electrostatic potentials (MEP) analysis. The study is carried out in nine complexes formed between benzene and simple lineal molecules, where hydrogen, fluorine and chlorine atoms act as bridge atoms. All the results are obtained at MP2 level theory using cc-pVTZ basis set. Attention is focused on topological features observed at the intermolecular region such as bond, ring and cage critical points of the electron density, as well as the bond path, the gradient of the density maps, molecular graphs and interatomic surfaces. The strength of the interaction increases in the following order: F⋅⋅⋅π < Cl⋅⋅⋅π < H⋅⋅⋅π. Our results show that the fluorine atom has the capability to interact with the π−cloud to form an aromatic halogen bond, as long as the donor group is highly electron withdrawing. The Laplacian topology allows us to state that the halogen atoms can act as nucleophiles as well as electrophiles, showing clearly their dual character.     

A theoretical study on unusual intermolecular T-shaped X–H...π interactions between the singlet state HB=BH and HF, HCl, HCN or H2C2

The unusual T-shaped X–H...π hydrogen bonds are found between the B=B double bond of the singlet state HB=BH and the acid hydrogen of HF, HCl, HCN and H₂C₂ using MP2 and B3LYP methods at 6-311++G(2df,2p) and aug-cc-pVTZ levels. The binding energies follow the order of HB=BH...HF>HB=BH...HCl>HB=BH...HCN>HB=BH...H₂C₂. The hydrogen-bonded interactions in HB=BH...HX are found to be stronger than those in H₂C=CH₂...HX and OCB≡BCO...HX. The analyses of natural bond orbital (NBO) and the electron density shifts reveal that the nature of the T-shaped X–H...π hydrogen-bonded interaction is that much of the lost density from the π-orbital of B=B bond is shifted toward the hydrogen atom of the proton donor, leading to the electron density accumulation and the formation of the hydrogen bond. The atoms in molecules (AIM) theory have also been applied to characterize bond critical points and confirm that the B=B double bond can be a potential proton acceptor. The unusual T-shaped X–H...π hydrogen bonds are found between the B=B double bond of the singlet state HB=BH and the acid hydrogen of HF, HCl, HCN and H₂C₂     

Structural stability studies in adhesion molecules—role of cation–<Emphasis Type="Italic">π</Emphasis> interactions

Cell adhesion molecules are important for their various roles in many cellular events and responses. In the present study, we have analyzed the roles played by cation–π interactions in the structural stability of adhesion molecules. These interactions are mainly formed by long-range contacts. The occurrence of arginine is higher than lysine to form cation–π interactions. The secondary structure preferences of interacting residues are independent of amino acid class. Cation–π interactions might stabilize the interface between the terminus and core in this class of proteins. The results obtained in the present study will be useful in understanding the contribution of cation–π interactions to the overall stability of adhesion proteins.     

Electronic and molecular structure of M-DNA fragments

To study M-DNA molecular structure (such DNA with transition metal ions placed between the nucleic bases is able to conduct the electric current) and its conductivity mechanisms, we carried out ab initio quantum-mechanical calculations of electronic and spatial structures, thermodynamic characteristics of adenine-thymine (АТ) and guanine-cytosine (GC) base pair complexes with Zn²⁺ and Ni²⁺. To take into account the influence of the alkaline environment, calculations for these complexes were also carried out with hydroxyl and two water molecules. Computations were performed at MP2 level of theory using 6–31+G* basis set. Analogous calculations were carried out for (AC)(TG) stacking dimer of nucleic acid base pairs with two Zn²⁺. The calculation of the interaction energy in complexes has shown the preference of locating the metal ion (instead of the imino proton) between bases in M-DNA. The electronic transition energy calculation has revealed the reduction of the first singlet transition energy in АТ and GC complexes with Ni²⁺ from 4.5 eV to 0.4 - 0.6 eV. Ni²⁺ orbitals take part in the formation of HOMO and LUMO on the complexes investigated. It was shown that charges of metal ions incorporated into complexes with nucleic bases and in dimer decrease significantly.     

A novel nanohybrid of daunomycin and single-walled carbon nanotubes: photophysical properties and enhanced electrochemical activity

A nanohybrid adduct of the widely used, functional dye, daunomycin (DM), with single-walled carbon nanotubes (SWNTs) was prepared. Ultraviolet-visible-near infrared and fluorescence spectroscopy and electrochemistry of DM-functionalized SWNTs reveal that DM interacts with SWNTs through strong π–π stacking and there is a significant photo-induced charge-transfer interaction between the two components. Importantly, the novel adduct modified the glassy carbon (GC) electrode to give a much enhanced electrochemical activity than those of DM adsorbed onto not only the bare GC electrode but also the SWNTs-modified GC electrode.     

Theoretical investigations of the H···π and X (X = F,Cl, Br,I)···π complexes between hypohalous acids and benzene

The H···π and X (X = F, Cl, Br, I)···π interactions between hypohalous acids and benzene are investigated at the MP2/6-311++G(2d,2p) level. Four hydrogen-bonded and three halogen-bonded complexes were obtained. Ab initio calculations indicate that the X···π interaction between HOX and C₆H₆ is mainly electrostatically driven, and there is nearly an equal contribution from both electrostatic and dispersive energies in the case of XOH–C₆H₆ complexes. Natural bond orbital (NBO) analysis reveals that there exists charge transfer from benzene to hypohalous acids. Atom in molecules (AIM) analysis locates bond critical points (BCP) linking the hydrogen or halogen atom and carbon atom in benzene.     

<Superscript>13</Superscript>C-detected NMR experiments for measuring chemical shifts and coupling constants in nucleic acid bases

The paper presents a set of two-dimensional experiments that utilize direct ¹³C detection to provide proton–carbon, carbon–carbon and carbon–nitrogen correlations in the bases of nucleic acids. The set includes a ¹³C-detected proton–carbon correlation experiment for the measurement of ¹³C–¹³C couplings, the CaCb experiment for correlating two quaternary carbons, the HCaCb experiment for the ¹³C–¹³C correlations in cases where one of the carbons has a proton attached, the HCC-TOCSY experiment for correlating a proton with a network of coupled carbons, and a ¹³C-detected ¹³C–¹⁵N correlation experiment for detecting the nitrogen nuclei that cannot be detected via protons. The IPAP procedure is used for extracting the carbon–carbon couplings and/or carbon decoupling in the direct dimension, while the S³E procedure is preferred in the indirect dimension of the carbon–nitrogen experiment to obtain the value of the coupling constant. The experiments supply accurate values of ¹³C and ¹⁵N chemical shifts and carbon–carbon and carbon–nitrogen coupling constants. These values can help to reveal structural features of nucleic acids either directly or via induced changes when the sample is dissolved in oriented media. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Computational (DFT and TD DFT) study of the electron structure of the tautomers/conformers of uridine and deoxyuridine and the processes of intramolecular proton transfers

Six uridine and six deoxyuridine isomers were studied at the B3LYP and TD B3LYP theoretical level and 6–31+G(d) basis function. The stability and the excited states of the isomers were studied in order to clarify some known experimental data. It was established that the rotation of the oxo uracil ring in uridine is energetically more likely to occur in the excited state than in the ground state, driven by the bright ¹ ππ* state and the dark charge transfer ¹nπ* state. Very high energy barriers (on the S_o) were found for thermal intramolecular proton transfer processes.     

Supramolecular synthon pattern in solid clioquinol and cloxiquine (APIs of antibacterial,antifungal, antiaging and antituberculosis drugs) studied by 35Cl NQR, 1H-17O and 1H-14N NQDR and DFT/QTAIM

The quinolinol derivatives clioquinol (5-chloro-7-iodo-8-quinolinol, Quinoform) and cloxiquine (5-chloro-8-quinolinol) were studied experimentally in the solid state via ³⁵Cl NQR, ¹H-¹⁷O and ¹H-¹⁴N NQDR spectroscopies, and theoretically by density functional theory (DFT). The supramolecular synthon pattern of O–H···N hydrogen bonds linking dimers and π–π stacking interactions were described within the QTAIM (quantum theory of atoms in molecules) /DFT (density functional theory) formalism. Both proton donor and acceptor sites in O–H···N bonds were characterized using ¹H-¹⁷O and ¹H-¹⁴N NQDR spectroscopies and QTAIM. The possibility of the existence of O–H···H–O dihydrogen bonds was excluded. The weak intermolecular interactions in the crystals of clioquinol and cloxiquine were detected and examined. The results obtained in this work suggest that considerable differences in the NQR parameters for the planar and twisted supramolecular synthons permit differentiation between specific polymorphic forms, and indicate that the more planar supramolecular synthons are accompanied by a greater number of weaker hydrogen bonds linking them and stronger π···π stacking interactions.     

<Emphasis Type="Italic">π</Emphasis>–<Emphasis Type="Italic">π</Emphasis> interaction between aromatic ring and copper-coordinated His81 imidazole regulates the blue copper active-site structure

Noncovalent weak interactions play important roles in biological systems. In particular, such interactions in the second coordination shell of metal ions in proteins may modulate the structure and reactivity of the metal ion site in functionally significant ways. Recently, π–π interactions between metal ion coordinated histidine imidazoles and aromatic amino acids have been recognized as potentially important contributors to the properties of metal ion sites. In this paper we demonstrate that in pseudoazurin (a blue copper protein) the π–π interaction between a coordinated histidine imidazole ring and the side chains of aromatic amino acids in the second coordination sphere, significantly influences the properties of the blue copper site. Electronic absorption and electron paramagnetic resonance spectra indicate that the blue copper electronic structure is perturbed, as is the redox potential, by the introduction of a second coordination shell π–π interaction. We suggest that the π–π interaction with the metal ion coordinated histidine imidazole ring modulates the electron delocalization in the active site, and that such interactions may be functionally important in refining the reactivity of blue copper sites. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

The effect of anchoring group number on molecular structures and absorption spectra of triphenylamine sensitizers: a computational study

The molecular structures and absorption spectra of triphenylamine dyes containing different numbers of anchoring groups (S1-S3) were investigated by density functional theory (DFT) and time-dependent DFT. The calculated geometries indicate that strong conjugation is formed in the dyes. The interfacial charge transfer between the TiO₂ electrode and S1-S3 are electron injection processes from the excited dyes to the semiconductor conduction band. The simulated absorption bands are assigned to π → π* transitions according to the qualitative agreement between the experimental and calculated results. The effect of anchoring group number on the molecular structures, absorption spectra and photovoltaic performance were comparatively discussed.     

XH/<Emphasis Type="Italic">π</Emphasis> interactions with the <Emphasis Type="Italic">π</Emphasis> system of porphyrin ring in porphyrin-containing proteins

Searching structures of porphyrin-containing proteins from the Protein Data Bank revealed that the π system of every porphyrin ring is involved in XH/π interactions, with most of the porphyrins having several interactions. Both five-membered pyrrole rings and six-membered chelate rings are involved in XH/π interactions; the number of interactions with five-membered rings is larger than the number of interactions with six-membered rings. We found interactions with C–H and N–H groups as hydrogen-atom donors; however, the number of CH/π interactions is much larger than the number of NH/π interactions. The amino acids involved in the interactions show a high conservation score. Our results that every porphyrin is involved in XH/π interactions and that amino acids involved in these interactions are highly conserved demonstrate that XH/π interactions play an important role in porphyrin–protein stability. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.