Designing a New Class of Bases for Nucleic Acid Quadruplexes and Quadruplex-Active Ligands

Variable temperature FTIR spectra of polycrystalline purine nucleobases and estimating strengths of individual hydrogen bonds

In the first part of this work, we report the FTIR spectra of pure NH and isotopically substituted ND (10-15% D and 80-90% D) polycrystalline hypoxanthine, xanthine, adenine and guanine recorded in the 400-4000 cm^-1 range, as a function of temperature (10-300 K). We provide assignments of the stretching and out-of-plane bending amine (NH₂) and imine (NH) bands to the distinct H-bonds present in the crystal, based on the temperature sensitivity and isotopic exchange behavior. Empirical correlations between spectral and thermodynamic or structural parameters enabled us to estimate the energies and lengths of H-bonds in the studied nucleobase crystals and to correlate them with literature data. The empirical H-bonding energies are compared with H-bonding and stacking energies computed for hypoxanthine. In the second part, strategies for using the empirical correlations together with information extracted from quantum mechanical data (in particular from the Bader's quantum theory of atoms in molecules, QTAIM) for the evaluation of hydrogen bonding properties are discussed, and their advantages and drawbacks pointed out. The justification for a cooperative use of quantum-mechanical calculations with empirical spectra-energy correlations is discussed.

Anatomy of Base Pairing in DNA by Interacting Quantum Atoms

The formation of purine and pyrimidine base pairs (BPs), which contributes to shaping of the canonical and noncanonical 3D structures of nucleic acids, is one the most investigated phenomena in chemistry and life sciences. In this contribution, the anatomy of the bond energy (BDE) of the base-pairing interaction in 39 different arrangements found experimentally or predicted for DNA structures containing the four common nucleobases (A, C, G, T) in their neutral or protonated forms is described in light of the theory of interacting quantum atoms within the context of the quantum theory of atoms in molecules. The interplay of individual energy components involved in the three stages of the bond formation process (structural deformation, electron-density promotion, and intermolecular interaction) is studied. We recognized that for the neutral BPs, variations in the kinetic and electrostatic contributions to the BDE are rather negligible, leaving the exchange-correlation energy as the main stabilizing component. It is shown that the contribution of the exchange-correlation term can be recovered by including atoms that are formally assumed to be hydrogen bonded (primary interaction). In contrast, to recover the electrostatic component of interaction, one must consider both the primary and secondary (formally nonbonded atoms) interatomic interactions. The results of our study were employed to design new types of BPs with altered bonding anatomy. We demonstrate that improving the electrostatic characteristics of the BPs does not necessarily result in greater interaction energies if weak secondary hydrogen bonding is destroyed. However, the main tuning factor for systems with conserved interacting faces (primary interactions) is the electrostatic component of the interaction energy resulting from the secondary atom-atom electrostatics.

The evaluation of 5-amino- and 5-hydroxyuracil derivatives as potential quadruplex-forming agents

5-Substituted uracils (NH₂ or OH groups in position 5) have been examined theoretically and experimentally as potential building blocks in quadruplex structures. Our high level Density Functional Theory (DFT) calculations showed that the tetramer formation and stacking energies for 5-substituted uracils are similar to the energies of purine-based xanthine (X) or guanine (G) structures. As tetrads of 5-substituted uracils cover almost exactly the same area as purine tetrads, mixed tetrads or quadruplex structures based on X or G and 5-substituted uracil motifs are possible. According to the calculations, 5-hydroxyuracil-based structures are the best candidates for experimental implementation which was corroborated by the existence of higher complexes in the mass spectra of 1-benzyl-5-hydroxyuracil. These pyrimidine-based molecules can be used as efficient building blocks in different applications including aptamers, bio-sensors or - taking into account the larger cavity in the central region of 5-hydroxyuracil structures - as an artificial ion channel.

In What Ways Do Synthetic Nucleotides and Natural Base Lesions Alter the Structural Stability of G-Quadruplex Nucleic Acids?

Synthetic analogs of natural nucleotides have long been utilized for structural studies of canonical and noncanonical nucleic acids, including the extensively investigated polymorphic G-quadruplexes (GQs). Dependence on the sequence and nucleotide modifications of the folding landscape of GQs has been reviewed by several recent studies. Here, an overview is compiled on the thermodynamic stability of the modified GQ folds and on how the stereochemical preferences of more than 70 synthetic and natural derivatives of nucleotides substituting for natural ones determine the stability as well as the conformation. Groups of nucleotide analogs only stabilize or only destabilize the GQ, while the majority of analogs alter the GQ stability in both ways. This depends on the preferred syn or anti N-glycosidic linkage of the modified building blocks, the position of substitution, and the folding architecture of the native GQ. Natural base lesions and epigenetic modifications of GQs explored so far also stabilize or destabilize the GQ assemblies. Learning the effect of synthetic nucleotide analogs on the stability of GQs can assist in engineering a required stable GQ topology, and exploring the in vitro action of the single and clustered natural base damage on GQ architectures may provide indications for the cellular events.

Effect of Monovalent Ion Parameters on Molecular Dynamics Simulations of G-Quadruplexes

G-quadruplexes (GQs) are key noncanonical DNA and RNA architectures stabilized by desolvated monovalent cations present in their central channels. We analyze extended atomistic molecular dynamics simulations (∼580 μs in total) of GQs with 11 monovalent cation parametrizations, assessing GQ overall structural stability, dynamics of internal cations, and distortions of the G-tetrad geometries. Majority of simulations were executed with the SPC/E water model; however, test simulations with TIP3P and OPC water models are also reported. The identity and parametrization of ions strongly affect behavior of a tetramolecular d[GGG]₄ GQ, which is unstable with several ion parametrizations. The remaining studied RNA and DNA GQs are structurally stable, though the G-tetrad geometries are always deformed by bifurcated H-bonding in a parametrization-specific manner. Thus, basic 10-μs-scale simulations of fully folded GQs can be safely done with a number of cation parametrizations. However, there are parametrization-specific differences and basic force-field errors affecting the quantitative description of ion-tetrad interactions, which may significantly affect studies of the ion-binding processes and description of the GQ folding landscape. Our d[GGG]₄ simulations indirectly suggest that such studies will also be sensitive to the water models. During exchanges with bulk water, the Na⁺ ions move inside the GQs in a concerted manner, while larger relocations of the K⁺ ions are typically separated. We suggest that the Joung-Cheatham SPC/E K⁺ parameters represent a safe choice in simulation studies of GQs, though variation of ion parameters can be used for specific simulation goals.

Channeling through Two Stacked Guanine Quartets of One and Two Alkali Cations in the Li+, Na+, K+, and Rb+ Series. Assessment of the Accuracy of the SIBFA Anisotropic Polarizable Molecular Mechanics Potential

Stacking of guanine quartets (GQs) can trigger the formation of DNA or RNA quadruple helices, which play numerous biochemical roles. The GQs are stabilized by alkali cations, mainly K⁺ and Na⁺, which can reside in, or channel through, the central axis of the GQ stems. Further, ion conduction through GQ wires can be leveraged for nanochemistry applications. G-quadruplex systems have been extensively studied by classical molecular dynamics (MD) simulations using pair-additive force fields or by quantum-chemical (QC) calculations. However, the non-polarizable force fields are very approximate, while QC calculations lack the necessary sampling. Thus, ultimate description of GQ systems would require long-enough simulations using advanced polarizable molecular mechanics (MM). However, to perform such calculations, it is first mandatory to evaluate the method's accuracy using benchmark QC. We report such an evaluation for SIBFA polarizable MM, bearing on the channeling (movement) of an alkali cation (Li⁺, Na⁺, K⁺, or Rb⁺) along the axis of two stacked G quartets interacting with either one or two ions. The QC energy profiles display markedly different features depending upon the cation but can be retrieved in the majority of cases by the SIBFA profiles. An appropriate balance of first-order (electrostatic and short-range repulsion) and second-order (polarization, charge-transfer, and dispersion) contributions within ΔE is mandatory. With two cations in the channel, the relative weights of the second-order contributions increase steadily upon increasing the ion size. In the G8 complexes with two K⁺ or two Rb⁺ cations, the sum of polarization and charge-transfer exceeds the first order terms for all ion positions.

Weak Supramolecular Interactions Governing Parallel and Antiparallel DNA Quadruplexes: Insights from Large-Scale Quantum Mechanics Analysis of Experimentally Derived Models

The topology and energetics of guanine (G) quadruplexes is governed by supramolecular interactions within their strands. In this work, an extensive quantum mechanical (QM) study has been performed to analyze supramolecular interactions that shape the stems of (4+0) parallel (P) and (2+2) antiparallel (AP) quadruplex systems. The large-scale (≈400 atoms) models of P and AP were constructed from high-quality experimental structures. The results provide evidence that each of the P and AP structures is shaped by a distinct network of supramolecular interactions. Analysis of electron topological characteristics of hydrogen bonds in P and AP systems indicates that the P model benefits from stronger intratetrad hydrogen bonding. For intertetrad stacking interactions, both noncovalent interaction plot and energy decomposition analysis approaches suggest that the stem of the P quadruplex benefits more from stacking than that of the AP stem; the difference in energetic stabilization for the two topologies is about 10 %. Stronger hydrogen-bonding and stacking interactions in the stem of the P quadruplex, relative to those in the AP system, can be an important indicator to explain the experimental observations that guanine-rich oligonucleotides tend to form all-parallel stems with an all-anti orientation of nucleobases. However, in addition to intrinsic stabilization, partial desolvation effects, which affect the energetics and dynamics of the G-quadruplex folding process, call for further investigations.

B-DNA model systems in non-terran bio-solvents: implications for structure, stability and replication

We have computationally analyzed a comprehensive series of Watson–Crick and mismatched B-DNA base pairs, in the gas phase and in several solvents, including toluene, chloroform, ammonia, methanol and water, using dispersion-corrected density functional theory and implicit solvation.

Modified Guanines as Constituents of Smart Ligands for Nucleic Acid Quadruplexes

Repetitive guanine-rich nucleic acid sequences play a crucial role in maintaining genome stability and the cell life cycle and represent potential targets for regulatory drugs. Recently, it has been demonstrated that guanine-based ligands with a porphyrin core can be used as markers of G-quadruplex assemblies in cell tissues. Herein, model systems of guanine-based ligands are explored by DFT methods. The energies of formation of modified guanine tetrads and those of modified tetrads stacked on the top of natural guanine tetrads have been calculated. The interaction energy has been decomposed into contributions from hydrogen bonding, stacking, and ion coordination and a twist-rise potential energy scan has been performed to find the individual local minima. Energy decomposition analysis reveals the impact of various substituents (F, Cl, Br, I, Me, NMe2 ) on individual energy terms. In addition, cooperative reinforcement in forming the modified and stacked tetrads, as well as the frontier orbitals participating in the hydrogen-bonding framework involving the HOMO-LUMO gap between the occupied σHOMO on the proton-accepting C=O and =N- groups and unoccupied σLUMO on the N-H groups, has been studied. The investigated systems are demonstrated to have a potential in ligand development, mainly due to stacking enhancement compared with natural guanine, which is used as a reference.

Nucleotides containing variously modified sugars: energetics, structure, and mechanical properties

The influence of various sugar residue modifications on intrinsic energetic, conformational, and mechanical properties of 2'-deoxyribonucleotide-5'-monophosphates (dNs) was comprehensively investigated using modern quantum chemical approaches. In total, fourteen sugar modifications, including double bonds and heteroatoms (S and N) inside the sugar ring, as well as fluorination in various positions, were analyzed. Among hundreds of possible conformational states of dNs, only two - AI and BI, corresponding to the most biologically significant forms of a double-helical DNA, were considered for each dN. It was established that the most of the studied modifications tend to strongly stabilize either AI or BI conformation of dNs both in the gas phase and in aqueous solution (modelled by implicit solvent models). Therefore, some of these modifications can be used as a tool for reducing structural polymorphism of nucleic acids in solution as well as for designing oligonucleotides with specific structural features. The evaluation of relaxed force constants (RFC) for glycosidic bonds suggests that the majority of the studied modifications of the sugar residue yield increased strengths of glycosidic bonds in dNs, and can therefore be used for designing modified nucleic acids with an increased resistance to abasic lesions. The most significant reinforcement of the glycosidic bond occurs in dNs containing the CF2 group instead of the O4' oxygen and the fluorine atom at the 2'-α-position. The calculation of the RFC and vibrational root-mean-square (VRMS) deviations for conformational degrees of freedom revealed a strong dependence between mechanical properties of dNs and their energetic characteristics. In particular, electronic energies of AI and BI conformers of dNs calculated in vacuo are closely connected with the values of relaxed force constants (RFC) for the δ angle: the higher RFC(δ) values correspond to more energetically favorable conformers.

Substituting CF2 for O4' in Components of Nucleic Acids: Towards Systems with Reduced Propensity to Form Abasic Lesions

Intrinsic structural features and energetics of nucleotides containing variously fluorinated sugars as potential building blocks of DNA duplexes and quadruplexes are explored systematically using the modern methods of density functional theory (DFT) and quantum chemical topology (QCT). Our results suggest that fluorination at the 2'-β or 2'-α,β positions somewhat stabilizes in vacuo the AI relative to the BI conformations. In contrast, substitution of the CF2 group for the O4' atom (O4'-CF2 modification) leads to a preference of the BI relative to AI DNA-like conformers. All the studied modifications result in a noticeable increase in the stability of the glycosidic bond [estimated by the relaxed force constants (RFC) approach], with particularly encouraging results for the O4'-CF2 derivative. Consequently, the O4'-CF2 modified systems are suggested and explored as promising scaffolds for the development of duplex and quadruplex structures with reduced propensity to form abasic lesions and to undergo DNA damage.