Tautomeric equilibrium, stability, and hydrogen bonding in 2'-deoxyguanosine monophosphate complexed with Mg2+

An insight to the molecular interactions of the FDA approved HIV PR drugs against L38L↑N↑L PR mutant

The aspartate protease of the human immune deficiency type-1 virus (HIV-1) has become a crucial antiviral target in which many useful antiretroviral inhibitors have been developed. However, it seems the emergence of new HIV-1 PR mutations enhances drug resistance, hence, the available FDA approved drugs show less activity towards the protease. A mutation and insertion designated L38L↑N↑L PR was recently reported from subtype of C-SA HIV-1. An integrated two-layered ONIOM (QM:MM) method was employed in this study to examine the binding affinities of the nine HIV PR inhibitors against this mutant. The computed binding free energies as well as experimental data revealed a reduced inhibitory activity towards the L38L↑N↑L PR in comparison with subtype C-SA HIV-1 PR. This observation suggests that the insertion and mutations significantly affect the binding affinities or characteristics of the HIV PIs and/or parent PR. The same trend for the computational binding free energies was observed for eight of the nine inhibitors with respect to the experimental binding free energies. The outcome of this study shows that ONIOM method can be used as a reliable computational approach to rationalize lead compounds against specific targets. The nature of the intermolecular interactions in terms of the host-guest hydrogen bond interactions is discussed using the atoms in molecules (AIM) analysis. Natural bond orbital analysis was also used to determine the extent of charge transfer between the QM region of the L38L↑N↑L PR enzyme and FDA approved drugs. AIM analysis showed that the interaction between the QM region of the L38L↑N↑L PR and FDA approved drugs are electrostatic dominant, the bond stability computed from the NBO analysis supports the results from the AIM application. Future studies will focus on the improvement of the computational model by considering explicit water molecules in the active pocket. We believe that this approach has the potential to provide information that will aid in the design of much improved HIV-1 PR antiviral drugs.

Protonation induces base rotation of purine nucleotides pdGuo and pGuo

Infrared multiple photon dissociation (IRMPD) action spectra of the protonated forms of 2'-deoxyguanosine-5'-monophosphate and guanosine-5'-monophosphate, [pdGuo+H](+) and [pGuo+H](+), are measured over the IR fingerprint and hydrogen-stretching regions using the FELIX free electron laser and an OPO/OPA laser system. Electronic structure calculations are performed to generate low-energy conformations of [pdGuo+H](+) and [pGuo+H](+) and determine their relative stabilities at the B3LYP/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) and MP2(full)/6-311+G(2d,2p)//B3LYP/6-311+G(d,p) levels of theory. Comparative analyses of the measured IRMPD action spectra and B3LYP/6-311+G(d,p) linear IR spectra computed for the low-energy conformers are performed to determine the most favorable site of protonation and the conformers present in the experiments. These comparisons and the computed energetics find that N7 protonation is considerably preferred over O6 and N3, and the N7 protonated ground-state conformers of [pdGuo+H](+) and [pGuo+H](+) are populated in the experiments. The 2'-hydroxyl substituent does not significantly impact the stable low-energy conformers of [pdGuo+H](+)vs. those of [pGuo+H](+). The effect of the 2'-hydroxyl substituent is primarily reflected in the relative intensities of the measured IRMPD bands, as the IRMPD profiles of [pdGuo+H](+) and [pGuo+H](+) are quite similar. Comparisons to previous IRMPD spectroscopy investigations of the protonated forms of the guanine nucleosides, [dGuo+H](+) and [Guo+H](+), and deprotonated forms of the guanine nucleotides, [pdGuo-H](-) and [pGuo-H](-), provide insight into the effects of the phosphate moiety and protonation on the conformational features of the nucleobase and sugar moieties. Protonation is found to induce base rotation of the guanine residue to an anti orientation vs. the syn orientation found for the deprotonated forms of the guanine nucleotides.

Molecular dynamics study on conformational differences between dGMP and 8-oxo-dGMP: Effects of metal ions

The modified nucleotide base 7,8-dihydro-8-oxo-guanine (8-oxo-G) is one of the major sources of spontaneous mutagenesis. Nucleotide-sanitizing enzymes, such as the MutT homolog-1 (MTH1) and nudix-type motif 5 (NUDT5), selectively remove 8-oxo-G from the cellular pool of nucleotides. Previous studies showed that, although the syn conformation generally predominates in purine nucleotides with a bulky substituent at the 8-position, 8-oxo-dGMP binds to both MTH1 and NUDT5 in the anti conformation. This study was initiated to investigate the possibility that 8-oxo-dGMP itself may adopt the anti conformation. Molecular dynamics simulations of mononucleotides (dGMP, 8-oxo-dGMP) in aqueous solution were performed. 8-oxo-dGMP adopted the anti conformation as well as the syn conformation, and the proportion of adopting the anti conformation increased in the presence of metal ions. When 8-oxo-dGMP was in the anti conformation, a metal ion was located between the oxygen atom of phosphate and the oxygen atom at the 8-position of 8-oxo-G. The types of stable anti conformations of 8-oxo-dGMP differed, depending on the ionic radii and charges of coexisting ions. These data suggested a role for metal ions, other than as cofactors for the hydrolysis of the di- and tri-phosphate forms of mononucleotides; that the metal ions help retain the anti conformation of the N-glycosidic torsion angle of 8-oxo-dGMP to promote the binding between the 8-oxo-G deoxynucleotide and the nucleotide-sanitizing enzymes.

B-DNA characteristics are preserved in double stranded d(A)3·d(T)3 and d(G)3·d(C)3 mini-helixes: conclusions from DFT/M06-2X study

We report the results of the first comprehensive DFT study on the d(A)3·d(T)3 and d(G)3·d(C)3 nucleic acid duplexes. The ability of mini-helixes to preserve the conformation of B-DNA in the gas phase and under the influence of such factors as: solvent, uncompensated charge, and counter-ions was evaluated using M06-2X functional with 6-31G(d,p) basis set. The accuracy of the models was ascertained based on their ability to reproduce key structural features of natural B-DNA. Analysis of the helicity suggests that the helical conformations adopt geometrical parameters which are close to those of the B-DNA form. The torsion angles fall somewhere between the values observed for BI/BII conformational classes. The comparative analysis of parameters of isolated Watson-Crick base pairs versus B-DNA-like conformations indicates the same tendency of base-pair polarization and hydration. Specifically, effects of polarization of nucleobases in continuum type dielectric medium mimicking water are stronger than those caused by the presence of backbone. Polar environment as well as the presence of counterions stabilizes duplexes, facilitating helix formation. Substantial conformational changes of nucleotides upon duplex formation decrease the binding energy. In spite of structural and energetic changes, the placement of a mini-helix into the gas phase does not lead to significant disruption of the structure. On the contrary, the duplex preserves its helicity and the strands remain bound.

Glycosidic bond cleavage in deoxynucleotides: effects of solvent and the DNA phosphate backbone in the computational model

Density functional theory (B3LYP) was employed to examine the hydrolysis of the canonical 2'-deoxynucleotides in varied environments (gas phase or water) using different computational models for the sugar residue (methyl or phosphate group at C5') and nucleophile (water activated through full or partial proton abstraction). Regardless of the degree of nucleophile activation, our results show that key geometrical parameters along the reaction pathway are notably altered upon direct inclusion of solvent effects in the optimization routine, which leads to significant changes in the reaction energetics and better agreement with experiment. Therefore, despite the wide use of gas-phase calculations in the literature, small model computational work, as well as large-scale enzyme models, that strive to understand nucleotide deglycosylation must adequately describe the environment. Alternatively, although inclusion of the phosphate group at C5' also affects the geometries of important stationary points, the effects cancel to yield unchanged deglycosylation barriers, and therefore smaller computational models can be used to estimate the energy associated with nucleotide deglycosylation, with the 5' phosphate group included if full (geometric) details of the reaction are desired. Hydrogen-bonding interactions with the nucleobase can significantly reduce the barrier to deglycosylation, which supports suggestions that discrete hydrogen-bonding interactions with active-site amino acid residues can play a significant role in enzyme-catalyzed nucleobase excision. Taken together with previous studies, the present work provides vital clues about the components that must be included in future studies of the deglycosylation of isolated noncanonical nucleotides, as well as the corresponding enzyme-catalyzed reactions.

Conformational flexibility of C8-phenoxylguanine adducts in deoxydinucleoside monophosphates

M06-2X/6-31G(d,p) is used to calculate the structure of all natural deoxydinucleoside monophosphates with G in the 5' or 3' position, the anti or syn conformation, and each natural (A, C, G, T) base in the corresponding flanking position. When the ortho or para C8-phenoxyl-2'-deoxyguanosine (C8-phenoxyl-dG) adduct replaces G in each model, there is little change in the relative base-base orientation or backbone conformation. However, the orientation of the C8-phenoxyl group can be characterized according to the position (5' versus 3'), conformation (anti versus syn), and isomer (ortho versus para) of damage. Although the degree of coplanarity between the phenoxyl ring and G base in the ortho adduct is highly affected by the sequence since the hydroxyl group can interact with neighboring bases, the para adduct generally does not exhibit discrete interactions with flanking bases. For both adducts, steric clashes between the phenoxyl group and the backbone or flanking base destabilize the anti conformation preferred by the natural nucleotide and thereby result in a clear preference for the syn conformation regardless of the sequence or position. This contrasts the conclusions drawn from smaller (nucleoside, nucleotide) models previously used in the literature, which stresses the importance of using models that address the steric constraints present due to the surrounding environment. Since replication errors for other C8-dG bulky adducts have been linked to a preference for the syn conformation, our findings provide insight into the possible mutagenicity of phenolic adducts.

Developing a computational model that accurately reproduces the structural features of a dinucleoside monophosphate unit within B-DNA

The ability of a dinucleoside monophosphate to mimic the conformation of B-DNA was tested using a combination of different phosphate models (anionic, neutral, counterion), environments (gas, water), and density functionals (B3LYP, MPWB1K, M06-2X) with the 6-31G(d,p) basis set. Three sequences (5'-GX(Py)-3', where X(Py) = T, U or (Br)U) were considered, which vary in the (natural or modified) 3' pyrimidine nucleobase (X(Py)). These bases were selected due to their presence in natural DNA, structural similarity to T and/or applications in radical-initiated anti-tumour therapies. The accuracy of each of the 54 model, method and sequence combinations was judged based on the ability to reproduce key structural features of natural B-DNA, including the stacked base-base orientation and important backbone torsion angles. B3LYP yields distorted or tilted relative base-base orientations, while many computational challenges were encountered for MPWB1K. Despite wide use in computational studies of DNA, the neutral (protonated) phosphate model could not consistently predict a stacked arrangement for all sequences. In contrast, stacked base-base arrangements were obtained for all sequences with M06-2X in conjunction with both the anionic and (sodium) counterion phosphate models. However, comparison of calculated and experimental backbone conformations reveals the charge-neutralized counterion phosphate model best mimics B-DNA. Structures optimized with implicit solvent (water) are comparable to gas-phase structures, which suggests similar results should be obtained in an environment of intermediate polarity. We recommend the use of either gas-phase or water M06-2X optimizations with the counterion phosphate model to study the structure and/or reactivity of natural or modified DNA with a dinucleoside monophosphate.

Effect of sodium and acetate ions on 8-hydroxyguanine formation in irradiated aqueous solutions of DNA and 2'-deoxyguanosine 5'-monophosphate

PURPOSE

The aim of this work was to study the combined effect of sodium and acetate ions on the radiation yield of 8-hydroxyguanine (8-OHG), one of the major DNA base lesions induced by free radicals.

MATERIALS AND METHODS

Aqueous solutions of DNA and 2'-deoxyguanosine 5'-monophosphate (dGMP) with various concentrations of sodium acetate and sodium perchlorate were γ-irradiated, enzymatically digested and analyzed by high-performance liquid chromatography (HPLC) methods.

RESULTS

It was found that both salts decrease the 8-OHG radiation yield in the concentration range studied for both DNA and dGMP, except in the case of dGMP wherein an increase in yield occurs in the concentration range from 0.1-1 mM. The dependence of the 8-hydroxy-2'-deoxyguanosine radiation yield on the concentration of both sodium acetate and sodium perchlorate have different shapes and have steeper slopes for the DNA compared with the dGMP solutions.

CONCLUSIONS

The observed decrease in the radiation yield of 8-OHG with increasing concentrations of sodium acetate is consistent with the hypothesis that sodium acetate produces two concentration-dependent effects in the DNA solutions: (1) A conformational change in the DNA caused by Na(+) counterions; and (2) free radical reactions related to the radiolysis of acetate ion.

Folding pathways of human telomeric type-1 and type-2 G-quadruplex structures

We have investigated new folding pathways of human telomeric type-1 and type-2 G-quadruplex conformations via intermediate hairpin and triplex structures. The stabilization energies calculated by ab initio methods evidenced the formation of a hairpin structure with Hoogsteen GG base pairs. Further calculations revealed that the G-triplet is more stable than the hairpin conformation and equally stable when compared to the G-tetrad. This indicated the possibility of a triplex intermediate. The overall folding is facilitated by K(+) association in each step, as it decreases the electrostatic repulsion. The K(+) binding site was identified by molecular dynamics simulations. We then focused on the syn/anti arrangement and found that the anti conformation of deoxyguanosine is more stable than the syn conformation, which indicated that folding would increase the number of anti conformations. The K(+) binding to a hairpin near the second lateral TTA loop was found to be preferable, considering entropic effects. Stacking of G-tetrads with the same conformation (anti/anti or syn/syn) is more stable than mixed stacking (anti/syn and vice versa). These results suggest the formation of type-1 and type-2 G-quadruplex structures with the possibility of hairpin and triplex intermediates.

Recognition and sensing of nucleoside monophosphates by a dicopper(II) cryptate

The dimetallic cryptate [Cu(2)(II)(1)](4+) selectively recognizes guanosine monophosphate with respect to other nucleoside monophosphates (NMPs) in a MeOH/water solution at pH 7. Recognition is efficiently signaled through the displacement of the indicator 6-carboxyfluorescein bound to the receptor, monitoring its yellow fluorescent emission. Titration experiments evidenced the occurrence of several simultaneous equilibria involving 1:1 and 2:1 receptor/NMP and receptor/indicator complexes. It was demonstrated that the added NMP displaces the indicator from the 2:1 receptor/indicator complex, forming the 1:1 receptor/analyte inclusion complex. Recognition selectivity is thus ascribed to the nature of nucleotide donor atoms involved in the coordination and their ability to encompass the Cu(II)-Cu(II) distance within the cryptate.

Electrically controlled DNA adhesion

The ability to control the interaction of polyelectrolytes, such as DNA or proteins, with charged surfaces is of pivotal importance for a multitude of biotechnological applications. Previously, we measured the desorption forces of single polymers on charged surfaces using an atomic force microscope. Here, we show that the adhesion of DNA on gold electrodes modified with self-assembled monolayers can be biased by the composition of the monolayer and externally controlled by means of the electrode potential. Positive potentials induced DNA adsorption onto OH-terminated electrodes with adhesion forces up to 25 pN (at +0.5 V versus Ag/AgCl), whereas negative potentials suppressed DNA adsorption. The measured contributions of the DNA backbone phosphate charges and the doubly charged terminal phosphate on adsorption agreed with a model based on the Gouy-Chapman theory. Experiments on an NH(2)-terminated electrode revealed a similar force modulation range of the coulomb component of the desorption force. These findings are important for the development of new DNA-based biochips or supramolecular structures.

Effect of a pH change on the conformational stability of the modified nucleotide queuosine monophosphate

The naturally occurring modified nucleotide queuosine 5'-monophosphate (QMP) related to biochemical regulatory pathways in the cell was investigated using quantum chemical approaches. The relative stability of biologically relevant conformations of QMP in solvent under a pH change was predicted at the BVP86/TZVP and MP2/TZVP levels of theory. Hydrogen bonding in QMP was studied using Bader's approach. The acidity constants of QMP were estimated using the COSMO-RS theory. It has been found that the neutral and anionic forms of QMP are the most stable in the physiological pH range. These forms correspond to the anti/north conformation and exist as zwitterionic tautomers having a negatively charged phosphate group (-1 for neutral and -2 for anionic) and a positively charged secondary amine group in the side chain. It was also found that QMP possesses the syn conformation in the cationic state at pH < 5.0 and undergoes syn to anti conformation transition when the pH increases, remaining in the anti conformation at the higher pH values. The marker IR bands specific for the anionic and neutral QMP forms in the 2300-2700 cm(-1) region were assigned to H-bonded NH groups of the QMP side chain. The bands between 800 and 1300 cm(-1) of the "fingerprint" (400-1500 cm(-1)) region were assigned to the vibrations of the ribose ring, the phosphate group and the side chain of QMP. The predicted IR spectra can be useful for the assignment of vibration bands in the experiential spectra of QMP or identification of the QMP forms. The revealed peculiarities of the QMP conformation sensitivity to a pH change as well as additional formed H-bonds could be responsible for specific nucleotide interactions with enzymes.

Effect of the ribose versus 2'-deoxyribose residue on the metal ion-binding properties of purine nucleotides

The interaction between metal ions and nucleotides is well characterized, as is their importance for metabolic processes, e.g. in the synthesis of nucleic acids. Hence, it is surprising to find that no detailed comparison is available of the metal ion-binding properties between nucleoside 5'-phosphates and 2'-deoxynucleoside 5'-phosphates. Therefore, we have measured here by potentiometric pH titrations the stabilities of several metal ion complexes formed with 2'-deoxyadenosine 5'-monophosphate (dAMP2-), 2'-deoxyadenosine 5'-diphosphate (dADP3-) and 2'-deoxyadenosine 5'-triphosphate (dATP4-). These results are compared with previous data measured under the same conditions and available in the literature for the adenosine 5'-phosphates, AMP(2-), ADP(3-) and ATP(4-), as well as guanosine 5'-monophosphate (GMP(2-)) and 2'-deoxyguanosine 5'-monophosphate (dGMP(2-)). Hence, in total four nucleotide pairs, GMP(2-)/dGMP(2-), AMP(2-)/dAMP(2-), ADP(3-)/dADP(3-) and ATP(4-)/dATP(4-) (= NP/dNP), could be compared for the four metal ions Mg2+, Ni2+, Cu2+ and Zn2+ (= M2+). The comparisons show that complex stability and extent of macrochelate formation between the phosphate-coordinated metal ion and N7 of the purine residue is very similar (or even identical) for the AMP(2-)/dAMP(2-) and ADP(3-)/dADP(3-) pairs. In the case of the complexes formed with ATP(4-)/dATP(4-) the 2'-deoxy complexes are somewhat more stable and show also a slightly enhanced tendency for macrochelate formation. This is different for guanine nucleotides: the stabilities of the M(dGMP) complexes are clearly higher, as are the formation degrees of their macrochelates, than is the case with the M(GMP) complexes. This enhanced complex stability and greater tendency to form macrochelates can be attributed to the enhanced basicity (DeltapKaca. 0.2) of N7 in the 2'-deoxy compound. These results allow general conclusions regarding nucleic acids to be made.