Modeling DNA polymerase μ motions: subtle transitions before chemistry

Time-lapse crystallography snapshots of a double-strand break repair polymerase in action

DNA polymerase (pol) μ is a DNA-dependent polymerase that incorporates nucleotides during gap-filling synthesis in the non-homologous end-joining pathway of double-strand break repair. Here we report time-lapse X-ray crystallography snapshots of catalytic events during gap-filling DNA synthesis by pol μ. Unique catalytic intermediates and active site conformational changes that underlie catalysis are uncovered, and a transient third (product) metal ion is observed in the product state. The product manganese coordinates phosphate oxygens of the inserted nucleotide and PPi. The product metal is not observed during DNA synthesis in the presence of magnesium. Kinetic analyses indicate that manganese increases the rate constant for deoxynucleoside 5'-triphosphate insertion compared to magnesium. The likely product stabilization role of the manganese product metal in pol μ is discussed. These observations provide insight on structural attributes of this X-family double-strand break repair polymerase that impact its biological function in genome maintenance.DNA polymerase (pol) μ functions in DNA double-strand break repair. Here the authors use time-lapse X-ray crystallography to capture the states of pol µ during the conversion from pre-catalytic to product complex and observe a third transiently bound metal ion in the product state.

Molecular Dynamics Simulations and Structural Analysis to Decipher Functional Impact of a Twenty Residue Insert in the Ternary Complex of Mus musculus TdT Isoform

Insertions/deletions are common evolutionary tools employed to alter the structural and functional repertoire of protein domains. An insert situated proximal to the active site or ligand binding site frequently impacts protein function; however, the effect of distal indels on protein activity and/or stability are often not studied. In this paper, we have investigated a distal insert, which influences the function and stability of a unique DNA polymerase, called terminal deoxynucleotidyl transferase (TdT). TdT (EC:2.7.7.31) is a monomeric 58 kDa protein belonging to family X of eukaryotic DNA polymerases and known for its role in V(D)J recombination as well as in non-homologous end-joining (NHEJ) pathways. Two murine isoforms of TdT, with a length difference of twenty residues and having different biochemical properties, have been studied. All-atom molecular dynamics simulations at different temperatures and interaction network analyses were performed on the short and long-length isoforms. We observed conformational changes in the regions distal to the insert position (thumb subdomain) in the longer isoform, which indirectly affects the activity and stability of the enzyme through a mediating loop (Loop1). A structural rationale could be provided to explain the reduced polymerization rate as well as increased thermosensitivity of the longer isoform caused by peripherally located length variations within a DNA polymerase. These observations increase our understanding of the roles of length variants in introducing functional diversity in protein families in general.

Structural basis for a novel mechanism of DNA bridging and alignment in eukaryotic DSB DNA repair

Eukaryotic DNA polymerase mu of the PolX family can promote the association of the two 3'-protruding ends of a DNA double-strand break (DSB) being repaired (DNA synapsis) even in the absence of the core non-homologous end-joining (NHEJ) machinery. Here, we show that terminal deoxynucleotidyltransferase (TdT), a closely related PolX involved in V(D)J recombination, has the same property. We solved its crystal structure with an annealed DNA synapsis containing one micro-homology (MH) base pair and one nascent base pair. This structure reveals how the N-terminal domain and Loop 1 of Tdt cooperate for bridging the two DNA ends, providing a templating base in trans and limiting the MH search region to only two base pairs. A network of ordered water molecules is proposed to assist the incorporation of any nucleotide independently of the in trans templating base. These data are consistent with a recent model that explains the statistics of sequences synthesized in vivo by Tdt based solely on this dinucleotide step. Site-directed mutagenesis and functional tests suggest that this structural model is also valid for Pol mu during NHEJ.

Molecular dynamics study of the opening mechanism for DNA polymerase I

During DNA replication, DNA polymerases follow an induced fit mechanism in order to rapidly distinguish between correct and incorrect dNTP substrates. The dynamics of this process are crucial to the overall effectiveness of catalysis. Although X-ray crystal structures of DNA polymerase I with substrate dNTPs have revealed key structural states along the catalytic pathway, solution fluorescence studies indicate that those key states are populated in the absence of substrate. Herein, we report the first atomistic simulations showing the conformational changes between the closed, open, and ajar conformations of DNA polymerase I in the binary (enzyme:DNA) state to better understand its dynamics. We have applied long time-scale, unbiased molecular dynamics to investigate the opening process of the fingers domain in the absence of substrate for B. stearothermophilis DNA polymerase in silico. These simulations are biologically and/or physiologically relevant as they shed light on the transitions between states in this important enzyme. All closed and ajar simulations successfully transitioned into the fully open conformation, which is known to be the dominant binary enzyme-DNA conformation from solution and crystallographic studies. Furthermore, we have detailed the key stages in the opening process starting from the open and ajar crystal structures, including the observation of a previously unknown key intermediate structure. Four backbone dihedrals were identified as important during the opening process, and their movements provide insight into the recognition of dNTP substrate molecules by the polymerase binary state. In addition to revealing the opening mechanism, this study also demonstrates our ability to study biological events of DNA polymerase using current computational methods without biasing the dynamics.

Decision-making during NHEJ: a network of interactions in human Polμ implicated in substrate recognition and end-bridging

Human Polμ is a DNA polymerase belonging to the X family that has been implicated in the non-homologous end-joining (NHEJ) pathway during repair of double-strand breaks in DNA. Loop1 is a flexible piece of Polμ which has a critical role during terminal transferase and end-joining activities: it acts as a pseudo-template when the template strand is discontinuous or unavailable, whilst diffusing away if present to avoid steric clashes. Mutational analysis and inspection of the 3D structures available allowed us to identify a network of residues in charge of sensing the presence or absence of discontinuities in the template strand, which will in turn determine the final position adopted by Loop1. This network is formed by the previously uncharacterized thumb mini-loop (NSH motif) and the positively charged helix N, which contribute to the correct positioning of Loop1 and to juxtapose the discontinuous template strand during NHEJ of incompatible ends. Accordingly, single mutation of specific conserved residues in these motifs, whilst irrelevant in most of the cases for gap filling, largely affected terminal transferase and end-joining activities. Other point mutations in the ‘hinges’ of Loop1, such as residues Phe³⁸⁵ or Phe³⁸⁹, corroborated the flexibility requirements of this motif.

C(α) torsion angles as a flexible criterion to extract secrets from a molecular dynamics simulation

Given the increasing complexity of simulated molecular systems, and the fact that simulation times have now reached milliseconds to seconds, immense amounts of data (in the gigabyte to terabyte range) are produced in current molecular dynamics simulations. Manual analysis of these data is a very time-consuming task, and important events that lead from one intermediate structure to another can become occluded in the noise resulting from random thermal fluctuations. To overcome these problems and facilitate a semi-automated data analysis, we introduce in this work a measure based on C(α) torsion angles: torsion angles formed by four consecutive C(α) atoms. This measure describes changes in the backbones of large systems on a residual length scale (i.e., a small number of residues at a time). Cluster analysis of individual C(α) torsion angles and its fuzzification led to continuous time patches representing (meta)stable conformations and to the identification of events acting as transitions between these conformations. The importance of a change in torsion angle to structural integrity is assessed by comparing this change to the average fluctuations in the same torsion angle over the complete simulation. Using this novel measure in combination with other measures such as the root mean square deviation (RMSD) and time series of distance measures, we performed an in-depth analysis of a simulation of the open form of DNA polymerase I. The times at which major conformational changes occur and the most important parts of the molecule and their interrelations were pinpointed in this analysis. The simultaneous determination of the time points and localizations of major events is a significant advantage of the new bottom-up approach presented here, as compared to many other (top-down) approaches in which only the similarity of the complete structure is analyzed.

Structures of the Leishmania infantum polymerase beta

Protozoans of the genus Leishmania, the pathogenic agent causing leishmaniasis, encode the family X DNA polymerase Li Pol β. Here, we report the first crystal structures of Li Pol β. Our pre- and post-catalytic structures show that the polymerase adopts the common family X DNA polymerase fold. However, in contrast to other family X DNA polymerases, the dNTP-induced conformational changes in Li Pol β are much more subtle. Moreover, pre- and post-catalytic structures reveal that Li Pol β interacts with the template strand through a nonconserved, variable region known as loop3. Li Pol β Δloop3 mutants display a higher catalytic rate, catalytic efficiency and overall error rates with respect to WT Li Pol β. These results further demonstrate the subtle structural variability that exists within this family of enzymes and provides insight into how this variability underlies the substantial functional differences among their members.

Sustained active site rigidity during synthesis by human DNA polymerase μ

DNA polymerase μ (Pol μ) is the only template-dependent human DNA polymerase capable of repairing double-strand DNA breaks (DSBs) with unpaired 3' ends in nonhomologous end joining (NHEJ). To probe this function, we structurally characterized Pol μ's catalytic cycle for single-nucleotide incorporation. These structures indicate that, unlike other template-dependent DNA polymerases, Pol μ shows no large-scale conformational changes in protein subdomains, amino acid side chains or DNA upon dNTP binding or catalysis. Instead, the only major conformational change is seen earlier in the catalytic cycle, when the flexible loop 1 region repositions upon DNA binding. Pol μ variants with changes in loop 1 have altered catalytic properties and are partially defective in NHEJ. The results indicate that specific loop 1 residues contribute to Pol μ's unique ability to catalyze template-dependent NHEJ of DSBs with unpaired 3' ends.

A specific N-terminal extension of the 8 kDa domain is required for DNA end-bridging by human Polμ and Polλ

Human DNA polymerases mu (Polµ) and lambda (Polλ) are X family members involved in the repair of double-strand breaks in DNA during non-homologous end joining. Crucial abilities of these enzymes include bridging of the two 3′ single-stranded overhangs and trans-polymerization using one 3′ end as primer and the other as template, to minimize sequence loss. In this context, we have studied the importance of a previously uncharacterised sequence (‘brooch’), located at the N-terminal boundary of the Polß-like polymerase core, and formed by Tyr¹⁴¹, Ala¹⁴², Cys¹⁴³, Gln¹⁴⁴ and Arg¹⁴⁵ in Polµ, and by Trp²³⁹, Val²⁴⁰, Cys²⁴¹, Ala²⁴² and Gln²⁴³ in Polλ. The brooch is potentially implicated in the maintenance of a closed conformation throughout the catalytic cycle, and our studies indicate that it could be a target of Cdk phosphorylation in Polµ. The brooch is irrelevant for 1 nt gap filling, but of specific importance during end joining: single mutations in the conserved residues reduced the formation of two ended synapses and strongly diminished the ability of Polµ and polymerase lambda to perform non-homologous end joining reactions in vitro.

Structures of intermediates along the catalytic cycle of terminal deoxynucleotidyltransferase: dynamical aspects of the two-metal ion mechanism

Terminal deoxynucleotidyltransferase (Tdt) is a non-templated eukaryotic DNA polymerase of the polX family that is responsible for the random addition of nucleotides at the V(D)J junctions of immunoglobulins and T-cell receptors. Here we describe a series of high-resolution X-ray structures that mimic the pre-catalytic state, the post-catalytic state and a competent state that can be transformed into the two other ones in crystallo via the addition of dAMPcPP and Zn(2+), respectively. We examined the effect of Mn(2+), Co(2+) and Zn(2+) because they all have a marked influence on the kinetics of the reaction. We demonstrate a dynamic role of divalent transition metal ions bound to site A: (i) Zn(2+) (or Co(2+)) in Metal A site changes coordination from octahedral to tetrahedral after the chemical step, which explains the known higher affinity of Tdt for the primer strand when these ions are present, and (ii) metal A has to leave to allow the translocation of the primer strand and to clear the active site, a typical feature for a ratchet-like mechanism. Except for Zn(2+), the sugar puckering of the primer strand 3' terminus changes from C2'-endo to C3'-endo during catalysis. In addition, our data are compatible with a scheme where metal A is the last component that binds to the active site to complete its productive assembly, as already inferred in human pol beta. The new structures have potential implications for modeling pol mu, a closely related polX implicated in the repair of DNA double-strand breaks, in a complex with a DNA synapsis.

"Gate-keeper" residues and active-site rearrangements in DNA polymerase μ help discriminate non-cognate nucleotides

Incorporating the cognate instead of non-cognate substrates is crucial for DNA polymerase function. Here we analyze molecular dynamics simulations of DNA polymerase μ (pol μ) bound to different non-cognate incoming nucleotides including A:dCTP, A:dGTP, A(syn):dGTP, A:dATP, A(syn):dATP, T:dCTP, and T:dGTP to study the structure-function relationships involved with aberrant base pairs in the conformational pathway; while a pol μ complex with the A:dTTP base pair is available, no solved non-cognate structures are available. We observe distinct differences of the non-cognate systems compared to the cognate system. Specifically, the motions of active-site residue His329 and Asp330 distort the active site, and Trp436, Gln440, Glu443 and Arg444 tend to tighten the nucleotide-binding pocket when non-cognate nucleotides are bound; the latter effect may further lead to an altered electrostatic potential within the active site. That most of these “gate-keeper” residues are located farther apart from the upstream primer in pol μ, compared to other X family members, also suggests an interesting relation to pol μ's ability to incorporate nucleotides when the upstream primer is not paired. By examining the correlated motions within pol μ complexes, we also observe different patterns of correlations between non-cognate systems and the cognate system, especially decreased interactions between the incoming nucleotides and the nucleotide-binding pocket. Altered correlated motions in non-cognate systems agree with our recently proposed hybrid conformational selection/induced-fit models. Taken together, our studies propose the following order for difficulty of non-cognate system insertions by pol μ: T:dGTP<A(syn):dATP<T:dCTP<A:dGTP<A(syn):dGTP<A:dCTP<A:dATP. This sequence agrees with available kinetic data for non-cognate nucleotide insertions, with the exception of A:dGTP, which may be more sensitive to the template sequence. The structures and conformational aspects predicted here are experimentally testable.

DNA polymerase μ (pol μ) is an enzyme that participates in DNA repair and thus has a central role in maintaining the integrity of genetic information. To efficiently repair the DNA, discriminating the cognate instead of non-cognate nucleotides (“fidelity-checking”) is required. Here we analyze molecular dynamics simulations of pol μ bound to different non-cognate nucleotides to study the structure-function relationships involved in the fidelity-checking mechanism of pol μ on the atomic level. Our results suggest that His329, Asp330, Trp436, Gln440, Glu443, and Arg444 are of great importance for pol μ's fidelity-checking mechanism. We also observe altered patterns of correlated motions within pol μ complex when non-cognate instead of cognate nucleotides are bound, which agrees with our recently proposed hybrid conformational selection/induced-fit models. Taken together, our studies help interpret the available kinetic data of various non-cognate nucleotide insertions by pol μ. We also suggest experimentally testable predictions; for example, a point mutation like E443M may reduce the ability of pol μ to insert the cognate more than of non-cognate nucleotides. Our studies suggest an interesting relation to pol μ's unique ability to incorporate nucleotides when the upstream primer is not paired.