Metal ion dependence of DNA cleavage by SepMI and EhoI restriction endonucleases

Is Metal Stabilization of the Leaving Group Required or Can Lysine Facilitate Phosphodiester Bond Cleavage in Nucleic Acids? A Computational Study of EndoV

Endonuclease V (EndoV) is a single-metal-dependent enzyme that repairs deaminated DNA nucleobases in cells by cleaving the phosphodiester bond, and this enzyme has proven to be a powerful tool in biotechnology and medicine. The catalytic mechanism used by EndoV must be understood to design new disease detection and therapeutic solutions and further exploit the enzyme in interdisciplinary applications. This study has used a mixed molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) approach to compare eight distinct catalytic pathways and provides the first proposed mechanism for bacterial EndoV. The calculations demonstrate that mechanisms involving either direct or indirect metal coordination to the leaving group of the substrate previously proposed for other nucleases are unlikely for EndoV, regardless of the general base (histidine, aspartate, and substrate phosphate moiety). Instead, distinct catalytic pathways are characterized for EndoV that involve K139 stabilizing the leaving group, a metal-coordinated water stabilizing the transition structure, and either H214 or a substrate phosphate group activating the water nucleophile. In silico K139A and H214A mutational results support the newly proposed roles of these residues. Although this is a previously unseen combination of general base, general acid, and metal-binding architecture for a one-metal-dependent endonuclease, our proposed catalytic mechanisms are fully consistent with experimental kinetic, structural, and mutational data. In addition to substantiating a growing body of literature, suggesting that one metal is enough to catalyze P-O bond cleavage in nucleic acids, this new fundamental understanding of the catalytic function will promote the exploration of new and improved applications of EndoV.

In-silicon studies on hydration in EcoRI-cognate DNA complex

Restriction endonucleases (REs) cleave DNA at specific site in presence of Mg2+ ion. Experiments further emphasize the role of hydration in metal ion specificity and sequence specificity of DNA cleavage. However, the relation between hydration and specificity has not been understood till date. This leads us to study via all-atom molecular dynamics (MD) simulations how the hydration around the scissile phosphate group changes in presence of Mg2+ and Ca2+ and depend on the DNA sequence. We observe the least number of hydrogen bonds around the scissile phosphate group in presence of Mg2+ ion. We further find that the hydrogen bonds decrease at the scissile phosphate on mutating one base pair in the cleavage region of the DNA in Mg2+ loaded EcoRI-DNA complex. We also perform steered MD simulations and observe that the rate of decrease of fraction of hydrogen bonds is slower in the mutated complex than the unmutated complex.

The Impact of DFT Functional, Cluster Model Size, and Implicit Solvation on the Structural Description of Single-Metal-Mediated DNA Phosphodiester Bond Cleavage: The Case Study of APE1

Phosphodiester bond hydrolysis in nucleic acids is a ubiquitous reaction that can be facilitated by enzymes called nucleases, which often use metal ions to achieve catalytic function. While a two-metal-mediated pathway has been well established for many enzymes, there is growing support that some enzymes require only one metal for the catalytic step. Using human apurinic/apyrimidinic endonuclease (APE1) as a prototypical example and cluster models, this study clarifies the impact of DFT functional, cluster model size, and implicit solvation on single-metal-mediated phosphodiester bond cleavage and provides insight into how to efficiently model this chemistry. Initially, a model containing 69 atoms built from a high-resolution X-ray crystal structure is used to explore the reaction pathway mapped by a range of DFT functionals and basis sets, which provides support for the use of standard functionals (M06-2X and B3LYP-D3) to study this reaction. Subsequently, systematically increasing the model size to 185 atoms by including additional amino acids and altering residue truncation points highlights that small models containing only a few amino acids or β carbon truncation points introduce model strains and lead to incorrect metal coordination. Indeed, a model that contains all key residues (general base and acid, residues that stabilize the substrate, and amino acids that maintain the metal coordination) is required for an accurate structural depiction of the one-metal-mediated phosphodiester bond hydrolysis by APE1, which results in 185 atoms. The additional inclusion of the broader enzyme environment through continuum solvation models has negligible effects. The insights gained in the present work can be used to direct future computational studies of other one-metal-dependent nucleases to provide a greater understanding of how nature achieves this difficult chemistry.

The metal dependence of single-metal mediated phosphodiester bond cleavage: a QM/MM study of a multifaceted human enzyme

Nucleases catalyze the cleavage of phosphodiester bonds in nucleic acids using a range of metal cofactors. Although it is well accepted that many nucleases rely on two metal ions, the one-metal mediated pathway is debated. Furthermore, one-metal mediated nucleases maintain activity in the presence of many different metals, but the underlying reasons for this broad metal specificity are unknown. The human apurinic/apyrimidinic endonuclease (APE1), which plays a key role in DNA repair, transcription regulation, and gene expression, is a prototypical example of a one-metal dependent nuclease. Although Mg²⁺ is the native metal cofactor, APE1 remains catalytically active in the presence of several metals, with the rate decreasing as Mg²⁺ > Mn²⁺ > Ni²⁺ > Zn²⁺, while Ca²⁺ completely abolished the activity. The present work uses quantum mechanics-molecular mechanics techniques to map APE1-facilitated phosphodiester bond hydrolysis in the presence of these metals. The structural differences in stationary points along the reaction pathway shed light on the interplay between several factors that allow APE1 to remain catalytically active for various metals, with the trend in the barrier heights correlating with the experimentally reported APE1 catalytic activity. In contrast, Ca²⁺ significantly changes the metal coordination and active site geometry, and thus completely inhibits catalysis. Our work thereby provides support for the controversial single-metal mediated phosphodiester bond cleavage and clarifies uncertainties regarding the role of the metal and metal identity in this important reaction. This information is key for future medicinal and biotechnological applications including disease diagnosis and treatment, and protein engineering.

Assembly of Cas7 subunits of Leptospira on the mature crRNA of CRISPR-Cas I-B is modulated by divalent ions

The spirochete Leptospira interrogans serovar Copenhageni harbors the genetic elements of the CRISPR-Cas type I-B system in its genome. CRISPR-Cas is a CRISPR RNA (crRNA) mediated adaptive immune system in most prokaryotes against mobile genetic elements (MGEs). To eliminate the intruding MGEs, CRISPR-Cas type I systems utilize a Cascade (CRISPR-associated complex for antiviral defense) complex composed of Cas5, Cas6, Cas7, and Cas8 bound with a crRNA. The Cas7 is essentially known to constitute the major component of the Cascade complex. The present study reports the biochemical characterization of the Cas7 (LinCas7) from the CRISPR-Cas type I-B system of L. interrogans serovar Copenhageni. The pure recombinant LinCas7 (rLinCas7) exists as a monomer in the solution by size exclusion chromatography. The rLinCas7 demonstrates an endoDNase activity dependent upon divalent Mg²⁺ ions, monovalent ions, pH, temperature, and substrate size. Analysis of ribonucleoprotein composite (rLinCas7-crRNA) by electron microscopy and native-PAGE demonstrated that rLinCas7 could oligomerize on the mature CRISPR RNA (crRNA) framework in the presence of Mg²⁺ ions. The ribonucleoprotein composite attains a helical shape similar to the backbone of the Cascade complex. However, in the absence of Mg²⁺ ions, rLinCas7 acts as an RNase. The fluorescence spectroscopy disclosed a weak interaction (K_d = 26.81 mM) between rLinCas7 and Mg²⁺ ions, leading to an overall conformational change in rLinCas7 that modulates the rLinCas7's activity on DNA and RNA substrates. The nuclease activity of LinCas7 characterized in this study aids to the functional divergences among proteins of the Cas7 family from different CRISPR-Cas systems in various organisms.

Microscopic insight to specificity of metal ion cofactor in DNA cleavage by restriction endonuclease EcoRV

Restriction endonucleases protect bacterial cells against bacteriophage infection by cleaving the incoming foreign DNA into fragments. In presence of Mg²⁺ ions, EcoRV is able to cleave the DNA but not in presence of Ca²⁺ , although the protein binds to DNA in presence of both metal ions. We make an attempt to understand this difference using conformational thermodynamics. We calculate the changes in conformational free energy and entropy of conformational degrees of freedom, like DNA base pair steps and dihedral angles of protein residues in Mg²⁺ (A)-EcoRV-DNA complex compared to Ca²⁺ (S)-EcoRV-DNA complex using all-atom molecular dynamics (MD) trajectories of the complexes. We find that despite conformational stability and order in both complexes, the individual degrees of freedom behave differently in the presence of two different metal ions. The base pairs in cleavage region are highly disordered in Ca²⁺ (S)-EcoRV-DNA compared to Mg²⁺ (A)-EcoRV-DNA. One of the acidic residues ASP90, coordinating to the metal ion in the vicinity of the cleavage site, is conformationally destabilized and disordered, while basic residue LYS92 gets conformational stability and order in Ca²⁺ (S) bound complex than in Mg²⁺ (A) bound complex. The enhanced fluctuations hinder placement of the metal ion in the vicinity of the scissile phosphate of DNA. Similar loss of conformational stability and order in the cleavage region is observed by the replacement of the metal ion. Considering the placement of the metal ion near scissile phosphate as requirement for cleavage action, our results suggest that the changes in conformational stability and order of the base pair steps and the protein residues lead to cofactor sensitivity of the enzyme. Our method based on fluctuations of microscopic conformational variables can be applied to understand enzyme activities in other protein-DNA systems.