DFT study of the interaction of cytidine and 2′-deoxycytidine with Li+, Na+, and K+: effects of metal cationization on sugar puckering and stability of the N-glycosidic bond

First principles modelling of the ion binding capacity of finger millet

Finger millet, a cereal grain widely consumed in India and Africa, has gained more attention in recent years due to its high dietary fibre (arabinoxylan) and trace mineral content, and its climate resilience. The aim of this study was to understand the interactions between potassium (K+), calcium (Ca2+) and zinc (Zn2+) ions and the arabinoxylan structure and determine its ion-binding capacity. Three variations of a proposed model of the arabinoxylan structure were constructed and first principles Density Functional Theory calculations were carried out to determine the cation-binding capacity of the arabinoxylan complexes. Zn2+-arabinoxylan complexes were highly unstable and thermodynamically unfavourable in all three models. Ca2+ and K+ ions, however, form thermodynamically stable complexes, particularly involving two glucuronic acid residues as a binding pocket. Glucuronic acid residues are found to play a key role in stabilising the cation-arabinoxylan complex, and steric effects are more important to the stability than charge density. Our results highlight the most important structural features of the millet fibre regarding ion-storage capacity, and provide valuable preliminary data for confirmatory experimental studies and for the planning of clinical trials where the bioavailability of bound ions following digestion may be tested.

DFT investigation of the regioselective allylation of pyrimidine 2'-deoxynucleosides

To understand the regioselectivity observed in the allylation of pyrimidine nucleosides and to identify the factors directing the reaction, a theoretical study of the regioselective allylation was carried out. Several key points were considered such as: the structure of the deprotonated nucleobase in the presence of Na+; the effect of the solvent on the dissociation and aggregation reactions of thymidine/Na+ ion pair; and the likely allylation reaction mechanisms involved. The results showed that the regioselectivity observed experimentally can be attributed to a greater stability of a dimeric form coupled to an increase of the reaction barrier in THF due to larger Na+ binding to the nucleobase.

IRMPD Action Spectroscopy, ER-CID Experiments, and Theoretical Studies of Sodium Cationized Thymidine and 5-Methyluridine: Kinetic Trapping During the ESI Desolvation Process Preserves the Solution Structure of [Thd+Na]<sup/>

Thymidine (dThd) is a fundamental building block of DNA nucleic acids, whereas 5-methyluridine (Thd) is a common modified nucleoside found in tRNA. In order to determine the conformations of the sodium cationized thymine nucleosides [dThd+Na]⁺ and [Thd+Na]⁺ produced by electrospray ionization, their infrared multiple photon dissociation (IRMPD) action spectra are measured. Complementary electronic structure calculations are performed to determine the stable low-energy conformations of these complexes. Geometry optimizations and frequency analyses are performed at the B3LYP/6-311+G(d,p) level of theory, whereas energies are calculated at the B3LYP/6-311+G(2d,2p) level of theory. As protonation preferentially stabilizes minor tautomers of dThd and Thd, tautomerization facilitated by Na⁺ binding is also considered. Comparisons of the measured IRMPD and computed IR spectra find that [dThd+Na]⁺ prefers tridentate (O2,O4',O5') coordination to the canonical 2,4-diketo form of dThd with thymine in a syn orientation. In contrast, [Thd+Na]⁺ prefers bidentate (O2,O2') coordination to the canonical 2,4-diketo tautomer of Thd with thymine in an anti orientation. Although 2,4-dihydroxy tautomers and O2 protonated thymine nucleosides coexist in the gas phase, no evidence for minor tautomers is observed for the sodium cationized species. Consistent with experimental observations, the computational results confirm that the sodium cationized thymine nucleosides exhibit a strong preference for the canonical form of the thymine nucleobase. Survival yield analyses based on energy-resolved collision-induced dissociation (ER-CID) experiments suggest that the relative stabilities of protonated and sodium cationized dThd and Thd follow the order [dThd+H]⁺ < [Thd+H]⁺ < [dThd+Na]⁺ < [Thd+Na]⁺. Graphical Abstract ᅟ.

Pyrrolyl-, 2-(2-thienyl)pyrrolyl- and 2,5-bis(2-thienyl)pyrrolyl-nucleosides: synthesis, molecular and electronic structure, and redox behaviour of C5-thymidine derivatives

A series of modified nucleosides based on thymidine have been prepared by Pd-catalysed cross-coupling between N-alkyl-alkynyl functionalised pyrrolyl- (py), 2-(2-thienyl)pyrrolyl- (tp) or 2,5-bis(2-thienyl)pyrrolyl (tpt) groups with 5-iodo-2'-deoxyuridine. The length of the alkyl chain linking the nucleoside and pyrrolyl-containing unit, N(CH(2))(n)C[triple bond, length as m-dash]C-nucleoside (where n = 1-3) was also varied. The compounds have been characterised by (1)H NMR, ES-MS, UV-vis, cyclic voltammetry (CV) and, in some cases, single-crystal X-ray diffraction. Cyclic voltammetry studies demonstrated that all the py-, tp- and tpt-alkynyl derivatives 1-7 can be electrochemically polymerised to form conductive materials. It was found that increasing the N-alkyl chain length in these cases resulted in only minor changes in the oxidation potential. The same behaviour was observed for the tp- and tpt-modified nucleosides 9-12; however, the py-derivative, 8, produced a poorly conducting material. DFT calculations on the one-electron oxidised cation of the modified nucleosides bearing tp or tpt showed that spin density is located on the pyrrolyl and thienyl units in all cases and that the coplanarity of adjacent rings increases upon oxidation. In contrast, in the corresponding pyrrolyl cases the spin density is distributed over the whole molecule, suggesting that polymerisation does not occur solely at the pyrrolyl-Cα position and the conjugation is interrupted.

Effects of nucleophile, oxidative damage, and nucleobase orientation on the glycosidic bond cleavage in deoxyguanosine

Deglycosylation of nucleotides occurs during many essential biological processes, including DNA repair, and is initiated by a variety of nucleophiles. In the present work, density functional theory (B3LYP) was used to investigate the thermodynamics and kinetics of the glycosidic bond cleavage reaction in the model nucleoside forms of guanine and its major oxidation product, 8-oxoguanine. Base excision facilitated by four different nucleophiles (hydroxyl anion (fully activated water), formate-water complex (partially activated water), lysine, and proline) was considered, which spans nucleophiles involved in a collection of spontaneous and enzyme-catalyzed processes. Because some enzymes that catalyze deglycosylation can accommodate more than one orientation of the base with respect to the sugar moiety, the effects of the (anti/syn) base orientation on the barrier height were also considered. We find that the nucleophile has a very large effect on the overall (gas-phase) reaction energetics. Although this effect decreases in different (polar) environments, the nucleophile has the greatest influence on the overall reaction as compared to whether the base is damaged or to the base orientation. Furthermore, the effects are significant in environments that most closely resemble (nonpolar) enzymatic active sites. Our results provide a greater understanding of the relative effects of the nucleophile, damage to the nucleobase, and the nucleobase orientation with respect to the sugar moiety on the deglycosylation pathway, which provide qualitative explanations for relative base excision rates observed in some biological systems.

Modeling the dissociative hydrolysis of the natural DNA nucleosides

Two-dimensional PCM-B3LYP/6-31+G(d) potential energy surfaces for the hydrolysis of the four natural 2'-deoxyribonucleosides (2'-deoxyadenosine, 2'-deoxyguanosine, 2'-deoxycytidine, and thymidine) are characterized using a model that includes both implicit (bulk) solvent effects and (three or four) explicit water molecules in the optimization routine. For the first time, the experimentally predicted dissociative (S(N)1) mechanism is found to be favored over the synchronous (S(N)2) pathway for all nucleosides studied. Due to the success of our model in stabilizing the charge-separated intermediates along the S(N)1 pathway, it is proposed that the new model presented here is the smallest system capable of generating the experimentally predicted oxacarbenium cation intermediate. We therefore stress that dissociative mechanisms should be studied with methodologies that account for the (bulk) environment in the optimization routine, where these effects are often only included as a correction to the energy in the current literature. In addition to accounting for charge stabilization through implicit solvation, nucleophile activation and leaving group stabilization should also be explicitly introduced into the model to further stabilize the system. Our work also emphasizes the importance of studying the Gibbs surface, which in some cases provides a better description of chemically important regions of the reaction surface or changes the calculated trend in the magnitude of dissociative barriers. In addition, it is proposed that the methodology presented in this study can be used to calculate uncatalyzed deglycosylation barriers for a range of DNA nucleosides, which when compared to the corresponding enzyme-catalyzed reactions, will allow the prediction of the rate enhancement (barrier reduction) due to the enzyme.

Structures of human MST3 kinase in complex with adenine, ADP and Mn2+

The MST family is a subclass of mammalian serine/threonine kinases that are related to the yeast sterile-20 protein and are implicated in regulating cell growth and transformation. The MST3 protein contains a 300-residue catalytic domain and a 130-residue regulatory domain, which can be cleaved by caspase and activated by autophosphorylation, promoting apoptosis. Here, five crystal structures of the catalytic domain of MST3 are presented, including a complex with ADP and manganese, a unique cofactor preferred by the enzyme, and a complex with adenine. Similar to other protein kinases, the catalytic domain of MST3 folds into two lobes: the smaller N lobe forms the nucleotide-binding site and the larger C lobe recognizes the polypeptide substrate. The bound ADP and Mn(2+) ions are covered by a glycine-rich loop and held in place by Asn149 and Asp162. A different orientation was observed for the ligand in the MST3-adenine complex. In the activation loop, the side chain of Thr178 is phosphorylated and is sandwiched by Arg143 and Arg176. Comparison of this structure with other similar kinase structures shows a 180 degrees rotation of the loop, leading to activation of the enzyme. The well defined protein-ligand interactions also provide useful information for the design of potent inhibitors.