Toward antituberculosis drugs: in silico screening of synthetic compounds against Mycobacterium tuberculosisl,d-transpeptidase 2

Abstract

Mycobacterium tuberculosis (Mtb) the main causative agent of tuberculosis, is the main reason why this disease continues to be a global public health threat. It is therefore imperative to find a novel antitubercular drug target that is unique to the structural machinery or is essential to the growth and survival of the bacterium. One such target is the enzyme l,d-transpeptidase 2, also known as Ldt_Mt2, a protein primarily responsible for the catalysis of 3→3 cross-linkages that make up the mycolyl–arabinogalactan–peptidoglycan complex of Mtb. In this study, structure-based pharmacophore screening, molecular docking, and in silico toxicity evaluations were employed in screening compounds from a database of synthetic compounds. Out of the 4.5 million database compounds, 18 structures were identified as high-scoring, high-binding hits with very satisfactory absorption, distribution, metabolism, excretion, and toxicity properties. Two out of the 18 compounds were further subjected to in vitro bioactivity assays, with one exhibiting a good inhibitory activity against the Mtb H37Ra strain.

Keywords:

• antituberculosis drug discovery
• virtual screening
• docking

Introduction

The cell wall of Mycobacterium tuberculosis (Mtb) is made up of a lipid membrane interspersed with cell wall proteins and complex layers of peptidoglycan, arabinogalactan, and mycolic acids known as the cell wall core or the mycolyl–arabinogalactan–peptidoglycan (mAGP) complex. While the part of the cell wall composed of lipids and proteins can be easily disrupted using solvents, the mAGP complex remains insoluble, making it indispensable for the survival of the pathogen and providing resistance to common antibiotics. Thus, the mAGP complex is a very attractive target for drug development against tuberculosis (TB).¹ The peptidoglycan layer of mycobacteria, as well as other β-lactam resistant bacteria, is predominantly made of 3→3 linkages instead of the more common 4→3 cross-linkages. l,d-transpeptidases catalyze the formation of 3→3 cross-linkages, one type of which is the l,d-transpeptidase 2 (Ldt_Mt2; MT2594, product of gene Rv2518c), a protein that was observed to be pertinent for Mtb’s virulence and growth during the chronic phase of the disease.^2,3 While the combination of clavulanate and meropenem shows potential in treating patients with TB by inhibiting Mtb’s mycobacterial β-lactamase (BlaC) and Ldt_Mt2, respectively, the short half-life of meropenem forestalls its use as an anti-TB agent.⁴ Furthermore, both imipenem, another carbapenem that has shown activity against Ldt_Mt2, and meropenem, are compounds that are currently being tested for the treatment of a wide range of infections other than TB.⁵ Nonetheless, this may later lead to resistance, as the drug is not specific for TB.

TB, an infectious disease caused by Mtb, is one of the major health concerns worldwide, resulting in almost 1.5 million deaths in the year 2014.⁶ Despite the presence of current drug treatments for TB, including the use of isoniazid, rifampicin, ethambutol, pyrazinamide, and streptomycin,⁷ the incidence of TB remains high, especially in developing countries, such as the Philippines.⁶ The alarming rise of multi- and extensively drug-resistant TB patients has become a serious global health threat. The emergence of resistant strains is due to problems, such as the arrival of human immunodeficiency virus⁸ and poor patient compliance because of the extensive treatment regimen.⁹ Thus, new anti-TB drugs that can shorten the treatment regimen and/or target the resistant TB strains are urgently needed.

Depending on the target disease and treatment method, the whole process of drug discovery and development can take up to 15 years and cost millions of dollars before the drug can actually reach the market. Out of the millions of candidate drug compounds initially screened for a certain disease, only a few make it to clinical trials, and even after that, <10% of those compounds from clinical trials successfully get the final approval.¹⁰ For TB, the initial phase of drug discovery aims to identify “leads” with anti-TB activity and desirable physicochemical, pharmacokinetic, and toxicity properties through classic wet laboratory testing.¹¹ However, this approach requires Biosafety Level 3 (BSL3) laboratory equipment and expert research skills in handling sensitive protocols used during drug screening.¹² Also, with the continuous advancements in technology and the availability of the Mtb genome and three-dimensional (3D) structure of potential enzyme targets of TB, virtual target-based screening can be utilized. The use of computational screening methods can potentially reduce the cost, time, and effort needed for the initial screening of candidate drug compounds with pharmacological activity against TB.¹³

With the use of a powerful computer, a computational software, and structural data of a protein target and compound libraries, initial screening of millions of compounds can be performed in less time than what is used in the classic drug discovery and development process, helping to prioritize compound testing and minimizing randomization in the laboratory. In this work, >4 million synthetic compounds were screened based on a pharmacophore that satisfies the electronic and structural requirements of the drug target’s binding site. The high-scoring hits were subsequently docked to the target and were rank-ordered based on their binding energies. The high-affinity hits were further evaluated in silico for their potential pharmacokinetics and pharmacodynamics properties.

Materials and methods

All computational work was performed using Accelrys Discovery Studio 4.0 (DS 4.0) on a Windows 7 Home Edition with an Intel^® Core™ i7-3770 3.40 GHz quad core processor, 4 GB RAM, and 64-bit operating system. Protein structures were downloaded from Research Collaboratory for Structural Bioinformatics protein databank, and imipenem and meropenem structures were taken from the National Center for Biotechnology Information website. Enamine real database containing compound structures was downloaded from the enamine website.¹⁴

Structure-based pharmacophore modeling

Preparation of 3D protein structure and library compounds

The 3D structure of Ldt_Mt2 complexed with a peptidoglycan fragment (PDB ID: 3TUR) solved at 1.72 Å resolution² was retrieved. The bound peptidoglycan fragment was removed, and the protein was prepared using the Prepare Protein protocol of DS 4.0 (BIOVIA, Tokyo, Japan) using the default parameters. The Prepare Protein protocol primes the protein for input into other protocols in DS 4.0 by inserting missing atoms in incomplete residues, optimizing side-chain conformation, modeling missing loop regions, removing alternate conformations, and protonating titratable residues at pH 7.4.¹⁵ The enamine compound database was downloaded and prepared using the Prepare Ligands protocol. The compounds in the enamine database were prepared using Prepare Ligands protocol.

Optimization of protein structure and root-mean-square deviation

Minimization protocol was used to optimize the protein structure for screening. The default algorithm parameter, Smart Minimizer, was used to minimize the structure by executing 1,000 steps of steepest descent using an RMS gradient acceptance of 3, followed by conjugate gradient minimization, which locates an unconstrained local minimum for the input structure.^15,16 The root-mean-square deviation (RMSD) of the prepared protein structure was then calculated against the original protein file using the Superimpose Proteins tool. The protein structures were superimposed based on Cα pairs.

Generation of structure-based pharmacophore model

The binding site of Ldt_Mt2 was identified based on literature’s description, that is, the site that contains the catalytic triad Cys354, His336, and Ser337.^2,17,18 After identification of the binding site, a binding sphere was generated using the Binding Site tool in DS 4.0 with a radius of 10 Å. The Interaction Generation tool of DS 4.0 was used to generate a pharmacophore model that complements the chemical features (hydrophobic, H-donor, and H-acceptor) in the protein’s active site. The Edit and Cluster Pharmacophore tool was used to cluster the common pharmacophore properties down to <30 features.

Virtual screening of compounds

Preparation of 3D compound libraries

Approximately 4.5 million database compounds were screened in this work. The test compounds, as well as imipenem and meropenem, were prepared using the Prepare Ligands protocol with default parameters. The Prepare Ligands protocol primes the ligands for use in other protocols by removing duplicate structures, generating isomers and tautomers, generating 3D conformations, and other functions specified by the user.¹⁵

Database building

The Build 3D Database protocol was used to create compound databases for easier screening. The compound database was built based on Catalyst algorithms, which create compact, indexed compound databases used for pharmacophore screening.¹⁵

Pharmacophore-based screening

The generated structure-based pharmacophore model was employed to screen the compound databases using the Screen Library protocol, which uses the flexible search method of Catalyst. The Screen Library protocol enumerates numerous possible structures from an input query pharmacophore model and screens a compound 3D database. Screening was carried out twice, one for rigid fitting method and another for flexible fitting method. The rigid fitting method screens the compounds without modifying their input conformation, whereas in the flexible fitting method, each ligand conformation was slightly modified to better fit the pharmacophore model.¹⁵ The base fit value for rigid screening was arbitrarily set to 2.5, while the base fit value for flexible screening was set to 3.0. Any compound with a fit value <2.5 and <3.0 was not chosen for subsequent molecular docking and in silico toxicity screening, respectively.

Molecular docking

The hit compounds from the pharmacophore screening were docked to the prepared Ldt_Mt2 active site using the CHARMm-based DOCKER (CDOCKER) docking protocol. CDOCKER is a grid-based molecular dynamics-simulated-annealing-based algorithm docking procedure that utilizes CHARMm force fields. It allows full ligand flexibility in the docking process by producing several ligand poses when the ligand is docked into the receptor’s binding site and by applying molecular dynamics-based simulated annealing and in situ minimization.¹⁹ Meropenem and imipenem, carbapenems known to inhibit Ldt_Mt2,^2,18 were also docked and compared to the original bound conformation to validate the docking method and compare binding affinities and active site interaction.

Calculation of binding energies

Calculate Binding Energies protocol of DS4.0 was used to calculate the binding affinity of all docked compounds. This protocol computes for the binding energy using the formula energy of binding = energy of complex − energy of ligand − energy of receptor.¹⁵ The calculated binding energies of meropenem and imipenem were used as the baseline comparison for the selection of compounds with the best binding affinity to Ldt_Mt2.

In silico absorption, distribution, metabolism, excretion, and toxicity and toxicity prediction by komputer assisted technology screening

The compounds with the best binding affinity were further screened in silico for their pharmacokinetics and pharmacodynamics properties using the absorption, distribution, metabolism, excretion, and toxicity (ADMET) and toxicity prediction by komputer assisted technology (TOPKAT) protocols of DS 4.0. The solubility, absorption, plasma protein binding, CYP2D6 inhibition, and hepatotoxicity of each compound were evaluated by the use of ADMET.¹⁵ The candidate compounds were also subjected to TOPKAT calculations to determine the probability for carcinogenicity, mutagenicity, and other toxicity measures.

Resazurin-based microtiter plate assay

Bacterial strains and culture conditions

Mtb H37Ra (ATCC 25177) glycerol stocks, provided by the Marine Natural Products Laboratory at the University of the Philippines Marine Science Institute, were first thawed out and subcultured on Middlebrook 7H11 agar supplemented with 10% oleic acid–albumin–dextrose–catalase (Titan Media, Delhi, India). Plates were then incubated at 37°C for 3–4 weeks and then subcultured on Middlebrook 7H9 broth supplemented with 10% albumin–dextrose–catalase (Titan Media). Broth tubes were then incubated in a shaking incubator at 37°C, 150 rpm for another 3–4 weeks.

Compound preparation

Enamine compounds were all procured from Enamine Ltd (Kiev, Ukraine) and then solubilized at a stock concentration of 2 mg/mL in dimethyl sulfoxide (DMSO). Rifampicin (Sigma-Aldrich Co., St Louis, MO, USA), the assay positive control, was also solubilized in DMSO at a stock concentration of 1 mg/mL. All the compound stock solutions were stored at −20°C.

Resazurin microtiter assay

Preparation of the assay inoculum was first done by adjusting the turbidity of the broth culture to match a McFarland No 1 standard (A625 nm ≈0.25). The adjusted culture was then diluted further to a 1:49 mixture of culture: M7H9 broth. This resulting inoculum was used for the assay proper within 30 minutes of preparation.

Stock solutions of enamine compounds and rifampicin were thawed and diluted to arrive at two concentrations: one at a high concentration of 2,500 μM, and the other at a low concentration of 20 μM. Since each test well has a final volume of 200 μL and only 2 μL of the compound is added per well, the final drug and rifampicin well concentrations were 25 μM and 0.2 μM. Each sample was tested in quadruplicates.

The resazurin assay proper was performed in 96-well, flat-bottomed microtiter plates. Four types of assay wells were prepared. First, sample wells consisted of 98 μL of Middlebrook 7H9 broth supplemented with 10% ADC, 100 μL of the H37Ra inoculum, and 2 μL of the prepared enamine compound. On the other hand, the rifampicin positive control wells were added with 98 μL of M7H9 broth with 10% ADC, 100 μL of the inoculum, and 2 μL of the prepared rifampicin standard. Sterility control wells were also prepared, which consisted of 2 μL of DMSO and 198 μL of M7H9 broth with 10% ADC. Moreover, growth control wells were prepared, which consisted of 98 μL of M7H9 broth with 10% ADC, 100 μL of H37Ra inoculum, and 2 μL of DMSO. Lastly, four media-only wells were also added with 200 μL of the M7H9 broth. Well plates were then incubated at 37°C, 150 rpm for 7 days in a shaking incubator. After 7 days, 20 μL of 0.02% resazurin (Sigma-Aldrich Co.) was added to all the wells and then incubated for an additional 24 hours. The fluorescence readings of all wells were then read at an excitation filter of 530 nm and an emission filter of 590 nm.

Ligand interaction analysis

Ligand interactions for both the originally bound and docked meropenem were analyzed with the aid of 2D protein–ligand interaction diagrams. Ligand interactions of the hit compounds were also analyzed to compare the interactions with that of the known inhibitor.

Results and discussion

Crystal structure of Ldt_Mt2

Most bacteria rely on the action of d,d-transpeptidases in catalyzing the biosynthesis of the peptidoglycan layer of their cell walls, particularly in the formation of 4→3 cross-linkages (Figure 1B). A structurally unrelated trans-peptidase also present in most bacteria, l,d-transpeptidase, only plays a minimal role in the organisms’ peptidoglycan synthesis.^2,20 Presently, β-lactams are the mostly widely used class of antibacterial drugs that are designed to target 4→3 transpeptide linkages by acting as suicide substrates of d,dD-transpeptidases.²¹ However, in the case of the nonreplicating Mtb, 3→3 linkages (Figure 1A) predominate the transpeptide network of their peptidoglycan layer, which in turn are primarily catalyzed by l,d-transpeptidases.^22,23

Figure 1 (A) 3→3 linkages and (B) 4→3 linkages of the bacterial peptidoglycan layer.

In 1974, Wietzerbin et al first discovered the existence of these 3→3 linkages in the peptidoglycan layer of mycobacteria.²² Nevertheless, their role and significance in the biochemistry pathway and physiology of mycobacteria remained unknown until recently, when a study conducted by Gupta et al³ showed that a specific Ldt_Mt2 is essential for the maintenance of the virulence and drug resistance of Mtb strains. An inactivation of Ldt_Mt2 resulted in an altered colony morphology, even though no difference in the lipid composition of the cell wall was detected.³ Ldt_Mt2 is the primary l,d-transpeptidase in Mtb; Ldt_Mt2 is expressed at least ten-fold higher than Ldt_Mt1 in all phases of growth of the bacteria, which is evidenced by the fact that continuous Ldt_Mt1 expression was not able to compensate the loss of Ldt_Mt2 activity.^3,18 Loss of Ldt_Mt2 was also found to compromise the bacteria’s ability to acclimate during chronic infection and has displayed an increased susceptibility to β-lactams.³ Therefore, the protein Ldt_Mt2 serves as an excellent target for novel drug discovery against the drug-resistant strains of Mtb.

β-lactams were formerly thought to be ineffectual against Mtb due to the presence of the endogenous BlaC.²⁴ However, a specific class of β-lactams called carbapenems were shown to have activity against l,d-transpeptidases from Enterococcus faecium (Ldt_fm) and Mtb (Ldt_Mt1), since they were poorly hydrolyzed by BlaC.² In addition to this, clavulanic acid is found to be effective in inhibiting Mtb BlaC and can be used in combination with these carbapenems.²⁴ Meropenem and imipenem, as known inhibitors of Ldt_Mt2, were used as the reference in searching for lead-like compounds in this work.

The crystal structure of Ldt_Mt2 with bound peptidoglycan fragment (PDB ID: 3TUR) was used in this study. The structural data of 3TUR contain two monomers of Ldt_Mt2. However, these molecules are not equivalent to a physiological dimer, since it was determined that Ldt_Mt2 behaves as a monomer in solution.² A more recent crystal structure of Ldt_Mt2 in complex with meropenem (PDB ID: 4GSU) that can be used to validate the docking results of meropenem became available.¹⁸ With the release of this new PDB data, the structural similarity of the protein file used in this study (ie, 3TUR) was calculated using 4GSU as a reference, yielding an RMSD value of 0.714 Å. This result indicates very tight structural resemblance between the two proteins. Subsequently, the 3D structure of the Ldt_Mt2 protein with bound peptidoglycan fragments (Figure 2A) was subjected to Prepare Protein and Minimization protocols to clean and fix the protein and find the most stable protein conformation, respectively. The conformation of the 3D structure changed to some extent after optimization procedures were performed, albeit the deviation was minuscule (RMSD =0.66 Å; Figure 2B).

Figure 2 (A) Three-dimensional structure of Ldt_Mt2 protein of Mtb (Ldt_Mt2, PDB ID: 3TUR). (B) Molecular overlay picture of the downloaded Ldt_Mt2 protein structure (green) and prepared Ldt_Mt2 protein structure (pink).

Abbreviations: Ldt_Mt2, l,d-transpeptidase 2; Mtb, Mycobacterium tuberculosis.

Pharmacophore screening, molecular docking, in silico ADMET calculations, and resazurin microtiter assay

As in l,d-transpeptidases from Bacillus subtilis and E. faecium, the YkuD domain of Ldt_Mt2 has a catalytic triad Cys354–His336–Ser337.^2,17,18 These residues were used to identify the active site of Ldt_Mt2, generate the binding sphere, and subsequently, the pharmacophore model. A pharmacophore is a set of steric and electronic characteristics utilized in screening of compounds to help prioritize compounds with potential pharmacological activity against the biological target structure to inhibit its activity.²⁵ A structure-based pharmacophore model identifies areas in a target protein that are hydrophobic, hydrogen donor, or hydrogen acceptor. The generated pharmacophore model for Ldt_Mt2 contains a total of 25 features: seven hydrophobic, ten donor, and eight acceptor. The model was employed to screen 4.5 million compounds from the enamine REAL database before docking to the prepared protein structure.

The binding energies for meropenem was determined to be −112.83 kcal/mol. Hit compounds with high fit values (data not shown) were docked individually to the target protein using CDOCKER, and their affinity was measured by calculating their binding energies. Ninety-four compounds showed superior binding affinity than meropenem (data not shown).

Meanwhile, conventional drug discovery entails the need of wet laboratory analysis with the use of high-throughput screening, after which preclinical tests are performed to determine pharmacokinetic and toxicity properties of the compounds.¹² Frequently, adverse findings are revealed at this late stage of drug discovery and development,²⁶ wasting effort, time, and resources. In silico ADMET evaluation is an expedient step for predictive quantitative structure–property studies that can be applied in drug discovery.²⁷ It reduces the need for expensive and expansive in vitro pharmacokinetics and toxicity screening.

ADMET and TOPKAT protocols of DS 4.0 were employed to furnish (absorption, distribution, metabolism, and excretion) and toxicity information for the hit compounds from the docking results. The ADMET protocol calculates for the probability of a query compound’s absorption, solubility, inhibition of CYP2D6, plasma protein binding, and hepatotoxicity. The TOPKAT module, on the other hand, computes for the probability of an input compound to be carcinogenic, mutagenic, ocular irritant, and other user-defined toxicity measures. In this study, TOPKAT was used to determine the probability of carcinogenicity, mutagenicity, developmental toxicity, and aerobic biodegradability.¹⁵ The 94 high-binding compounds were subjected to in silico ADMET screening, yielding 18 compounds that had the most satisfactory results as shown in Table 1.

Table 1 ADMET and TOPKAT values of meropenem and the top 18 hit compounds

Compound	WOE	AM	DTP	AB	Absorption	Solubility	Hepatotoxicity	CYP2D6 inhibitor	PPB
Meropenem	+	+	+	−	Moderate	Optimal	No	No	<90%
1	−	−	−a	−b	Good	Good	No	No	<90%
2	−	−	−	−b	Good	Good	Yes	No	<90%
3	−	Indb	−a	−b	Good	Good	No	No	>90%
4	+b	−	−a	−b	Good	Good	No	No	<90%
5	−a	−	−	−b	Good	Good	No	Yes	<90%
6	−	−	−	−b	Good	Good	No	No	<90%
7	−	−	−a	−b	Good	Good	No	No	>90%
8	−	−	−a	−	Good	Good	Yes	No	<90%
9	+b	−	−	Indb	Good	Good	Yes	No	<90%
10	−	−	+b	−b	Good	Optimal	No	No	<90%
11	−b	Ind	+b	−b	Good	Low	No	No	<90%
12	−	−	−b	−b	Good	Good	No	No	>90%
13	−	−a	−a	−b	Good	Good	Yes	No	<90%
14	−a	−	−b	−b	Good	Good	No	Yes	<90%
15	+b	−b	−b	−b	Good	Good	No	No	<90%
16	−	−	−a	−b	Good	Good	No	No	<90%
17	−	−	−a	−b	Good	Good	No	No	>90%
18	−	−	−	−b	Good	Good	Yes	No	<90%

Notes: (−): P=0–0.29; Ind: P=0.30–0.69; and (+): P=0.70–1.00. The OPS is a unique multivariate descriptor space in which the model is applicable. Assessment of this is needed to determine if the chemical structure being examined is within a model’s OPS. The probability results may be accepted with confidence, subjected to the results obtained from hypothesis testing.

a Outside of OPS but within OPS limits and

b outside of OPS and OPS limit.

Abbreviations: ADMET, absorption, distribution, metabolism, excretion, and toxicity; TOPKAT, toxicity prediction by komputer assisted technology; WOE, weight of evidence carcinogenicity; AM, Ames mutagenicity; DTP, developmental toxicity potential; AB, aerobic biodegradability; PPB, plasma protein binding; OPS, optimum prediction space.

Among these 18 compounds, two were available for acquisition from Enamine Ltd. These two compounds, namely compound 1 and compound 2, consequently underwent in vitro antimycobacterial screening using the resazurin microtiter assay (REMA). However, due to limitations in biosafety facilities and resources, we were constrained to use only Mtb H37Ra, an avirulent strain of the mycobacterium. Results of the REMA screening showed that compound 2 demonstrated excellent activity against H37Ra with percentage growth inhibition of ~94% at a concentration of 25 μM and 89% at 0.2 μM (Figure 3). Moreover, this remarkable action of compound 2 was sustained for a period of at least 3 days (Figure 4).

Figure 3 H37Ra growth inhibition (%) by compound 1 and compound 2 in comparison to the positive control, rifampicin.

Note: Error bars are the standard error of the mean.

Figure 4 H37Ra growth inhibition (%) by compound 2 at 25 μM and 0.2 μM after 24 hours, 48 hours, and 72 hours.

Note: Error bars are the standard error of the mean.

Ligand interaction analysis

In silico docking of meropenem, a known inhibitor into the active site of Ldt_Mt2 has shown characteristic protein–ligand interactions. A number of interactions with this complex was also observed in the Compound2-Ldt_Mt2 complex (Table 2).

Table 2 2D structures and ligand interactions of meropenem and compound 2 in Ldt_Mt2’s active site

Compound	2D structure	Binding energy (kcal/mol)	van der Waals	Conventional hydrogen bond	Carbon hydrogen bond	Attractive charge	Π–sulfur	Π–alkyl	Alkyl	Π–lone pair	Π–Π T-shaped	Π–sigma	Unfavorable donor–donor
Meropenem (docked)		148.884	Val322,	Trp340,	Thr320,	His336	Trp340,	His352,	His352,	–	–	–	–
			Met303,	Thr320	Arg319,		His336	Trp340	Trp340
			Tyr318,		Cys354,
			His352,		Gly332
			Thr350,
			Ser341,
			Ser351,
			Asn356,
			Lys282,
			Gly353,
			Val333,
			Phe334
Compound 2		164.695	Lys282,	Cys354,	–	–	Trp340,	–		Thr350	His352	Trp340	Lys282
			Thr350,	His336,			His352,
			Arg319,	Asn356,			Cys354
			Ser341,	Trp340,
			Ser351,	His352
			Thr320,
			Tyr318,
			Val322,
			Met303,
			His336,
			Asn356

Notes: Specific amino acids involved for each type of ligand interaction are enumerated for each compound. Calculated binding energies in kilocalorie per mole are also shown.

Abbreviations: 2D, two dimensional; Ldt_Mt2, l,d-transpeptidase 2.

Optimization of the protein model involved docking the meropenem ligand into the Ldt_Mt2 binding site. The docked Ldt_Mt2–meropenem complex, as shown in Figure 5, displayed hydrogen bonding with Cys354 (d₁=3.02 Å, d₂=3.74 Å) and Gly353 (d=3.38 Å). In Ldt_Mt2, His352, Gly353, and Cys354 form an oxyanion cavity, wherein the carbonyl group of the opened meropenem β-lactam ring is hydrogen bound to the mentioned residues. The docking data of meropenem depicted H-bonds with His352 and Cys354, which were proposed as two of the three key catalytic residues in Ldt_Mt2.^17,18,28 The H-bond with Cys354 could be the representation of the thioester bond between meropenem and Ldt_Mt2, since covalent bonding cannot be computed via Accelrys DS. In both originally bound and docked meropenem, a number of similar contacts were identified: van der Waals interaction with Val333; polar interaction with Gly353, Cys354, Tyr318, and His352; and side-chain H-bond with Cys354.

Figure 5 2D ligand interaction diagram of meropenem with Ldt_Mt2’s active site.

Abbreviations: 2D, two dimensional; Ldt_Mt2, l,d-transpeptidase 2.

Compound 2, which is the compound observed to be active in vitro against the TB strain, when docked within the Ldt_Mt2 protein, was predicted to have pi interaction with the sulfur of Cys354, which could be parallel to the meropenem–Cys354 thioester bond.^17,18 Of the catalytic triad, His336 and Cys354 were predicted to have H-bonding with meropenem, while other conventional H-bonds observed were with residues Trp340, His352, and Asn356. A Π–sulfur contact between meropenem and Trp340, a T-shaped Π–Π interaction with His352, and a Π–lone pair with Thr350 were observed. A number of van der Waals interactions were also predicted (Figure 6). A significant donor–donor clash was predicted between the meropenem pyrroline and Lys282. It is suggested that eliminating Lys282 or substituting the residue with another amino acid may further improve the stability of the complex.

Figure 6 2D ligand interaction diagram of compound 2 with the Ldt_Mt2 active site.

Abbreviations: 2D, two dimensional; Ldt_Mt2, l,d-transpeptidase 2.

Conclusion

Pharmacophore-based virtual screening, molecular docking, as well as in silico ADMET evaluation of compounds from the enamine database were performed in order to identify a new class of potential antitubercular lead compounds. Out of the 4.5 million compounds screened, 18 compounds were found to have better binding energies than meropenem and with satisfactory in silico ADMET properties. Out of the two of the 18 compounds that were tested in vitro, compound 2 was found to have an excellent bioactivity against Mtb H37Ra. Consequently, this lead compound may lead to a novel class of anti-TB drugs in the future.

Acknowledgments

This work was funded by the UP System Emerging Interdisciplinary Research Program (OVPAA-EIDR 12-001-121102).

Disclosure

The authors report no conflicts of interest in this work.