Base Pair Complexes Research Articles

Electronic Properties of Artificial Metal‐DNA Base Pair Complexes Formed from Hydroxypyridone Base

AbstractThe metal‐mediated hydroxypyridone base pairs containing Ni2+, Cu2+, Zn2+, Pd2+, Cd2+, Pt2+, and Hg2+ ions were designed theoretically using by density functional theory (DFT) method in Gaussian 09 program. The calculations which were performed in the gas phase, methanol, water, and formamide were done at B3LYP functional with LANL2DZ and 6–311++G (d, p) basis sets which were used metals and rest atoms, respectively. The Gibbs free energies of the complexes were calculated and their relative stabilities were determined in the mentioned fields. Furthermore, geometric parameters, HOMO, LUMO, energy gap energies of the compounds were computed in all media and effect on these parameters of solvents were determined. Based on these parameters, the relative conductivities of the investigated molecules were predicted. On the other hand, using Atomistix ToolKit (ATK‐VNL) 2016 program, nonequilibrium Green's function (NEGF) calculations based on the combination DFT of the mentioned complexes were carried out. The transport properties of the molecules were found in detailed via NEGF‐DFT method.

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.

Exploration and characterization of the memcapacitor and memristor properties of Ni\u2013DNA nanowire devices

A 2-μm-long Ni ion-chelated DNA molecule (Ni–DNA) was found for the first time to possess both memcapacitor and memristor properties; this Ni–DNA molecule is known as a dual memory circuit element (memelement). As a memelement, the state of impedance on Ni–DNA is proportional to the unit number of Ni ions containing a base pair complex (Ni–bp), which is determined by the previously applied external voltage. Interestingly, the impedances of Ni–DNA change in response to a change in the sweeping frequencies of the external bias. Our simulation results also indicate that changes in the effective resistance and capacitance of Ni–bp may be attributed to changes in the Ni ion redox species in the Ni–bp of a Ni–DNA nanowire. Therefore, the working mechanism of a nanowire-type memcapacitor and memristor is revealed. In summary, the Ni–DNA nanowire is shown to be a multi-dimensional memory device, whose memory state depends on the length of DNA and applied external voltages/frequencies.

DFT study on metal-mediated uracil base pair complexes

The most stable of metal-mediated uracil base pair complexes were determined. Method was used density functional theory, B3LYP. The calculations of systems containing C, H, N, O were described by 6-311++G(d,p) and cc-PVTZ basis sets and LANL2DZ and SDD basis sets was used for transition metals. Then Egap values of complexes were calculated and the electrical conductivity of the complexes for single nanowires was studied by band theory. Metal-mediated uracil base pair complexes which will be used as conductive wires in nanotechnology were predicted. In nanoworld, this study is expected to show a way for practical applications.

The solvent (water) and metal effects on HOMO-LUMO gaps of guanine base pair: A computational study

The most stable structures of guanine dimer and metal-mediated guanine base pair complexes were determined both in vacuum and solvent (water). Density functional theory (DFT) method is generally used in the calculations. The calculations of systems containing C, H, N, O is described by 6-311++G(d,p) basis set and LANL2DZ basis set is used for transition metals. Some geometrical parameters, the LUMO, HOMO energy levels and energy gaps of compounds are calculated. Hydrogen bond stabilization energies in dimer and charge transfer in metal complexes are also determined using NBO analysis. The solvent and metal effects on energy gaps of molecules are found. Cu-mediated guanine complex can be preferred for nanotechnological applications due to its good electronic properties.

A multi-step model for facilitated unwinding of the yeast U4/U6 RNA duplex.

The small nuclear RNA (snRNA) components of the spliceosome undergo many conformational rearrangements during its assembly, catalytic activation and disassembly. The U4 and U6 snRNAs are incorporated into the spliceosome as a base-paired complex within the U4/U6.U5 small nuclear ribonucleoprotein (tri-snRNP). U4 and U6 are then unwound in order for U6 to pair with U2 to form the spliceosome's active site. After splicing, U2/U6 is unwound and U6 annealed to U4 to reassemble the tri-snRNP. U6 rearrangements are crucial for spliceosome formation but are poorly understood. We have used single-molecule Förster resonance energy transfer and unwinding assays to identify interactions that promote U4/U6 unwinding and have studied their impact in yeast. We find that U4/U6 is efficiently unwound using DNA oligonucleotides by coupling unwinding of U4/U6 stem II with strand invasion of stem I. Unwinding is stimulated by the U6 telestem, which transiently forms in the intact U4/U6 RNA complex. Stabilization of the telestem in vivo results in accumulation of U4/U6 di-snRNP and impairs yeast growth. Our data reveal conserved mechanisms for U4/U6 unwinding and indicate telestem dynamics are critical for tri-snRNP assembly and stability.

Conformational stabilities of iminoallantoin and its base pairs in DNA: implications for mutagenicity.

The types of mutations induced by oxidatively damaged products of DNA are continuously in debate. For example, some biochemical studies have proposed that guanidinohydantoin (Gh) would induce exclusively G to C mutations, while other studies have predicted a mixture of various mutations including G to C, G to T and G to A. In addition to the nature of mutations, the exact reasons of these mutations are also not properly understood. It is suggested that Gh can easily isomerize to iminoallantoin (Ia) in a pH-dependent manner and the transition becomes complete at pH > 8. In order to understand Gh/Ia-induced mutations, we have here studied the role of the most stable tautomer of Ia in the R- and S-enantiomeric configurations in promoting mismatch base pair complexes in DNA by employing a density functional theoretical (DFT) approach. It is found that Ia can have 39 different possible tautomeric forms each in the R- and S-enantiomeric configurations, out of which the most stable tautomer would involve the deprotonation of the N1 atom and protonation of the N3 atom. The most stable tautomer of Ia can adopt three different rotameric conformations (Ia1, Ia2, and Ia3) of comparable stabilities. It is further revealed that these rotamers of Ia can interact with different bases of DNA in 88 different possible ways. However, the interaction of G with Ia3 in both the anti- and syn-conformations would be the most stable. It is further revealed that the base pairing patterns, binding energies and electronic environments of anti-Ia3:G and G:T complexes are similar. In addition to this, it is also found that the binding patterns and energies of Gh1:G and Ia3:G complexes are similar. Based on these results, it is proposed that under physiological conditions, Gh1 may be responsible for the observed G to C mutations in DNA, while in an acidic environment Ia3 may be responsible for the same mutations. This study has led to a solid foundation for further high resolution structural studies to completely unravel Ia-induced mutagenicity in DNA.

Modeling Systems with π–π Interactions Using the Hartree–Fock Method with an Empirical Dispersion Correction

The accuracy of the Hartree-Fock method with an empirical dispersion correction, HF-D3, to model interaction energies and locate constrained minimum geometries is tested against more conventional correlation methods, such as second-order perturbation theory and coupled cluster theory, and against the sophisticated effective fragment potential model. HF-D3 was applied to substituted-benzene dimers in both sandwich and T-shaped configurations and to DNA base pair complexes in both hydrogen bonded and stacked geometries. Overall, HF-D3 is found to be a plausible and cost efficient substitute for higher levels of electronic structure theory, such as MP2, in systems with π-π interactions.

The R- and S-diastereoisomeric effects on the guanidinohydantoin-induced mutations in DNA.

Direct and indirect oxidation of guanine in DNA produces guanidinohydantoin (Gh), which is capable of inhibiting replication and inducing mutations during cellular activities. Although some biochemical studies have proposed that Gh may induce exclusively G to C mutations in DNA, other studies have predicted the occurrence of both G to C and G to T mutations. However, the exact reasons for these mutations and the dubious character of Gh in this context are not yet understood. Further, due to insufficient structural data, the electronic structure of Gh that can participate in the formation of different base pair complexes in DNA is also not known. Here, density functional theory (DFT) is used to find the most stable tautomers of Gh at the base level out of a total 112 possible tautomers and their involvement in mutagenesis is investigated by computing structures, energies and electronic properties of different base pair complexes formed between the syn- and anti-conformations of the most stable tautomer of Gh (aGh) and all the bases of DNA. It is found that aGh can coexist in R- and S-diastereoisomeric configurations. Due to the flexible guanidinium group, it can rotate about the N3-C4 bond in each of the above diastereoisomers to form two different stable conformations (aGh1 and aGh2). It is further shown that among the different base pair complexes involving aGh1, syn-aGh1:G is the most stable. This indicates that G would be easily incorporated against syn-aGh1 giving rise to G to C mutations in DNA. However, in the case of aGh2, G is the preferred base pair partner of syn-aGh2 and T is the preferred base pair partner of anti-aGh2. This implies that in addition to G to C mutations, the occurrence of aGh2 in DNA may also induce G to A mutations. Further, due to similarities between base pairing patterns and binding energies of syn-aGh1:A and syn-aGh2:A complexes with those of the T:A complex, DNA polymerases may mistakenly insert A opposite aGh1 or aGh2 by misrecognizing the latter as T. This may ultimately induce G to T mutations in DNA. However, as the constraints imposed by the DNA backbones and stacking interactions were not considered here, the possibilities of aGh2:T and aGh2:A base pairs need to be investigated experimentally. It is further found that the mutagenic character of aGh in the R- and S-diastereoisomeric forms is similar.

Application of long-range corrected density-functional theory to excess electron attachment to biomolecules

The knowledge of excess electron binding mechanisms in biomolecules is very important for molecular-level understanding of DNA damage by low-energy electrons. We here focus on two different excess electron binding mechanisms from the theoretical side; diffuse dipole-bound and valence-bound anionic states. We have found that long-range corrected (LC) density-functional theory (DFT) methods combined with a set of diffuse basis functions can reasonably describe electronic transformation processes between these two different anionic states. The applicability of the LC-DFT calculations is demonstrated for uracil, guanine–cytosine base pair, fructose, and uracil-iodide anion complex.

How DNA Polymerase X Preferentially Accommodates Incoming dATP Opposite 8-Oxoguanine on the Template

The modified base 8-oxo-7,8-dihydro-2′-deoxyguanosine (oxoG) is a common DNA adduct produced by the oxidation of DNA by reactive oxygen species. Kinetic data reveal that DNA polymerase X (pol X) from the African swine fever virus incorporates adenine (dATP) opposite to oxoG with higher efficiency than the non-damaged G:C basepair. To help interpret the kinetic data, we perform molecular dynamics simulations of pol X/DNA complexes, in which the template base opposite to the incoming dNTP (dCTP, dATP, dGTP) is oxoG. Our results suggest that pol X accommodates the oxoGsyn:A mispair by sampling closed active conformations that mirror those observed in traditional Watson-Crick complexes. Moreover, for both the oxoGsyn:A and oxoG:C ternary complexes, conformational sampling of the polymerase follows previously described large subdomain movements, local residue motions, and active site reorganization. Interestingly, the oxoGsyn:A system exhibits superior active site geometry in comparison to the oxoG:C system. Simulations for the other mismatch basepair complexes reveal large protein subdomain movement for all systems, except for oxoG:G, which samples conformations close to the open state. In addition, active site geometry and basepairing of the template base with the incoming nucleotide, reveal distortions and misalignments that range from moderate (i.e., oxoG:Asyn) to extreme (i.e., oxoGanti/syn:G). These results agree with the available kinetic data for pol X and provide structural insights regarding the mechanism by which this polymerase can accommodate incoming nucleotides opposite oxoG. Our simulations also support the notion that α-helix E is involved both in DNA binding and active site stabilization. Our proposed mechanism by which pol X can preferentially accommodate dATP opposite template oxoG further underscores the role that enzyme dynamics and conformational sampling operate in polymerase fidelity and function.

Accurate prediction of polarised high order electrostatic interactions for hydrogen bonded complexes using the machine learning method kriging

As intermolecular interactions such as the hydrogen bond are electrostatic in origin, rigorous treatment of this term within force field methodologies should be mandatory. We present a method able of accurately reproducing such interactions for seven van der Waals complexes. It uses atomic multipole moments up to hexadecupole moment mapped to the positions of the nuclear coordinates by the machine learning method kriging. Models were built at three levels of theory: HF/6-31G**, B3LYP/aug-cc-pVDZ and M06-2X/aug-cc-pVDZ. The quality of the kriging models was measured by their ability to predict the electrostatic interaction energy between atoms in external test examples for which the true energies are known. At all levels of theory, >90% of test cases for small van der Waals complexes were predicted within 1kJmol−1, decreasing to 60–70% of test cases for larger base pair complexes. Models built on moments obtained at B3LYP and M06-2X level generally outperformed those at HF level. For all systems the individual interactions were predicted with a mean unsigned error of less than 1kJmol−1.

Caffeine as base analogue of adenine or guanine: A theoretical study

The results of quantum mechanical calculations, including binding energies and results of the population analysis show that the GC and AT base pair complexes are more stable than the CAF-X ones (where CAF is caffeine and X=adenine (A), thymine (T), cytosine (C) and guanine (G)). Structural similarity between the CAF molecule and purine bases (G and A) provides the possibility of incorporation of the CAF molecule into the DNA macromolecule. By comparing the CAF-A and CAF-T complexes with the AT base pair, and the CAF-G and CAF-C complexes with the GC base pair, it was found that the formation of the CAF-T complex is more probable than the other complexes. Thus, the CAF molecule acts as an analogue base of A and can be incorporated into the DNA macromolecule and paired with T and C in normal and rare state, respectively. Indeed, the results show that formation of the CAF-C complex is less probable than the CAF-T one and an AT to GC conversion is rarely occurred in the next DNA replication, so the CAF molecule may be considered as a weak mutagenic compound. To examine solvent effect, the binding energies have been calculated in solvent for the most important structures of the CAF-G, CAF-T, CAF-A and CAF-C complexes. The results in solvent are in agreement with those in the gas phase.

Computational Studies on Effects of MDMA as an Anticancer Drug on DNA

This research is designed to further understand the effects of the novel drug MDMA on biologic receptor of DNA. The ultimate goal is to design drugs that have higher affinity with DNA. Understanding the physicochemical properties of the drug as well as the mechanism by which it interacts with DNA should ultimately enable the rational design of novel anticancer or antiviral drugs. Molecular modeling on the complex formed between MDMA and DNA presented this complex to be fully capable of participating in the formation of a stable intercalation site. Furthermore, the molecular geometries of MDMA and DNA bases (Adenine, Guanine, Cytosine, and Thymine) were optimized with the aid of the B3LYP/6-31G* method. The properties of the isolated intercalator and its stacking interactions with adenine···thymine (AT) and guanine···cytosine (GC) nucleic acid base pairs were studied with the DFTB method. DFTB method is an approximate version of the DFT method that was extended to cover the London dispersion energy. The B3LYP/6-31G* stabilization energies of the intercalator···base pair complexes were found to be -9.40 and -12.57 kcal/mol for AT···MDMA and GC···MDMA, respectively. Results from comparison of the DFTB method and HF method conclude close results and support each other.

Metal binding and substrate positioning by evolutionarily invariant U6 sequences in catalytically active protein-free snRNAs

We have previously shown that a base-paired complex formed by two of the spliceosomal RNA components, U6 and U2 small nuclear RNAs (snRNAs), can catalyze a two-step splicing reaction that depended on an evolutionarily invariant region in U6, the ACAGAGA box. Here we further analyze this RNA-catalyzed reaction and show that while the 5' and 3' splice site substrates are juxtaposed and positioned near the ACAGAGA sequence in U6, the role of the snRNAs in the reaction is beyond mere juxtaposition of the substrates and likely involves the formation of a sophisticated active site. Interestingly, the snRNA-catalyzed reaction is metal dependent, as is the case with other known splicing RNA enzymes, and terbium(III) cleavage reactions indicate metal binding by the U6/U2 complex within the evolutionarily conserved regions of U6. The above results, combined with the structural similarities between U6 and catalytically critical domains in group II self-splicing introns, suggest that the base-paired complex of U6 and U2 snRNAs is a vestigial ribozyme and a likely descendant of a group II-like self-splicing intron.

Electronic and molecular structure of M-DNA fragments

To study M-DNA molecular structure (such DNA with transition metal ions placed between the nucleic bases is able to conduct the electric current) and its conductivity mechanisms, we carried out ab initio quantum-mechanical calculations of electronic and spatial structures, thermodynamic characteristics of adenine-thymine (АТ) and guanine-cytosine (GC) base pair complexes with Zn(2+) and Ni(2+). To take into account the influence of the alkaline environment, calculations for these complexes were also carried out with hydroxyl and two water molecules. Computations were performed at MP2 level of theory using 6-31+G* basis set. Analogous calculations were carried out for (AC)(TG) stacking dimer of nucleic acid base pairs with two Zn(2+). The calculation of the interaction energy in complexes has shown the preference of locating the metal ion (instead of the imino proton) between bases in M-DNA. The electronic transition energy calculation has revealed the reduction of the first singlet transition energy in АТ and GC complexes with Ni(2+) from 4.5eV to 0.4 - 0.6eV. Ni(2+) orbitals take part in the formation of HOMO and LUMO on the complexes investigated. It was shown that charges of metal ions incorporated into complexes with nucleic bases and in dimer decrease significantly.

Protein-free small nuclear RNAs catalyze a two-step splicing reaction

Pre-mRNA splicing is a crucial step in eukaryotic gene expression and is carried out by a highly complex ribonucleoprotein assembly, the spliceosome. Many fundamental aspects of spliceosomal function, including the identity of catalytic domains, remain unknown. We show that a base-paired complex of U6 and U2 small nuclear RNAs, in the absence of the approximately 200 other spliceosomal components, performs a two-step reaction with two short RNA oligonucleotides as substrates that results in the formation of a linear RNA product containing portions of both oligonucleotides. This reaction, which is chemically identical to splicing, is dependent on and occurs in proximity of sequences known to be critical for splicing in vivo. These results prove that the complex formed by U6 and U2 RNAs is a ribozyme and can potentially carry out RNA-based catalysis in the spliceosome.

Configurations of base-pair complexes in solutions

A theoretical search for the most stable conformations (i.e., stacked or hydrogen bonded) of the base pairs A-U and G-C in water, CCl4, and CHCl3 solutions is presented. The calculations of free energies indicate a significant role of the solvent in determining the conformations of the base-pair complexes. The application of the continuum method yields preferred conformations in good agreement with experiment. Results of the calculations with this method emphasize the importance of both the electrostatic interactions between the two bases in a complex, and the dipolar interaction of the complex with the entire medium. In calculations with the solvation shell method, the last term, i.e., dipolar interaction of the complex with the entire medium, was added. With this modification the prediction of the solvation shell model agrees both with the continuum model and with experiment, i.e., in water the stacked conformation of the bases is preferred.

Theoretical investigation of interaction between Gatifloxacin and DNA: Implications for anticancer drug design

This study is an attempt to have a better understanding of the physicochemical interaction between a novel anticancer drug, Gatifloxacin (GTFX) and its biological receptor, DNA. The eventual purpose is to design drugs which interact more with DNA. Considering the physicochemical properties of the drug as well as the mechanism by which it interacts with DNA, it should eventually allow the rational design of novel anticancer or antiviral drugs. Molecular modeling on the complex formed between Gatifloxacin and DNA presented the full ability of the drug for participating in the formation of a stable intercalation site. Besides, the molecular geometries of Gatifloxacin (GTFX) and the DNA bases (Adenine, Guanine, Cytosine and Thymine) were optimized with the aid of B3LYP/6-31G⁎ method. The properties of the isolated intercalator and its stacking interactions with the adenine⋯thymine (AT) and guanine⋯cytosine (GC) nucleic acid base pairs, were studied using DFTB method, (an approximate version of the DFT method) that was extended to cover the London dispersion energy. The B3LYP/6-31G⁎ stabilization energies of the intercalator⋯base pair complexes were found to be − 26.99 and − 37.62 kcal/mol for AT⋯GTFX and GC⋯GTFX, respectively. It was completed that the dispersion energy and the electrostatic interaction contributed to the stability of the intercalator⋯DNA base pair complexes.

Prediction of Interaction Energies of Substituted Hydrogen-Bonded Watson−Crick Cytosine:Guanine8XBase Pairs

We investigated the variation in the interaction energy between the Watson-Crick hydrogen-bonded DNA base pairs guanine and cytosine (G(8X):C), where guanine is substituted in the C8 position by 37 different functional groups. Base pairs were optimized at the B3LYP/6-311+G(2d,p) level. A base pair complex containing a more strongly electron-withdrawing group remarkably forms a more stable base pair with C. Multivariate linear regression provided a quantitative relationship between the interaction energies and descriptors generated by the quantum chemical topology (QCT) approach. The descriptors were sampled from the monomers only, not the supermolecular base pair complexes. A model with r2 = 0.96 and a root-mean-square (rms) value of 0.6 kJ/mol was obtained for a training set of 28 base pair complexes. The model was tested by an external test set of 9 complexes, yielding r2 = 0.99 and an rms value of 0.2 kJ/mol. The results indicated that the bonds C6=O6 and N2-H2 at the hydrogen-bonded frontier of the guanine derivatives play an important role in transmitting the substituent effects. A linear correlation between substitution energies and Hammett constants (sigma(m)) was also obtained for all 37 substituents, yielding r2 = 0.82 and an rms value of 1.2 kJ/mol. The model based on QCT descriptors can therefore be used for the prediction of the interaction energy of the base pair G(8x):C, strictly based on data for the G(8x) monomers only.