Hydration at the TD Damaged Site of DNA and its Role in the Formation of Complex with T4 Endonuclease V

Abstract

An analysis of the distribution of water around DNA surface focusing on the role of the distribution of water molecules in the proper recognition of damaged site by repair enzyme T4 Endonuclease V was performed. The native DNA dodecamer, dodecamer with the thymine dimer (TD) and complex of DNA and part of repair enzyme T4 Endonuclease V were examined throughout the 500 ps of molecular dynamics simulation. During simulation the number of water molecules close to the DNA atoms and the residence time were calculated. There is an increase in number of water molecules lying in the close vicinity to TD if compared with those lying close to two native thymines (TT). Densely populated area with water molecules around TD is one of the factors detected by enzyme during scanning process. The residence time was found higher for molecule of the complex and the six water molecules were found occupying the stabile positions between the TD and catalytic center close to atoms P, C3′ and N3. These molecules originate water mediated hydrogen bond network that contribute to the stability of complex required for the onset of repair process.     

Homology Modeling and Molecular Dynamics Simulations of PBCV-1 Glycosylase Complexed with UV-damaged DNA

The UV-light damage specific DNA glycosylase from Chlorella virus strain PBCV-1 (pyrimidine dimer glycosylase; PDG) incises DNA at sites containing UV-induced thymidine dimers by catalyzing the breakage of the N-C-1 glycosyl bond. As the amino acid sequence of PDG is 41% identical to that of T4 endonuclease V (Endo V), and potential key active site residues are conserved, we used coordinates from a crystal structure of an Endo V complexed with DNA containing a cis-syn cyclobutane thymidine-dimer as a template to model a similar complex of PDG. Quantum mechanical calculations of the damaged base pair and the distance geometry based program DIAMOD were used to generate a PDG/DNA model whose backbone root mean square deviation (RMSD) to the Endo V/DNA structure was 0.5 Å, 0.5 Å, and 0.8 Å for DNA, protein, and the whole complex, respectively. To better understand structural details that could account for differences in activity of the two enzymes, molecular dynamics simulations were used to follow protein-DNA interactions in an aqueous environment. The simulations of the Endo V/DNA complex indicate new roles for Arg22 and Arg26 in the active site in recognizing irregular pairing and maintaining the strand separation needed for incision of the damaged bases. The model for the PDG/DNA complex and simulations thereof indicate a similar mechanism for DNA binding by this enzyme despite significant differences in residues maintaining the flipped-out adenine and strand separation in the area of damage. According to our model, PDGs increased affinity for substrate is probably due to a higher surface charge. Further, reduced packing density in the active site could account for PDGs activity on trans-syn II cyclobutane dimers.Electronic Supplementary Material available.     

Endonuclease V of Escherichia coli.

A small endodeoxyribonuclease )2.3 S) that is active on single-stranded DNA has been extensively purified from Escherichia coli so as to be free of other known DNases. It has an alkaline pH optimum (9.5), requires Mg2+, and makes 3'-hydroxy and 5'-phosphate termini. The nuclease nicks duplex DNA, particularly if treated with OsO4, irradiated with ultraviolet light, or exposed to pH 5. The uracil-containing duplex DNA from the Bacillus subtilis phage PBS-2 is an especially good substrate; it is made acid-soluble by levels of the enzyme which fail to produce any acid-soluble material in other single-stranded or duplex DNAs. Neither RNA nor RNA-DNA hybrid are degraded by the enzyme. The enzyme specificity suggests that it might act at abnormal regions in DNA, so that its in vivo function could be to initiate an excision repair sequence. Its high activity on uracil-containing DNA could imply that the enzyme provides an alternative mechanism for excising uracil residues from DNA to the pathway utilizing uracil-DNA N-glycosidase. We suggest that this enzyme be designated as endonuclease V of E. coli.     

Structure of Escherichia coli AlkA in Complex with Undamaged DNA

Because DNA damage is so rare, DNA glycosylases interact for the most part with undamaged DNA. Whereas the structural basis for recognition of DNA lesions by glycosylases has been studied extensively, less is known about the nature of the interaction between these proteins and undamaged DNA. Here we report the crystal structures of the DNA glycosylase AlkA in complex with undamaged DNA. The structures revealed a recognition mode in which the DNA is nearly straight, with no amino acid side chains inserted into the duplex, and the target base pair is fully intrahelical. A comparison of the present structures with that of AlkA recognizing an extrahelical lesion revealed conformational changes in both the DNA and protein as the glycosylase transitions from the interrogation of undamaged DNA to catalysis of nucleobase excision. Modeling studies with the cytotoxic lesion 3-methyladenine and accompanying biochemical experiments suggested that AlkA actively interrogates the minor groove of the DNA while probing for the presence of lesions.     

Mechanisms for defining supercoiling set point of DNA gyrase orthologs: I. A nonconserved acidic C-terminal tail modulates Escherichia coli gyrase activity

DNA topoisomerases manage chromosome supercoiling and organization in all cells. Gyrase, a prokaryotic type IIA topoisomerase, consumes ATP to introduce negative supercoils through a strand passage mechanism. All type IIA topoisomerases employ a similar set of catalytic domains for function; however, the activity and specificity of gyrase are augmented by a specialized DNA binding and wrapping element, termed the C-terminal domain (CTD), which is appended to its GyrA subunit. We have discovered that a nonconserved, acidic tail at the extreme C terminus of the Escherichia coli GyrA CTD has a dramatic and unexpected impact on gyrase function. Removal of the CTD tail enables GyrA to introduce writhe into DNA in the absence of GyrB, an activity exhibited by other GyrA orthologs, but not by wild-type E. coli GyrA. Strikingly, a "tail-less" gyrase holoenzyme is markedly impaired for DNA supercoiling capacity, but displays normal ATPase function. Our findings reveal that the E. coli GyrA tail regulates DNA wrapping by the CTD to increase the coupling efficiency between ATP turnover and supercoiling, demonstrating that CTD functions can be fine-tuned to control gyrase activity in a highly sophisticated manner.     

Conformational Dynamics of DNA Repair by Escherichia coli Endonuclease III

Escherichia coli endonuclease III (Endo III or Nth) is a DNA glycosylase with a broad substrate specificity for oxidized or reduced pyrimidine bases. Endo III possesses two types of activities: N-glycosylase (hydrolysis of the N-glycosidic bond) and AP lyase (elimination of the 3′-phosphate of the AP-site). We report a pre-steady-state kinetic analysis of structural rearrangements of the DNA substrates and uncleavable ligands during their interaction with Endo III. Oligonucleotide duplexes containing 5,6-dihydrouracil, a natural abasic site, its tetrahydrofuran analog, and undamaged duplexes carried fluorescent DNA base analogs 2-aminopurine and 1,3-diaza-2-oxophenoxazine as environment-sensitive reporter groups. The results suggest that Endo III induces several fast sequential conformational changes in DNA during binding, lesion recognition, and adjustment to a catalytically competent conformation. A comparison of two fluorophores allowed us to distinguish between the events occurring in the damaged and undamaged DNA strand. Combining our data with the available structures of Endo III, we conclude that this glycosylase uses a multistep mechanism of damage recognition, which likely involves Gln⁴¹ and Leu⁸¹ as DNA lesion sensors.     

Three-dimensional structure of the KdpFABC complex of Escherichia coli by electron tomography of two-dimensional crystals

The KdpFABC complex (Kdp) functions as a K⁺ pump in Escherichia coli and is a member of the family of P-type ATPases. Unlike other family members, Kdp has a unique oligomeric composition and is notable for segregating K⁺ transport and ATP hydrolysis onto separate subunits (KdpA and KdpB, respectively). We have produced two-dimensional crystals of the KdpFABC complex within reconstituted lipid bilayers and determined its three-dimensional structure from negatively stained samples using a combination of electron tomography and real-space averaging. The resulting map is at a resolution of 2.4 nm and reveals a dimer of Kdp molecules as the asymmetric unit; however, only the cytoplasmic domains are visible due to the lack of stain penetration within the lipid bilayer. The sizes of these cytoplasmic domains are consistent with Kdp and, using a pseudo-atomic model, we have described the subunit interactions that stabilize the Kdp dimer within the larger crystallographic array. These results illustrate the utility of electron tomography in structure determination of ordered assemblies, especially when disorder is severe enough to hamper conventional crystallographic analysis.     

Modeling of mammalian olfactory receptors and docking of odorants

Olfactory receptors (ORs) belong to the superfamily of G protein-coupled receptors (GPCRs), the second largest class of genes after those related to immunity, and account for about 3 % of mammalian genomes. ORs are present in all multicellular organisms and represent more than half the GPCRs in mammalian species (e.g., the mouse OR repertoire contains >1,000 functional genes). ORs are mainly expressed in the olfactory epithelium where they detect odorant molecules, but they are also expressed in a number of other cells, such as sperm cells, although their functions in these cells remain mostly unknown. It has recently been reported that ORs are present in tumoral tissues where they are expressed at different levels than in healthy tissues. A specific OR is over-expressed in prostate cancer cells, and activation of this OR has been shown to inhibit the proliferation of these cells. Odorant stimulation of some of these receptors results in inhibition of cell proliferation. Even though their biological role has not yet been elucidated, these receptors might constitute new targets for diagnosis and therapeutics. It is important to understand the activation mechanism of these receptors at the molecular level, in particular to be able to predict which ligands are likely to activate a particular receptor (‘deorphanization’) or to design antagonists for a given receptor. In this review, we describe the in silico methodologies used to model the three-dimensional (3D) structure of ORs (in the more general framework of GPCR modeling) and to dock ligands into these 3D structures.     

Molecular Modeling of Mycobacterium Tuberculosis DNA Gyrase and its Molecular Docking Study with Gatifloxacin Inhibitors

Abstract

Mycobacterium tuberculosis (Mt) is a leading cause of infectious disease in the world today. This outlook is aggravated by a growing number of Mt infections in individuals who are immunocompromised as a result of HIV infections. Thus, new and more potent anti-tuberculosis agents are necessary. Therefore, DNA gyrase was selected as a target enzyme to combat Mt. In this work, the first three-dimensional molecular model of the hypothetical structures for the Mycobacterium tuberculosis DNA gyrase (mtDNAg) was elucidated by a homology modeling method. In addition, the orientations and binding affinities of some gatifloxacin analogs with those new structures were investigated. Our findings could be helpful for the design of new more potent gatifloxacin analogs.     

Modulation of RAG/DNA complex by HSP70 in V(D)J recombination

V(D)J recombination, a site-specific gene rearrangement process, requires two RAG1 and RAG2 proteins specifically recognizing recombination signal sequences and forming DNA double-strand breaks. The broken DNA ends tightly bound to RAG proteins are joined by repair proteins. Here, we found that heat shock protein 70 was associated with RAG2 following two-step affinity chromatography purification. It was also co-immunoprecipitated with RAG2 in pro-B cells. Purified HSP70 protein disrupted RAG/DNA complexes assembled in vitro and also inhibited the V(D)J cleavage (both nick and hairpin formation) in a dose-dependent manner. This HSP70 action required ATP energy. These data suggest that HSP70 might play a crucial role in disassembling RAG/DNA complexes stably formed during V(D)J recombination.     

脱氮硫杆菌硫化合物载体SoxYZ蛋白的同源建模和结构分析

脱氮硫杆菌(Thiobacillus denitrificans)中的Sox蛋白在硫代谢过程中起着至关重要的作用,硫化合物需先与硫氧化基因族(sox)编码的蛋白质Sox YZ二聚体共价连接后才能与其他酶发生相互作用。利用同源建模法构建硫化合物载体Sox YZ蛋白的二聚体结构并验证了其合理性。二聚体相互作用分析发现Sox YZ蛋白的溶剂可及表面积(solvent accessible surface,SAS)为10 922.92,疏水率为50.85%;亚基Sox Y和Sox Z界面处共含有12个氢键和1个Π键来维持其三维结构的稳定性;二聚体表面呈现明显的正负电势互补,两亚基界面处氨基酸残基的VDW作用能和静电作用能分别为-80.925 13kcal/mol和-323.856 57kcal/mol,这说明静电作用是二聚体形成的主要驱动力;Sox Z亚基的残基Thr28、Arg31、Lys32、Ser64、Gly65、Val66、Ser67对Sox Y亚基活性位点构象的稳定有重要作用。     

Differences in Binding Sites of Two Melatonin Receptors Help to Explain Their Selectivity to Some Melatonin Analogs: A Molecular Modeling Study

Abstract

Numerous diseases have been linked to the malfunction of G-protein coupled receptors (GP-CRs). Their adequate treatment requires rational design of new high-affinity and high-selectivity drugs targeting these receptors. In this work, we report three-dimensional models of the human MT₁ and MT₂ melatonin receptors, members of the GPCR family. The models are based on the X-ray structure of bovine rhodopsin. The computational approach employs an original procedure for optimization of receptor-ligand structures. It includes rotation of one of the transmembrane α-helices around its axis with simultaneous assessment of quality of the resulting complexes according to a number of criteria we have developed for this purpose. The optimal geometry of the receptor-ligand binding is selected based on the analysis of complementarity of hydrophobic/hydrophilic properties between the ligand and its protein environment in the binding site. The elaborated “optimized” models are employed to explore the details of protein-ligand interactions for melatonin and a number of its analogs with known affinity to MT₁ and MT₂ receptors. The models permit rationalization of experimental data, including those that were not used in model building. The perspectives opened by the constructed models and by the optimization procedure in the design of new drugs are discussed.     

Molecular Modeling and docking of Wheat Hydroquinone Glucosyl transferase by using Hydroquinone,Phenyl phosphorodiamate and n-(n butyl) Phosphorothiocic Triamide as Inhibitors

In agriculture high urease activity during urea fertilization causes substantial environmental and economical problems by releasingabnormally large amount of ammonia into the atmosphere which leads to plant damage as well as ammonia toxicity. All over theworld, urea is the most widely applied nitrogen fertilizer. Due to the action of enzyme urease; urea nitrogen is lost as volatileammonia. For efficient use of nitrogen fertilizer, urease inhibitor along with the urea fertilizer is one of the best promisingstrategies. Urease inhibitors also provide an insight in understanding the mechanism of enzyme catalyzed reaction, the role ofvarious amino acids in catalytic activity present at the active site of enzyme and the importance of nickel to this metallo enzyme. Bykeeping it in view, the present study was designed to dock three urease inhibitors namely Hydroquinone (HQ), PhenylPhosphorodiamate (PPD) and N-(n-butyl) Phosphorothiocic triamide (NBPT) against Hydroquinone glucosyltransferase usingmolecular docking approach. The 3D structure of Hydroquinone glucosyltransferase was predicted using homology modelingapproach and quality of the structure was assured using Ramachandran plot. This study revealed important interactions amongthe urease inhibitors and Hydroquinone glucosyltransferase. Thus, it can be inferred that these inhibitors may serve as future antitoxic constituent against plant toxins.     

Solution structure of the IIAChitobiose-HPr complex of the N,N'-diacetylchitobiose branch of the Escherichia coli phosphotransferase system

The solution structure of the complex of enzyme IIA of the N,N'-diacetylchitobiose (Chb) transporter with the histidine phosphocarrier protein HPr has been solved by NMR. The IIA(Chb)-HPr complex completes the structure elucidation of representative cytoplasmic complexes for all four sugar branches of the bacterial phosphoryl transfer system (PTS). The active site His-89 of IIA(Chb) was mutated to Glu to mimic the phosphorylated state. IIA(Chb)(H89E) and HPr form a weak complex with a K(D) of ~0.7 mM. The interacting binding surfaces, concave for IIA(Chb) and convex for HPr, complement each other in terms of shape, residue type, and charge distribution, with predominantly hydrophobic residues, interspersed by some uncharged polar residues, located centrally, and polar and charged residues at the periphery. The active site histidine of HPr, His-15, is buried within the active site cleft of IIA(Chb) formed at the interface of two adjacent subunits of the IIA(Chb) trimer, thereby coming into close proximity with the active site residue, H89E, of IIA(Chb). A His89-P-His-15 pentacoordinate phosphoryl transition state can readily be modeled without necessitating any significant conformational changes, thereby facilitating rapid phosphoryl transfer. Comparison of the IIA(Chb)-HPr complex with the IIA(Chb)-IIB(Chb) complex, as well as with other cytoplasmic complexes of the PTS, highlights a unifying mechanism for recognition of structurally diverse partners. This involves generating similar binding surfaces from entirely different underlying structural elements, large interaction surfaces coupled with extensive redundancy, and side chain conformational plasticity to optimize diverse sets of intermolecular interactions.     

Interactions of replication versus repair DNA substrates with the Pol I DNA polymerases from Escherichia coli and Thermus aquaticus

Different DNA polymerases partition differently between replication and repair pathways. In this study we examine if two Pol I family polymerases from evolutionarily distant organisms also differ in their preferences for replication versus repair substrates. The DNA binding preferences of Klenow and Klentaq DNA polymerases, from Escherichia coli and Thermus aquaticus respectively, have been studied using a fluorescence competition binding assay. Klenow polymerase binds primed-template DNA (the replication substrate) with up to 50× higher affinity than it binds to nicked DNA, DNA with a 2 base single-stranded gap, blunt-ended DNA, or to a DNA end with a 3′ overhang. In contrast, Klentaq binds all of these DNAs almost identically, indicating that Klenow has a stronger ability to discriminate between replication and repair substrates than Klentaq. In contrast, both polymerases bind mismatched primed-template and blunt-ended DNA tighter than they bind matched primed-template DNA, suggesting that these two proteins may share a similar mechanism to identify mismatched DNA, despite the fact that Klentaq has no proofreading ability. In addition, the presence or absence of 5′- or 3′-phosphates has slightly different effects on DNA binding by the two polymerases, but again reinforce Klenow's more effective substrate discrimination capability.     

Mass spectrometric study of the Escherichia coli repressor proteins, IcIR and GcIR, and their complexes with DNA

In Escherichia coli, the IclR protein regulates both the aceBAK operon and its own synthesis. Database homology searches have identified many IclR-like proteins, now known as the IclR family, which can be identified by a conserved C-terminal region. We have cloned and purified one of these proteins, which we have named GclR (glyoxylate carboligase repressor). Although purification is straightforward, both the IclR and GclR proteins are difficult to manipulate, requiring high salt (up to 0.6 M KCl) for solubility. With the advent of nanospray ionization, we could transfer the proteins into much higher concentrations of volatile buffer than had been practical with ordinary electrospray. In 0.5 M ammonium bicarbonate buffer, both proteins were stable as tetramers, with a small amount of dimer. In a separate experiment, we found that IclR protein selected from a random pool a sequence which matched exactly that of the presumed binding region of the GclR protein, although IclR does not regulate the gcl gene. We designed a 29 bp synthetic DNA to which IclR and GclR bind, and with which we were able to form noncovalent DNA-protein complexes for further mass spectrometry analysis. These complexes were far more stable than the proteins alone, and we have evidence of a stoichiometry which has not been described previously with (protein monomer : dsDNA) = (4 : 1).