Amino Proton Exchange Processes in Mononucleosides

The catalytic mechanism of the RNA methyltransferase METTL3

The complex of methyltransferase-like proteins 3 and 14 (METTL3-14) is the major enzyme that deposits N6-methyladenosine (m6A) modifications on messenger RNA (mRNA) in humans. METTL3-14 plays key roles in various biological processes through its methyltransferase (MTase) activity. However, little is known about its substrate recognition and methyl transfer mechanism from its cofactor and methyl donor S-adenosylmethionine (SAM). Here, we study the MTase mechanism of METTL3-14 by a combined experimental and multiscale simulation approach using bisubstrate analogues (BAs), conjugates of a SAM-like moiety connected to the N6-atom of adenosine. Molecular dynamics simulations based on crystal structures of METTL3-14 with BAs suggest that the Y406 side chain of METTL3 is involved in the recruitment of adenosine and release of m6A. A crystal structure with a BA representing the transition state of methyl transfer shows a direct involvement of the METTL3 side chains E481 and K513 in adenosine binding which is supported by mutational analysis. Quantum mechanics/molecular mechanics (QM/MM) free energy calculations indicate that methyl transfer occurs without prior deprotonation of adenosine-N6. Furthermore, the QM/MM calculations provide further support for the role of electrostatic contributions of E481 and K513 to catalysis. The multidisciplinary approach used here sheds light on the (co)substrate binding mechanism, catalytic step, and (co)product release, and suggests that the latter step is rate-limiting for METTL3. The atomistic information on the substrate binding and methyl transfer reaction of METTL3 can be useful for understanding the mechanisms of other RNA MTases and for the design of transition state analogues as their inhibitors.

The catalytic mechanism of the RNA methyltransferase METTL3

The complex of methyltransferase-like proteins 3 and 14 (METTL3-14) is the major enzyme that deposits N6-methyladenosine (m6A) modifications on mRNA in humans. METTL3-14 plays key roles in various biological processes through its methyltransferase (MTase) activity. However, little is known about its substrate recognition and methyl transfer mechanism from its cofactor and methyl donor S-adenosylmethionine (SAM). Here, we study the MTase mechanism of METTL3-14 by a combined experimental and multiscale simulation approach using bisubstrate analogues (BAs), conjugates of a SAM-like moiety connected to the N6-atom of adenosine. Molecular dynamics simulations based on crystal structures of METTL3-14 with BAs suggest that the Y406 side chain of METTL3 is involved in the recruitment of adenosine and release of m6A. A crystal structure with a bisubstrate analogue representing the transition state of methyl transfer shows a direct involvement of the METTL3 side chains E481 and K513 in adenosine binding which is supported by mutational analysis. Quantum mechanics/molecular mechanics (QM/MM) free energy calculations indicate that methyl transfer occurs without prior deprotonation of adenosine-N6. Furthermore, the QM/MM calculations provide further support for the role of electrostatic contributions of E481 and K513 to catalysis. The multidisciplinary approach used here sheds light on the (co)substrate binding mechanism, catalytic step, and (co)product release catalysed by METTL3, and suggests that the latter step is rate-limiting. The atomistic information on the substrate binding and methyl transfer reaction of METTL3 can be useful for understanding the mechanisms of other RNA MTases and for the design of transition state analogues as their inhibitors.

Mapping H+ in the Nanoscale (A2C4)2-Ag8 Fluorophore

When strands of DNA encapsulate silver clusters, supramolecular optical chromophores develop. However, how a particular structure endows a specific spectrum remains poorly understood. Here, we used neutron diffraction to map protonation in (A₂C₄)₂-Ag₈, a green-emitting fluorophore with a "Big Dipper" arrangement of silvers. The DNA host has two substructures with distinct protonation patterns. Three cytosines from each strand collectively chelate handle-like array of three silvers, and calorimetry studies suggest Ag⁺ cross-links. The twisted cytosines are further joined by hydrogen bonds from fully protonated amines. The adenines and their neighboring cytosine from each strand anchor a dipper-like group of five silvers via their deprotonated endo- and exocyclic nitrogens. Typically, exocyclic amines are strongly basic, so their acidification and deprotonation in (A₂C₄)₂-Ag₈ suggest that silvers perturb the electron distribution in the aromatic nucleobases. The different protonation states in (A₂C₄)₂-Ag₈ suggest that atomic level structures can pinpoint how to control and tune the electronic spectra of these nanoscale chromophores.

Exciton decay mechanism in DNA single strands: back-electron transfer and ultrafast base motions

The photochemistry of DNA systems is characterized by the ultraviolet (UV) absorption of π-stacked nucleobases, resulting in exciton states delocalized over several bases. As their relaxation sensitively depends on local stacking conformations, disentangling the ensuing electronic and structural dynamics has remained an experimental challenge, despite their fundamental role in protecting the genome from potentially harmful UV radiation. Here we use transient absorption and transient absorption anisotropy spectroscopy with broadband femtosecond deep-UV pulses (250–360 nm) to resolve the exciton dynamics of UV-excited adenosine single strands under physiological conditions. Due to the exceptional deep-UV bandwidth and polarization sensitivity of our experimental approach, we simultaneously resolve the population dynamics, charge-transfer (CT) character and conformational changes encoded in the UV transition dipoles of the π-stacked nucleotides. Whilst UV excitation forms fully charge-separated CT excitons in less than 0.3 ps, we find that most decay back to the ground state via a back-electron transfer. Based on the anisotropy measurements, we propose that this mechanism is accompanied by a structural relaxation of the photoexcited base-stack, involving an inter-base rotation of the nucleotides. Our results finally complete the exciton relaxation mechanism for adenosine single strands and offer a direct view into the coupling of electronic and structural dynamics in aggregated photochemical systems.

Despite its key role in DNA photochemistry, the decay mechanism of excitons in stacked bases has remained difficult to resolve. Ultrafast polarization spectroscopy now reveals a back-electron transfer and ultrafast base motions in adenosine strands.

Structural and functional characterization of TrmM in m6 A modification of bacterial tRNA

N6 -methyladenosine (m6 A), widely distributed in both coding and noncoding RNAs, regulates the epigenetic signals and RNA metabolism in eukaryotes. Although this posttranscriptional modification is frequently observed in messenger and ribosomal RNA, it is relatively rare in transfer RNA. In Escherichia coli, TrmM encoded by yfiC is the tRNA-specific N6 methyltransferase, which modifies the A37 residue of tRNAVal (cmo5 UAC) using S-adenosyl-l-methionine as a methyl donor. However, the structure-function relationship of this enzyme is not completely understood. In this report, we determined two x-ray crystal structures of Mycoplasma capricolum TrmM with and without S-adenosyl-l-homocysteine, which is a reaction product. We also demonstrated the cellular and in vitro activities of this enzyme in the m6 A modification of tRNA and the requirement of a divalent metal ion for its function, which is unprecedented in other RNA N6 methyltransferases, including the E. coli TrmM. Our results reveal that the dimeric form of M. capricolum TrmM is important for efficient tRNA binding and catalysis, thereby offering insights into the distinct substrate specificity of the monomeric E. coli homolog.

More than Proton Detection-New Avenues for NMR Spectroscopy of RNA

Ribonucleic acid oligonucleotides (RNAs) play pivotal roles in cellular function (riboswitches), chemical biology applications (SELEX-derived aptamers), cell biology and biomedical applications (transcriptomics). Furthermore, a growing number of RNA forms (long non-coding RNAs, circular RNAs) but also RNA modifications are identified, showing the ever increasing functional diversity of RNAs. To describe and understand this functional diversity, structural studies of RNA are increasingly important. However, they are often more challenging than protein structural studies as RNAs are substantially more dynamic and their function is often linked to their structural transitions between alternative conformations. NMR is a prime technique to characterize these structural dynamics with atomic resolution. To extend the NMR size limitation and to characterize large RNAs and their complexes above 200 nucleotides, new NMR techniques have been developed. This Minireview reports on the development of NMR methods that utilize detection on low-γ nuclei (heteronuclei like 13 C or 15 N with lower gyromagnetic ratio than 1 H) to obtain unique structural and dynamic information for large RNA molecules in solution. Experiments involve through-bond correlations of nucleobases and the phosphodiester backbone of RNA for chemical shift assignment and make information on hydrogen bonding uniquely accessible. Previously unobservable NMR resonances of amino groups in RNA nucleobases are now detected in experiments involving conformational exchange-resistant double-quantum 1 H coherences, detected by 13 C NMR spectroscopy. Furthermore, 13 C and 15 N chemical shifts provide valuable information on conformations. All the covered aspects point to the advantages of low-γ nuclei detection experiments in RNA.

Novel 13 C-detected NMR Experiments for the Precise Detection of RNA Structure

Up to now, NMR spectroscopic investigations of RNA have utilized imino proton resonances as reporters for base pairing and RNA structure. The nucleobase amino groups are often neglected, since most of their resonances are broadened beyond detection due to rotational motion around the C-NH₂ bond. Here, we present ¹³ C-detected NMR experiments for the characterization of all RNA amino groups irrespective of their motional behavior. We have developed a C(N)H-HDQC experiment that enables the observation of a complete set of sharp amino resonances through the detection of proton-NH₂ double quantum coherences. Further, we present an "amino"-NOESY experiment to detect NOEs to amino protons, which are undetectable by any other conventional NOESY experiment. Together, these experiments allow the exploration of additional chemical shift information and inter-residual proton distances important for high-resolution RNA secondary and tertiary structure determination.

Insights into Catalytic and tRNA Recognition Mechanism of the Dual-Specific tRNA Methyltransferase from Thermococcus kodakarensis

The tRNA methyltransferase Trm10, conserved throughout Eukarya and Archaea, catalyzes N1-methylation of purine residues at position 9 using S-adenosyl methionine as the methyl donor. The Trm10 family exhibits diverse target nucleotide specificity, with some homologs that are obligate m¹G₉ or m¹A₉-specific enzymes, while others are bifunctional enzymes catalyzing both m¹G₉ and m¹A₉. This variability is particularly intriguing given different chemical properties of the target N1 atom of guanine and adenine. Here we performed an extensive kinetic and mutational analysis of the m¹G₉ and m¹A₉-catalyzing Trm10 from Thermococcus kodakarensis to gain insight into the active site that facilitates this unique bifunctionality. These results suggest that the rate-determining step for catalysis likely involves a conformational change to correctly position the substrate tRNA in the active site. In this model, kinetic preferences for certain tRNA can be explained by variations in the overall stability of the folded substrate tRNA, consistent with tRNA-specific differences in metal ion dependence. Together, these results provide new insight into the substrate recognition, active site and catalytic mechanism of m¹G/m¹A catalyzing bifunctional enzymes.

Rotation of Guanine Amino Groups in G-Quadruplexes: A Probe for Local Structure and Ligand Binding

Nucleic acids are dynamic molecules whose functions may depend on their conformational fluctuations and local motions. In particular, amino groups are dynamic components of nucleic acids that participate in the formation of various secondary structures such as G-quadruplexes. Here, we present a cost-efficient NMR method to quantify the rotational dynamics of guanine amino groups in G-quadruplex nucleic acids. An isolated spectrum of amino protons from a specific tetrad-bound guanine can be extracted from the nuclear Overhauser effect spectroscopy spectrum based on the close proximity between the intra-residue imino and amino protons. We apply the method in different structural contexts of G-quadruplexes and their complexes. Our results highlight the role of stacking and hydrogen-bond interactions in restraining amino-group rotation. The measurement of the rotation rate of individual amino groups could give insight into the dynamic processes occurring at specific locations within G-quadruplex nucleic acids, providing valuable probes for local structure, dynamics, and ligand binding.

The m1A(58) modification in eubacterial tRNA: An overview of tRNA recognition and mechanism of catalysis by TrmI

The enzymes of the TrmI family catalyze the formation of the m(1)A58 modification in tRNA. We previously solved the crystal structure of the Thermus thermophilus enzyme and conducted a biophysical study to characterize the interaction between TrmI and tRNA. TrmI enzymes are active as a tetramer and up to two tRNAs can bind to TrmI simultaneously. In this paper, we present the structures of two TrmI mutants (D170A and Y78A). These residues are conserved in the active site of TrmIs and their mutations result in a dramatic alteration of TrmI activity. Both structures of TrmI mutants revealed the flexibility of the N-terminal domain that is probably important to bind tRNA. The structure of TrmI Y78A catalytic domain is unmodified regarding the binding of the SAM co-factor and the conformation of residues potentially interacting with the substrate adenine. This structure reinforces the previously proposed role of Y78, i.e. stabilize the conformation of the A58 ribose needed to hold the adenosine in the active site. The structure of the D170A mutant shows a flexible active site with one loop occupying in part the place of the co-factor and the second loop moving at the entrance to the active site. This structure and recent data confirms the central role of D170 residue binding the amino moiety of SAM and the exocyclic amino group of adenine. Possible mechanisms for methyl transfer are then discussed.

Fluctuations in the DNA double helix: A critical review

A critical overview of the extensive literature on fluctuations in the DNA double helix is presented. Both theory and experiment are comprehensively reviewed and analyzed. Fluctuations, which open up the DNA double helix making bases accessible for hydrogen exchange and chemical modification, are the main focus of the review. Theoretical descriptions of the DNA fluctuations are discussed with special emphasis on most popular among them: the nonlinear-dynamic Peyrard-Bishop-Dauxois (PBD) model and the empirical two-state (or helix-coil) model. The experimental data on the issue are comprehensibly overviewed in the historical retrospective with main emphasis on the hydrogen exchange data and formaldehyde kinetics. The theoretical descriptions are critically evaluated from the viewpoint of their applicability to describe DNA in water environment and from the viewpoint of agreement of their predictions with the reliable experimental data. The presented analysis makes it possible to conclude that, while the two-state model is most adequate from theoretical viewpoint and its predictions, based on an empirical parametrization, agree with experimental data very well, the PBD model is inapplicable to DNA in water from theoretical viewpoint on one hand and it makes predictions totally incompatible with reliable experimental data on the other. In particular, it is argued that any oscillation movements of nucleotides, assumed by the PBD model, are severely damped in water, that no "bubbles", which the PBD model predicts, exist in reality in linear DNA well below the melting range and the lifetime of an open state in DNA is actually 5 orders of magnitude longer than the value predicted by the PBD model.

Methylation stabilizes the imino tautomer of dAMP and amino tautomer of dCMP in solution

Alkylating agents cause methylation of adenosine and cytidine in DNA to generate 1-methyladenosine and 3-methylcytidine. These modified nucleosides can serve as regulators of cells or can act as agents of mutagenesis depending on the context and the partner enzymes. Solution structures and the chemical interactions with enzymes that lead to their recognition are of inherent interest. At physiological pH, 1-methyladenosine and 3-methylcytidine are presumed to be in the protonated amino forms in the literature. We report the structures, ionization states, and UV resonance Raman spectra of both substrates over a range of pH (2.5-11.0). The Raman excitation wavelength was tuned to selectively enhance Raman scattering from the nucleobase (260 nm) and further specifically from the imino form (210 nm) of 1-me-dAMP. We find that contrary to the general assumption, 1-me-dAMP is present in its neutral imino form at physiological pH and 3-me-dCMP is in the amino form.

Crystal structure of Thermus thermophilus tRNA m1A58 methyltransferase and biophysical characterization of its interaction with tRNA

Methyltransferases from the m(1)A(58) tRNA methyltransferase (TrmI) family catalyze the S-adenosyl-l-methionine-dependent N(1)-methylation of tRNA adenosine 58. The crystal structure of Thermus thermophilus TrmI, in complex with S-adenosyl-l-homocysteine, was determined at 1.7 A resolution. This structure is closely related to that of Mycobacterium tuberculosis TrmI, and their comparison enabled us to enlighten two grooves in the TrmI structure that are large enough and electrostatically compatible to accommodate one tRNA per face of TrmI tetramer. We have then conducted a biophysical study based on electrospray ionization mass spectrometry, site-directed mutagenesis, and molecular docking. First, we confirmed the tetrameric oligomerization state of TrmI, and we showed that this protein remains tetrameric upon tRNA binding, with formation of complexes involving one to two molecules of tRNA per TrmI tetramer. Second, three key residues for the methylation reaction were identified: the universally conserved D170 and two conserved aromatic residues Y78 and Y194. We then used molecular docking to position a N(9)-methyladenine in the active site of TrmI. The N(9)-methyladenine snugly fits into the catalytic cleft, where the side chain of D170 acts as a bidentate ligand binding the amino moiety of S-adenosyl-l-methionine and the exocyclic amino group of the adenosine. Y194 interacts with the N(9)-methyladenine ring, whereas Y78 can stabilize the sugar ring. From our results, we propose that the conserved residues that form the catalytic cavity (D170, Y78, and Y194) are essential for fashioning an optimized shape of the catalytic pocket.

Not All G-Quadruplexes Exhibit Ion-Channel-like Properties: NMR Study of Ammonium Ion (Non)movement within the d(G3T4G4)2 Quadruplex

A solution-state NMR study on 15NH4(+) ion movement within d(G(3)T(4)G(4))(2), a dimeric G-quadruplex consisting of three G-quartets and two T(4) loops, rather unexpectedly demonstrated the absence of 15NH4(+) ion movement between the binding sites U and L along the central axis of the G-quadruplex. Distinct temperature dependences of autocorrelation signals for U and L binding sites have been observed in 15N-1H NzExHSQC spectra which correlate with the local stiffness of the G-quadruplex. The volumes of the cross-peaks, which are the result of 15NH4(+) ion movement, have been interpreted in terms of rate constants, T(1) relaxation, and proton exchange. 15NH4(+) ion movements from the binding sites U and L into the bulk solution are characterized by lifetimes of 139 ms and 1.7 s at 298 K, respectively. The 12 times faster movement from the binding site U demonstrates that 15NH4(+) ion movement is controlled by the structure of T4 loop residues, which through diagonal- vs edge-type orientations impose distinct steric restraints for cations to leave or enter the G-quadruplex. Arrhenius-type analysis has afforded an activation energy of 66 kJ mol(-)1 for the UB process, while it could not be determined for the LB process due to slow rates at temperatures below 298 K. We further the use of the 15NH4(+) ion as an NMR probe to gain insight into the occupancy of binding sites by cations and kinetics of ion movement which are intrinsically correlated with the structural details, dynamic fluctuations, and local flexibility of the DNA structure.

Mutagenesis, genotoxicity, and repair of 1-methyladenine, 3-alkylcytosines, 1-methylguanine, and 3-methylthymine in alkB Escherichia coli

AlkB repairs 1-alkyladenine and 3-methylcytosine lesions in DNA by directly reversing the base damage. Although repair studies with randomly alkylated substrates have been performed, the miscoding nature of these and related individually alkylated bases and the suppression of mutagenesis by AlkB within cells have not yet been explored. Here, we address the miscoding potential of 1-methyldeoxyadenosine (m1A), 3-methyldeoxycytidine (m3C), 3-ethyldeoxycytidine (e3C), 1-methyldeoxyguanosine (m1G), and 3-methyldeoxythymidine (m3T) by synthesizing single-stranded vectors containing each alkylated base, followed by vector passage through Escherichia coli. In SOS(-), AlkB-deficient cells, m1A was only 1% mutagenic; however, m3C and e3C were 30% mutagenic, rising to 70% in SOS(+) cells. In contrast, the mutagenicity of m1G and m3T in AlkB(-) cells dropped slightly when SOS polymerases were expressed (m1G from 80% to 66% and m3T from 60% to 53%). Mutagenicity was abrogated for m1A, m3C, and e3C in wild-type (AlkB(+)) cells, whereas m3T mutagenicity was only partially reduced. Remarkably, m1G mutagenicity was also eliminated in AlkB(+) cells, establishing it as a natural AlkB substrate. All lesions were blocks to replication in AlkB-deficient cells. The m1A, m3C, and e3C blockades were completely removed in wild-type cells; the m1G blockade was partially removed and that for m3T was unaffected by the presence of AlkB. All lesions demonstrated enhanced bypass when SOS polymerases were induced. This work provides direct evidence that AlkB suppresses both genotoxicity and mutagenesis by physiologically realistic low doses of 1-alkylpurine and 3-alkylpyrimidine DNA damage in vivo.

The RNA i-motif

Oligodeoxynucleotides with stretches of cytidine residues associate into a four-stranded structure, the i-motif, in which two head-to-tail, intercalated, parallel-stranded duplexes are held together by hemiprotonated C.C+ pairs. We have investigated the possibility of forming an i-motif structure with C-rich ribonucleic acids. The four C-rich RNAs studied, r(UC5), r(C5), r(C5U) and r(UC3), associate into multiple intercalated structures at acidic pH. r(UC5) forms two i-motif structures that differ by their intercalation topologies. We report on a structural study of the main form and we analyze the small conformational differences found by comparison with the DNA i-motif. The stacking topology of the main structure avoids one of the six 2'-OH/2'-OH repulsive contacts expected in a fully intercalated structure. The C3'-endo pucker of the RNA sugars and the orientation of the intercalated C.C+ pairs result in a modest widening of the narrow grooves at the steps where the hydroxyl groups are in close contact. The free energy of the RNA i-motif, on average -4 kJ mol(-1) per C.C+ pair, is half of the value found in DNA i-motif structures.

Determination of the residence time of water molecules hydrating B'- DNA and B-DNA, by one-dimensional zero-enhancement nuclear Overhauser effect spectroscopy

The residence time of water in the minor groove of the d(CGCGAATTCGCG) duplex has been determined by a recent measurement combining nuclear Overhauser enhancements (NOE, ROE) and 17O relaxation dispersion. The time is in the range of nanoseconds, so that it may be measured by a rather simple method proposed here, namely the choice of conditions such that the NOE between the observed DNA proton and a nearby water proton is zero. This condition is realized when the residence time of the water molecule is 0.178 times the nuclear magnetic resonance period (e.g. 0.297 ns at 600 MHz). It may be achieved by varying the magnetic field and/or the temperature. The zero-NOE measurement may be performed by one-dimensional NMR, and has therefore good sensitivity. We have developed excitation sequences which suppress two spurious contributions to the NOE: from neighboring exchangeable protons and from H3' protons whose chemical shift is close to that of water. The method is applied here to the comparison of residence times of water next to B-DNA and next to B'-DNA, the latter corresponding to better stacked, propeller-twisted base-pairs and a correspondingly narrower minor groove. In the minor groove of [d(CGCGAATTCGCG)]2, a B'-DNA duplex, the residence time of the water molecule next to H2 of adenine(6) (underlined), is 0.6 ns at 10 degreesC, in good agreement with the value obtained previously. The residence time is slightly but distinctly shorter for the water next to A5, suggesting non-cooperative departure of these two molecules which are presumed to be part of the hydration spine. Near A5 and A4 of [d(AAAAATTTTT)]2, another B'-DNA duplex, the residence times are approximately twice as long, but the activation enthalpies are about the same, ca. 38 kJ/mol. The residence time in the minor groove of the regular B-DNA sequence d(CGCGATCGCG) was 0.3 ns at 10 degreesC, shorter than in the case of the B'-DNA sequences by factors of 2 and 4, respectively. The temperature dependence is less, with an activation enthalpy of 27 kJ/mol. The major groove residence times are comparable for the three sequences, and a few times shorter than those of minor groove water. A value of 0.36 ns, or even more in case of rotation of water, is obtained around -8 degreesC. The most striking aspect of these results is the relatively small difference in the residence times of reputedly fast and slow-exchanging water molecules bound to DNA in biological conditions. This suggests that the spine of hydration is perhaps not a major stabilizer of the B'-DNA structure as compared with B-DNA.