Modeling of flap endonuclease interactions with DNA substrate

Implicit Solvent with Explicit Ions Generalized Born Model in Molecular Dynamics: Application to DNA

The ion atmosphere surrounding highly charged biomolecules, such as nucleic acids, is crucial for their dynamics, structure, and interactions. Here, we develop an approach for the explicit treatment of ions within an implicit solvent framework suitable for atomistic simulations of biomolecules. The proposed implicit solvent/explicit ions model, GBION, is based on a modified generalized Born (GB) model; it includes separate, modified GB terms for solute-ion and ion-ion interactions. The model is implemented in the AMBER package (version 24), and its performance is thoroughly investigated in atomistic molecular dynamics (MD) simulations of double-stranded DNA on a microsecond time scale. The aggregate characteristics of monovalent (Na+ and K+) and trivalent (Cobalt Hexammine, CoHex3+) counterion distributions around double-stranded DNA predicted by the model are in reasonable agreement with the experiment (where available), all-atom explicit water MD simulations, and the expectation from the Manning condensation theory. The radial distributions of monovalent cations around DNA are reasonably close to the ones obtained using the explicit water model: expressed in units of energy, the maximum deviations of local ion concentrations from the explicit solvent reference are within 1 kBT, comparable to the corresponding deviations expected between different established explicit water models. The proposed GBION model is able to simulate DNA fragments in a large volume of solvent with explicit ions with little additional computational overhead compared with the fully implicit GB treatment of ions. Ions simulated using the developed model explore conformational space at least 2 orders of magnitude faster than in the explicit solvent. These advantages allowed us to observe and explore an unexpected "stacking" mode of DNA condensation in the presence of trivalent counterions (CoHex3+) that was revealed by recent experiments.

Biochemical characterization and mutational analysis of a novel flap endonuclease 1 from Thermococcus barophilus Ch5

Flap endonuclease 1 (FEN1) plays important roles in DNA replication, repair and recombination. Herein, we report biochemical characteristics and catalytic mechanism of a novel FEN1 from the hyperthermophilic euryarchaeon Thermococcus barophilus Ch5 (Tb-FEN1). As expected, the recombinant Tb-FEN1 can cleave 5'-flap DNA. However, the enzyme has no activity on cleaving pseudo Y DNA, which sharply contrasts with other archaeal and eukaryotic FEN1 homologs. Tb-FEN1 retains 24% relative activity after heating at 100 °C for 20 min, demonstrating that it is the most thermostable among all reported FEN1 proteins. The enzyme displays maximal activity in a wide range of pH from 7.0 to 9.5. The Tb-FEN1 activity is dependent on a divalent metal ion, among which Mg²⁺ and Mn²⁺ are optimal. Enzyme activity is inhibited by NaCl. Kinetic analyzes estimated that an activation energy for removal of 5'-flap from DNA by Tb-FEN1 was 35.7 ± 4.3 kcal/mol, which is the first report on energy barrier for excising 5'-flap from DNA by a FEN1 enzyme. Mutational studies demonstrate that the K87A, R94A and E154A amino acid substitutions abolish cleavage activity and reduce 5'-flap DNA binding efficiencies, suggesting that residues K87, R94, and E154 in Tb-FEN1 are essential for catalysis and DNA binding as well. Overall, Tb-FEN1 is an extremely thermostable endonuclease with unusual features.

Gene cloning and characterization of Tk1281, a flap endonuclease 1 from Thermococcus kodakarensis

Flap endonuclease is a structure-specific nuclease which cleaves 5'-flap of bifurcated DNA substrates. Genome sequence of Thermococcus kodakarensis harbors an open reading frame, Tk1281, exhibiting high homology with archaeal flap endonucleases 1. The corresponding gene was cloned and expressed in Escherichia coli, and the gene product was purified to apparent homogeneity. Tk1281 was a monomer of 38 kDa and catalyzed the cleavage of 5'-flap from double-stranded DNA substrate containing single-stranded DNA flap. The highest cleavage activity was observed at 80 °C and pH 7.5. Under optimal conditions, Tk1281 exhibited apparent V_max and K_m values of 278 nmol/min/mg and 37 μM, respectively, against a 54-nucleotide double-stranded substrate containing a single-stranded 5'-flap of 27 nucleotides. A unique feature of Tk1281 is its highest activation in the presence of Co²⁺ and no activation with Mn²⁺. To the best of our knowledge, this is the first cloning and characterization of a flap endonuclease from the genus Thermococcus.

Positioning the 5'-flap junction in the active site controls the rate of flap endonuclease-1-catalyzed DNA cleavage

Flap endonucleases catalyze cleavage of single-stranded DNA flaps formed during replication, repair, and recombination and are therefore essential for genome processing and stability. Recent crystal structures of DNA-bound human flap endonuclease (hFEN1) offer new insights into how conformational changes in the DNA and hFEN1 may facilitate the reaction mechanism. For example, previous biochemical studies of DNA conformation performed under non-catalytic conditions with Ca²⁺ have suggested that base unpairing at the 5'-flap:template junction is an important step in the reaction, but the new structural data suggest otherwise. To clarify the role of DNA changes in the kinetic mechanism, we measured a series of transient steps, from substrate binding to product release, during the hFEN1-catalyzed reaction in the presence of Mg²⁺ We found that whereas hFEN1 binds and bends DNA at a fast, diffusion-limited rate, much slower Mg²⁺-dependent conformational changes in DNA around the active site are subsequently necessary and rate-limiting for 5'-flap cleavage. These changes are reported overall by fluorescence of 2-aminopurine at the 5'-flap:template junction, indicating that local DNA distortion (e.g. disruption of base stacking observed in structures), associated with positioning the 5'-flap scissile phosphodiester bond in the hFEN1 active site, controls catalysis. hFEN1 residues with distinct roles in the catalytic mechanism, including those binding metal ions (Asp-34 and Asp-181), steering the 5'-flap through the active site and binding the scissile phosphate (Lys-93 and Arg-100), and stacking against the base 5' to the scissile phosphate (Tyr-40), all contribute to these rate-limiting conformational changes, ensuring efficient and specific cleavage of 5'-flaps.

DNA and Protein Requirements for Substrate Conformational Changes Necessary for Human Flap Endonuclease-1-catalyzed Reaction

Human flap endonuclease-1 (hFEN1) catalyzes the essential removal of single-stranded flaps arising at DNA junctions during replication and repair processes. hFEN1 biological function must be precisely controlled, and consequently, the protein relies on a combination of protein and substrate conformational changes as a prerequisite for reaction. These include substrate bending at the duplex-duplex junction and transfer of unpaired reacting duplex end into the active site. When present, 5'-flaps are thought to thread under the helical cap, limiting reaction to flaps with free 5'-terminiin vivo Here we monitored DNA bending by FRET and DNA unpairing using 2-aminopurine exciton pair CD to determine the DNA and protein requirements for these substrate conformational changes. Binding of DNA to hFEN1 in a bent conformation occurred independently of 5'-flap accommodation and did not require active site metal ions or the presence of conserved active site residues. More stringent requirements exist for transfer of the substrate to the active site. Placement of the scissile phosphate diester in the active site required the presence of divalent metal ions, a free 5'-flap (if present), a Watson-Crick base pair at the terminus of the reacting duplex, and the intact secondary structure of the enzyme helical cap. Optimal positioning of the scissile phosphate additionally required active site conserved residues Tyr(40), Asp(181), and Arg(100)and a reacting duplex 5'-phosphate. These studies suggest a FEN1 reaction mechanism where junctions are bound and 5'-flaps are threaded (when present), and finally the substrate is transferred onto active site metals initiating cleavage.

Real-time PCR probe optimization using design of experiments approach

Primer and probe sequence designs are among the most critical input factors in real-time polymerase chain reaction (PCR) assay optimization. In this study, we present the use of statistical design of experiments (DOE) approach as a general guideline for probe optimization and more specifically focus on design optimization of label-free hydrolysis probes that are designated as mediator probes (MPs), which are used in reverse transcription MP PCR (RT-MP PCR). The effect of three input factors on assay performance was investigated: distance between primer and mediator probe cleavage site; dimer stability of MP and target sequence (influenza B virus); and dimer stability of the mediator and universal reporter (UR). The results indicated that the latter dimer stability had the greatest influence on assay performance, with RT-MP PCR efficiency increased by up to 10% with changes to this input factor. With an optimal design configuration, a detection limit of 3-14 target copies/10 μl reaction could be achieved. This improved detection limit was confirmed for another UR design and for a second target sequence, human metapneumovirus, with 7-11 copies/10 μl reaction detected in an optimum case. The DOE approach for improving oligonucleotide designs for real-time PCR not only produces excellent results but may also reduce the number of experiments that need to be performed, thus reducing costs and experimental times.

Rate-determining Step of Flap Endonuclease 1 (FEN1) Reflects a Kinetic Bias against Long Flaps and Trinucleotide Repeat Sequences

Flap endonuclease 1 (FEN1) is a structure-specific nuclease responsible for removing 5'-flaps formed during Okazaki fragment maturation and long patch base excision repair. In this work, we use rapid quench flow techniques to examine the rates of 5'-flap removal on DNA substrates of varying length and sequence. Of particular interest are flaps containing trinucleotide repeats (TNR), which have been proposed to affect FEN1 activity and cause genetic instability. We report that FEN1 processes substrates containing flaps of 30 nucleotides or fewer at comparable single-turnover rates. However, for flaps longer than 30 nucleotides, FEN1 kinetically discriminates substrates based on flap length and flap sequence. In particular, FEN1 removes flaps containing TNR sequences at a rate slower than mixed sequence flaps of the same length. Furthermore, multiple-turnover kinetic analysis reveals that the rate-determining step of FEN1 switches as a function of flap length from product release to chemistry (or a step prior to chemistry). These results provide a kinetic perspective on the role of FEN1 in DNA replication and repair and contribute to our understanding of FEN1 in mediating genetic instability of TNR sequences.

Introducing Charge Hydration Asymmetry into the Generalized Born Model

The effect of charge hydration asymmetry (CHA)-non-invariance of solvation free energy upon solute charge inversion-is missing from the standard linear response continuum electrostatics. The proposed charge hydration asymmetric-generalized Born (CHA-GB) approximation introduces this effect into the popular generalized Born (GB) model. The CHA is added to the GB equation via an analytical correction that quantifies the specific propensity of CHA of a given water model; the latter is determined by the charge distribution within the water model. Significant variations in CHA seen in explicit water (TIP3P, TIP4P-Ew, and TIP5P-E) free energy calculations on charge-inverted "molecular bracelets" are closely reproduced by CHA-GB, with the accuracy similar to models such as SEA and 3D-RISM that go beyond the linear response. Compared against reference explicit (TIP3P) electrostatic solvation free energies, CHA-GB shows about a 40% improvement in accuracy over the canonical GB, tested on a diverse set of 248 rigid small neutral molecules (root mean square error, rmse = 0.88 kcal/mol for CHA-GB vs 1.24 kcal/mol for GB) and 48 conformations of amino acid analogs (rmse = 0.81 kcal/mol vs 1.26 kcal/mol). CHA-GB employs a novel definition of the dielectric boundary that does not subsume the CHA effects into the intrinsic atomic radii. The strategy leads to finding a new set of intrinsic atomic radii optimized for CHA-GB; these radii show physically meaningful variation with the atom type, in contrast to the radii set optimized for GB. Compared to several popular radii sets used with the original GB model, the new radii set shows better transferability between different classes of molecules.

VERIFICATION OF THE GENERALIZED BORN MODEL AT SHORT DISTANCES

The generalized Born (GB) model, one of the implicit solvent models, is widely applied in molecular dynamics (MD) simulations as a simple description of the solvation effect. In the GB model, an empirical function called the Still's formula, with the algorithmic simplicity, is utilized to calculate the solvation energy due to the polarization, termed as ΔG pol . Applications of the GB model have exhibited reasonable accuracy and high computational efficiency. However, there is still room for improvements. Most of the attempts to improve the GB model focus on optimizing effective Born radii. Contrarily, limited researches have been performed to improve the feasibility of the Still's formula. In this paper, analytical methods was applied to investigate the validity of the Still's formula at short distance. Taking advantage of the toroidal coordinates and Mehler–Fock transform, the analytical solutions of the GB model at short distances was derived explicitly for the first time. Additionally, the solvation energy was numerically computed using proper algorithms based on the analytical solutions and compared with ΔG pol calculated in the GB model. With the analysis on the deficiencies of the Still's formula at short distances, potential methods to improve the validity of the GB model were discussed.

Flap endonuclease activity of gene 6 exonuclease of bacteriophage T7

Flap endonucleases remove flap structures generated during DNA replication. Gene 6 protein of bacteriophage T7 is a 5'-3'-exonuclease specific for dsDNA. Here we show that gene 6 protein also possesses a structure-specific endonuclease activity similar to known flap endonucleases. The flap endonuclease activity is less active relative to its exonuclease activity. The major cleavage by the endonuclease activity occurs at a position one nucleotide into the duplex region adjacent to a dsDNA-ssDNA junction. The efficiency of cleavage of the flap decreases with increasing length of the 5'-overhang. A 3'-single-stranded tail arising from the same end of the duplex as the 5'-tail inhibits gene 6 protein flap endonuclease activity. The released flap is not degraded further, but the exonuclease activity then proceeds to hydrolyze the 5'-terminal strand of the duplex. T7 gene 2.5 single-stranded DNA-binding protein stimulates the exonuclease and also the endonuclease activity. This stimulation is attributed to a specific interaction between the two proteins because Escherichia coli single-stranded DNA binding protein does not produce this stimulatory effect. The ability of gene 6 protein to remove 5'-terminal overhangs as well as to remove nucleotides from the 5'-termini enables it to effectively process the 5'-termini of Okazaki fragments before they are ligated.

Heat conductivity of DNA double helix

Thermal conductivity of isolated single molecule DNA fragments is of importance for nanotechnology, but has not yet been measured experimentally. Theoretical estimates based on simplified (1D) models predict anomalously high thermal conductivity. To investigate thermal properties of single molecule DNA we have developed a 3D coarse-grained (CG) model that retains the realism of the full all-atom description, but is significantly more efficient. Within the proposed model each nucleotide is represented by 6 particles or grains; the grains interact via effective potentials inferred from classical molecular dynamics (MD) trajectories based on a well-established all-atom potential function. Comparisons of 10 ns long MD trajectories between the CG and the corresponding all-atom model show similar root-mean-square deviations from the canonical B-form DNA, and similar structural fluctuations. At the same time, the CG model is 10 to 100 times faster depending on the length of the DNA fragment in the simulation. Analysis of dispersion curves derived from the CG model yields longitudinal sound velocity and torsional stiffness in close agreement with existing experiments. The computational efficiency of the CG model makes it possible to calculate thermal conductivity of a single DNA molecule not yet available experimentally. For a uniform (polyG-polyC) DNA, the estimated conductivity coefficient is 0.3 W/mK which is half the value of thermal conductivity for water. This result is in stark contrast with estimates of thermal conductivity for simplified, effectively 1D chains ("beads on a spring") that predict anomalous (infinite) thermal conductivity. Thus, full 3D character of DNA double-helix retained in the proposed model appears to be essential for describing its thermal properties at a single molecule level.

Structures of human exonuclease 1 DNA complexes suggest a unified mechanism for nuclease family

Human exonuclease 1 (hExo1) plays important roles in DNA repair and recombination processes that maintain genomic integrity. It is a member of the 5' structure-specific nuclease family of exonucleases and endonucleases that includes FEN-1, XPG, and GEN1. We present structures of hExo1 in complex with a DNA substrate, followed by mutagenesis studies, and propose a common mechanism by which this nuclease family recognizes and processes diverse DNA structures. hExo1 induces a sharp bend in the DNA at nicks or gaps. Frayed 5' ends of nicked duplexes resemble flap junctions, unifying the mechanisms of endo- and exonucleolytic processing. Conformational control of a mobile region in the catalytic site suggests a mechanism for allosteric regulation by binding to protein partners. The relative arrangement of substrate binding sites in these enzymes provides an elegant solution to a complex geometrical puzzle of substrate recognition and processing.

Structure of flap endonuclease 1 from the hyperthermophilic archaeon Desulfurococcus amylolyticus

Flap endonuclease 1 (FEN1) is a key enzyme in DNA repair and DNA replication. It is a structure-specific nuclease that removes 5'-overhanging flaps and the RNA/DNA primer during maturation of the Okazaki fragment. Homologues of FEN1 exist in a wide range of bacteria, archaea and eukaryotes. In order to further understand the structural basis of the DNA recognition, binding and cleavage mechanism of FEN1, the structure of FEN1 from the hyperthermophilic archaeon Desulfurococcus amylolyticus (DaFEN1) was determined at 2.00 Å resolution. The overall fold of DaFEN1 was similar to those of other archaeal FEN1 proteins; however, the helical clamp and the flexible loop exhibited a putative substrate-binding pocket with a unique conformation.

Substrate recognition and catalysis by flap endonucleases and related enzymes

FENs (flap endonucleases) and related FEN-like enzymes [EXO-1 (exonuclease-1), GEN-1 (gap endonuclease 1) and XPG (xeroderma pigmentosum complementation group G)] are a family of bivalent-metal-ion-dependent nucleases that catalyse structure-specific hydrolysis of DNA duplex-containing nucleic acid structures during DNA replication, repair and recombination. In the case of FENs, the ability to catalyse reactions on a variety of substrates has been rationalized as a result of combined functional and structural studies. Analyses of FENs also exemplify controversies regarding the two-metal-ion mechanism. However, kinetic studies of T5FEN (bacteriophage T5 FEN) reveal that a two-metal-ion-like mechanism for chemical catalysis is plausible. Consideration of the metallobiochemistry and the positioning of substrate in metal-free structures has led to the proposal that the duplex termini of substrates are unpaired in the catalytically active form and that FENs and related enzymes may recognize breathing duplex termini within more complex structures. An outstanding issue in FEN catalysis is the role played by the intermediate (I) domain arch or clamp. It has been proposed that FENs thread the 5′-portion of their substrates through this arch, which is wide enough to accommodate single-stranded, but not double-stranded, DNA. However, FENs exhibit gap endonuclease activity acting upon substrates that have a region of 5′-duplex. Moreover, the action of other FEN family members such as GEN-1, proposed to target Holliday junctions without termini, appears incompatible with a threading mechanism. An alterative is that the I domain is used as a clamp. A future challenge is to clarify the role of this domain in FENs and related enzymes.

New approach to real-time nucleic acids detection: folding polymerase chain reaction amplicons into a secondary structure to improve cleavage of Forster resonance energy transfer probes in 5'-nuclease assays

The article describes a new technology for real-time polymerase chain reaction (PCR) detection of nucleic acids. Similar to Taqman, this new method, named Snake, utilizes the 5′-nuclease activity of Thermus aquaticus (Taq) DNA polymerase that cleaves dual-labeled Förster resonance energy transfer (FRET) probes and generates a fluorescent signal during PCR. However, the mechanism of the probe cleavage in Snake is different. In this assay, PCR amplicons fold into stem–loop secondary structures. Hybridization of FRET probes to one of these structures leads to the formation of optimal substrates for the 5′-nuclease activity of Taq. The stem–loop structures in the Snake amplicons are introduced by the unique design of one of the PCR primers, which carries a special 5′-flap sequence. It was found that at a certain length of these 5′-flap sequences the folded Snake amplicons have very little, if any, effect on PCR yield but benefit many aspects of the detection process, particularly the signal productivity. Unlike Taqman, the Snake system favors the use of short FRET probes with improved fluorescence background. The head-to-head comparison study of Snake and Taqman revealed that these two technologies have more differences than similarities with respect to their responses to changes in PCR protocol, e.g. the variations in primer concentration, annealing time, PCR asymmetry. The optimal PCR protocol for Snake has been identified. The technology’s real-time performance was compared to a number of conventional assays including Taqman, 3′-MGB-Taqman, Molecular Beacon and Scorpion primers. The test trial showed that Snake supersedes the conventional assays in the signal productivity and detection of sequence variations as small as single nucleotide polymorphisms. Due to the assay’s cost-effectiveness and simplicity of design, the technology is anticipated to quickly replace all known conventional methods currently used for real-time nucleic acid detection.

Active site substitutions delineate distinct classes of eubacterial flap endonuclease

FENs (flap endonucleases) play essential roles in DNA replication, pivotally in the resolution of Okazaki fragments. In eubacteria, DNA PolI (polymerase I) contains a flap processing domain, the N-terminal 5′→3′ exonuclease. We present evidence of paralogous FEN-encoding genes present in many eubacteria. Two distinct classes of these independent FEN-encoding genes exist with four groups of eubacteria, being identified based on the number and type of FEN gene encoded. The respective proteins possess distinct motifs hallmarking their differentiation. Crucially, based on primary sequence and predicted secondary structural motifs, we reveal key differences at their active sites. These results are supported by biochemical characterization of two family members - ExoIX (exonuclease IX) from Escherichia coli and SaFEN (Staphylococcus aureus FEN). These proteins displayed marked differences in their ability to process a range of branched and linear DNA structures. On bifurcated substrates, SaFEN exhibited similar substrate specificity to previously characterized FENs. In quantitative exonuclease assays, SaFEN maintained a comparable activity with that reported for PolI. However, ExoIX showed no observable enzymatic activity. A threaded model is presented for SaFEN, demonstrating the probable interaction of this newly identified class of FEN with divalent metal ions and a branched DNA substrate. The results from the present study provide an intriguing model for the cellular role of these FEN sub-classes and illustrate the evolutionary importance of processing aberrant DNA, which has led to their maintenance alongside DNA PolI in many eubacteria.

Crystallization and preliminary crystallographic analysis of the catalytic domain of human flap endonuclease 1 in complex with a nicked DNA product: use of a DPCS kit for efficient protein-DNA complex crystallization

Flap endonuclease 1 (FEN1) is a structure-specific nuclease that removes the RNA/DNA primer associated with Okazaki fragments in DNA replication. Here, crystals of the complex between the catalytic domain of human FEN1 and a DNA product have been obtained. For efficient crystallization screening, a DNA-protein complex crystallization screening (DPCS) kit was designed based on commercial crystallization kits. The crystal was found to belong to space group P2(1), with unit-cell parameters a = 61.0, b = 101.3, c = 106.4 A, beta = 106.4 degrees. The asymmetric unit is predicted to contain two complexes in the crystallographic asymmetric unit. A diffraction data set was collected to a resolution of 2.75 A.

Crystal structure of bacteriophage T4 5' nuclease in complex with a branched DNA reveals how flap endonuclease-1 family nucleases bind their substrates

Bacteriophage T4 RNase H, a flap endonuclease-1 family nuclease, removes RNA primers from lagging strand fragments. It has both 5' nuclease and flap endonuclease activities. Our previous structure of native T4 RNase H (PDB code 1TFR) revealed an active site composed of highly conserved Asp residues and two bound hydrated magnesium ions. Here, we report the crystal structure of T4 RNase H in complex with a fork DNA substrate bound in its active site. This is the first structure of a flap endonuclease-1 family protein with its complete branched substrate. The fork duplex interacts with an extended loop of the helix-hairpin-helix motif class 2. The 5' arm crosses over the active site, extending below the bridge (helical arch) region. Cleavage assays of this DNA substrate identify a primary cut site 7-bases in from the 5' arm. The scissile phosphate, the first bond in the duplex DNA adjacent to the 5' arm, lies above a magnesium binding site. The less ordered 3' arm reaches toward the C and N termini of the enzyme, which are binding sites for T4 32 protein and T4 45 clamp, respectively. In the crystal structure, the scissile bond is located within the double-stranded DNA, between the first two duplex nucleotides next to the 5' arm, and lies above a magnesium binding site. This complex provides important insight into substrate recognition and specificity of the flap endonuclease-1 enzymes.

Comparison of the catalytic parameters and reaction specificities of a phage and an archaeal flap endonuclease

Flap endonucleases (FENs) catalyse the exonucleolytic hydrolysis of blunt-ended duplex DNA substrates and the endonucleolytic cleavage of 5'-bifurcated nucleic acids at the junction formed between single and double-stranded DNA. The specificity and catalytic parameters of FENs derived from T5 bacteriophage and Archaeoglobus fulgidus were studied with a range of single oligonucleotide DNA substrates. These substrates contained one or more hairpin turns and mimic duplex, 5'-overhanging duplex, pseudo-Y, nicked DNA, and flap structures. The FEN-catalysed reaction properties of nicked DNA and flap structures possessing an extrahelical 3'-nucleotide (nt) were also characterised. The phage enzyme produced multiple reaction products of differing length with all the substrates tested, except when the length of duplex DNA downstream of the reaction site was truncated. Only larger DNAs containing two duplex regions are effective substrates for the archaeal enzyme and undergo reaction at multiple sites when they lack a 3'-extrahelical nucleotide. However, a single product corresponding to reaction 1 nt into the double-stranded region occurred with A. fulgidus FEN when substrates possessed a 3'-extrahelical nt. Steady-state and pre-steady-state catalytic parameters reveal that the phage enzyme is rate-limited by product release with all the substrates tested. Single-turnover maximal rates of reaction are similar with most substrates. In contrast, turnover numbers for T5FEN decrease as the size of the DNA substrate is increased. Comparison of the catalytic parameters of the A. fulgidus FEN employing flap and double-flap substrates indicates that binding interactions with the 3'-extrahelical nucleotide stabilise the ground state FEN-DNA interaction, leading to stimulation of comparative reactions at DNA concentrations below saturation with the single flap substrate. Maximal multiple turnover rates of the archaeal enzyme with flap and double flap substrates are similar. A model is proposed to account for the varying specificities of the two enzymes with regard to cleavage patterns and substrate preferences.

Examining interactions of HIV-1 reverse transcriptase with single-stranded template nucleotides by nucleoside analog interference

Crystallographic studies have implicated several residues of the p66 fingers subdomain of human immunodeficiency virus type-1 reverse transcriptase in contacting the single-stranded template overhang immediately ahead of the DNA polymerase catalytic center. This interaction presumably assists in inducing the appropriate geometry on the template base for efficient and accurate incorporation of the incoming dNTP. To investigate this, we introduced nucleoside analogs either individually or in tandem into the DNA template ahead of the catalytic center and investigated whether they induce pausing of the replication machinery before serving as the template base. Analogs included abasic tetrahydrofuran linkages, neutralizing methylphosphonate linkages, and conformationally locked nucleosides. In addition, several Phe-61 mutants were included in our analysis, based on previous data indicating that altering this residue affects both strand displacement synthesis and the fidelity of DNA synthesis. We demonstrate here that altering the topology of the template strand two nucleotides ahead of the catalytic center can interrupt DNA synthesis. Mutating Phe-61 to either Ala or Leu accentuates this defect, whereas replacement with an aromatic residue (Trp) allows the mutant enzyme to bypass the template analogs with relative ease.

The DNA-protein interaction modes of FEN-1 with gap substrates and their implication in preventing duplication mutations

Flap endonuclease-1 (FEN-1) is a structure-specific nuclease best known for its involvement in RNA primer removal and long-patch base excision repair. This enzyme is known to possess 5'-flap endo- (FEN) and 5'-3' exo- (EXO) nuclease activities. Recently, FEN-1 has been reported to also possess a gap endonuclease (GEN) activity, which is possibly involved in apoptotic DNA fragmentation and the resolution of stalled DNA replication forks. In the current study, we compare the kinetics of these activities to shed light on the aspects of DNA structure and FEN-1 DNA-binding elements that affect substrate cleavage. By using DNA binding deficient mutants of FEN-1, we determine that the GEN activity is analogous to FEN activity in that the single-stranded DNA region of DNA substrates interacts with the clamp region of FEN-1. In addition, we show that the C-terminal extension of human FEN-1 likely interacts with the downstream duplex portion of all substrates. Taken together, a substrate-binding model that explains how FEN-1, which has a single active center, can have seemingly different activities is proposed. Furthermore, based on the evidence that GEN activity in complex with WRN protein cleaves hairpin and internal loop substrates, we suggest that the GEN activity may prevent repeat expansions and duplication mutations.

Flap endonuclease 1: a central component of DNA metabolism

One strand of cellular DNA is generated as RNA-initiated discontinuous segments called Okazaki fragments that later are joined. The RNA terminated region is displaced into a 5' single-stranded flap, which is removed by the structure-specific flap endonuclease 1 (FEN1), leaving a nick for ligation. Similarly, in long-patch base excision repair, a damaged nucleotide is displaced into a flap and removed by FEN1. FEN1 is a genome stabilization factor that prevents flaps from equilibrating into structures that lead to duplications and deletions. As an endonuclease, FEN1 enters the flap from the 5' end and then tracks to cleave the flap base. Cleavage is oriented by the formation of a double flap. Analyses of FEN1 crystal structures suggest mechanisms for tracking and cleavage. Some flaps can form self-annealed and template bubble structures that interfere with FEN1. FEN1 interacts with other nucleases and helicases that allow it to act efficiently on structured flaps. Genetic and biochemical analyses continue to reveal many roles of FEN1.

Structure-specific DNA-induced conformational changes in Taq polymerase revealed by small angle neutron scattering

The DNA polymerase I from Thermus aquaticus (Taq polymerase) performs lagging-strand DNA synthesis and DNA repair. Taq polymerase contains a polymerase domain for synthesizing a new DNA strand and a 5'-nuclease domain for cleaving RNA primers or damaged DNA strands. The extended crystal structure of Taq polymerase poses a puzzle on how this enzyme coordinates its polymerase and the nuclease activities to generate only a nick. Using contrast variation solution small angle neutron scattering, we have examined the conformational changes that occur in Taq polymerase upon binding "overlap flap" DNA, a structure-specific DNA substrate that mimics the substrate in strand replacement reactions. In solution, apoTaq polymerase has an overall expanded equilibrium conformation similar to that in the crystal structure. Upon binding to the DNA substrate, both the polymerase and the nuclease domains adopt more compact overall conformations, but these changes are not enough to bring the two active sites close enough to generate a nick. Reconstruction of the three-dimensional molecular envelope from small angle neutron scattering data shows that in the DNA-bound form, the nuclease domain is lifted up relative to its position in the non-DNA-bound form so as to be in closer contact with the thumb and palm subdomains of the polymerase domain. The results suggest that a form of structure sensing is responsible for the coordination of the polymerase and nuclease activities in nick generation. However, interactions between the polymerase and the nuclease domains can assist in the transfer of the DNA substrate from one active site to the other.

Simulation and modeling of nucleic acid structure, dynamics and interactions

In moving towards the simulation of larger nucleic acid assemblies over longer timescales that include more accurate representations of the environment, we are nearing the end of an era characterized by single nanosecond molecular dynamics simulation of nucleic acids. We are excited by the promise and predictability of the modeling methods, yet remain prudently cautious of sampling and force field limitations. Highlights include the accurate representation of subtle drug-DNA interactions, the detailed study of modified and unusual nucleic acid structures, insight into the influence of dynamics on the structure of DNA, and exploration of the interaction of solvent and ions with nucleic acids.

Interaction interface of human flap endonuclease-1 with its DNA substrates

Flap endonuclease-1 or FEN-1 is a structure-specific and multifunctional nuclease critical for DNA replication, repair, and recombination; however, its interaction with DNA substrates has not been fully understood. In the current study, we have defined the borders of the interaction between the FEN-1 protein and its DNA substrates and identified six clusters of conserved positively charged amino acid residues, which are in direct contact with DNA substrate. To map further the corresponding interactions between FEN-1 residues and DNA substrates, we performed biochemical assays employing a series of flap DNA substrates lacking some structural components and a series of binding-deficient point mutants of FEN-1. It was revealed that Arg(47), Arg(70), and Lys(326)-Arg(327) of FEN-1 interact with the upstream duplex of DNA substrates, whereas Lys(244)-Arg(245) interact with the downstream duplex. This result indicates the orientation of the FEN-1-DNA interaction. Moreover, Arg(70) and Arg(47) were determined to interact with the sites around the 2nd nucleotide (Arg(70)) or the 5th/6th nucleotide (Arg(47)) of the template strand in the upstream duplex portion counting from the nick point of the flap substrate. Together with previously published data and the crystallographic ainformation from the FEN-1.DNA complex that we published recently (Chapados, B. R., Hosfield, D. J., Han, S., Qiu, J., Yelent, B., Shen, B., Tainer, J. A. (2004) Cell 116, 39-50) we are able to propose a reasonable model for how the human FEN-1 protein interacts with its DNA substrates.

Structural basis for FEN-1 substrate specificity and PCNA-mediated activation in DNA replication and repair

Flap EndoNuclease-1 (FEN-1) and the processivity factor proliferating cell nuclear antigen (PCNA) are central to DNA replication and repair. To clarify the molecular basis of FEN-1 specificity and PCNA activation, we report here structures of FEN-1:DNA and PCNA:FEN-1-peptide complexes, along with fluorescence resonance energy transfer (FRET) and mutational results. FEN-1 binds the unpaired 3' DNA end (3' flap), opens and kinks the DNA, and promotes conformational closing of a flexible helical clamp to facilitate 5' cleavage specificity. Ordering of unstructured C-terminal regions in FEN-1 and PCNA creates an intermolecular beta sheet interface that directly links adjacent PCNA and DNA binding regions of FEN-1 and suggests how PCNA stimulates FEN-1 activity. The DNA and protein conformational changes, composite complex structures, FRET, and mutational results support enzyme-PCNA alignments and a kinked DNA pivot point that appear suitable to coordinate rotary handoffs of kinked DNA intermediates among enzymes localized by the three PCNA binding sites.

DNA-thumb interactions and processivity of T7 DNA polymerase in comparison to yeast polymerase eta

The replicative polymerase of bacteriophage T7 is structurally and mechanistically well characterized. The crystal structure of T7 DNA polymerase or gene 5 protein complexed to its processivity factor, Escherichia coli thioredoxin, a primer-template, and a dideoxynucleotide reveals how this enzyme interacts with the 3'-end of the primer-template, but does not show how thioredoxin confers processivity to the polymerase. In the crystal structure highly conserved amino acids Asn(335) and Ser(338) of the thumb subdomain of T7 DNA polymerase are seen to interact with phosphates 7 and 8 of the DNA template strand. Results with a mutant T7 DNA polymerase in which aliphatic residues are substituted for these amino acids and experiments with different length and methylphosphonate-modified primer-templates demonstrate that these interactions are essential for processive synthesis and d(A.T)(n) tract bypass. Our data with methylphosphonate-modified DNA suggests that thioredoxin confers processivity to T7 DNA polymerase in part by causing an interaction with the phosphate backbone or minor groove of DNA. Residues Asn(335) and Ser(338) may also function with a nearby helix-loop-helix motif located at residues 339-372 to enclose the DNA during processive synthesis. Our results suggest that this structure must be held close to the DNA by ionic interactions to function. These interactions also allow for DNA sliding but physically block the passage of a 3T bulge in the template. In contrast, yeast polymerase eta, a polymerase that non-mutagenically repairs cis-syn thymidine dimers, allows the same bulge to slide past its thumb subdomain during synthesis. A relaxed thumb interaction with the DNA could account for the notably low processivity of polymerase eta.

Aromatic residues located close to the active center are essential for the catalytic reaction of flap endonuclease-1 from hyperthermophilic archaeon Pyrococcus horikoshii

Flap endonuclease-1 (FEN-1) possessing 5'-flap endonuclease and 5'-->3' exonuclease activity plays important roles in DNA replication and repair. In this study, the kinetic parameters of mutants at highly conserved aromatic residues, Tyr33, Phe35, Phe79, and Phe278-Phe279, in the vicinity of the catalytic centers of FEN-1 were examined. The substitution of these aromatic residues with alanine led to a large reduction in kcat values, although these mutants retained Km values similar to that of the wild-type enzyme. Notably, the kcat of Y33A and F79A decreased 333-fold and 71-fold, respectively, compared with that of the wild-type enzyme. The aromatic residues Tyr33 and Phe79, and the aromatic cluster Phe278-Phe279 mainly contributed to the recognition of the substrates without the 3' projection of the upstream strand (the nick, 5'-recess-end, single-flap, and pseudo-Y substrates) for the both exo- and endo-activities, but played minor roles in recognizing the substrates with the 3' projection (the double flap substrate and the nick substrate with the 3' projection). The replacement of Tyr33, Phe79, and Phe278-Phe279, with non-charged aromatic residues, but not with aliphatic hydrophobic residues, recovered the kcat values almost fully for the substrates without the 3' projection of the upstream strand, suggesting that the aromatic groups of Tyr33, Phe79, and Phe278-Phe279 might be involved in the catalytic reaction, probably via multiple stacking interactions with nucleotide bases. The stacking interactions of Tyr33 and Phe79 might play important roles in fixing the template strand and the downstream strand, respectively, in close proximity to the active center to achieve the productive transient state leading to the hydrolysis.