Molecular dynamics simulations reveal that apo‐HisJ can sample a closed conformation

The Escherichia coli histidine binding protein HisJ is a type II periplasmic binding protein (PBP) that preferentially binds histidine and interacts with its cytoplasmic membrane ABC transporter, HisQMP₂, to initiate histidine transport. HisJ is a bilobal protein where the larger Domain 1 is connected to the smaller Domain 2 via two linking strands. Type II PBPs are thought to undergo “Venus flytrap” movements where the protein is able to reversibly transition from an open to a closed conformation. To explore the accessibility of the closed conformation to the apo state of the protein, we performed a set of all‐atom molecular dynamics simulations of HisJ starting from four different conformations: apo‐open, apo‐closed, apo‐semiopen, and holo‐closed. The simulations reveal that the closed conformation is less dynamic than the open one. HisJ experienced closing motions and explored semiopen conformations that reverted to closed forms resembling those found in the holo‐closed state. Essential dynamics analysis of the simulations identified domain closing/opening and twisting as main motions. The formation of specific inter‐hinge strand and interdomain polar interactions contributed to the adoption of the closed apo‐conformations although they are up to 2.5‐fold less prevalent compared with the holo‐closed simulations. The overall sampling of the closed form by apo‐HisJ provides a rationale for the binding of unliganded PBPs with their cytoplasmic membrane ABC transporters. Proteins 2014; 82:386–398. © 2013 Wiley Periodicals, Inc.     

Oligomeric states of bacteriophage T7 gene 4 primase/helicase

Electron microscopic and crystallographic data have shown that the gene 4 primase/helicase encoded by bacteriophage T7 can form both hexamers and heptamers. After cross-linking with glutaraldehyde to stabilize the oligomeric protein, hexamers and heptamers can be distinguished either by negative stain electron microscopy or electrophoretic analysis using polyacrylamide gels. We find that hexamers predominate in the presence of either dTTP or beta,gamma-methylene dTTP whereas the ratio between hexamers and heptamers is nearly the converse in the presence of dTDP. When formed, heptamers are unable to efficiently bind either single-stranded DNA or double-stranded DNA. We postulate that a switch between heptamer to hexamer may provide a ring-opening mechanism for the single-stranded DNA binding pathway. Accordingly, we observe that in the presence of both nucleoside di- and triphosphates the gene 4 protein exists as a hexamer when bound to single-stranded DNA and as a mixture of heptamer and hexamer when not bound to single-stranded DNA. Furthermore, altering regions of the gene 4 protein postulated to be conformational switches for dTTP-dependent helicase activity leads to modulation of the heptamer to hexamer ratio.     

Molecular cloning of a splicing variant of human RECQL helicase

A cDNA was isolated from the human heart library. This cDNA variant was produced by the deletion of 176 bases at the 5(') end of human RecQL gene, presumably by an alternative mRNA splicing. The cDNA contains two open reading frames and so may encode two isoforms of human RECQL. The first isoform is a 105 amino acid protein with the first 53 N-terminal amino acids identical to the known sequence of RECQL protein and followed by 52 amino acids introduced by in-frame premature stop codon. The second isoform is a 537 amino acid protein that has the same sequence as the published human RECQL helicase, except the first 112 amino acids at the N-terminal end were absent. The possible roles of both of these proteins are discussed.     

荧光偏振技术研究Bloom解旋酶催化核心与双链DNA的相互作用

Bloom 综合症(BLM)解旋酶是RecQ家族DNA解旋酶中的一个重要成员,参与了DNA复制、修复、转录、重组以及端粒的维持等细胞代谢过程,在维持染色体的稳定性中具有重要的作用．BLM解旋酶的突变可导致Bloom综合症,患者遗传不稳定易患多种类型癌症．本研究运用荧光偏振技术研究BLM解旋酶催化核心(BLM642-1290)与双链DNA(dsDNA)的相互作用,分析其相关特征参数,了解BLM642-1290解旋酶与dsDNA的结合和解链特性．结果表明,BLM642-1290解旋酶与dsDNA的结合和解链和dsDNA3’末端的单链DNA(ssDNA)长度有关;解旋酶优先结合于dsDNA底物的ssDNA末端,且每分子解旋酶可结合9.6 nt的ssDNA;dsDNA3’末端ssDNA的长度为9.6 nt时,解旋酶的解链效率达到最大且不再随其长度而变化．另外,BLM642-1290解旋酶也能够结合和解链钝末端dsDNA,但其结合亲和力和解链效率低于有3’末端ssDNA的dsDNA．推测BLM642-1290解旋酶在与dsDNA底物结合和解链时是单体形式,可能以尺蠖的形式解开dsDNA．这些结果可为进一步研究BLM解旋酶的功能特征提供理论基础．     

Partial purification and characterization of a DNA helicase from chloroplasts of Glycine max

A DNA helicase activity was detected in extracts of purified chloroplasts from the SB-1 cell line of Glycine max and partially purified by column chromatography on DEAE cellulose, phosphocellulose, and single-stranded DNA cellulose. The chloroplast helicase has a DNA-dependent ATPase activity, and its strand displacement activity is strictly dependent upon the presence of a nucleoside triphosphate and Mg²⁺ or Mn²⁺. Strand displacement activity does not require a free unannealed single-strand or replication fork-like structure.     

Bacillus subtilis bacteriophage SPP1 G40P helicase lacking the n-terminal domain unwinds DNA bidirectionally

Bacillus subtilis bacteriophage SPP1 G40P hexameric replicative DNA helicase unidirectionally translocates with a 5'-->3' polarity while separating the DNA strands. A G40P mutant derivative lacking the N-terminal domain (containing amino acid residues 110-442 from G40P, G40PDeltaN109) was purified and characterized. G40PDeltaN109 showed an ATPase activity that was dependent on the presence of single-stranded (ss) DNA. Unlike G40P, G40PDeltaN109 was shown to bind with similar affinity both ssDNA arms of forked structures by nuclease protection assays. In a pH-dependent manner, G40PDeltaN109 unwound a branched double-arm substrate preferentially with a 3'-->5' polarity. Our results show that the linker region and the C-terminal domain of G40P are sufficient to render an enzyme capable of encircling the ssDNA tails of the forked DNA and to unwind DNA with both 5'-->3' and 3'-->5' polarity. The presence of the N-terminal domain, which does not play an essential role in helicase action, might be required indirectly for strand discrimination and polarity of translocation.     

Coupling translocation with nucleic acid unwinding by NS3 helicase

We present a semiquantitative model for translocation and unwinding activities of monomeric nonstructural protein 3 (NS3) helicase. The model is based on structural, biochemical, and single-molecule measurements. The model predicts that the NS3 helicase actively unwinds duplex by reducing more than 50% the free energy that stabilizes base pairing/stacking. The unwinding activity slows the movement of the helicase in a sequence-dependent manner, lowering the average unwinding efficiency to less than 1 bp per ATP cycle. When bound with ATP, the NS3 helicase can display significant translocational diffusion. This increases displacement fluctuations of the helicase, decreases the average unwinding efficiency, and enhances the sequence dependence. Also, interactions between the helicase and the duplex stabilize the helicase at the junction, facilitating the helicase's unwinding activity while preventing it from dissociating. In the presence of translocational diffusion during active unwinding, the dissociation rate of the helicase also exhibits sequence dependence. Based on unwinding velocity fluctuations measured from single-molecule experiments, we estimate the diffusion rate to be on the order of 10 s^− 1 . The generic features of coupling single-stranded nucleic acid translocation with duplex unwinding presented in this work may apply generally to a class of helicases.     

Molecular dynamics simulations of xDNA

xDNA is a modified DNA, which contains natural as well as expanded bases. Expanded bases are generated by the addition of a benzene spacer to the natural bases. A set of AMBER force‐field parameters were derived for the expanded bases and the structural dynamics of the xDNA decamer ( xT5 ′ G xT A xC xG C xA xG T3′ ) · ( xA5′ C T xG C G xT A xC A3′) was explored using a 22 ns molecular dynamics simulation in explicit solvent. During the simulation, the duplex retained its Watson‐Crick base‐pairing and double helical structure, with deviations from the starting B‐form geometry towards A‐form; the deviations are mainly in the backbone torsion angles and in the helical parameters. The sugar pucker of the residues were distributed among a variety of modes; C2′ endo, C1′ exo, O4′ endo, C4′ exo, C2′ exo, and C3′ endo. The enhanced stacking interactions on account of the modification in the bases could help to retain the duplex nature of the helix with minor deviations from the ideal geometry. In our simulation, the xDNA showed a reduced minor groove width and an enlarged major groove width in comparison with the NMR structure. Both the grooves are larger than that of standard B‐DNA, but major groove width is larger than that of A‐DNA with almost equal minor groove width. The enlarged groove widths and the possibility of additional hydration in the grooves makes xDNA a potential molecule for various applications. © 2009 Wiley Periodicals, Inc. Biopolymers 91: 351–360, 2009. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Identification of the PcrA DNA helicase reaction pathway by applying advanced targeted molecular dynamic simulations

PcrA DNA helicase uses the free energy of hydrolysis and binding of ATP to unwind double-stranded DNA (ds-DNA). There are two states of PcrA, termed the substrate and product complexes and, through the conformational changes between these two states, PcrA moves along ds-DNA and separates the two strands. In this study, two different methods, namely chain minimisation (CM, less reliable method) and auto targeted molecular dynamic (TMD) simulation (more reliable), were performed to generate two different initial reaction pathways between these two states, and then fixed root mean square distance (RMSD) TMD simulation was performed to optimise these two initial pathways. In general, the two optimised pathways share very similar major conformational changes, but are different in the minor motions. The potential energy profiles of the two improved pathways are generally similar, but the one generated by the improved TMD path is slightly lower. Considering the poor reliability of the initial path generated by CM and insignificant improvements of the auto-TMD path, our study suggests that fixed RMSD TMD simulation can generate reliable reaction pathways, but the different initial paths still have some influence on the detailed conformational analysis.     

Optimization of a novel potent and selective bacterial DNA helicase inhibitor scaffold from a high throughput screening hit

Benzobisthiazole derivatives were identified as novel helicase inhibitors through high throughput screening against purified Staphylococcus aureus (Sa) and Bacillus anthracis (Ba) replicative helicases. Chemical optimization has produced compound 59 with nanomolar potency against the DNA duplex strand unwinding activities of both B. anthracis and S. aureus helicases. Selectivity index (SI = CC₅₀/IC₅₀) values for 59 were greater than 500. Kinetic studies demonstrated that the benzobisthiazole-based bacterial helicase inhibitors act competitively with the DNA substrate. Therefore, benzobisthiazole helicase inhibitors represent a promising new scaffold for evaluation as antibacterial agents.     

DbpA is a region‐specific RNA helicase

DbpA is a DEAD‐box RNA helicase implicated in RNA structural rearrangements in the peptidyl transferase center. DbpA contains an RNA binding domain, responsible for tight binding of DbpA to hairpin 92 of 23S ribosomal RNA, and a RecA‐like catalytic core responsible for double‐helix unwinding. It is not known if DbpA unwinds only the RNA helices that are part of a specific RNA structure, or if DbpA unwinds any RNA helices within the catalytic core's grasp. In other words, it is not known if DbpA is a site‐specific enzyme or region‐specific enzyme. In this study, we used protein and RNA engineering to investigate if DbpA is a region‐specific or a site‐specific enzyme. Our data suggest that DbpA is a region‐specific enzyme. This conclusion has an important implication for the physiological role of DbpA. It suggests that during ribosome assembly, DbpA could bind with its C‐terminal RNA binding domain to hairpin 92, while its catalytic core may unwind any double‐helices in its vicinity. The only requirement for a double‐helix to serve as a DbpA substrate is for the double‐helix to be positioned within the catalytic core's grasp.     

Mutant bovine odorant-binding protein: Temperature affects the protein stability and dynamics as revealed by infrared spectroscopy and molecular dynamics simulations

Bovine odorant-binding protein (bOBP), a member of the lipocalin family, presents the so-called 3D "domain-swapped" protein structure. In fact, in solution, it appears as a dimer in which each monomer is composed by the classical lipocalin fold, with a central beta-barrel followed by a stretch of residues and the alpha-helix domain protruding out of the barrel and crossing the dimer interface. Recently, a deswapped mutant form of bOBP was obtained, in which a Gly residue was inserted after position 121 and the two residues in position 64 and 156 were replaced by Cys residues for restoring the disulfide bridge common to the lipocalin family. In this work, we used Fourier transform infrared spectroscopy and molecular dynamics simulations to investigate the effect of temperature on the structural stability and conformational dynamics of the mutant bOBP. The spectroscopic and molecular simulation data pointed out that the hydrophobic regions of the protein matrix appear to be an important factor for the protein stability and integrity. In addition, it was also found that the mutant bOBP is significantly stabilized by the binding of the ligand, which may have an impact on the biological function of bOBP. The obtained results will allow for a better use of this protein as probe for the design of advanced protein-based biosensors for the detection of compounds used in the fabrication of explosive powders.     

Molecular dynamics simulations of HPr under hydrostatic pressure

The histidine-containing protein (HPr) plays an important role in the phosphotransferase system (PTS). The deformations induced on the protein structure at high hydrostatic pressure values (4, 50, 100, 150, and 200 MPa) were previously (H. Kalbitzer, A. G?rler, H. Li, P. Dubovskii, A. Hengstenberg, C. Kowolik, H. Yamada, and K. Akasaka, Protein Science 2000, Vol. 9, pp. 693-703) analyzed by NMR experiments: the nonlinear variations of the amide chemical shifts at high pressure values were supposed to arise from induced shifts in the protein conformational equilibrium. Molecular dynamics (MD) simulations are here performed, to analyze the protein internal mobility at 0.1 MPa, and to relate the nonlinear variations of chemical shifts observed at high pressure, to variations in conformational equilibrium. The global features of the protein structure are only slightly modified along the pressure. Nevertheless, the values of the Voronoi residues volumes show that the residues of alpha-helices are more compressed that those belonging to the beta-sheet. The alpha-helices are also displaying the largest internal mobility and deformation in the simulations. The nonlinearity of the 1H chemical shifts, computed from the MD simulation snapshots, is in qualitative agreement with the nonlinearity of the experimentally observed chemical shifts.     

Double-stranded DNA-induced localized unfolding of HCV NS3 helicase subdomain 2

The NS3 helicase of the hepatitis C virus (HCV) unwinds double-stranded (ds) nucleic acid (NA) in an NTP-dependent fashion. Mechanistic details of this process are, however, largely unknown for the HCV helicase. We have studied the binding of dsDNA to an engineered version of subdomain 2 of the HCV helicase (d(2Delta)NS3h) by NMR and circular dichroism. Binding of dsDNA to d(2Delta)NS3h induces a local unfolding of helix (alpha(3)), which includes residues of conserved helicase motif VI (Q(460)RxxRxxR(467)), and strands (beta(1) and beta(8)) from the central beta-sheet. This also occurs upon lowering the pH (4.4) and introducing an R461A point mutation, which disrupt salt bridges with Asp 412 and Asp 427 in the protein structure. NMR studies on d(2Delta)NS3h in the partially unfolded state at low pH map the dsDNA binding site to residues previously shown to be involved in single-stranded DNA binding. Sequence alignment and structural comparison suggest that these Arg-Asp interactions are highly conserved in SF2 DEx(D/H) proteins. Thus, modulation of these interactions by dsNA may allow SF2 helicases to switch between conformations required for helicase function.     

T7 DNA helicase: a molecular motor that processively and unidirectionally translocates along single-stranded DNA

DNA helicases are molecular motors that use the energy from NTP hydrolysis to drive the process of duplex DNA strand separation. Here, we measure the translocation and energy coupling efficiency of a replicative DNA helicase from bacteriophage T7 that is a member of a class of helicases that assembles into ring-shaped hexamers. Presteady state kinetics of DNA-stimulated dTTP hydrolysis activity of T7 helicase were measured using a real time assay as a function of ssDNA length, which provided evidence for unidirectional translocation of T7 helicase along ssDNA. Global fitting of the kinetic data provided an average translocation rate of 132 bases per second per hexamer at 18 degrees C. While translocating along ssDNA, T7 helicase hydrolyzes dTTP at a rate of 49 dTTP per second per hexamer, which indicates that the energy from hydrolysis of one dTTP drives unidirectional movement of T7 helicase along two to three bases of ssDNA. One of the features that distinguishes this ring helicase is its processivity, which was determined to be 0.99996, which indicated that T7 helicase travels on an average about 75kb of ssDNA before dissociating. We propose that the ability of T7 helicase to translocate unidirectionally along ssDNA in an efficient manner plays a crucial role in DNA unwinding.