True stabilization energies for the optimal planar hydrogen-bonded and stacked structures of guanine...cytosine, adenine...thymine, and their 9- and 1-methyl derivatives: complete basis set calculations at the MP2 and CCSD(T) levels and comparison with experiment

Ab initio Modeling of Hydrogen Bonding of Remdesivir and Adenosine with Uridine

Remdesivir (RDV) emerged as an effective drug against the SARS-CoV-2 virus pandemic. One of the crucial steps in the mechanism of action of RDV is its incorporation into the growing RNA strand. RDV, an adenosine analogue, forms Watson-Crick (WC) type hydrogen bonds with uridine in the complementary strand and the strength of this interaction will control efficacy of RDV. While there is a plethora of structural and energetic information available about WC H-bonds in natural base pairs, the interaction of RDV with uridine has not been studied yet at the atomic level. In this article, we aim to bridge this gap, to understand RDV and its hydrogen bonding interactions, by employing density functional theory (DFT) at the M06-2X/cc-pVDZ level. The interaction energy, QTAIM analysis, NBO and SAPT2 are performed for RDV, adenosine, and their complex with uridine to gain insights into the nature of hydrogen bonding. The computations show that RDV has similar geometry, energetic, molecular orbitals, and aromaticity as adenosine, suggesting that RDV is an effective adenosine analogue. The important geometrical parameters, such as bond distances and red-shift in the stretching vibrational modes of adenosine, RDV and uridine identify two WC-type H-bonds. The relative strength of these two H-bonds is computed using QTAIM parameters and the computed hydrogen bond energy. Finally, the SAPT2 study is performed at the minima and at non-equilibrium base pair distances to understand the dominant intermolecular physical force. This study, based on a thorough analysis of a variety of computations, suggests that both adenosine and RDV have similar structure, energetic, and hydrogen bonding behaviour.

Tryptophan-PNA gc Conjugates Self-Assemble to Form Fibers

Peptide nucleic acids and their conjugates to peptides can self-assemble and generate complex architectures. In this work, we explored the self-assembly of PNA dimers conjugated to the dipeptide WW. Our studies suggest that the indole ring of tryptophan promotes aggregation of the conjugates. The onset of fluorescence is observed upon self-assembly. The structure of self-assembled WWgc is concentration-dependent, being spherical at low concentrations and fibrous at high concentrations. As suggested by molecular modeling studies, fibers are stabilized by stacking interactions between tryptophans and Watson-Crick hydrogen bonds between nucleobases.

High-Order Quantum-Mechanical Analysis of Hydrogen Bonding in Hachimoji and Natural DNA Base Pairs

High-order quantum chemistry is applied to hydrogen-bonded natural DNA nucleobase pairs [adenine:thymine (A:T) and guanine:cytosine (G:C)] and non-natural Hachimoji nucleobase pairs [isoguanine:1-methylcytosine (B:S) and 2-aminoimidazo[1,2a][1,3,5]triazin-4(1H)-one:6-amino-5-nitropyridin-2-one (P:Z)] to see how the intermolecular interaction energies and their energetic components (electrostatics, exchange-repulsion, induction/polarization, and London dispersion interactions) vary among the base pairs. We examined the Hoogsteen (HG) geometries in addition to the traditional Watson-Crick (WC) geometries. Coupled-cluster theory through perturbative triples [CCSD(T)] extrapolated to the complete basis set (CBS) limit and high-order symmetry-adapted perturbation theory (SAPT) at the SAPT2+(3)(CCD)δMP2/aug-cc-pVTZ level are used to estimate highly accurate noncovalent interaction energies. Electrostatic interactions are the most attractive component of the interaction energies, but the sum of induction/polarization and London dispersion is nearly as large, for all base pairs and geometries considered. Interestingly, the non-natural Hachimoji base pairs interact more strongly than the corresponding natural base pairs, by -21.8 (B:S) and -0.3 (P:Z) kcal mol^-1 in the WC geometries, according to CCSD(T)/CBS. This is consistent with the H-bond distances being generally shorter in the non-natural base pairs. The natural base pairs are energetically more stabilized in their Hoogsteen geometries than in their WC geometries. The Hoogsteen geometry makes the A:T base pair slightly more stable, by -0.8 kcal mol^-1, and it greatly stabilizes the G:C⁺ base pair, by -15.3 kcal mol^-1. The G:C⁺ stabilization is mainly due to the fact that C has typically added a proton when found in Hoogsteen geometries. By contrast, Hoogsteen geometries are substantially less favorable than WC geometries for non-natural Hachimoji base pairs, by 17.3 (B:S) and 13.8 (P:Z) kcal mol^-1.

Photoionization Dynamics and Proton Transfer within the Adenine-Thymine Nucleobase Pair

Studying the stability of hydrogen-bonded nucleobase pairs, at the heart of the genetic code, is of utmost importance for an in-depth understanding of basic mechanisms of life and biomolecular evolution. We present here a VUV single photon ionization dynamic study of the nucleobase pair adenine-thymine (AT), revealing its ionization and dissociative ionization thresholds via double imaging electron/ion coincidence spectroscopy. The experimental data, consisting of cluster mass-resolved threshold photoelectron spectra and photon energy-dependent ion kinetic energy release distributions, allow the unambiguous distinction of the dissociation of AT into protonated adenine AH⁺ and a dehydrogenated thymine radical T(-H) from dissociative ionization processes of other nucleobase clusters. Comparison to high-level ab initio calculations indicates that our experimental observations can be explained by a single hydrogen-bonded conformer present in our molecular beam and allows the estimation of an upper limit of the barrier of the proton transfer in the ionized AT pair.

Aggregation of nucleobases and metabolites: Adenine-theobromine trimers

The selection of cytosine, guanine, thymine, and adenine as components of the information biopolymers was a complex process influenced by several factors. Among them, the intermolecular interactions may have played a determinant role. Thus, a deep understanding of the intermolecular interactions between nucleobases and other prebiotic molecules may help understand the first instants of chemical evolution. Following this hypothesis, we present here a combined spectroscopic and computational study of theobromine₂-adenine and thebromine-adenine₂ trimers. While adenine is a nucleobase, theobromine was probably part of the prebiotic chemistry. The trimers were formed in jets and probed by a combination of UV and IR spectroscopic techniques. The spectra were interpreted in light of the predictions obtained using density-functional methods. The results suggest the existence of a subtle balance between formation of hydrogen bonds and π-π interactions. Thus, while theobromine₂-adenine tends to form complex in stacked structures, theobromine-adenine₂ prefers formation of planar structures, maximizing the interaction by hydrogen bonds. The small energy difference between planar and stacked structures highlights the importance of accurately modeling the dispersion forces in the functionals to produce reliable predictions.

The Structure of DNA Fragments: Quantum-Chemical Modelling

In this review, we analyze and systematize our computational studies of the nucleic acid duplex formations and thermodynamic stability under the different factors of investigation. The proposed structural models of mini-helix contains N nucleobase pairs (N = 3-5); QM structural data suggest that the helical conformations of mini-helix adopt geometrical parameters comparable to those of natural A- and B-DNA forms under specific conditions as micro hydration and charge compensation. The gas-phase models adopt non regular conformations between the helical form and a ladder form.. The natural helical shape of DNA mini-helix is stabilized by the presence of counterions or by explicit micro-hydration of the major and minor groves. The presence of aqueous solution is shown as a minor factor for the helical shape formation. The studies are performed at the level of density functional theory.

The Importance of Charge Transfer and Solvent Screening in the Interactions of Backbones and Functional Groups in Amino Acid Residues and Nucleotides

Quantum mechanical (QM) calculations at the level of density-functional tight-binding are applied to a protein–DNA complex (PDB: 2o8b) consisting of 3763 atoms, averaging 100 snapshots from molecular dynamics simulations. A detailed comparison of QM and force field (Amber) results is presented. It is shown that, when solvent screening is taken into account, the contributions of the backbones are small, and the binding of nucleotides in the double helix is governed by the base–base interactions. On the other hand, the backbones can make a substantial contribution to the binding of amino acid residues to nucleotides and other residues. The effect of charge transfer on the interactions is also analyzed, revealing that the actual charge of nucleotides and amino acid residues can differ by as much as 6 and 8% from the formal integer charge, respectively. The effect of interactions on topological models (protein -residue networks) is elucidated.

Recognition in the Domain of Molecular Chirality: From Noncovalent Interactions to Separation of Enantiomers

It is not a coincidence that both chirality and noncovalent interactions are ubiquitous in nature and synthetic molecular systems. Noncovalent interactivity between chiral molecules underlies enantioselective recognition as a fundamental phenomenon regulating life and human activities. Thus, noncovalent interactions represent the narrative thread of a fascinating story which goes across several disciplines of medical, chemical, physical, biological, and other natural sciences. This review has been conceived with the awareness that a modern attitude toward molecular chirality and its consequences needs to be founded on multidisciplinary approaches to disclose the molecular basis of essential enantioselective phenomena in the domain of chemical, physical, and life sciences. With the primary aim of discussing this topic in an integrated way, a comprehensive pool of rational and systematic multidisciplinary information is provided, which concerns the fundamentals of chirality, a description of noncovalent interactions, and their implications in enantioselective processes occurring in different contexts. A specific focus is devoted to enantioselection in chromatography and electromigration techniques because of their unique feature as "multistep" processes. A second motivation for writing this review is to make a clear statement about the state of the art, the tools we have at our disposal, and what is still missing to fully understand the mechanisms underlying enantioselective recognition.

Structure-activity features of purines and their receptors: implications in cell physiopathology

The purine molecular structure consists of fused pyrimidine and imidazole rings. Purines are main pieces that conform the structure of nucleic acids which rule the inheritance processes. Purines also work as metabolic intermediates in different cell functions and as messengers in the signaling pathways throughout cellular communication. Purines, mainly ATP and adenosine (ADO), perform their functional and pharmacological properties because of their structural/chemical characteristics that make them either targets of mutagenesis, mother frameworks for designing molecules with controlled effects (e.g. anti-cancer), or chemical donors (e.g., of methyl groups, which represent a potential chemoprotective action against cancer). Purines functions also come from their effect on specific receptors, channel-linked and G-protein coupled for ATP, and exclusively G-coupled receptors for ADO (also known as ADORAs), which are involved in cell signaling pathways, there, purines work as chemical messengers with autocrine, paracrine, and endocrine actions that regulate cell metabolism and immune response in tumor progression which depends on the receptor types involved in these signals. Purines also have antioxidant and anti-inflammatory properties and participate in the cell energy homeostasis. Therefore, purine physiology is important for a variety of functions relevant to cellular health; thus, when these molecules present a homeostatic imbalance, the stability and survival of the cellular systems become compromised.

Exploring the Influence of Intermolecular Interactions in Prebiotic Chemistry Using Laser Spectroscopy and Calculations

One of the most fascinating questions in chemistry is why nature chose CGAT as the alphabet of life. Very likely, such selection was the result of multiple factors and a long period of refinement. Here, we explore how the intermolecular interactions influenced such process, by characterizing the formation of dimers between adenine, theobromine and 4-aminopyrimidine. Using a combination of mass-resolved excitation spectroscopy and DFT calculations, we determined the structure of adenine-theobromine and 4-aminopyrimidine-theobromine dimers. The binding energy of these dimers is very close to the canonical adenine-thymine nucleobases. Likewise, the dimers are able to adopt Watson-Crick conformations. These findings seem to indicate that there were many options available to build the first versions of the informational polymers, which also had to compete with other molecules, such as 4-aminopyrimidine, which does not have a valid attaching point for a saccharide. For some reason, nature did not select the most strongly-bonded partners or if it did, such proto-bases were later replaced by the nowadays canonical CGAT.

Noncovalent Interactions by the Quantum Monte Carlo Method: Strong Influence of Isotropic Jastrow Cutoff Radii

We present a paradigmatic example of a strong effect of Jastrow cutoff radii setup on the accuracy of noncovalent interaction energy differences within one-determinant Slater-Jastrow fixed-node diffusion Monte Carlo (1FNDMC) simulations using isotropic Jastrow terms and effective-core potentials. Analysis of total energies, absolute and relative errors, and local energy variance of energy differences vs the reference results suggests a simple procedure to marginalize the related biases. The presented data showcase improvements in dispersion-bounded systems within such a 1FNDMC method.

Arene-Perfluoroarene Interactions in Solution

A systematic study of arene-perfluoroarene interactions in solution is presented. Using a combination of NMR titration experiments, X-ray crystallography, and computational analysis, we analyze the effects of fluorination, substituents, ring size, and solvation on the arene-perfluoroarene interaction. We find that fluorination, extension of the π systems, and enhancement of solvent polarity greatly stabilize the stacking energy up to 3 orders of magnitude (K_a = <1 to 6000 M^-1), with the highest K_a achieved for the interaction of water-soluble variants of perfluoronaphthalene and anthracene in buffered D₂O (pD = 12). Combining computational and experimental results, we conclude that this impressive binding energy is a result of enthalpically favorable electrostatic and dispersion interactions as well as the entropically driven hydrophobic effect. The enhanced understanding of arene-perfluoroarene interactions in aqueous solution sets the stage for the implementation of this abiotic intermolecular interaction in biology and medicine.

W-RESP: Well-Restrained Electrostatic Potential-Derived Charges. Revisiting the Charge Derivation Model

Representation of electrostatic interactions by a Coulombic pairwise potential between atom-centered partial charges is a fundamental and crucial part of empirical force fields used in classical molecular dynamics simulations. The broad success of the AMBER force-field family originates mainly from the restrained electrostatic potential (RESP) charge model, which derives partial charges to reproduce the electrostatic field around the molecules. However, the description of the electrostatic potential around molecules by standard RESP may be biased for some types of molecules. In this study, we modified the RESP charge derivation model to improve its description of the electrostatic potential around molecules and thus electrostatic interactions in the force field. In particular, we reoptimized the atomic radii for definition of the grid points around the molecule, redesigned the restraining scheme, and included extra point (EP) charges. The RESP fitting was significantly improved for aromatic heterocyclic molecules. Thus, the suggested W-RESP(-EP) charge derivation model shows some potential for improving the performance of the nucleic acid force fields, for which the poor description of nonbonded interactions, such as the underestimated stability of base pairing, is well-established. We also report some preliminary simulation tests (around 1 ms of simulation data) on A-RNA duplexes, tetranucleotides, and tetraloops. The simulations reveal no adverse effects, while the description of base-pairing interactions might be improved. The new charges can thus be used in future attempts to improve the nucleic acid simulation force fields, in combination with reparametrization of the other terms.

Replacing thymine with a strongly pairing fifth Base: A combined quantum mechanics and molecular dynamics study

The non-natural ethynylmethylpyridone C-nucleoside (W), a thymidine (T) analogue that can be incorporated in oligonucleotides by automated synthesis, has recently been reported to form a high fidelity base pair with adenosine (A) and to be well accommodated in B-DNA duplexes. The enhanced binding affinity for A of W, as compared to T, makes it an ideal modification for biotechnological applications, such as efficient probe hybridization for the parallel detection of multiple DNA strands. In order to complement the experimental study and rationalize the impact of the non-natural W nucleoside on the structure, stability and dynamics of DNA structures, we performed quantum mechanics (QM) calculations along with molecular dynamics (MD) simulations. Consistently with the experimental study, our QM calculations show that the A:W base pair has an increased stability as compared to the natural A:T pair, due to an additional CH-π interaction. Furthermore, we show that mispairing between W and guanine (G) causes a distortion in the planarity of the base pair, thus explaining the destabilization of DNA duplexes featuring a G:W pair. MD simulations show that incorporation of single or multiple consecutive A:W pairs in DNA duplexes causes minor changes to the intra- and inter-base geometrical parameters, while a moderate widening/shrinking of the major/minor groove of the duplexes is observed. QM calculations applied to selected stacks from the MD simulations also show an increased stacking energy for W, over T, with the neighboring bases.

Intrinsic dynamic and static nature of each HB in the multi-HBs between nucleobase pairs and its behavior, elucidated with QTAIM dual functional analysis and QC calculations

The intrinsic dynamic and static nature of each HB in the multi-HBs between nucleobase pairs (Nu-Nu') is elucidated with QTAIM dual functional analysis (QTAIM-DFA). Perturbed structures generated using coordinates derived from the compliance constants (C ii ) are employed for QTAIM-DFA. The method is called CIV. Two, three, or four HBs are detected for Nu-Nu'. Each HB in Nu-Nu' is predicted to have the nature of CT-TBP (trigonal bipyramidal adduct formation through charge transfer (CT)), CT-MC (molecular complex formation through CT), or t-HBwc (typical HB with covalency), while the vdW nature is predicted for the C-H⋯X interactions, for example. Energies for the formation of the pairs (ΔE) are linearly correlated with the total values of C ii -1 in Nu-Nu'. The total C ii -1 values are obtained by summing each C ii -1 value, similarly to the case of Ohm's law for the parallel connection in the electric resistance. The total ΔE value for a nucleobase pair could be fractionalized to each HB, based on each C ii -1 value. The perturbed structures generated with CIV are very close to those generated with the partial optimization method, when the changes in the interaction distances are very small. The results provide useful insights for better understanding DNA processes, although they are highly enzymatic.

DNA-specific selectivity in pairing of model nucleobases in the solid state

Solid-state may serve as the reaction medium for selective recognition between model nucleobases.

Binding of metal ions and water molecules to nucleic acid bases: the influence of water molecule coordination to a metal ion on water-nucleic acid base hydrogen bonds

The interactions of nucleic acid bases with non-coordinated and coordinated water molecules were studied by analyzing data in the Protein Data Bank (PDB) and by quantum chemical calculations. The analysis of the data in the crystal structures from the PDB indicates that hydrogen bonds involving oxygen or nitrogen atoms of nucleic acid bases and water molecules are shorter when water is bonded to a metal ion. These results are in agreement with the quantum chemical calculations on geometries and interaction energies of hydrogen bonds; the calculations on model systems show that hydrogen bonds of nucleic acid bases with water bonded to a metal ion are stronger than hydrogen bonds with non-coordinated water. These calculated values are similar to the strength of hydrogen bonds between nucleic acid bases. The results presented in this paper may be relevant to understand the role of water molecules and metal ions in the process of replication and stabilization of nucleic acids and also to understand the possible toxicity of metal ion interactions with nucleic acids.

Synthesis, characterization and crystal structures of novel fluorinated di(thiazolyl)benzene derivatives

A series of 11 novel fluorinated and non-fluorinated di(thiazolyl)benzenes have been synthesized via microwave assisted Stille coupling and characterized using X-ray crystallography.

Does the stability of the stacking motif surpass the planar motif in 2-amino-4-nitrophenol? - a CCSD(T) analysis

In this work we analyzed O-H...O, O-H...N, and N-H...O contacts existing in the 2-amino-4-nitrophenol structure engaged in ANP molecules through quantum chemical methods. Furthermore, the above contacts were favored to comprehensively understand the stability of noncovalent interactions, π stacking and hydrogen bonding, surviving in 2-amino-4-nitrophenol. The geometries of π stacking and hydrogen bond interactions between two 2-amino-4-nitrophenols were optimized at BLYP-D3/def2-QZVP with dispersion 3 and MP2/cc-pVTZ levels of theory, and their stability was compared using the CCSD(T) interaction energies. The analyses predicted a particularly strong π stacking interaction of 2-amino-4-nitrophenol with hydrogen bond due to the narrow equivalent configuration of NO2 interactions with the other 2-amino-4-nitrophenols. Furthermore, this work focused on analyzing the stability of the individual hydrogen bonds existing in planar and stacked arrangements. Graphical abstract Stacked and planar motif in 2-amino-4-nitrophenol.

Association of Nucleobases in Hydrated Ionic Liquid from Biased Molecular Dynamics Simulations

We employed metadynamics-based classical molecular dynamics simulations to methylated adenine-thymine (mA-mT) and guanine-cytosine (mG-mC) base pairs to see favorable conformations in various concentrations of hydrated 1-ethyl, 3-methyl imidazolium acetate. We investigated various stacked and hydrogen-bonded conformations of association of base pairs through appropriately chosen collective variables. Stacked conformations more favored in water for both base pairs, whereas Watson-Crick (WC) hydrogen-bonding conformations are favored in pure and hydrated ionic liquids (ILs) except for 0.75 mol fraction IL. We observe that EMIm cations surround the base pairs in WC conformations creating a kind of hydrophobic cavity and protect the hydrogen bonds between base pairs. However, the five-membered heteroaromatic rings of cations stack with the nucleobases in the cation-base-cation (π-π-π) model, which resembles the base-base-base stacking in a DNA duplex. Interestingly, from additional simulations of 0.5 mol fraction hydrated choline dihydrogen phosphate IL, we observe that the stacked conformations become more favored than the WC conformation due to the absence of π-bonds in cations. The calculated values of relative solubility of base pairs in pure and hydrated ionic liquids compared to those in pure water correlate well with the free energy values of WC and stacked conformations.

Analogues of P and Z as Efficient Artificially Expanded Genetic Information System

To artificially expand the genetic information system and to realize artificial life, it is necessary to discover new functional DNA bases that can form stable duplex DNA and participate in error-free replication. It is recently proposed that the 2-amino-imidazo[1,2- a]-1,3,5-triazin-4(8 H)one (P) and 6-amino-5-nitro-2(1 H)-pyridone (Z) would form a base pair complex, which is more stable than that of the normal G-C base pair and would produce an unperturbed duplex DNA. Here, by using quantum chemical calculations in aqueous medium, it is shown that the P and Z molecules can be modified with the help of electron-withdrawing and -donating substituents mainly found in B-DNA to generate new bases that can produce even more stable base pairs. Among the various bases studied, P3, P4, Z3, and Z5 are found to produce base pairs, which are about 2-15 kcal/mol more stable than the P-Z base pair. It is further shown that these base pairs can be stacked onto the G-C and A-T base pairs to produce stable dimers. The consecutive stacking of these base pairs is found to yield even more stable dimers. The influence of charge penetration effects and backbone atoms in stabilizing these dimers are also discussed. It is thus proposed that the P3, P4, Z3, and Z5 would form promiscuous artificial genetic information system and can be used for different biological applications. However, the evaluations of the dynamical effects of these bases in DNA-containing several nucleotides and the efficacy of DNA polymerases to replicate these bases would provide more insights.

The influence of substituents and the environment on the NMR shielding constants of supramolecular complexes based on A-T and A-U base pairs

In the present study, we have theoretically analyzed supramolecular complexes based on the Watson-Crick A-T and A-U base pairs using dispersion-corrected density functional theory (DFT). Hydrogen atoms H8 and/or H6 in the natural adenine and thymine/uracil bases were replaced, respectively, by substituents X8, Y6 = NH^-, NH₂, NH₃⁺ (N series), O^-, OH, OH₂⁺ (O series), F, Cl or Br (halogen series). We examined the effect of the substituents on the hydrogen-bond lengths, strength and bonding mechanism, and the NMR shielding constants of the C2-adenine and C2-thymine/uracil atoms in the base pairs. The general belief in the literature that there is a direct connection between changes in the hydrogen-bond strength and the C2-adenine shielding constant is conclusively rejected by our computations.

Deprotonated guanine·cytosine and 9-methylguanine·cytosine base pairs and their "non-statistical" kinetics: a combined guided-ion beam and computational study

We report a guided-ion beam mass spectrometric study on collision-induced dissociation (CID) of deprotonated guanine(G)·cytosine(C) base pairs and their 9-methylguanine (9MG) analogue with Xe, including measurements of product cross sections as a function of collision energy and determination of dissociation thresholds. DFT, RI-MP2 and DLPNO-CCSD(T) calculations and Rice-Ramsperger-Kassel-Marcus (RRKM) modeling were performed to elucidate structures and kinetics. The experiment and theoretical study have provided considerable insight into tautomerization, intra-base-pair proton transfer and dissociation of deprotonated G·C and 9MG·C. In contrast to the previously reported lowest-energy deprotonated base pair structure G·[C-H1]^- that consists of H-bonded neutral guanine and N1-deprotonated cytosine, we found that proton transfer from guanine N1 to cytosine N3 within G·[C-H1]^- (or 9MG·[C-H1]^-) leads to another slightly more stable conformer denoted as G·[C-H1]^-_PT1 (or 9MG·[C-H1]^-_PT1). The conventional (non-proton-transferred) and the proton-transferred conformers are close in energy and interconvert quickly, but they can be distinguished by dissociation products. The conventional structure dissociates into deprotonated cytosine and neutral guanine, while the other dissociates into deprotonated guanine and neutral cytosine. The two dissociation asymptotes have similar threshold energies, but surprisingly the CID product mass spectra of deprotonated G·C and 9MG·C are both overwhelmingly dominated by deprotonated G or 9MG, with their branching ratios greater than RRKM predictions by one to two orders of magnitude. The proton-transferred structures of deprotonated base pairs and the "unexpected" non-statistical kinetics provide new leads for understanding purine-pyrimidine interactions, forming rare nucleobase tautomers, and base pair opening.

S⋯S and S⋯P chalcogen bonding in solution: a cryospectroscopic study of the complexes of 2,2,4,4-tetrafluoro-1,3-dithietane with dimethyl sulfide and trimethylphosphine

Experimental evidence on the formation of S⋯S and P⋯S chalcogen bonded complexes between 2,2,4,4-tetrafluoro-1,3-dithiethane and the Lewis bases dimethyl sulfide and trimethylphosphine is obtained using infrared spectroscopy of solutions in liquid krypton.

Structure and Electronic Properties of Unnatural Base Pairs: The Role of Dispersion Interactions

Recent reports of the successful incorporation of unnatural base pairs (UBPs), such as d5SICS-dNaM, in the gene sequence and replication with DNA is an important milestone in synthetic biology. Followed by this, several other UBPs, such as dTPT3-dNaM, dTPT3-dFIMO, dTPT3-IMO, dTPT3-FEMO, FTPT3-NaM, FTPT3-FIMO, FTPT3-IMO, and FTPT3-FEMO, have demonstrated similar or better retention and fidelity inside cells. Of these base pairs, dNaM-dTPT3 has been optimized to be a better fit inside a pAIO plasmid. Based on both implicit and explicit dispersion-corrected density functional theory (DFT) calculations, we show that although this set of UBPs is significantly diverse in elemental and structural configuration, the members do share a common trait of favoring a slipped parallel stacked dimer arrangement. Unlike the natural bases (A, T, G, C, and U), this set of UBPs has a negligible affinity for a Watson-Crick (WC)-type planar structure because they are invariably more stable within slipped parallel stacked orientations. We also observed that all the UBPs have either similar or higher binding energies with the natural bases in similar stacked orientations. When arranged between two natural base pairs, the UBPs exhibited a binding energy similar to that of three-base sequences of natural bases. Our computational data show that the most promising base pairs are 5SICS-NaM, TPT3-NaM, and TPT3-FEMO. These results are consistent with recent progress on experimental research into UBPs along with our previous calculations on the d5SICS-dNaM pair and, therefore, strengthen the hypothesis that hydrogen bonding might not be absolutely essential and that interbase stacking dispersion interactions play a key role in the stabilization of genetic materials.

From an Easily Accessible Pentacarbonylcobalt(I) Salt to Piano-Stool Cations [(arene)Co(CO)2 ]. Chemistry 2017;23:14658-14664. [PMID: 28796933 DOI: 10.1002/chem.201703589]

The facile synthesis of a pentacarbonyl cobalt(I) salt without the need for a superacid as solvent is presented. This salt, [Co(CO)₅ ]⁺ [Al(OR^F )₄ ]^- {R^F =C(CF₃ )₃ }, readily accessible on a multigram scale, undergoes substitution reactions with arenes yielding the hitherto unknown class of two-legged cobalt piano-stool complexes [(arene)Co(CO)₂ ]⁺ with four different arene ligands. Such a substitution chemistry would have been impossible in superacid solution, as the arenes used would have been oxidized and/or protonated. Thus, the general approach described herein may have a wide synthetic use. Additionally, the thermochemistry of the piano-stool complexes is shown to be not easy to describe computationally and most of the established DFT methods overestimate the reaction energies. Only CCSD(T) calculations close to the basis set limit gave energies fully agreeing with the experiment.

Conformational flexibility and base-pairing tendency of the tobacco carcinogen O6-[4-oxo-4-(3-pyridyl)butyl]guanine

The present work uses DFT calculations to characterize the conformational and hydrogen-bonding properties of O6-[4-oxo-4-(3-pyridyl)butyl]guanine (POB-G), a DNA adduct caused by tobacco. POB-G is found to adopt many isoenergetic conformations that allow for discrete interactions between the bulky moiety and the adducted G and/or pairing base. The calculated structure and stability of the hydrogen-bonded pairs between the Watson-Crick or Hoogsteen face of POB-G and the canonical DNA nucleobases fully rationalize the previously reported mutational spectra. Specifically, the stable, non-distorted pseudo-Watson-Crick POB-G:T pair explains the predominant G➔A mutations, while the stable, yet marginally distorted pairs between the Watson-Crick face of POB-G and A or C clarify the G➔T mutations and non-mutagenic replication. Finally, the stable, yet highly distorted Hoogsteen POB-G:G pair rationalizes the experimentally-observed insertion but lack of persistence of G opposite POB-G in DNA. Overall, these structural insights are critical for guiding future studies that strive to fully understand the adduct mutagenicity, including the accessible conformations and the replication of POB-G-adducted DNA.

Toward an Expanded Genome: Structural and Computational Characterization of an Artificially Expanded Genetic Information System

Although the fundamental properties of DNA as first proposed by Watson and Crick in 1953 provided a basic understanding of how duplex DNA was organized and might be replicated, it was not until the first crystal structures of DNA (Z-DNA in 1979, B-DNA in 1980, and A-DNA in 1982) that the true complexity of the molecule began to be appreciated. Many crystal structures of oligonucleotides have since shed light on the helical forms that "Watson-Crick" DNA can adopt, their associated groove widths, and the properties of the nucleobase pairs and their interactions in all three helical forms. Additional understanding of the properties of Watson-Crick DNA has been provided by computational studies employing a variety of theoretical methods. Together with these studies devoted to understanding Watson-Crick DNA, recent efforts to expand the genetic alphabet have founded a new field in synthetic biology. One of these efforts, the artificially expanded genetic information system (AEGIS) developed by Steven Benner and co-workers, takes advantage of orthogonal hydrogen bonding to produce DNA comprised of six nucleobase pairs, of which the most extensively studied is referred to as P:Z with P being 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one) and Z being 6-amino-5-nitro-2(1H)-pyridone. P:Z forms three edge-on hydrogen bonds that differ from standard Watson-Crick pairs in the arrangement of acceptors and donor groups; P presents acceptor, acceptor, donor, and Z presents donor, donor, acceptor. Z is unique among the AEGIS nucleobases in having a nitro group present in the major groove. PZ-containing DNA has been exploited in a number of clinical applications and is being used to develop receptors and catalysts. Ultimately, the grand challenge will be to create a semisynthetic organism with an expanded genome. Furthermore, just as our understanding of the properties of natural DNA have benefited from structural and computational characterization, so too will our understanding of artificial DNA. This Account focuses on the structural and biophysical properties of AEGIS DNA containing P:Z pairs. We begin with the fundamental properties of P:Z nucleobase pairs, including their electrostatic potential and hydrogen-bonding energies, as elucidated by quantum mechanical calculations. We then examine the impact of including multiple consecutive P:Z pairs into duplex DNA providing an opportunity to investigate stacking interactions between P:Z pairs. The self-complementary 5'-CTTATPPTAZZATAAG was crystallized in B-form using the host-guest system along with analogous natural sequences including Gs or As. Use of the host-guest system to characterize B-DNA obviates a number of limitations on the structural characterization of sequences of interest; these include the ability to crystallize the desired sequences and to distinguish structural effects imparted by the lattice constraints from those inherent in the sequence itself. On the other hand, 3/6ZP, 5'-CTTATPPPZZZATAAG, was crystallized in A-form in a DNA-only lattice allowing a comparative analysis of P:Z pairs in two of the biologically relevant helical forms: A- and B-DNA. Computational studies on the 3/6ZP sequence starting in A-form provide additional evidence for a more energetically favorable stacking interaction, which we term the "slide" conformer, observed in the A-form crystal structure; this unusual stacking interaction plays a major role in altering the conformational dynamics observed for the PZ-containing duplex as compared to a GC-containing "control" duplex in long time scale molecular dynamics simulations. This combined use of structural and computational strategies paves the way for obtaining a detailed description of artificial DNA, both in how it differs from Watson-Crick DNA and in the rational discovery of proteins, such as endonucleases, transcription factors, and polymerases, which can specifically manipulate DNA containing AEGIS nucleobase pairs.

Effect of Fluorination on the Competition of Halogen Bonding and Hydrogen Bonding: Complexes of Fluoroiodomethane with Dimethyl Ether and Trimethylamine

To further rationalize the competition between halogen and hydrogen bonding, a combined experimental and theoretical study on the weakly bound molecular complexes formed between the combined halogen bond/hydrogen bond donor fluoroiodomethane and the Lewis bases dimethyl ether and trimethylamine (in standard and fully deuterated form) is presented. The experimental data are obtained by recording infrared and Raman spectra of mixtures of the compounds in liquid krypton, at temperatures between 120 and 156 K. The experiments are supported by ab initio calculations at the MP2/aug-cc-pVDZ-PP level, statistical thermodynamics and Monte Carlo free energy perturbation calculations. For the mixtures containing fluoroiodomethane and dimethyl ether a hydrogen-bonded complex with an experimental complexation enthalpy of -7.0(2) kJ mol^-1 is identified. Only a single weak spectral feature is observed which can be tentatively assigned to the halogen-bonded complex. For the mixtures involving trimethylamine, both halogen- and hydrogen-bonded complexes are observed, the experimental complexation enthalpies being -12.5(1) and -9.6(2) kJ mol^-1 respectively. To evaluate the influence of fluorination on the competition between halogen and hydrogen bonding, the results obtained for fluoroiodomethane are compared with those of a previous study involving difluoroiodomethane.

Consecutive non-natural PZ nucleobase pairs in DNA impact helical structure as seen in 50 μs molecular dynamics simulations

Z

Little is known about the influence of multiple consecutive 'non-standard' ( , 6-amino-5-nitro-2(1H)-pyridone, and , 2-amino-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one) nucleobase pairs on the structural parameters of duplex DNA. nucleobase pairs follow standard rules for Watson-Crick base pairing but have rearranged hydrogen bonding donor and acceptor groups. Using the X-ray crystal structure as a starting point, we have modeled the motions of a DNA duplex built from a self-complementary oligonucleotide (5΄-CTTATPPPZZZATAAG-3΄) in water over a period of 50 μs and calculated DNA local parameters, step parameters, helix parameters, and major/minor groove widths to examine how the presence of multiple, consecutive nucleobase pairs might impact helical structure. In these simulations, the -containing DNA duplex exhibits a significantly wider major groove and greater average values of stagger, slide, rise, twist and h-rise than observed for a 'control' oligonucleotide in which nucleobase pairs are replaced by . The molecular origins of these structural changes are likely associated with at least two differences between and . First, the electrostatic properties of differ from in terms of density distribution and dipole moment. Second, differences are seen in the base stacking of pairs in dinucleotide steps, arising from energetically favorable stacking of the nitro group in with π-electrons of the adjacent base.

Taking the halogen bonding–hydrogen bonding competition one step further: complexes of difluoroiodomethane with trimethylphosphine, dimethyl sulfide and chloromethane

To rationalize the driving factors in the competition of halogen bonding and hydrogen bonding, the complexes of the combined halogen-/hydrogen-bond donor difluoroiodomethane with the Lewis bases trimethylphosphine, dimethyl sulfide and chloromethane are studied. For all Lewis bases,ab initiocalculations lead to halogen- and hydrogen-bonded complexes. Fourier transform–IR experiments involving solutions of mixtures of difluoroiodomethane with trimethylphosphine(-d9) or dimethyl sulfide(-d6) in liquid krypton confirm the coexistence of a halogen-bonded and hydrogen-bonded complex. Also for solutions containing chloromethane, evidence of the formation of binary associations is found, but no definitive assignment of the multiple complex bands could be made. Using van't Hoff plots, the experimental complexation enthalpies for the halogen- and hydrogen-bonded complex of difluoroiodomethane with trimethylphosphine are determined to be −15.4 (4) and −10.5 (3) kJ mol−1, respectively, while for the halogen- and hydrogen-bonded complexes with dimethyl sulfide, the values are −11.3 (5) and −7.7 (6) kJ mol−1, respectively. The experimental observation that for both trimethylphospine and dimethyl sulfide the halogen-bonded complex is more stable than the hydrogen-bonded complex supports the finding that softer Lewis bases tend to favor iodine halogen bonding over hydrogen bonding.

Competition between stacked and hydrogen bonded structures of cytosine aggregates

The four bases of DNA constitute what is known as the “alphabet of life”.

Solvent effects on the excited state characteristics of adenine–thymine base pairs

A systematic analysis of the excited state characteristics of the DNA base pair adenine–thymine in stacked and Watson–Crick hydrogen bonded configurations has been carried out in this study.

Thymine cocrystals based on DNA-inspired binding motifs

Cocrystal design through DNA-type base pairing between a nucleic acid and complementary heterocycles.