Spectroscopy ofd-wave superconductors using DNA as a probing tip

Andreev reflection in DNA molecules terminated by ad-wave superconductor is investigated for demonstrating advantages in using DNA as the probe for the spectroscopy of the superconductor. DNA molecules are incorporated in the simulations using a two-leg ladder model with a simplification as constructed by homopolymers. The increase of the Andreev reflection probability at zero bias originating from the midgap surface states ofd-wave superconductors appears even when the DNA molecule is coupled strongly. The zero-bias peak is enhanced by orders of magnitude when the coupling is weakened. The one-dimensional transport in DNA strands gives rise to the remarkable sensitivity in the spectroscopy, where the changes of the reflection probability caused by the midgap states are also in orders of magnitude when the voltage bias is varied and when the orientation of thed-wave symmetry is inclined with respect to the superconductor surface. The quantum interference of the transport between the two strands in DNA modifies the zero-bias increase. A narrow dip occurs at zero bias with a plateau sandwiched by two peaks in the immediate vicinity of the zero bias. The characteristics of these transmission resonances are dependent on the parameters describing the model molecules, and so the width of the resonance peaks, for instance, enables us to evaluate the strength of the inter-strand coupling.

Quantum Effects Explain the Twist Angle in the Helical Structure of DNA

Why are DNA bases stacked in a double helix structure? We combined three theoretical approaches to demonstrate how one core concept derived from quantum mechanics (Pauli repulsion) annihilates the contribution of dispersion to the π-π stacking. The helical architecture is governed by a combination of exchange and electrostatic forces, a result that is interpreted from both a computational and a biological perspective.

The Influence of 2'-Deoxyguanosine Lesions on the Electronic Properties of OXOG:::C Base Pairs in Ds-DNA: A Comparative Analysis of Theoretical Studies

DNA is continuously exposed to a variety of harmful factors, which, on the one hand, can force undesirable processes such as ageing, carcinogenesis and mutagenesis, while on the other hand, can accelerate evolutionary changes. Of all the canonical nucleosides, 2'-deoxyguanosine (dG) exhibits the lowest ionization potential, making it particularly prone to the one-electron oxidizing process. The most abundant type of nucleobase damage is constituted by 7,8-dihydro-8-oxo-2'-deoxyguanosine (OXOdG), with an oxidation potential that is 0.56 V lower than that of canonical dG. All this has led to OXOdG, as an isolated lesion, being perceived as a sink for radical cations in the genome. In this paper, a comparative analysis of the electronic properties of an OXOGC base pair within the context of a clustered DNA lesion (CDL) has been conducted. It is based on previous DFT studies that were carried out at the M06-2x/6-31++G** level of theory in non-equilibrated and equilibrated condensed phases. The results of the comparative analysis presented here reveal the following: (A) The ionization potentials of OXOG4C2 were largely unaffected by a second lesion. (B) The positive charge and spin were found predominantly on the OXOG4C2 moiety. (C) The electron-hole transfers A3T3→G4C2 and G4C2←A5T1 were found in the Marcus inverted region and were resistant to the presence of a second DNA lesion in close proximity. It can therefore be reasonably postulated that OXOGC becomes the sink for a radical cation migrating through the double helix, irrespective of the presence of other 2'-deoxyguanosine lesions in the CDL structure.

Repurposing the antihistamine drug bilastine as an anti-cancer metallic drug entity: synthesis and single-crystal X-ray structure analysis of metal-based bilastine and phen [Co(II), Cu(II) and Zn(II)] tailored anticancer chemotherapeutic agents against resistant cancer cells

Bilastine (BLA), 2-(4-(2-(4-(1-(2-ethoxyethyl)-1H-benzo[d]imidazole-2-yl)-piperidin-1-yl)-ethyl)-phenyl)-2-methylpropanoic acid, is an active antihistamine drug. With the idea of repurposing drugs from the existing pool of 'active' pharmaceutical ingredients, the therapeutic potency of bilastine as an anticancer agent was investigated via the tailored synthesis of a metal-based anticancer drug formulation of the type [BLA(phen)2M(II)]+·X-, where M = Co, Cu, and Zn and X- = NO3 and ClO4. The synthesized metal-based chemotherapeutics derived from the bilastine drug that acts as a ligand were thoroughly characterized using spectroscopic techniques, namely, UV-vis, FT-IR, and EPR (in the case of 1 and 2); 1H-NMR and 13C-NMR (in the case of 3); ESI-MS and single-crystal X-ray diffraction studies. Comprehensive biological studies (DNA binding, cleavage, and cytotoxic activity) using various biophysical and gel electrophoretic methods were carried out to validate their potential as anticancer agents. The cytotoxic activity of 'therapeutically promising' copper(II)-based drug candidate 2 was evaluated against MCF-7, MBA-MD-231, HeLa, HepG2, and Mia-PaCa-2 cancer cells via an SRB assay, and the results demonstrated 2 as a potent anticancer agent at low nanomolar concentrations against all tested cancer cells, preferably with a much superior anticancer efficacy against human pancreatic cancer cells.

Effects of Chalcogen Atoms on Excited-State Double-Proton Transfer Behavior for 3,6-bis(4,5-Dihydroxyoxazo-2-yl)benzene-1,2-diol Derivatives: A Computational Investigation

The impact of the chalcogen atomic electronegativity (O, S, and Se atoms) of new organic molecules on excited-state dynamical reactions is self-evident. Inspired by this kind of distinguished photochemical characteristic, in this work, we performed a computational investigation of chalcogen-substituted 3,6-bis(4,5-dihydroxyoxazo-2-yl)benzene-1,2-diol (BDYBD) derivatives (i.e., BDYBD-O, BDYBD-S, and BDYBD-Se). In this paper, we pay close attention to characteristic BDYBD derivatives that contain intramolecular double hydrogen bonds (O1-H2···N3 and O4-H5···N6). The main goal of this study was to explore how changes in atomic electronegativity affect the way hydrogen bonds interact and how excited molecules affect transfer protons. We go into further detail in the main text of the paper. By fixing our attention to geometrical variations and infrared (IR) vibrational spectra between the S0 and S1 states, exploring hydrogen bonding behaviors using the core-valence bifurcation (CVB) index, and simulating hydrogen bonding energy (EHB) via the atom in molecule (AIM) method, we clarified the photo-induced strengthened dual hydrogen bonding interactions that facilitate the excited-state dual-proton transfer (ESDPT) behavior of BDYBD derivatives. The reorganization of charge stemming from photoexcitation further verifies the tendencies of ESDPT reactions. We relied on constructing potential energy surfaces (PESs) by adopting a restrictive optimization approach, and herein, we finally clarify the gradual ESDPT mechanism of BDYBD derivatives. Particularly, we confirm that the variation in chalcogen atomic electronegativity has a regulatory effect on the ESDPT behavior of BDYBD derivatives; that is, the lower the atomic electronegativity, the more favorable it is for the gradual ESDPT reaction.

Why we are made of proteins and nucleic acids: Structural biology views on extraterrestrial life

Is it a miracle that life exists on the Earth, or is it a common phenomenon in the universe? If extraterrestrial organisms exist, what are they like? To answer these questions, we must understand what kinds of molecules could evolve into life, or in other words, what properties are generally required to perform biological functions and store genetic information. This review summarizes recent findings on simple ancestral proteins, outlines the basic knowledge in textbooks, and discusses the generally required properties for biological molecules from structural biology viewpoints (e.g., restriction of shapes, and types of intra- and intermolecular interactions), leading to the conclusion that proteins and nucleic acids are at least one of the simplest (and perhaps very common) forms of catalytic and genetic biopolymers in the universe. This review article is an extended version of the Japanese article, On the Origin of Life: Coevolution between RNA and Peptide, published in SEIBUTSU BUTSURI Vol. 61, p. 232-235 (2021).

Complementary base pair interactions between different rare tautomers of the second-generation artificial genetic alphabets

The functionality of a semisynthetic DNA in the biological environment will depend on the base pair nature of its complementary base pairs. To understand this, base pair interactions between complementary bases of recently proposed eight second-generation artificial nucleobases are studied herein by considering their rare tautomeric conformations and a dispersion-corrected density functional theoretic method. It is found that the binding energies of two hydrogen-bonded complementary base pairs are more negative than those of the three hydrogen-bonded base pairs. However, as the former base pairs are endothermic, the semisynthetic duplex DNA would involve the latter base pairs.

Interaction of an anticancer benzopyrane derivative with DNA: Biophysical, biochemical, and molecular modeling studies

BACKGROUND

SIMR1281 is a potent anticancer lead candidate with multi- target activity against several proteins; however, its mechanism of action at the molecular level is not fully understood. Revealing the mechanism and the origin of multitarget activity is important for the rational identification and optimization of multitarget drugs.

METHODS

We have used a variety of biophysical (circular dichroism, isothermal titration calorimetry, viscosity, and UV DNA melting), biochemical (topoisomerase I & II assays) and computational (molecular docking and MD simulations) methods to study the interaction of SIMR1281 with duplex DNA structures.

RESULTS

The biophysical results revealed that SIMR1281 binds to dsDNA via an intercalation-binding mode with an average binding constant of 3.1 × 10⁶ M^-1. This binding mode was confirmed by the topoisomerases' inhibition assays and molecular modeling simulations, which showed the intercalation of the benzopyrane moiety between DNA base pairs, while the remaining moieties (thiazole and phenyl rings) sit in the minor groove and interact with the flanking base pairs adjacent to the intercalation site.

CONCLUSIONS

The DNA binding characteristics of SIMR1281, which can disrupt/inhibit DNA function as confirmed by the topoisomerases' inhibition assays, indicate that the observed multi-target activity might originate from ligand intervention at nucleic acids level rather than due to direct interactions with multiple biological targets at the protein level.

GENERAL SIGNIFICANCE

The findings of this study could be helpful to guide future optimization of benzopyrane-based ligands for therapeutic purposes.

Activating charge-transfer state formation in strongly-coupled dimers using DNA scaffolds

Strongly-coupled multichromophoric assemblies orchestrate the absorption, transport, and conversion of photonic energy in natural and synthetic systems. Programming these functionalities involves the production of materials in which chromophore placement is precisely controlled. DNA nanomaterials have emerged as a programmable scaffold that introduces the control necessary to select desired excitonic properties. While the ability to control photophysical processes, such as energy transport, has been established, similar control over photochemical processes, such as interchromophore charge transfer, has not been demonstrated in DNA. In particular, charge transfer requires the presence of close-range interchromophoric interactions, which have a particularly steep distance dependence, but are required for eventual energy conversion. Here, we report a DNA-chromophore platform in which long-range excitonic couplings and short-range charge-transfer couplings can be tailored. Using combinatorial screening, we discovered chromophore geometries that enhance or suppress photochemistry. We combined spectroscopic and computational results to establish the presence of symmetry-breaking charge transfer in DNA-scaffolded squaraines, which had not been previously achieved in these chromophores. Our results demonstrate that the geometric control introduced through the DNA can access otherwise inaccessible processes and program the evolution of excitonic states of molecular chromophores, opening up opportunities for designer photoactive materials for light harvesting and computation.

RNABPDB: Molecular Modeling of RNA Structure-From Base Pair Analysis in Crystals to Structure Prediction

The stable three-dimensional structure of RNA is known to play several important biochemical roles, from post-transcriptional gene regulation to enzymatic action. These structures contain double-helical regions, which often have different types of non-canonical base pairs in addition to Watson-Crick base pairs. Hence, it is important to study their structures from experimentally obtained or even predicted ones, to understand their role, or to develop a drug against the potential targets. Molecular Modeling of RNA double helices containing non-canonical base pairs is a difficult process, particularly due to the unavailability of structural features of non-Watson-Crick base pairs. Here we show a composite web-server with an associated database that allows one to generate the structure of RNA double helix containing non-canonical base pairs using consensus parameters obtained from the database. The database classification is followed by an evaluation of the central tendency of the structural parameters as well as a quantitative estimation of interaction strengths. These parameters are used to construct three-dimensional structures of double helices composed of Watson-Crick and/or non-canonical base pairs. Our benchmark study to regenerate double-helical fragments of many experimentally derived RNA structures indicate very high accuracy. This composite server is expected to be highly useful in understanding functions of various pre-miRNA by modeling structures of the molecules and estimating binding efficiency. The database can be accessed from http://hdrnas.saha.ac.in/rnabpdb .

Solvent effects on the NMR shieldings of stacked DNA base pairs

Stacking effects are among the most important effects in DNA. We have recently studied their influence in fragments of DNA through the analysis of NMR magnetic shieldings, firstly in vacuo. As a continuation of this line of research we show here the influence of solvent effects on the shieldings through the application of both explicit and implicit models. We found that the explicit solvent model is more appropriate for consideration due to the results matching better in general with experiments, as well as providing clear knowledge of the electronic origin of the value of the shieldings. Our study is grounded on a recently developed theoretical model of our own, by which we are able to learn about the magnetic effects of given fragments of DNA molecules on selected base pairs. We use the shieldings of the atoms of a central base pair (guanine-cytosine) of a selected fragment of DNA molecules as descriptors of physical effects, like π-stacking and solvent effects. They can be taken separately and altogether. The effect of π-stacking is introduced through the addition of some pairs above and below of the central base pair, and now, the solvent effect is considered including a network of water molecules that consist of two solvation layers, which were fixed in the calculations performed in all fragments. We show that the solvent effects enhance the stacking effects on the magnetic shieldings of atoms that belong to the external N-H bonds. The net effect is of deshielding on both atoms. There is also a deshielding effect on the carbon atoms that belong to CO bonds, for which the oxygen atom has an explicit hydrogen bond (HB) with a solvent water molecule. Solvent effects are found to be no higher than a few percent of the total value of the shieldings (between 1% and 5%) for most atoms, although there are few for which such an effect can be higher. There is one nitrogen atom, the acceptor of the HB between guanine and cytosine, that is more highly shielded (around 15 ppm or 10%) when the explicit solvent is considered. In a similar manner, the most external nitrogen atom of cytosine and the hydrogen atom that is bonded to it are highly deshielded (around 10 ppm for nitrogen and around 3 ppm for hydrogen).

Affinity and Correlation in DNA

A statistical analysis of important DNA sequences and related proteins has been performed to study the relationships between monomers, and some general considerations about these macromolecules can be provided from the results. First, the most important relationship between sites in all the DNA sequences examined is that between two consecutive base pairs. This is an indication of an energetic stabilization due to the stacking interaction of these couples of base pairs. Secondly, the difference between human chromosome sequences and their coding parts is relevant both in the relationships between sites and in some specific compositional rules, such as the second Chargaff rule. Third, the evidence of the relationship in two successive triplets of DNA coding sequences generates a relationship between two successive amino acids in the proteins. This is obviously impossible if all the relationships between the sites are statistical evidence and do not involve causes; therefore, in this article, due to stacking interactions and this relationship in coding sequences, we will divide the concept of the relationship between sites into two concepts: affinity and correlation, the first with physical causes and the second without. Finally, from the statistical analyses carried out, it will emerge that the human genome is uniform, with the only significant exception being the Y chromosome.

Determination of two-photon absorption in nucleobase analogues: a QR-DFT perspective

With the prevalence of fluorescence spectroscopy in biological systems, and the benefits of two-photon absorption techniques, presented here is an assessment of the two-photon accessibility of modern fluorescent nucleobase analogues utilising quadratic response DFT. Due to the complex environment experienced by these nucleobases, the two-photon spectra of each analogue has been assessed in the presence of both [Formula: see text]-stacked and hydrogen-bonding interactions involving the canonical nucleobases. Findings suggest that the [Formula: see text]-stacking environment provides a more significant effect on the spectra of the analogues studies than a hydrogen-bonding environment; analogue structures presenting high two-photon cross-section values for one or more states coincide with polycyclic extensions to preserved canonical base structure, as observed in the qA family of analogues, while analogue structures more closely resembling the structure of the base in question present a much more muted spectra in comparison. Results from this investigation have also allowed for the derivation of a number of design rules for the development of potential, two-photon specific, analogues for future use in both imaging and potential photochemical activation.

B-DNA Structure and Stability: The Role of Nucleotide Composition and Order

We have quantum chemically analyzed the influence of nucleotide composition and sequence (that is, order) on the stability of double-stranded B-DNA triplets in aqueous solution. To this end, we have investigated the structure and bonding of all 32 possible DNA duplexes with Watson-Crick base pairing, using dispersion-corrected DFT at the BLYP-D3(BJ)/TZ2P level and COSMO for simulating aqueous solvation. We find enhanced stabilities for duplexes possessing a higher GC base pair content. Our activation strain analyses unexpectedly identify the loss of stacking interactions within individual strands as a destabilizing factor in the duplex formation, in addition to the better-known effects of partial desolvation. Furthermore, we show that the sequence-dependent differences in the interaction energy for duplexes of the same overall base pair composition result from the so-called "diagonal interactions" or "cross terms". Whether cross terms are stabilizing or destabilizing depends on the nature of the electrostatic interaction between polar functional groups in the pertinent nucleobases.

Orientation Dependence of DNA Blunt-End Stacking Studied by Free-Energy Simulations

DNA blunt ends can associate mediated by stacking interactions between the terminal base pairs that form blunt ends. The blunt end association plays a role in DNA repair and recombination processes and can also be of importance for the design of DNA-based nano-materials. Its function depends on the sequence and on the geometric arrangement that leads to stable interaction. For a stacked state, the relative orientation (twisting) of the base pairs is important. Molecular dynamics and advanced sampling simulations were used to calculate free energy change associated with twist changes of the stacked blunt-end base pairs. The calculations reproduce blunt stacking arrangements found in crystal structures of DNA oligonucleotides as free energy minima. To elucidate the physical origin of the stabilization of certain angular arrangements, the interactions between backbone atoms in the blunt-end stack were switched off in additional free energy calculations. It allows us to decipher the contributions to stacking stabilization due to the nucleobases and the backbone and to analyze the sequence dependence of the angular stacking preferences. Good qualitative agreement was also found for the comparison with quantum mechanical calculations. The results may help in the design of novel DNA-based materials.

Unveiling the complex pattern of intermolecular interactions responsible for the stability of the DNA duplex

Herein, we provide new insights into the intermolecular interactions responsible for the intrinsic stability of the duplex structure of a large portion of human B-DNA by using advanced quantum mechanical methods. Our results indicate that (i) the effect of non-neighboring bases on the inter-strand interaction is negligibly small, (ii) London dispersion effects are essential for the stability of the duplex structure, (iii) the largest contribution to the stability of the duplex structure is the Watson-Crick base pairing - consistent with previous computational investigations, (iv) the effect of stacking between adjacent bases is relatively small but still essential for the duplex structure stability and (v) there are no cooperativity effects between intra-strand stacking and inter-strand base pairing interactions. These results are consistent with atomic force microscope measurements and provide the first theoretical validation of nearest neighbor approaches for predicting thermodynamic data of arbitrary DNA sequences.

Pyrrole and Pyridine in the Water Environment-Effect of Discrete and Continuum Solvation Models

Properties of pyrrole and pyridine molecules upon different hydrations were investigated through density functional theory. Complexes of studied molecules with a cluster of 50 water molecules were considered, and the polarizable continuum model of solvation (PCM) was also taken into account. For comparative purposes, all mentioned calculations were repeated for single pyrrole and pyridine molecules and their complexes with one water molecule. Aromaticities of solvated pyrrole and pyridine rings were studied using several geometric- and electronic-based aromaticity criteria. Special attention was paid to studying the properties of formed hydrogen bonds between pyrrole or pyridine and either a single water molecule or several water molecules of the cluster. Overall, a comprehensive description of two very important heterocyclic compounds, that is, pyrrole and pyridine, in both discrete and continuum water solutions, is extensively presented.

The influence of oxoG on the electronic properties of ds-DNA. Damage versus mismatch: A theoretical approach

The seed of life is concealed in the base sequence in DNA. This macromolecule is continuously exposed to harmful factors which can cause it damage. The stability of genetic information depends on the protein efficiency of repair systems. Glycosylases are the scouts which recognize and remove damaged bases. Their efficiency depends on how rapidly they recognize DNA lesions. One theory states that charge transfer is involved in protein cross talking through ds-DNA. For these reasons a comparative analysis of ds-oligo containing a mismatched base pair dA:::dG and a damaged dA::dG^OXO is proposed. Additionally, the electronic properties of the short ds-oligo in the context of non-equilibrated and equilibrated solvent modes were taken into theoretical consideration. All energetic calculations were performed at the M062x/6-31++G** level of theory, while for geometry optimized ONIOM methodology was used. The lowest adiabatic ionization potential was assigned for DNA containing a dA:dG^OXO pair. Moreover, the adiabatic electron affinity was assigned at the same level for the mismatched and lesioned ds-oligo. Surprisingly, in the non-equilibrated mode, a significantly higher vertical electro affinity was found for lesioned DNA. The higher VEA in a non-equilibrated solvent state supported faster recognition in the A:G^OXO base pair than A:G by MutY glycosylases under electron transfer mechanism.

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.

Cross-Talk between Overlap Interactions in Biomolecules: A Case Study of the β-Turn Motif

Noncovalent interactions play a pivotal role in regulating protein conformation, stability and dynamics. Among the quantum mechanical (QM) overlap-based noncovalent interactions, n→π* is the best understood with studies ranging from small molecules to β-turns of model proteins such as GB1. However, these investigations do not explore the interplay between multiple overlap interactions in contributing to local structure and stability. In this work, we identify and characterize all noncovalent overlap interactions in the β-turn, an important secondary structural element that facilitates the folding of a polypeptide chain. Invoking a QM framework of natural bond orbitals, we demonstrate the role of several additional interactions such as n→σ* and π→π* that are energetically comparable to or larger than n→π*. We find that these interactions are sensitive to changes in the side chain of the residues in the β-turn of GB1, suggesting that the n→π* may not be the only component in dictating β-turn conformation and stability. Furthermore, a database search of n→σ* and π→π* in the PDB reveals that they are prevalent in most proteins and have significant interaction energies (∼1 kcal/mol). This indicates that all overlap interactions must be taken into account to obtain a comprehensive picture of their contributions to protein structure and energetics. Lastly, based on the extent of QM overlaps and interaction energies, we propose geometric criteria using which these additional interactions can be efficiently tracked in broad database searches.

DNA Ring-Opening Polymerization Driven by Base Stacking

Supramolecular polymers, relying on reversible intermolecular interactions, promise a wide range of applications, including optoelectronic materials, self-healing materials, and biomedical delivery materials. Among potential molecular candidates, DNA strands act as an excellent platform. DNA has a well-established secondary structure (double helix), and its intermolecular interactions can be readily thermodynamically engineered and kinetically controlled. Extensive studies have demonstrated that various DNA motifs can polymerize/assemble into large polymers with different topology, geometry, and dimensionalities. Most of the reported polymerization is driven by hybridization of DNA strands. Herein, we report a novel system of DNA supramolecular polymerization that is driven by DNA base stacking. The polymerization has been confirmed by native polyacrylamide gel electrophoresis (PAGE) and atomic force microscopy (AFM). We believe that this work will expand the toolbox for DNA supramolecular polymerization and would, with further development, increase further control of DNA supramolecular polymerization.

Understanding alkali metal cation affinities of multi-layer guanine quadruplex DNA

To gain better understanding of the stabilizing interactions between metal ions and DNA quadruplexes, dispersion-corrected density functional theory (DFT-D) based calculations were performed on double-, triple- and four-layer guanine tetrads interacting with alkali metal cations. All computations were performed in aqueous solution that mimics artificial supramolecular conditions where guanine bases assemble into stacked quartets as well as biological environments in which telomeric quadruplexes are formed. To facilitate the computations on these significant larger systems, optimization of the DFT description was performed first by evaluating the performance of partial reduced basis sets. Analysis of the stabilizing interactions between alkali cations and the DNA bases in double and triple-layer guanine quadruplex DNA reproduced the experimental affinity trend of the order Li⁺< Rb⁺ < Na⁺ < K⁺. The desolvation and the size of alkali metal cations are thought to be responsible for the order of affinity. Nevertheless, for the alkali metal cation species individually, the magnitude of the bond energy stays equal for binding as first, second or third cation in double, triple and four-layer guanine quadruplexes, respectively. This is the result of an interplay between a decreasingly stabilizing interaction energy and increasingly stabilizing solvation effects, along the consecutive binding events. This diminished interaction energy is the result of destabilizing electrostatic repulsion between the hosted alkali metal cations. This work emphasizes the stabilizing effect of aqueous solvent on large highly charged biomolecules.

Stabilizing solvent effects and electrostatic repulsion are responsible for the constant alkali metal cation affinity in multi-layer guanine quadruplexes.

Enzymatic Formation of an Artificial Base Pair Using a Modified Purine Nucleoside Triphosphate

The expansion of the genetic alphabet with additional, unnatural base pairs (UBPs) is an important and long-standing goal in synthetic biology. Nucleotides acting as ligands for the coordination of metal cations have advanced as promising candidates for such an expansion of the genetic alphabet. However, the inclusion of artificial metal base pairs in nucleic acids mainly relies on solid-phase synthesis approaches, and very little is known about polymerase-mediated synthesis. Herein, we report the selective and high yielding enzymatic construction of a silver-mediated base pair (dIm^C-Ag^I-dPur^P) as well as a two-step protocol for the synthesis of DNA duplexes containing such an artificial metal base pair. Guided by DFT calculations, we also shed light into the mechanism of formation of this artificial base pair as well as into the structural and energetic preferences. The enzymatic synthesis of the dIm^C-Ag^I-dPur^P artificial metal base pair provides valuable insights for the design of future, more potent systems aiming at expanding the genetic alphabet.

Impact of 5-formylcytosine on the melting kinetics of DNA by 1H NMR chemical exchange

5-Formylcytosine (5fC) is a chemically edited, naturally occurring nucleobase which appears in the context of modified DNA strands. The understanding of the impact of 5fC on dsDNA physical properties is to date limited. In this work, we applied temperature-dependent 1H Chemical Exchange Saturation Transfer (CEST) NMR experiments to non-invasively and site-specifically measure the thermodynamic and kinetic influence of formylated cytosine nucleobase on the melting process involving dsDNA. Incorporation of 5fC within symmetrically positioned CpG sites destabilizes the whole dsDNA structure-as witnessed from the ∼2°C decrease in the melting temperature and 5-10 kJ mol-1 decrease in ΔG°-and affects the kinetic rates of association and dissociation. We observed an up to ∼5-fold enhancement of the dsDNA dissociation and an up to ∼3-fold reduction in ssDNA association rate constants, over multiple temperatures and for several proton reporters. Eyring and van't Hoff analysis proved that the destabilization is not localized, instead all base-pairs are affected and the transition states resembles the single-stranded conformation. These results advance our knowledge about the role of 5fC as a semi-permanent epigenetic modification and assist in the understanding of its interactions with reader proteins.

Mutual Relations between Substituent Effect, Hydrogen Bonding, and Aromaticity in Adenine-Uracil and Adenine-Adenine Base Pairs

The electronic structure of substituted molecules is governed, to a significant extent, by the substituent effect (SE). In this paper, SEs in selected nucleic acid base pairs (Watson-Crick, Hoogsteen, adenine-adenine) are analyzed, with special emphasis on their influence on intramolecular interactions, aromaticity, and base pair hydrogen bonding. Quantum chemistry methods-DFT calculations, the natural bond orbital (NBO) approach, the Harmonic Oscillator Model of Aromaticity (HOMA) index, the charge of the substituent active region (cSAR) model, and the quantum theory of atoms in molecules (QTAIM)-are used to compare SEs acting on adenine moiety and H-bonds from various substitution positions. Comparisons of classical SEs in adenine with those observed in para- and meta-substituted benzenes allow for the better interpretation of the obtained results. Hydrogen bond stability and its other characteristics (e.g., covalency) can be significantly changed as a result of the SE, and its consequences are dependent on the substitution position. These changes allow us to investigate specific relations between H-bond parameters, leading to conclusions concerning the nature of hydrogen bonding in adenine dimers-e.g., H-bonds formed by five-membered ring nitrogen acceptor atoms have an inferior, less pronounced covalent nature as compared to those formed by six-membered ring nitrogen. The energies of individual H-bonds (obtained by the NBO method) are analyzed and compared to those predicted by the Espinosa-Molins-Lecomte (EML) model. Moreover, both SE and H-bonds can significantly affect the aromaticity of adenine rings; long-distance SEs on π-electron delocalization are also documented.

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.

Revisiting the Potential Energy Surface of the Stacked Cytosine Dimer: FNO-CCSD(T) Interaction Energies, SAPT Decompositions, and Benchmarking

Nucleobase stacking interactions are crucial for the stability of nucleic acids. This study investigates base stacking energies of the cytosine homodimer in different configurations, including intermolecular separation plots, detailed twist dependence, and displaced structures. Highly accurate ab initio quantum chemical single point energies using an energy function based on MP2 complete basis set extrapolation ([6 → 7]ZaPa-NR) and a CCSD(T)/cc-pVTZ-F12 high-level correction are presented as new reference data, providing the most accurate stacking energies of nucleobase dimers currently available. Accurate SAPT2+(3)δMP2 energy decomposition is used to obtain detailed insights into the nature of base stacking interactions at varying vertical distances and twist values. The ab initio symmetry adapted perturbation theory (SAPT) energy decomposition suggests that the base stacking originates from an intricate interplay between dispersion attraction, short-range exchange-repulsion, and Coulomb interaction. The interpretation of the SAPT data is a complex issue as key energy terms vary substantially in the region of optimal (low energy) base stacking geometries. Thus, attempts to highlight one leading stabilizing SAPT base stacking term may be misleading and the outcome strongly depends on the used geometries within the range of geometries sampled in nucleic acids upon thermal fluctuations. Modern dispersion-corrected density functional theory (among them DSD-BLYP-D3, ωB97M-V, and ωB97M-D3BJ) is benchmarked and often reaches up to spectroscopic accuracy (below 1 kJ/mol). The classical AMBER force field is benchmarked with multiple different sets of point-charges (e.g. HF, DFT, and MP2-based) and is found to produce reasonable agreement with the benchmark data.

Molecular Mechanisms of DNA Replication and Repair Machinery: Insights from Microscopic Simulations

Reproduction, the hallmark of biological activity, requires making an accurate copy of the genetic material to allow the progeny to inherit parental traits. In all living cells, the process of DNA replication is carried out by a concerted action of multiple protein species forming a loose protein-nucleic acid complex, the replisome. Proofreading and error correction generally accompany replication but also occur independently, safeguarding genetic information through all phases of the cell cycle. Advances in biochemical characterization of intracellular processes, proteomics and the advent of single-molecule biophysics have brought about a treasure trove of information awaiting to be assembled into an accurate mechanistic model of the DNA replication process. In this review, we describe recent efforts to model elements of DNA replication and repair processes using computer simulations, an approach that has gained immense popularity in many areas of molecular biophysics but has yet to become mainstream in the DNA metabolism community. We highlight the use of diverse computational methods to address specific problems of the fields and discuss unexplored possibilities that lie ahead for the computational approaches in these areas.

Intermolecular magnetic interactions in stacked DNA base pairs

The influence of pi-stacking on the magnetic properties of atoms that belong to adenine-thymine and guanine-cytosine pairs in sequences of three and five layers of DNA base pairs was analysed. As probes we used NMR spectroscopic parameters, which are among the most useful tools to learn about the transmission of magnetic interactions in molecules. Four DFT functionals were employed: B3LYP, BHANDLYP, KT2 and KT3, together with the SOPPA method. Besides, given that the number of non-hydrogen atoms of the supramolecular systems studied here is larger than 50 we applied a locally dense basis set scheme. Our results show that the piling up of a few Watson-Crick base pairs above and below a given pair modifies its NMR spectroscopic parameters by an amount that may be measurable and the percentage of variation does not depend on dispersion. We found that magnetic shieldings are more sensitive than J-couplings, and also that some atoms are more sensitive than others. Stacking affects the shielding of non-hydrogen atoms like nitrogens, that are donors in hydrogen bonds, HBs, and the carbons bonded to them. The amount of variation of these shieldings was found to be from 2% to 5% when the pairs are considered first as isolated, and then, placed in the middle of a sequence of three layers of base pairs. Such a variation becomes vanishingly small when the sequence contains more than three layers, showing that the stacking effect on NMR spectroscopic parameters has a local nature. We have also found a pattern for shieldings. First, equivalent atoms of similar monomers (thymine and adenine, or guanine and cytosine) have similar values of absolute shieldings in isolated pairs, and the amount of variation from isolated pairs to aggregates of a few pairs is also similar, meaning that equivalent atoms are affected in a similar manner by pi-stacking. Second, the hydrogen atoms which belong to hydrogen bonds are more sensitive to the piling up than the non-hydrogen atoms.

Unveiled electric profiles within hydrogen bonds suggest DNA base pairs with similar bond strengths

Electrical forces are the background of all the interactions occurring in biochemical systems. From here and by using a combination of ab-initio and ad-hoc models, we introduce the first description of electric field profiles with intrabond resolution to support a characterization of single bond forces attending to its electrical origin. This fundamental issue has eluded a physical description so far. Our method is applied to describe hydrogen bonds (HB) in DNA base pairs. Numerical results reveal that base pairs in DNA could be equivalent considering HB strength contributions, which challenges previous interpretations of thermodynamic properties of DNA based on the assumption that Adenine/Thymine pairs are weaker than Guanine/Cytosine pairs due to the sole difference in the number of HB. Thus, our methodology provides solid foundations to support the development of extended models intended to go deeper into the molecular mechanisms of DNA functioning.

DNA Polymerase λ Active Site Favors a Mutagenic Mispair between the Enol Form of Deoxyguanosine Triphosphate Substrate and the Keto Form of Thymidine Template: A Free Energy Perturbation Study

Human DNA polymerase λ is an intermediate fidelity member of the X family, which plays a role in DNA repair. Recent X-ray diffraction structures of a ternary complex of a loop-deletion mutant of polymerase λ, a deoxyguanosine triphosphate analogue, and a gapped DNA show that guanine and thymine form a mutagenic mispair with an unexpected Watson-Crick-like geometry rather than a wobble geometry. Hence, there is an intriguing possibility that either thymine in the DNA or guanine in the deoxyguanosine triphosphate analogue may spend a substantial fraction of time in a deprotonated or enol form (both are minor species in aqueous solution) in the active site of the polymerase λ mutant. The experiments do not determine particular forms of the nucleobases that contribute to this mutagenic mispair. Thus, we investigate the thermodynamics of formation of various mispairs between guanine and thymine in the ternary complex at a neutral pH using classical molecular dynamics simulations and the free energy perturbation method. Our free energy calculations, as well as a comparison of the experimental and computed structures of mispairs, indicate that the Watson-Crick-like mispair between the enol tautomer of guanine and the keto tautomer of thymine is dominant. The wobble mispair between the keto forms of guanine and thymine and the Watson-Crick-like mispair between the keto tautomer of guanine and the enol tautomer of thymine are less prevalent, and mispairs that involve deprotonated guanine or thymine are thermodynamically unlikely. These findings are consistent with the experiment and relevant for understanding mechanisms of spontaneous mutagenesis.

Effects of oxygenation on the intercalation of 1,10-phenanthroline-5,6/4,7-dione between DNA base pairs: a computational study

The effects of oxygen in positions 4,7 and 5,6 of phenanthroline have been studied computationally when this ligand intercalates between DNA base pairs. Our results indicate that solvation energy could be the driving force of the process and thus, it can be also related with the cytotoxicity of the drug.

Stacking interactions involving non-Watson–Crick basepairs: dispersion corrected density functional theory studies

Stacking interactions between a non Watson–Crick G:A S:HT basepair and C:G basepair is predicted in terms of roll, twist and slide basepair step parameters using DFT-D augmented with coarse-grain energy penalty for sugar–phosphate backbone.

How the electron-deficient cavity of heterocalixarenes recognizes anions: insights from computation

The nature of bridging heteroatoms in a heterocalixarene structure has a crucial role in anion recognition.

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.

The role of alkali metal cations in the stabilization of guanine quadruplexes: why K+ is the best

The desolvation and size of monovalent alkali metal ions are of equal importance for the cation affinity of guanine quadruplexes.