Apo and nickel-bound forms of the Pyrococcus horikoshii species of the metalloregulatory protein: NikR characterized by molecular dynamics simulations

Dinickel enzyme evolved to metabolize the pharmaceutical metformin and its implications for wastewater and human microbiomes

Metformin is the first-line treatment for type II diabetes patients and a pervasive pollutant with more than 180 million kg ingested globally and entering wastewater. The drug's direct mode of action is currently unknown but is linked to effects on gut microbiomes and may involve specific gut microbial reactions to the drug. In wastewater treatment plants, metformin is known to be transformed by microbes to guanylurea, although genes encoding this metabolism had not been elucidated. In the present study, we revealed the function of two genes responsible for metformin decomposition (mfmA and mfmB) found in isolated bacteria from activated sludge. MfmA and MfmB form an active heterocomplex (MfmAB) and are members of the ureohydrolase protein superfamily with binuclear metal-dependent activity. MfmAB is nickel-dependent and catalyzes the hydrolysis of metformin to dimethylamine and guanylurea with a catalytic efficiency (kcat/KM) of 9.6 × 103 M-1s-1 and KM for metformin of 0.82 mM. MfmAB shows preferential activity for metformin, being able to discriminate other close substrates by several orders of magnitude. Crystal structures of MfmAB show coordination of binuclear nickel bound in the active site of the MfmA subunit but not MfmB subunits, indicating that MfmA is the active site for the MfmAB complex. Mutagenesis of residues conserved in the MfmA active site revealed those critical to metformin hydrolase activity and its small substrate binding pocket allowed for modeling of bound metformin. This study characterizes the products of the mfmAB genes identified in wastewater treatment plants on three continents, suggesting that metformin hydrolase is widespread globally in wastewater.

Antibacterial properties and urease suppression ability of Lactobacillus inhibit the development of infectious urinary stones caused by Proteus mirabilis

Infectious urolithiasis is a type of urolithiasis, that is caused by infections of the urinary tract by bacteria producing urease such as Proteus mirabilis. Lactobacillus spp. have an antagonistic effect against many pathogens by secreting molecules, including organic acids. The aim of the study was to determine the impact of Lactobacillus strains isolated from human urine on crystallization of urine components caused by P. mirabilis by measuring bacterial viability (CFU/mL), pH, ammonia release, concentration of crystallized salts and by observing crystals by phase contrast microscopy. Moreover, the effect of lactic acid on the activity of urease was examined by the kinetic method and in silico study. In the presence of selected Lactobacillus strains, the crystallization process was inhibited. The results indicate that one of the mechanisms of this action was the antibacterial effect of Lactobacillus, especially in the presence of L. gasseri, where ten times less P. mirabilis bacteria was observed, compared to the control. It was also demonstrated that lactic acid inhibited urease activity by a competitive mechanism and had a higher binding affinity to the enzyme than urea. These results demonstrate that Lactobacillus and lactic acid have a great impact on the urinary stones development, which in the future may help to support the treatment of this health problem.

In silico profiling and structural insights of zinc metal ion on O6-methylguanine methyl transferase and its interactions using molecular dynamics approach

O6-methylguanine DNA methyl transferase (MGMT) is a metalloenzyme participating in the repair of alkylated DNA. In this research, we performed a comparative study for evaluating the impact of zinc metal ion on the behavior and interactions of MGMT in the both enzymatic forms of apo MGMT and holo MGMT. DNA and proliferating cell nuclear antigen (PCNA), as partners of MGMT, were utilized to evaluate molecular interactions by virtual microscopy of molecular dynamics simulation. The stability and conformational alterations of each forms (apo and holo) MGMT-PCNA, and (apo and holo) MGMT-DNA complexes were calculated by MM/PBSA method. A total of seven systems including apo MGMT, holo MGMT, free PCNA, apo MGMT-PCNA, holo MGMT-PCNA, apo MGMT-DNA, and holo MGMT-DNA complexes were simulated. In this study, we found that holo MGMT was more stable and had better folding and functional properties than that of apo MGMT. Simulation analysis of (apo and holo) MGMT-PCNA complexes displayed that the sequences of the amino acids involved in the interactions were different in the two forms of MGMT. The important amino acids of holo MGMT involved in its interaction with PCNA included E92, K101, A119, G122, N123, P124, and K125, whereas the important amino acids of apo MGMT included R128, R135, S152, N157, Y158, and L162. Virtual microscopy of molecular dynamics simulation showed that the R128 and its surrounding residues were important amino acids involved in the interaction of holo MGMT with DNA that was exactly consistent with X-ray crystallography structure. In the apo form of the protein, the N157 and its surrounding residues were important amino acids involved in the interaction with DNA. The binding free energies of - 387.976, - 396.226, - 622.227, and - 617.333 kcal/mol were obtained for holo MGMT-PCNA, apo MGMT-PCNA, holo MGMT-DNA, and apo MGMT-DNA complexes, respectively. The principle result of this research was that the area of molecular interactions differed between the two states of MGMT. Therefore, in investigations of metalloproteins, the metal ion must be preserved in their structures. Finally, it is recommended to use the holo form of metalloproteins in in vitro and in silico researches.

Metal binding mediated conformational change of XPA protein:a potential cytotoxic mechanism of nickel in the nucleotide excision repair

Nucleotide excision repair (NER) is a pivotal life process for repairing DNA nucleotide mismatch caused by chemicals, metal ions, radiation, and other factors. As the initiation step of NER, the xeroderma pigmentosum complementation group A protein (XPA) recognizes damaged DNA molecules, and recruits the replication protein A (RPA), another important player in the NER process. The stability of the Zn(2+)-chelated Zn-finger domain of XPA center core portion (i.e., XPA98-210) is the foundation of its biological functionality, while the displacement of the Zn(2+) by toxic metal ions (such as Ni(2+), a known human carcinogen and allergen) may impair the effectiveness of NER and hence elevate the chance of carcinogenesis. In this study, we first calculated the force field parameters for the bonded model in the metal center of the XPA98-210 system, showing that the calculated results, including charges, bonds, angles etc., are congruent with previously reported results measured by spectrometry experiments and quantum chemistry computation. Then, comparative molecular dynamics simulations using these parameters revealed the changes in the conformation and motion mode of XPA98-210 Zn-finger after the substitution of Zn(2+) by Ni(2+). The results showed that Ni(2+) dramatically disrupted the relative positions of the four Cys residues in the Zn-finger structure, forcing them to collapse from a tetrahedron into an almost planar structure. Finally, we acquired the binding mode of XPA98-210 with its ligands RPA70N and DNA based on molecular docking and structural alignment. We found that XPA98-210's Zn-finger domain primarily binds to a V-shaped cleft in RPA70N, while the cationic band in its C-terminal subdomain participates in the recognition of damaged DNA. In addition, this article sheds light on the multi-component interaction pattern among XPA, DNA, and other NER-related proteins (i.e., RPA70N, RPA70A, RPA70B, RPA70C, RPA32, and RPA14) based on previously reported structural biology information. Thus, we derived a putative cytotoxic mechanism associated with the nickel ion, where the Ni(2+) disrupts the conformation of the XPA Zn-finger, directly weakening its interaction with RPA70N, and thus lowering the effectiveness of the NER process. In sum, this work not only provides a theoretical insight into the multi-protein interactions involved in the NER process and potential cytotoxic mechanism associated with Ni(2+) binding in XPA, but may also facilitate rational anti-cancer drug design based on the NER mechanism.

Genome-wide transcriptional response of the archaeon Thermococcus gammatolerans to cadmium

Thermococcus gammatolerans, the most radioresistant archaeon known to date, is an anaerobic and hyperthermophilic sulfur-reducing organism living in deep-sea hydrothermal vents. Knowledge of mechanisms underlying archaeal metal tolerance in such metal-rich ecosystem is still poorly documented. We showed that T. gammatolerans exhibits high resistance to cadmium (Cd), cobalt (Co) and zinc (Zn), a weaker tolerance to nickel (Ni), copper (Cu) and arsenate (AsO₄) and that cells exposed to 1 mM Cd exhibit a cellular Cd concentration of 67 µM. A time-dependent transcriptomic analysis using microarrays was performed at a non-toxic (100 µM) and a toxic (1 mM) Cd dose. The reliability of microarray data was strengthened by real time RT-PCR validations. Altogether, 114 Cd responsive genes were revealed and a substantial subset of genes is related to metal homeostasis, drug detoxification, re-oxidization of cofactors and ATP production. This first genome-wide expression profiling study of archaeal cells challenged with Cd showed that T. gammatolerans withstands induced stress through pathways observed in both prokaryotes and eukaryotes but also through new and original strategies. T. gammatolerans cells challenged with 1 mM Cd basically promote: 1) the induction of several transporter/permease encoding genes, probably to detoxify the cell; 2) the upregulation of Fe transporters encoding genes to likely compensate Cd damages in iron-containing proteins; 3) the induction of membrane-bound hydrogenase (Mbh) and membrane-bound hydrogenlyase (Mhy2) subunits encoding genes involved in recycling reduced cofactors and/or in proton translocation for energy production. By contrast to other organisms, redox homeostasis genes appear constitutively expressed and only a few genes encoding DNA repair proteins are regulated. We compared the expression of 27 Cd responsive genes in other stress conditions (Zn, Ni, heat shock, γ-rays), and showed that the Cd transcriptional pattern is comparable to other metal stress transcriptional responses (Cd, Zn, Ni) but not to a general stress response.

Metal site occupancy and allosteric switching in bacterial metal sensor proteins

All prokaryotes encode a panel of metal sensor or metalloregulatory proteins that govern the expression of genes that allows an organism to quickly adapt to toxicity or deprivation of both biologically essential transition metal ions, e.g., Zn, Cu, Fe, and heavy metal pollutants. As such, metal sensor proteins can be considered arbiters of intracellular transition metal bioavailability and thus potentially control the metallation state of the metalloproteins in the cell. Metal sensor proteins are specialized allosteric proteins that regulate transcription as a result direct binding of one or two cognate metal ions, to the exclusion of all others. In most cases, the binding of the cognate metal ion induces a structural change in a protein oligomer that either activates or inhibits operator DNA binding. A quantitative measure of the degree to which a particular metal drives metalloregulation of operator DNA-binding is the allosteric coupling free energy, ΔGc. In this review, we summarize recent work directed toward understanding metal occupancy and metal selectivity of these allosteric switches in selected families of metal sensor proteins and examine the structural origins of ΔGc in the functional context a thermodynamic "set-point" model of intracellular metal homeostasis.

Helicobacter pylori NikR protein exhibits distinct conformations when bound to different promoters

Helicobacter pylori NikR (HpNikR) is a ribbon-helix-helix (RHH) DNA-binding protein that binds to several different promoter regions. The binding site sequences are not absolutely conserved. The ability of HpNikR to discriminate specific DNA sites resides partly in its nine-amino acid N-terminal arm. Previously, indirect evidence indicated that the arm exists in different conformations when HpNikR is bound to the nixA and ureA promoters. Here, we directly examined HpNikR conformation when it was bound to nixA and ureA DNA fragments by tethering (S)-1{[bis(carboxymethyl)amino]methyl}-2-{4-[(2-bromoacetyl)amino]phenylethyl}(carboxymethyl)amino]acetic acid, iron(III) to different positions in the N-terminal arm and RHH DNA binding domain. Different cleavage patterns at each promoter directly demonstrated that both the RHH domain and the arm adopt different conformations on the nixA and ureA promoters. Additionally, the two RHH domain dimers of the HpNikR tetramer are in distinct conformations at ureA. Site-directed mutagenesis identified an interchain salt bridge (Lys(48)-Glu(47')) in the RHH domain remote from the DNA binding interface that is required for high affinity binding to ureA but not nixA. Finally, DNA affinity measurements of wild-type HpNikR and a salt bridge mutant (K48A) to hybrid nixA-ureA promoters demonstrated that inverted repeat half-sites, spacers, and flanking DNA are all required for sequence-specific DNA binding by HpNikR. Notably, the spacer region made the largest contribution to DNA affinity. HpNikR exhibits a substantially expanded regulon compared with other NikR proteins. The results presented here provide a molecular basis for understanding regulatory network expansion by NikR as well as other prokaryotic regulatory proteins.

Holo-Ni(II)HpNikR is an asymmetric tetramer containing two different nickel-binding sites

The metalloregulatory protein NikR from Helicobacter pylori (HpNikR) is a master regulator of gene expression which both activates and represses specific genes in response to nickel availability. Here, we report the first crystal structure (at 2.37 Å resolution) of Ni(II)HpNikR prepared directly from the holo protein. The protein contains four nickel ions located in two distinct coordination environments. Two nickel ions are bound to sites in a four-coordinate square-planar geometry as predicted on the basis of the structures of NikR from Escherichia coli and Pyrococcus horikoshii . The remaining two nickel ions are bound to sites with unexpected 5- or 6-coordination geometries which were previously thought to be involved in nickel incorporation into the protein. The nickel with 5-/6-coordination geometry utilizes three histidines from two separate monomeric HpNikR units along with two or three water molecules as ligands. The spatial location of the nickel in the 5-/6-coordinate site is within approximately 5 Å of the expected site if a 4-coordinate square-planar geometry occurred. Two of the histidines that participate as ligands in the 5-/6-coordinate site would also participate as ligands if the 4-coordinate site was occupied, making it impossible for both sites to be occupied simultaneously. DFT calculations show that the 5-/6-coordinate geometries are energetically favorable when the local protein environment is included in the calculations. The presence of two distinct coordination environments in HpNikR is suggested to be related to the specificity and binding affinity of this transcription factor for DNA.

Structural Survey of Zinc Containing Proteins and the Development of the Zinc AMBER Force Field (ZAFF)

Currently the Protein Data Bank (PDB) contains over 18,000 structures that contain a metal ion including Na, Mg, K, Ca, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, Pd, Ag, Cd, Ir, Pt, Au, and Hg. In general, carrying out classical molecular dynamics (MD) simulations of metalloproteins is a convoluted and time consuming process. Herein, we describe MCPB (Metal Center Parameter Builder), which allows one, to conveniently and rapidly incorporate metal ions using the bonded plus electrostatics model (Hoops et al., J. Am. Chem. Soc. 1991, 113, 8262-8270) into the AMBER Force Field (FF). MCPB was used to develop a Zinc FF, ZAFF, which is compatible with the existing AMBER FFs. The PDB was mined for all Zn containing structures with most being tetrahedrally bound. The most abundant primary shell ligand combinations were extracted and FFs were created. These include Zn bound to CCCC, CCCH, CCHH, CHHH, HHHH, HHHO, HHOO, HOOO, HHHD, and HHDD (O = water and the remaining are 1 letter amino acid codes). Bond and angle force constants and RESP charges were obtained from B3LYP/6-31G* calculations of model structures from the various primary shell combinations. MCPB and ZAFF can be used to create FFs for MD simulations of metalloproteins to study enzyme catalysis, drug design and metalloprotein crystal refinement.