Conformational Dynamics of Thermus aquaticus DNA Polymerase I during Catalysis

Despite the fact that DNA polymerases have been investigated for many years and are commonly used as tools in a number of molecular biology assays, many details of the kinetic mechanism they use to catalyze DNA synthesis remain unclear. Structural and kinetic studies have characterized a rapid, pre-catalytic open-to-close conformational change of the Finger domain during nucleotide binding for many DNA polymerases including Thermus aquaticus DNA polymerase I (Taq Pol), a thermostable enzyme commonly used for DNA amplification in PCR. However, little has been performed to characterize the motions of other structural domains of Taq Pol or any other DNA polymerase during catalysis. Here, we used stopped-flow Förster resonance energy transfer to investigate the conformational dynamics of all five structural domains of the full-length Taq Pol relative to the DNA substrate during nucleotide binding and incorporation. Our study provides evidence for a rapid conformational change step induced by dNTP binding and a subsequent global conformational transition involving all domains of Taq Pol during catalysis. Additionally, our study shows that the rate of the global transition was greatly increased with the truncated form of Taq Pol lacking the N-terminal domain. Finally, we utilized a mutant of Taq Pol containing a de novo disulfide bond to demonstrate that limiting protein conformational flexibility greatly reduced the polymerization activity of Taq Pol.     

Dynamics of Loading the Escherichia coli DNA Polymerase Processivity Clamp

ABSTRACT

Sliding clamps and clamp loaders are processivity factors required for efficient DNA replication. Sliding clamps are ring-shaped complexes that tether DNA polymerases to DNA to increase the processivity of synthesis. Clamp loaders assemble these ring-shaped clamps onto DNA in an ATP-dependent reaction. The overall process of clamp loading is dynamic in that protein–protein and protein–DNA interactions must actively change in a coordinated fashion to complete the mechanical clamp-loading reaction cycle. The clamp loader must initially have a high affinity for both the clamp and DNA to bring these macromolecules together, but then must release the clamp on DNA for synthesis to begin. Evidence is presented for a mechanism in which the clamp-loading reaction comprises a series of binding reactions to ATP, the clamp, DNA, and ADP, each of which promotes some change in the conformation of the clamp loader that alters interactions with the next component of the pathway. These changes in interactions must be rapid enough to allow the clamp loader to keep pace with replication fork movement. This review focuses on the measurement of dynamic and transient interactions required to assemble the Escherichia coli sliding clamp on DNA.     

Molecular Dynamics of Channelrhodopsin at the Early Stages of Channel Opening

Channelrhodopsin (ChR) is a light-gated cation channel that responds to blue light. Since ChR can be readily expressed in specific neurons to precisely control their activities by light, it has become a powerful tool in neuroscience. Although the recently solved crystal structure of a chimeric ChR, C1C2, provided the structural basis for ChR, our understanding of the molecular mechanism of ChR still remains limited. Here we performed electrophysiological analyses and all-atom molecular dynamics (MD) simulations, to investigate the importance of the intracellular and central constrictions of the ion conducting pore observed in the crystal structure of C1C2. Our electrophysiological analysis revealed that two glutamate residues, Glu122 and Glu129, in the intracellular and central constrictions, respectively, should be deprotonated in the photocycle. The simulation results suggested that the deprotonation of Glu129 in the central constriction leads to ion leakage in the ground state, and implied that the protonation of Glu129 is important for preventing ion leakage in the ground state. Moreover, we modeled the 13-cis retinal bound; i.e., activated C1C2, and performed MD simulations to investigate the conformational changes in the early stage of the photocycle. Our simulations suggested that retinal photoisomerization induces the conformational change toward channel opening, including the movements of TM6, TM7 and TM2. These insights into the dynamics of the ground states and the early photocycle stages enhance our understanding of the channel function of ChR.     

Global Conformational Dynamics of a Y-Family DNA Polymerase during Catalysis

Replicative DNA polymerases are stalled by damaged DNA while the newly discovered Y-family DNA polymerases are recruited to rescue these stalled replication forks, thereby enhancing cell survival. The Y-family DNA polymerases, characterized by low fidelity and processivity, are able to bypass different classes of DNA lesions. A variety of kinetic and structural studies have established a minimal reaction pathway common to all DNA polymerases, although the conformational intermediates are not well defined. Furthermore, the identification of the rate-limiting step of nucleotide incorporation catalyzed by any DNA polymerase has been a matter of long debate. By monitoring time-dependent fluorescence resonance energy transfer (FRET) signal changes at multiple sites in each domain and DNA during catalysis, we present here a real-time picture of the global conformational transitions of a model Y-family enzyme: DNA polymerase IV (Dpo4) from Sulfolobus solfataricus. Our results provide evidence for a hypothetical DNA translocation event followed by a rapid protein conformational change prior to catalysis and a subsequent slow, post-chemistry protein conformational change. Surprisingly, the DNA translocation step was induced by the binding of a correct nucleotide. Moreover, we have determined the directions, rates, and activation energy barriers of the protein conformational transitions, which indicated that the four domains of Dpo4 moved in a synchronized manner. These results showed conclusively that a pre-chemistry conformational change associated with domain movements was too fast to be the rate-limiting step. Rather, the rearrangement of active site residues limited the rate of correct nucleotide incorporation. Collectively, the conformational dynamics of Dpo4 offer insights into how the inter-domain movements are related to enzymatic function and their concerted interactions with other proteins at the replication fork.     

Mechanism of Replicative DNA Polymerase Delta Pausing and a Potential Role for DNA Polymerase Kappa in Common Fragile Site Replication

Common fragile sites (CFSs) are hot spots of chromosomal breakage, and CFS breakage models involve perturbations of DNA replication. Here, we analyzed the contribution of specific repetitive DNA sequence elements within CFSs to the inhibition of DNA synthesis by replicative and specialized DNA polymerases (Pols). The efficiency of in vitro DNA synthesis was quantitated using templates corresponding to regions within FRA16D and FRA3B harboring AT-rich microsatellite and quasi-palindrome (QP) sequences. QPs were predicted to form stems of ~ 75–100% self-homology, separated by 3–9 bases of intervening sequences. Analysis of DNA synthesis progression by human Pol δ demonstrated significant synthesis perturbation both at [A]_n and [TA]_n repeats in a length-dependent manner and at short (< 40 base pairs) QP sequences. DNA synthesis by the Y-family polymerase κ was significantly more efficient than Pol δ through both types of repetitive elements. Using DNA trap experiments, we show that Pol δ pauses within CFS sequences are sites of enzyme dissociation, and dissociation was observed in the presence of RFC-loaded PCNA. We propose that enrichment of microsatellite and QP elements at CFS regions contributes to fragility by perturbing replication through multiple mechanisms, including replicative Pol pausing and dissociation. Our finding that Pol δ dissociates at specific CFS sequences is significant, since dissociation of the replication machinery and inability to efficiently recover the replication fork can lead to fork collapse and/or formation of double-strand breaks in vivo. Our biochemical studies also extend the potential involvement of Y-family polymerases in CFS maintenance to include polymerase κ.     

Single-Stranded DNA within Nanopores: Conformational Dynamics and Implications for Sequencing; a Molecular Dynamics Simulation Study

Engineered protein nanopores, such as those based on α-hemolysin from Staphylococcus aureus have shown great promise as components of next-generation DNA sequencing devices. However, before such protein nanopores can be used to their full potential, the conformational dynamics and translocation pathway of the DNA within them must be characterized at the individual molecule level. Here, we employ atomistic molecular dynamics simulations of single-stranded DNA movement through a model α-hemolysin pore under an applied electric field. The simulations enable characterization of the conformations adopted by single-stranded DNA, and allow exploration of how the conformations may impact on translocation within the wild-type model pore and a number of mutants. Our results show that specific interactions between the protein nanopore and the DNA can have a significant impact on the DNA conformation often leading to localized coiling, which in turn, can alter the order in which the DNA bases exit the nanopore. Thus, our simulations show that strategies to control the conformation of DNA within a protein nanopore would be a distinct advantage for the purposes of DNA sequencing.     

The Role of Molecular Dynamics Simulations for the Study of Slow Dynamics

Abstract

A two step strategy is proposed to study dynamical properties of a physical system much slower than the time scales accessible by molecular dynamics simulations. The strategy is applied to investigate the slow dynamics of supercooled liquids.     

Mechanism of Suppression of Chromosomal Instability by DNA Polymerase POLQ

Although a defect in the DNA polymerase POLQ leads to ionizing radiation sensitivity in mammalian cells, the relevant enzymatic pathway has not been identified. Here we define the specific mechanism by which POLQ restricts harmful DNA instability. Our experiments show that Polq-null murine cells are selectively hypersensitive to DNA strand breaking agents, and that damage resistance requires the DNA polymerase activity of POLQ. Using a DNA break end joining assay in cells, we monitored repair of DNA ends with long 3′ single-stranded overhangs. End joining events retaining much of the overhang were dependent on POLQ, and independent of Ku70. To analyze the repair function in more detail, we examined immunoglobulin class switch joining between DNA segments in antibody genes. POLQ participates in end joining of a DNA break during immunoglobulin class-switching, producing insertions of base pairs at the joins with homology to IgH switch-region sequences. Biochemical experiments with purified human POLQ protein revealed the mechanism generating the insertions during DNA end joining, relying on the unique ability of POLQ to extend DNA from minimally paired primers. DNA breaks at the IgH locus can sometimes join with breaks in Myc, creating a chromosome translocation. We found a marked increase in Myc/IgH translocations in Polq-defective mice, showing that POLQ suppresses genomic instability and genome rearrangements originating at DNA double-strand breaks. This work clearly defines a role and mechanism for mammalian POLQ in an alternative end joining pathway that suppresses the formation of chromosomal translocations. Our findings depart from the prevailing view that alternative end joining processes are generically translocation-prone.     

大肠杆菌DNA聚合酶Ⅰ的“碱基选择”作用

用双脱氧核苷酸研究了大肠杆菌ＤＮＡ聚合酶Ⅰ的Ｋｌｅｎｏｗ片段的“碱基选择”作用．当引物末端为３’ＯＨ时,Ｍｎ２＋代替Ｍｇ２＋对ｄＮＴＰ的Ｋｍ影响不大,而使ｄｄＮＴＰ的Ｋｉ大大降低,非配对的ｄｄＮＴＰ仍无抑制作用。当引物末端为ｄｄＮＭＰ且与模板配对时是ＤＮＡ合成的抑制剂,与模板配对的下一个ｄＮＴＰ可增强这一抑制作用,Ｍｎ２＋代替Ｍｇ２＋使ｄＮＴＰ的诱导作用大大增强而ｄｄＮＴＰ无诱导作用。这些结果表明引物末端对于酶选择核苷酸起着重要作用,ＤＮＡ聚合酶Ⅰ选择核苷酸的过程应是个有序的构象变化的过程。下一互补核苷酸的结合抑制３＇→５＇外切活性对引物末端ｄｄＮＭＰ的水解作用为此提供了又一证据。     

Taq DNA聚合酶功能区域的定位

通过参Ｕ法定点突变产生了ＴａｑＤＮＡ聚合酶Ｎ端分别缺失３个,２３５个,２８７个和４４３个氨基酸的４个缺失体,利用Ｂａｌ－３１连续缺失法产生了ＴａｑＤＮＡ聚合酶的Ｃ端分别缺失了２个、１６个、２９个、３２个、３４个氨基酸的５个缺失体．经ＤＮＡ聚合酶活性测定表明Ｎ端缺失３个,２３５个,２８７个氨基酸后活力和完整的Ｔａｑ相近,而缺失４４３个氨基酸后则失去了ＤＮＡ聚合酶活力;Ｃ端的５个缺失体都失去了ＤＮＡ聚合酶活性．据此ＴａｑＤＮＡ聚合酶的功能区域被定位在２８７～８３２氨基酸之间．     

On the Use of Molecular Dynamics Simulations for Probing Allostery through DNA

A recent study described an allosteric effect in which the binding of a protein to DNA is influenced by another protein bound nearby. The effect shows a periodicity of ∼10 basepairs and decays with increasing protein-protein distance. As a mechanistic explanation, the authors reported a similar periodic, decaying pattern of the correlation coefficient between major groove widths inferred from a shorter molecular dynamics simulation. Here we show that in a state-of-the-art, microsecond-long simulation of the same DNA sequence, the periodicity of the correlation coefficient is not observed. To study the problem further, we extend an earlier mechanical model of DNA allostery based on constrained minimization of effective quadratic deformation energy of the DNA. We demonstrate that, if the constraints mimicking the bound proteins are properly applied, the periodicity in the binding energy is indeed recovered.     

A Combined Quantum/Classical Molecular Dynamics Study of the Catalytic Mechanism of HIV Protease

Based on available three-dimensional structures of enzyme-inhibitor complexes, the mechanism of the reaction catalysed by HIV protease is studied using molecular dynamics simulations with molecular mechanics and combined quantum-mechanics/molecular-mechanics potential energy functions. The results support the general acid/general base catalysis mechanism, with Asp25′ protonated in the enzyme-substrate complex. In the enzyme-substrate complex, the lytic water molecule binds at a position different from the positions of the hydroxyl groups in various aspartic protease-inhibitor complexes. The carboxyl groups at the active site also adopt a different orientation. However, when the lytic water molecule approaches the scissile peptide, the reaction centre changes gradually to a conformation close to that derived from X-ray diffraction studies of various enzyme-inhibitor complexes. The proton transfer processes can take place only after the lytic water molecule has approached the scissile peptide bond to a certain degree. Qualitatively, the free-energy barrier associated with the nucleophilic attack step, which takes place at physiological pH, is comparable with the acid or base-catalysed reactions of model systems. The structure of the tetrahedral intermediate resulting from the nucleophilic attack step also indicates a straightforward pathway of the next reaction step, i.e. the breaking of the C-N bond.     

A Molecular Dynamics Simulation of a Polyamine-Induced Conformational Change of DNA. A Possible Mechanism for the B to Z Transition

Abstract

A 75ps molecular dynamics simulation has been performed on a fully solvated complex of spermine with the B DNA decamer (dGdC)₅ · (dGdC)₅. The simulation indicates a possible mechanism by which polyamines might induce the formation of a left-handed helix, the B to Z transition. Spermine was initially located in the major groove, hydrogen bonded to the helix. During the simulation the ligand migrates deeper into the DNA, maintaining strong hydrogen bonding to the central guanine bases and destroying the Watson-Crick base pairing with their respective cytosines. Significant rotation of these and other cytosine bases was observed, in part due to interactions of the helix with the aminopropyl chains of spermine. An intermediate B_II conformation might be of importance in this process.     

Molecular Dynamics Simulation Study of DNA Triplex Formed by Mixed Sequences in Solution

Abstract

The unrestrained molecular dynamics simulation of the triple helical DNA with mix sequences d(GACTGGTGAC)?d(CTGACCACTG)*d (GACTGGTGAC), using the particle mesh Ewald sum, is presented here. The Ewald summation method effectively eliminates the usual “cut-off” of the long—range interactions and allowed us to evaluate the full effect of the electrostatic forces. The AMBER5.0 force field has been used during the simulation in solvent. The MD results support a dynamically stable model of DNA triplex over the entire length of the trajectory. The duplex structure assumes the conformation, which is very close to B-DNA. In mixed sequences the purine bases occurs in both strand of DNA duplex. The bases of third strand do not favor the Hoogsteen or/and reverse Hoogsteen type of Hydrogen bonding but they form hydrogen bonds with the bases of both the strand of DNA duplex. The orientation of the third strand is parallel to one of the strand of duplex and all nucleotides (C, A, G & T) show isomorphic behavior with respect to the DNA duplex. The conformation of all the three strands is almost same except few exceptions. Due to interaction of third strand the conformational change in the duplex structure and a finite amount of displacement in the W-C base pairs have been observed. The conformational variation of the back bone torsion angles and helicoidal parameters, groove widths have been discussed. The sequence—dependent effects on local conformation, helicoidal and morphological structure, width of the grooves of DNA helix may have important implication for understanding the functional energetics and specificity of interactions of DNA and its triplexes with proteins, pharmaceutical agents and other legends.     

持续性催化酶类机理研究及其分子动力学模拟

酶分子催化机理研究是生命科学研究领域一个重要的问题.近80年来,过渡态理论在解释酶催化机理问题上占据了主导地位,结合热力学循环、锁钥学说、诱导契合学说以及酶活性中心柔性学说等理论,可以很好地解释多种酶分子的催化过程.近年来,随着蛋白质结构解析方法、单分子分析检测技术及计算机模拟技术的发展,人们对酶分子催化机理的认识愈加深刻.但持续性催化酶类的催化动力研究表明,过渡态理论的解释并不充分.本文对酶催化机理研究的相关进展进行了综述,并针对持续性酶类催化动态过程的特点提出了可能的研究方向及可行的研究方法.