How proton transfer impacts hachimoji DNA

Hachimoji DNA is a synthetic nucleic acid extension of DNA, formed by an additional four bases, Z, P, S, and B, that can encode information and sustain Darwinian evolution. In this paper, we aim to look into the properties of hachimoji DNA and investigate the probability of proton transfer between the bases, resulting in base mismatch under replication. First, we present a proton transfer mechanism for hachimoji DNA, analogous to the one presented by Löwdin years prior. Then, we use density functional theory to calculate proton transfer rates, tunnelling factors and the kinetic isotope effect in hachimoji DNA. We determined that the reaction barriers are sufficiently low that proton transfer is likely to occur even at biological temperatures. Furthermore, the rates of proton transfer of hachimoji DNA are much faster than in Watson-Crick DNA due to the barrier for Z-P and S-B being 30% lower than in G-C and A-T. Suggesting that proton transfer occurs more frequently in hachimoji DNA than canonical DNA, potentially leading to a higher mutation rate.

From structure and dynamics to biomolecular functions: The ubiquitous role of solvent in biology

Biological activity requires a solvent that can provide a suitable environment, which satisfies the twin need for stability and the ability to change. Among all the solvents water plays the most important role. We review, analyze, and comment on recent works on the structure and dynamics of water around biomolecules and their role in specific biological functions. While studies in the past have focused on understanding the biomolecule-water interactions through a hydration layer; recently the attention has shifted towards understanding functions at a molecular level. Such a microscopic understanding clearly requires elucidation of detailed dynamical processes where solvent molecules play an important role. Finally, we comment on the advances made in understanding the role of water inside a biological cell.

Assessing the sequence dependence of pyrimidine-pyrimidone (6-4) photoproduct in a duplex double-stranded DNA: A pitfall for microsecond range simulation

Sequence dependence of the (6-4) photoproduct conformational landscape when embedded in six 25-bp duplexes is evaluated along extensive unbiased and enhanced (replica exchange with solute tempering, REST2) molecular dynamics simulations. The structural reorganization as the central pyrimidines become covalently tethered is traced back in terms of non-covalent interactions, DNA bending, and extrusion of adenines of the opposite strands. The close sequence pattern impacts the conformational landscape around the lesion, inducing different upstream and downstream flexibilities. Moreover, REST2 simulations allow us to probe structures possibly important for damaged DNA recognition.

Purine tautomeric preferences and bond-length alternation in relation with protonation-deprotonation and alkali metal cationization

Quantum chemical calculations were carried out for deprotonated (P-) and protonated purine (PH+) and for adducts with one alkali metal cation (P-M+ and PM+, where M+ is Li+ or Na+) in the gas phase {B3LYP/6-311+G(d,p)}, a model of perfectly apolar environment, and for selected structures in aqueous solution {PCM(water)//B3LYP/6-311+G(d,p)}, a reference polar medium for biological studies. All potential isomers of purine derivatives were considered, the favored structures indicated, and the preferred sites for protonation/deprotonation and cationization reactions determined. Proton and metal cation basicities of purine in the gas phase were discussed and compared with those of imidazole and pyrimidine. Bond-length alternations in the P, PH+, P-M+, and PM+ forms were quantitatively measured using the harmonic oscillator model of electron delocalization (HOMED) indices and compared with those for P. Variations of the HOMED values when proceeding from the purine structural building blocks, pyrimidine and imidazole, to the bicyclic purine system were also examined. Generally, the isolated NH isomers exhibit a strongly delocalized π-system (HOMED > 0.8). Deprotonation slightly increases the HOMED values, whereas protonation and cationization change the HOMED indices in different way. For bidentate M+-adducts, the HOMED values are larger than 0.9 like for the largely delocalized P-. The HOMED values correlate well in a comprehensive relationship with the relative Gibbs energies (ΔG) calculated for individual isomers whatever the purine form is, neutral, protonated, or cationized. When PCM-DFT model was utilized for P-, PH+, PM+, and P-M+ (M+ = Li+) both electron delocalization and relative stability are different from those for the molecules in vacuo. The solvation effects cause a slight increase in HOMEDs, whereas the ΔEs decrease, but in different ways. Hence, contribution of particular isomers in the isomeric mixtures of PH+, PM+, and P-M+ also varies. HOMED variations for the favored neutral, deprotonated, protonated, and lithiated forms of purine in the gas phase and aqueous solution.

Non-poissonian Distribution of Point Mutations in DNA

In general, for chemical reactions occurring in systems, where fluctuations are not negligibly small, it is necessary to introduce a master equation for distribution of probability of fluctuations. It has been established that the monomolecular reactions of a type as A ↔ X are described by the master equation, which leads to a Poisson distribution with the variance equal to the average value N 0. However, the consideration of the Löwdin mechanism as autocatalytic non-linear chemical reactions such as A + X ↔ 2X and the corresponding master equation lead to a non-Poissonian probability distribution of fluctuations. In the presented work, first-order autocatalysis has been applied to the Löwdin's mechanism of spontaneous mutations in DNA. Describing double proton transfers between complimentary nucleotide bases along the chain by first-order autocatalytic reactions, the corresponding master equation for protons in tautomeric states becomes non-linear, and at non-equilibrium conditions this leads to the non-Poissonian distribution of spontaneous mutations in DNA. It is also suggested that the accumulation of large fluctuations of successive cooperative concerted protons along the chain may produce higher non-linearities which could have a significant impact on some biochemical processes, occurring in DNA.

Double Proton Transfer in the Dimer of Formic Acid: An Efficient Quantum Mechanical Scheme

Double proton transfer plays an important role in biology and chemistry, such as with DNA base pairs, proteins and molecular clusters, and direct information about these processes can be obtained from tunneling splittings. Carboxylic acid dimers are prototypes for multiple proton transfer, of which the formic acid dimer is the simplest one. Here, we present efficient quantum dynamics calculations of ground-state and fundamental excitation tunneling splittings in the formic acid dimer and its deuterium isotopologues. These are achieved with a multidimensional scheme developed by us, in which the saddle-point normal coordinates are chosen, the basis functions are customized for the proton transfer process, and the preconditioned inexact spectral transform method is used to solve the resultant eigenvalue problem. Our computational results are in excellent agreement with the most recent experiments (Zhang et al., 2017; Li et al., 2019).

Theoretical Study of Adenine to Guanine Transition Assisted by Water and Formic Acid Using Steered Molecular Dynamic Simulations

The free energy profile of the adenine to guanine transition in the gas and aqueous phases was obtained by applying steered molecular dynamic (SMD) simulations. Three processes were considered to explain the mechanism assisted by water and formic acid molecules. The first process is hydrolytic deamination of adenine, then oxidation of the hypoxanthine previously formed, and finally, the animation from xanthine to guanine. In the gas phase these processes indicate a slow and not spontaneous conversion (ΔG g = 4.07 kcal·mol-1, k = 5.59·10-40 s-1), and a lifetime for guanine of τ = 7.75·10+22 s. The presence of solvent makes the transition more difficult by increasing the reaction energy to 26.90 kcal·mol-1 and decreasing the speed of the process to 1.63·10-55 s-1. However, it decreases the energy of the deamination process to -9.63 kcal·mol-1 and the lifetime of guanine base to τ = 6.85·10+17 s when the surrounding medium used in the transition process is aqueous. The results show that the guanine could participate in genetic mutations based on the lifetimes obtained. Transition states and intermediates structures were analyzed at the molecular dynamic level. This allows to follow the mechanism over time and to calculate thermodynamic and kinetic properties.

Extended Efficacy of Low-Dose Valganciclovir for Prevention of Cytomegalovirus Disease in Intermediate-Risk Kidney Transplant Recipients: Two-Year Follow-Up

OBJECTIVES

In a previous study, we evaluated 1-year outcomes of using low-dose valganciclovir prophylaxis for cytomegalovirus infection in intermediate-risk kidney transplant recipients. Whether this effect persists in the long term is unknown. We aimed to evaluate the 2-year follow up of such adopted prophylaxis.

MATERIALS AND METHODS

We randomized 2 matched groups of kidney transplant recipients (1:1) to receive valganciclovir as 450 mg daily (group 1) or 900 mg daily (group 2) for the first 6 months after kidney transplant. The final analysis included 196 patients as intermediate-risk patients (98 in each treatment group) after exclusion of 5 high-risk patients. Serologically, all patients were at moderate risk for cytomegalovirus infection. Long-term outcomes including cytomegalovirus disease, acute rejection, new-onset diabetes after transplant, graft loss, and patient survival were assessed.

RESULTS

Through year 2 of follow-up, cytomegalovirus infection was reported in only 1 patient in group 1 (at month 13) and 1 patient in group 2 (at month 19) (not significant). Biopsy-proven acute rejection episodes were not statistically different between the groups (2 episodes in group 1 and 6 in group 2; P = .431). New-onset diabetes posttransplant was reported in 8.1% in group 1 and 13.2% in group 2 (P = .535). Graft failure was equal in both groups (1 in each group) at 2 years of follow up (not significant). Patient survival was comparable in both groups (100% in group 1 versus 97.9% in group 2; P = .661). The total number of cytomegalovirus infections at 2 years was numerically less in group 1 (P = .128).

CONCLUSIONS

Low-dose valganciclovir prophylaxis for 6 months was associated with sustained reduction of cytomegalovirus infection up to 2 years after kidney transplant without significant impact on the acute rejection, new-onset diabetes posttransplant, or patient and graft outcomes.

Effects of Positive and Negative Ionization on Prototropy in Pyrimidine Bases: An Unusual Case of Isocytosine

Intramolecular proton-transfers (prototropic conversions) have been studied for the guanine building block isocytosine (iC), and effects of positive ionization, called one-electron oxidation (iC - e → iC^+•), and negative ionization, called one-electron reduction (iC + e → iC^-•), on tautomeric conversions when proceeding from neutral to ionized isocytosine have been discussed. Although radical cations and radical anions are very short-lived species, the ionization effects could be investigated by quantum-chemical methods. Such kind of studies gives some information about the labile protons and the most basic positions in the neutral and radical forms of the tautomeric system. For investigations, the complete isomeric mixture of isocytosine has been considered and calculations performed in two extreme environments, apolar {DFT(B3LYP)/6-311+G(d,p)} and polar {PCM(water)//DFT(B3LYP)/6-311+G(d,p)}. For selected isomers, the G4 theory has also been applied. There are no good relations for energetic parameters of neutral and ionized forms. Ionization energies depend on localization of labile protons. Tautomeric equilibria for neutral and ionized isocytosine, favored sites of protonation and deprotonation, and favored structures of protonated and deprotonated forms strongly depend on environment. Acidity of iC^+• is close to that of the iC conjugate acid, and basicity of iC^-• is close to that of the iC conjugate base. This increase of acid-base properties of charged radicals explains the proton-transfer in ionized pairs of nucleobases. When compared to other pyrimidine bases such as uracil (U) and cytosine (C), which exhibit analogous tautomeric equilibria between nine prototropic tautomers as isocytosine, the tautomeric preferences for iC, iC^+•, iC^-•, U, U^+•, U^-•, C, C^+•, and C^-• are completely different. The differences suggest that acid-base properties of functional groups, their stabilities, and ionization energies play a principal role in proton-transfers for pyrimidine bases and influence compositions of tautomeric mixtures.

Adenine–thymine tautomerization under the influence of strong homogeneous electric fields

This paper reveals the key energy correction terms contributing to the low probability of AT tautomerization under electric fields.

Theoretical study of mechanisms for the hydrolytic deamination of cytosine via steered molecular dynamic simulations

Gibbs free energy profiles of the cytosine deamination assisted by a water molecule in a discrete aqueous medium were obtained by the application of Steered Molecular Dynamic (SMD) simulations. Two pathways were considered to explain the mechanism of this process, where the water molecule attacks the C–N bond to give an intermediate (an amino–hydroxy–oxo structure in the A-path, and a hydroxy–oxo in the B-path) as the determinant step of reaction. Stationary structures along both energy profiles were analyzed at molecular dynamics level, obtaining states with higher free energies than those from electronic calculations in the gas phase and in solution described as a continuous medium. From the results obtained, the more complex A-pathway, with five steps, was kinetically the most favorable (with an endergonic reaction energy of 7.41 kcal mol⁻¹, a high barrier of 67.53 kcal mol⁻¹, and a small velocity constant k₂ = 1.80 × 10⁻³⁷ s⁻¹), concluding that the uracil base can participate in a spontaneous genetic mutation since the uracil–ammonia complex has a long lifetime of 6.10 × 10²⁷ s. This process turns out exergonic and faster when carried out in gas phase simulation or electronic calculation with a continuous medium, due to the disappearance of explicit water molecules that can compete with the assistant molecule.

Gibbs free energy profiles of the cytosine deamination assisted by a water molecule in a discrete aqueous medium were obtained by the application of Steered Molecular Dynamic (SMD) simulations.

Influence of Coupling and Embedding Schemes on QM Size Convergence in QM/MM Approaches for the Example of a Proton Transfer in DNA

The influence of embedding and coupling schemes on the convergence of the QM size in the QM/MM approach is investigated for the transfer of a proton in a DNA base pair. We find that the embedding scheme (mechanical or electrostatic) has a much greater impact on the convergence behavior than the coupling scheme (additive QM/MM or subtractive ONIOM). To achieve size convergence, QM regions with up to 6000 atoms are necessary for pure QM or mechanical embedding. In contrast, electrostatic embedding converges faster: for the example of the transfer of a proton between DNA base pairs, we recommend including at least five base pairs and 5 Å of solvent (including counterions) into the QM region, i.e., a total of 1150 atoms.

Exposing the G-quadruplex to electric fields: the role played by telomeres in the propagation of DNA errors

We use quantum calculations to assess the impact of external electric fields on the stability of G-quadruplex, a key structure in telomere functionality.

Molecular phylogenetic study on few morphotypes of a patellogastropod Cellana karachiensis from northern Arabian Sea reveals unexpected genetic diversity

A group of limpets, Cellana karachiensis, exhibiting phenotypic plasticity were examined from Gujarat coastline India, using molecular phylogeny. Previous examination of the COI genes established the presence of three different haplotypes X, Y and Z, while present study showed three more haplotypes X₁, X₂ and Z₁. Thus, a total of six COI gene haplotypes, having 99.23% to 99.85% sequence similarity, were observed with variations at six sites. Bayesian phylogenetic analysis shows divergence of lineages X-Y, X₁-X₂ and Z-Z₁. Careful observation of nucleotide alterations showed a nonrandom mutation with more A↔G and C↔T transitions between closely related species of the genus Cellana. A stretch of 17 base pair, within COI region, was marked as region with high degree of variability between species of Cellana. Results suggest that this could be the beginning of speciation, with partial or complete reproductive barrier or these are already distinct species in different stages of evolution.

Next generation sequencing of benzo(a)pyrene-induced lacZ mutants identifies a germ cell-specific mutation spectrum

De novo mutations are implicated in a variety of genetic diseases and arise primarily in the male germline. We investigated whether male germ cells have unique mechanisms for spontaneous or chemically-induced mutation relative to somatic cells using the MutaMouse model. We recovered lacZ transgenes from sperm 42 days after a 28-day exposure to benzo(a)pyrene (BaP, 100 mg/kg/day) to assess mutations arising in dividing spermatogonia. BaP caused a 3.4-fold increase in lacZ mutant frequency over controls which increased to 4.1-fold after clonal correction. We then used next generation sequencing to compare the spontaneous and BaP-induced mutation spectra in sperm and bone marrow. The spontaneous spectrum in sperm had significantly more G:C to A:T transitions and fewer mutations at A:T basepairs than bone marrow. BaP predominantly induced G:C to T:A transversions in both cell types, and both were enriched for mutations at CpG dinucleotides. However, BaP induced significantly more deletions in sperm, but more G:C to A:T transitions and G:C to C:G transversions in bone marrow. Differences in error-prone translesion DNA synthesis polymerases may underlie the observed spectrum differences between sperm and bone marrow. These findings suggest that mutations in sperm can arise via mechanisms that are unique to male germ cells.

Effects of ionization on stability of 1-methylcytosine - DFT and PCM studies

Consequences of ionization were studied by quantum-chemical methods (DFT and PCM) for 1-methylcytosine (MC)—a model of the nucleobase cytosine (C) connected with sugar in DNA. For calculations, three prototropic tautomers (one amino and two imino forms) and two imino zwitterions were considered, including conformational or configurational isomerism of exo heterogroups. Ionization and interactions between neighboring groups affect intramolecular proton-transfers, geometric and thermodynamic parameters, and electron delocalization for individual isomers. We discovered that an imino isomer is present in the isomeric mixture in the highest amount for positively ionized MC. Its contribution in neutral and negatively ionized MC is considerably smaller. Acid-base parameters for selected radical ions were estimated in the gas phase and compared to those of neutral MC. Gas-phase acidity of radical cations is close to that of the conjugate acid of MC, and gas-phase basicity of radical anions is close to that of the conjugate base of MC. Various routes of amino-imino conversion between neutral and ionized isomers were considered. Energetic-barrier for intramolecular proton-transfer in MC is close to that in the parent system—formamidine.

A novel conception for spontaneous transversions caused by homo-pyrimidine DNA mismatches: a QM/QTAIM highlight

We have firstly shown that the T·T(w) and C·C(w) DNA mismatches with wobble (w) geometry stay in slow tautomeric equilibrium with short T·T*(WC) and C·C*(WC) Watson-Crick (WC) mispairs. These non-dissociative tautomeric rearrangements are controlled by the plane-symmetric, highly stable, highly polar and zwitterionic transition states. The obtained results allow us to understand in what way the T·T(w) and C·C(w) mismatches acquire enzymatically competent T·T*(WC) and C·C*(WC) conformations directly in the hydrophobic recognition pocket of a high-fidelity DNA-polymerase, thereby producing thermodynamically non-equilibrium spontaneous transversions. The simplest numerical estimation of the frequency ratio of the TT to CC spontaneous transversions satisfactorily agrees with experimental data.