Wing 1 of protein HOP2 is as important as helix 3 in DNA binding by MD simulation

The repair of programmed DNA double-strand breaks through recombination is required for proper association and disjunction of the meiotic homologous chromosomes. Meiosis-specific protein HOP2 plays essential roles in recombination by promoting recombinase activities. The N-terminal domain of HOP2 interacts with DNA through helix 3 (H3) and wing 1 (W1). Mutations in wing 1 (Y65A/K67A/Q68A) slightly weakened the binding but mutations in helices 2 and 3 (Q30A/K44A/K49A) nearly abolished the binding. To better understand such differential effects at atomic level, molecular dynamics simulations were employed. Despite losing some hydrogen bonds, the W1-mutant DNA complex was rescued by stronger hydrophobic interactions. For the wild type and W1-mutant, the protein was found to slide along the DNA grooves as the DNA rolls along its double-helix axis. This motion could be functionally important to facilitate the precise positioning of the single-stranded DNA with the homologous double-stranded DNA. The sliding motion was reduced in the W1-mutant. The H-mutant nearly lost all intermolecular interactions. Moreover, an additional mutation in wing 1 (Y65A/K67A/Q68A/K69A) also caused complete complex dissociation. Therefore, both wing 1 and helix 3 make important contribution to the DNA binding, which could be important to the strand invasion function of HOP2 homodimer and HOP2-MND1 heterodimer. Similar to cocking a medieval crossbow with the archer’s foot placed in the stirrup, wing 1 may push the minor groove to cause distortion while helix 3 grabs the major groove.     

R248Q mutation—Beyond p53‐DNA binding

R248 in the DNA binding domain (DBD) of p53 interacts directly with the minor groove of DNA. Earlier nuclear magnetic resonance (NMR) studies indicated that the R248Q mutation resulted in conformation changes in parts of DBD far from the mutation site. However, how information propagates from the mutation site to the rest of the DBD is still not well understood. We performed a series of all‐atom molecular dynamics (MD) simulations to dissect sterics and charge effects of R248 on p53‐DBD conformation: (i) wild‐type p53 DBD; (ii) p53 DBD with an electrically neutral arginine side‐chain; (iii) p53 DBD with R248A; (iv) p53 DBD with R248W; and (v) p53 DBD with R248Q. Our results agree well with experimental observations of global conformational changes induced by the R248Q mutation. Our simulations suggest that both charge‐ and sterics are important in the dynamics of the loop (L3) where the mutation resides. We show that helix 2 (H2) dynamics is altered as a result of a change in the hydrogen bonding partner of D281. In turn, neighboring L1 dynamics is altered: in mutants, L1 predominantly adopts the recessed conformation and is unable to interact with the major groove of DNA. We focused our attention the R248Q mutant that is commonly found in a wide range of cancer and observed changes at the zinc‐binding pocket that might account for the dominant negative effects of R248Q. Furthermore, in our simulations, the S6/S7 turn was more frequently solvent exposed in R248Q, suggesting that there is a greater tendency of R248Q to partially unfold and possibly lead to an increased aggregation propensity. Finally, based on the observations made in our simulations, we propose strategies for the rescue of R248Q mutants. Proteins 2015; 83:2240–2250. © 2015 Wiley Periodicals, Inc.     

Molecular and structural insight into plasmodium falciparum RIO2 kinase

Among approximately 65 kinases of the malarial genome, RIO2 (right open reading frame) kinase belonging to the atypical class of kinase is unique because along with a kinase domain, it has a highly conserved N-terminal winged helix (wHTH) domain. The wHTH domain resembles the wing like domain found in DNA binding proteins and is situated near to the kinase domain. Ligand binding to this domain may reposition the kinase domain leading to inhibition of enzyme function and could be utilized as a novel allosteric site to design inhibitor. In the present study, we have generated a model of RIO2 kinase from Plasmodium falciparum utilizing multiple modeling, simulation approach. A novel putative DNA-binding site is identified for the first time in PfRIO2 kinase to understand the DNA binding events involving wHTH domain and flexible loop. Induced fit DNA docking followed by minimization, molecular dynamics simulation, energetic scoring and binding mode studies are used to reveal the structural basis of PfRIO2-ATP-DNA complex. Ser105 as a potential site of phosphorylation is revealed through the structural studies of ATP binding in PfRIO2. Overall the present study discloses the structural facets of unknown PfRIO2 complex and opens an avenue toward exploration of novel drug target.     

Characterization of DNA Binding Property of the HIV-1 Host Factor and Tumor Suppressor Protein Integrase Interactor 1 (INI1/hSNF5)

Integrase Interactor 1 (INI1/hSNF5) is a component of the hSWI/SNF chromatin remodeling complex. The INI1 gene is either deleted or mutated in rhabdoid cancers like ATRT (Atypical terratoid and rhabdoid tumor). INI1 is also a host factor for HIV-1 replication. INI1 binds DNA non-specifically. However, the mechanism of DNA binding and its biological role are unknown. From agarose gel retardation assay (AGRA), Ni-NTA pull-down and atomic force microscopy (AFM) studies we show that amino acids 105–183 of INI1 comprise the minimal DNA binding domain (DBD). The INI1 DBD is absent in plants and in yeast SNF5. It is present in Caenorhabditis elegans SNF5, Drosophila melanogaster homologue SNR1 and is a highly conserved domain in vertebrates. The DNA binding property of this domain in SNR1, that is only 58% identical to INI1/hSNF5, is conserved. Analytical ultracentrifugation studies of INI1 DBD and INI1 DBD:DNA complexes at different concentrations show that the DBD exists as a monomer at low protein concentration and two molecules of monomer binds one molecule of DNA. At high protein concentration, it exists as a dimer and binds two DNA molecules. Furthermore, isothermal calorimetry (ITC) experiments demonstrate that the DBD monomer binds DNA with a stoichiometry (N) of ∼0.5 and K_d  = 0.94 µM whereas the DBD dimer binds two DNA molecules sequentially with K’_d1 = 222 µM and K’_d2 = 1.16 µM. Monomeric DBD binding to DNA is enthalpy driven (ΔH = –29.9 KJ/mole). Dimeric DBD binding to DNA is sequential with the first binding event driven by positive entropy (ΔH’₁ = 115.7 KJ/mole, TΔS’₁ = 136.8 KJ/mole) and the second binding event driven by negative enthalpy (ΔH’₂ = –106.3 KJ/mole, TΔS’₂ = –75.7 KJ/mole). Our model for INI1 DBD binding to DNA provides new insights into the mechanism of DNA binding by INI1.     

Molecular Dynamics Simulation in Solvent of the Estrogen Receptor Protein DNA Binding Domain in Complex with a Non-Consensus Estrogen Response Element DNA Sequence

Abstract

We investigated protein/DNA interactions, using molecular dynamics simulations computed between a 10 Angstom water layer model of the estrogen receptor (ER) protein DNA binding domain (DBD) amino acids and DNA of a non-consensus estrogen response element (ERE) consisting of 29 nucleotide base pairs. This ERE nucleotide sequence occurs naturally upstream of the Xenopus laevis Vitelligenin AI gene. The ER DBD is encoded by three exons. Namely, exons 2 and 3 which encode the two zinc binding motifs and a sequence of exon 4 which encodes a predicted alpha helix. We generated a computer model of the ER DBD using atomic coordinates derived from the average of 30 nuclear magnetic resonance (NMR) spectroscopy coordinate sets. Amino acids on the carboxyl end of the ER DBD were disordered in both X-ray crystallography and NMR determinations and no coordinates were reported. This disordered region includes 10 amino acids of a predicted alpha helix encoded in exon 4 at the exon 3/4 splice junction. These amino acids are known to be important in DNA binding and are also believed to function as a nuclear translocation signal sequence for the ER protein. We generated a computer model of the predicted alpha helix consisting of the 10 amino acids encoded in exon 4 and attached this helix to the carboxyl end of the ER DBD at the exon 3/4 splice junction site. We docked the ER DBD model within the DNA major groove halfsites of the 29 base pair non-consensus ERE and flanking nucleotides. We constructed a solvated model with the ER DBD/ERE complex surrounded by a ten Angstrom water layer and conducted molecular dynamics simulations. Hydrogen bonding interactions were monitored. In addition, van der Waals and electrostatic interaction energies were calculated. Amino acids of the ER DBD DNA recognition helix formed both direct and water mediated hydrogen bonds at cognate codon-anticodon nucleotide base and backbone sites within the ERE DNA right major groove halfsite. Amino acids of the ER DBD exon 4 encoded predicted alpha helix formed direct and water mediated H-bonds with base and backbone sites of their cognate codon-anticodon nucleotides within the minor grooves flanking the ERE DNA major groove halfsites. These interactions together induced bending of the DNA into the protein.     

Improving dideoxynucleotide-triphosphate utilisation by the hyper-thermophilic DNA polymerase from the archaeon Pyrococcus furiosus

Polymerases from the Pol-I family which are able to efficiently use ddNTPs have demonstrated a much improved performance when used to sequence DNA. A number of mutations have been made to the gene coding for the Pol-II family DNA polymerase from the archaeon Pyrococcus furiosus with the aim of improving ddNTP utilisation. ‘Rational’ alterations to amino acids likely to be near the dNTP binding site (based on sequence homologies and structural information) did not yield the desired level of selectivity for ddNTPs. However, alteration at four positions (Q472, A486, L490 and Y497) gave rise to variants which incorporated ddNTPs better than the wild type, allowing sequencing reactions to be carried out at lowered ddNTP:dNTP ratios. Wild-type Pfu–Pol required a ddNTP:dNTP ratio of 30:1; values of 5:1 (Q472H), 1:3 (L490W), 1:5 (A486Y) and 5:1 (Y497A) were found with the four mutants; A486Y representing a 150-fold improvement over the wild type. A486, L490 and Y497 are on an α-helix that lines the dNTP binding groove, but the side chains of the three amino acids point away from this groove; Q472 is in a loop that connects this α-helix to a second long helix. None of the four amino acids can contact the dNTP directly. Therefore, the increased selectivity for ddNTPs is likely to arise from two factors: (i) small overall changes in conformation that subtly alter the nucleotide triphosphate binding site such that ddNTPs become favoured; (ii) interference with a conformational change that may be critical both for the polymerisation step and discrimination between different nucleotide triphosphates.     

Structural basis of replication origin recognition by the DnaA protein

Escherichia coli DnaA binds to 9 bp sequences (DnaA boxes) in the replication origin, oriC, to form a complex initiating chromosomal DNA replication. In the present study, we determined the crystal structure of its DNA-binding domain (domain IV) complexed with a DnaA box at 2.1 Å resolution. DnaA domain IV contains a helix–turn–helix motif for DNA binding. One helix and a loop of the helix– turn–helix motif are inserted into the major groove and 5 bp (3′ two-thirds of the DnaA box sequence) are recognized through base-specific hydrogen bonds and van der Waals contacts with the C5-methyl groups of thymines. In the minor groove, Arg399, located in the loop adjacent to the motif, recognizes three more base pairs (5′ one-third of the DnaA box sequence) by base-specific hydrogen bonds. DNA bending by ~28° was also observed in the complex. These base-specific interactions explain how DnaA exhibits higher affinity for the strong DnaA boxes (R1, R2 and R4) than the weak DnaA boxes (R3 and M) in the replication origin.     

Mechanism of high-mobility group protein B enhancement of progesterone receptor sequence-specific DNA binding

The DNA-binding domain (DBD) of progesterone receptor (PR) is bipartite containing a zinc module core that interacts with progesterone response elements (PRE), and a short flexible carboxyl terminal extension (CTE) that interacts with the minor groove flanking the PRE. The chromosomal high-mobility group B proteins (HMGB), defined as DNA architectural proteins capable of bending DNA, also function as auxiliary factors that increase the DNA-binding affinity of PR and other steroid receptors by mechanisms that are not well defined. Here we show that the CTE of PR contains a specific binding site for HMGB that is required for stimulation of PR-PRE binding, whereas the DNA architectural properties of HMGB are dispensable. Specific PRE DNA inhibited HMGB binding to the CTE, indicating that DNA and HMGB–CTE interactions are mutually exclusive. Exogenous CTE peptide increased PR-binding affinity for PRE as did deletion of the CTE. In a PR-binding site selection assay, A/T sequences flanking the PRE were enriched by HMGB, indicating that PR DNA-binding specificity is also altered by HMGB. We conclude that a transient HMGB–CTE interaction alters a repressive conformation of the flexible CTE enabling it to bind to preferred sequences flanking the PRE.     

Molecular dynamics simulations of the glucocorticoid receptor DNA-binding domain in complex with DNA and free in solution.

Molecular dynamics simulations have been performed on the glucocorticoid receptor DNA binding domain (GR DBD) in aqueous solution as a dimer in complex with DNA and as a free monomer. In the simulated complex, we find a slightly increased bending of the DNA helix axis compared with the crystal structure in the spacer region of DNA between the two half-sites that are recognized by GR DBD. The bend is mainly caused by an increased number of interactions between DNA and the N-terminal extended region of the sequence specifically bound monomer. The recognition helices of GR DBD are pulled further into the DNA major groove leading to a weakening of the intrahelical hydrogen bonds in the middle of the helices. Many ordered water molecules with long residence times are found at the intermolecular interfaces of the complex. The hydrogen-bonding networks (including water bridges) on either side of the DNA major groove involve residues that are highly conserved within the family of nuclear receptors. Very similar hydrogen-bonding networks are found in the estrogen receptor (ER) DBD in complex with DNA, which suggests that this is a common feature for proper positioning of the recognition helix in ER DBD and GR DBD.     

Triple helical DNA in a duplex context and base pair opening

It is fundamental to explore in atomic detail the behavior of DNA triple helices as a means to understand the role they might play in vivo and to better engineer their use in genetic technologies, such as antigene therapy. To this aim we have performed atomistic simulations of a purine-rich antiparallel triple helix stretch of 10 base triplets flanked by canonical Watson–Crick double helices. At the same time we have explored the thermodynamic behavior of a flipping Watson–Crick base pair in the context of the triple and double helix. The third strand can be accommodated in a B-like duplex conformation. Upon binding, the double helix changes shape, and becomes more rigid. The triple-helical region increases its major groove width mainly by oversliding in the negative direction. The resulting conformations are somewhere between the A and B conformations with base pairs remaining almost perpendicular to the helical axis. The neighboring duplex regions maintain a B DNA conformation. Base pair opening in the duplex regions is more probable than in the triplex and binding of the Hoogsteen strand does not influence base pair breathing in the neighboring duplex region.     

The DNA sequence specificity of HMG boxes lies in the minor wing of the structure.

To establish the basis of sequence-specific DNA recognition by HMG boxes we separately transferred the minor and major wings from the sequence-specific HMG box of TCF1 alpha into their equivalent position in the non-sequence-specific box 2 of HMG1. Thus chimera THT1 contains the minor wing (of 11 N-terminal and 25 C-terminal residues) from the HMG box of TCF1 alpha and the major wing (the 45 residue central section) from HMG1 box 2, whilst the situation is reversed in chimera HTH1. The structural integrity of the two chimeric proteins was established by CD, NMR and their binding to four-way junction DNA. Gel retardation and circular permutation assays showed that only chimera THT1, containing the TCF1 alpha minor wing, formed a sequence-specific complex and bent the DNA. The bend angle was estimated to be 59 degrees for chimera THT1 and 52 degrees for the HMG box of TCF1 alpha. Our results, in combination with mutagenesis and other data, suggests a model for the DNA binding of HMG boxes in which the N-terminal residues and part of helix 1 contact the minor groove on the outside of a bent DNA duplex.