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

Thermodynamic Reversal and Structural Correlation of 24-Crown-8/Protonated Tryptophan and 24-Crown 8/Protonated Serine Noncovalent Complexes in the Gas Phase vs in Solution: Quantum Chemical Analysis

We investigate the structures of 24-crown-8/H+/l-tryptophan (CR/TrpH+) and 24-crown-8/H+/l-serine (CR/SerH+) noncovalent host-guest complex both in the gas phase and in an aqueous solution by quantum chemical methods. The Gibbs free energies of the complex in the two phases are calculated to determine the thermodynamically most favorable conformer in each phase. Our predictions indicate that both the carboxyl and the ammonium in CR/TrpH+ and the ammonium in the CR/SerH+ complexes in the lowest Gibbs free energy configurations form hydrogen bonds (H-bonds) with the CR host in the gas phase, while the conformer with the "naked" (devoid of H-bond with the CR host) -CO2H (and/or -OH) is much less favorable (Gibbs free energy higher by >3.6 kcal/mol). In the solution phase, however, a "thermodynamic reversal" occurs, making the higher Gibbs free energy gas-phase CR/TrpH+ and CR/SerH+ conformers thermodynamically more favorable under the influence of solvent molecules. Consequently, the global minimum Gibbs free energy structure in solution is structurally correlated with the thermodynamically much less gas-phase conformer. Discussions are provided concerning the possibility of elucidating host-guest-solvent interactions in solution from the gas-phase host-guest configurations in molecular detail.

Study of protonated dimers of cytosine, cytidine, and deoxycytidine using survival yield method and quantum mechanics calculations

RATIONALE

Cytosine and its conjugates are prone to form protonated, triply-bonded dimers. Therefore, the nucleic-acid cytosine-rich sequence forms the four-stranded noncanonical secondary structure known as the intercalated motif (i-motif). This process has resulted in studies on cytosine protonated dimers. This communication focuses on the protonated dimers of cytosine and its nucleoside using the survival yield (SY) method and quantum mechanics calculations.

METHODS

To obtain the precursor ion fragmentation curve, the plot of SY against Ecomδ , the product ion spectra of the protonated dimers were obtained using a Waters/Micromass Q-TOF Premier mass spectrometer. Quantum mechanics calculations were performed using GAUSSIAN 16, and full geometry optimizations and energy calculations were performed within the density functional theory framework at B3LYP/6-31G(d,p).

RESULTS

The precursor ion fragmentation curve allowed the rating of the gas-phase stabilities of the analyzed protonated dimers. Substitution of sugar moiety at N1 cytosine atom decreased the gas-phase stabilities of the protonated dimers. The deoxycytidine dimer was found to be more stable than the cytidine dimer and cytidine-deoxycytidine dimer. Quantum chemical calculations indicated that cytosine aminohydroxy tautomer may be involved in the formation of protonated cytosine-cytosine nucleoside dimers but not in the formation of cytosine dimers.

CONCLUSIONS

The results obtained for nucleoside dimers indicated that the SY method may reflect the i-motif stabilities observed under physiological conditions. Therefore, the analysis of other protonated dimers of variously substituted cytosine-cytosine nucleoside using the SY method may be important to study the effect of cytosine substitution on the i-motif stabilities. Cytosine tautomer containing C2-OH… N(2H)-C4 moiety may be involved in the formation of protonated cytosine-cytosine nucleoside dimers.

Distinguishing gas phase lactose and lactulose complexed with sodiated L-arginine by IRMPD spectroscopy and DFT calculations

We present the origin of the observed differentiation of lactose and lactulose achieved by complexation with sodiated L-arginine (ArgNa+). We find that the infrared multiphoton dissociation (IRMPD) bands in 3600-3650 and >3650 cm-1 regimes for gas phase lactose and lactulose, respectively, vanish when forming host-guest complexes with ArgNa+. We interpret these differences in the IRMPD spectra by scrutinizing the interactions between the functional groups (guanidium, -CO2-Na+) in ArgNa+ and -OHs in lactose/lactulose. Our calculated structures and infrared spectra of lactose/ArgNa+ and lactulose/ArgNa+ host-guest pairs indicate that the functional groups interact with the low- and high-frequency -OH stretch modes of lactose and lactulose, respectively, in the 3600-3720 cm-1 window.

Unveiling host-guest-solvent interactions in solution by identifying highly unstable host-guest configurations in thermal non-equilibrium gas phase

We propose a novel scheme of examining the host-guest-solvent interactions in solution from their gas phase structures. By adopting the permethylated β-cyclodextrin (perm β-CD)-protonated L-Lysine non-covalent complex as a prototypical system, we present the infrared multiple photon dissociation (IRMPD) spectrum of the gas phase complex produced by electrospray ionization technique. In order to elucidate the structure of perm β-CD)/LysH⁺ complex in the gas phase, we carry out quantum chemical calculations to assign the two strong peaks at 3,340 and 3,560 cm^-1 in the IRMPD spectrum, finding that the carboxyl forms loose hydrogen bonding with the perm β-CD, whereas the ammonium group of L-Lysine is away from the perm β-CD unit. By simulating the structures of perm β-CD/H⁺/L-Lysine complex in solution using the supramolecule/continuum model, we find that the extremely unstable gas phase structure corresponds to the most stable conformer in solution.

Hydrogen Bonding in Natural and Unnatural Base Pairs-A Local Vibrational Mode Study

In this work hydrogen bonding in a diverse set of 36 unnatural and the three natural Watson Crick base pairs adenine (A)-thymine (T), adenine (A)-uracil (U) and guanine (G)-cytosine (C) was assessed utilizing local vibrational force constants derived from the local mode analysis, originally introduced by Konkoli and Cremer as a unique bond strength measure based on vibrational spectroscopy. The local mode analysis was complemented by the topological analysis of the electronic density and the natural bond orbital analysis. The most interesting findings of our study are that (i) hydrogen bonding in Watson Crick base pairs is not exceptionally strong and (ii) the N-H⋯N is the most favorable hydrogen bond in both unnatural and natural base pairs while O-H⋯N/O bonds are the less favorable in unnatural base pairs and not found at all in natural base pairs. In addition, the important role of non-classical C-H⋯N/O bonds for the stabilization of base pairs was revealed, especially the role of C-H⋯O bonds in Watson Crick base pairs. Hydrogen bonding in Watson Crick base pairs modeled in the DNA via a QM/MM approach showed that the DNA environment increases the strength of the central N-H⋯N bond and the C-H⋯O bonds, and at the same time decreases the strength of the N-H⋯O bond. However, the general trends observed in the gas phase calculations remain unchanged. The new methodology presented and tested in this work provides the bioengineering community with an efficient design tool to assess and predict the type and strength of hydrogen bonding in artificial base pairs.

A vibrational spectroscopic and computational study of gaseous protonated and alkali metal cationized G-C base pairs

The structures and properties of metal cationized complexes of 9-ethylguanine (9eG) and 1-methylcytosine (1mC), (9eG:1mC)M⁺, where M⁺ = Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺ as well as the protonated complex, (9eG:1mC)H⁺, have been studied using a combination of IRMPD spectroscopy and computational methods. For (9eG:1mC)H⁺, the dominant structure is a Hoogsteen type complex with the proton covalently bound to N3 of 1mC despite this being the third best protonation site of the two bases; based on proton affinities N7 of 9eG should be protonated. However, this structural oddity can be explained considering both the number of hydrogen bonds that can be formed when N3 of 1mC is protonated as well as the strong ion-induced dipole interaction that exists between an N3 protonated 1mC and 9eG due to the higher polarizability of 9eG. The anomalous dissociation of (9eG:1mC)H⁺, forming much more (1mC)H⁺ than would be predicted based on the computed thermochemistry, can be explained as being due to the structural oddity of the protonation site and that the barrier to proton transfer from N3 of 1mC to N7 of 9eG grows dramatically as the base pair begins to dissociate. For the (9eG:1mC)M⁺; M = Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺ complexes, single unique structures could not be assigned. However, the experimental spectra were consistent with the computed spectra. For (9eG:1mC)Li⁺, the lowest energy structure is one in which Li⁺ is bound to O6 of 9eG and both O2 and N3 of 1mC; there is also an interbase hydrogen bond from the amine of 1mC to N7 of 9eG. For Na⁺, K⁺, and Rb⁺, similar binding of the metal cation to 1mC is calculated but, unlike Li⁺, the lowest energy structure is one in which the metal cation is bound to N7 of 9eG; there is also an interbase hydrogen bond between the amine of 1mC and the carbonyl of 9eG. The lowest energy structure for the Cs complex is the Watson-Crick type base pairing with Cs⁺ binding only to 9eG through O6 and N7 and with three hydrogen bonds between 9eG and 1mC. It also interesting to note that the Watson-Crick base pairing structure gets lower in Gibbs energy relative to the lowest energy complexes as the metal gets larger. This indicates that the smaller, more densely charged cations have a greater propensity to interfere with Watson-Crick base pairing than do the larger, less densely charged metal cations.

An IRMPD spectroscopic and computational study of protonated guanine-containing mismatched base pairs in the gas phase

Infrared multiple photon dissociation spectroscopy has been used to probe the structures of the three protonated base-pair mismatches containing 9-ethylguanine (9eG) in the gas phase. Some of these protonated base-pairs have been identified in RNA.