Nina Liu

Drug Discovery Using Chemical Systems Biology: Repositioning the Safe Medicine Comtan to Treat Multi-Drug and Extensively Drug Resistant Tuberculosis

The rise of multi-drug resistant (MDR) and extensively drug resistant (XDR) tuberculosis around the world, including in industrialized nations, poses a great threat to human health and defines a need to develop new, effective and inexpensive anti-tubercular agents. Previously we developed a chemical systems biology approach to identify off-targets of major pharmaceuticals on a proteome-wide scale. In this paper we further demonstrate the value of this approach through the discovery that existing commercially available drugs, prescribed for the treatment of Parkinsons disease, have the potential to treat MDR and XDR tuberculosis. These drugs, entacapone and tolcapone, are predicted to bind to the enzyme InhA and directly inhibit substrate binding. The prediction is validated by in vitro and InhA kinetic assays using tablets of Comtan, whose active component is entacapone. The minimal inhibition concentration (MIC99) of entacapone for Mycobacterium tuberculosis (M.tuberculosis) is approximately 260.0 µM, well below the toxicity concentration determined by an in vitro cytotoxicity model using a human neuroblastoma cell line. Moreover, kinetic assays indicate that Comtan inhibits InhA activity by 47.0% at an entacapone concentration of approximately 80 µM. Thus the active component in Comtan represents a promising lead compound for developing a new class of anti-tubercular therapeutics with excellent safety profiles. More generally, the protocol described in this paper can be included in a drug discovery pipeline in an effort to discover novel drug leads with desired safety profiles, and therefore accelerate the development of new drugs.

A Slow, Tight Binding Inhibitor of InhA, the Enoyl-Acyl Carrier Protein Reductase from Mycobacterium tuberculosis

InhA, the enoyl-ACP reductase in Mycobacterium tuberculosis is an attractive target for the development of novel drugs against tuberculosis, a disease that kills more than two million people each year. InhA is the target of the current first line drug isoniazid for the treatment of tuberculosis infections. Compounds that directly target InhA and do not require activation by the mycobacterial catalase-peroxidase KatG are promising candidates for treating infections caused by isoniazid-resistant strains. Previously we reported the synthesis of several diphenyl ethers with nanomolar affinity for InhA. However, these compounds are rapid reversible inhibitors of the enzyme, and based on the knowledge that long drug target residence times are an important factor for in vivo drug activity, we set out to generate a slow onset inhibitor of InhA using structure-based drug design. 2-(o-Tolyloxy)-5-hexylphenol (PT70) is a slow, tight binding inhibitor of InhA with a K1 value of 22 pm. PT70 binds preferentially to the InhA·NAD+ complex and has a residence time of 24 min on the target, which is 14,000 times longer than that of the rapid reversible inhibitor from which it is derived. The 1.8 Å crystal structure of the ternary complex between InhA, NAD+, and PT70 reveals the molecular details of enzyme-inhibitor recognition and supports the hypothesis that slow onset inhibition is coupled to ordering of an active site loop, which leads to the closure of the substrate-binding pocket.

A machine learning-based method to improve docking scoring functions and its application to drug repurposing.

Docking scoring functions are notoriously weak predictors of binding affinity. They typically assign a common set of weights to the individual energy terms that contribute to the overall energy score; however, these weights should be gene family dependent. In addition, they incorrectly assume that individual interactions contribute toward the total binding affinity in an additive manner. In reality, noncovalent interactions often depend on one another in a nonlinear manner. In this paper, we show how the use of support vector machines (SVMs), trained by associating sets of individual energy terms retrieved from molecular docking with the known binding affinity of each compound from high-throughput screening experiments, can be used to improve the correlation between known binding affinities and those predicted by the docking program eHiTS. We construct two prediction models: a regression model trained using IC(50) values from BindingDB, and a classification model trained using active and decoy compounds from the Directory of Useful Decoys (DUD). Moreover, to address the issue of overrepresentation of negative data in high-throughput screening data sets, we have designed a multiple-planar SVM training procedure for the classification model. The increased performance that both SVMs give when compared with the original eHiTS scoring function highlights the potential for using nonlinear methods when deriving overall energy scores from their individual components. We apply the above methodology to train a new scoring function for direct inhibitors of Mycobacterium tuberculosis (M.tb) InhA. By combining ligand binding site comparison with the new scoring function, we propose that phosphodiesterase inhibitors can potentially be repurposed to target M.tb InhA. Our methodology may be applied to other gene families for which target structures and activity data are available, as demonstrated in the work presented here.

Synthesis and in vitro antimycobacterial activity of B-ring modified diaryl ether InhA inhibitors

Previous structure-based design studies resulted in the discovery of alkyl substituted diphenyl ether inhibitors of InhA, the enoyl reductase from Mycobacterium tuberculosis. Compounds such as 5-hexyl-2-phenoxyphenol 19 are nM inhibitors of InhA and inhibit the growth of both sensitive and isoniazid-resistant strains of Mycobacterium tuberculosis with MIC(90) values of 1-2 microg/mL. However, despite their promising in vitro activity, these compounds have ClogP values of over 5. In efforts to reduce the lipophilicity of the compounds, and potentially enhance compound bioavailability, a series of B ring analogues of 19 were synthesized that contained either heterocylic nitrogen rings or phenyl rings having amino, nitro, amide, or piperazine functionalities. Compounds 3c, 3e, and 14a show comparable MIC(90) values to that of 19, but have improved ClogP values.

Synthesis and SAR studies of 1,4-benzoxazine MenB inhibitors: novel antibacterial agents against Mycobacterium tuberculosis.

Menaquinone is an essential component of the electron transport chain in many pathogens and consequently enzymes in the menaquinone biosynthesis pathway are potential drug targets for the development of novel antibacterial agents. In order to identify leads that target MenB, the 1,4-dihydroxy-2-naphthoyl-CoA synthase from Mycobacterium tuberculosis, a high-throughput screen was performed. Several 1,4-benzoxazines were identified in this screen and subsequent SAR studies resulted in the discovery of compounds with excellent antibacterial activity against M. tuberculosis H37Rv with MIC values as low as 0.6μg/ml. The 1,4-benzoxazine scaffold is thus a promising foundation for the development of antitubercular agents.

A Structural and Energetic Model for the Slow-Onset Inhibition of the Mycobacterium tuberculosis Enoyl-ACP Reductase InhA

Slow-onset enzyme inhibitors are of great interest for drug discovery programs since the slow dissociation of the inhibitor from the drug–target complex results in sustained target occupancy leading to improved pharmacodynamics. However, the structural basis for slow-onset inhibition is often not fully understood, hindering the development of structure-kinetic relationships and the rational optimization of drug-target residence time. Previously we demonstrated that slow-onset inhibition of the Mycobacterium tuberculosis enoyl-ACP reductase InhA correlated with motions of a substrate-binding loop (SBL) near the active site. In the present work, X-ray crystallography and molecular dynamics simulations have been used to map the structural and energetic changes of the SBL that occur upon enzyme inhibition. Helix-6 within the SBL adopts an open conformation when the inhibitor structure or binding kinetics is substrate-like. In contrast, slow-onset inhibition results in large-scale local refolding in which helix-6 adopts a closed conformation not normally populated during substrate turnover. The open and closed conformations of helix-6 are hypothesized to represent the EI and EI* states on the two-step induced-fit reaction coordinate for enzyme inhibition. These two states were used as the end points for nudged elastic band molecular dynamics simulations resulting in two-dimensional potential energy profiles that reveal the barrier between EI and EI*, thus rationalizing the binding kinetics observed with different inhibitors. Our findings indicate that the structural basis for slow-onset kinetics can be understood once the structures of both EI and EI* have been identified, thus providing a starting point for the rational control of enzyme–inhibitor binding kinetics.

Time‐Dependent Diaryl Ether Inhibitors of InhA: Structure–Activity Relationship Studies of Enzyme Inhibition, Antibacterial Activity, and in vivo Efficacy

The diaryl ethers are a novel class of antituberculosis drug candidates that inhibit InhA, the enoyl‐ACP reductase involved in the fatty acid biosynthesis (FASII) pathway, and have antibacterial activity against both drug‐sensitive and drug‐resistant strains of Mycobacterium tuberculosis. In the present work, we demonstrate that two time‐dependent B‐ring modified diaryl ether InhA inhibitors have antibacterial activity in a mouse model of TB infection when delivered by intraperitoneal injection. We propose that the efficacy of these compounds is related to their residence time on the enzyme, and to identify structural features that modulate drug–target residence time in this system, we have explored the inhibition of InhA by a series of B‐ring modified analogues. Seven ortho‐substituted compounds were found to be time‐dependent inhibitors of InhA, where the slow step leading to the final enzyme–inhibitor complex (EI*) is thought to correlate with closure and ordering of the InhA substrate binding loop. A detailed mechanistic understanding of the molecular basis for residence time in this system will facilitate the development of InhA inhibitors with improved in vivo activity.

Mechanism and inhibition of the FabI enoyl-ACP reductase from Burkholderia pseudomallei

OBJECTIVES As an initial step in developing novel antibacterials against Burkholderia pseudomallei, we have characterized the FabI enoyl-ACP reductase homologues in the type II fatty acid biosynthesis pathway from this organism and performed an initial enzyme inhibition study. METHODS A BLAST analysis identified two FabI enoyl-ACP reductase homologues, bpmFabI-1 and bpmFabI-2, in the B. pseudomallei genome, which were cloned, overexpressed in Escherichia coli and purified. Steady-state kinetics was used to determine the reaction mechanism and the sensitivity of bpmFabI-1 to four diphenyl ether FabI inhibitors. The antibacterial activity of the inhibitors was assessed using a wild-type strain of Burkholderia thailandensis (E264) and an efflux pump mutant (Bt38). RESULTS Consistent with its annotation as an enoyl-ACP reductase, bpmFabI-1 catalysed the NADH-dependent reduction of 2-trans-dodecenoyl-CoA via a sequential Bi Bi mechanism. In contrast, bpmFabI-2 was inactive with all substrates tested and only bpmfabI-1 was transcriptionally active under the growth conditions employed. The sensitivity of bpmFabI-1 to four diphenyl ethers was evaluated and in each case the compounds were slow-onset inhibitors with K(i) values of 0.5-2 nM. In addition, triclosan and PT01 had MIC values of 30 and 70 mg/L for B. pseudomallei as well as a wild-type strain of B. thailandensis (E264), but MIC values of <1 mg/L for the efflux pump mutant Bt38. A reduction in MIC values was also observed for the pump mutant strain with the other diphenyl ethers. CONCLUSIONS Provided that efflux can be circumvented, bpmFabI-1 is a suitable target for drug discovery.

Mechanism of the intramolecular Claisen condensation reaction catalyzed by MenB, a crotonase superfamily member.

MenB, the 1,4-dihydroxy-2-naphthoyl-CoA synthase from the bacterial menaquinone biosynthesis pathway, catalyzes an intramolecular Claisen condensation (Dieckmann reaction) in which the electrophile is an unactivated carboxylic acid. Mechanistic studies on this crotonase family member have been hindered by partial active site disorder in existing MenB X-ray structures. In the current work the 2.0 Å structure of O-succinylbenzoyl-aminoCoA (OSB-NCoA) bound to the MenB from Escherichia coli provides important insight into the catalytic mechanism by revealing the position of all active site residues. This has been accomplished by the use of a stable analogue of the O-succinylbenzoyl-CoA (OSB-CoA) substrate in which the CoA thiol has been replaced by an amine. The resulting OSB-NCoA is stable, and the X-ray structure of this molecule bound to MenB reveals the structure of the enzyme-substrate complex poised for carbon-carbon bond formation. The structural data support a mechanism in which two conserved active site Tyr residues, Y97 and Y258, participate directly in the intramolecular transfer of the substrate α-proton to the benzylic carboxylate of the substrate, leading to protonation of the electrophile and formation of the required carbanion. Y97 and Y258 are also ideally positioned to function as the second oxyanion hole required for stabilization of the tetrahedral intermediate formed during carbon-carbon bond formation. In contrast, D163, which is structurally homologous to the acid-base catalyst E144 in crotonase (enoyl-CoA hydratase), is not directly involved in carbanion formation and may instead play a structural role by stabilizing the loop that carries Y97. When similar studies were performed on the MenB from Mycobacterium tuberculosis, a twisted hexamer was unexpectedly observed, demonstrating the flexibility of the interfacial loops that are involved in the generation of the novel tertiary and quaternary structures found in the crotonase superfamily. This work reinforces the utility of using a stable substrate analogue as a mechanistic probe in which only one atom has been altered leading to a decrease in α-proton acidity.

Structure of the Yersinia pestis FabV enoyl-ACP reductase and its interaction with two 2-pyridone inhibitors

The recently discovered FabV enoyl-ACP reductase, which catalyzes the last step of the bacterial fatty acid biosynthesis (FAS-II) pathway, is a promising but unexploited drug target against the reemerging pathogen Yersinia pestis. The structure of Y. pestis FabV in complex with its cofactor reveals that the enzyme features the common architecture of the short-chain dehydrogenase reductase superfamily, but contains additional structural elements that are mostly folded around the usually flexible substrate-binding loop, thereby stabilizing it in a very tight conformation that seals the active site. The structures of FabV in complex with NADH and two newly developed 2-pyridone inhibitors provide insights for the development of new lead compounds, and suggest a mechanism by which the substrate-binding loop opens to admit the inhibitor, a motion that could also be coupled to the interaction of FabV with the acyl-carrier protein substrate.