Can DNA-binding proteins of replisome tautomerize nucleotide bases? Ab initio model study

Quasi-intrinsic fluorescent probes for detecting the DNA adduct (ABPdG) based on an excited-state intermolecular charge transfer mechanism

N'-(2'-Deoxyguanosin-8-yl)-4-aminobiphenyl (^ABPdG) is one of the most representative carcinogenic DNA adducts formed by human exposure to 4-aminobiphenyl (4-ABP) during dye production, rubber-manufacturing processes and cigarette smoke. Accordingly, the ultrasensitive detection of ABP-derived adducts in DNA with minimal interference to the native structures becomes key for elucidating carcinogenesis mechanisms and mitigating the risk of cancer. In view of the lack of efficient optical emission in ^ABPG, we report a theoretical study on the photophysical properties of a set of quasi-intrinsic fluorescent C-analogues, which can form stable W-C base pairs with ^ABPG. It is found that fluorophore replacement and ring-expansion can bring a red-shifted absorption and bright photoluminescence due to additional π-conjugation. In particular, because the tricyclic cytosine analogue 1,3-diaza-2-oxophenoxazine (tC^O) possesses distinct optical properties, it is proposed as a biosensor to identify ^ABPG. The TDDFT-calculated absorption maximum of tC^O is red-shifted by 97 nm in comparison with that of the native C base, which contributes to selective excitation after incorporating into the nucleic acids. Although the fluorescence is insensitive to base pairing with natural guanine, the excited state intermolecular charge transfer (ESICT)-governed "OFF-ON" signal can be observed in the presence and absence of ^ABPG. Moreover, to evaluate the direct availability of the bright C-analogues with high selectivity for the deoxyguanosine adduct ^ABPG in DNA, we further investigated thoroughly the effects of its linking to deoxyribose on its absorption and emission, which shows little difference from that of experiment.

Photophysical properties of fluorescent nucleobase P-analogues expected to monitor DNA replication

In this work, we computationally design a series of fluorescent purine analogues based on the 2-amino-8-(1'-β-D-2'-deoxyribofuranosyl)-imidazo[1,2-a]-1,3,5-triazin-4(8H)-one (P) to monitor the DNA replication process with merely a minimal perturbation to the natural structure of nucleic acid. The P-modified fluorescent probes present red-shifted absorption spectra and enhanced photoluminescence due to the additional π-conjugation resulting from the fluorophore modification and the ring-expansion. Efficient fluorescence quenching of P-analogues occurs upon pairing with the complementary 6-amino-5-nitro-3-(1'-β-D-2'-deoxyribofuranosyl)-2(1H)-pyridone (Z) due to the nonradiative relaxation from the low-lying dark excited state to the ground state of Z moiety. Especially, the P3 and the P7, which have high fluorescence intensity in both gas and liquid phases, are proposed as the sensors for studying conformational switching in the presence and absence of a complementary sequence. Also examined are the influences of hydration and the linking to deoxyribose on absorption and emission processes. Besides, the potential phosphorescence emission of these modified base pairs is taken into account by constructing the relaxed potential energy curves of S₀, T₁ and S₁ states.

Novel Tautomerisation Mechanisms of the Biologically Important Conformers of the Reverse Löwdin, Hoogsteen, and Reverse Hoogsteen G*·C* DNA Base Pairs via Proton Transfer: A Quantum-Mechanical Survey

For the first time, in this study with the use of QM/QTAIM methods we have exhaustively investigated the tautomerization of the biologically-important conformers of the G*·C* DNA base pair-reverse Löwdin G*·C*(rWC), Hoogsteen G*'·C*(H), and reverse Hoogsteen G*'·C*(rH) DNA base pairs-via the single (SPT) or double (DPT) proton transfer along the neighboring intermolecular H-bonds. These tautomeric reactions finally lead to the formation of the novel G· C O 2 * (rWC), G N 2 * · C(rWC), G*'N2·C(rWC), G N 7 * · C(H), and G*'N7·C(rH) DNA base mispairs. Gibbs free energies of activation for these reactions are within the range 3.64-31.65 kcal·mol-1 in vacuum under normal conditions. All TSs are planar structures (Cs symmetry) with a single exception-the essentially non-planar transition state TSG*·C*(rWC)↔G+·C-(rWC) (C1 symmetry). Analysis of the kinetic parameters of the considered tautomerization reactions indicates that in reality only the reverse Hoogsteen G*'·C*(rH) base pair undergoes tautomerization. However, the population of its tautomerised state G*'N7·C(rH) amounts to an insignificant value-2.3·10-17. So, the G*·C*(rWC), G*'·C*(H), and G*'·C*(rH) base pairs possess a permanent tautomeric status, which does not depend on proton mobility along the neighboring H-bonds. The investigated tautomerization processes were analyzed in details by applying the author's unique methodology-sweeps of the main physical and chemical parameters along the intrinsic reaction coordinate (IRC). In general, the obtained data demonstrate the tautomeric mobility and diversity of the G*·C* DNA base pair.

The Role of Proton Transfer on Mutations

Hydrogen bonds play a critical role in nucleobase studies as they encode genes, map protein structures, provide stability to the base pairs, and are involved in spontaneous and induced mutations. Proton transfer mechanism is a critical phenomenon that is related to the acid-base characteristics of the nucleobases in Watson-Crick base pairs. The energetic and dynamical behavior of the proton can be depicted from these characteristics and their adjustment to the water molecules or the surrounding ions. Further, new pathways open up in which protonated nucleobases are generated by proton transfer from the ionized water molecules and elimination of a hydroxyl radical in this review, the analysis will be focused on understanding the mechanism of untargeted mutations in canonical, wobble, Hoogsteen pairs, and mutagenic tautomers through the non-covalent interactions. Further, rare tautomer formation through the single proton transfer (SPT) and the double proton transfer (DPT), quantum tunneling in nucleobases, radiation-induced bystander effects, role of water in proton transfer (PT) reactions, PT in anticancer drugs-DNA interaction, displacement and oriental polarization, possible models for mutations in DNA, genome instability, and role of proton transfer using kinetic parameters for RNA will be discussed.

Drug–DNA interaction, a joint DFT-D3/MD study on safranal as an anticancer and DNA nanostructure model

In this research, using a combination of quantum mechanics and molecular dynamic (MD) simulations, the interaction of safranal (2,6,6-trimethylcyclohexa-1,3-dien-1-carboxaldehyde) as an anti-cancer drug and Dickerson B-DNA was studied. MD simulations were executed for 35 ns in water. Binding energy analysis in three definite parts of the B-DNA and comparison between different contributions of the binding energy shows that the van der Waals energy part of the interaction is impressive among the standard molecular mechanic energy terms. On the basis of Gibbs energies, it is confirmed that the most important interactions in the safranal complex are related to the A–T and C–G rich regions, which is in agreement with the experimental data. Quantum theory of atoms in molecules and natural bond orbital analyses were applied. A diminution in the electronic chemical potential of the safranal–DNA complex in comparison with the isolated DNA, 0.026 and 0.022 au for the S1 region and 0.012 and 0.017 au for the S2 region, was obtained in the gas phase and water, respectively, which increases the complex stability. An enhancement in the electrophilicity character, during the complexation process, shows the electron charge flux between the safranal and DNA, especially in water. The strengths of the CH⋯O bonds at the center of safranal–DNA interaction were also evaluated. A mean value of 0.06 au for the electron density of the bond critical point of the H⋯O in the complex confirms the H-bond formation during the complexation.

Formic acid catalyzed isomerization of protonated cytosine: a lower barrier reaction for tautomer production of potential biological importance

Tautomerism in nucleotide bases is one of the possible mechanisms of DNA mutation. In spite of numerous studies on the structure and energy of protonated cytosine tautomers, little information is available on the process of their intra- and intermolecular tautomerizations. The catalytic ability of H₂O, HCOOH, and the HCOOHH₂O group to facilitate the tautomerism of the Cyt2t⁺ to CytN3⁺ isomer has been studied. It is shown that the activation free energies of tautomerism in the gas phase are 161.17, 58.96, 26.06, and 15.69 kJ mol^-1, respectively, when the reaction is carried out in the absence and presence of H₂O, HCOOH, or the HCOOHH₂O group. The formation of a doubly hydrogen bonded transition state is central to lowering the activation free energy and facilitating the intramolecular hydrogen atom transfer that is required for isomerization. In the aqueous phase, although the solvent effects of water significantly decrease the activation free energy of intramolecular tautomerization, the isomerization of the Cyt2t⁺ to CytN3⁺ isomer remains unfavorable, and the HCOOH and HCOOHH₂O group mediated mechanisms are still more favorable. Meanwhile, conventional transition state theory (CTST) followed by Wigner tunneling correction is then applied to estimate the rate constants. The rate constant with Wigner tunneling correction for direct tautomerization is obviously smaller than that of HCOOH-mediated tautomerization, which is the most plausible mechanism. Finally, another important finding is that the product complex (CytN3⁺HCOOH) is in the rapid tautomeric equilibrium with the reaction complex (Cyt2t⁺HCOOH) (τ_99.9% = 3.84 × 10^-12 s), which is implemented by the mechanism of the concerted synchronous double proton transfer. Its lifetime of the formed CytN3⁺HCOOH complex (τ = 8.33 × 10^-9 s) is almost one order of magnitude larger than the time required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication (several ns), which is further dissociated into the CytN3⁺ and HCOOH monomers. The results of the present study demonstrate the feasibility of acid catalysis for DNA base isomerization reactions that would otherwise be forbidden.

The A·T(rWC)/A·T(H)/A·T(rH) ↔ A·T*(rwWC)/A·T*(wH)/A·T*(rwH) mutagenic tautomerization via sequential proton transfer: a QM/QTAIM study

In this study for the first time we have revealed by QM and QTAIM calculations at the MP2/aug-cc-pVDZ//B3LYP/6-311++G(d,p) level of QM theory the novel routes of the mutagenic tautomerization of three biologically important A·T DNA base pairs – reverse Watson–Crick A·T(rWC), Hoogsteen A·T(H) and reverse Hoogsteen A·T(rH) – followed by their rebuilding into the wobble (w) A·T*(rw_WC), A·T*(w_H) and A·T*(rw_H) base mispairs by the participation of the mutagenic tautomers of the DNA bases (denoted by asterisk) and vice versa, thus complementing the physico-chemical property of the canonical A·T(WC) Watson–Crick DNA base pair reported earlier (Brovarets' et al., RSC Adv., 2015, 5, 99594–99605). These non-dissociative tautomeric transformations in the classical A·T(rWC), A·T(H) and A·T(rH) DNA base pairs proceed similarly to the canonical A·T(WC) DNA base pair via the intrapair sequential proton transfer with shifting towards major or minor grooves of DNA followed by further double proton transfer along the intermolecular H-bonds and are controlled by the plane symmetric and highly stable transition states – tight ion pairs formed by the A⁺ nucleobase, protonated by the N1/N7 nitrogen atoms, and T⁻ nucleobase, deprotonated by the N3H imino group. Comparison of the estimated populations of the tautomerised states (10⁻²¹ to 10⁻¹⁴) with similar characteristics for the canonical A·T(WC) DNA base pair (10⁻⁸ to 10⁻⁷) leads authors to the conclusion, that only a base pair with WC architecture can be a building block of the DNA macromolecule as a genetic material, which is able for the evolutionary self-development. Among all four classical DNA base pairs, only A·T(WC) DNA base pair can ensure the proper rate of the spontaneous point errors of replication in DNA.

We discovered tautomeric wobbling of the classical A·T DNA base pairs. This data evidence, that only a base pair with Watson–Crick architecture can be a building block of the DNA macromolecule as a genetic material, which is able for the evolutionary self-development.

Molecular dynamic simulation and DFT study on the Drug-DNA interaction; Crocetin as an anti-cancer and DNA nanostructure model

In this research, the interaction of Crocetin as an anti-cancer drug and a Dickerson DNA has been investigated. 25 ns molecular dynamic simulations of Crocetin and DNA composed of 12 base pairs and a sequence of d(CGCGAATTCGCG)₂ were done in water. Three definite parts of the B-DNA were considered in analyzing the best interactive site from the thermodynamic point of view. Binding energy analysis showed that van der Waals interaction is the most important part related to the reciprocal O and H atoms of the Crocetin and DNA. Stabilizing interactions, obtained by ΔG calculations, showed that maximum and minimum interactions are related to the S1 and S3 regions, respectively. This means that the most probable van der Waals interaction site of the Dickerson B-DNA and Crocetin is located in the minor groove of DNA. Two sharp peaks at 2.55 and 1.75 Å in radial distribution functions of the PO⋯HO and NH⋯OC parts are related to new hydrogen bonds between the Crocetin and DNA in the complex which can be considered as the driving force of the anti-cancer mechanism of the Crocetin. Average values of 0.3 au and zero for the electron densities of the H⋯O bonds for DNA and complex, obtained by Quantum theory of atoms in molecules (QTAIM), means that the origin of DNA instability after complexation may be related to the H-bond denaturation by Crocetin. Finally, the evaluation of the dispersion interactions using the dispersion functional, -148.76 kcal.mol^-1, confirmed the importance of the dispersion interaction in drug-DNA complex.

Tautomeric transition between wobble A·C DNA base mispair and Watson-Crick-like A·C* mismatch: microstructural mechanism and biological significance

Here, we use MP2/DFT quantum-chemical methods combined with Quantum Theory of Atoms in Molecules to study the tautomeric transition between wobble A·C(w) mismatch and Watson-Crick-like A·C*(WC) base mispair, proceeding non-dissociatively via sequential proton transfer between bases through the planar, highly stable and zwitterionic TS(A∙C-)(A∙C(W)<-->A∙C&(WC)) transition state joined by the participation of (A)N6(+)H∙∙∙N4(-)(C), (A)N1(+)H∙∙∙N4(-)(C) and (A)C2(+)H∙∙∙N3(-)(C) H-bonds. Notably, the A·C(w) ↔ A·C*(WC) tautomerization reaction is accompanied by 10 unique patterns of the specific intermolecular interactions that consistently replace each other. Our data suggest that biologically significant A·C(w) → A·C*(WC) tautomerization is a kinetically controlled pathway for formation of the enzymatically competent Watson-Crick-like A·C*(WC) DNA base mispair in the essentially hydrophobic recognition pocket of the high-fidelity DNA-polymerase, responsible for the occurrence of spontaneous point AC/CA incorporation errors during DNA biosynthesis.

In silico studies with substituted adenines to achieve a remarkable stability of mispairs with thymine nucleobase

The modified adenine and thymine mispair achieves a remarkable stability, which can presumably help the DNA lesions to be less cytotoxic.

Wobble↔Watson-Crick tautomeric transitions in the homo-purine DNA mismatches: a key to the intimate mechanisms of the spontaneous transversions

The intrinsic capability of the homo-purine DNA base mispairs to perform wobble↔Watson-Crick/Topal-Fresco tautomeric transitions via the sequential intrapair double proton transfer was discovered for the first time using QM (MP2/DFT) and QTAIM methodologies that are crucial for understanding the microstructural mechanisms of the spontaneous transversions.

Proton tunneling in the A∙T Watson-Crick DNA base pair: myth or reality?

The results and conclusions reached by Godbeer et al. in their recent work, that proton tunneling in the A∙T(WC) Watson-Crick (WC) DNA base pair occurs according to the Löwdin's (L) model, but with a small (~10(-9)) probability were critically analyzed. Here, it was shown that this finding overestimates the possibility of the proton tunneling at the A∙T(WC)↔A*∙T*(L) tautomerization, because this process cannot be implemented as a chemical reaction. Furthermore, it was outlined those biologically important nucleobase mispairs (A∙A*↔A*∙A, G∙G*↔G*∙G, T∙T*↔T*∙T, C∙C*↔C*∙C, H∙H*↔H*∙H (H - hypoxanthine)) - the players in the field of the spontaneous point mutagenesis - where the tunneling of protons is expected and for which the application of the model proposed by Godbeer et al. can be productive.

Structure and conformational analysis of the anti-HIV reverse transcriptase inhibitor AZT using MP2 and DFT methods. Differences with the natural nucleoside thymidine. Simulation of the 1st phosphorylation step with ATP

A comprehensive quantum-chemical investigation of the conformational landscape of the HIV-1 reverse transcriptase inhibitor AZT (3'-azido-3'-deoxythymidine) nucleoside analogue was carried out. The whole conformational parameters (χ, γ, β, δ, ϕ, P, νmax) were analysed as well as the NBO charges. The search located at least 55 stable structures, 9 of which were by MP2 within a 1 kcal mol(-1) electronic energy range of the global minimum. Most conformers were anti or high-anti around the glycoside bond and with North sugar ring puckering angles. The distribution of all the conformers according to the ranges of stability of the characteristic torsional angles was established. The results obtained were in accordance with those found in related anti-HIV nucleoside analogues. The best conformer in the anti form corresponded to the calculated values by MP2 of χ = -126.9°, β = 176.4° and γ = 49.1°. An analysis of the lowest vibrations in conformer C1 was carried out. The first hydration shell was simulated and the structural differences with the natural nucleoside deoxythymidine (dT) were determined. The first phosphorylation step was simulated by interacting ATP with the best hydrated clusters of AZT and dT. The Na cations act as a bridge between the phosphate moieties of ATP making it easy for -P3O3 to receive the H5' proton from AZT or dT. A proton-transfer mechanism is proposed through the water molecules. When the number of the water molecules surrounding AZT is lower than 8, the first phosphorylation step of AZT can be carried out. However, the appropriate orientation of the O5'-H in dT avoids this limitation and it can be performed with large numbers of water molecules.

How many tautomerization pathways connect Watson-Crick-like G*·T DNA base mispair and wobble mismatches?

In this study, we have theoretically demonstrated the intrinsic ability of the wobble G·T(w)/G*·T*(w)/G·T(w1)/G·T(w2) and Watson-Crick-like G*·T(WC) DNA base mispairs to interconvert into each other via the DPT tautomerization. We have established that among all these transitions, only one single G·T(w) ↔ G*·T(WC) pathway is eligible from a biological perspective. It involves short-lived intermediate - the G·T*(WC) base mispair - and is governed by the planar, highly stable, and zwitterionic [Formula: see text] transition state stabilized by the participation of the unique pattern of the five intermolecular O6(+)H⋯O4(-), O6(+)H⋯N3(-), N1(+)H⋯N3(-), N1(+)H⋯O2(-), and N2(+)H⋯O2(-) H-bonds. This non-dissociative G·T(w) ↔ G*·T(WC) tautomerization occurs without opening of the pair: Bases within mispair remain connected by 14 different patterns of the specific intermolecular interactions that successively change each other along the IRC. Novel kinetically controlled mechanism of the thermodynamically non-equilibrium spontaneous point GT/TG incorporation errors has been suggested. The mutagenic effect of the analogues of the nucleotide bases, in particular 5-bromouracil, can be attributed to the decreasing of the barrier of the acquisition by the wobble pair containing these compounds of the enzymatically competent Watson-Crick's geometry via the intrapair mutagenic tautomerization directly in the essentially hydrophobic recognition pocket of the replication DNA-polymerase machinery. Proposed approaches are able to explain experimental data, namely growth of the rate of the spontaneous point incorporation errors during DNA biosynthesis with increasing temperature.

How does the long G·G* Watson-Crick DNA base mispair comprising keto and enol tautomers of the guanine tautomerise? The results of a QM/QTAIM investigation

The double proton transfer (DPT) in the long G·G* Watson-Crick base mispair (|C6N1(G*)N1C6(G)| = 36.4°; C1 symmetry), involving keto and enol tautomers of the guanine (G) nucleobase, along two intermolecular neighboring O6H···O6 (8.39) and N1···HN1 (6.14 kcal mol(-1)) H-bonds that were established to be slightly anti-cooperative, leads to its transformation into the G*·G base mispair through a single transition state (|C6N1N1C6| = 37.1°; C1), namely to the interconversion into itself. It was shown that the G·G* ↔ G*·G tautomerisation via the DPT is assisted by the third specific contact, that sequentially switches along the intrinsic reaction coordinate (IRC) in an original way: (G)N2H···N2(G*) H-bond (-25.13 to -10.37) → N2···N2 van der Waals contact (-10.37 to -9.23) → (G)N2···HN2(G*) H-bond (-9.23 to 0.79) → (G*)N2···HN2(G) H-bond (0.79 to 7.35 Bohr). The DPT tautomerisation was found to proceed through the asynchronous concerted mechanism by employing the QM/QTAIM approach and the methodology of the scans of the geometric, electron-topological, energetic, polar and NBO properties along the IRC. Nine key points, that can be considered as part of the tautomerisation repertoire, have been established and analyzed in detail. Furthermore, it was shown that the G·G* or G*·G base mispair is a thermodynamically and dynamically stable structure with a lifetime of 8.22 × 10(-10) s and all 6 low-frequency intermolecular vibrations are able to develop during this time span. Lastly, our results highlight the importance of the G·G* ↔ G*·G DPT tautomerisation, which can have implications for biological and chemical sensing applications.

Molecular docking and dynamics simulations on the interaction of cationic porphyrin-anthraquinone hybrids with DNA G-quadruplexes

A series of cationic porphyrin-anthraquinone hybrids bearing either pyridine, imidazole, or pyrazole rings at the meso-positions have been investigated for their interaction with DNA G-quadruplexes by employing molecular docking and molecular dynamics simulations. Three types of DNA G-quadruplexes were utilized, which comprise parallel, antiparallel, and mixed hybrid topologies. The porphyrin hybrids have a preference to bind with parallel and mixed hybrid structures compared to the antiparallel structure. This preference arises from the end stacking of porphyrin moiety following G-stem and loop binding of anthraquinone tail, which is not found in the antiparallel due to the presence of diagonal and lateral loops that crowd the G-quartet. The binding to the antiparallel, instead, occurred with poorer affinity through both the loop and wide groove. All sites of porphyrin binding were confirmed by 6 ns molecular dynamics simulation, as well as by the negative value of the total binding free energies that were calculated using the MMPBSA method. Free energy analysis shows that the favorable contribution came from the electrostatic term, which supposedly originated from the interaction of either cationic pyridinium, pyrazole, or imidazole groups and the anionic phosphate backbone, and also from the van der Waals energy, which primarily contributed through end stacking interaction.

DPT tautomerisation of the G·A(syn) and A*·G*(syn) DNA mismatches: a QM/QTAIM combined atomistic investigation

By applying a combined QM and QTAIM atomistic computational approach we have established for the first time that the G·A(syn) and A*·G*(syn) DNA mismatches (rare tautomers are marked with an asterisk), causing spontaneous transversions with substantially various probabilities, radically differ from each other in their ability to tautomerise through the double proton transfer (DPT). The A*·G*(syn) mismatch tautomerises quite easily (ΔΔG(TS) ≈ 4·kT at room temperature) into the A·G*(syn) mismatch through the asynchronous concerted mechanism, whereas the G·A(syn) base mispair does not tautomerise via the DPT at all, since there is no local minimum corresponding to the tautomerised G*·A*(syn) mismatch on the potential energy surface. It was established that the A·G*(syn) base mispair is a dynamically unstable H-bonded complex with an extremely short lifetime of 2.17 × 10(-13) s. Consequently, the obtained results allow us to believe that spontaneous or forced dissociation of both the G·A(syn) and A*·G*(syn) DNA mismatches by the DNA-polymerase occurs with the preservation of the tautomeric status of the bases.

Theoretical insight into the structural mechanism for the binding of vinblastine with tubulin

Vinblastine (VLB) is one of vinca alkaloids with high cytotoxicity toward cancer cells approved for clinical use. However, because of drug resistance, toxicity, and other side effects caused from the use of VLB, new vinca alkaloids with higher cytotoxicity toward cancer cells and other good qualities need to develop. One strategy is to further study and better understand the essence why VLB possesses the high cytotoxicity toward cancer cells. In present work, by using molecular simulation, molecular docking, density functional calculation, and the crystal structure of α,β-tubulin complex, we find two modes labeled in catharanthine moiety (CM) and vindoline moiety (VM) modes of VLB bound with the interface of α,β-tubulin to probe the essence why VLB has the high cytotoxicity toward cancer cells. In the CM mode, nine key residues B-Ser178, B-Asp179, B-Glu183, B-Tyr210, B-Asp226, C-Lys326, C-Asp327, C-Lys336, and C-Lys352 from the α,β-tubulin complex are determined as the active sites for the interaction of VLB with α,β-tubulin. Some of them such as B-Ser178, B-Glu183, B-Tyr210, B-Asp226, C-Lys326, C-Asp327, and C-Lys336 are newly identified as the active sites in present work. The affinity between VLB and the active pocket within the interface of α,β-tubulin is -60.8 kJ mol(-1) in the CM mode. In the VM mode, that is a new mode established in present paper, nine similar key residues B-Lys176, B-Ser178, B-Asp179, B-Glu183, B-Tyr210, B-Asp226, C-Lys326, C-Asp327, and C-Lys336 from the α,β-tubulin complex are found as the active sites for the interaction with VLB. The difference is from one key residue C-Lys352 in the CM mode changed to the key residue B-Lys176 in the VM mode. The affinity between VLB and the active pocket within the interface of α,β-tubulin is -96.3 kJ mol(-1) in the VM mode. Based on the results obtained in present work, and because VLB looks like two faces, composed of CM and VM both to have similar polar active groups, to interact with the active sites, we suggest double-faces sticking mechanism for the binding of VLB to the interface of α,β-tubulin. The double-faces sticking mechanism can be used to qualitatively explain high cytotoxicity toward cancer cells of vinca alkaloids including vinblastine, vincristine, vindestine, and vinorelbine approved for clinical use and vinflunine still in a phase III clinical trial. Furthermore, this mechanism will be applied to develop novel vinca alkaloids with much higher cytotoxicity toward cancer cells.

Genetic algorithms coupled with quantum mechanics for refinement of force fields for RNA simulation: a case study of glycosidic torsions in the canonical ribonucleosides

We report the use of genetic algorithms (GA) as a method to refine force field parameters in order to determine RNA energy.

New structural hypostases of the A·T and G·C Watson–Crick DNA base pairs caused by their mutagenic tautomerisation in a wobble manner: a QM/QTAIM prediction

Our investigation reveals the hitherto unknown ability of the canonical Watson–Crick DNA base pairs to switch into wobble mismatches with mutagenic tautomers, clarifying the nature of genome instability.

FT-IR and FT-Raman spectra of 6-chlorouracil: molecular structure, tautomerism and solid state simulation. A comparison between 5-chlorouracil and 6-chlorouracil

A Raman and IR study of the biomolecule 6-chlorouracil was carried out in the solid state. The unit cell found in the crystal was simulated as a tetramer form by density functional calculations. Specific scale factors and scaling equations deduced from uracil molecule were employed in the predicted wavenumbers of 6-chlorouracil. The scaled wavenumbers were used in the reassignment of the IR and Raman experimental bands. Good reproduction of the experimental wavenumbers is obtained and the % error is very small in the majority of cases. A comparison between the molecular structure and charge distribution of 6-chlorouracil and 5-chlorouracil molecules was presented. The effect of the hydration with the PCM model in the molecular structure and charges was discussed. The optimum tautomers of 6-chlorouracil were optimized and analyzed. Six of them were related to those of uracil molecule. The effect of the halogen substitution in the sixth position of the pyrimidine ring in the stability of the different tautomers was evaluated. HOMO and LUMO orbital energy analysis were carried out.

The nature of the transition mismatches with Watson-Crick architecture: the G*·T or G·T* DNA base mispair or both? A QM/QTAIM perspective for the biological problem

This study provides the first accurate investigation of the tautomerization of the biologically important guanine*·thymine (G*·T) DNA base mispair with Watson-Crick geometry, involving the enol mutagenic tautomer of the G and the keto tautomer of the T, into the G·T* mispair (∆G = .99 kcal mol(-1), population = 15.8% obtained at the MP2 level of quantum-mechanical theory in the continuum with ε = 4), formed by the keto tautomer of the G and the enol mutagenic tautomer of the T base, using DFT and MP2 methods in vacuum and in the weakly polar medium (ε = 4), characteristic for the hydrophobic interfaces of specific protein-nucleic acid interactions. We were first able to show that the G*·T↔G·T* tautomerization occurs through the asynchronous concerted double proton transfer along two antiparallel O6H···O4 and N1···HN3 H-bonds and is assisted by the third N2H···O2 H-bond, that exists along the entire reaction pathway. The obtained results indicate that the G·T* base mispair is stable from the thermodynamic point of view complex, while it is dynamically unstable structure in vacuum and dynamically stable structure in the continuum with ε = 4 with lifetime of 6.4·10(-12) s, that, on the one side, makes it possible to develop all six low-frequency intermolecular vibrations, but, on the other side, it is by three orders less than the time (several ns) required for the replication machinery to forcibly dissociate a base pair into the monomers during DNA replication. One of the more significant findings to emerge from this study is that the short-lived G·T* base mispair, which electronic interaction energy between the bases (-23.76 kcal mol(-1)) exceeds the analogical value for the G·C Watson-Crick nucleobase pair (-20.38 kcal mol(-1)), "escapes from the hands" of the DNA replication machinery by fast transforming into the G*·T mismatch playing an indirect role of its supplier during the DNA replication. So, exactly the G*·T mismatch was established to play the crucial role in the spontaneous point mutagenesis.

Atomistic nature of the DPT tautomerisation of the biologically important C·C* DNA base mispair containing amino and imino tautomers of cytosine: a QM and QTAIM approach

A theoretical study of tautomerisation of the biologically important cytosine·cytosine* (C·C*) DNA mismatch with a propeller-like structure (|C4N3N3C4| = 32.4°; C1 symmetry) and cis-oriented N1H glycosidic bonds, formed by the amino and imino tautomers of the C nucleobase, via the asynchronous concerted double proton transfer (DPT) along two H-bonds through the transition state (TSC·C*↔C*·C) (|C4N3N3C4| = 48.5°; C1 symmetry) into the C*·C mispair was carried out for the first time. It was established that the C·C*/C*·C DNA base mispair is associated by the antiparallel N4H···N4 (6.66 kcal mol(-1)), N3H···N3 (6.47 kcal mol(-1)) H-bonds and the O2···O2 van der Waals (vdW) contact (0.33 kcal mol(-1)), while the zwitterionic TSC·C*↔C*·C is stabilized by the parallel N4(+)H···N4(-) (13.55 kcal mol(-1)), N3(+)H···N3(-) (13.20 kcal mol(-1)) H-bonds and the O2(+)···O2(-) vdW contact (0.60 kcal mol(-1)). It was shown that the C·C* ↔ C*·C tautomerisation via the DPT is assisted by the O2···O2 vdW contact, that in contrast to the two others N4H···N4 and N3H···N3 H-bonds exists along the entire intrinsic reaction coordinate (IRC) range. The positive values of the Grunenberg's compliance constants (30.919 and 21.384 Å mdyn(-1) for C·C*/C*·C and TSC·C*↔C*·C, respectively) indicate that the O2···O2 vdW contact is a stabilizing closed-shell interaction. It was found that the middle N3H···N3 H-bond is anti-cooperative with the upper N4H···N4 H-bond and cooperative with the lower O2···O2 vdW contact. The 9 key points, which can be considered as electron-topological "fingerprints" of the asynchronous concerted C·C* ↔ C*·C tautomerisation process via the DPT were revealed along the IRC and examined in detail. It was shown that the C·C*/C*·C base mispair is a thermodynamically and dynamically stable structure. Its lifetime is equal to 1.53 × 10(-7) s at the MP2/cc-pVQZ//B3LYP/6-311++G(d,p) level of theory in vacuum. All 6 low-frequency intermolecular vibrations are able to develop during this time span.

Conformational analysis, molecular structure and solid state simulation of the antiviral drug acyclovir (zovirax) using density functional theory methods

The five tautomers of the drug acyclovir (ACV) were determined and optimised at the MP2 and B3LYP quantum chemical levels of theory. The stability of the tautomers was correlated with different parameters. On the most stable tautomer N1 was carried out a comprehensive conformational analysis, and the whole conformational parameters (R, β, Φ, φ1, φ2, φ3, φ4, φ5) were studied as well as the NBO Natural atomic charges. The calculations were carried out with full relaxation of all geometrical parameters. The search located at least 78 stable structures within 8.5 kcal/mol electronic energy range of the global minimum, and classified in two groups according to the positive or negative value of the torsional angle j1. In the nitrogen atoms and in the O2' and O5' oxygen atoms of the most stable conformer appear a higher reactivity than in the natural nucleoside deoxyguanosine. The solid state was simulated through a dimer and tetramer forms and the structural parameters were compared with the X-ray crystal data available. Several general conclusions were emphasized.

DPT tautomerisation of the wobble guanine·thymine DNA base mispair is not mutagenic: QM and QTAIM arguments

We have shown for the first time, connecting QM methods with QTAIM analysis and using the methodology of the sweeps of the energetical, electron-topological and geometrical parameters, that the tautomerisation of the wobble guanine·thymine (wG·T) DNA base mispair into the wG(*)·T(*) base mispair induced by the double proton transfer (DPT), which undergoes a concerted asynchronous pathway, is not mutagenic. The wG·T → wG(*)·T(*) DPT tautomerisation does not result in the transition of the G base into its mutagenic tautomeric form G(*) able to mispair with the T base within the Watson-Crick base pairing scheme. This observation is explained by the so-called quantum protection of the wG·T DNA base mispair from its mutagenic tautomerisation - the dynamical non-stability of the tautomerised wG(*)·T(*) base mispair and significantly negative value of the Gibbs free energy of activation for the reverse reaction of the wG·T → wG(*)·T(*) DPT tautomerisation.

The binding modes of cationic porphyrin-anthraquinone hybrids to DNA duplexes: in silico study

Cationic porphyrin-anthraquinone hybrids bearing peripheral substituents, either pyridine, imidazole, or pyrazole rings have been investigated for their binding mode to DNA duplexes. The four kinds of DNA duplexes were used, which represent intercalation and groove binding modes. AutoDock 4.2 was used to dock nine hybrid compounds to four DNA duplexes, while monitoring of conformational changes of four best hybrid-DNA complexes during 2 ns was performed by Amber9 molecular dynamics package. The binding energy calculation of best four complexes was then carried out using MMPBSA method. The hybrid compounds interacted to DNA duplexes through intercalation and groove binding modes. The minor groove binding of DNA was energetically preferred by cationic porphyrin hybrids due to favorable electrostatic and van der Waals interactions. Both electrostatic and van der Waals contributions were able to distinguish the binding mode of porphyrin hybrid to DNA duplexes.

[Structural and energetic properties of the four configurations of the A.T and G.C DNA base pairs]

Using the methods of non-empirical quantum chemistry at the MP2/6-311++G(2df,pd)// B3LYP/6-311++G(d,p) level of theory it was established for the first time, that Hoogsteen, reverse Hoogsteen, Watson-Crick and reverse Watson-Crick configurations of the A.T and G.C DNA base pairs are isoelectronic and isomorphic structures with similar dynamic properties. Based on these results, non-ionisation mechanism of the Hoogsteen <"breathing" of the G*.C* DNA base pair, namely transformation of the tautomerised (Lowdin's) G-C base pair with Watson-Crick geometry into the Hoogsteen electroneutral G*.C* H base pair stabilized by the three O6H...N4, N3H...N7 and C8H...02 H-bonds, was postulated. It is suggested that such scenario activates only in those cases, when DNA is not located in aqueous solution, but works together with proteins and cytosine protonation at the N3 atom is precluded.

Is the DPT tautomerization of the long A·G Watson-Crick DNA base mispair a source of the adenine and guanine mutagenic tautomers? A QM and QTAIM response to the biologically important question

Herein, we first address the question posed in the title by establishing the tautomerization trajectory via the double proton transfer of the adenine·guanine (A·G) DNA base mispair formed by the canonical tautomers of the A and G bases into the A*·G* DNA base mispair, involving mutagenic tautomers, with the use of the quantum-mechanical calculations and quantum theory of atoms in molecules (QTAIM). It was detected that the A·G ↔ A*·G* tautomerization proceeds through the asynchronous concerted mechanism. It was revealed that the A·G base mispair is stabilized by the N6H···O6 (5.68) and N1H···N1 (6.51) hydrogen bonds (H-bonds) and the N2H···HC2 dihydrogen bond (DH-bond) (0.68 kcal·mol(-1) ), whereas the A*·G* base mispair-by the O6H···N6 (10.88), N1H···N1 (7.01) and C2H···N2 H-bonds (0.42 kcal·mol(-1) ). The N2H···HC2 DH-bond smoothly and without bifurcation transforms into the C2H···N2 H-bond at the IRC = -10.07 Bohr in the course of the A·G ↔ A*·G* tautomerization. Using the sweeps of the energies of the intermolecular H-bonds, it was observed that the N6H···O6 H-bond is anticooperative to the two others-N1H···N1 and N2H···HC2 in the A·G base mispair, while the latters are significantly cooperative, mutually strengthening each other. In opposite, all three O6H···N6, N1H···N1, and C2H···N2 H-bonds are cooperative in the A*·G* base mispair. All in all, we established the dynamical instability of the А*·G* base mispair with a short lifetime (4.83·10(-14) s), enabling it not to be deemed feasible source of the A* and G* mutagenic tautomers of the DNA bases. The small lifetime of the А*·G* base mispair is predetermined by the negative value of the Gibbs free energy for the A*·G* → A·G transition. Moreover, all of the six low-frequency intermolecular vibrations cannot develop during this lifetime that additionally confirms the aforementioned results. Thus, the A*·G* base mispair cannot be considered as a source of the mutagenic tautomers of the DNA bases, as the A·G base mispair dissociates during DNA replication exceptionally into the A and G monomers in the canonical tautomeric form.

A theoretical investigation on the conformation and the interaction of CHF₂OCF₂CHF₂ (desflurane II) with one water molecule

The conformation and the interaction of CHF₂OCF₂CHF₂ (desflurane II) with one water molecule is investigated theoretically using the ab initio MP2/aug-cc-pvdz and DFT-based M062X/6-311++G(d,p) methods. The calculations include the optimized geometries, the harmonic frequencies of relevant vibrational modes along with a natural bond orbital (NBO) analysis including the NBO charges, the hybridization of the C atom and the intra- and intermolecular hyperconjugation energies. In the two most stable conformers, the CH bond of the F2HCO- group occupies the gauche position. The hyperconjugation energies are about the same for both conformers and the conformational preference depends on the interaction between the non-bonded F and H atoms. The deprotonation enthalpies of the CH bonds are about the same for both conformers, the proton affinity of the less stable conformer being 3 kcal mol−1 higher. Both conformers of desflurane II interact with water forming cyclic complexes characterized by CH…O and OH…F hydrogen bonds. The binding energies are moderate, ranging from −2.4 to −3.2 kcal mol−1 at the MP2 level. The origin of the blue shifts of the ν(CH) vibrations is analyzed. In three of the complexes, the water molecule acts as an electron donor. Interestingly, in these cases a charge transfer is also directed to the non bonded OH group of the water molecule. This effect seems to be a property of polyfluorinated ethers.

Configuration of wobble base pairs having pyrimidines as anticodon wobble bases: significance for codon degeneracy

Degeneracy of the genetic code was attributed by Crick to imprecise hydrogen-bonded base-pairing at the wobble position during codon-anticodon pairing. The Crick wobble rules define but do not explain the RNA base pair combinations allowed at this position. We select six pyrimidine bases functioning as anticodon wobble bases (AWBs) to study their H-bonded pairing properties with the four major RNA bases using density functional theory at the B3LYP/6-31G(d,p) level. This is done to assess the extent to which the configuration of a solitary RNA wobble base pair may in itself determine specificity and degeneracy of the genetic code by allowing or disallowing the given base pair during codon-anticodon pairing. Calculated values of select configuration markers for the base pairs screen well between allowed and disallowed base pairs for most cases examined here, where the base pair width emerges as an important factor. A few allowed wobble pairs invoke the involvement of RNA nucleoside conformation, as well as involvement of the exocyclic substituent in H-bonding. This study, however, cannot explain the disallowed status of the Ura⋯Gua wobble pair on the basis of configuration alone. Explanation of the allowed status of the V⋯Ura pair requires further study on the mediatory role of water molecules. Apart from these two cases, these computational results are sufficient, on the basis of base pair configuration alone, to account for the specificity and degeneracy of the genetic code for all known cases of codon-anticodon pairing which involve the pyrimidine AWBs studied here.