Parthiban et al.
acid. Crystal structure of MTB LAT is available with its
substrate lysine (10), and we have designed novel inhibi-
tors using e-pharmacophore approach. One of the most
active lead was derivatised to study SAR and its activity
against nutrient-starved MTB model.
erated and further used for virtual screening of a commer-
cial database Asinex containing five lakh molecules and
our in-house database of 3000 molecules, with must
match five features of six-point pharmacophore features.
Compounds retrieved by e-pharmacophore filter using
Phase with a fit value above 1.0 were regarded as poten-
tial hits and were carried forward for high-throughput vir-
tual screening. Compounds resulting with a score of
≥ ꢀ4.0 kcal/mol and with two or more hydrogen bonds
were subjected to another round of docking by Glide XP.
Glide XP scores combine accurately physics-based scor-
ing terms and thorough sampling and result in compounds
with docking scores between ꢀ7.649 and ꢀ6.222 kcal/
mol. The final short listing of the molecules was based on
the protein–ligand interaction in the active site cleft through
hydrogen bonding with Arg170, Gln274, Ala129, Gly128,
and p–p stacking interactions with Phe167. The 2D repre-
sentation of top hits identified from Asinex database has
been represented in Figure 2. Finally, top seven hits were
selected from the Asinex and our in-house database and
were evaluated biologically against MTB LAT enzyme inhi-
bition studies at various concentrations from 100 to 10 lM.
In this study, our effort is to identify novel small molecule
inhibitors that bind to the active site; this is the place
where the substrate binds and reaction occurs of MTB
LAT protein started with e-pharmacophore approach,
which combined both the aspects of structure-based and
ligand-based techniques using the crystal structure of LAT
from MTB in external aldimine form in complex with its
ꢀ
substrate lysine (PDB code: 2CJD) with 2.00 A resolution.
Pharmacophore hypotheses based on mapping of the
energetic terms from the extra precision (XP) Glide scoring
function onto atom centers were obtained using PHASE
(
PHASE, v3.3; Schr o€ dinger, LLC, New York, NY, USA) with
the default set of six chemical features: hydrogen bond
acceptor (A), hydrogen bond donor (D), hydrophobic (H),
negative ionizable (N), positive ionizable (P), and aromatic
ring (R). The Glide XP scoring terms were computed, and
the energies were mapped onto atoms. The pharma-
cophore sites were generated, and the Glide XP energies
from the atoms that comprised each pharmacophore sites
were summed up. These sites were then ranked based on
the individual energies, and the most favorable sites were
selected for the pharmacophore hypothesis. This pharma-
cophore model was then used as query for virtual screen-
ing. A six-point e-pharmacophore model was generated
with MTB LAT protein domain. The pharmacophoric sites
established included three hydrogen bond acceptor (A1,
A3, and A4), one ring aromatic (R14), and two negative
ionizable group (N9 and N11) (Figure 1), which were gen-
Among the seven compounds tested against MTB LAT,
four molecules showed IC50 of ≤ 25 lM. The most promis-
ing compounds from the in vitro MTB LAT assay, Lead 1
having IC50 of 18.06 lM, were then evaluated in more
detail to provide a structural basis for the further optimiza-
tion by chemical synthesis. In silico molecular docking
studies were carried out with the crystal structure of MTB
LAT in external aldimine form in complex with its substrate
L-lysine with bound PMP [PDB ID: 2CJD]. The reference
ligand (L-lysine with bound PMP) was re-docked with the
active site residues of the MTB LAT to validate the active
site cavity. The ligand exhibited
a Glide score of
ꢀ
5.81 kcal/mol and was found in the vicinity of Arg170,
Glu243, Arg422, Val63, Gly169, Gly127, Gly128, Ala129,
Val132, Phe167, Asn271, Val273, Glu238, Gly242, and
Gln274 amino acid residues. Re-docking results showed
that the compound exhibited similar interactions as that of
ꢀ
the original crystal structure with RMSD of 0.97 A sug-
gesting reliability of the docking method (Figure 3).
The docking orientation of the Lead 1, when compared
with the crystal ligand, the lead compound, revealed
important hydrogen bonding interaction with Arg170 amino
acid residue. The docking score of lead was found to be
ꢀ
7.06 kcal/mol, and a strong metal co-ordination bond
with the oxygen atom of Ash271 was observed as shown
in Figure 4. In addition to hydrogen bonding interactions,
the phenyl ring interacted via p–p interaction with phenyl
ring of Phe167 amino acid residue, analogues to the one
observed with the crystal ligand L-lysine. The orientation of
lead compound was found to be closely similar to that of
crystal ligand exhibiting a very good fitness in the active
site pocket, making this compound the most active
(IC50 = 18.06 lM) as shown in Table 1.
Figure 1: The energy based pharmacophoric features of the
reference ligand bound to the protein.
2
66
Chem Biol Drug Des 2016; 87: 265–274