Biochemical and mechanistic analysis of the cleavage of branched DNA by human ANKLE1

ANKLE1 is a nuclease that provides a final opportunity to process unresolved junctions in DNA that would otherwise create chromosomal linkages blocking cell division. It is a GIY-YIG nuclease. We have expressed an active domain of human ANKLE1 containing the GIY-YIG nuclease domain in bacteria, that is monomeric in solution and when bound to a DNA Y-junction, and unilaterally cleaves a cruciform junction. Using an AlphaFold model of the enzyme we identify the key active residues, and show that mutation of each leads to impairment of activity. There are two components in the catalytic mechanism. Cleavage rate is pH dependent, corresponding to a pKa of 6.9, suggesting an involvement of the conserved histidine in proton transfer. The reaction rate depends on the nature of the divalent cation, likely bound by glutamate and asparagine side chains, and is log-linear with the metal ion pKa. We propose that the reaction is subject to general acid-base catalysis, using a combination of tyrosine and histidine acting as general base and water directly coordinated to the metal ion as general acid. The reaction is temperature dependent; activation energy Ea = 37 kcal mol-1, suggesting that cleavage is coupled to opening of DNA in the transition state.

Catalytic mechanism and pH dependence of a methyltransferase ribozyme (MTR1) from computational enzymology

A methyltransferase ribozyme (MTR1) was selected in vitro to catalyze alkyl transfer from exogenous O6-methylguanine (O6mG) to a target adenine N1, and recently, high-resolution crystal structures have become available. We use a combination of classical molecular dynamics, ab initio quantum mechanical/molecular mechanical (QM/MM) and alchemical free energy (AFE) simulations to elucidate the atomic-level solution mechanism of MTR1. Simulations identify an active reactant state involving protonation of C10 that hydrogen bonds with O6mG:N1. The deduced mechanism involves a stepwise mechanism with two transition states corresponding to proton transfer from C10:N3 to O6mG:N1 and rate-controlling methyl transfer (19.4  kcal·mol-1 barrier). AFE simulations predict the pKa for C10 to be 6.3, close to the experimental apparent pKa of 6.2, further implicating it as a critical general acid. The intrinsic rate derived from QM/MM simulations, together with pKa calculations, enables us to predict an activity-pH profile that agrees well with experiment. The insights gained provide further support for a putative RNA world and establish new design principles for RNA-based biochemical tools.

Ribocentre: a database of ribozymes

Ribozymes are excellent systems in which to study 'sequence - structure - function' relationships in RNA molecules. Understanding these relationships may greatly help structural modeling and design of functional RNA structures and some functional structural modules could be repurposed in molecular design. At present, there is no comprehensive database summarising all the natural ribozyme families. We have therefore created Ribocentre, a database that collects together sequence, structure and mechanistic data on 21 ribozyme families. This includes available information on timelines, sequence families, secondary and tertiary structures, catalytic mechanisms, applications of the ribozymes together with key publications. The database is publicly available at https://www.ribocentre.org.