Na+, Mg2+, and Zn2+ Binding to All Tautomers of Adenine, Cytosine, and Thymine and the Eight Most Stable Keto/Enol Tautomers of Guanine: A Correlated ab Initio Quantum Chemical Study

Cryogenic Ion Spectroscopy of Protonated and Sodiated Methyladenine Derivatives

Ultraviolet photodissociation (UVPD) spectra of protonated 9-methyladenine (H⁺9MA), protonated 7-methyl adenine (H⁺7MA), protonated 3-methyladenine (H⁺3MA), and sodiated 7-methyladenine (Na⁺7MA) near the origin bands of the S₀-S₁ transition were obtained using cryogenic ion spectroscopy. The UV-UV hole burning, infrared (IR) ion-dip, and IR-UV double resonance spectra showed that all the ions were present as single isomers in a cryogenic ion trap. The UVPD spectrum of H⁺9MA exhibited only a broad absorption band, whereas the spectra of H⁺7MA, H⁺3MA, and Na⁺7MA displayed moderately or well-resolved vibronic bands. Potential energy profiles were computed to understand the reason for the different bandwidths of the vibronic bands in the spectra. The broadening of the bands was correlated with the slopes between the Franck-Condon point and the conical intersection between the S₁ and S₀ states in the potential energy profiles, thus reflecting the deactivation rates in the S₁ state.

Cryogenic ion spectroscopy of adenine complexes containing alkali metal cations

Cryogenic ion spectroscopy was used to characterize adenine complexes containing alkali metal cations (M⁺A, M = Cs, Rb, K, Na, and Li) produced by electrospray ionization. The ultraviolet (UV) photodissociation spectra of the complexes stored in a cryogenic ion trap exhibited well-resolved vibronic bands near their origin bands of the S₀-S₁ transition. The UV-UV hole-burning and infrared ion-dip spectra showed that all the M⁺A ions in the ion trap were single isomers of M⁺A7a, where the M⁺ ion was not bound to canonical 9H-adenine (A9) but bound to a rare tautomer, 7H-adenine (A7). Density functional theory calculations showed lower tautomerization barriers for M⁺A9 than for bare A9 in aqueous solution. We suggest that M⁺ ions not only play a catalytic role in the tautomerization of A9 to A7 but also increase the tautomerization yield by forming stable M⁺A7a isomers.

A TD-DFT study of the excited dissociative electronic states of the DNA nucleobases bound to Li. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2021;244:118862. [PMID: 32927359 DOI: 10.1016/j.saa.2020.118862]

In this work, the first and second dissociative potential curves of adenineLi⁺ (ADLi⁺), guanineLi⁺ (GUALi⁺), cytosineLi⁺(CYTLi⁺), and thymineLi⁺ (THYLi⁺) complexes, related to the dissociation of their LiO and LiN bonds, have been calculated in the gas phase and water, separately. For this purpose, the fifteen excited potential curves, in order of increasing energy, were calculated for each complex and the dissociative potential curves were distinguished from them considering the conical intersection points. The time-dependent density functional theory (TD-DFT) method employing the M06-2X functional was used for the calculations. It was observed that the electron transfer from the DNA base to the Li⁺ took place during the dissociation of complexes in the gas phase. The electrostatic field of water blocked this charge transfer and led to the excited DNA base and Li⁺ in its ground state. The vertical excitation energy for the desorption of the Li fragment as a neutral and cation species from the DNA bases was determined. The effect of the interaction site of Li⁺ on the dissociative potential curves was also investigated.

A theoretical study on the role of stability of cytosine and its tautomers in DNA (deoxyribonucleic acid), and investigation of interactions of Na+, K+, Mg2+, Ca2+, Zn2+ metal ions and OH radical with cytosine tautomers

In the present study, 21 cytosine tautomers were investigated so that some tautomers were reported for the first time in the gas phase and aqueous solution. C₃ tautomer was the most stable tautomer in gas phase but C₁ was the most stable structure in aqueous solution. The potential energy surface of all trajectories was determined for 21 tautomers and 22 transition states. Also, interactions of cytosine tautomers with Na⁺, K⁺, Mg²⁺, Ca²⁺ and Zn²⁺ metal ions were studied in gas phase and aqueous solution. Three types of interactions among metal ions and (N1 and O10), (N3 and O10) and (N3 and N9) of cytosine tautomers were investigated. The study of interaction energies of all complexes showed the stability of complexes in which interactions among Mg²⁺ and Zn²⁺ with tautomers were stronger than interactions among Ca²⁺, Na⁺ and K⁺ with tautomers, respectively. Some interactions of metal ions with cytosine tautomers made the most stable tautomers. So, the stability of rare tutomeric forms had a significant effect on stabilization of anomalous DNA (deoxyribonucleic acid) double helix and spontaneous mutations. Also, one of the most important causes of mutations in DNA (deoxyribonucleic acid) was the reaction of OH radical with nucleotide bases. So, interactions of OH radical with cytosine and its tautomers were investigated in gas phase and aqueous solution.Communicated by Ramaswamy H. Sarma.

Effect of Alkali Metal Cations on Length and Strength of Hydrogen Bonds in DNA Base Pairs

For many years, non-covalently bonded complexes of nucleobases have attracted considerable interest. However, there is a lack of information about the nature of hydrogen bonding between nucleobases when the bonding is affected by metal coordination to one of the nucleobases, and how the individual hydrogen bonds and aromaticity of nucleobases respond to the presence of the metal cation. Here we report a DFT computational study of nucleobase pairs interacting with alkali metal cations. The metal cations contribute to the stabilization of the base pairs to varying degrees depending on their position. The energy decomposition analysis revealed that the nature of bonding between nucleobases does not change much upon metal coordination. The effect of the cations on individual hydrogen bonds were described by changes in VDD charges on frontier atoms, H-bond length, bond energy from NBO analysis, and the delocalization index from QTAIM calculations. The aromaticity changes were determined by a HOMA index.

Role of backbones on the interaction of metal ions with deoxyribonucleic acid and peptide nucleic acid: A DFT study

Metal ion interaction with deoxyribonucleic acid and peptide nucleic acid were studied using B3LYP-D3/6-311++g(d,p)//B3LYP/6-31 + G(d) level of theory in aqueous phase employing polarized continuum (PCM) model. This study reports the role of backbones on deoxyribonucleic acid and peptide nucleic acid for complexation with different metal ions. The systematic study performed with DFT calculations reveals that central binding (Type-4) shows the strongest binding compared to the other binding modes because of the involvement of the backbone as well as the nitrogenous bases. The charged backbone of DNA nucleotides contributes significantly towards binding with the metal ions. The deoxyguanosine monophosphate (dGMP) clearly indicates the strongest binding upon complexation with Mg²⁺ (-49.6 kcal/mol), Zn²⁺ (-45.3 kcal/mol) and Cu²⁺ (-148.4 kcal/mol), respectively. The neutral backbone of PNA also assists to complex the metal ions with PNA nucleotides. The Mg²⁺ and Cu²⁺ prefer to bind with the PNA-Cytosine (-32.9 kcal/mol & -132.9 kcal/mol) in central binding mode (type-4). PNA-Adenine-Zn²⁺ (-29.1 kcal/mol) is the preferred binding mode (type-4) compared to other modes of interaction for this metal ion with PNA-Adenine nucleotide. The Cu²⁺ ion showed the superior complexation ability with deoxyribonucleic acid and peptide nucleic acid compared to Mg²⁺ and Zn²⁺ ions. The cation-π complexation with the bases of nucleotides was also obtained with Cu²⁺ ion. The AIM (atoms in molecule) theory has been applied to examine the nature of the interaction of Mg²⁺, Zn²⁺, and Cu²⁺ ion to the deoxyribonucleic acid and peptide nucleic acid. The alkaline earth metal, Mg²⁺ ion shows electrostatic nature while interaction with deoxyribonucleic acid and peptide nucleic acid, however, the transition metal ions (Zn²⁺, Cu²⁺) showed partly covalent nature as well with deoxyribonucleic acid and peptide nucleic acid. The optical properties calculated for the binding of metal ions with deoxyribonucleic acid and peptide nucleic acid showed a diagnostic signature to ascertain the interaction of metal ions with such nucleotides. Cu²⁺ ion showed larger red shifts in the absorption spectrum values upon complexation with the DNAs and PNAs. The calculated results suggest that such metal ions would prefer to bind with the DNA compared to PNA in DNA-PNA duplexes. The preference for the binding of metal ions with DNA nucleotides is largely attributed to the contribution of charged backbones compared to the neutral PNA backbones.

Hydrogen bonding in alkali metal cation-bound i-motif-like dimers of 1-methyl cytosine: an IRMPD spectroscopic and computational study

The structures of alkali metal cation bound 1-methylcytosine (1-mCyt) dimers were explored using vibrational spectroscopy in the form of infrared multiple photon dissociation (IRMPD) spectroscopy and by computational methods. For the smaller alkali metal cations, Li⁺ and Na⁺, only non-hydrogen bonded symmetric anti-parallel structures were observed in agreement with the lowest energy computed structures. For K⁺, Rb⁺, and Cs⁺ the vibrational spectra in the N-H stretch region showed strong evidence for hydrogen bonding in agreement with the lowest energy structures which contained hydrogen bonding interactions between the amine group of one cytosine and the carbonyl oxygen of the other cytosine. The lowest energy structures for these complexes were compared to previously studied cytosine complexes [(Cyt)₂M]⁺ where M = Li, Na, and K. The calculations are in agreement that only the non-hydrogen bonded structures would be observed for these cytosine complexes.

Two-dimensional self-assembled nanostructures of nucleobases and their related derivatives on Au(111)

The construction of two-dimensional (2D) self-assembled nanostructures has been one of the considerably interesting areas of on-surface chemistry in the past few decades, and has benefited from the rapid development and improvement of scanning probe microscopy techniques. In this research field, many attempts have been made in the controllable fabrication of well-ordered and multifunctional surface nanostructures, which attracted interest because of the prospect for artificial design of functional molecular nanodevices. DNA and RNA are considered to be programmable self-assembly systems and it is possible to use their base sequences to encode instructions for assembly in a predetermined fashion at the nanometer scale. As important constituents of nucleic acids, nucleobases, with intrinsic functional groups for hydrogen bonding, coordination bonding, and electrostatic interactions, can be employed as a potential system for the versatile construction of various biomolecular nanostructures, which may be used to structure the self-assembly of DNA-based artificial molecular constructions and play an important role in novel biosensors based on surface functionalization. In this article, we will review the recent progress of on-surface self-assembly of nucleobases and their derivatives together with different reactants (e.g., metals, halogens, salts and water), and as a result, various 2D surface nanostructures are summarized.

The influence of isolated and penta-hydrated Zn2+ on some of the intramolecular proton-transfer processes of thymine: a quantum chemical study

Zinc cation (Zn2+) plays an important role in the chemistry of DNA base pairs. In this work, the influence of isolated and penta-hydrated Zn2+ on some of the intramolecular proton-transfer processes of thymine (T) is investigated by the density functional theory method. It is shown that the calculated binding energies between Zn2+ and T are exothermic in vacuum. Compared to T, Zn2+ increases the stability of tautomer T' by 28.7 kcal mol-1, promoting the intramolecular proton transfer of T. But in a micro-water environment, the attachment processes of Zn2+ to T hydrates, penta-hydrated Zn2+ to T, and penta-hydrated Zn2+ to T hydrates lead to the rearrangement of molecules and the redistribution of charges. The conventional T is still the most stable form and the influence of Zn2+ is much reduced and the proton transfer is thermodynamically unfavored. The detailed characterization is helpful to understand the genotoxicity of zinc ions.

Self-assembled uracil complexes containing tautomeric uracils: an IRMPD spectroscopic and computation study of the structures of gaseous uracilnCa2+ (n = 4, 5, or 6) complexes

The structures of doubly-charged uracil (U) complexes with Ca²⁺, U_nCa²⁺ (n = 4, 5, 6), were studied by infrared multiphoton dissociation (IRMPD) spectroscopy and computational methods.

Do the fragments from decomposed ZIF-8 greatly affect some of the intramolecular proton-transfer of thymine? A quantum chemical study

The intramolecular proton-transfer processes of thymine were investigated by the density functional theory method. It is shown that the mutation from keto (T) to enol (T′) form is affected by zeolitic imidazolate framework-8 (ZIF-8) fragments such as single 2-methylimidazole neutral crystals (M), and negatively charged 2-methylimidazole ligands (M⁻). Results show that with the number (n) of water (w) molecules that assist proton-transfer increasing from 1 to 4, the order of the tautomeric energy barriers (in kcal mol⁻¹) is T-2w (16.3) < T-1w (17.6) < T-3w (17.8) < T-4w (20.5). In the presence of M, the order of energy barrier is MT-2w (16.6) < MT-1w (17.7) < MT-3w (18.9) < MT-4w (20.8). M⁻ has a catalysis effect on the energy barrier and the order is M⁻T-2w (14.4) < M⁻T-3w (15.2) < M⁻T-1w (16.3) < M⁻T-4w (16.8). The attachment of the M⁻ fragment slightly promotes the proton-transfer processes in some instances. The characterization of the proton-transfer processes is helpful to understand the genotoxicity of ZIF-8 during drug delivery applications.

Investigations on whether fragments from decomposed ZIF-8 would affect the intramolecular proton-transfer of thymine by DFT modeling.

IRMPD Action Spectroscopy, ER-CID Experiments, and Theoretical Studies of Sodium Cationized Thymidine and 5-Methyluridine: Kinetic Trapping During the ESI Desolvation Process Preserves the Solution Structure of [Thd+Na]<sup/>

Thymidine (dThd) is a fundamental building block of DNA nucleic acids, whereas 5-methyluridine (Thd) is a common modified nucleoside found in tRNA. In order to determine the conformations of the sodium cationized thymine nucleosides [dThd+Na]⁺ and [Thd+Na]⁺ produced by electrospray ionization, their infrared multiple photon dissociation (IRMPD) action spectra are measured. Complementary electronic structure calculations are performed to determine the stable low-energy conformations of these complexes. Geometry optimizations and frequency analyses are performed at the B3LYP/6-311+G(d,p) level of theory, whereas energies are calculated at the B3LYP/6-311+G(2d,2p) level of theory. As protonation preferentially stabilizes minor tautomers of dThd and Thd, tautomerization facilitated by Na⁺ binding is also considered. Comparisons of the measured IRMPD and computed IR spectra find that [dThd+Na]⁺ prefers tridentate (O2,O4',O5') coordination to the canonical 2,4-diketo form of dThd with thymine in a syn orientation. In contrast, [Thd+Na]⁺ prefers bidentate (O2,O2') coordination to the canonical 2,4-diketo tautomer of Thd with thymine in an anti orientation. Although 2,4-dihydroxy tautomers and O2 protonated thymine nucleosides coexist in the gas phase, no evidence for minor tautomers is observed for the sodium cationized species. Consistent with experimental observations, the computational results confirm that the sodium cationized thymine nucleosides exhibit a strong preference for the canonical form of the thymine nucleobase. Survival yield analyses based on energy-resolved collision-induced dissociation (ER-CID) experiments suggest that the relative stabilities of protonated and sodium cationized dThd and Thd follow the order [dThd+H]⁺ < [Thd+H]⁺ < [dThd+Na]⁺ < [Thd+Na]⁺. Graphical Abstract ᅟ.

Influence of Sodium Cationization versus Protonation on the Gas-Phase Conformations and Glycosidic Bond Stabilities of 2'-Deoxyadenosine and Adenosine

The influence of noncovalent interactions with a sodium cation on the gas-phase structures and N-glycosidic bond stabilities of 2'-deoxyadenosine (dAdo) and adenosine (Ado), [dAdo+Na](+) and [Ado+Na](+), are probed via infrared multiple photon dissociation (IRMPD) action spectroscopy and energy-resolved collision-induced dissociation (ER-CID) experiments. ER-CID experiments are also performed on the protonated forms of these nucleosides, [dAdo+H](+) and [Ado+H](+), for comparison purposes. Complementary electronic structure calculations are performed to determine the structures and relative stabilities of the stable low-energy conformations of the sodium cationized nucleoside complexes and to predict their IR spectra. Comparison between the measured IRMPD action spectra and calculated IR spectra enables the conformations of the sodium cationized nucleosides present in the experiments to be elucidated. The influence of sodium cationization versus protonation on the structures and IR spectra is elucidated by comparison to IRMPD and theoretical results previously reported for the protonated forms of these nucleosides. The influence of sodium cationization versus protonation on the glycosidic bond stability of the adenine nucleosides is determined by comparison of the ER-CID behavior of these systems. All structures present in the experiments are found to involve tridentate binding of Na(+) to the N3, O4', and O5' atoms forming favorable 5- and 6-membered chelation rings, which requires that adenine rotate to a syn configuration. This mode of sodium cation binding results in moderate flexibility of the sugar moiety such that the sugar puckering of the conformations present varies between C2'-endo and O4'-endo. Sodium cationization is found to be less effective toward activating the N-glycosidic bond than protonation for both dAdo and Ado. Both the IRMPD yields and ER-CID behavior indicate that the 2'-hydroxyl substituent of Ado stabilizes the N-glycosidic bond relative to that of dAdo.

Interaction of Cu(+) with cytosine and formation of i-motif-like C-M(+)-C complexes: alkali versus coinage metals

The Watson-Crick structure of DNA is among the most well-known molecular structures of our time. However, alternative base-pairing motifs are also known to occur, often depending on base sequence, pH, or the presence of cations. Pairing of cytosine (C) bases induced by the sharing of a single proton (C-H(+)-C) may give rise to the so-called i-motif, which occurs primarily in expanded trinucleotide repeats and the telomeric region of DNA, particularly at low pH. At physiological pH, silver cations were recently found to stabilize C dimers in a C-Ag(+)-C structure analogous to the hemiprotonated C-dimer. Here we use infrared ion spectroscopy in combination with density functional theory calculations at the B3LYP/6-311G+(2df,2p) level to show that copper in the 1+ oxidation state induces an analogous formation of C-Cu(+)-C structures. In contrast to protons and these transition metal ions, alkali metal ions induce a different dimer structure, where each ligand coordinates the alkali metal ion in a bidentate fashion in which the N3 and O2 atoms of both cytosine ligands coordinate to the metal ion, sacrificing hydrogen-bonding interactions between the ligands for improved chelation of the metal cation.

Zn2+ blocks annealing of complementary single-stranded DNA in a sequence-selective manner

Zinc is the second most abundant trace element essential for all living organisms. In human body, 30–40% of the total zinc ion (Zn²⁺) is localized in the nucleus. Intranuclear free Zn²⁺ sparks caused by reactive oxygen species have been observed in eukaryotic cells, but question if these free Zn²⁺ outrages could have affected annealing of complementary single-stranded (ss) DNA, a crucial step in DNA synthesis, repair and recombination, has never been raised. Here the author reports that Zn²⁺ blocks annealing of complementary ssDNA in a sequence-selective manner under near-physiological conditions as demonstrated in vitro using a low-temperature EDTA-free agarose gel electrophoresis (LTEAGE) procedure. Specifically, it is shown that Zn²⁺ does not block annealing of repetitive DNA sequences lacking CG/GC sites that are the major components of junk DNA. It is also demonstrated that Zn²⁺ blocks end-joining of double-stranded (ds) DNA fragments with 3′ overhangs mimicking double-strand breaks, and prevents renaturation of long stretches (>1 kb) of denatured dsDNA, in which Zn²⁺-tolerant intronic DNA provides annealing protection on otherwise Zn²⁺-sensitive coding DNA. These findings raise a challenging hypothesis that Zn²⁺-ssDNA interaction might be among natural forces driving eukaryotic genomes to maintain the Zn²⁺-tolerant repetitive DNA for adapting to the Zn²⁺-rich nucleus.

An experimental and theoretical vibrational study of interaction of adenine and thymine with artificial seawaters: A prebiotic chemistry experiment

Nucleic acid bases play important roles in living beings. Thus, their interaction with salts the prebiotic Earth could be an important issue for the understanding of origin of life. In this study, the effect of pH and artificial seawaters on the structure of adenine and thymine was studied via parallel determinations using FT-IR, Raman spectroscopy and theoretical calculations. Thymine and adenine lyophilized in solutions at basic and acidic conditions showed characteristic bands of the enol-imino tautomer due to the deprotonation and the hydrochloride form due to protonation, respectively. The interaction of thymine and adenine with different seawaters representative of different geological periods on Earth was also studied. In the case of thymine a strong interaction with Sr(2+) promoted changes in the Raman and infrared spectra. For adenine changes in infrared and Raman spectra were observed in the presence of salts from all seawaters tested. The experimental results were compared to theoretical calculations, which showed structural changes due to the presence of ions Na(+), Mg(2+), Ca(2+) and Sr(2+) of artificial seawaters. For thymine the bands arising from C4=C5 and C6=O stretching were shifted to lower values, and for adenine, a new band at 1310cm(-1) was observed. The reactivity of adenine and thymine was studied by comparing changes in nucleophilicity and energy of the HOMO orbital.

Atomic-scale investigation on the facilitation and inhibition of guanine tautomerization at Au(111) surface

Nucleobase tautomerization might induce mismatch of base pairing. Metals, involved in many important biophysical processes, have been theoretically proven to be capable of affecting tautomeric equilibria and stabilities of different nucleobase tautomers. However, direct real-space evidence on demonstrating different nucleobase tautomers and further revealing the effect of metals on their tautomerization at surfaces has not been reported to date. From the interplay of high-resolution STM imaging and DFT calculations, we show for the first time that tautomerization of guanine from G/9H to G/7H is facilitated on Au(111) by heating, whereas such tautomerization process is effectively inhibited by introducing Ni atoms due to its preferential coordination at the N7 site of G/9H. These findings may help to elucidate possible influence of metals on nucleobase tautomerization and provide from a molecular level some theoretical basis on metal-based drug design.

Tautomerisation of thymine acts against the Hückel 4N + 2 rule. The effect of metal ions and H-bond complexations on the electronic structure of thymine

Not necessarily the π-electron delocalization is responsible for the stability of thymine tautomers.

INFLUENCE OF CATION-HETEROATOM (Li+, Na+, AND K+) INTERACTION ON THE STRUCTURAL AND THERMOCHEMICAL PROPERTIES OF 2′-DEOXYTHYMIDINE NUCLEOSIDE: QTAIM AND NBO ANALYZES

Density functional theory (DFT) method and B3LYP/6-311++G(d,p) basis set were used to determine coordination geometries, binding strength, and metal ion affinity (MIA) for interaction of 2′-deoxythymidine (dT) with alkali metal cations including Li+, Na+ , and K+ . Calculations demonstrated that the interaction of dT with these cations is tri-coordinated η (O2, O4′, O5′). Among these cations, Li+ cation exhibited the most tendency for interaction with dT. Cations via their interaction with dT can affect the N-glycosidic bond length, the values of pseudorotation of the sugar ring, the orientation of base unit with respect to sugar ring and the acidity of O5′H, O3′H, and N3H groups in 2′-dT nucleoside. Natural bond orbital (NBO) analysis was performed to calculate the charge transfer and natural population analysis of the complexes. Quantum theory of atoms in molecules (QTAIM) was also applied to determine the nature of interactions. It was shown that in these complexes, (dT- Li+ , dT- Na+ , and dT- K+ ), the bonds are an electrostatic (closed-shell) interaction in the nature.

Rigorous incorporation of tautomers, ionization species, and different binding modes into ligand-based and receptor-based 3D-QSAR methods

Speciation of drug candidates and receptors caused by ionization, tautomerism, and/or covalent hydration complicates ligandand receptor-based predictions of binding affinities by 3-dimensional structure-activity relationships (3D-QSAR). The speciation problem is exacerbated by tendency of tautomers to bind in multiple conformations or orientations (modes) in the same binding site. New forms of the 3D-QSAR correlation equations, capable of capturing this complexity, can be developed using the time hierarchy of all steps that lie behind the monitored biological process - binding, enzyme inhibition or receptor activity. In most cases, reversible interconversions of individual ligand and receptor species can be treated as quickly established equilibria because they are finished in a small fraction of the exposure time that is used to determine biological effects. The speciation equilibria are satisfactorily approximated by invariant fractions of individual ligand and receptor species for buffered experimental or in vivo conditions. For such situations, the observed drug-receptor association constant of a ligand is expressed as the sum of products, for each ligand and receptor species pair, of the association microconstant and the fractions of involved species. For multiple binding modes, each microconstant is expressed as the sum of microconstants of individual modes. This master equation leads to new 3D-QSAR correlation equations integrating the results of all molecular simulations or calculations, which are run for each ligand-receptor species pair separately. The multispecies, multimode 3D-QSAR approach is illustrated by a ligand-based correlation of transthyretin binding of thyroxine analogs and by a receptor-based correlation of inhibition of MK2 by benzothiophenes and pyrrolopyrimidines.

Gas-phase hydration thermochemistry of sodiated and potassiated nucleic acid bases

Hydration reactions of sodiated and potassiated nucleic acid bases (uracil, thymine, cytosine, and adenine) produced by electrospray have been studied in a gas phase using the pulsed ion-beam high-pressure mass spectrometer. The thermochemical properties, ΔH(o)(n), ΔS(o)(n), and ΔG(o)(n), for the hydrated systems were obtained from hydration equilibrium measurement. The structural aspects of the hydrated complexes are discussed in conjunction with available literature data. The correlation between water binding energies in the hydrated complexes and the corresponding metal ion affinities of nucleobases suggests that a significant (if not dominant) amount of the canonical structure of cytosine undergoes tautomerization during electrospray ionization, and the thermochemical values for cationized cytosine probably correspond to a mixture of tautomeric complexes.

Structures and physical properties of gaseous metal cationized biological ions

Metal chelation can alter the activity of free biomolecules by modifying their structures or stabilizing higher energy tautomers. In recent years, mass spectrometric techniques have been used to investigate the effects of metal complexation with proteins, nucleobases and nucleotides, where small conformational changes can have significant physiological consequences. In particular, infrared multiple photon dissociation spectroscopy has emerged as an important tool for determining the structure and reactivity of gas-phase ions. Unlike other mass spectrometric approaches, this method is able to directly resolve structural isomers using characteristic vibrational signatures. Other activation and dissociation methods, such as blackbody infrared radiative dissociation or collision-induced dissociation can also reveal information about the thermochemistry and dissociative pathways of these biological ions. This information can then be used to provide information about the structures of the ionic complexes under study. In this article, we review the use of gas-phase techniques in characterizing metal-bound biomolecules. Particular attention will be given to our own contributions, which detail the ability of metal cations to disrupt nucleobase pairs, direct the self-assembly of nucleobase clusters and stabilize non-canonical isomers of amino acids.

Theoretical studies of the in-solution isomeric protonation of non-aromatic six-member rings with two nitrogens

For exploring the preferred site for hydrogen bond formation, theoretical calculations have been performed for a number of six-member, nonaromatic rings allowing for alternative protonation on the ring nitrogens. Gas-phase protonation studies for test molecules indicate that the B3LYP/aug-cc-pvtz and QCISD(T)(CBS) calculations approach the experimental values within about 1 kcal/mol with considerable improvement for relative enthalpies and free energies. Relative free energies calculated at the IEF-PCM/B3LYP/aug-cc-pvtz level predict favorable protonation on the tertiary rather than on the secondary nitrogen both in aqueous solution and in a dichloromethane solvent for saturated rings. Protonation on a nitrogen atom next to a C═C bond is disfavored due to a large increase in internal energy. Monte Carlo simulations considering a counterion and Ewald summation for the long-range electrostatic effects for a 0.1 molar model system predict ΔG(solv)/MC values generally less negative than from the IEF-PCM calculations. These results make the protonation on the tertiary nitrogen even more favored. The solute-solvent pair-energy distribution depends sensitively on the applied model. In conclusion, the freely moving anion has been considered as the most relevant model with overall neutrality for the system and applying the least restrictions.

Structures of alkali metal ion-adenine complexes and hydrated complexes by IRMPD spectroscopy and electronic structure calculations

Complexes between adenine and the alkali metal ions Li(+), Na(+), K(+), and Cs(+) have been investigated by infrared multiple photon dissociation (IRMPD) spectroscopy between 2800 and 3900 cm(-1), as have some singly hydrated complexes. The IRMPD spectra clearly show the N-H stretching and the NH(2) symmetric and asymmetric stretching vibrations of adenine; and for the solvated ions, the O-H stretching vibrations are observed. These experimental spectra were compared with those for a variety of possible structures, including both A9 (A9 refers to the tautomer where hydrogen is on the nitrogen in position 9 of adenine, see Scheme 1) and A7 adenine tautomers, computed using B3-LYP/6-31+G(d,p). By comparing the experimental and the simulated spectra it is possible to rule out various structures and to further assign structures to the species probed in these experiments. Single-point calculations on the B3-LYP/6-31+G(d,p) geometries have been performed at MP2/6-311++G(2d, p) to obtain good estimates of the relative thermochemistries for the different structures. In all cases the computed IR spectrum for the lowest energy structure is consistent with the experimental IRMPD spectrum, but in some cases structural assignment cannot be confirmed based solely upon comparison with the experimental spectra so computed thermochemistries can be used to rule out high-energy structures. On the basis of the IRMPD spectra and the energy calculations, all adenine-M(+) and adenine-M(+)-H(2)O are concluded to be composed of the A7 tautomer of adenine, which is bound to the cations in a bidentate fashion through N3 and N9 (see Scheme 1 for numbering convention). For the hydrated ions water binds directly to the metal ion through oxygen, as would be expected since the metal contains most positive charge density. For the hydrated lithium cation-bound adenine dimer, the water molecule is concluded to be hydrogen bonded to a free basic site of one of the adenine monomers, which is also bound to the lithium cation. Experimental and theoretical results on adenine-Li(+)-H(2)O suggest that the electrosprayed adenine-Li(+) resembles the lowest-energy solution phase ion rather than the lowest-energy gas-phase ion, which is the imine form.