skip to main content


Title: Computational Investigation of a Series of Small Molecules as Lead Compounds for Lysyl hydroxylase-2 (LH2) Inhibition

The catalytic function of Lysyl hydroxylase-2 (LH2), a member of the Fe(II)/αKG-dependent oxygenase superfamily, is to catalyze the hydroxylation of lysine to hydroxylysine in collagen, resulting in stable hydroxylysine aldehyde-derived collagen cross-links (HLCCs). Reports show that high amounts of LH2 lead to the accumulation of HLCCs, causing fibrosis and specific types of cancer metastasis. Some members of the Fe(II)/αKG-dependent family have also been reported to have intramolecular O2 tunnels, which aid in transporting one of the required co-substrates into the active site. While LH2 can be a promising target to combat these diseases, efficacious inhibitors are still lacking. We have used computational simulations to investigate a series of forty-four small molecules as lead compounds for LH2 inhibition. Tunneling analyses indicate the existence of several intra-molecular tunnels. The lengths of the calculated O2-transporting tunnels in holoenzymes are relatively longer than the apoenzyme suggesting that the ligands may affect the enzyme's structure and possibly block (at least partially) the tunnels. The sequence alignment analysis between LH enzymes from different organisms shows that all the amino acid residues with the highest occurrence rate in the oxygen tunnels are conserved. Our results suggest that the enolate form of diketone compounds establishes stronger interactions with the Fe(II) in the active site. Branching the enolate compounds with functional groups such as phenyl and pyridinyl enhances the interaction with various residues around the active site. Our results provide information about possible leads for further LH2 inhibition design and development.

 
more » « less
Award ID(s):
2117247
PAR ID:
10558183
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ; ; ; ;
Publisher / Repository:
ACS
Date Published:
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Oxalate decarboxylase is an Mn- and O2-dependent enzyme in the bicupin superfamily that catalyzes the redox-neutral disproportionation of the oxalate monoanion to form carbon dioxide and formate. Its best-studied isozyme is from Bacillus subtilis where it is stress-induced under low pH conditions. Current mechanistic schemes assume a monodentate binding mode of the substrate to the N-terminal active site Mn ion to make space for a presumed O2 molecule, despite the fact that oxalate generally prefers to bind bidentate to Mn. We report on X-band 13C-electron nuclear double resonance (ENDOR) experiments on 13C-labeled oxalate bound to the active-site Mn(II) in wild-type oxalate decarboxylase at high pH, the catalytically impaired W96F mutant enzyme at low pH, and Mn(II) in aqueous solution. The ENDOR spectra of these samples are practically identical, which shows that the substrate binds bidentate (κO, κO’) to the active site Mn(II) ion. Domain-based local pair natural orbital coupled cluster singles and doubles (DLPNO-CCSD) calculations of the expected 13C hyperfine coupling constants for bidentate bound oxalate predict ENDOR spectra in good agreement with the experiment, supporting bidentate bound substrate. Geometry optimization of a substrate-bound minimal active site model by density functional theory shows two possible substrate coordination geometries, bidentate and monodentate. The bidentate structure is energetically preferred by ~4.7 kcal/mol. Our results revise a long-standing hypothesis regarding substrate binding in the enzyme and suggest that dioxygen does not bind to the active site Mn ion after substrate binds. The results are in agreement with our recent mechanistic hypothesis of substrate activation via a long-range electron transfer process involving the C-terminal Mn ion.

     
    more » « less
  2. Cleavage of aromatic carbon–chlorine bonds is critical for the degradation of toxic industrial compounds. Here, we solved the X-ray crystal structure of chlorothalonil dehalogenase (Chd) from Pseudomonas sp. CTN-3, with 15 of its N-terminal residues truncated (Chd T ), using single-wavelength anomalous dispersion refined to 1.96 Å resolution. Chd has low sequence identity (<15%) compared with all other proteins whose structures are currently available, and to the best of our knowledge, we present the first structure of a Zn(II)-dependent aromatic dehalogenase that does not require a coenzyme. Chd T forms a “head-to-tail” homodimer, formed between two α-helices from each monomer, with three Zn(II)-binding sites, two of which occupy the active sites, whereas the third anchors a structural site at the homodimer interface. The catalytic Zn(II) ions are solvent-accessible via a large hydrophobic (8.5 × 17.8 Å) opening to bulk solvent and two hydrophilic branched channels. Each active-site Zn(II) ion resides in a distorted trigonal bipyramid geometry with His 117 , His 257 , Asp 116 , Asn 216 , and a water/hydroxide as ligands. A conserved His residue, His 114 , is hydrogen-bonded to the Zn(II)-bound water/hydroxide and likely functions as the general acid-base. We examined substrate binding by docking chlorothalonil (2,4,5,6-tetrachloroisophtalonitrile, TPN) into the hydrophobic channel and observed that the most energetically favorable pose includes a TPN orientation that coordinates to the active-site Zn(II) ions via a CN and that maximizes a π–π interaction with Trp 227 . On the basis of these results, along with previously reported kinetics data, we propose a refined catalytic mechanism for Chd-mediated TPN dehalogenation. 
    more » « less
  3. Small molecule inhibitors of snake venom metalloproteinases (SVMPs) could provide a means to rapidly halt the progression of local tissue damage following viperid snake envenomations. In this study, we examine the ability of candidate compounds based on a pentacyclic triterpene skeleton to inhibit SVMPs. We leverage molecular dynamics simulations to estimate the free energies of the candidate compounds for binding to BaP1, a P-I type SVMP, and compare these results with experimental assays of proteolytic activity inhibition in a homologous enzyme (Batx-I). Both simulation and experiment suggest that betulinic acid is the most active candidate, with the simulations predicting a standard binding free energy of Δ G ∘ = − 11.0 ± 1.4 kcal/mol. The simulations also reveal the atomic interactions that underlie binding between the triterpenic acids and BaP1, most notably the electrostatic interaction between carboxylate groups of the compounds and the zinc cofactor of BaP1. Together, our simulations and experiments suggest that occlusion of the S1 ′ subsite is essential for inhibition of proteolytic activity. While all active compounds make hydrophobic contacts in the S1 ′ site, β -boswellic acid, with its distinct carboxylate position, does not occlude the S1 ′ site in simulation and exhibits negligible activity in experiment. 
    more » « less
  4. Hydroxyl radical (•OH) is produced in soils from oxidation of reduced iron (Fe(II)) by dissolved oxygen (O2) and can oxidize dissolved organic carbon (DOC) to carbon dioxide (CO2). Understanding the role of •OH on CO2 production in soils requires knowing whether Fe(II) production or O2 supply to soils limits •OH production. To test the relative importance of Fe(II) production versus O2 supply, we measured changes in Fe(II) and O2 and in situ •OH production during simulated precipitation events and during common, waterlogged conditions in mesocosms from two landscape ages and the two dominant vegetation types of the Arctic. The balance of Fe(II) production and consumption controlled •OH production during precipitation events that supplied O2 to the soils. During static, waterlogged conditions, •OH production was controlled by O2 supply because Fe(II) production was higher than its consumption (oxidation) by O2. An average precipitation event (4 mm) resulted in 200 µmol •OH m−2 per day produced compared to 60 µmol •OH m−2 per day produced during waterlogged conditions. These findings suggest that the oxidation of DOC to CO2 by •OH in arctic soils, a process potentially as important as microbial respiration of DOC in arctic surface waters, will depend on the patterns and amounts of rainfall that oxygenate the soil. 
    more » « less
  5. The iron-containing heterodimeric MbnBC enzyme complex plays a central role in the biosynthesis of methanobactins (Mbns), ribosomally synthesized, posttranslationally modified natural products that bind copper with high affinity. MbnBC catalyzes a four-electron oxidation of a cysteine residue in its precursor-peptide substrate, MbnA, to an oxazolone ring and an adjacent thioamide group. Initial studies of MbnBC indicated the presence of both diiron and triiron species, complicating identification of the catalytically active species. Here, we present evidence through activity assays combined with electron paramagnetic resonance (EPR) and Mössbauer spectroscopic analysis that the active species is a mixed-valent, antiferromagnetically coupled Fe(II)Fe(III) center. Consistent with this assignment, heterologous expression of the MbnBC complex in culture medium containing less iron yielded purified protein with less bound iron but greater activity in vitro. The maximally activated MbnBC prepared in this manner could modify both cysteine residues in MbnA, in contrast to prior findings that only the first cysteine could be processed. Site-directed mutagenesis and multiple crystal structures clearly identify the two essential Fe ions in the active cluster as well as the location of the previously detected third Fe site. Moreover, structural modeling indicates a role for MbnC in recognition of the MbnA leader peptide. These results add a biosynthetic oxidative rearrangement reaction to the repertoire of nonheme diiron enzymes and provide a foundation for elucidating the MbnBC mechanism. 
    more » « less