Vibrational spectroscopy simulation of solvation effects on a G-quadruplex

It is commonly believed that solvation effects on the vibrational properties of a solute are easily accounted for by simple rules of thumbs, that is, solvating a polar molecule in a polar medium has the only effect of red shifting all its spectroscopical features and, similarly, solvating a polar molecule in a nonpolar medium has the opposite effect. In this work, we use theoretical vibrational spectroscopy at quasi-classical and quantum approximate semiclassical level to gain atomistic insights about solvent-solute interactions for 2'-deoxyguanosine and the G-quadruplex. We employ the quasi-classical trajectory method to include full anharmonicity into our calculated spectra, and then introduce quantum nuclear effects by means of divide-and-conquer semiclassical spectroscopy calculations. Solvation is treated explicitly leading to a good reproducibility of the available experimental data and reliable predictions when an experimental reference is missing.Communicated by Ramaswamy H. Sarma.

Mass Spectrometry of Nucleic Acid Noncovalent Complexes

Nucleic acids have been among the first targets for antitumor drugs and antibiotics. With the unveiling of new biological roles in regulation of gene expression, specific DNA and RNA structures have become very attractive targets, especially when the corresponding proteins are undruggable. Biophysical assays to assess target structure as well as ligand binding stoichiometry, affinity, specificity, and binding modes are part of the drug development process. Mass spectrometry offers unique advantages as a biophysical method owing to its ability to distinguish each stoichiometry present in a mixture. In addition, advanced mass spectrometry approaches (reactive probing, fragmentation techniques, ion mobility spectrometry, ion spectroscopy) provide more detailed information on the complexes. Here, we review the fundamentals of mass spectrometry and all its particularities when studying noncovalent nucleic acid structures, and then review what has been learned thanks to mass spectrometry on nucleic acid structures, self-assemblies (e.g., duplexes or G-quadruplexes), and their complexes with ligands.

Elucidating Noncovalent Reaction Mechanisms: G-Quartet as an Intermediate in G-Quadruplex Assembly

In analogy to covalent reactions, the understanding of noncovalent association pathways is fundamental to influence and control any supramolecular process. Following an approach that is reminiscent of covalent methodologies, we study here, for the first time, the mechanism of G-quadruplex formation in organic solvents. Our results support a reaction pathway in which the cation shifts the equilibrium towards a G-quartet transient intermediate, which then acts as a template in the formation of the G-quadruplex product.

Hydrogen bonding influences collision-induced dissociation of Na+ -bound guanine tetrads

Na+ -bound guanine (G)-tetrads possess square planar structures formed solely by noncovalent interactions including multiple hydrogen bonds. Unlike G-tetrads facilitated by other alkali metal ions, an intriguing behavior in collision-induced dissociation (CID) has been observed in Na+ -bound G-tetrads, which features a preferential, simultaneous loss of two G ligands in the low energy regime. To understand this unique behavior, we investigated the CID of Na+ -bound G-tetrads with mixed ligands of G and 9-methylguanine (9mG), [Na·Gm ·9mGn ]+ (m + n = 4), and [Li·9mG4 ]+ for comparison. In the CID experiments, the simultaneous losses of two ligands were by far more pronounced than the loss of a single ligand for all five Na+ -bound G-tetrads. However, it appeared that the CID of [Li·9mG4 ]+ prefers to lose single ligands sequentially. An analysis of the fragment abundances suggested that the generation of Na+ -bound dimeric fragments might have occurred with two adjacent ligands. This theoretical study predicted for [Li·9mG4 ]+ that the loss of a single ligand is more energetically favorable than the production of neutral hydrogen-bonded fragments by 35.5 kJ/mol (ΔG). This contradicts our previous calculations for [Na·9mG4 ]+ that a neutral loss of hydrogen-bonded dimers provides the lowest energy product state of Na+ -bound dimeric fragments, which is lower than that of Na+ -bound trimeric fragments by 15.6 kJ/mol. From the results, this comparative study suggests that the pronounced generation of Na+ -bound dimeric fragments in CID of the G-tetrads is likely promoted by the dissociation pathway associated with neutral loss of hydrogen-bonded dimers. It thus demonstrates that multiple hydrogen bonding participating in formation of Na+ -bound G-tetrads may also strongly influence the fate of dissociating complexes of G-tetrads.

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.

Guidelines for the assembly of hydrogen-bonded macrocycles

This article highlights selected examples on the synthesis of hydrogen-bonded macrocycles from ditopic molecules and analyze the main factors, often interrelated, that influence the equilibrium between ring and chain species.

The K2(9-ethylguanine)122+ quadruplex is more stable to unimolecular dissociation than the K(9-ethylguanine)8+ quadruplex in the gas phase: a BIRD, energy resolved SORI-CID, IRMPD spectroscopic, and computational study

A combination of experimental trapped-ion mass spectrometric studies and computational chemistry has been used to assess the intrinsic properties of the potassiated 9-ethylguanine (9eG) self-assembled quadruplex, K₂(9eG)₁₂²⁺, in the gas phase.

Quantum Mechanical Investigation of the G-Quadruplex Systems of Human Telomere

The three G-quadruplexes involved in the human telomere have been studied with an accurate quantum mechanical approach, and the possibility of reducing them to a simpler model has been tested. The similarities and the differences of these three systems are shown and discussed. Each system has been analyzed through different properties and compared to the others. In particular, we have considered: (1) the shape of the cavity and the atomic charges around it; (2) the electric field in and out of the cavity; (3) the stabilization energy due to the stacking of G-tetrads, to the H-bonds and to the ion interactions; and, finally, (4) to study the mechanism of the process of the ion inclusion in the cavity, the curves of potential energy due to the movement of the Na+ and K+ ions toward the cavity. The results suggest that a detailed study is essential in order to obtain the quantitative properties of these complex systems, but also that some qualitative behaviors can be schematized. Our study makes it clear that the entry of an ion in the cavity of these systems is a complex process, where it is possible to find stable structures with the ion out and in the cavity. Moreover, it is possible that more than one diabatic state is involved in this process.

Quantum mechanical investigation of G-quartet systems of DNA

Minima of the electric field and positions of K⁺ and Na⁺ (zero of the x-coordinate is the center of the cavity).

The intrinsic stabilities and structures of alkali metal cationized guanine quadruplexes

The structures and stabilities of self-assembled guanine quadruplexes, M(9eG)₈⁺ (M = Na, K, Rb, Cs; 9eG = 9-ethylguanine), have been studied in the gas phase by blackbody infrared radiative dissociation kinetics to determine the effect the metal cations have on the decomposition energies and reactions of the quadruplex.

Hydrogen-Bonded Macrocyclic Supramolecular Systems in Solution and on Surfaces

Cyclization into closed assemblies is the most recurrent approach to realize the noncovalent synthesis of discrete, well-defined nanostructures. This review article particularly focuses on the noncovalent synthesis of monocyclic hydrogen-bonded systems that are self-assembled from a single molecule with two binding-sites. Taking advantage of intramolecular binding events, which are favored with respect to intermolecular binding in solution, can afford quantitative amounts of a given supramolecular species under thermodynamic control. The size of the assembly depends on geometric issues such as the monomer structure and the directionality of the binding interaction, whereas the fidelity achieved relies largely on structural preorganization, low degrees of conformational flexibility, and templating effects. Here, we discuss several examples described in the literature in which cycles of different sizes, from dimers to hexamers, are studied by diverse solution or surface characterization techniques.

Guanine tetrads: an IRMPD spectroscopy, energy resolved SORI-CID, and computational study of M(9-ethylguanine)4+ (M = Li, Na, K, Rb, Cs) in the gas phase

The electrostatic potential map of the neutral guanine tetrad showing the negative hole in the centre.