ACS Publications. Most Trusted. Most Cited. Most Read
Recruiting Mechanism and Functional Role of a Third Metal Ion in the Enzymatic Activity of 5′ Structure-Specific Nucleases
My Activity

Figure 1Loading Img
  • Open Access
Article

Recruiting Mechanism and Functional Role of a Third Metal Ion in the Enzymatic Activity of 5′ Structure-Specific Nucleases
Click to copy article linkArticle link copied!

Open PDFSupporting Information (2)

Journal of the American Chemical Society

Cite this: J. Am. Chem. Soc. 2020, 142, 6, 2823–2834
Click to copy citationCitation copied!
https://doi.org/10.1021/jacs.9b10656
Published January 15, 2020

Copyright © 2020 American Chemical Society. This publication is licensed under

CC-BY 4.0 .

Abstract

Click to copy section linkSection link copied!

Enzymes of the 5′ structure-specific nuclease family are crucial for DNA repair, replication, and recombination. One such enzyme is the human exonuclease 1 (hExo1) metalloenzyme, which cleaves DNA strands, acting primarily as a processive 5′-3′ exonuclease and secondarily as a 5′-flap endonuclease. Recently, in crystallo reaction intermediates have elucidated how hExo1 exerts hydrolysis of DNA phosphodiester bonds. These hExo1 structures show a third metal ion intermittently bound close to the two-metal-ion active site, to which recessed ends or 5′-flap substrates bind. Evidence of this third ion has been observed in several nucleic-acid-processing metalloenzymes. However, there is still debate over what triggers the (un)binding of this transient third ion during catalysis and whether this ion has a catalytic function. Using extended molecular dynamics and enhanced sampling free-energy simulations, we observed that the carboxyl side chain of Glu89 (located along the arch motif in hExo1) flips frequently from the reactant state to the product state. The conformational flipping of Glu89 allows one metal ion to be recruited from the bulk and promptly positioned near the catalytic center. This is in line with the structural evidence. Additionally, our simulations show that the third metal ion assists the departure, through the mobile arch, of the nucleotide monophosphate product from the catalytic site. Structural comparisons of nuclease enzymes suggest that this Glu(Asp)-mediated mechanism for third ion recruitment and nucleic acid hydrolysis may be shared by other 5′ structure-specific nucleases.

This publication is licensed under

CC-BY 4.0 .
  • cc licence
  • by licence
Copyright © 2020 American Chemical Society

Introduction

Click to copy section linkSection link copied!

Recent structural data have shown the recurring presence of a third metal ion close to the two-metal-ion center of nucleic-acid-processing enzymes. (1−4) This third ion has been captured during different stages of catalysis of vital enzymatic reactions involved in DNA repair, recombination, and replication processes. (5−9) These reactions are often related to cancer progression. (10−13) Indeed, over the past few years, a third ion has been observed, or hypothesized, close to the two-metal-ion catalytic site of polymerases, (4) nucleases, (14−18) and topoisomerases. (19−21) This suggests that the third metal ion may be actively involved in catalysis. (1,4,22) However, there is debate over how this ion is recruited from the bulk and transiently binds the enzyme and how it could play a role in catalysis.
In this context, recent time-resolved in crystallo reaction intermediates (23) have elucidated how human exonuclease 1 (hExo1) exerts its catalytic function, with sequential structures showing how the enzyme/DNA complex evolves during catalysis. hExo1 is an essential hydrolytic enzyme for genome maintenance. Belonging to the RAD2/XPG family, (24−30) hExo1 is a 5′ structure-specific metallonuclease, which carries out a primary exonucleolytic activity on the 5′ recessed end and a secondary endonucleolytic cleavage on 5′-flap of the substrate DNA strand. (31,32)
The structures show the enzymatic mechanism for DNA hydrolysis in hExo1, which starts in the precatalytic state with the intact double-strand DNA (dsDNA) recognized and bound (tethered) to the helix-two-turn-helix (H2TH) motif and to a monovalent (K+/Na+) ion state in hExo1 (Figure 1). Then catalysis begins with formation of the assembled active site, where the dsDNA bifurcates into the 5′ and 3′ single strands (i.e., the dsDNA “junction”). At this point, the scissile phosphate of the processed single 5′ strand is properly located at the reactive metal center at the N-terminal domain. In this state, the catalytic residues Lys85 and Arg92 interact with the scissile phosphate, (33) after a rotation (clamped conformation) of the mobile helical arch formed by two α-helices (α4-α5) located near the junction. Here, the side chain of the guide residues Tyr32 and His36 are also rotated (Figure 1). These structural motifs contribute to the “threading mechanism”, whereby the 5′-flap DNA passes through the helical arch. (34) In this way, basic residues steer the phosphate of the 5′ strand, promoting the proper location for hydrolysis of the scissile phosphate on top of the two catalytic ions, as expected for the recognized two-metal-ion mechanism. (34−40)

Figure 1

Figure 1. Catalytic domain of hExo1 in complex with DNA substrate and the two catalytic metal ions (PDB ID 5V06). (Left) hExo1 in cartoon and with colors for different structural motifs. (Right) Closer view of the active site. Three metal ions (MgA, MgB, MgC) are in orange, nucleophilic water molecule (Nu) is in red, two guide residues (Tyr32, His36) in yellow, and residues of the catalytic pocket (Gly2, Asp30, Asp78, Asp152, Asp171, Asp173) are in cyan. Scissile phosphate is correctly positioned for the nucleophilic attack, and MgC is coordinated by the 5′ terminal phosphate.

At this point, hydrolysis of the 5′ recessed end or the 5′-flap DNA substrate occurs. (41−43) In hExo1, this is proposed to be favored by a structured network of interactions involving Arg95, Arg96, Arg121, and Asn124, which are located along the mobile arch. In particular, a key role in phosphate steering is proposed for the Arg96 and Asn124, which interact with the phosphate next to the scissile one (i.e., the terminal 5′ phosphate in the 5′ recessed-end substrate). This is similar to what has been observed in the enzyme hFEN1. (34) After DNA hydrolysis, the nucleotide monophosphate group can leave the active site, with the “free” enzyme that now has Tyr32 and His36 back in their initial conformation.
Remarkably, these structural data show a transient third metal ion that is intermittently located close to the catalytic site during exonuclease catalysis (Scheme 1). This suggests that the transient third ion may play a role in substrate hydrolysis and/or leaving group departure. (44) Indeed, during hExo1 catalysis, four different structures of the assembled active site were solved in the presence of a second-shell and solvent-exposed third metal ion, preserved close to the two-metal ion center (Figure S1). This ion is not found in the structure of the cleaved product, demonstrating its transient nature during catalysis. (1,23)

Scheme 1

Scheme 1. Schematic Representation of the Motion of the Transient Third Metal Ion, Which We Found to Be Intermittently Recruited/Released by Glu89 (switching its inner/outer conformations) during Exonuclease Catalysisa

aThis structural evidence suggests that the transient third ion may be crucial for substrate hydrolysis and/or leaving group departure.

Here, we used force-field-based molecular dynamic (MD) simulations coupled to enhanced sampling free-energy calculations to investigate the role of the third ion during hExo1 catalysis. We compared several systems of wild-type and mutated hExo1 from the reactant state to the product state. We found that the second-shell and conserved Glu89 residue selects, recruits, and places the third ion close to the two-metal-ion catalytic site. We show that this negatively charged residue is functional and conserved among Exo1 belonging to different organisms and that this enzymatic mechanism is likely shared by other nucleases.

Results

Click to copy section linkSection link copied!

Glu89 Selects, Recruits, and Places a Third Ion Close to the Catalytic Site of hExo1

First, we ran multiple unbiased force-field-based molecular dynamic (MD) simulations (∼700 ns in total) of the wild-type (wt) reactant state (Figure S2). We considered the enzyme in the reactive state, replacing the nonreactive Mn ions in the prereactive crystal structure (PDB ID 5V06) (23) with native Mg ions. (31) Notably, this structure features a third Mg ion (MgC) close to the two-metal-ion catalytic site. This ion is bound to the 5′ phosphate of the processed strand.
Residues Tyr32 and His36 are thought to guide substrate binding through interactions with their flexible side chain. (23) Located above the reaction center, these residues maintain the crystallographic conformation, in which their side chain points “down” toward the two catalytic ions MgA and MgB beneath it (Figure S3). These ions maintain an octahedral coordination (45) throughout the simulated time scale (i.e., 360 ns). As a result, the nucleophilic water molecule remains optimally positioned for nucleophilic attack, sitting on top of MgA, in front of the substrate’s scissile phosphodiester bond, at 3.88 ± 0.11 Å (Figure S3).
Despite the overall stability of the protein–DNA complex, the side chains at the base of the mobile arch (i.e., at the gateway) sample different conformations. In particular, the initial clamped conformation shows some flexibility over time, with Lys85 moving slightly further away from the scissile phosphate (∼6 Å vs 4 Å in PDB ID 5V06 see Figure S4). Moreover, we observed an amplified mobility of MgC, reflected in its enhanced RMSD of 1.63 ± 0.42 Å, as compared to the two catalytic ions and its anchor 5′ nucleotide, which are highly stable with an RMSD of 0.54 ± 0.17 and 0.82 ± 0.18 Å, respectively. Indeed, MgC moves from its initial position, forming new interactions with the flexible Glu89 carboxylate along the α-helix of the gateway. As a result, the mobile MgC alternates between a bidentate and a monodentate binding coordination with the 5′ phosphate and Glu89 (Figure S5), always maintaining an octahedral shell. In this regard, the motions of the Glu89 side chain are described by the pseudodihedral angle ϕ (taken along the N, Cα, Cδ, Cγ bonds; see Figure S5), which oscillates from ∼120° to ∼−10° (Figure S5), never reaching a conformation of the product state (−36°, PDB ID 5V0A). (23) Concomitantly, we observed that Arg96, which is located along the mobile arch, shortens its distance from the 5′ phosphate, reaching a value of ∼4.5 ± 0.19 Å (compared to an initial distance of 6 Å in PDB ID 5V06). This distance corresponds well to the value in the crystal structure of the product state (5.1 Å, PDB ID 5V0A) (Figure S6). Taken together, these results support the hypothesis of a gradual removal and possible departure of the third ion during catalysis, as suggested by comparing the structural data of the reactant (MgC present, PDB ID 5V06) and product states (MgC missing, PDB ID 5V0A).
To further characterize the structural impact of MgC bound close to the two-metal-ion catalytic site in the reactant state, we manually removed it and ran multiple simulations of the solvated system (∼1 μs, in total). This protein–DNA complex showed no major difference in the overall backbone stability (see Figure S2) compared to the three-ion reactant state (see above). However, in these replications, the catalytic residues Lys85 and Arg92 maintained their native H-bond pattern with the scissile phosphate (as in the crystallographic structure) for the whole simulation. Thus, the catalytic residue Lys85 behaved differently to the three-ion system. Also, in these simulations, at times Glu89 interacted transiently with Arg93, located along the mobile arch. This interaction was observed only when Glu89 adopted the outer conformations, in the absence of MgC. Concurrently, Arg96 maintained its starting orientation and never interacted with 5′ phosphate (Figure S6). Again, this differs from the observations in the presence of MgC (see previous paragraph).
Intriguingly, during the initial equilibration phase (∼10 ns), the side chain of Glu89 undergoes a marked rearrangement from the initial inner conformation to an outer conformation toward the bulk water. This rotation is captured well by the Glu89 pseudodihedral angle ϕ, which changes from positive values of ∼+100° (inner) to negative values of ∼−95° (outer) (Figure S7). Importantly, in this new solvent-exposed conformation, we observed that Glu89 carboxylate transiently recruits and binds monovalent ions from the bulk (either K+ or Na+, freely diffusing in solution). This result is also confirmed by the radial distribution function, g(r), calculated as the variation of the density of the ions from the center of mass of the 5′ phosphate group (Figure 2), which displays two peaks at ∼2.7 and ∼3.3 Å.

Figure 2

Figure 2. Radial distribution function, g(r), calculated for ions around 10 Å from the center of mass of the 5′ phosphate group. Plot shows the presence of ions ∼3 Å from the 5′ phosphate group for the RS2M system. In this system, a K+ ion approached the negatively charged group. For the RSGlu89Ala and PSGlu89Ala systems, there are no ions within ∼5.5 Å of the 5′ phosphate, as indicated by the g(r) values of ∼0. In the upper right corner, the 5′ phosphate group and Glu89Ala residues are shown in licorice (taken from the PSGlu89Ala simulations).

After metal binding (∼10 ns), Glu89 flips back into its initial inner conformation, carrying the coordinated metal closer to the 5′ phosphate, at ∼3.25 Å. This metal ion is thus brought into a very similar location compared to MgC in the crystal structure of the prereactive state (PDB ID 5V06). After a few hundred nanoseconds (e.g., ∼200 ns for K+ ion), the third metal ion departs spontaneously from the catalytic site. Glu89 then flips again into its outer conformation toward the bulk. These unprompted metal binding and release events, synchronized with the flipping of the Glu89 side chain, were observed multiple times in our extended simulations (Figure S7). This indicates that Glu89 may recruit a third metal ion, bringing it transiently closer to the catalytic site. (46)
To further test Glu89’s role as metal-ion recruiter, we ran multiple MD simulations (∼1 μs in total), inserting the Glu89Ala mutation in the reactant state in the absence of MgC. The overall stability of the enzyme–DNA complex was maintained, with a low RMSD value of 1.27 ± 0.15 Å (see Figure S2). The prereactive state at the active site was also maintained throughout the simulations, that is, the two catalytic metal ions, MgA and MgB, stably maintained their internuclear distance. The nucleophilic water molecule remained properly positioned in front of the scissile phosphodiester bond and the catalytic residues. Finally, Lys85 and Arg92 maintained their initial interaction network with the scissile phosphate. Notably, in this mutated system, we did not see any ion approaching the terminal 5′ phosphate from the bulk solvent. This result is supported by g(r), which confirms that no ion is located within ∼6 Å of the center of mass of the 5′ phosphate group (Figure 2), further suggesting that Glu89 recruits a third metal ion from the bulk.

Third Ion Promotes Leaving Group Departure after DNA Hydrolysis

Here, we used MD simulations (∼1 μs in total) of the products of the wild-type (wt) native state. Thus, we inserted a native aspartate at the Asp225Ala mutation and replaced Mn with native Mg ions in the postreactive crystal structure (PDB ID 5V0A). Notably, at this catalytic stage, the DNA’s processed strand is enzymatically cleaved, with the consequent generation of the leaving group, i.e., the adenosine 5′-monophosphate (AMP) nucleotide, which is now detached from the newly formed 5′ recessed-end substrate. Importantly, the third ion is not present at the catalytic active site in the crystallographic structure.
In the postreactive crystal structure, the Glu89 side chain adopted an intermediate conformation (ϕ of Glu89 was −36°, Figure S1) between the inner (∼+100°) and the outer (∼−95°) conformations. Then during the MD simulations, the Glu89 side chain stably adopted an outer conformation (ϕ of Glu89 becomes ∼−100°). However, after ∼50 ns, we observed the unprompted approach of a third Mg2+ ion from the bulk water (Mgbulk), which came close to the Glu89 side chain. This transient ion thus reached a position close to the 5′-monophosphate of the AMP at a distance of ∼3.3 Å, which was equivalent to that in the prereactive simulations and crystal structure (PDB ID 5V06). In this position, the third metal ion interacted with the 5′ phosphate and Glu89 for the remaining simulation time (Figure 3A). During this event, the two catalytic metal ions moved apart slowly, reaching a distance of ∼5.5 Å (compared to 3.9 Å in the starting model PDB ID 5V0A). The drifting of the internuclear two-metal-ion distance was coupled to a shift in the leaving AMP. This is described well by the collective variable CV1, which measures the distance between the center of mass (COM) of the heavy atoms of AMP and the COM of the Cα of the aspartates in the first coordination shell of the two-metal-ion center (i.e., Asp152, Asp171, Asp173) (Figure S8A). During our simulations, CV1 increased by ∼2 Å, from 9 to 11 Å, reflecting the partial exit of the leaving AMP (Figure S8B and S8C). Moreover, His36, which was initially in the down orientation, immediately rotated into the up conformation (Figure S9), forming a π–π interaction with AMP. This interaction helps the initial displacement of the leaving AMP. Notably, the up conformation of His36 was found in the structure of the enzyme after AMP departure (PDB ID 5V0B), (23) which further suggests the need of this rotation during leaving group release. We also noted a gradual and slight opening in the gateway region at the bottom of the mobile arch, which however maintained an ordered secondary structure (Figure S10). This event is described well by the increase of ∼1 Å of the two distances d1 and d2, which reflect the opening of the α4/α5 interhelix passage (calculated using the Cα of Glu89 and Arg92 along the α4 helix and the Cα of Asn124 and Ile125, located at the bottom of α5 helix; see Figure 4). The probability density function of ϕ, calculated for two states (Mgbulk > 4 Å or Mgbulk < 4 Å from the 5′ phosphate), shows the relative peaks of the two conformations assumed by Glu89, i.e., inner and outer (Figure 3B).

Figure 3

Figure 3. (A) Distance (dMG in yellow) between the third Mg2+ ion, from the bulk (Mgbulk), and the phosphorus of the 5′ phosphate group of AMP. Inset, representation (snapshot from the PS2M simulations) of Mgbulk approaching the terminal 5′ phosphate as well as the pseudo-dihedral angle ϕ of Glu89 side chain (defined by the N–Cα–Cδ−-Cγ atoms). (B) Probability density of the pseudo-dihedral angle ϕ in Glu89 during the simulation: (blue) probability density for dMG values > 4 Å shows the outer conformation as the most populated; (red) probability density for dMG values < 4 Å shows the inner conformation is the most populated.

Figure 4

Figure 4. (A) Graphic representation of d1 and d2 distances (PDB ID 5V0A). (B) Probability density of the distances d1 and d2 calculated during simulations of the systems PS2M (in red), PS3M (in orange), and PSGlu89Ala (in green).

In the products, Arg96 invariantly interacts with the 5′ phosphate in AMP (Figure S6), as also reported for the reactant state simulations in the presence of the MgC. Moreover, after ∼450 ns, we noted a second arginine residue located along the α4 helix, Arg93, which approached the same 5′ phosphate of the AMP, at ∼5 Å. These interactions thus favor the positional shift and partial departure of the AMP from the catalytic site, as indicated by CV1, with consequent destabilization of the two-metal-ion site. We compared these results with an additional postreactive model, which initially contained a third ion at the catalytic site (∼1 μs in total, see SI). These simulations confirmed the enhanced instability of the prone-to-escape leaving group with a partial shift from its starting position during the simulations (Figure S8). Taken together, these results suggest that complete leaving group departure is eventually expected, although this would require longer simulations and would also likely implicate the overtaking of an energetic barrier. (47)
As with the reactant state (see above), we also ran simulations (∼1 μs in total) in the product state of a system with the Glu89Ala mutation in the absence of MgC. The overall stability of the enzyme–DNA complex was maintained (Figure S2). Importantly, in the absence of the Glu89 recruiter, no ion approached the 5′ phosphate of the leaving group, as shown by g(r) (Figure 2). As a result, the leaving group also showed higher stability in its position, as highlighted by the low RMSD value of 1.71 ± 0.18 Å. Moreover, Arg93 almost never interacted with the 5′-phosphate of AMP, and Arg96 stably maintained its initial interaction with AMP. In addition, we did not observe any opening in the gateway region, with the distances d1 and d2 remaining unchanged during the simulations (Figure 4). These results support the hypothesis that Glu89 recruits MgC before (or during) DNA cleavage. In return, MgC seems to promote the release of the leaving group after the chemical step for phosphodiester bond hydrolysis, acting as a shuttle for the AMP departure. (47)

Energetics of the Glu89 Flipping and Leaving Group Departure via Metadynamics Simulations

To sample and determine the semiquantitative energetics of the inner ↔ outer conformational switch of Glu89, we used the pseudodihedral angle ϕ as the collective variable to run multiple metadynamics simulations, with and without the third metal ion at the catalytic site, in the reactant and product states (for a total of ∼920 ns).
In the reactant state with the third ion, Glu89 tended to adopt inner conformations located in an energy minimum at ϕ ≈ 70°, while outer conformations were not visited due to their high energy (Figure 5, red profile). In the inner conformation, MgC stayed close to the reactive center. However, in the absence of the third metal ion, Glu89 could be found in two isoenergetic minima, i.e., inner and outer conformations, separated by a barrier of only ∼3 kcal mol–1 (Figure 5, blue profile). This explains the fact that both Glu89 conformations were similarly populated in our unbiased MD simulations of the system without the third metal ion.

Figure 5

Figure 5. (Bottom) Free energy surface obtained through well-tempered metadynamics simulations for RS3M (red), RS2M (blue), PS2M (green), and PS3M (light purple) systems. Results show three conformations (outer, intermediate, inner). (Top) Graphic representations, taken from PS3M simulaitons, of the three conformations are shown in licorice.

In the product state, regardless of the presence or absence of the third ion, Glu89 visited both the inner and the outer conformations (Figure 5, green and light purple profiles). The conformational switch showed a barrier of ∼4.5 kcal mol–1, with or without MgC. We also located a metastable conformation of Glu89 bound to MgC, at ϕ ≈ −40°, in which the glutamate’s side chain adopted an intermediate orientation between the two minima (inner and outer). Interestingly, this metastable state corresponds well to the crystallographic conformation (ϕ = −36°, PDB ID 5V0A) in which MgC is missing. This is likely because, at this point, Glu89 is already solvent exposed.
Then we evaluated possible pathways and energetics for the release of the leaving group from hExo1 in the presence and absence of MgC. We used confined metadynamics, (48) which enhances the sampling of transversal and often slow degrees of freedom of complex (rare) events, such as the exit of the adenosine monophosphate (AMP) nucleotide from the catalytic site (see SI for further information). Here, the collective variable was CV1 (Figure 6), which captures the degree of departure of the leaving AMP from the reactive site (see definition of CV1, above).

Figure 6

Figure 6. Free energy surface obtained through confined well-tempered metadynamics simulations for PS2M (blue) and PS3M (yellow) systems. (Top) Schematic representation of the CV1, exemplified using a snapshot from PS3M simulations. It represents the distance between the center of mass (COM) of the heavy atoms of the nucleotide leaving group and the COM of the Cα of the aspartates (Asp152, Asp171, Asp173) in the first coordination shell of MgA, MgB. Two different minima, at 8.8 and 12.4 Å, agree with the MD results, in which a partial exit of the leaving group (CV1 ≈ 12.4 Å) was seen only in the presence of MgC (Figure S8).

In the presence of MgC, AMP fell into a minimum at CV1 = ∼12 Å, showing that the leaving group is already shifted out of the catalytic site (CV1 = ∼8.5 Å in the uncleaved prereactive state). This is in line with the plain MD simulations, where it was only in the presence of the third metal ion that the leaving group partially exited from the active site, increasing the CV1 value of ∼2 Å, reaching a value of ∼12 Å (Figure S8). At this point, to allow the full departure of the leaving group, the freed AMP stayed complexed with MgC. In this way, the AMP/MgC complex exited from the catalytic site, passing through the aperture under the mobile arch formed only when MgC is present (due to MgC-mediated AMP drifting out from the catalytic center). At this point, the Glu89 sampled the inner and outer conformations until the leaving group overcame the gateway region. At this point, Glu89 stably adopted the inner conformation, in agreement with the crystallographic structure of the complex after release of the leaving group (PDB ID 5V0B). The physical step for AMP/MgC unbinding occurred with a barrier of ∼16.5 kcal/mol (Figure 6).
In the absence of the third ion, the energy barrier for the overall unbinding process was much higher at ∼35 kcal/mol. Indeed, the system behaved quite differently. In the initial configuration, the leaving AMP fell into an energy minimum where it was poorly solvated, at CV1= ∼9 Å (compared to the case where it was solvated and complexed with MgC, CV1= ∼12 Å, Figure 6). From this state, the exit of the AMP alone had to overcome a first barrier of ∼20 kcal/mol, which is already higher than the barrier in the presence of MgC. Then the exit path showed metastable states (relative energy minima at CV1 = ∼13.5 Å and CV1 = ∼17 Å, Figure 6) where the leaving AMP seemed to be transiently trapped by the formation of short-lived interactions with the enzyme. This may slow the AMP unbinding kinetics. Notably, these transient interactions were not formed for the AMP/MgC leaving complex. It is also worth noticing that in both systems PS3M and PS2M we observed the exit of MgB from the catalytic pocket, which occurred concertedly with the exit of the AMP leaving (see Movie in SI). In detail, the MgB catalytic ion remained coordinated to the oxygen of the OH in C3 position of the sugar of the AMP leaving group. Interestingly, the concerted departure of MgB, MgC, and the leaving group AMP agrees with previous studies showing that MgB dissociates from the catalytic site ∼200 times faster compared MgA. (49) Thus, the concerted exit of both MgB and MgC is found here to promote AMP departure by stabilizing the newly formed negative charge on the leaving group after substrate hydrolysis. (47)

Discussion

Click to copy section linkSection link copied!

Recently, a time series of structural intermediates captured during human exonuclease1 (hExo1) catalysis has revealed the presence of a third metal ion (MgC) close to the active site. (23) Intriguingly, a transient metal ion was recently observed in a few DNA/RNA-processing enzymes. (1,2,4,50) The role of this third additional metal ion at the catalytic center is still unclear. (44) Here, we used force-field-based molecular dynamics (MD) simulations and free-energy calculations to investigate the recruiting mechanism and functional role of a third metal for hExo1 catalysis. We simulated and compared several model systems built with recent hExo1 structures of the wild-type (wt) hExo1/DNA complex in the reactant and product states with and without MgC.
During our multiple and extended MD simulations (∼6 μs in total) we first observed that Glu89 sometimes oscillated but mostly maintained its starting conformation. In this inner conformation, the Glu89 carboxylate group points toward the 5′ phosphate. However, at times and only in the absence of MgC the Glu89 carboxylate group switched its orientation, adopting outer conformations that pointed toward the bulk solvent. Free-energy calculations confirmed that in the reactant state and in the presence of MgC Glu89 tended to adopt inner conformations, which are located in an energy minimum at ϕ ≈ 70°. Outer conformations are not visited due to their high energy. However, in the absence of MgC, the inner and outer conformations become isoenergetic, with a barrier of only ∼3 kcal/mol in between.
Importantly, following this conformational switch in equilibrium MD simulations in the absence of MgC, we observed that transient monovalent ions were freely recruited from the bulk by the outer conformation of Glu89. This conformation therefore seems to act as an anchor point for (third)metal–enzyme complexation. Then Glu89 could switch back and adopt the inner conformation, bringing the bound metal ion (either K+ or Na+, from these simulations) close to the terminal 5′ phosphate (∼3.5 Å) (Figure S7). This third metal was spontaneously recruited and located in the same position as the third ion captured in the prereactive crystal structure (PDB ID 5V06). This Glu89-mediated mechanism for metal recruitment was further validated by simulations of mutated Glu89Ala systems. These simulations confirmed that in the absence of Glu89 no metal ion from the bulk was spontaneously recruited close to the catalytic center. Interestingly, the role of Glu89 in hExo1 is similar to the role previously proposed for Glu188 in Bacillus halodurans ribonuclease H (BhRNase H), where MD simulations suggested that this residue attracted a transient third ion. (46) Intriguingly, a transient third solvent-exposed cation was found close to the two-metal-ion active site of D. mobilis homing endonuclease, I-DmoI. (51)
In the unbiased MD simulations of the product state we observed the unprompted entry of the third metal ion MgC from the bulk, reconstituting the three-metal-ion system (Figure 3). This happened concomitantly to the rotation of the Glu89 side chain from outer to inner, thus destabilizing the geometry of the catalytic active site. Notably, during this process with the double-strand DNA already bound to hExo1, the enzyme maintained an ordered secondary structure of the mobile arch. After the initial DNA binding and subsequent DNA hydrolysis, the enzyme will have to position the “next” scissile bond into the active site. An ordered structure may thus favor the processivity of the exonuclease activity of this enzyme. At this stage the system evolved toward the final catalytic step, i.e., the exit of the leaving group from the catalytic site. In this system, Glu89 was free to populate the inner and outer conformations, overcoming an energy barrier of ∼4.5 kcal/mol, calculated from metadynamics simulations. In addition, we computed the energetics for the full release of the leaving group in the presence or absence of MgC using confined well-tempered metadynamic simulations. (48) The energetic barrier for AMP departure was ∼16.5 kcal/mol in the presence of MgC and ∼35 kcal/mol in its absence.
The Gibbs free energy (ΔG) for the overall catalytic process of hExo1 is 19.6 kcal/mol, computed using the experimental kcat for hExo1 (see SI for further information). (31) This energy value corresponds fairly well to our estimation of the free-energy barrier for the unbinding process of the leaving group in the presence of the third ion, i.e., ∼16.5 kcal/mol. The leaving group departure may therefore be rate limiting for the exonuclease catalytic process in hExo1, as already proposed for other metallonucleases (e.g., FENs, APE1, PvuII, MunI, NaeI, SfiI, EcoRI, EcoRV). (52−59)
These results suggest a mechanism where Glu89 recruits a third metal ion in the reactant state. Clearly, quantum calculations are needed to evaluate the mechanistic implications of this additional ion for the chemical step of DNA hydrolysis. (60−63) However, from these classical MD simulations, it emerges that the third ion promotes leaving group departure, acting as a shuttle for the exit of the nucleotide monophosphate product from the catalytic site. Notably, this result is in line with evidence of a third-ion-mediated leaving mechanism for pyrophosphate departure in polymerase enzymes. (47)
To further test this mechanistic hypothesis and investigate whether this enzymatic strategy is shared by other nucleases, we performed sequence alignments via the Needleman–Wunsch algorithm (64) using 10 different eukaryotic species of Exo1 (see SI for more information). We found that the Glu89 is fully conserved among these enzymes, as for those residues forming the reaction center, and second-shell residues like Lys85 and Arg92, and the guide residues Tyr32 and His36 (Figure S11). This suggests that Glu89 is an integral part of the enzymatic machinery for efficient catalysis in Exo1.
We also performed structural comparisons using recent crystallographic structures of additional nucleases and identified a shared spatial localization in these enzymes of an acidic residue (Glu/Asp), in analogy to Glu89 in hExo1. In such enzymes, in fact, we always identified the presence of a Glu/Asp residue located in a second-shell sphere cantered on the two-metal-ion active site. This acidic residue is always situated in a solvent-accessible position (thus able to recruit ions from the bulk), being strategically located on the side of the expected exit path for leaving group departure, with respect to the catalytic center. For example, human ExoG, (65,66) which is 5′ metallo-exonuclease enzyme cocrystallized in complex with the DNA substrate, has a glutamate (Glu317) residue located near the terminal 5′ phosphate. Here, Glu317 resides in a solvent-exposed area. Notably, Glu317 can assume different orientations in the available crystals (see PDB ID 5T5C vs PDB ID 5T40), (65) which suggests that this glutamate may act as a recruiter of metal ions in the same way as for Glu89 in hExo1 (Figure 7A). Another case is the human λ-Exonuclease, (67,68) where Glu36 is located close to the 5′ phosphate of the DNA substrate (Figure 7B). Here, too, it has been hypothesized that a third metal ion may transiently bind close to the two-metal-ion site, likely aiding the leaving group departure. (49) A further example is RecJ nuclease, (69) where Asp158 is solvent exposed and close to the active site in a similar position as Glu89 in hExo1 (Figure 7C).

Figure 7

Figure 7. Close views of the active site of 5′ metallonuclease members that possess an analogous acid residue (light green) close to the two-metal-ion center (MA, MB, in orange), the active site (in cyan), and the leaving group (indicated by a dashed line). (A) Human ExoG, in which Glu317 is pointing in the inner (PDB ID 5T5C) and (A′) outer conformations (PDB ID 5T40, merged with the DNA substrate from PDB ID 5T5C). (B) Escherichia phage T5Fen (PDB ID 5HNK). (C) Human λ-Exonuclease X (PDB ID 3SM4). (D) D. radiodurans RecJ (PDB ID 5F55). (E) Sequence alignment of M. smegmatis FenA and E. phage T5Fen. Conserved acid residue (Glu/Asp) is indicated in orange.

Then we looked at the hExo1 family member bacteriophage T5 flap endonuclease (T5 Fen, PDB ID 5HNK). (16) We identified Glu83, which is located along the mobile arch, in analogy to Glu89 in hExo1. Finally, we considered the recent high-resolution X-ray structure of M. smegmatis FenA, (17) which is a 5′ structure-specific nuclease and close homologue of phage T5Fen, and its sequence alignment with T5Fen. It suggests that Glu83 in T5Fen may correspond to Asp85 in FenA (Figure 7D and 7E). Notably, holo forms of T5Fen and FenA were crystallized in complex with three metal ions in the active site. (16,17) This further corroborates the idea that multimetal-ion catalytic sites may be necessary for nucleic-acid processing in these enzymes. Indeed, while the presence of a third metal ion, located in the vicinity of the two-metal-ion active site, is a novel aspect in polymerases and nucleases, the exact position of such additional metal ion with respect to the catalytic center can vary. (4) For example, T5FEN and FenA enzymes were solved with a third ion in a different relative position, although always in close proximity of the reactive two-metal-ion center. It is thus plausible that the third transient ion may play different roles during catalysis, according to its specific location at the catalytic center.
Another intriguing aspect is the recurring presence of second-shell positively charged residues that surround the metal-aided catalytic site in nucleic-acid-processing enzyme. (3) During our extended MD simulations, we observed the interaction of Arg93, along the helical arch, with the terminal phosphate of the substrate. We structurally aligned hExo1 (PDB ID 5V0E) (23) with hFEN1 (PDB ID 5KSE), (34) and we noted that the terminal guanidine, amide, and amine groups of Arg93, Asn124 (in hExo1), and Lys132 (in hFEN1) residues were all within a sphere of ∼3 Å and close to the 5′ phosphate in the leaving group (Figure S12). This result further suggests that Arg93, together with Arg96, may have a crucial role in 5′ phosphate steering. (34) In this respect, Arg93 is an important point of control in Exo1 poly(ADP-ribose) binding, as proved by in vitro and in vivo assays on natural Arg93Gly mutation. (70) Moreover, a common feature of different nucleic-acid-processing enzymes is a solvent-exposed and positively charged residue that interacts with the negatively charged moiety in the leaving group (e.g., Arg96–5′ phosphate). (71) Interestingly, this Arg96–5′ phosphate interaction comes in addition to an already complex architecture characterized by a number of positively charged residues surrounding the active site. Indeed, hExo1 is one of a large set of nucleic-acid-processing enzymes characterized by the recurring presence of positively charged elements in the vicinity of the reactive site. (3) These elements in hExo1 comprise the catalytic Lys85 and Arg92. These positively charged residues are thought to play a role in the phosphate steering during the threading mechanism. These residues therefore seem crucial for substrate recognition, binding, catalysis, translocation, and initial product release. However, the Glu89-mediated recruitment and binding of a third ion from the bulk also seems necessary for the full departure of the leaving group from the enzyme’s catalytic site in hExo1 and likely in several other nuclease enzymes.

Conclusions

Click to copy section linkSection link copied!

Our results provide new insights into the functional role of a third metal ion, which was recently found transiently located at the catalytic site of nuclease enzymes during catalysis. Using molecular dynamics and free-energy simulations applied to multiple systems, we considered the conformational switch of the side chain of a specific residue, Glu89, which is located near the active site in hExo1. We noted that this conformational switch favors the recruitment of a third metal ion from the bulk. The third metal ion is thus promptly positioned near the catalytic center, in accordance with the structural evidence. Our simulations also indicate that this ion serves as an exit shuttle for the leaving group departure from the catalytic site after DNA hydrolysis. The exit mechanism is also favored by the initial involvement of positively charged residues, which are located in an extended and highly structured second-shell area at the two-metal-ion active site. (3,72,73) Finally, our structural analyses of nuclease enzymes show that such a negatively charged residue (Glu/Asp) is persistently found in a similar, structurally conserved, and strategic position in several other 5′ structure-specific nucleases, which seem to share this enzymatic mechanism to promote DNA hydrolysis. These findings may have an implication for de novo enzyme engineering and structure-based drug design. (74,75)

Methods

Click to copy section linkSection link copied!

Structural Models

We used six different model systems: (i) the wild-type (wt) reactant state (RS3M), based on the recent time-resolved X-ray structure of the complex of the hExo1(PDB ID 5V06), which includes the 5′ recessed-end DNA substrate (in this structure, the third metal ion is located close to the catalytic site); (ii) the same wt reactant state with the third metal ion manually removed (RS2M); (iii) the mutated Glu89Ala reactant state (RSGlu89Ala), modeled on the same X-ray structure (PDB ID 5V06); (iv) the product state (PS2M), based on the recent time-resolved X-ray structure of the ternary complex of hExo1(PDB ID 5V0A), characterized by the newly formed 5′ recessed-end DNA substrate and the leaving adenosine monophosphate (AMP); (v) the product state with the third metal ion close to the 5′ phosphate of the AMP (PS3M) (this system is modeled on the X-ray structure of the prereactive complex, in which we manually cleaved the scissile phosphate); (vi) the mutated Glu89Ala product state (PSGlu89Ala), modeled on the same X-ray structure (PDB ID 5V0A). (23) Protein coordinates and relevant trajectory files are available from the authors upon request.

Classical Molecular Dynamics Simulations

To investigate the functional dynamics of the hExo1/DNA complex, we used extensive force-field-based MD simulations, which are highly informative for complex enzyme/nucleic acid assemblies. (76−80) Here, the AMBER/ff14SB (81) and OL15 (82,83) force fields were used to treat the hExo1 enzyme and the DNA, respectively. The terminal 5′ thymine monophosphate, the 5′ adenosine monophosphate, and the terminal 5′ guanidine monophosphate were treated with the general Amber force field (GAFF). (84) The atomic charges were derived by fitting the electrostatic potential according to the Merz–Singh–Kollman scheme, (85) the RESP fitting procedure. (86) The length of all bonds involving hydrogen atoms was constrained to the equilibrium using the P-LINCS algorithm, (87) and a time integration step of 2 fs was used. All simulations were performed using GROMACS 5.1 code. (88) Long-range electrostatic interactions were calculated with the particle mesh Ewald method with a Fourier grid spacing of 1.6 Å. (89,90) Periodic boundary conditions in the three directions of Cartesian space were applied. The magnesium ions were treated with a nonbonded approach based on the “atoms in molecules” theory partitioning scheme. (91,92) The systems were solvated with TIP3P water molecules (93) and neutralized adding Mg2+, Na+, K+, and Cl ions, as indicated in the crystallization procedure. (23) The total number of atoms was ∼60 000 for each system (see SI for more information). We followed a two-step procedure for the MD simulations: first, the equilibration phase in which we followed different procedures depending on the system (see SI). Then the production phase was carried out in an NPT ensemble, a constant temperature of 310 K imposed using the velocity-rescaling thermostat, (94) and a constant pressure of 1 bar maintained with a Parrinello–Rahman barostat. (95) We collected MD simulations of ∼0.7 μs for RS3M and ∼1 μs for each of RS2M, RSGlu89Ala, PS2M, PS3M, and PSGlu89Ala for a total of ∼6 μs of MD.

Free-Energy Calculations

We used well-tempered metadynamics (96) to characterize and estimate the free-energy landscape of Glu89 conformational flexibility. We selected a collective variable (CV) that distinguished between the inner, the outer, and the intermediate conformations adopted by Glu89 during the MD simulations. Thus, the selected CV was the pseudodihedral angle ϕ defined by the N, Cα, Cδ, and Cγ atoms on Glu89 (see SI). In particular, based on our MD simulations, inner conformations are adopted at ∼70° < ϕ < ∼100°, intermediate conformations at ∼−40° < ϕ < ∼−10°, and outer conformations at ∼−150° < ϕ < ∼−100°. We performed well-tempered metadynamics by biasing the CV using an initial hill height of 0.02 kcal mol–1, a hill width of 0.35 rad, a fictitious CV temperature of 1550 K, and a deposition rate of 1 ps–1. The simulations were conducted until convergence (see SI for more information).
Well-tempered metadynamics was also used to evaluate possible pathways and the semiquantitative energetics for the release of the leaving group (i.e., adenosine 5′-monophosphate, AMP) from the active site. We used a confined metadynamics approach, (48) which excludes regions of the conformational space that are not relevant to the chemical event under investigation. The selected CV was the distance between the center of mass (COM) of the heavy atoms of AMP and the COM of the Cα of the aspartates (Asp152, Asp171, Asp173) in the first coordination shell of the two catalytic metal ions (see SI). This CV indicates the degree of departure of AMP from the active site in our MD simulations. We used an initial hill height of 0.29 kcal mol–1, a hill width of 0.6 Å, a fictitious CV temperature of 3720 K, and a deposition rate of 1 ps–1. The simulations were conducted until convergence (see SI for more information).

Supporting Information

Click to copy section linkSection link copied!

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/jacs.9b10656.

  • Setup protocols of the force-field-based MD simulations; details on enhanced sampling free energy calculations and additional analyses; PDB analyses, structure, and sequence alignments (PDF)

  • Illustrative movie of the overall catalytic process of hExo1, created by merging PDB structures with fragments of trajectories of equilibrium and metadynamics simulations (MP4)

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

Click to copy section linkSection link copied!

  • Corresponding Author
  • Authors
  • Notes
    The authors declare no competing financial interest.

Acknowledgments

Click to copy section linkSection link copied!

M.D.V. thanks the Italian Association for Cancer Research (AIRC) for financial support (IG 23679). V.G. thanks the European Molecular Biology Organization (EMBO) for financial support (ALTF 103-2018). We thank Giuseppina La Sala for helpful discussions. We thank Grace Fox for her proofreading and copyediting.

References

Click to copy section linkSection link copied!

This article references 96 other publications.

  1. 1
    Yang, W.; Weng, P. J.; Gao, Y. A New Paradigm of DNA Synthesis: Three-Metal-Ion Catalysis. Cell Biosci. 2016, 6 (1), 51,  DOI: 10.1186/s13578-016-0118-2
  2. 2
    Gao, Y.; Yang, W. Capture of a Third Mg2+ Is Essential for Catalyzing DNA Synthesis. Science 2016, 352 (6291), 13341337,  DOI: 10.1126/science.aad9633
  3. 3
    Genna, V.; Colombo, M.; De Vivo, M.; Marcia, M. Second-Shell Basic Residues Expand the Two-Metal-Ion Architecture of DNA and RNA Processing Enzymes. Structure 2018, 26 (1), 4050,  DOI: 10.1016/j.str.2017.11.008
  4. 4
    Genna, V.; Donati, E.; De Vivo, M. The Catalytic Mechanism of DNA and RNA Polymerases. ACS Catal. 2018, 8 (12), 1110311118,  DOI: 10.1021/acscatal.8b03363
  5. 5
    Pavlov, Y. I.; Shcherbakova, P. V.; Rogozin, I. B. Roles of DNA Polymerases in Replication, Repair, and Recombination in Eukaryotes. Int. Rev. Cytol. 2006, 255, 41132,  DOI: 10.1016/S0074-7696(06)55002-8
  6. 6
    Nishino, T.; Morikawa, K. Structure and Function of Nucleases in DNA Repair: Shape, Grip and Blade of the DNA Scissors. Oncogene 2002, 21 (58), 90229032,  DOI: 10.1038/sj.onc.1206135
  7. 7
    Marti, T. M.; Fleck, O. DNA Repair Nucleases. Cell. Mol. Life Sci. 2004, 61 (3), 336354,  DOI: 10.1007/s00018-003-3223-4
  8. 8
    Potapov, V.; Fu, X.; Dai, N.; Corrêa, I. R.; Tanner, N. A.; Ong, J. L. Base Modifications Affecting RNA Polymerase and Reverse Transcriptase Fidelity. Nucleic Acids Res. 2018, 46 (11), 57535763,  DOI: 10.1093/nar/gky341
  9. 9
    Nimonkar, A. V.; Ozsoy, A. Z.; Genschel, J.; Modrich, P.; Kowalczykowski, S. C. Human Exonuclease 1 and BLM Helicase Interact to Resect DNA and Initiate DNA Repair. Proc. Natl. Acad. Sci. U. S. A. 2008, 105 (44), 1690616911,  DOI: 10.1073/pnas.0809380105
  10. 10
    Loeb, L. A.; Monnat, R. J. DNA Polymerases and Human Disease. Nat. Rev. Genet. 2008, 9 (8), 594604,  DOI: 10.1038/nrg2345
  11. 11
    Zheng, L.; Jia, J.; Finger, L. D.; Guo, Z.; Zer, C.; Shen, B. Functional Regulation of FEN1 Nuclease and Its Link to Cancer. Nucleic Acids Res. 2011, 39 (3), 781794,  DOI: 10.1093/nar/gkq884
  12. 12
    Peltomäki, P. Role of DNA Mismatch Repair Defects in the Pathogenesis of Human Cancer. J. Clin. Oncol. 2003, 21 (6), 11741179,  DOI: 10.1200/JCO.2003.04.060
  13. 13
    Dai, Y.; Tang, Z.; Yang, Z.; Zhang, L.; Deng, Q.; Zhang, X.; Yu, Y.; Liu, X.; Zhu, J. EXO1 Overexpression Is Associated with Poor Prognosis of Hepatocellular Carcinoma Patients. Cell Cycle 2018, 17 (19–20), 23862397,  DOI: 10.1080/15384101.2018.1534511
  14. 14
    Ivanov, I.; Tainer, J. A.; McCammon, J. A. Unraveling the Three-Metal-Ion Catalytic Mechanism of the DNA Repair Enzyme Endonuclease IV. Proc. Natl. Acad. Sci. U. S. A. 2007, 104 (5), 14651470,  DOI: 10.1073/pnas.0603468104
  15. 15
    Prieto, J.; Redondo, P.; Merino, N.; Villate, M.; Montoya, G.; Blanco, F. J.; Molina, R. Structure of the I-SceI Nuclease Complexed with Its DsDNA Target and Three Catalytic Metal Ions. Acta Crystallogr., Sect. F: Struct. Biol. Commun. 2016, 72 (6), 473479,  DOI: 10.1107/S2053230X16007512
  16. 16
    AlMalki, F. A.; Flemming, C. S.; Zhang, J.; Feng, M.; Sedelnikova, S. E.; Ceska, T.; Rafferty, J. B.; Sayers, J. R.; Artymiuk, P. J. Direct Observation of DNA Threading in Flap Endonuclease Complexes. Nat. Struct. Mol. Biol. 2016, 23 (7), 640646,  DOI: 10.1038/nsmb.3241
  17. 17
    Uson, M. L.; Carl, A.; Goldgur, Y.; Shuman, S. Crystal Structure and Mutational Analysis of Mycobacterium Smegmatis FenA Highlight Active Site Amino Acids and Three Metal Ions Essential for Flap Endonuclease and 5 Exonuclease Activities. Nucleic Acids Res. 2018, 46 (8), 41644175,  DOI: 10.1093/nar/gky238
  18. 18
    Yang, W. Nucleases: Diversity of Structure, Function and Mechanism. Q. Rev. Biophys. 2011, 44 (1), 193,  DOI: 10.1017/S0033583510000181
  19. 19
    Palermo, G.; Cavalli, A.; Klein, M. L.; Alfonso-Prieto, M.; Dal Peraro, M.; De Vivo, M. Catalytic Metal Ions and Enzymatic Processing of DNA and RNA. Acc. Chem. Res. 2015, 48 (2), 220228,  DOI: 10.1021/ar500314j
  20. 20
    Palermo, G.; Stenta, M.; Cavalli, A.; Dal Peraro, M.; De Vivo, M. Molecular Simulations Highlight the Role of Metals in Catalysis and Inhibition of Type II Topoisomerase. J. Chem. Theory Comput. 2013, 9 (2), 857862,  DOI: 10.1021/ct300691u
  21. 21
    Schmidt, B. H.; Burgin, A. B.; Deweese, J. E.; Osheroff, N.; Berger, J. M. A Novel and Unified Two-Metal Mechanism for DNA Cleavage by Type II and IA Topoisomerases. Nature 2010, 465 (7298), 641644,  DOI: 10.1038/nature08974
  22. 22
    Perera, L.; Freudenthal, B. D.; Beard, W. A.; Shock, D. D.; Pedersen, L. G.; Wilson, S. H. Requirement for Transient Metal Ions Revealed through Computational Analysis for DNA Polymerase Going in Reverse. Proc. Natl. Acad. Sci. U. S. A. 2015, 112 (38), E5228E5236,  DOI: 10.1073/pnas.1511207112
  23. 23
    Shi, Y.; Hellinga, H. W.; Beese, L. S. Interplay of Catalysis, Fidelity, Threading, and Processivity in the Exo- and Endonucleolytic Reactions of Human Exonuclease I. Proc. Natl. Acad. Sci. U. S. A. 2017, 114 (23), 60106015,  DOI: 10.1073/pnas.1704845114
  24. 24
    Jimeno, S.; Herrera-Moyano, E.; Ortega, P.; Aguilera, A. Differential Effect of the Overexpression of Rad2/XPG Family Endonucleases on Genome Integrity in Yeast and Human Cells. DNA Repair 2017, 57, 6675,  DOI: 10.1016/j.dnarep.2017.06.030
  25. 25
    Emmert, S.; Schneider, T. D.; Khan, S. G.; Kraemer, K. H. The Human XPG Gene: Gene Architecture, Alternative Splicing and Single Nucleotide Polymorphisms. Nucleic Acids Res. 2001, 29 (7), 14431452,  DOI: 10.1093/nar/29.7.1443
  26. 26
    Qiu, J.; Qian, Y.; Chen, V.; Guan, M.-X.; Shen, B. Human Exonuclease 1 Functionally Complements Its Yeast Homologues in DNA Recombination, RNA Primer Removal, and Mutation Avoidance. J. Biol. Chem. 1999, 274 (25), 1789317900,  DOI: 10.1074/jbc.274.25.17893
  27. 27
    Genschel, J.; Bazemore, L. R.; Modrich, P. Human Exonuclease I Is Required for 5′ and 3′ Mismatch Repair. J. Biol. Chem. 2002, 277 (15), 1330213311,  DOI: 10.1074/jbc.M111854200
  28. 28
    Wei, K.; Clark, A. B.; Wong, E.; Kane, M. F.; Mazur, D. J.; Parris, T.; Kolas, N. K.; Russell, R.; Hou, H.; Kneitz, B.; Yang, G.; Kunkel, T. A.; Kolodner, R. D.; Cohen, P. E.; Edelmann, W. Inactivation of Exonuclease I in Mice Results in DNA Mismatch Repair Defects, Increased Cancer Susceptibility, and Male and Female Sterility. Genes Dev. 2003, 17 (5), 603614,  DOI: 10.1101/gad.1060603
  29. 29
    Zhu, Z.; Chung, W. H.; Shim, E. Y.; Lee, S. E.; Ira, G. Sgs1 Helicase and Two Nucleases Dna2 and Exo1 Resect DNA Double-Strand Break Ends. Cell 2008, 134 (6), 981994,  DOI: 10.1016/j.cell.2008.08.037
  30. 30
    Vallur, A. C.; Maizels, N. Complementary Roles for Exonuclease 1 and Flap Endonuclease 1 in Maintenance of Triplet Repeats. J. Biol. Chem. 2010, 285 (37), 2851428519,  DOI: 10.1074/jbc.M110.132738
  31. 31
    Lee, B. I.; Wilson, D. M. The RAD2 Domain of Human Exonuclease 1 Exhibits 5′ to 3′ Exonuclease and Flap Structure-Specific Endonuclease Activities. J. Biol. Chem. 1999, 274 (53), 3776337769,  DOI: 10.1074/jbc.274.53.37763
  32. 32
    Keijzers, G.; Bohr, V. A.; Rasmussen, L. J. Human Exonuclease 1 (EXO1) Activity Characterization and Its Function on FLAP Structures. Biosci. Rep. 2015, 35 (3), e00206  DOI: 10.1042/BSR20150058
  33. 33
    Orans, J.; McSweeney, E. A.; Iyer, R. R.; Hast, M. A.; Hellinga, H. W.; Modrich, P.; Beese, L. S. Structures of Human Exonuclease 1 DNA Complexes Suggest a Unified Mechanism for Nuclease Family. Cell 2011, 145 (2), 212223,  DOI: 10.1016/j.cell.2011.03.005
  34. 34
    Tsutakawa, S. E.; Thompson, M. J.; Arvai, A. S.; Neil, A. J.; Shaw, S. J.; Algasaier, S. I.; Kim, J. C.; Finger, L. D.; Jardine, E.; Gotham, V. J. B.; Sarker, A. H.; Her, M. Z.; Rashid, F.; Hamdan, S. M.; Mirkin, S. M.; Grasby, J. A.; Tainer, J. A. Phosphate Steering by Flap Endonuclease 1 Promotes 5′-Flap Specificity and Incision to Prevent Genome Instability. Nat. Commun. 2017, 8, 15855,  DOI: 10.1038/ncomms15855
  35. 35
    Steitz, T. A.; Steitz, J. A. A General Two-Metal-Ion Mechanism for Catalytic RNA. Proc. Natl. Acad. Sci. U. S. A. 1993, 90 (14), 64986502,  DOI: 10.1073/pnas.90.14.6498
  36. 36
    Yang, W.; Lee, J. Y.; Nowotny, M. Making and Breaking Nucleic Acids: Two-Mg2+-Ion Catalysis and Substrate Specificity. Mol. Cell 2006, 22 (1), 513,  DOI: 10.1016/j.molcel.2006.03.013
  37. 37
    Dupureur, C. M. One Is Enough: Insights into the Two-Metal Ion Nuclease Mechanism from Global Analysis and Computational Studies. Metallomics. 2010, 2 (9), 609620,  DOI: 10.1039/c0mt00013b
  38. 38
    Beese, L. S.; Steitz, T. A. Structural basis for the 3'-5' exonuclease activity of Escherichia coli DNA polymerase 1: a two metal ion mechanism. EMBO J. 1991, 10 (1), 2533,  DOI: 10.1002/j.1460-2075.1991.tb07917.x
  39. 39
    De Vivo, M.; Dal Peraro, M.; Klein, M. L. Phosphodiester Cleavage in Ribonuclease H Occurs via an Associative Two-Metal-Aided Catalytic Mechanism. J. Am. Chem. Soc. 2008, 130 (33), 1095510962,  DOI: 10.1021/ja8005786
  40. 40
    Genna, V.; Vidossich, P.; Ippoliti, E.; Carloni, P.; De Vivo, M. A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases. J. Am. Chem. Soc. 2016, 138 (44), 1459214598,  DOI: 10.1021/jacs.6b05475
  41. 41
    Shaw, S. J.; Finger, L. D.; Grasby, J. A. Human Exonuclease 1 Threads 5′-Flap Substrates through Its Helical Arch. Biochemistry 2017, 56 (29), 37043707,  DOI: 10.1021/acs.biochem.7b00507
  42. 42
    Tomlinson, C. G.; Atack, J. M.; Chapados, B.; Tainer, J. A.; Grasby, J. A. Substrate Recognition and Catalysis by Flap Endonucleases and Related Enzymes. Biochem. Soc. Trans. 2010, 38 (2), 433437,  DOI: 10.1042/BST0380433
  43. 43
    Bennet, I. A; Finger, L D.; Baxter, N. J; Ambrose, B.; Hounslow, A. M; Thompson, M. J; Exell, J. C; Shahari, N. N. B M.; Craggs, T. D; Waltho, J. P; Grasby, J. A Regional Conformational Flexibility Couples Substrate Specificity and Scissile Phosphate Diester Selectivity in Human Flap Endonuclease 1. Nucleic Acids Res. 2018, 46 (11), 56185633,  DOI: 10.1093/nar/gky293
  44. 44
    Raper, A. T.; Reed, A. J.; Suo, Z. Kinetic Mechanism of DNA Polymerases: Contributions of Conformational Dynamics and a Third Divalent Metal Ion. Chem. Rev. 2018, 118 (2), 60006025,  DOI: 10.1021/acs.chemrev.7b00685
  45. 45
    Black, C. B.; Huang, H. W.; Cowan, J. A. Biological Coordination Chemistry of Magnesium, Sodium, and Potassium Ions. Protein and Nucleotide Binding Sites. Coord. Chem. Rev. 1994, 135–136 (C), 165202,  DOI: 10.1016/0010-8545(94)80068-5
  46. 46
    Ho, M.-H.; De Vivo, M.; Dal Peraro, M.; Klein, M. L. Understanding the Effect of Magnesium Ion Concentration on the Catalytic Activity of Ribonuclease H through Computation: Does a Third Metal Binding Site Modulate Endonuclease Activity?. J. Am. Chem. Soc. 2010, 132 (39), 1370213712,  DOI: 10.1021/ja102933y
  47. 47
    Genna, V.; Gaspari, R.; Dal Peraro, M.; De Vivo, M. Cooperative Motion of a Key Positively Charged Residue and Metal Ions for DNA Replication Catalyzed by Human DNA Polymerase-η. Nucleic Acids Res. 2016, 44 (6), 28272836,  DOI: 10.1093/nar/gkw128
  48. 48
    La Sala, G.; Riccardi, L.; Gaspari, R.; Cavalli, A.; Hantschel, O.; De Vivo, M. HRD Motif as the Central Hub of the Signaling Network for Activation Loop Autophosphorylation in Abl Kinase. J. Chem. Theory Comput. 2016, 12 (11), 55635574,  DOI: 10.1021/acs.jctc.6b00600
  49. 49
    Hwang, W.; Yoo, J.; Lee, Y.; Park, S.; Hoang, P. L.; Cho, H.; Yu, J.; Hoa Vo, T. M.; Shin, M.; Jin, M. S.; Park, D.; Hyeon, C.; Lee, G. Dynamic Coordination of Two-Metal-Ions Orchestrates λ-Exonuclease Catalysis. Nat. Commun. 2018, 9 (1), 4404,  DOI: 10.1038/s41467-018-06750-9
  50. 50
    Samara, N. L.; Yang, W. Cation Trafficking Propels RNA Hydrolysis. Nat. Struct. Mol. Biol. 2018, 25 (8), 715721,  DOI: 10.1038/s41594-018-0099-4
  51. 51
    Molina, R.; Stella, S.; Redondo, P.; Gomez, H.; Marcaida, M. J.; Orozco, M.; Prieto, J.; Montoya, G. Visualizing Phosphodiester-Bond Hydrolysis by an Endonuclease. Nat. Struct. Mol. Biol. 2015, 22 (1), 6572,  DOI: 10.1038/nsmb.2932
  52. 52
    Sasnauskas, G.; Jeltsch, A.; Pingoud, A.; Siksnys, V. Plasmid DNA Cleavage by MunI Restriction Enzyme: Single-Turnover and Steady-State Kinetic Analysis. Biochemistry 1999, 38 (13), 40284036,  DOI: 10.1021/bi982456n
  53. 53
    Xie, F.; Dupureur, C. M. Kinetic Analysis of Product Release and Metal Ions in a Metallonuclease. Arch. Biochem. Biophys. 2009, 483 (1), 19,  DOI: 10.1016/j.abb.2009.01.001
  54. 54
    Sengerová, B.; Tomlinson, C.; Atack, J. M.; Williams, R.; Sayers, J. R.; Williams, N. H.; Grasby, J. A. Brønsted Analysis and Rate-Limiting Steps for the T5 Flap Endonuclease Catalyzed Hydrolysis of Exonucleolytic Substrates. Biochemistry 2010, 49 (37), 80858093,  DOI: 10.1021/bi100895j
  55. 55
    Taylor, J. D.; Halford, S. E. Discrimination between DNA Sequences by the EcoRV Restriction Endonuclease. Biochemistry 1989, 28 (15), 61986207,  DOI: 10.1021/bi00441a011
  56. 56
    Yang, C. C.; Baxter, B. K.; Topal, M. D. DNA Cleavage by NaeI: Protein Purification, Rate-Limiting Step, and Accuracy. Biochemistry 1994, 33 (49), 1491814925,  DOI: 10.1021/bi00253a031
  57. 57
    Williams, R.; Sengerova, B.; Osborne, S.; Syson, K.; Ault, S.; Kilgour, A.; Chapados, B. R.; Tainer, J. A.; Sayers, J. R.; Grasby, J. A. Comparison of the Catalytic Parameters and Reaction Specificities of a Phage and an Archaeal Flap Endonuclease. J. Mol. Biol. 2007, 371 (1), 3448,  DOI: 10.1016/j.jmb.2007.04.063
  58. 58
    Finger, D. L.; Blanchard, M. S.; Theimer, C. A.; Sengerová, B.; Singh, P. S.; Chavez, V.; Liu, F.; Grasby, J. A.; Shen, B. The 3′-Flap Pocket of Human Flap Endonuclease 1 Is Critical for Substrate Binding and Catalysis. J. Biol. Chem. 2009, 284 (33), 2218422194,  DOI: 10.1074/jbc.M109.015065
  59. 59
    Kuznetsova, A. A.; Fedorova, O. S.; Kuznetsov, N. A. Kinetic Features of 3′-5′ Exonuclease Activity of Human AP-Endonuclease APE1. Molecules 2018, 23 (9), 2101,  DOI: 10.3390/molecules23092101
  60. 60
    Warshel, A.; Levitt, M. Theoretical Studies of Enzymic Reactions: Dielectric, Electrostatic and Steric Stabilization of the Carbonium Ion in the Reaction of Lysozyme. J. Mol. Biol. 1976, 103 (2), 227249,  DOI: 10.1016/0022-2836(76)90311-9
  61. 61
    Brunk, E.; Rothlisberger, U. Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. Chem. Rev. 2015, 115 (12), 62176263,  DOI: 10.1021/cr500628b
  62. 62
    Stevens, D. R.; Hammes-Schiffer, S. Exploring the Role of the Third Active Site Metal Ion in DNA Polymerase η with QM/MM Free Energy Simulations. J. Am. Chem. Soc. 2018, 140 (28), 89658969,  DOI: 10.1021/jacs.8b05177
  63. 63
    Field, M. J.; Bash, P. A.; Karplus, M. A Combined Quantum Mechanical and Molecular Mechanical Potential for Molecular Dynamics Simulations. J. Comput. Chem. 1990, 11 (6), 700733,  DOI: 10.1002/jcc.540110605
  64. 64
    Needleman, S. B.; Wunsch, C. D. A General Method Applicable to Search for Similarities in Amino Acid Sequence of Two Proteins. J. Mol. Biol. 1970, 48 (3), 443453,  DOI: 10.1016/0022-2836(70)90057-4
  65. 65
    Szymanski, M. R.; Yu, W.; Gmyrek, A. M.; White, M. A.; Molineux, I. J.; Lee, J. C.; Yin, Y. W. A Domain in Human EXOG Converts Apoptotic Endonuclease to DNA-Repair Exonuclease. Nat. Commun. 2017, 8, 14959,  DOI: 10.1038/ncomms14959
  66. 66
    Wu, C.-C.; Lin, J. L. J.; Yang-Yen, H.-F.; Yuan, H. S. A Unique Exonuclease ExoG Cleaves between RNA and DNA in Mitochondrial DNA Replication. Nucleic Acids Res. 2019, 47 (10), 54055419,  DOI: 10.1093/nar/gkz241
  67. 67
    Zhang, J.; Mccabe, K. A.; Bell, C. E. Crystal Structures of λ Exonuclease in Complex with DNA Suggest an Electrostatic Ratchet Mechanism for Processivity. Proc. Natl. Acad. Sci. U. S. A. 2011, 108 (29), 1187211877,  DOI: 10.1073/pnas.1103467108
  68. 68
    Zhang, J.; Pan, X.; Bell, C. E. Crystal Structure of λ Exonuclease in Complex with DNA and Ca2+. Biochemistry 2014, 53 (47), 74157425,  DOI: 10.1021/bi501155q
  69. 69
    Cheng, K.; Xu, H.; Chen, X.; Wang, L.; Tian, B.; Zhao, Y.; Hua, Y. Structural Basis for DNA 5′-End Resection by RecJ. eLife 2016, 5, e14294  DOI: 10.7554/eLife.14294
  70. 70
    Zhang, F.; Shi, J.; Chen, S. H.; Bian, C.; Yu, X. The PIN Domain of EXO1 Recognizes Poly(ADP-Ribose) in DNA Damage Response. Nucleic Acids Res. 2015, 43 (22), 1078210794,  DOI: 10.1093/nar/gkv939
  71. 71
    Genna, V.; Carloni, P.; De Vivo, M. A Strategically Located Arg/Lys Residue Promotes Correct Base Paring during Nucleic Acid Biosynthesis in Polymerases. J. Am. Chem. Soc. 2018, 140 (9), 33123321,  DOI: 10.1021/jacs.7b12446
  72. 72
    Palermo, G.; Miao, Y.; Walker, R. C.; Jinek, M.; McCammon, J. A. CRISPR-Cas9 Conformational Activation as Elucidated from Enhanced Molecular Simulations. Proc. Natl. Acad. Sci. U. S. A. 2017, 114 (28), 72607265,  DOI: 10.1073/pnas.1707645114
  73. 73
    Mulholland, A. J.; Roitberg, A. E.; Tuñón, I. Enzyme Dynamics and Catalysis in the Mechanism of DNA Polymerase. Theor. Chem. Acc. 2012, 131, 1286,  DOI: 10.1007/s00214-012-1286-8
  74. 74
    Riccardi, L.; Genna, V.; De Vivo, M. Metal–Ligand Interactions in Drug Design. Nat. Rev. Chem. 2018, 2 (7), 100112,  DOI: 10.1038/s41570-018-0018-6
  75. 75
    De Vivo, M. Bridging Quantum Mechanics and Structure-Based Drug Design. Front. Biosci., Landmark Ed. 2011, 16, 16191633,  DOI: 10.2741/3809
  76. 76
    Genna, V.; Marcia, M.; De Vivo, M. A Transient and Flexible Cation−π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases. J. Am. Chem. Soc. 2019, 141 (27), 1077010776,  DOI: 10.1021/jacs.9b03663
  77. 77
    Yan, C.; Dodd, T.; He, Y.; Tainer, J. A.; Tsutakawa, S. E.; Ivanov, I. Transcription Preinitiation Complex Structure and Dynamics Provide Insight into Genetic Diseases. Nat. Struct. Mol. Biol. 2019, 26 (6), 397406,  DOI: 10.1038/s41594-019-0220-3
  78. 78
    Casalino, L.; Palermo, G.; Spinello, A.; Rothlisberger, U.; Magistrato, A. All-Atom Simulations Disentangle the Functional Dynamics Underlying Gene Maturation in the Intron Lariat Spliceosome. Proc. Natl. Acad. Sci. U. S. A. 2018, 115 (26), 65846589,  DOI: 10.1073/pnas.1802963115
  79. 79
    Orellana, L.; Thorne, A. H.; Lema, R.; Gustavsson, J.; Parisian, A. D.; Hospital, A.; Cordeiro, T. N.; Bernadó, P.; Scott, A. M.; Brun-Heath, I.; Lindahl, E.; Cavenee, W. K.; Furnari, F. B.; Orozco, M. Oncogenic Mutations at the EGFR Ectodomain Structurally Converge to Remove a Steric Hindrance on a Kinase-Coupled Cryptic Epitope. Proc. Natl. Acad. Sci. U. S. A. 2019, 116 (20), 1000910018,  DOI: 10.1073/pnas.1821442116
  80. 80
    Palermo, G.; Casalino, L.; Magistrato, A.; McCammon, A. J. Understanding the Mechanistic Basis of Non-Coding RNA through Molecular Dynamics Simulations. J. Struct. Biol. 2019, 206 (3), 267279,  DOI: 10.1016/j.jsb.2019.03.004
  81. 81
    Maier, J. A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K. E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11 (8), 36963713,  DOI: 10.1021/acs.jctc.5b00255
  82. 82
    Galindo-Murillo, R.; Robertson, J. C.; Zgarbová, M.; Šponer, J.; Otyepka, M.; Jurečka, P.; Cheatham, T. E. Assessing the Current State of Amber Force Field Modifications for DNA. J. Chem. Theory Comput. 2016, 12 (8), 41144127,  DOI: 10.1021/acs.jctc.6b00186
  83. 83
    Dans, P. D.; Ivani, I.; Hospital, A.; Portella, G.; González, C.; Orozco, M. How Accurate Are Accurate Force-Fields for B-DNA?. Nucleic Acids Res. 2017, 45 (7), 42174230,  DOI: 10.1093/nar/gkw1355
  84. 84
    Wang, J.; Wolf, R. M.; Caldwell, J. W.; Kollman, P. A.; Case, D. A. Development and Testing of a General Amber Force Field. J. Comput. Chem. 2004, 25 (9), 11571174,  DOI: 10.1002/jcc.20035
  85. 85
    Singh, U. C.; Kollman, P. A. An Approach to Computing Electrostatic Charges for Molecules. J. Comput. Chem. 1984, 5 (2), 129145,  DOI: 10.1002/jcc.540050204
  86. 86
    Bayly, C. I.; Cieplak, P.; Cornell, W.; Kollman, P. A. A Well-Behaved Electrostatic Potential Based Method Using Charge Restraints for Deriving Atomic Charges: The RESP Model. J. Phys. Chem. 1993, 97 (40), 1026910280,  DOI: 10.1021/j100142a004
  87. 87
    Hess, B. P-LINCS: A Parallel Linear Constraint Solver for Molecular Simulation. J. Chem. Theory Comput. 2008, 4 (1), 116122,  DOI: 10.1021/ct700200b
  88. 88
    Abraham, M. J.; Murtola, T.; Schulz, R.; Pall, S.; Smith, J. C.; Hess, B.; Lindahl, E. Gromacs: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1–2, 1925,  DOI: 10.1016/j.softx.2015.06.001
  89. 89
    Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N·log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98 (12), 1008910092,  DOI: 10.1063/1.464397
  90. 90
    Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103 (19), 85778593,  DOI: 10.1063/1.470117
  91. 91
    Matta, C. F.; Bader, R. F. W. Atoms-in-molecules study of the genetically encoded amino acids. III. Bond and atomic properties and their correlations with experiment including mutation-induced changes in protein stability and genetic coding. Proteins: Struct., Funct., Genet. 2003, 52 (3), 360399,  DOI: 10.1002/prot.10414
  92. 92
    Peraro, M. D.; Spiegel, K.; Lamoureux, G.; De Vivo, M.; DeGrado, W. F.; Klein, M. L. Modeling the Charge Distribution at Metal Sites in Proteins for Molecular Dynamics Simulations. J. Struct. Biol. 2007, 157 (3), 444453,  DOI: 10.1016/j.jsb.2006.10.019
  93. 93
    Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79 (2), 926935,  DOI: 10.1063/1.445869
  94. 94
    Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126 (1), 014101,  DOI: 10.1063/1.2408420
  95. 95
    Parrinello, M.; Rahman, A. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52 (12), 71827190,  DOI: 10.1063/1.328693
  96. 96
    Barducci, A.; Bussi, G.; Parrinello, M. Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method. Phys. Rev. Lett. 2008, 100 (2), 020603  DOI: 10.1103/PhysRevLett.100.020603

Cited By

Click to copy section linkSection link copied!
Citation Statements
Explore this article's citation statements on scite.ai

This article is cited by 14 publications.

  1. Dylan J. Nikkel, Rajwinder Kaur, Stacey D. Wetmore. How Can One Metal Power Nucleic Acid Phosphodiester Bond Cleavage by a Nuclease? Multiscale Computational Studies Highlight a Diverse Mechanistic Landscape. The Journal of Physical Chemistry B 2025, 129 (1) , 3-18. https://doi.org/10.1021/acs.jpcb.4c05875
  2. Rajwinder Kaur, Dylan J. Nikkel, Stacey D. Wetmore. Mechanism of Nucleic Acid Phosphodiester Bond Cleavage by Human Endonuclease V: MD and QM/MM Calculations Reveal a Versatile Metal Dependence. The Journal of Physical Chemistry B 2024, 128 (39) , 9455-9469. https://doi.org/10.1021/acs.jpcb.4c05846
  3. Dylan J. Nikkel, Stacey D. Wetmore. Distinctive Formation of a DNA–Protein Cross-Link during the Repair of DNA Oxidative Damage: Insights into Human Disease from MD Simulations and QM/MM Calculations. Journal of the American Chemical Society 2023, 145 (24) , 13114-13125. https://doi.org/10.1021/jacs.3c01773
  4. Elisa Donati, Pietro Vidossich, Marco De Vivo. Molecular Mechanism of Phosphate Steering for DNA Binding, Cleavage Localization, and Substrate Release in Nucleases. ACS Catalysis 2021, 11 (21) , 13244-13254. https://doi.org/10.1021/acscatal.1c03346
  5. Jacopo Manigrasso, Marco De Vivo, Giulia Palermo. Controlled Trafficking of Multiple and Diverse Cations Prompts Nucleic Acid Hydrolysis. ACS Catalysis 2021, 11 (14) , 8786-8797. https://doi.org/10.1021/acscatal.1c01825
  6. Marco Marcia, Jacopo Manigrasso, Marco De Vivo. Finding the Ion in the RNA-Stack: Can Computational Models Accurately Predict Key Functional Elements in Large Macromolecular Complexes?. Journal of Chemical Information and Modeling 2021, 61 (6) , 2511-2515. https://doi.org/10.1021/acs.jcim.1c00572
  7. Pietro Vidossich, Luis Eduardo Castañeda Moreno, Cristiano Mota, Daniele de Sanctis, Gian Pietro Miscione, Marco De Vivo. Functional Implications of Second-Shell Basic Residues for dUTPase DR2231 Enzymatic Specificity. ACS Catalysis 2020, 10 (23) , 13825-13833. https://doi.org/10.1021/acscatal.0c04148
  8. Yanbin Xu, Anlian Huang, Wei Yi, Guosheng Chen, Siming Huang, Gangfeng Ouyang. Nanozyme engineering in structurally explicit framework: Design mechanisms and biosensing applications. Coordination Chemistry Reviews 2024, 500 , 215517. https://doi.org/10.1016/j.ccr.2023.215517
  9. Vito Genna, Javier Iglesias-Fernández, Laura Reyes-Fraile, Nuria Villegas, Kevin Guckian, Punit Seth, Brad Wan, Cristina Cabrero, Montserrat Terrazas, Isabelle Brun-Heath, Carlos González, Simone Sciabola, Anabella Villalobos, Modesto Orozco. Controlled sulfur-based engineering confers mouldability to phosphorothioate antisense oligonucleotides. Nucleic Acids Research 2023, 51 (10) , 4713-4725. https://doi.org/10.1093/nar/gkad309
  10. T. Kubař, M. Elstner, Q. Cui. Hybrid Quantum Mechanical/Molecular Mechanical Methods For Studying Energy Transduction in Biomolecular Machines. Annual Review of Biophysics 2023, 52 (1) , 525-551. https://doi.org/10.1146/annurev-biophys-111622-091140
  11. Caleb Chang, Christie Lee Luo, Yang Gao. In crystallo observation of three metal ion promoted DNA polymerase misincorporation. Nature Communications 2022, 13 (1) https://doi.org/10.1038/s41467-022-30005-3
  12. Inacrist Geronimo, Pietro Vidossich, Elisa Donati, Marco De Vivo. Computational investigations of polymerase enzymes: Structure, function, inhibition, and biotechnology. WIREs Computational Molecular Science 2021, 11 (6) https://doi.org/10.1002/wcms.1534
  13. Łukasz Nierzwicki, Pablo Ricardo Arantes, Aakash Saha, Giulia Palermo. Establishing the allosteric mechanism in CRISPR‐Cas9. WIREs Computational Molecular Science 2021, 11 (3) https://doi.org/10.1002/wcms.1503
  14. Concetta Cozza, Françisco M. Raymo, Adriana Pietropaolo. Shape factors in the binding of soft fluorescent nanoshuttles with target receptors. Molecular Systems Design & Engineering 2021, 6 (4) , 281-285. https://doi.org/10.1039/D0ME00164C

Journal of the American Chemical Society

Cite this: J. Am. Chem. Soc. 2020, 142, 6, 2823–2834
Click to copy citationCitation copied!
https://doi.org/10.1021/jacs.9b10656
Published January 15, 2020

Copyright © 2020 American Chemical Society. This publication is licensed under

CC-BY 4.0 .

Article Views

2381

Altmetric

-

Citations

Learn about these metrics

Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

  • Abstract

    Figure 1

    Figure 1. Catalytic domain of hExo1 in complex with DNA substrate and the two catalytic metal ions (PDB ID 5V06). (Left) hExo1 in cartoon and with colors for different structural motifs. (Right) Closer view of the active site. Three metal ions (MgA, MgB, MgC) are in orange, nucleophilic water molecule (Nu) is in red, two guide residues (Tyr32, His36) in yellow, and residues of the catalytic pocket (Gly2, Asp30, Asp78, Asp152, Asp171, Asp173) are in cyan. Scissile phosphate is correctly positioned for the nucleophilic attack, and MgC is coordinated by the 5′ terminal phosphate.

    Scheme 1

    Scheme 1. Schematic Representation of the Motion of the Transient Third Metal Ion, Which We Found to Be Intermittently Recruited/Released by Glu89 (switching its inner/outer conformations) during Exonuclease Catalysisa

    aThis structural evidence suggests that the transient third ion may be crucial for substrate hydrolysis and/or leaving group departure.

    Figure 2

    Figure 2. Radial distribution function, g(r), calculated for ions around 10 Å from the center of mass of the 5′ phosphate group. Plot shows the presence of ions ∼3 Å from the 5′ phosphate group for the RS2M system. In this system, a K+ ion approached the negatively charged group. For the RSGlu89Ala and PSGlu89Ala systems, there are no ions within ∼5.5 Å of the 5′ phosphate, as indicated by the g(r) values of ∼0. In the upper right corner, the 5′ phosphate group and Glu89Ala residues are shown in licorice (taken from the PSGlu89Ala simulations).

    Figure 3

    Figure 3. (A) Distance (dMG in yellow) between the third Mg2+ ion, from the bulk (Mgbulk), and the phosphorus of the 5′ phosphate group of AMP. Inset, representation (snapshot from the PS2M simulations) of Mgbulk approaching the terminal 5′ phosphate as well as the pseudo-dihedral angle ϕ of Glu89 side chain (defined by the N–Cα–Cδ−-Cγ atoms). (B) Probability density of the pseudo-dihedral angle ϕ in Glu89 during the simulation: (blue) probability density for dMG values > 4 Å shows the outer conformation as the most populated; (red) probability density for dMG values < 4 Å shows the inner conformation is the most populated.

    Figure 4

    Figure 4. (A) Graphic representation of d1 and d2 distances (PDB ID 5V0A). (B) Probability density of the distances d1 and d2 calculated during simulations of the systems PS2M (in red), PS3M (in orange), and PSGlu89Ala (in green).

    Figure 5

    Figure 5. (Bottom) Free energy surface obtained through well-tempered metadynamics simulations for RS3M (red), RS2M (blue), PS2M (green), and PS3M (light purple) systems. Results show three conformations (outer, intermediate, inner). (Top) Graphic representations, taken from PS3M simulaitons, of the three conformations are shown in licorice.

    Figure 6

    Figure 6. Free energy surface obtained through confined well-tempered metadynamics simulations for PS2M (blue) and PS3M (yellow) systems. (Top) Schematic representation of the CV1, exemplified using a snapshot from PS3M simulations. It represents the distance between the center of mass (COM) of the heavy atoms of the nucleotide leaving group and the COM of the Cα of the aspartates (Asp152, Asp171, Asp173) in the first coordination shell of MgA, MgB. Two different minima, at 8.8 and 12.4 Å, agree with the MD results, in which a partial exit of the leaving group (CV1 ≈ 12.4 Å) was seen only in the presence of MgC (Figure S8).

    Figure 7

    Figure 7. Close views of the active site of 5′ metallonuclease members that possess an analogous acid residue (light green) close to the two-metal-ion center (MA, MB, in orange), the active site (in cyan), and the leaving group (indicated by a dashed line). (A) Human ExoG, in which Glu317 is pointing in the inner (PDB ID 5T5C) and (A′) outer conformations (PDB ID 5T40, merged with the DNA substrate from PDB ID 5T5C). (B) Escherichia phage T5Fen (PDB ID 5HNK). (C) Human λ-Exonuclease X (PDB ID 3SM4). (D) D. radiodurans RecJ (PDB ID 5F55). (E) Sequence alignment of M. smegmatis FenA and E. phage T5Fen. Conserved acid residue (Glu/Asp) is indicated in orange.

  • References


    This article references 96 other publications.

    1. 1
      Yang, W.; Weng, P. J.; Gao, Y. A New Paradigm of DNA Synthesis: Three-Metal-Ion Catalysis. Cell Biosci. 2016, 6 (1), 51,  DOI: 10.1186/s13578-016-0118-2
    2. 2
      Gao, Y.; Yang, W. Capture of a Third Mg2+ Is Essential for Catalyzing DNA Synthesis. Science 2016, 352 (6291), 13341337,  DOI: 10.1126/science.aad9633
    3. 3
      Genna, V.; Colombo, M.; De Vivo, M.; Marcia, M. Second-Shell Basic Residues Expand the Two-Metal-Ion Architecture of DNA and RNA Processing Enzymes. Structure 2018, 26 (1), 4050,  DOI: 10.1016/j.str.2017.11.008
    4. 4
      Genna, V.; Donati, E.; De Vivo, M. The Catalytic Mechanism of DNA and RNA Polymerases. ACS Catal. 2018, 8 (12), 1110311118,  DOI: 10.1021/acscatal.8b03363
    5. 5
      Pavlov, Y. I.; Shcherbakova, P. V.; Rogozin, I. B. Roles of DNA Polymerases in Replication, Repair, and Recombination in Eukaryotes. Int. Rev. Cytol. 2006, 255, 41132,  DOI: 10.1016/S0074-7696(06)55002-8
    6. 6
      Nishino, T.; Morikawa, K. Structure and Function of Nucleases in DNA Repair: Shape, Grip and Blade of the DNA Scissors. Oncogene 2002, 21 (58), 90229032,  DOI: 10.1038/sj.onc.1206135
    7. 7
      Marti, T. M.; Fleck, O. DNA Repair Nucleases. Cell. Mol. Life Sci. 2004, 61 (3), 336354,  DOI: 10.1007/s00018-003-3223-4
    8. 8
      Potapov, V.; Fu, X.; Dai, N.; Corrêa, I. R.; Tanner, N. A.; Ong, J. L. Base Modifications Affecting RNA Polymerase and Reverse Transcriptase Fidelity. Nucleic Acids Res. 2018, 46 (11), 57535763,  DOI: 10.1093/nar/gky341
    9. 9
      Nimonkar, A. V.; Ozsoy, A. Z.; Genschel, J.; Modrich, P.; Kowalczykowski, S. C. Human Exonuclease 1 and BLM Helicase Interact to Resect DNA and Initiate DNA Repair. Proc. Natl. Acad. Sci. U. S. A. 2008, 105 (44), 1690616911,  DOI: 10.1073/pnas.0809380105
    10. 10
      Loeb, L. A.; Monnat, R. J. DNA Polymerases and Human Disease. Nat. Rev. Genet. 2008, 9 (8), 594604,  DOI: 10.1038/nrg2345
    11. 11
      Zheng, L.; Jia, J.; Finger, L. D.; Guo, Z.; Zer, C.; Shen, B. Functional Regulation of FEN1 Nuclease and Its Link to Cancer. Nucleic Acids Res. 2011, 39 (3), 781794,  DOI: 10.1093/nar/gkq884
    12. 12
      Peltomäki, P. Role of DNA Mismatch Repair Defects in the Pathogenesis of Human Cancer. J. Clin. Oncol. 2003, 21 (6), 11741179,  DOI: 10.1200/JCO.2003.04.060
    13. 13
      Dai, Y.; Tang, Z.; Yang, Z.; Zhang, L.; Deng, Q.; Zhang, X.; Yu, Y.; Liu, X.; Zhu, J. EXO1 Overexpression Is Associated with Poor Prognosis of Hepatocellular Carcinoma Patients. Cell Cycle 2018, 17 (19–20), 23862397,  DOI: 10.1080/15384101.2018.1534511
    14. 14
      Ivanov, I.; Tainer, J. A.; McCammon, J. A. Unraveling the Three-Metal-Ion Catalytic Mechanism of the DNA Repair Enzyme Endonuclease IV. Proc. Natl. Acad. Sci. U. S. A. 2007, 104 (5), 14651470,  DOI: 10.1073/pnas.0603468104
    15. 15
      Prieto, J.; Redondo, P.; Merino, N.; Villate, M.; Montoya, G.; Blanco, F. J.; Molina, R. Structure of the I-SceI Nuclease Complexed with Its DsDNA Target and Three Catalytic Metal Ions. Acta Crystallogr., Sect. F: Struct. Biol. Commun. 2016, 72 (6), 473479,  DOI: 10.1107/S2053230X16007512
    16. 16
      AlMalki, F. A.; Flemming, C. S.; Zhang, J.; Feng, M.; Sedelnikova, S. E.; Ceska, T.; Rafferty, J. B.; Sayers, J. R.; Artymiuk, P. J. Direct Observation of DNA Threading in Flap Endonuclease Complexes. Nat. Struct. Mol. Biol. 2016, 23 (7), 640646,  DOI: 10.1038/nsmb.3241
    17. 17
      Uson, M. L.; Carl, A.; Goldgur, Y.; Shuman, S. Crystal Structure and Mutational Analysis of Mycobacterium Smegmatis FenA Highlight Active Site Amino Acids and Three Metal Ions Essential for Flap Endonuclease and 5 Exonuclease Activities. Nucleic Acids Res. 2018, 46 (8), 41644175,  DOI: 10.1093/nar/gky238
    18. 18
      Yang, W. Nucleases: Diversity of Structure, Function and Mechanism. Q. Rev. Biophys. 2011, 44 (1), 193,  DOI: 10.1017/S0033583510000181
    19. 19
      Palermo, G.; Cavalli, A.; Klein, M. L.; Alfonso-Prieto, M.; Dal Peraro, M.; De Vivo, M. Catalytic Metal Ions and Enzymatic Processing of DNA and RNA. Acc. Chem. Res. 2015, 48 (2), 220228,  DOI: 10.1021/ar500314j
    20. 20
      Palermo, G.; Stenta, M.; Cavalli, A.; Dal Peraro, M.; De Vivo, M. Molecular Simulations Highlight the Role of Metals in Catalysis and Inhibition of Type II Topoisomerase. J. Chem. Theory Comput. 2013, 9 (2), 857862,  DOI: 10.1021/ct300691u
    21. 21
      Schmidt, B. H.; Burgin, A. B.; Deweese, J. E.; Osheroff, N.; Berger, J. M. A Novel and Unified Two-Metal Mechanism for DNA Cleavage by Type II and IA Topoisomerases. Nature 2010, 465 (7298), 641644,  DOI: 10.1038/nature08974
    22. 22
      Perera, L.; Freudenthal, B. D.; Beard, W. A.; Shock, D. D.; Pedersen, L. G.; Wilson, S. H. Requirement for Transient Metal Ions Revealed through Computational Analysis for DNA Polymerase Going in Reverse. Proc. Natl. Acad. Sci. U. S. A. 2015, 112 (38), E5228E5236,  DOI: 10.1073/pnas.1511207112
    23. 23
      Shi, Y.; Hellinga, H. W.; Beese, L. S. Interplay of Catalysis, Fidelity, Threading, and Processivity in the Exo- and Endonucleolytic Reactions of Human Exonuclease I. Proc. Natl. Acad. Sci. U. S. A. 2017, 114 (23), 60106015,  DOI: 10.1073/pnas.1704845114
    24. 24
      Jimeno, S.; Herrera-Moyano, E.; Ortega, P.; Aguilera, A. Differential Effect of the Overexpression of Rad2/XPG Family Endonucleases on Genome Integrity in Yeast and Human Cells. DNA Repair 2017, 57, 6675,  DOI: 10.1016/j.dnarep.2017.06.030
    25. 25
      Emmert, S.; Schneider, T. D.; Khan, S. G.; Kraemer, K. H. The Human XPG Gene: Gene Architecture, Alternative Splicing and Single Nucleotide Polymorphisms. Nucleic Acids Res. 2001, 29 (7), 14431452,  DOI: 10.1093/nar/29.7.1443
    26. 26
      Qiu, J.; Qian, Y.; Chen, V.; Guan, M.-X.; Shen, B. Human Exonuclease 1 Functionally Complements Its Yeast Homologues in DNA Recombination, RNA Primer Removal, and Mutation Avoidance. J. Biol. Chem. 1999, 274 (25), 1789317900,  DOI: 10.1074/jbc.274.25.17893
    27. 27
      Genschel, J.; Bazemore, L. R.; Modrich, P. Human Exonuclease I Is Required for 5′ and 3′ Mismatch Repair. J. Biol. Chem. 2002, 277 (15), 1330213311,  DOI: 10.1074/jbc.M111854200
    28. 28
      Wei, K.; Clark, A. B.; Wong, E.; Kane, M. F.; Mazur, D. J.; Parris, T.; Kolas, N. K.; Russell, R.; Hou, H.; Kneitz, B.; Yang, G.; Kunkel, T. A.; Kolodner, R. D.; Cohen, P. E.; Edelmann, W. Inactivation of Exonuclease I in Mice Results in DNA Mismatch Repair Defects, Increased Cancer Susceptibility, and Male and Female Sterility. Genes Dev. 2003, 17 (5), 603614,  DOI: 10.1101/gad.1060603
    29. 29
      Zhu, Z.; Chung, W. H.; Shim, E. Y.; Lee, S. E.; Ira, G. Sgs1 Helicase and Two Nucleases Dna2 and Exo1 Resect DNA Double-Strand Break Ends. Cell 2008, 134 (6), 981994,  DOI: 10.1016/j.cell.2008.08.037
    30. 30
      Vallur, A. C.; Maizels, N. Complementary Roles for Exonuclease 1 and Flap Endonuclease 1 in Maintenance of Triplet Repeats. J. Biol. Chem. 2010, 285 (37), 2851428519,  DOI: 10.1074/jbc.M110.132738
    31. 31
      Lee, B. I.; Wilson, D. M. The RAD2 Domain of Human Exonuclease 1 Exhibits 5′ to 3′ Exonuclease and Flap Structure-Specific Endonuclease Activities. J. Biol. Chem. 1999, 274 (53), 3776337769,  DOI: 10.1074/jbc.274.53.37763
    32. 32
      Keijzers, G.; Bohr, V. A.; Rasmussen, L. J. Human Exonuclease 1 (EXO1) Activity Characterization and Its Function on FLAP Structures. Biosci. Rep. 2015, 35 (3), e00206  DOI: 10.1042/BSR20150058
    33. 33
      Orans, J.; McSweeney, E. A.; Iyer, R. R.; Hast, M. A.; Hellinga, H. W.; Modrich, P.; Beese, L. S. Structures of Human Exonuclease 1 DNA Complexes Suggest a Unified Mechanism for Nuclease Family. Cell 2011, 145 (2), 212223,  DOI: 10.1016/j.cell.2011.03.005
    34. 34
      Tsutakawa, S. E.; Thompson, M. J.; Arvai, A. S.; Neil, A. J.; Shaw, S. J.; Algasaier, S. I.; Kim, J. C.; Finger, L. D.; Jardine, E.; Gotham, V. J. B.; Sarker, A. H.; Her, M. Z.; Rashid, F.; Hamdan, S. M.; Mirkin, S. M.; Grasby, J. A.; Tainer, J. A. Phosphate Steering by Flap Endonuclease 1 Promotes 5′-Flap Specificity and Incision to Prevent Genome Instability. Nat. Commun. 2017, 8, 15855,  DOI: 10.1038/ncomms15855
    35. 35
      Steitz, T. A.; Steitz, J. A. A General Two-Metal-Ion Mechanism for Catalytic RNA. Proc. Natl. Acad. Sci. U. S. A. 1993, 90 (14), 64986502,  DOI: 10.1073/pnas.90.14.6498
    36. 36
      Yang, W.; Lee, J. Y.; Nowotny, M. Making and Breaking Nucleic Acids: Two-Mg2+-Ion Catalysis and Substrate Specificity. Mol. Cell 2006, 22 (1), 513,  DOI: 10.1016/j.molcel.2006.03.013
    37. 37
      Dupureur, C. M. One Is Enough: Insights into the Two-Metal Ion Nuclease Mechanism from Global Analysis and Computational Studies. Metallomics. 2010, 2 (9), 609620,  DOI: 10.1039/c0mt00013b
    38. 38
      Beese, L. S.; Steitz, T. A. Structural basis for the 3'-5' exonuclease activity of Escherichia coli DNA polymerase 1: a two metal ion mechanism. EMBO J. 1991, 10 (1), 2533,  DOI: 10.1002/j.1460-2075.1991.tb07917.x
    39. 39
      De Vivo, M.; Dal Peraro, M.; Klein, M. L. Phosphodiester Cleavage in Ribonuclease H Occurs via an Associative Two-Metal-Aided Catalytic Mechanism. J. Am. Chem. Soc. 2008, 130 (33), 1095510962,  DOI: 10.1021/ja8005786
    40. 40
      Genna, V.; Vidossich, P.; Ippoliti, E.; Carloni, P.; De Vivo, M. A Self-Activated Mechanism for Nucleic Acid Polymerization Catalyzed by DNA/RNA Polymerases. J. Am. Chem. Soc. 2016, 138 (44), 1459214598,  DOI: 10.1021/jacs.6b05475
    41. 41
      Shaw, S. J.; Finger, L. D.; Grasby, J. A. Human Exonuclease 1 Threads 5′-Flap Substrates through Its Helical Arch. Biochemistry 2017, 56 (29), 37043707,  DOI: 10.1021/acs.biochem.7b00507
    42. 42
      Tomlinson, C. G.; Atack, J. M.; Chapados, B.; Tainer, J. A.; Grasby, J. A. Substrate Recognition and Catalysis by Flap Endonucleases and Related Enzymes. Biochem. Soc. Trans. 2010, 38 (2), 433437,  DOI: 10.1042/BST0380433
    43. 43
      Bennet, I. A; Finger, L D.; Baxter, N. J; Ambrose, B.; Hounslow, A. M; Thompson, M. J; Exell, J. C; Shahari, N. N. B M.; Craggs, T. D; Waltho, J. P; Grasby, J. A Regional Conformational Flexibility Couples Substrate Specificity and Scissile Phosphate Diester Selectivity in Human Flap Endonuclease 1. Nucleic Acids Res. 2018, 46 (11), 56185633,  DOI: 10.1093/nar/gky293
    44. 44
      Raper, A. T.; Reed, A. J.; Suo, Z. Kinetic Mechanism of DNA Polymerases: Contributions of Conformational Dynamics and a Third Divalent Metal Ion. Chem. Rev. 2018, 118 (2), 60006025,  DOI: 10.1021/acs.chemrev.7b00685
    45. 45
      Black, C. B.; Huang, H. W.; Cowan, J. A. Biological Coordination Chemistry of Magnesium, Sodium, and Potassium Ions. Protein and Nucleotide Binding Sites. Coord. Chem. Rev. 1994, 135–136 (C), 165202,  DOI: 10.1016/0010-8545(94)80068-5
    46. 46
      Ho, M.-H.; De Vivo, M.; Dal Peraro, M.; Klein, M. L. Understanding the Effect of Magnesium Ion Concentration on the Catalytic Activity of Ribonuclease H through Computation: Does a Third Metal Binding Site Modulate Endonuclease Activity?. J. Am. Chem. Soc. 2010, 132 (39), 1370213712,  DOI: 10.1021/ja102933y
    47. 47
      Genna, V.; Gaspari, R.; Dal Peraro, M.; De Vivo, M. Cooperative Motion of a Key Positively Charged Residue and Metal Ions for DNA Replication Catalyzed by Human DNA Polymerase-η. Nucleic Acids Res. 2016, 44 (6), 28272836,  DOI: 10.1093/nar/gkw128
    48. 48
      La Sala, G.; Riccardi, L.; Gaspari, R.; Cavalli, A.; Hantschel, O.; De Vivo, M. HRD Motif as the Central Hub of the Signaling Network for Activation Loop Autophosphorylation in Abl Kinase. J. Chem. Theory Comput. 2016, 12 (11), 55635574,  DOI: 10.1021/acs.jctc.6b00600
    49. 49
      Hwang, W.; Yoo, J.; Lee, Y.; Park, S.; Hoang, P. L.; Cho, H.; Yu, J.; Hoa Vo, T. M.; Shin, M.; Jin, M. S.; Park, D.; Hyeon, C.; Lee, G. Dynamic Coordination of Two-Metal-Ions Orchestrates λ-Exonuclease Catalysis. Nat. Commun. 2018, 9 (1), 4404,  DOI: 10.1038/s41467-018-06750-9
    50. 50
      Samara, N. L.; Yang, W. Cation Trafficking Propels RNA Hydrolysis. Nat. Struct. Mol. Biol. 2018, 25 (8), 715721,  DOI: 10.1038/s41594-018-0099-4
    51. 51
      Molina, R.; Stella, S.; Redondo, P.; Gomez, H.; Marcaida, M. J.; Orozco, M.; Prieto, J.; Montoya, G. Visualizing Phosphodiester-Bond Hydrolysis by an Endonuclease. Nat. Struct. Mol. Biol. 2015, 22 (1), 6572,  DOI: 10.1038/nsmb.2932
    52. 52
      Sasnauskas, G.; Jeltsch, A.; Pingoud, A.; Siksnys, V. Plasmid DNA Cleavage by MunI Restriction Enzyme: Single-Turnover and Steady-State Kinetic Analysis. Biochemistry 1999, 38 (13), 40284036,  DOI: 10.1021/bi982456n
    53. 53
      Xie, F.; Dupureur, C. M. Kinetic Analysis of Product Release and Metal Ions in a Metallonuclease. Arch. Biochem. Biophys. 2009, 483 (1), 19,  DOI: 10.1016/j.abb.2009.01.001
    54. 54
      Sengerová, B.; Tomlinson, C.; Atack, J. M.; Williams, R.; Sayers, J. R.; Williams, N. H.; Grasby, J. A. Brønsted Analysis and Rate-Limiting Steps for the T5 Flap Endonuclease Catalyzed Hydrolysis of Exonucleolytic Substrates. Biochemistry 2010, 49 (37), 80858093,  DOI: 10.1021/bi100895j
    55. 55
      Taylor, J. D.; Halford, S. E. Discrimination between DNA Sequences by the EcoRV Restriction Endonuclease. Biochemistry 1989, 28 (15), 61986207,  DOI: 10.1021/bi00441a011
    56. 56
      Yang, C. C.; Baxter, B. K.; Topal, M. D. DNA Cleavage by NaeI: Protein Purification, Rate-Limiting Step, and Accuracy. Biochemistry 1994, 33 (49), 1491814925,  DOI: 10.1021/bi00253a031
    57. 57
      Williams, R.; Sengerova, B.; Osborne, S.; Syson, K.; Ault, S.; Kilgour, A.; Chapados, B. R.; Tainer, J. A.; Sayers, J. R.; Grasby, J. A. Comparison of the Catalytic Parameters and Reaction Specificities of a Phage and an Archaeal Flap Endonuclease. J. Mol. Biol. 2007, 371 (1), 3448,  DOI: 10.1016/j.jmb.2007.04.063
    58. 58
      Finger, D. L.; Blanchard, M. S.; Theimer, C. A.; Sengerová, B.; Singh, P. S.; Chavez, V.; Liu, F.; Grasby, J. A.; Shen, B. The 3′-Flap Pocket of Human Flap Endonuclease 1 Is Critical for Substrate Binding and Catalysis. J. Biol. Chem. 2009, 284 (33), 2218422194,  DOI: 10.1074/jbc.M109.015065
    59. 59
      Kuznetsova, A. A.; Fedorova, O. S.; Kuznetsov, N. A. Kinetic Features of 3′-5′ Exonuclease Activity of Human AP-Endonuclease APE1. Molecules 2018, 23 (9), 2101,  DOI: 10.3390/molecules23092101
    60. 60
      Warshel, A.; Levitt, M. Theoretical Studies of Enzymic Reactions: Dielectric, Electrostatic and Steric Stabilization of the Carbonium Ion in the Reaction of Lysozyme. J. Mol. Biol. 1976, 103 (2), 227249,  DOI: 10.1016/0022-2836(76)90311-9
    61. 61
      Brunk, E.; Rothlisberger, U. Mixed Quantum Mechanical/Molecular Mechanical Molecular Dynamics Simulations of Biological Systems in Ground and Electronically Excited States. Chem. Rev. 2015, 115 (12), 62176263,  DOI: 10.1021/cr500628b
    62. 62
      Stevens, D. R.; Hammes-Schiffer, S. Exploring the Role of the Third Active Site Metal Ion in DNA Polymerase η with QM/MM Free Energy Simulations. J. Am. Chem. Soc. 2018, 140 (28), 89658969,  DOI: 10.1021/jacs.8b05177
    63. 63
      Field, M. J.; Bash, P. A.; Karplus, M. A Combined Quantum Mechanical and Molecular Mechanical Potential for Molecular Dynamics Simulations. J. Comput. Chem. 1990, 11 (6), 700733,  DOI: 10.1002/jcc.540110605
    64. 64
      Needleman, S. B.; Wunsch, C. D. A General Method Applicable to Search for Similarities in Amino Acid Sequence of Two Proteins. J. Mol. Biol. 1970, 48 (3), 443453,  DOI: 10.1016/0022-2836(70)90057-4
    65. 65
      Szymanski, M. R.; Yu, W.; Gmyrek, A. M.; White, M. A.; Molineux, I. J.; Lee, J. C.; Yin, Y. W. A Domain in Human EXOG Converts Apoptotic Endonuclease to DNA-Repair Exonuclease. Nat. Commun. 2017, 8, 14959,  DOI: 10.1038/ncomms14959
    66. 66
      Wu, C.-C.; Lin, J. L. J.; Yang-Yen, H.-F.; Yuan, H. S. A Unique Exonuclease ExoG Cleaves between RNA and DNA in Mitochondrial DNA Replication. Nucleic Acids Res. 2019, 47 (10), 54055419,  DOI: 10.1093/nar/gkz241
    67. 67
      Zhang, J.; Mccabe, K. A.; Bell, C. E. Crystal Structures of λ Exonuclease in Complex with DNA Suggest an Electrostatic Ratchet Mechanism for Processivity. Proc. Natl. Acad. Sci. U. S. A. 2011, 108 (29), 1187211877,  DOI: 10.1073/pnas.1103467108
    68. 68
      Zhang, J.; Pan, X.; Bell, C. E. Crystal Structure of λ Exonuclease in Complex with DNA and Ca2+. Biochemistry 2014, 53 (47), 74157425,  DOI: 10.1021/bi501155q
    69. 69
      Cheng, K.; Xu, H.; Chen, X.; Wang, L.; Tian, B.; Zhao, Y.; Hua, Y. Structural Basis for DNA 5′-End Resection by RecJ. eLife 2016, 5, e14294  DOI: 10.7554/eLife.14294
    70. 70
      Zhang, F.; Shi, J.; Chen, S. H.; Bian, C.; Yu, X. The PIN Domain of EXO1 Recognizes Poly(ADP-Ribose) in DNA Damage Response. Nucleic Acids Res. 2015, 43 (22), 1078210794,  DOI: 10.1093/nar/gkv939
    71. 71
      Genna, V.; Carloni, P.; De Vivo, M. A Strategically Located Arg/Lys Residue Promotes Correct Base Paring during Nucleic Acid Biosynthesis in Polymerases. J. Am. Chem. Soc. 2018, 140 (9), 33123321,  DOI: 10.1021/jacs.7b12446
    72. 72
      Palermo, G.; Miao, Y.; Walker, R. C.; Jinek, M.; McCammon, J. A. CRISPR-Cas9 Conformational Activation as Elucidated from Enhanced Molecular Simulations. Proc. Natl. Acad. Sci. U. S. A. 2017, 114 (28), 72607265,  DOI: 10.1073/pnas.1707645114
    73. 73
      Mulholland, A. J.; Roitberg, A. E.; Tuñón, I. Enzyme Dynamics and Catalysis in the Mechanism of DNA Polymerase. Theor. Chem. Acc. 2012, 131, 1286,  DOI: 10.1007/s00214-012-1286-8
    74. 74
      Riccardi, L.; Genna, V.; De Vivo, M. Metal–Ligand Interactions in Drug Design. Nat. Rev. Chem. 2018, 2 (7), 100112,  DOI: 10.1038/s41570-018-0018-6
    75. 75
      De Vivo, M. Bridging Quantum Mechanics and Structure-Based Drug Design. Front. Biosci., Landmark Ed. 2011, 16, 16191633,  DOI: 10.2741/3809
    76. 76
      Genna, V.; Marcia, M.; De Vivo, M. A Transient and Flexible Cation−π Interaction Promotes Hydrolysis of Nucleic Acids in DNA and RNA Nucleases. J. Am. Chem. Soc. 2019, 141 (27), 1077010776,  DOI: 10.1021/jacs.9b03663
    77. 77
      Yan, C.; Dodd, T.; He, Y.; Tainer, J. A.; Tsutakawa, S. E.; Ivanov, I. Transcription Preinitiation Complex Structure and Dynamics Provide Insight into Genetic Diseases. Nat. Struct. Mol. Biol. 2019, 26 (6), 397406,  DOI: 10.1038/s41594-019-0220-3
    78. 78
      Casalino, L.; Palermo, G.; Spinello, A.; Rothlisberger, U.; Magistrato, A. All-Atom Simulations Disentangle the Functional Dynamics Underlying Gene Maturation in the Intron Lariat Spliceosome. Proc. Natl. Acad. Sci. U. S. A. 2018, 115 (26), 65846589,  DOI: 10.1073/pnas.1802963115
    79. 79
      Orellana, L.; Thorne, A. H.; Lema, R.; Gustavsson, J.; Parisian, A. D.; Hospital, A.; Cordeiro, T. N.; Bernadó, P.; Scott, A. M.; Brun-Heath, I.; Lindahl, E.; Cavenee, W. K.; Furnari, F. B.; Orozco, M. Oncogenic Mutations at the EGFR Ectodomain Structurally Converge to Remove a Steric Hindrance on a Kinase-Coupled Cryptic Epitope. Proc. Natl. Acad. Sci. U. S. A. 2019, 116 (20), 1000910018,  DOI: 10.1073/pnas.1821442116
    80. 80
      Palermo, G.; Casalino, L.; Magistrato, A.; McCammon, A. J. Understanding the Mechanistic Basis of Non-Coding RNA through Molecular Dynamics Simulations. J. Struct. Biol. 2019, 206 (3), 267279,  DOI: 10.1016/j.jsb.2019.03.004
    81. 81
      Maier, J. A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K. E.; Simmerling, C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J. Chem. Theory Comput. 2015, 11 (8), 36963713,  DOI: 10.1021/acs.jctc.5b00255
    82. 82
      Galindo-Murillo, R.; Robertson, J. C.; Zgarbová, M.; Šponer, J.; Otyepka, M.; Jurečka, P.; Cheatham, T. E. Assessing the Current State of Amber Force Field Modifications for DNA. J. Chem. Theory Comput. 2016, 12 (8), 41144127,  DOI: 10.1021/acs.jctc.6b00186
    83. 83
      Dans, P. D.; Ivani, I.; Hospital, A.; Portella, G.; González, C.; Orozco, M. How Accurate Are Accurate Force-Fields for B-DNA?. Nucleic Acids Res. 2017, 45 (7), 42174230,  DOI: 10.1093/nar/gkw1355
    84. 84
      Wang, J.; Wolf, R. M.; Caldwell, J. W.; Kollman, P. A.; Case, D. A. Development and Testing of a General Amber Force Field. J. Comput. Chem. 2004, 25 (9), 11571174,  DOI: 10.1002/jcc.20035
    85. 85
      Singh, U. C.; Kollman, P. A. An Approach to Computing Electrostatic Charges for Molecules. J. Comput. Chem. 1984, 5 (2), 129145,  DOI: 10.1002/jcc.540050204
    86. 86
      Bayly, C. I.; Cieplak, P.; Cornell, W.; Kollman, P. A. A Well-Behaved Electrostatic Potential Based Method Using Charge Restraints for Deriving Atomic Charges: The RESP Model. J. Phys. Chem. 1993, 97 (40), 1026910280,  DOI: 10.1021/j100142a004
    87. 87
      Hess, B. P-LINCS: A Parallel Linear Constraint Solver for Molecular Simulation. J. Chem. Theory Comput. 2008, 4 (1), 116122,  DOI: 10.1021/ct700200b
    88. 88
      Abraham, M. J.; Murtola, T.; Schulz, R.; Pall, S.; Smith, J. C.; Hess, B.; Lindahl, E. Gromacs: High Performance Molecular Simulations through Multi-Level Parallelism from Laptops to Supercomputers. SoftwareX 2015, 1–2, 1925,  DOI: 10.1016/j.softx.2015.06.001
    89. 89
      Darden, T.; York, D.; Pedersen, L. Particle Mesh Ewald: An N·log(N) Method for Ewald Sums in Large Systems. J. Chem. Phys. 1993, 98 (12), 1008910092,  DOI: 10.1063/1.464397
    90. 90
      Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A Smooth Particle Mesh Ewald Method. J. Chem. Phys. 1995, 103 (19), 85778593,  DOI: 10.1063/1.470117
    91. 91
      Matta, C. F.; Bader, R. F. W. Atoms-in-molecules study of the genetically encoded amino acids. III. Bond and atomic properties and their correlations with experiment including mutation-induced changes in protein stability and genetic coding. Proteins: Struct., Funct., Genet. 2003, 52 (3), 360399,  DOI: 10.1002/prot.10414
    92. 92
      Peraro, M. D.; Spiegel, K.; Lamoureux, G.; De Vivo, M.; DeGrado, W. F.; Klein, M. L. Modeling the Charge Distribution at Metal Sites in Proteins for Molecular Dynamics Simulations. J. Struct. Biol. 2007, 157 (3), 444453,  DOI: 10.1016/j.jsb.2006.10.019
    93. 93
      Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of Simple Potential Functions for Simulating Liquid Water. J. Chem. Phys. 1983, 79 (2), 926935,  DOI: 10.1063/1.445869
    94. 94
      Bussi, G.; Donadio, D.; Parrinello, M. Canonical Sampling through Velocity Rescaling. J. Chem. Phys. 2007, 126 (1), 014101,  DOI: 10.1063/1.2408420
    95. 95
      Parrinello, M.; Rahman, A. Polymorphic Transitions in Single Crystals: A New Molecular Dynamics Method. J. Appl. Phys. 1981, 52 (12), 71827190,  DOI: 10.1063/1.328693
    96. 96
      Barducci, A.; Bussi, G.; Parrinello, M. Well-Tempered Metadynamics: A Smoothly Converging and Tunable Free-Energy Method. Phys. Rev. Lett. 2008, 100 (2), 020603  DOI: 10.1103/PhysRevLett.100.020603
  • Supporting Information

    Supporting Information


    The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/jacs.9b10656.

    • Setup protocols of the force-field-based MD simulations; details on enhanced sampling free energy calculations and additional analyses; PDB analyses, structure, and sequence alignments (PDF)

    • Illustrative movie of the overall catalytic process of hExo1, created by merging PDB structures with fragments of trajectories of equilibrium and metadynamics simulations (MP4)


    Terms & Conditions

    Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.