Ana Trapero1, Angela Pacitto2, Vinayak Singh3, Mohamad Sabbah1, Anthony G Coyne1, Valerie Mizrahi3, Tom L Blundell2, David B Ascher2,4, Chris Abell1. 1. Department of Chemistry , University of Cambridge , Lensfield Road , Cambridge CB2 1EW , United Kingdom. 2. Department of Biochemistry , University of Cambridge , 80 Tennis Court Road , Cambridge CB2 1GA , United Kingdom. 3. MRC/NHLS/UCT Molecular Mycobacteriology Research Unit & DST/NRF Centre of Excellence for Biomedical TB Research, Institute of Infectious Disease and Molecular Medicine and Division of Medical Microbiology, Faculty of Health Sciences , University of Cape Town , Rondebosch 7701 , Cape Town , South Africa. 4. Department of Biochemistry and Molecular Biology, Bio21 Institute , University of Melbourne , 30 Flemington Road , Parkville , Victoria 3052 , Australia.
Abstract
Tuberculosis (TB) remains a major cause of mortality worldwide, and improved treatments are needed to combat emergence of drug resistance. Inosine 5'-monophosphate dehydrogenase (IMPDH), a crucial enzyme required for de novo synthesis of guanine nucleotides, is an attractive TB drug target. Herein, we describe the identification of potent IMPDH inhibitors using fragment-based screening and structure-based design techniques. Screening of a fragment library for Mycobacterium thermoresistible ( Mth) IMPDH ΔCBS inhibitors identified a low affinity phenylimidazole derivative. X-ray crystallography of the Mth IMPDH ΔCBS-IMP-inhibitor complex revealed that two molecules of the fragment were bound in the NAD binding pocket of IMPDH. Linking the two molecules of the fragment afforded compounds with more than 1000-fold improvement in IMPDH affinity over the initial fragment hit.
Tuberculosis (TB) remains a major cause of mortality worldwide, and improved treatments are needed to combat emergence of drug resistance. Inosine 5'-monophosphate dehydrogenase (IMPDH), a crucial enzyme required for de novo synthesis of guanine nucleotides, is an attractive TB drug target. Herein, we describe the identification of potent IMPDH inhibitors using fragment-based screening and structure-based design techniques. Screening of a fragment library for Mycobacterium thermoresistible ( Mth) IMPDH ΔCBS inhibitors identified a low affinity phenylimidazole derivative. X-ray crystallography of the Mth IMPDH ΔCBS-IMP-inhibitor complex revealed that two molecules of the fragment were bound in the NAD binding pocket of IMPDH. Linking the two molecules of the fragment afforded compounds with more than 1000-fold improvement in IMPDH affinity over the initial fragment hit.
Tuberculosis
(TB) is a contagious infectious disease caused by Mycobacterium
tuberculosis (Mtb), which
can be transmitted through the air as droplets. The infection predominantly
affects the lungs, but it can spread to other parts of the body, especially
in patients with a suppressed immune system.The World Health
Organization (WHO) has estimated that nearly one-third of the world’s population
is infected with Mtb, leading to 1.8 million TB deaths
in 2015.[1] Although there has been a slow
decline in new TB cases and TB-related deaths in recent years, the
emergence and spread of multidrug resistant (MDR) and extensively
drug resistant (XDR) strains of Mtb has increased
the threat that this disease poses for global public health. According
to the WHO, approximately 480,000 cases of MDR-TB emerged in 2015,
and the cure rate of those patients was only 50%.[1]Current TB treatments require combinations of four
first-line drugs,
isoniazid, rifampicin, ethambutol, pyrazinamide, and streptomycin,
which must be taken for six months or longer.[2] Resistant strains are not susceptible to the standard drugs, and
although MDR-TB is treatable using second-line drugs, such treatments
have a number severe side effects.[3] Consequently,
there is an urgent need for the development of novel and more effective
drugs for the treatment of drug resistant TB.Inosine-5′-monophosphate
dehydrogenase (IMPDH, E.C. 1.1.1.205)
has received considerable interest in recent years as an important
target enzyme for immunosuppressive,[4] anticancer,[5,6] and antiviral drugs.[7] Most recently,
IMPDH has emerged as a promising antimicrobial drug target.[8−11]IMPDH catalyzes the first unique step in the de novo synthesis of guanine nucleotides, the oxidation of inosine 5′-monophosphate
(IMP) to xanthosine 5′-monophosphate (XMP) with the concomitant
reduction of the cofactor nicotinamide adenine dinucleotide (NAD+) to NADH (Figure ).[12] XMP is then subsequently converted
to guanosine 5′-monophosphate (GMP) by a GMP synthetase.[13]
Figure 1
Purine nucleotide biosynthesis. The commonly occurring
guanine
nucleotide biosynthetic and salvage reactions are shown, as is the
adenine nucleotide biosynthetic pathway. The IMPDH reaction is boxed.
NK, nucleoside kinase; HPRT, hypoxanthine phosphoribosyl transferase;
XPRT, xanthine phosphoribosyl transferase; GPRT, guanine phosphoribosyl
transferase; GMPS, guanosine 5′-monophosphate synthetase; GMPR,
guanosine 5′-monophosphate reductase; ADSS, adenylosuccinate
synthetase; ADSL, adenylosuccinate lyase.
Purine nucleotide biosynthesis. The commonly occurring
guanine
nucleotide biosynthetic and salvage reactions are shown, as is the
adenine nucleotide biosynthetic pathway. The IMPDH reaction is boxed.
NK, nucleoside kinase; HPRT, hypoxanthine phosphoribosyl transferase;
XPRT, xanthine phosphoribosyl transferase; GPRT, guanine phosphoribosyl
transferase; GMPS, guanosine 5′-monophosphate synthetase; GMPR,
guanosine 5′-monophosphate reductase; ADSS, adenylosuccinate
synthetase; ADSL, adenylosuccinate lyase.IMPDH has been deemed essential in every pathogen analyzed
to date,
including Mtb, Staphylococcus aureus, and Pseudomonas aeruginosa, which are three of
the most serious bacterial threats.[14−16] However, this has been
somewhat contradictory,[17] in comparing
cell versus animal work. IMPDH is an ubiquitous enzyme present in
several eukaryotes, bacteria, and protozoa.[18] The IMPDH reaction involves two chemical transformations. The first
step of the IMPDH catalyzed reaction involves the attack of catalytic
Cys on substrate IMP followed by hydride transfer to NAD+, forming the covalent enzyme intermediate E-XMP*. In the second
step, E-XMP* is hydrolyzed to XMP.[11] The
enzyme exists in two different conformations, an open form that accommodates
both the substrate and cofactor during the first step and a closed
form where the active site flap moves into the NAD+-binding
site for the E-XMP* hydrolysis.[19]In recent years, there has been considerable effort aimed at identifying
small molecule inhibitors of IMPDH as potential antitubercular agents.
X-ray crystal structures of a truncated form of the Mtb enzyme in complex with some of these compounds have been reported.[20−26]In antibacterial drug discovery, and especially in TB drug
discovery,
high-throughput screening (HTS) typically identifies a number of leads
that show high potency in vitro, but most did not
show any translation to an in vivo effect. It is
also inevitable that the HTS libraries represent only a small fraction
of possible chemical space and so limit confidence in finding a good
starting point for subsequent development. Phenotypic screens can
potentially lead to the identification of a molecule that modifies
a disease phenotype by acting on a previously undescribed target or
by acting simultaneously on more than one target.[27] However, for many of these hits the relevant target or
targets has not yet been identified, thus preventing further target-based
optimization of the compounds.[28,29]The previously
reported IMPDH inhibitors, such as compounds 7759844
(1), MAD1, P41, VCC234718, and DDD00079282 (Figure ), were identified
by phenotypic screening or target based HTS of compound libraries.[21,23−25]
Figure 2
Structures of previously reported IMPDH inhibitors (α
Ki values against IMP). All compounds showed uncompetitive
inhibition
with respect to IMP.
Structures of previously reported IMPDH inhibitors (α
Ki values against IMP). All compounds showed uncompetitive
inhibition
with respect to IMP.We have sought to develop IMPDH inhibitors using a fragment-based
approach. Fragment-based drug discovery (FBDD) is now established
in both industry and academia as an alternative approach to high-throughput
screening for the generation of hits or chemical tools for drug targets.[30]We have previously reported the discovery
of several series of
novel and potent inhibitors using FBDD to target Mtb. Previously we have reported the fragment elaboration strategies
that we have applied which have included fragment growing, merging,
and linking. Although fragment linking is conceptually the most appealing
strategy for fragment elaboration, in practice, this strategy can
be challenging where the choice of the optimal fragment linker can
be crucial.[31,32]In the fragment-based approach,
biophysical techniques are usually
used to identify small chemical compounds (fragments) that bind with
low affinity to the drug target. X-ray crystallography is then usually
employed to establish the binding mode of the fragment and to facilitate
the design of elaborated fragments. The availability of high-resolution
X-ray crystal structures of a truncated form of the IMPDH,[21,24,25] in both the substrate-free and
substrate/ligand-bound forms, makes this enzyme attractive for a fragment-based
approach.In this work, the discovery of a new class of potent
nM inhibitors
of IMPDH using a FBDD approach is reported. A library of 960 fragments
was screened against Mth IMPDH ΔCBS using a
biochemical assay. The fragment hits from this assay were examined
using X-ray crystallography, and an X-ray crystal structure of one
of the fragment complexes was solved to a resolution of 1.45 Å.
Examination of the X-ray crystal structure suggested a strategy of
fragment-linking for optimization of this fragment hit.
Results and Discussion
Fragment
Screening
An in-house fragment library composed
of 960 fragments was screened using a biochemical assay against Mth IMPDH ΔCBS. Mth IMPDH, which
shares 85% sequence identity with Mtb IMPDH and is
100% identical in the active site,[24,25] was chosen
for the fragment screening and structural studies because it gave
higher protein expression yields and better diffracting crystals than
the Mtb orthologue. IMPDH activity was monitored
spectrophotometrically by measuring the formation of NADH at 340 nm.
The biochemical assay was performed at a fragment concentration of
1 mM, and hits were retested in triplicate. Compound 1 (7759844) previously reported as IMPDH inhibitor was used as a positive
control in assays (Table ).[23]
Table 1
Structures
and Activities for Compound 1 (7759844) and the Most
Potent Fragment Hits Found in the
Screen against Mth IMPDH ΔCBS
% Inhibition
at 10 μM.
% Inhibition
at 10 μM.The screen
resulted in 18 hits (1.9% hit rate), where a hit was
defined as a compound that gave greater than 50% inhibition at a concentration
of 1 mM. A complete list of fragment hits identified is included in Table . A number of common
scaffolds were observed, in particular phenylimidazole (fragments 2–4), aminopyrazole (fragments 5–6), and naphthol (fragments 7–11), and the remaining compounds contained a substituted phenyl
or a heterocyclic five membered ring (fragments 12–19).The IC50 values of six of the most active
fragments
were measured and the IC50 and ligand efficiency (LE) values
of these fragments are summarized in Table . The fragment screen provided an array of
hits with IC50 ranging from 325 μM to 674 μM
and ligand efficiencies from 0.31 to 0.42. The inhibition constant
[Ki] with respect to both substrates IMP
and NAD+ was determined by assaying various concentrations
of each inhibitor with five different concentrations of substrate
and a fixed saturating concentration of the cosubstrate. The inhibition
data for these fragments are summarized in Table . All compounds yielded an uncompetitive
inhibition pattern with respect to NAD+ with Ki values ranging from 262 to 525 μM. Fragments 14, 16, and 19 yielded a mixed inhibition
with respect to IMP with Ki values ranging
from 126 to 398 μM, and compounds 1, 2, 11, and 18 yielded an uncompetitive inhibition
with Ki values ranging from 361 to 609
μM.
Table 2
IC50, Ligand Efficiency,
and Alpha Ki Values of Fragment Hits and
Compound 1 (7759844) against Mth IMPDH
ΔCBS
Compd
IC50 (μM) (LE)a
IMP αKi (μM)
NAD αKi (μM)
1
0.77 ± 0.06 (0.40)
0.86 ± 0.03 (UC)b
0.55 ± 0.02 (UC)
2
674 ± 53 (0.36)
609 ± 3 (UC)
512 ± 23 (UC)
11
336 ± 6 (0.39)
361 ± 53 (UC)
310 ± 30 (UC)
14
400 ± 7 (0.42)
238 ± 6 (Mixed)
262 ± 39 (UC)
16
433 ± 60 (0.31)
398 ± 81 (Mixed)
525 ± 64 (UC)
18
512 ± 37 (0.37)
554 ± 70 (UC)
452 ± 26 (UC)
19
325 ± 9 (0.37)
126 ± 34 (Mixed)
355 ± 47 (UC)
Ligand efficiency
was calculated
using the equation LE = (1.37 × pIC50)/HA, where HA
means heavy atom, i.e., a non-hydrogen atom.
UC: Uncompetitive inhibition.
Ligand efficiency
was calculated
using the equation LE = (1.37 × pIC50)/HA, where HA
means heavy atom, i.e., a non-hydrogen atom.UC: Uncompetitive inhibition.Inhibition constants of compound 2 toward full-length Mtb IMPDH were also
determined. Compound 2 inhibited
full-length Mtb IMPDH enzyme with a Ki value of 572 ± 14 μM with IMP as the substrate
and a Ki value of 534 ± 18 μM
with NAD as the substrate, which are similar to the Ki values observed for the Mth IMPDH ΔCBS
enzyme.
X-ray Structure of Compounds 1 (7759844) and 2
Compound 1 and the six fragment hits
(2, 11, 14, 16, 18, and 19) were selected for structural
characterization using X-ray crystallography by soaking into preformed
crystals of Mth IMPDH ΔCBS as previously described.[20] After molecular replacement, clear electron
density was observed in the 2F0 – Fc difference map (σ = 3.0) for IMP, and
in addition density was observed for one molecule of compound 1 (Figure A) and two molecules of compound 2 (Figure B), which were partially occupying
the NAD+ binding site. None of the other fragments (11, 14, 16, 18, and 19) showed any electron density in the X-ray crystal structures.
Although the kinetic studies of compounds 1 and 2 suggested that these two compounds are uncompetitive with
respect to NAD+, the binding mode of compounds 1 and 2 closely resembles that of other previously reported
uncompetitive inhibitors of Mtb IMPDH.[21,24,25] It has been proposed that the
uncompetitive mode of inhibition of IMPDH inhibitors with respect
to NAD+ is consistent with their binding preferentially
to the covalent IMPDH-XMP* intermediate after NADH has been released.[24,25,33]
Figure 3
X-ray crystal structure of compounds 1 and 2 bound to Mth IMPDH ΔCBS.
Ligand interactions
are represented as dotted lines; hydrogen bonds are represented in
red, polar interactions in orange, ionic interactions in yellow, and
aromatic and π interactions in green dotted lines. Protein–ligand
interactions were analyzed using Arpeggio. (a) Interactions made by 1 (green) in the X-ray crystal structure of the complex of
IMPDH with IMP (orange). (b) Interactions made by 2 (yellow)
in the X-ray crystal structure of the complex of IMPDH with IMP. (c)
Structural alignment of the IMPDH crystal structures of 1 (green) and 2 (yellow).
X-ray crystal structure of compounds 1 and 2 bound to Mth IMPDH ΔCBS.
Ligand interactions
are represented as dotted lines; hydrogen bonds are represented in
red, polar interactions in orange, ionic interactions in yellow, and
aromatic and π interactions in green dotted lines. Protein–ligand
interactions were analyzed using Arpeggio. (a) Interactions made by 1 (green) in the X-ray crystal structure of the complex of
IMPDH with IMP (orange). (b) Interactions made by 2 (yellow)
in the X-ray crystal structure of the complex of IMPDH with IMP. (c)
Structural alignment of the IMPDH crystal structures of 1 (green) and 2 (yellow).The structure of Mth IMPDH ΔCBS with
compound 1 showed that the inhibitor binds in the NAD
pocket in a near
identical manner to our recently described IMPDH inhibitors.[24,25] Compound 1 formed strong π interactions with
the hypoxanthine group of IMP, P45′, Y471′, and polar
interactions with G409, E442, P45′, and G470′. The electron
density revealed two molecules of compound 2 within the
NAD binding pocket of IMPDH. One molecule of compound 2 stacked with IMP, forming extensive π interactions with the
hypoxanthine group of IMP. This fragment was further stabilized through
polar interactions, hydrogen bonds, and π interactions to surrounding
residues in the active site pocket, including A269, G318, and E442.
The other molecule of compound 2 sits closer to the opening
of the active site, making polar interactions with N273 and E442,
and π interactions with H270 and Y471′.A comparison
of the structure of Mth IMPDH ΔCBS
with compound 1 with the fragment 2 structure
shows that the two molecules of 2 mimic the position
of the larger inhibitor 1 (Figure C).
Fragment Elaboration
Fragment 2 was selected
as the starting point for exploration because of the ease of synthetic
modification and the availability of a X-ray crystal structure to
guide the optimization. For chemical elaboration of 2, fragment linking as well as fragment growing were considered. As
the two molecules of fragment 2 are found to bind in
adjacent regions of the target protein, the fragment-linking approach
was the more attractive option. However, before fragment linking,
fragment 2 was further optimized with the aim of improving
the binding affinity. The structures and inhibitory activities of
these compounds against Mth IMPDH ΔCBS are
summarized in Tables and 4. The corresponding data for fragment 2 have also been included for comparative purposes.
Table 3
Structures and Activities for Fragment 2 and Compounds 20–28 against Mth IMPDH ΔCBS
Table 4
Structures and Activities for Fragment 2 and Compounds 29–35 against Mth IMPDH ΔCBS
Ligand efficiency
was calculated
using the equation LE = (1.37 × pIC50)/HA, where HA
means heavy atom, i.e., a non-hydrogen atom.
ND: not determined. IC50 values were
determined for compounds that showed >50% inhibition. %Inhibition at 10 μM.
Ligand efficiency
was calculated
using the equation LE = (1.37 × pIC50)/HA, where HA
means heavy atom, i.e., a non-hydrogen atom.ND: not determined. IC50 values were
determined for compounds that showed >50% inhibition. %Inhibition at 10 μM.All compounds were evaluated at a concentration of
100 μM
with Mth IMPDH ΔCBS.
Fragment Growing
The fragment-growing strategy involves
using structure-based drug design to form additional interactions
by growing out from the starting fragment. Fragment 2 was modified at the 2-position of the imidazole ring to explore
the introduction of various aromatic rings linked by a thioacetamide
(20–22) to form π interactions
with the hypoxanthine group of IMP (Figure S1). Such modifications gave compounds with improved Mth IMPDH ΔCBS inhibition. The phenyl and benzofuran derivatives
(20 and 21) showed 13 and 31% inhibition,
respectively. Mth IMPDH ΔCBS inhibition was
shown to be sensitive to minor modifications of the phenyl substituent
groups; for example, the 4-iodo substituted 22 showed
30% inhibition at 100 μM, whereas the nonsubstituted compound, 20, showed 13% inhibition at the same concentration. The effect
of the removal of the 4-bromo group was investigated and compounds 23–25 were synthesized. Removal of the
bromo substituent in compounds 20–22 (13–31% inhibition at 100 μM) was tolerated (23–25, 3–46% inhibition at 100
μM). The importance of the aromatic amide linked by a thioacetamide
was subsequently examined. Replacing the phenyl with isopropyl (26) resulted in complete loss of activity. Substitution of
the thioacetamide by a thioacetic acid also led to a complete loss
in activity (compounds 27 and 28).
Fragment
Linking and SAR
Examination of the X-ray crystal
structure of the previously reported inhibitor 1 when
overlaid with fragment 2 revealed that the distance between
the 4-position of the phenyl ring of fragment 2 and the
2-position of the imidazole ring represents the closest approach of
the molecules (Figure S1). On the basis
of this structural information, three different linkers were designed
to connect the two copies of the fragment 2 at these
positions (compounds 29–31). Initially,
a thioacetamide and urea linker moieties were examined. Compounds 29 and 30 showed 20% and 56% of Mth IMPDH ΔCBS inhibition, respectively, at 100 μM (Table ). Interestingly,
compound 30 showed a 12-fold improvement in Mth IMPDH ΔCBS inhibitory activity with an IC50 of
58 μM, compared to the fragment 2.The lactate
linker was then used, but all attempts to link 4-phenylimidazole with
4-(4-bromophenyl)-1H-imidazole were unsuccessful.
We therefore decided to synthesize compound 31, which
contains 1-methyl-4-phenyl-1H-imidazole and 4-(4-bromophenyl)-1H-imidazole linked with a lactate analogue, as in compound 1. Compound 31 (Table ) showed markedly improved Mth IMPDH ΔCBS inhibition with a LE of 0.29 and IC50 of 0.52 μM, which is 1300-fold more potent compared to the
fragment 2. Although compound 31 binds in
the cofactor site, the mechanism of inhibition can vary depending
on its relative affinities for the E·IMP and E-XMP* complexes.[8,11] Kinetic evaluation of compound 31 showed the mode of Mth IMPDH ΔCBS inhibition was uncompetitive with respect
to both IMP and the NAD+ cofactor (see Figure S2, Supporting Information) with a Ki value of 0.30 ± 0.02 μM with IMP as the substrate
and a Ki value of 0.20 ± 0.01 with
NAD as the substrate.Inhibition constants of compound 31 toward full-length Mtb IMPDH were also
determined. Compound 31 inhibited full-length Mtb IMPDH enzyme with a Ki value
of 0.61 ± 0.05 μM with IMP
as the substrate and a Ki value of 0.39
± 0.02 μM with NAD as the substrate. The inhibition constants
were consistent with the data using the Mth IMPDH
ΔCBS enzyme.Removal of the bromo substituent in compound 29 to
give compound 32 was well tolerated (Table ). The importance of the imidazole
group for the inhibitory activity against Mth IMPDH
ΔCBS was confirmed by replacing of the 4-(4-bromophenyl)-1H-imidazole substituent of compound 31 with
4-(4-bromophenyl)oxazole (33) which resulted in a 4-fold
loss of activity (Table ). Replacing the 4-(4-bromophenyl)-1H-imidazole
of 31 by a phenyl (34) resulted in complete
loss of activity (Table ). However, the 4-iodophenyl derivative 35 demonstrated
slightly improved activity (IC50 = 0.47 μM, LE =
0.34) compared to compounds 31 (IC50 = 0.52
μM, LE = 0.29) and 1 (IC50 = 0.77 μM,
LE = 0.40). It is noteworthy that LE of compound 35 was
comparable to that of the original fragment hit 2 (LE
= 0.36) and other reported IMPDH inhibitors.[34]The importance of the 4-bromo substituent on the phenyl ring
in
compound 31 was also explored (compounds 36–39, Table ). Removal of the bromo substituent (36) resulted in a 5-fold loss in activity, whereas replacing this group
with an iodine (39) or a morpholine ring (37) resulted in less than 1.3-fold loss in activity. The electronic
nature of the substituent in this position had little effect on inhibitory
activity. For example, an electron-donating methoxy (38) retained activity comparable to that of the bromine derivative
(31). Among them, imidazoles 31, 37–38, and 39 were found to be potent
inhibitors of Mth IMPDH ΔCBS with IC50 values ranging from 520 to 690 nM.
Table 5
Structures
and Activities for Compounds 31, 36–39, and (S)-31 against Mth IMPDH ΔCBS
Ligand efficiency was calculated
using the equation LE = (1.37 × pIC50)/HA, where HA
means heavy atom, i.e., a non-hydrogen atom.
Ligand efficiency was calculated
using the equation LE = (1.37 × pIC50)/HA, where HA
means heavy atom, i.e., a non-hydrogen atom.The (S)-isomer of 31 was found to
bind preferentially (Table ), with the racemate 31 having approximately
half the Mth IMPDH ΔCBS inhibition of the (S)-isomer 31. This observation accords with
the results previously reported for other series of IMPDH inhibitors.[21,35]
X-ray Structure of Compound 31
From crystals
soaked with compound 31, the 2F0 – Fc difference map (σ
= 3.0) revealed strong density for the inhibitor. The structure of
compound 31 (Figure A) showed that it bound in a nearly identical manner
to compounds 1 and 2 in the NAD binding
pocking (Figure B),
stacking with IMP, and maintaining the interactions with H270, N273,
E442, and P45 and Y471 from the neighboring subunit. Compound 31 made additional interactions in the binding pocket, including
polar interactions with D267 and N297.
Figure 4
X-ray crystal structure
of compound 31 bound to Mth IMPDH ΔCBS.
Protein–ligand interactions
were analyzed using Arpeggio. (a) Interactions made by 31 in the X-ray crystal structure of the complex of IMPDH with IMP
(orange). Ligand interactions are represented as dotted lines; hydrogen
bonds are represented in red, polar interactions in orange, ionic
interactions in yellow, and aromatic and π interactions in green
dotted lines. (b) Structural alignment of the IMPDH crystal structures
of 1 (lilac), 2 (yellow), and 31 (green).
X-ray crystal structure
of compound 31 bound to Mth IMPDH ΔCBS.
Protein–ligand interactions
were analyzed using Arpeggio. (a) Interactions made by 31 in the X-ray crystal structure of the complex of IMPDH with IMP
(orange). Ligand interactions are represented as dotted lines; hydrogen
bonds are represented in red, polar interactions in orange, ionic
interactions in yellow, and aromatic and π interactions in green
dotted lines. (b) Structural alignment of the IMPDH crystal structures
of 1 (lilac), 2 (yellow), and 31 (green).
Whole-Cell Activity against Mtb
The
whole-cell activity of the most potent analogues 31, 33–39, and (S)-31in vitro was determined against Mtb H37Rv (see Table S1, Supporting
Information). No significant inhibition of bacterial growth was detected
for any of the compounds (MIC90 ≥ 50 μM) over
the tested concentration range (0–100 μM). There are
currently ongoing efforts to explain the lack of efficacy of the potent Mth IMPDH ΔCBS inhibitors described which could be
caused by low membrane permeability, poor metabolic stability, and/or
drug-efflux mechanisms.
Synthetic Chemistry
Compounds 20–29 and 32 were synthesized
from 4-phenyl-1H-imidazole-2-thiol or 4-(4-bromophenyl)-1H-imidazole-2-thiol according to the sequence described
in Scheme . Thioacetic
derivatives 27 and 28 were prepared by treatment
of 4-phenyl-1H-imidazole-2-thiol or 4-(4-bromophenyl)-1H-imidazole-2-thiol with 2-chloroacetic acid in the presence
of NaOH
followed by neutralization with hydrochloric acid.
Scheme 1
Synthesis of Thioacetamide
and Thioacetic Derivatives 20–29 and 32
Synthesis of Thioacetamide
and Thioacetic Derivatives 20–29 and 32
Reagent and conditions: (a)
2-chloroacetic acid, NaOH, EtOH, 80 °C, 69–72%; (b) isopropylamine,
HATU, DIPEA, EtOAc, rt, 73%; (c) chloroacetyl chloride, TEA, DCM,
rt, 89–99%; (d) NaOH, MeOH, H2O, 70 °C, 69–96%.Thioacetamide 26 was synthesized
by amide coupling
between thioacetic derivative 28 and isopropylamine.
Similarly, thioacetamides 20–25, 29, and 32 were prepared by cross-linking 4-phenyl-1H-imidazole-2-thiol or 4-(4-bromophenyl)-1H-imidazole-2-thiol with α-chloroacetamides 40–43, which were obtained by acylation of anilines with various
substituents with chloroacetyl chloride.The synthesis of urea 30 was achieved as shown in Scheme by coupling amine 44 and N-(4-(1H-imidazol-4-yl)phenyl)-1H-imidazole-1-carboxamide, which was obtained from commercially
available 4-(1H-imidazol-4-yl)aniline as a crude
intermediated, in the presence of N,N-diisopropylethylamine.
Amine 44 was synthesized from benzyl ((5-(4-bromophenyl)-1H-imidazol-2-yl)methyl)carbamate,[36] followed by the deprotection of the benzyloxycarbonyl group under
acidic condition.
Scheme 2
Synthesis of Urea 30
Reagent and conditions: (a)
CDI, THF, rt, quant.; (b) HCl (4 M), dioxane, 100 °C, quant.;
(c) DIPEA, DMF, rt, 26%.
Synthesis of Urea 30
Reagent and conditions: (a)
CDI, THF, rt, quant.; (b) HCl (4 M), dioxane, 100 °C, quant.;
(c) DIPEA, DMF, rt, 26%.Imidazole derivatives 31, 36–38, and 39 were prepared following the synthetic
procedure outlined in Scheme . 2-Aminoimidazoles 50–54 were synthesized according to a published microwave-assisted protocol.[37] In brief, 2-aminoimidazoles 50–54 were prepared by reaction of the commercially available
α-haloketones and N-acetylguanidine, followed
by deacetylation (Scheme ). Acid 59 was synthesized starting with imidazole 55, which was prepared by reaction of 2-bromo-4′-hydroxyacetophenone
with formamide as reported previously.[38] The phenol 57 was synthesized by alkylation of imidazole 55, followed by deprotection of the methyl ether with BBr3. Substituted phenol 57 was converted to the
ether 58 upon treatment with methyl 2-bromopropionate
in the presence of Cs2CO3. Enantiomerically
pure phenyl ethers were synthesized by using Mitsunobu reaction conditions
with ethyl d-lactate (Scheme ). After the hydrolysis of the ester group, the resulting
carboxylic acid 59 was treated with thionyl chloride
to give the acid chloride 60, which was reacted with
2-aminoimidazoles 50–54 to afford
imidazole derivatives 31, 36–38, and 39.
Scheme 3
Synthesis of 2-Acylaminoimidazole
Derivatives 31, 36–38, and 39
Reagent and conditions: (a)
acetylguanidine, CH3CN, μW, 100 °C, 71–88%;
(b) (1) 20% H2SO4, MeOH:H2O (1:1
v/v), μW, 100 °C; (2) aq Na2CO3,
79–95%; (c) formamide, 150 °C, 99%; (d) CH3I, Cs2CO3, DMF, rt, 88%; (e) BBr3 (1 M in DCM), −78 °C to rt, DCM, 96%; (f) Cs2CO3, methyl 2-bromopropionate, DMF, 60 °C, 92%; (g)
NaOH, THF:H2O (2:1 v/v), 80 °C, 65%; (h) SOCl2, 80 °C, quant.; (i) TEA, DCM, 40 °C, 61–72%.
Scheme 4
Synthesis of the S-Enantiomer of 31
Reagent and conditions: (a)
(1) Ethyl d-lactate, PPh3, DEAD, THF, rt; (2)
NaOH, THF:H2O (2:1 v/v), rt, 56%; (b) SOCl2,
80 °C, quant.; (c) 50, TEA, DCM, 40 °C, 69%.
Synthesis of 2-Acylaminoimidazole
Derivatives 31, 36–38, and 39
Reagent and conditions: (a)
acetylguanidine, CH3CN, μW, 100 °C, 71–88%;
(b) (1) 20% H2SO4, MeOH:H2O (1:1
v/v), μW, 100 °C; (2) aq Na2CO3,
79–95%; (c) formamide, 150 °C, 99%; (d) CH3I, Cs2CO3, DMF, rt, 88%; (e) BBr3 (1 M in DCM), −78 °C to rt, DCM, 96%; (f) Cs2CO3, methyl 2-bromopropionate, DMF, 60 °C, 92%; (g)
NaOH, THF:H2O (2:1 v/v), 80 °C, 65%; (h) SOCl2, 80 °C, quant.; (i) TEA, DCM, 40 °C, 61–72%.
Synthesis of the S-Enantiomer of 31
Reagent and conditions: (a)
(1) Ethyl d-lactate, PPh3, DEAD, THF, rt; (2)
NaOH, THF:H2O (2:1 v/v), rt, 56%; (b) SOCl2,
80 °C, quant.; (c) 50, TEA, DCM, 40 °C, 69%.The syntheses of 2-acylaminooxazole 33 and amides 34–35 were achieved
as shown in Scheme . 2-Acylaminooxazole
derivative 33 was obtained by coupling the acid chloride
derivative 60 with 2-aminooxazole derivative 61, which was prepared by reaction of 2,4′-dibromoacetophenone
with urea. Similarly, amides 34–35 were prepared by coupling the corresponding anilines with the acid
chloride derivative 60.
Scheme 5
Synthesis of 2-Acylaminooxazole
Derivative 33 and Amides 34–35
Reagent and conditions: (a)
urea, CH3CN, 80 °C, 78%; (b) TEA, DCM, rt (for 34 and 35) or 40 °C (for 33), 64–76%.
Synthesis of 2-Acylaminooxazole
Derivative 33 and Amides 34–35
Reagent and conditions: (a)
urea, CH3CN, 80 °C, 78%; (b) TEA, DCM, rt (for 34 and 35) or 40 °C (for 33), 64–76%.
Conclusions
FBDD
has emerged as a robust approach to identify small molecules
that bind to a wide range of therapeutic targets. Fragment elaboration
strategies have resulted in the development of a number of compounds
that have progressed into clinical trials. Within the area of TB drug
discovery, a number of HTS and phenotypic screens have been performed
during the past decade. Although HTS identified a number of leads
that show high potency in vitro, the translation
to an in vivo effect has proven challenging.This study illustrates the successful application of a fragment-based
approach followed by fragment optimization to obtain nanomolar affinity
ligands of IMPDH. A library of 960 fragments were screened against Mth IMPDH ΔCBS, and from the screen the phenylimidazole
fragment hit 2 (IC50 = 674 μM) was identified.
Kinetic experiments showed that 2 was an uncompetitive
inhibitor of Mth IMPDH ΔCBS with respect to
NAD+ and IMP. Two molecules of the fragment 2 were shown to bind at the NAD binding site of the enzyme. The X-ray
crystal structure also revealed that one molecule of fragment 2 makes π interactions with IMP and the other molecule
sits closer to the opening of the active site, making polar interactions
with N273 and E442 and π interactions with H270 and Y471′.
This provides potential for further optimization of fragment 2. To explore better the possibilities given by fragment 2, fragment-linking and fragment-growing strategies were employed,
resulting in low micromolar to nanomolar affinity compounds. Among
them, compounds 31, 35, and 37– 39 were the most potent IMDPH inhibitors of
the series described in this work with IC50 values between
0.47 and 0.69 μM, which represent >1000-fold improvement
in Mth IMPDH ΔCBS potency over the initial
fragment hit.
Compound 31 was shown to bind at the NAD binding site
of the enzyme, and the X-ray crystal structure also revealed that
it makes π interactions with IMP, maintaining the interactions
with H270, N273, E442, and P45 and Y471 from the neighboring subunit.
Moreover, compound 31 made additional interactions in
the binding pocket, including polar interactions with D267 and N297.
A comparison of this structure with the fragment 2 structure
shows that the two molecules of 2 mimic the position
of the larger inhibitor 31. This is the first example
of utilizing the fragment-based approach specifically to identify
new potent inhibitors of IMPDH. Further structural optimization to
improve the cellular activity of these analogues is ongoing with the
aim of developing novel classes of anti-TB agents.
Experimental Section
Chemistry
General Experimental Methods
Solvents
were distilled prior to use and dried by standard methods. Unless
otherwise stated, 1H and 13C NMR spectra were
obtained in CDCl3, MeOD, or DMSO-d6 solutions using either a Bruker 400 MHz AVANCE III HD Smart
Probe, 400 MHz QNP cryoprobe, or 500 MHz DCH cryoprobe spectrometer.
Chemical shifts (δ) are given in ppm relative to the residual
solvent peak (CDCl3: 1H, δ = 7.26 ppm; 13C, δ = 77.16 ppm), and the coupling constants (J) are reported in hertz (Hz). Optical rotations were measured
on a PerkinElmer Polarimeter 343 at 589 nm (Na D-line), and specific
rotations are reported in 10–1 deg cm2 g–1. Microwave reactions were performed using
a Biotage Initiator system under reaction conditions as indicated
for each individual reaction.Reactions were monitored by TLC
and LCMS to determine consumption of starting materials. Flash column
chromatography was performed using an Isolera Spektra One/Four purification
system and the appropriately sized Biotage SNAP column containing
KP-silica gel (50 μm). Solvents are reported as volume/volume
eluent mixture where applicable.High resolution mass spectra
(HRMS) were recorded using a Waters
LCT Premier Time of Flight (TOF) mass spectrometer or a Micromass
Quadrapole-Time of Flight (Q-TOF) spectrometer.Liquid chromatography
mass spectrometry (LCMS) was carried out
using an Ultra Performance Liquid Chromatographic system (UPLC) Waters
Acquity H-class coupled to a Waters SQ Mass Spectrometer detector.
Samples were detected using a Waters Acquity TUV detector at 2 wavelengths
(254 and 280 nm). Samples were run using an Acquity UPLC HSS column
and a flow rate of 0.8 mL/min. The eluent consisted of 0.1% formic
acid in water (A) and acetonitrile (B); gradient, from 95% A to 5%
A over a period of 4 or 7 min.All final compounds had a purity
greater than 95% by LCMS analysis.
General Method A: Synthesis
of Thioacetamides 20–25, 29, and 32
To a solution of the 2-chloroacetamide
derivative (0.17 mmol) in
MeOH (8 mL) was added 4-(4-bromophenyl)-1H-imidazole-2-thiol
or 4-phenyl-1H-imidazole-2-thiol (0.17 mmol), followed
by a solution of NaOH (0.67 mmol) in H2O (2.5 mL). The
reaction mixture was stirred at 70 °C for 3 h. After cooling
to rt, the solvents were removed in vacuo and the
resulting residue was taken up in 30 mL of EtOAc and washed with H2O (15 mL). The aqueous layer was extracted with EtOAc (3 ×
30 mL) and the combined organic layers were washed with brine and
dried over MgSO4. Filtration and evaporation afforded crude
products, which were purified as indicated below.
General Method
B: Synthesis of Thioacetic Derivatives 27 and 28
A solution of 4-phenyl-1H-imidazole-2-thiol
or 4-(4-bromophenyl)-1H-imidazole-2-thiol (2.80 mmol)
and NaOH (5.0 mmol) in EtOH (5 mL)
was reflux for 1 h. After cooling to rt, a solution of 2-chloroacetic
acid (2.80 mmol) in EtOH (2 mL) was added. The reaction was stirred
at reflux for an additional 3 h and then cooled to 0 °C. The
reaction mixture was diluted with cold water (5 mL) and acidified
with 1 M HCl. The precipitated product was collected by filtration
and washed with DCM (2 × 2 mL).
General Method C: Synthesis
of Compounds 31, 33, and 36–39
A
mixture of acid 59 (0.20 mmol) and SOCl2 (2
mL) was heated at 80 °C for 2 h. The solvent was removed under
reduced pressure to give the acid chloride 60 as a white
solid. The resulting solid was immediately dissolved in anhydrous
DCM, and the resulting solution was added slowly dropwise at 0 °C
to a solution of the corresponding substituted 2-aminoimidazoles (0.20
mmol) and TEA (0.80 mmol) in anhydrous DCM (5 mL). The reaction mixture
was stirred at 40 °C for 36 h and then diluted with DCM (20 mL)
and washed with saturated aqueous NaHCO3. The aqueous phase
layer was then extracted with DCM (2 × 20 mL), and the combined
organic layers were dried over MgSO4 and filtered, and
the solvent was removed under reduced pressure to afford a yellow
oil, which was purified by flash chromatography, eluting with the
solvent system specified.
General Method D: Synthesis of Amides 34–35
A mixture of acid 59 (0.20 mmol)
and SOCl2 (2 mL) was heated at 80 °C for 2 h. The
solvent was removed under reduced pressure to give 60 as a white solid. The resulting solid was immediately dissolved
in anhydrous DCM, and the resulting solution was added slowly dropwise
at 0 °C to a solution of the corresponding aniline (0.20 mmol)
and triethylamine (0.80 mmol) in anhydrous DCM (5 mL). The reaction
mixture was stirred at rt for 4 h and then diluted with DCM (20 mL)
and washed with saturated aqueous NaHCO3. The aqueous phase
layer was then extracted with DCM (2 × 20 mL), and the combined
organic layers were dried over MgSO4 and filtered, and
the solvent was removed under reduced pressure to afford a yellow
oil, which was purified by flash chromatography, eluting with the
solvent system specified.
General Method E: Synthesis of α-Chloroacetamides 40–43
Et3N (4.77 mmol)
followed by a solution of chloroacetyl chloride (4.77 mmol) in DCM
(3 mL) were added to a stirred solution of the corresponding aniline
(4.38 mmol) in DCM (5 mL) at rt. The reaction mixture was stirred
at rt for 2–4 h. The reaction was then diluted with DCM (20
mL) and washed with saturated aqueous NaHCO3, 1 M HCl,
and brine and dried over anhydrous MgSO4, and the solvent
was removed under reduced pressure. Compound 43 was purified
by flash chromatography eluting with the solvent system specified,
although other analogues were used in subsequent reactions without
further purification.
General Method F: Synthesis of Substituted N-(1H-Imidazol-2-yl)acetamides 45–49
A mixture of the corresponding 2-bromoacetophenone
derivative (0.38 mmol) and acetylguanidine (1.13 mmol) in anhydrous
acetonitrile (3 mL) was heated at 100 °C using microwave irradiation
for 15 min. The solvent was removed, and the residue was taken in
H2O (3 mL) and filtered, and the solid was washed with
H2O (2 mL × 2) and DCM (2 mL). The solid obtained
was used in the next step without further purification.
General Method
G: Synthesis of Substituted 2-Aminoimidazoles 50–54
To a solution of the corresponding
substituted N-(1H-imidazol-2-yl)acetamides
(0.31 mmol) in a 1:1 v/v mixture of MeOH and H2O (2.4 mL)
was added concentrated H2SO4 (0.6 mL), and the
reaction mixture was heated at 100 °C under microwave irradiation
for 15–30 min. The reaction mixture was concentrated, and the
resulting residue was resuspended in H2O (5 mL), and a
saturated aqueous Na2CO3 was added until pH
8. The product was extracted into EtOAc (3 × 40 mL). The combined
organic fractions were dried over MgSO4, and the solvent
was removed under reduced pressure. The resulting solid was used in
the next reaction without further purification.
To
a solution of 4-(1H-imidazol-4-yl)aniline (50 mg,
0.31 mmol) in THF (5 mL) was added 1,1-carbonyldiimidazole (76 mg,
0.47 mmol). The mixture was stirred at rt overnight, and then it was
filtered. The filter was washed with THF (2 × 3 mL) and the resulting
solid (30 mg, 0.12 mmol) was dissolved in DMF (3 mL), and 44 (42 mg, 0.12 mmol) and N,N-diisopropylethylamine
(74 μL, 0.48 mmol) were added. The reaction mixture was stirred
at rt for 14 h, and then it was poured into water, extracted with
EtOAc, dried over Na2SO4, and concentrated.
After cooling to 0 °C, a 1:5 v/v mixture of MeOH and DCM (12
mL) was added and the suspended solid was collected by filtration
and dried at vacuum to yield 30 (35 mg, 0.08 mmol, 26%
yield) as a white solid. 1H NMR (500 MHz, DMSO-d6) δ 12.03 (br s, 2H), 8.64 (br s, 1H),
7.73–7.68 (m, 2H), 7.64–7.59 (m, 2H), 7.57 (d, J = 1.9 Hz, 1H), 7.52–7.47 (m, 2H), 7.47–7.39
(m, 2H), 7.39–7.32 (m, 2H), 6.58 (br s, 1H), 4.33 (d, J = 5.5 Hz, 2H) ppm. 13C NMR (125 MHz, DMSO-d6) δ 155.5, 146.9, 140.6, 138.9, 138.9,
135.9, 134.6, 131.7, 128.8, 126.6, 125.1, 125.0, 119.0, 118.4, 118.2,
113.7, 111.7, 37.6 ppm. LCMS (ESI−) m/z 437.0 [M – H]−, retention time
5.07 min (97%). HRMS (ESI+): m/z calculated for C20H18BrN6O [M +
H]+: 437.0725. Found: 437.0725.
4-(4-Iodophenyl)-1H-imidazol-2-amine 54 (46 mg, 0.16 mmol) was reacted with the acid chloride 60 (43 mg, 0.16 mmol) and Et3N (90 μL, 0.65
mmol)
according to general method C. The crude product was purified by flash
chromatography (1–20% v/v MeOH in DCM) to give compound 39 as a yellow solid (50 mg, 0.10 mmol, 61% yield). 1H NMR (500 MHz, CDCl3) δ 10.75 (br s, 1H), 9.42
(br s, 1H), 7.69 (t, J = 8.3 Hz, 4H), 7.53–7.33
(m, 3H), 7.10 (d, J = 8.3 Hz, 2H), 6.95 (d, J = 8.9 Hz, 2H), 4.88 (q, J = 6.7 Hz, 1H),
3.71 (s, 3H), 1.66 (d, J = 6.7 Hz, 3H) ppm. 13C NMR (125 MHz, CDCl3) δ 171.5, 155.2, 141.7,
138.0, 137.7, 129.1, 126.5, 126.3, 115.8, 115.3, 108.1, 91.8, 74.5,
33.5, 18.5 ppm. LCMS (ESI+) m/z 514.2
[M + H]+, retention time 4.61 min (95%). HRMS (ESI+): m/z calculated for C22H21IN5O2 [M + H]+: 514.0740.
Found: 514.0740.
2-Chloro-N-phenylacetamide
(40)
Aniline (400 μL, 4.38 mmol) was reacted
with chloroacetyl
chloride (380 μL, 4.77 mmol) and Et3N (670 μL,
4.77 mmol) according to general method E, and used in subsequent reactions
without further purification. Compound 40 was obtained
(720 mg, 4.24 mmol, 96% yield) as a green solid and was used without
further purification. 1H NMR (400 MHz, CDCl3) δ 8.27 (br s, 1H), 7.58 (d, J = 7.7 Hz,
2H), 7.39 (t, J = 7.9 Hz, 2H), 7.20 (t, J = 7.4 Hz, 1H), 4.22 (s, 2H) ppm. 13C NMR (100 MHz, CDCl3) δ 163.9, 136.8, 129.3, 125.4, 120.2, 43.0 ppm. LCMS
(ESI+) m/z 170.1 [M + H]+, retention time 1.58 min (100%). HRMS (ESI+): m/z calculated for C8H9ClNO
[M + H]+: 170.0373. Found: 170.0366. NMR data is in accordance
with literature values.[39]
N-(Benzofuran-5-yl)-2-chloroacetamide (41)
Benzofuran-5-amine (140 mg, 1.05 mmol) was reacted
with chloroacetyl chloride (93 μL, 1.16 mmol) and Et3N (162 μL, 1.16 mmol) according to general method E, and it
was used in subsequent reactions without further purification. Compound 41 was obtained (219 mg, 1.04 mmol, 99% yield) as a green
solid and was used without further purification. 1H NMR
(400 MHz, CDCl3) δ 8.41 (br s, 1H), 7.94 (d, J = 2.2 Hz, 1H), 7.64 (d, J = 2.2 Hz, 1H),
7.47 (dd, J = 8.8, 0.8 Hz, 1H), 7.32 (dd, J = 8.8, 2.2 Hz, 1H), 6.76 (dd, J = 2.2,
0.9 Hz, 1H), 4.23 (s, 2H) ppm. 13C NMR (100 MHz, CDCl3) δ 163.9, 152.4, 146.1, 131.9, 127.9, 117.7, 113.3,
111.6, 106.8, 43.0 ppm. LCMS (ESI+) m/z 210.1 [M + H]+, retention time 1.71 min (100%). HRMS
(ESI+): m/z calculated for C11H9ClNO2 [M + H]+: 210.0322.
Found: 210.0319.
2-Chloro-N-(4-iodophenyl)acetamide
(42)
4-Iodoaniline (500 mg, 2.28 mmol) was reacted
with chloroacetyl
chloride (200 μL, 2.51 mmol) and Et3N (350 μL,
2.51 mmol) according to general method E. The crude product was purified
by flash chromatography (10–100% v/v EtOAc in petroleum ether)
to give compound 42 as a brown solid (600 mg, 2.03 mmol,
89% yield). 1H NMR (400 MHz, CDCl3) δ
8.32–8.11 (br s, 1H), 7.69 (d, J = 8.8 Hz,
2H), 7.36 (d, J = 8.7 Hz, 2H), 4.20 (s, 2H) ppm. 13C NMR (100 MHz, CDCl3) δ 163.8, 138.1, 136.4,
121.8, 88.7, 42.9 ppm. LCMS (ESI−) m/z 293.7 [M – H]−, retention time
2.38 min (100%). HRMS (ESI−): m/z calculated for C8H6ClINO [M – H]−: 293.9183. Found: 293.9187. NMR data is in accordance
with literature values.[40]
4 M HCl (2 mL) was added
slowly to a stirred solution of benzyl ((5-(4-bromophenyl)-1H-imidazol-2-yl)methyl)carbamate[36] (105 mg, 0.26 mmol) in dioxane (1 mL). The mixture was stirred at
100 °C for 4 h, and then the solvents were removed under reduced
pressure to give the desired product as a white solid (86 mg, 0.26
mmol, 100% yield). 1H NMR (500 MHz, MeOD) δ 8.00
(s, 1H), 7.75–7.67 (m, 4H), 4.55 (s, 2H) ppm. 13C NMR (125 MHz, MeOD) δ 140.8, 136.0, 133.7, 128.4, 127.6,
124.7, 117.8, 111.4, 34.8 ppm. LCMS (ESI−) m/z 250.1 [M – H]−, retention
time 2.13 min (100%). HRMS (ESI−): m/z calculated for C10H9BrN3 [M + H]+: 249.9980. Found: 249.9987.
Following general method F, from 2-bromo-1-(4-iodophenyl)ethan-1-one
(214 mg, 0.66 mmol) was obtained 49 (190 mg, 0.58 mmol,
88% yield) as a green solid. 1H NMR (500 MHz, DMSO-d6) δ 11.66 (br s, 1H), 11.21 (br s, 1H),
7.64 (d, J = 8.1 Hz, 2H), 7.51 (d, J = 8.1 Hz, 2H), 7.30 (s, 1H), 2.05 (s, 3H) ppm. 13C NMR
(125 MHz, DMSO-d6) δ 168.5, 141.4,
137.1, 135.1, 134.3, 126.2, 109.9, 90.9, 22.8. LCMS (ESI−) m/z 326.0 [M – H]−, retention time 2.43 min (100%). HRMS (ESI+): m/z calculated for C11H11IN3O [M + H]+: 327.9947. Found: 327.9945.
4-(4-Bromophenyl)-1H-imidazol-2-amine (50)
Following
general method G, from N-(4-(4-bromophenyl)-1H-imidazol-2-yl)acetamide 45 (88 mg, 0.31 mmol)
was obtained 50 (70 mg,
0.29 mmol, 95% yield) as a brown solid.1H NMR (500 MHz,
DMSO-d6) δ 7.52 (d, J = 8.1 Hz, 2H), 7.42 (d, J = 8.6 Hz, 2H), 7.04 (s,
1H), 5.40 (br s, 2H) ppm. 13C NMR (125 MHz, DMSO-d6) δ 150.4, 134.0, 131.2, 126.0, 125.4,
117.5, 110.5 ppm. LCMS (ESI+) m/z 240.1 [M + H]+, retention time 1.50 min (97%). HRMS (ESI+): m/z calculated for C9H9BrN3 [M + H]+: 237.9980. Found: 237.9980. NMR
data is in accordance with literature values.[37]
4-Phenyl-1H-imidazol-2-amine (51)
Following general method G, from N-(4-phenyl-1H-imidazol-2-yl)acetamide 46 (70 mg, 0.35 mmol)
was obtained 51 (47 mg, 0.29 mmol, 84% yield) as a red
solid. 1H NMR (400 MHz, MeOD) δ 7.62–7.46
(m, 2H), 7.31 (t, J = 7.7 Hz, 2H), 7.16 (td, J = 7.5, 1.2 Hz, 1H), 6.91 (s, 1H) ppm.13C NMR
(125 MHz, MeOD) δ 150.3, 132.9, 128.4, 128.1, 125.6, 123.5,
111.2 ppm. LCMS (ESI+) m/z 160.1
[M + H]+, retention time 1.45 min (100%). HRMS (ESI+): m/z calculated for C9H10BrN3 [M + H]+: 160.0875. Found: 160.0873. NMR
data is in accordance with literature values.[37]
4-(4-Morpholinophenyl)-1H-imidazol-2-amine
(52)
Following general method G, from N-(4-(4-morpholinophenyl)-1H-imidazol-2-yl)acetamide 47 (85 mg, 0.30 mmol) was obtained 52 (60 mg,
0.24 mmol, 82% yield) as a brown solid. 1H NMR (400 MHz,
MeOD) δ 7.47 (d, J = 8.8 Hz, 2H), 6.96 (d, J = 8.8 Hz, 2H), 6.78 (s, 1H), 3.91–3.64 (m, 4H),
3.19–3.02 (m, 4H) ppm. 13C NMR (125 MHz, MeOD) δ
149.9, 149.8, 129.3, 125.0, 124.5, 115.8, 113.0, 66.6, 49.4 ppm. LCMS
(ESI+) m/z 245.1 [M + H]+, retention time 1.50 min (100%). HRMS (ESI+): m/z calculated for C13H17N4O [M + H]+: 245.1402. Found: 245.1390.
4-(4-Methoxyphenyl)-1H-imidazol-2-amine (53)
Following
general method G, from N-(4-(4-methoxyphenyl)-1H-imidazol-2-yl)acetamide 48 (83 mg, 0.36 mmol)
was obtained 53 (65 mg,
0.34 mmol, 95% yield) as a brown solid. 1H NMR (400 MHz,
MeOD) δ 7.54–7.39 (m, 2H), 7.04–6.85 (m, 2H),
6.78 (s, 1H), 3.81 (s, 3H) ppm. 13C NMR (100 MHz, MeOD)
δ 158.3, 150.0, 125.9, 125.0, 124.8, 113.6, 109.8, 54.3 ppm.
LCMS (ESI+) m/z 190.2 [M + H]+, retention time 1.61 min (100%). HRMS (ESI+): m/z calculated for C10H12N3O [M + H]+: 190.0980. Found: 190.0979. NMR data
is in accordance with literature values.[37]
4-(4-Iodophenyl)-1H-imidazol-2-amine (54)
Following general method G, from N-(4-(4-iodophenyl)-1H-imidazol-2-yl)acetamide 49 (150 mg, 0.46 mmol) was obtained 54 (103 mg,
0.36 mmol, 79% yield) as a red solid. 1H NMR (500 MHz,
DMSO-d6) δ 7.58 (d, J = 8.3 Hz, 2H), 7.38 (d, J = 8.2 Hz, 2H), 7.03 (s,
1H), 5.35 (br s, 2H) ppm. 13C NMR (125 MHz, DMSO-d6) δ 150.8, 137.4, 134.6, 133.4, 126.0,
111.1, 90.2 ppm. LCMS (ESI+) m/z 286.1 [M + H]+, retention time 2.09 min (100%). HRMS
(ESI+): m/z calculated for C9H9IN3 [M + H]+: 285.9841.
Found: 285.9846.
4-(1H-Imidazol-4-yl)phenol
(55)
2-Bromo-1-(4-hydroxyphenyl)ethan-1-one
(1.0 g, 4.67 mmol)
was dissolved in formamide (5 mL), and the reaction mixture was heated
at 150 °C for 24 h. After cooling to rt, the resulting mixture
was diluted with EtOAc (40 mL) and washed with saturated aqueous NaHCO3 (40 mL). The aqueous layer was extracted with EtOAc (2 ×
40 mL), and the combined organic layers were dried over MgSO4, filtered, and concentrated to give 55 (0.75 g, 4.67
mmol, 99% yield) as a brown oil. 1H NMR (400 MHz, MeOD)
δ 8.17 (br s, 1H), 7.68 (d, J = 1.2 Hz, 1H),
7.64–7.44 (m, 2H), 7.25 (d, J = 1.1 Hz, 1H),
6.82 (d, J = 8.6 Hz, 2H) ppm.13C NMR (100
MHz, MeOD) δ 156.3, 135.1, 125.8, 124.3, 115.1 ppm. LCMS (ESI+) m/z 161.2 [M + H]+, retention
time 0.54 min (100%). HRMS (ESI+): m/z calculated for C9H9N2O [M + H]+: 161.0715. Found: 161.0720. NMR data is in accordance with
literature values.[38]
4-(4-Methoxyphenyl)-1-methyl-1H-imidazole (56)
To a solution of 55 (2.77 g, 17
mmol) in DMF (20 mL) containing cesium carbonate (12.4 g, 38 mmol)
at rt was added iodomethane (2.4 mL, 38 mmol). The reaction mixture
was stirred at rt for 8 h. The reaction mixture was then cooled to
0 °C and quenched with H2O (70 mL) and extracted with
EtOAc (3 × 60 mL). The combined organic layers were washed with
brine (3 × 100 mL), dried over MgSO4, filtered, and
concentrated to give 56 (2.8 g, 15 mmol, 88% yield) as
a brown solid. 1H NMR (400 MHz, CDCl3) δ
7.71 (d, J = 8.8 Hz, 2H), 7.46 (d, J = 1.3 Hz, 1H), 7.10 (d, J = 1.4 Hz, 1H), 6.94 (d, J = 8.9 Hz, 2H), 3.85 (s, 3H), 3.73 (s, 3H) ppm. 13C NMR (100 MHz, CDCl3) δ 158.6, 142.4, 137.8, 127.2,
126.0, 114.8, 114.0, 55.3, 33.5. LCMS (ESI+) m/z 189.2 [M + H]+, retention time 1.38 min (100%).
HRMS (ESI+): m/z calculated for
C11H13N2O [M + H]+: 189.1028.
Found: 189.1031. NMR data is in accordance with literature values.[41]
4-(1-Methyl-1H-imidazol-4-yl)phenol
(57)
A solution of 4-(4-methoxyphenyl)-1-methyl-1H-imidazole 56 (1.21 g, 6.42 mmol) in anhydrous
DCM (30 mL) at −78 °C was treated with BBr3 (16 mL, 16.05 mmol, 1 M in DCM). The reaction mixture was stirred
at −78 °C for 10 min and 2 h at rt. The reaction mixture
was then cooled to −78 °C and quenched with MeOH (4 mL).
The solvents were then removed under reduced pressure, and the crude
product was redissolved in EtOAc (50 mL) and washed with saturated
NaHCO3 solution. The aqueous phase was extracted with EtOAc
(2 × 50 mL), and the combined organic fractions were dried with
anhydrous MgSO4 and filtered off, and the solvent was removed
under vacuum to give 57 (1.30 g, 6.17 mmol, 96% yield)
as a brown solid, which was used in the next step without further
purification. 1H NMR (400 MHz, MeOD) δ 8.99–8.80
(m, 1H), 7.77 (d, J = 1.6 Hz, 1H), 7.55 (d, J = 8.7 Hz, 2H), 6.93 (d, J = 8.8 Hz, 2H),
3.99 (s, 3H) ppm. 13C NMR (100 MHz, MeOD) δ 159.1,
135.1, 134.6, 126.9, 117.3, 117.2, 115.8, 35.0 ppm. LCMS (ESI+) m/z 175.1 [M + H]+, retention
time 1.09 min (100%). HRMS (ESI+): m/z calculated for C10H11N2O [M + H]+: 175.0871. Found: 175.0870.
Cs2CO3 (3.46 g, 10.62
mmol) and methyl 2-bromopropanoate (1.29 g, 6.81 mmol) were added
to a solution of phenol derivative 57 (1.12 g, 5.32 mmol)
in anhydrous DMF (10 mL), and the mixture was stirred at rt for 1
h and at 60 °C for 5 h. After cooling to rt, the reaction was
diluted with water (80 mL) and extracted with EtOAc (2 × 80 mL).
The combined organic layers were washed with brine (3 × 100 mL),
dried over MgSO4, filtered, and concentrated. The crude
product was purified by flash chromatography (0–20% v/v MeOH
in DCM) to give 58 (1.30 g, 4.99 mmol, 92% yield) as
an orange oil. 1H NMR (400 MHz, CDCl3) δ
7.74–7.61 (m, 2H), 7.45 (d, J = 1.3 Hz, 1H),
7.08 (d, J = 1.4 Hz, 1H), 6.96–6.81 (m, 2H),
4.81 (q, J = 6.8 Hz, 1H), 3.77 (s, 3H), 3.71 (s,
3H), 1.64 (d, J = 6.8 Hz, 3H) ppm. 13C
NMR (100 MHz, CDCl3) δ 172.8, 156.5, 142.1, 137.8,
128.1, 126.0, 115.3, 115.1, 72.7, 52.3, 33.5, 18.6 ppm. LCMS (ESI+) m/z 261.3 [M + H]+, retention
time 1.50 min (100%). HRMS (ESI+): m/z calculated for C14H17N2O3 [M + H]+: 261.1239. Found: 261.1232.
Methyl 2-(4-(1-methyl-1H-imidazol-4-yl)phenoxy)propanoate 58 (0.53 g, 2.04 mmol)
was dissolved in a 2:1 v/v mixture of THF and H2O (7.5
mL), and then NaOH (0.16 g, 4.08 mmol) was added. The reaction mixture
was heated at 80 °C for 40 min. After the reaction mixture was
allowed to cool to rt, the solvents were removed. The mixture was
then diluted with H2O (5 mL) and acidified to pH 7–8
using 1 N HCl. The mixture was cooled, and the resulting precipitate
was collected by filtration, washed with water (2 × 3 mL) and
petroleum ether (2 × 3 mL), and dried in vacuo to afford acid 59 (0.33 g, 1.33 mmol, 65% yield) as
a white solid. 1H NMR (400 MHz, DMSO-d6) δ 7.69–7.55 (m, 3H), 7.46 (d, J = 1.3 Hz, 1H), 6.84 (d, J = 8.8 Hz, 2H),
4.81 (q, J = 6.8 Hz, 1H), 3.66 (s, 3H), 1.49 (d, J = 6.7 Hz, 3H). 13C NMR (100 MHz, DMSO-d6) δ 173.3, 156.1, 140.4, 138.2, 127.8,
125.4, 115.8, 114.8, 71.7, 33.1, 18.4. LCMS (ESI+) m/z 247.2 [M + H]+, retention time 1.43
min (100%). HRMS (ESI+): m/z calculated
for C13H15N2O3 [M + H]+: 247.1083. Found: 247.1080.
To a solution of 57 (212 mg, 1.21 mmol) in
anhydrous THF (5 mL) was added ethyl d-lactate (215 mg, 1.82
mmol). After the reaction mixture was cooled to 0 °C, PPh3 (478 mg, 1.82 mmol) and DEAD (286 μL, 1.82 mmol) were
added, and the reaction mixture was stirred for 16 h at rt. The reaction
mixture was poured into ice–water (40 mL) and extracted with
DCM (3 × 50 mL). The combined organic layers were washed with
brine, dried over anhydrous MgSO4, filtered, and concentrated.
The crude product was purified by flash chromatography (0–15%
v/v MeOH in DCM) to give ethyl (S)-2-(4-(1-methyl-1H-imidazol-4-yl)phenoxy)propanoate as an orange solid. LCMS
(ESI+) m/z 275.3 [M + H]+, retention time 2.06 min (94%). Ethyl (S)-2-(4-(1-methyl-1H-imidazol-4-yl)phenoxy)propanoate (200 mg, 0.73 mmol) was
dissolved in a 2:1 v/v mixture of THF and H2O (7.5 mL),
and then NaOH (58 mg, 1.46 mmol) was added. The reaction mixture was
stirred at rt for 3 h, and the solvents were removed. The residue
was dissolved in EtOAc (40 mL) and then extracted with H2O (40 mL). The aqueous layer was acidified to pH 7–8 using
1 N HCl. The mixture was cooled, and the resulting precipitate was
collected by filtration, washed with water (2 × 3 mL) and petroleum
ether (2 × 3 mL), and dried in vacuo to afford
acid ( (100 mg, 0.40 mmol, 56%
yield) as a white solid. [α]D25 −51.6 (c 1.0, CHCl3). LCMS (ESI+) m/z 247.2
[M + H]+, retention time 1.43 min (100%). HRMS (ESI+): m/z calculated for C13H15N2O3 [M + H]+: 247.1083.
Found: 247.1079.
4-(4-Bromophenyl)oxazol-2-amine (61)
A
mixture of 2-bromo-1-(4-bromophenyl)ethanone (200 mg, 0.72 mmol) and
urea (435 mg, 7.2 mmol) in anhydrous acetonitrile (5 mL) was heated
at 80 °C for 18 h. The solvent was removed, and the residue was
purified by flash chromatography (0–10% v/v MeOH in DCM) to
afford 61 (134 mg, 0.56 mmol, 78% yield) as a yellow
solid. 1H NMR (400 MHz, MeOD) δ 7.68 (s, 1H), 7.59–7.42
(m, 4H) ppm. 13C NMR (125 MHz, MeOD) δ 163.9, 139.8,
132.7, 132.1, 128.8, 127.9, 122.1 ppm. LCMS (ESI+) m/z 238.8 [M + H]+, retention time 1.95
min (100%). HRMS (ESI+): m/z calculated
for C9H8BrN2O [M + H]+: 238.9820. Found: 238.9816.
Enzyme
Assay
The activity of Mth IMPDH ΔCBS
was determined using a plate reader
by monitoring the production of NADH in absorbance at 340 nm and corrected
for noncatalyzed chemical reactions in the absence of Mth IMPDH ΔCBS. All the measurements were done in the assay buffer
(50 mM Tris HCl pH 8, 1 mM DTT, 1 mM EDTA, and 100 mM KCl) at 37 °C
with 20 nM Mth IMPDH ΔCBS, 2.8 mM NAD+, and 1 mM IMP in a total of 150 μL volume in a 96 well plate-based
format, and data were collected for 32 min. The reaction was initiated
by the addition of the substrate, IMP, at a concentration of 1 mM.
All reactions were performed in triplicate. Prior to reaction initiation,
the compounds were preincubated in a buffer with enzyme for 5 min.
The inhibitors were dissolved in DMSO-d6 and diluted to a final concentration of 1% v/v in experimental reactions.IC50 values were calculated by plotting the percentage
of inhibition against the logarithm of inhibitor concentration, and
dose–response curves were fitted using Prism software (GraphPad).The Ki value for NAD+ was
determined at a constant saturating IMP concentration (1 mM) and five
different concentrations of NAD+ (0.35, 0.70, 1.0, 1.4,
and 2.8 mM) in the presence of increasing concentrations of inhibitor.
The value of Ki for IMP was determined
at fixed saturating concentration of NAD+ (2.8 mM) and
different concentrations of IMP (0.12, 0.18, 0.25, 0.5, and 1 mM)
and inhibitor.The initial velocities at various inhibitor concentrations
were
determined based on the slope in the linear part of each reaction
containing the inhibitor and the uninhibited reaction. To determine
the inhibition constant (Ki values), the
initial rate data versus substrate concentration at different inhibitor
concentrations were fit using Prism software (GraphPad) to equations
for uncompetitive or mixed inhibition. For each inhibitor concentration,
the reciprocal of enzyme reaction velocity versus the reciprocal of
the substrate concentration was plotted in a Lineweaver–Burk
plot to determine the pattern of inhibition.
Protein
Purification, Crystallization, and Data
Collection of Mth IMPDH ΔCBS
Mth IMPDH was expressed, purified, and crystallized as previously
described.[24,25] Briefly, hexahistidine tagged Mth IMPDH ΔCBS in pHat2 was expressed overnight in
BL21 DE3 (NEB) cells at 18 °C by addition of 500 μM IPTG.
Cells were lysed in 50 mM Hepes, pH 8.0, 500 mM NaCl, 5% glycerol,
10 mM β-mercaptoethanol, and 20 mM imidazole, and the recombinant
protein purified using a Hi-Trap IMAC FF column (GE Healthcare) charged
with nickel and an elution gradient of up to 300 mM imidazole. The
hexahistidine tag was cleaved by TEV protease, and the purified Mth IMPDH ΔCBS was obtained by negative nickel gravity-flow
purification[42] and size exclusion chromatography
on a Superdex 200 gel filtration column equilibrated in 20 mM Hepes
pH 8.0, 500 mM NaCl, 5% glycerol, and 1 mM TCEP step. The recombinant Mth IMPDH ΔCBS was then concentrated to 12.5 mg/mL
for crystallization.Mth IMPDH ΔCBS protein
crystallized in 1 μL + 1 μL hanging drops with 100 mM
sodium acetate, pH 5.5, 200 mM calcium chloride, and 8–14%
isopropanol. Crystals were soaked overnight in drops of well solution
+ 5 mM IMP and either 5 mM Fragment 2 or Compound 1 or 1 mM Compound 31. Crystals were cryoprotected
by passing through drops containing well solution + 25% glycerol and
flash-frozen in liquid nitrogen. Data were collected from the crystals
at Diamond Light Source beamline.
Structure
Solution, Ligand Fitting, and Refinement
Data were processed
using XDS[43] and
Pointless (CCP4). To solve the structure, molecular replacement was
performed with Phenix Phaser[44] using a
previously solved IMP-bound Mth IMPDH ΔCBS
structure as a probe (PDB IDs: 5J5R; 5K4X; 5K4Z). Refinement was performed using Phenix.refine
and manually in Coot.[45] IMP and the inhibitors
were sequentially fitted into the density using the LigandFit function
of Phenix, and the structures were manually refined further using
Coot. Information regarding the crystallographic statistics can be
found in Table S2. Protein–ligand
interactions were analyzed using Arpeggio[46] and CSM-Lig.[47,48] Compound properties were evaluated
using pkCSM.[49] All figures made using Pymol
(Schrodinger).
Drug Susceptibility Testing
against Mtb
An Alamar Blue fluorescence-based
broth microdilution
assay was used to assess the minimum inhibitory concentration (MIC)
of compounds against Mtb H37Rv, as described previously.[25,50] Briefly, Mtb H37Rv was grown in standard Middlebrook
7H9 broth (BD) supplemented with OADC (BD), 0.2% glycerol, and 0.05%
Tween-80 to midexponential phase. Compounds dissolved in DMSO (1%)
were tested in clear-bottomed, round-well 96-well microtiter plates
at eight different concentrations using the standard anti-TB drugs,
rifampin and isoniazid, as positive controls. An inoculum of ∼105 bacteria was added to each well, and the plates were incubated
at 37 °C for 7 d. On day 7, 10 μL of Alamar Blue (Invitrogen)
was added to each well, and plates were further incubated at 37 °C
for 24 h. The fluorescence (excitation 544 nm; emission 590 nm) was
measured in a FLUOstar OPTIMA plate reader (BMG LABTECH, Offenberg,
Germany). Data were normalized to the minimum and maximum inhibition
controls to generate a dose response curve (% inhibition) from which
the MIC90 was determined.
Authors: Philip Onyebujoh; Alimuddin Zumla; Isabella Ribeiro; Roxana Rustomjee; Peter Mwaba; Melba Gomes; John M Grange Journal: Bull World Health Organ Date: 2005-11-10 Impact factor: 9.408
Authors: Evan W Orenstein; Sanjay Basu; N Sarita Shah; Jason R Andrews; Gerald H Friedland; Anthony P Moll; Neel R Gandhi; Alison P Galvani Journal: Lancet Infect Dis Date: 2009-03 Impact factor: 25.071
Authors: Yumi Park; Angela Pacitto; Tracy Bayliss; Laura A T Cleghorn; Zhe Wang; Travis Hartman; Kriti Arora; Thomas R Ioerger; Jim Sacchettini; Menico Rizzi; Stefano Donini; Tom L Blundell; David B Ascher; Kyu Rhee; Ardala Breda; Nian Zhou; Veronique Dartois; Surendranadha Reddy Jonnala; Laura E Via; Valerie Mizrahi; Ola Epemolu; Laste Stojanovski; Fred Simeons; Maria Osuna-Cabello; Lucy Ellis; Claire J MacKenzie; Alasdair R C Smith; Susan H Davis; Dinakaran Murugesan; Kirsteen I Buchanan; Penelope A Turner; Margaret Huggett; Fabio Zuccotto; Maria Jose Rebollo-Lopez; Maria Jose Lafuente-Monasterio; Olalla Sanz; Gracia Santos Diaz; Joël Lelièvre; Lluis Ballell; Carolyn Selenski; Matthew Axtman; Sonja Ghidelli-Disse; Hannah Pflaumer; Markus Bösche; Gerard Drewes; Gail M Freiberg; Matthew D Kurnick; Myron Srikumaran; Dale J Kempf; Simon R Green; Peter C Ray; Kevin Read; Paul Wyatt; Clifton E Barry; Helena I Boshoff Journal: ACS Infect Dis Date: 2016-10-17 Impact factor: 5.084
Authors: Jonathan A G Cox; Grace Mugumbate; Laura Vela-Glez Del Peral; Monika Jankute; Katherine A Abrahams; Peter Jervis; Stefan Jackenkroll; Arancha Perez; Carlos Alemparte; Jorge Esquivias; Joël Lelièvre; Fernando Ramon; David Barros; Lluis Ballell; Gurdyal S Besra Journal: Sci Rep Date: 2016-12-16 Impact factor: 4.379
Authors: Michael J Pepi; Shibin Chacko; Gary M Marqus; Vinayak Singh; Zhe Wang; Kyle Planck; Ryan T Cullinane; Penchala N Meka; Deviprasad R Gollapalli; Thomas R Ioerger; Kyu Y Rhee; Gregory D Cuny; Helena I M Boshoff; Lizbeth Hedstrom Journal: ACS Infect Dis Date: 2022-01-11 Impact factor: 5.578
Authors: Gyan Modi; Gary M Marqus; Mohana Rao Vippila; Deviprasad R Gollapalli; Youngchang Kim; Adhar C Manna; Shibin Chacko; Natalia Maltseva; Xingyou Wang; Ryan T Cullinane; Yubo Zhang; Judy L M Kotler; Petr Kuzmic; Minjia Zhang; Ann P Lawson; Andrzej Joachimiak; Ambrose Cheung; Barry B Snider; David M Rothstein; Gregory D Cuny; Lizbeth Hedstrom Journal: ACS Infect Dis Date: 2021-09-30 Impact factor: 5.578
Authors: Malancha Karmakar; Carlos H M Rodrigues; Kathryn E Holt; Sarah J Dunstan; Justin Denholm; David B Ascher Journal: PLoS One Date: 2019-05-29 Impact factor: 3.240