Role of magnesium ions in DNA recognition by the EcoRV restriction endonuclease

Docking Ligands into Flexible and Solvated Macromolecules. 8. Forming New Bonds─Challenges and Opportunities

Over the years, structure-based design programs and specifically docking small molecules to proteins have become prominent in drug discovery. However, many of these computational tools have been developed to primarily dock enzyme inhibitors (and ligands to other protein classes) relying heavily on hydrogen bonds and electrostatic and hydrophobic interactions. In reality, many drug targets either feature metal ions, can be targeted covalently, or are simply not even proteins (e.g., nucleic acids). Herein, we describe several new features that we have implemented into Fitted to broaden its applicability to a wide range of covalent enzyme inhibitors and to metalloenzymes, where metal coordination is essential for drug binding. This updated version of our docking program was tested for its ability to predict the correct binding mode of drug-sized molecules in a large variety of proteins. We also report new datasets that were essential to demonstrate areas of success and those where additional efforts are required. This resource could be used by other program developers to assess their own software.

On the Role of Magnesium Ions in the DNA-Scissoring Activity of the Restriction Endonuclease ApaI: Stochastic Kinetics from a Single Molecule to Mesoscopic Paradigm

Biochemical reactions occurring inside cells have significant stochastic signatures due to the low copy number of reacting species. Kinetics of DNA cleavage by restriction endonucleases are no exception as established by single-molecule experiments. Here, we propose a simple reaction scheme to understand the role of the cofactor magnesium ion in the action of the endonuclease ApaI. The methodology is based on the waiting time distribution of cleavage product formation that enables us to determine the corresponding rate both analytically and numerically. The theory is developed at the single-molecule level and then generalized to the biologically relevant case of a population of DNA-endonuclease complexes present inside a cell. The theoretical rate versus cofactor concentration curve is matched with relevant single-molecule experimental data that reveals positive cooperativity of cofactor binding and provides a reliable estimate of model parameters. Furthermore, a parameter range is identified where the dispersion of the waiting time, measured using the coefficient of variation, is significantly lower than the Poisson limit and becomes minimum at the in vivo magnesium ion concentration level. Such low dispersion can play a role in the robust DNA-scissoring activity of ApaI under in vivo conditions.

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.

Molecular movie of nucleotide binding to a motor protein

BACKGROUND

The SecA DEAD (Asp-Glu-Ala-Asp) motor protein uses binding and hydrolysis of adenosine triphosphate (ATP) to push secretory proteins across the plasma membrane of bacteria. The reaction coordinate of nucleotide exchange is unclear at the atomic level of detail.

METHODS

We performed multiple atomistic computations of the DEAD motor domain of SecA with different occupancies of the nucleotide and magnesium ion sites, for a total of ~1.7 μs simulation time. To characterize dynamics at the active site we analyzed hydrogen-bond networks.

RESULTS

ATP and ADP can bind spontaneously at the interface between the nucleotide binding domains, albeit at an intermediate binding site distinct from the native site. Binding of the nucleotide is facilitated by the presence of a magnesium ion close to the glutamic group of the conserved DEAD motif. In the absence of the magnesium ion, protein interactions of the ADP molecule are perturbed.

CONCLUSIONS

A protein hydrogen-bond network whose dynamics couples to the occupancy of the magnesium ion site helps guide the nucleotide along the nucleotide exchange path. In SecA, release of magnesium might be required to destabilize the ADP binding site prior to release of the nucleotide.

GENERAL SIGNIFICANCE

We identified dynamic hydrogen-bond networks that help control nucleotide exchange in SecA, and stabilize ADP at an intermediate site that could explain slow release. The reaction coordinate of the protein motor involves complex rearrangements of a hydrogen-bond network at the active site, with perturbation of the magnesium ion site likely occurring prior to the release of ADP.

Role of Mg2+ Ions in DNA Hydrolysis by EcoRV, Studied by the 3D-Reference Interaction Site Model and Molecular Dynamics

The role of Mg²⁺ ions during precursor formation in DNA hydrolysis by the homodimeric restriction enzyme EcoRV was elucidated based on the 3D-reference interaction site model (RISM) theory and the molecular dynamics (MD) simulation. From an analysis of the spatial distribution of Mg²⁺ in an active site using 3D-RISM, we identified a new position for Mg²⁺ in the X-ray EcoRV-DNA complex structure ( 1rvb ). We refer to the position as site IV^†. Site IV^† is almost at the same position as that of a Ca²⁺ ion in the superimposed X-ray crystallographic active-site structure of the PvuII-DNA complex ( 1f0o ). 3D-RISM was also used to locate the position of water molecules, including the water nucleophile at the active site. MD simulations were carried out with the initial structure having two Mg²⁺ ions at site IV^† and at site I*, experimentally identified by Horton et al., to find a stable complex structure in which the DNA fragment was rearranged to orient the scissile bond direction toward the water nucleophile. The equilibrium active-site structure of the EcoRV-DNA complex obtained from the MD simulation was similar to the superimposed X-ray crystallographic structure of the BamHI-DNA complex ( 2bam ). In the active-site structure, two metal ions have almost the same position (≤1.0 Å) as that of 2bam , and the scissile phosphate is twisted to orient the scissile bond toward the water nucleophile, as is the case in 2bam . We propose the equilibrium active-site structure obtained in this study as a precursor for the hydrolysis reaction of EcoRV.

Phosphorothioated Primers Lead to Loop-Mediated Isothermal Amplification at Low Temperatures

Loop-mediated isothermal amplification (LAMP) is an extremely powerful tool for the detection of nucleic acids with high sensitivity and specificity. However, LAMP shows optimal performance at around 65 °C, which limits applications in point-of-care-testing (POCT). Here, we have developed a version of LAMP that uses phosphorothioated primers (PS-LAMP) to enable more efficient hairpin formation and extension at the termini of growing concatamers, and that therefore works at much lower temperatures. By including additional factors such as chaotropes (urea) and single-stranded DNA binding protein (SSB), the sensitivities and selectivities for amplicon detection with PS-LAMP at 40 °C were comparable with a regular LAMP reaction at 65 °C.

Role of magnesium ions in the reaction mechanism at the interface between Tm1631 protein and its DNA ligand

A protein, Tm1631 from the hyperthermophilic organism Thermotoga maritima belongs to a domain of unknown function protein family. It was predicted that Tm1631 binds with the DNA and that the Tm1631–DNA complex is an endonuclease repair system with a DNA repair function (Konc et al. PLoS Comput Biol 9(11): e1003341, 2013). We observed that the severely bent, strained DNA binds to the protein for the entire 90 ns of classical molecular dynamics (MD) performed; we could observe no significant changes in the most distorted region of the DNA, where the cleavage of phosphodiester bond occurs. In this article, we modeled the reaction mechanism at the interface between Tm1631 and its proposed ligand, the DNA molecule, focusing on cleavage of the phosphodiester bond. After addition of two Mg²⁺ ions to the reaction center and extension of classical MD by 50 ns (totaling 140 ns), the DNA ligand stayed bolted to the protein. Results from density functional theory quantum mechanics/molecular mechanics (QM/MM) calculations suggest that the reaction is analogous to known endonuclease mechanisms: an enzyme reaction mechanism with two Mg²⁺ ions in the reaction center and a pentacovalent intermediate. The minimum energy pathway profile shows that the phosphodiester bond cleavage step of the reaction is kinetically controlled and not thermodynamically because of a lack of any energy barrier above the accuracy of the energy profile calculation. The role of ions is shown by comparing the results with the reaction mechanisms in the absence of the Mg²⁺ ions where there is a significantly higher reaction barrier than in the presence of the Mg²⁺ ions.Graphical abstract

Using single-turnover kinetics with osmotic stress to characterize the EcoRV cleavage reaction

Type II restriction endonucleases require metal ions to specifically cleave DNA at canonical sites. Despite the wealth of structural and biochemical information, the number of Mg(2+) ions used for cleavage by EcoRV, in particular, at physiological divalent ion concentrations has not been established. In this work, we employ a single-turnover technique that uses osmotic stress to probe reaction kinetics between an initial specific EcoRV-DNA complex formed in the absence of Mg(2+) and the final cleavage step. With osmotic stress, complex dissociation before cleavage is minimized and the reaction rates are slowed to a convenient time scale of minutes to hours. We find that cleavage occurs by a two-step mechanism that can be characterized by two rate constants. The dependence of these rate constants on Mg(2+) concentration and osmotic pressure gives the number of Mg(2+) ions and water molecules coupled to each kinetic step of the EcoRV cleavage reaction. Each kinetic step is coupled to the binding 1.5-2.5 Mg(2+) ions, the uptake of ∼30 water molecules, and the cleavage of a DNA single strand. We suggest that each kinetic step reflects an independent, rate-limiting conformational change of each monomer of the dimeric enzyme that allows Mg(2+) ion binding. This modified single-turnover protocol has general applicability for metalloenzymes.

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

Most of type II restriction endonucleases show an absolute requirement for divalent metal ions as cofactors for DNA cleavage. While Mg(2+) is the natural cofactor other metal ions can substitute it and mediate the catalysis, however Ca(2+) (alone) only supports DNA binding. To investigate the role of Mg(2+) in DNA cleavage by restriction endonucleases, we have studied the Mg(2+) and Mn(2+) concentration dependence of DNA cleavage by SepMI and EhoI. Digestion reactions were carried out at different Mg(2+) and Mn(2+) concentrations at constant ionic strength. These enzymes showed different behavior regarding the ions requirement, SepMI reached near maximal level of activity between 10 and 20mM while no activity was detected in the presence of Mn(2+) and in the presence of Ca(2+) cleavage activity was significantly decreased. However, EhoI was more highly active in the presence of Mn(2+) than in the presence of Mg(2+) and can be activated by Ca(2+). Our results propose the two-metal ion mechanism for EhoI and the one-metal ion mechanism for SepMI restriction endonuclease. The analysis of the kinetic parameters under steady state conditions showed that SepMI had a K(m) value for pTrcHisB DNA of 6.15 nM and a V(max) of 1.79×10(-2)nM min(-1), while EhoI had a K(m) for pUC19 plasmid of 8.66 nM and a V(max) of 2×10(-2)nM min(-1).