Next Article in Journal
Preparation of Biochar with Developed Mesoporous Structure from Poplar Leaf Activated by KHCO3 and Its Efficient Adsorption of Oxytetracycline Hydrochloride
Next Article in Special Issue
A QSAR Study for Antileishmanial 2-Phenyl-2,3-dihydrobenzofurans †
Previous Article in Journal
Characterization of Cell Death Induced by Imine Analogs of Trans-Resveratrol: Induction of Mitochondrial Dysfunction and Overproduction of Reactive Oxygen Species Leading to, or Not, Apoptosis without the Increase in the S-Phase of the Cell Cycle
Previous Article in Special Issue
eEF2K Inhibitor Design: The Progression of Exemplary Structure-Based Drug Design
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 7
10.3390/molecules28073186
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Targeting Toxoplasma gondii ME49 TgAPN2: A Bioinformatics Approach for Antiparasitic Drug Discovery

by
Ali Altharawi
Department of Pharmaceutical Chemistry, College of Pharmacy, Prince Sattam Bin Abdulaziz University, Al-Kharj 11942, Saudi Arabia
Molecules 2023, 28(7), 3186; https://doi.org/10.3390/molecules28073186
Submission received: 7 March 2023 / Revised: 29 March 2023 / Accepted: 29 March 2023 / Published: 3 April 2023
(This article belongs to the Special Issue Computational Approaches: Drug Discovery and Design in Medicinal Chemistry and Bioinformatics II)
Browse Figures
Versions Notes

Abstract

:
As fewer therapeutic options are available for treating toxoplasmosis, newer antiparasitic drugs that can block TgAPN2 M1 aminopeptidase are of significant value. Herein, we employed several computer-aided drug-design approaches with the objective of identifying drug molecules from the Asinex library with stable conformation and binding energy scores. By a structure-based virtual screening process, three molecules—LAS_52160953, LAS_51177972, and LAS_52506311—were identified as promising candidates, with binding affinity scores of −8.6 kcal/mol, −8.5 kcal/mol, and −8.3 kcal/mol, respectively. The compounds produced balanced interacting networks of hydrophilic and hydrophobic interactions, vital for holding the compounds at the docked cavity and stable binding conformation. The docked compound complexes with TgAPN2 were further subjected to molecular dynamic simulations that revealed mean RMSD for the LAS_52160953 complex of 1.45 Å), LAS_51177972 complex 1.02 Å, and LAS_52506311 complex 1.087 Å. Another round of binding free energy validation by MM-GBSA/MM-PBSA was done to confirm docking and simulation findings. The analysis predicted average MM-GBSA value of <−36 kcal/mol and <−35 kcal/mol by MM-PBSA. The compounds were further classified as appropriate candidates to be used as drug-like molecules and showed favorable pharmacokinetics. The shortlisted compounds showed promising biological potency against the TgAPN2 enzyme and may be used in experimental validation. They may also serve as parent structures to design novel derivatives with enhanced biological potency.

1. Introduction

The M1 aminopeptidases (also termed aminopeptidase N) are found in living systems and are anchored in cell membrane [1]. These enzymes are vital from a functionality perspective and are important in cell growth, development, maintenance, and defense [2]. They are documented as potential candidates for novel drugs and regarded as potential immunoprotected targets [3]. In the Plasmodium falciparum malaria parasite, a single M1 aminopeptidase enzyme has been reported, classified as PfA-M1, that plays a prime role in hemoglobin digestion to ensure pathogen growth and development in the blood stage of the parasite [4]. Targeted inhibition of PfA-M1 has been unveiled to stop P. falciparum inhibition in both laboratory and animal murine models [5]. Toxoplasma gondii is a causative agent of toxoplasmosis, which can cause serious complications in pregnant women [2,6,7,8,9]. The congenital pathogen infection results in mental retardation, miscarriage, and hearing and vision problems [6]. Toxoplasmosis management relies on chemotherapy comprising pyrimethamine with either clindamycin or sulfadiazine [10]. Drugs that can treat both active and latent infections are scarce, and thus new biomolecular targets and their inhibitors are needed for drug discovery [11].
The M1 aminopeptidases are of three types in the T. gondii ME49 strain—TgAPN1, TgAPN2, and TgAPN3—and encoded by genes present on different chromosomes [2]. TgAPN1 is experimentally active and shown as immunogenic. TgAPN3 is a metalloexopeptidase enzyme and a functional peptide, despite the N-terminal transmembrane anchor. TgAPN2 has localization within the parasite cytosol and is expressed at all stages of the parasite life cycle. The biological function of TgAPN2 is still not clear; however, it is predicted to play a role in tachyzoite growth. Its upregulation in bradyzoites also makes it an important target for drug development [2].
In contrast to conventional drug discovery, computer-aided drug design (CADD) is considered attractive, as it is more cost-effective and consumes less time and little human resources [12,13,14]. CADD is a multidisciplinary approach and involves structure-based and ligand-based virtual screening of drug libraries against given biological targets [15,16]. The aim of the study was to identify the best binding molecules of T. gondii ME49 TgAPN2 that show stable binding mode associated with lowest binding energy score. The aim was split into several objectives. First, screening of the Asinex drug library was done by using blind screening to identify novel binding sites for inhibitors [17]. The docking analysis includes both identification of compound binding mode as well as binding affinity score [18,19]. As the docking and virtual screening calculations are associated with high false-positive rates, the study was extended by molecular dynamic simulation analysis, which aids in validating hit compound binding modes and interaction stability along the length of simulation time [20]. Molecular dynamic simulation also helps in pointing to protein regions that are more flexible and stable in the presence of compounds and residues that are vital in compound binding and biomolecule catalytic mechanisms. Dynamic simulation analysis is key in determining the binding and chemical interactions of lead molecules to the receptors and deciphering conformation stability during the simulation [21,22,23,24,25,26,27,28,29,30,31,32]. Binding free energy analysis was conducted to revalidate the docking and simulation findings [20,21]. Binding free energy methods rely on simulation trajectories and play a key role in compound structure optimization for lead discovery [25,26,27,28,29,30,31,32]. These methods have gained considerable appreciation in recent years, as they have proved vital in estimation of relative binding free energies [33,34]. Compared to molecular docking- or virtual screening-based approaches, these methods are more significant in predicting compound binding [35,36,37].The outcomes of this study may hasten antiparasitic drug development for controlling T. gondii ME49 diseases. Furthermore, new chemical structures may be designed from the compound’s parent structure in order to achieve good biological potency and a safer pharmacokinetic profile.

2. Results and Discussion

2.1. Structure-Based Virtual Screening

Virtual screening was conducted to identify compounds that show stable binding at the TgAPN2 binding pocket. This was achieved by screening the Asinex library, which contains diverse chemical scaffolds and ready-to use-products for experimental investigation. The virtual screening identified 10 compounds as promising leads (Table 1). Among these, only three were complexed with TgAPN2 for docking conformation and interaction analysis.

2.2. Top Three Compounds’ Docking Analysis

The top three compounds—LAS_52160953, LAS_51177972, and LAS_52506311—were chosen based on the lowest binding energy value in the virtual screening process, with binding energy values of −8.6 kcal/mol, −8.5 kcal/mol and −8.3 kcal/mol, respectively. Chemically, the LAS_52160953, LAS_51177972 and LAS_52506311 are 1-(3-(3-(2-methyl-2,3-dihydro-1H-imidazol-1-yl)propoxy)benzyl)-4-((3,4,5-trimethylphenoxy)methyl)piperidin-4-ol, 5-carboxy-2-(2-(2,5-dimethoxyphenyl)-1H-pyrrol-1-yl)-4-((3,4,5-trimethylphenoxy)methyl)pyrimidine-1-ium and 3-(phenoxymethyl)-1-(4-(2-(piperidin-1-yl)ethoxy)benzyl)piperidin-3-ol. All three compounds were noticed to achieve deep binding inside the active pocket and produced rich hydrophilic and hydrophobic contacts (Figure 1). The LAS_52160953 1-methyl-4-((3,4,5-trimethylphenoxy)methyl)piperidin-4-ol chemical moiety was seen docked deep inside the pocket facing the pocket bottom and formed a hydrogen bond contact with Tyr914 at distance of 2.3 Å. The opposite 2-methyl-1-(3-phenoxypropyl)-2,3-dihydro-1H-imidazole ring is placed at a pocket-out position and formed multiple weak hydrophobic contact. The compound interacts with several important residues, such as Asp920, Arg889, Gly1369, Leu1372, Thr915, Arg1366, Glu836 and Gly799. The LAS_51177972 was among the most stable compound due to multiple hydrogen bonds contact with the enzyme. The compound 5-methoxy-1,2,3-trimethylbenzene ring tend to interact with enzyme base residues while the opposite 5-carboxy-2-(2-(2,5-dimethoxyphenyl)-1H-pyrrol-1-yl)pyrimidine-1-ium interact with the active pocket walls. The compounds formed hydrogen bonds with Arg1366, Tyr919, Ala916 and Asp920 at distance of 2.3 Å, 3.1 Å, 2.6 Å, and 2.2 Å, respectively. LAS_52506311 3-(phenoxymethyl)piperidin-3-ol is involved in hydrogen bond contact Asp865 while the 1-(2-phenoxyethyl)piperidine ring interactions are dominated by van der Waals involving residues of Asp920, Val862, Gly836, Phe885, Thr915, Arg1366, Arg889, Ala916, and Leu1372 etc. The control compound 1 binding energy score was −8.52 kcal/mol. The control mainly formed hydrogen bonds with Asn831, His835 and Glu924. The van der Waals contacts were Val883, Arg879, Phe885, Gln866, Lys878, Leu886, Val832, Val862, and Tyr919 (Figure S1).

2.3. Molecular Dynamic Simulations

Molecular dynamic simulations were conducted to investigate the physical movements of TgAPN2 in the presence of lead molecules. This helped in understanding the long-term compounds conformational stability with the enzyme and shed light on key molecular features of both receptor and compounds vital in intermolecular interactions. The analysis includes root mean square deviation (RMSD) [38,39], root mean square fluctuation (RMSF) [31,39] and radius of gyration (RoG) [4,5] (Figure 2). The analyses were done by mean of carbon alpha atoms. The RMSD describes the structure variation by measuring distance among the superimposed docked intermolecular conformation over the initial docked conformation. The RMSD demonstrated all the docked systems stable with mean RMSD value as LAS_52160953 complex (1.45 Å), LAS_51177972 complex (1.02 Å) and LAS_52506311 complex (1.087 Å). Among the systems, the first had more structural variations than the last two; however, it was still in a very stable range. All these values determine the complexes as highly stable from the docked conformation perspective. It further emphasizes that the compounds’ docked conformation with the enzyme is solid and no freedom of energy is available in the systems to detach the ligands. The control compound 1 RMSD is given in Figure S2, which shows mean RMSD of 2.4 Å. The reflects that the control docked complex with the enzyme showed more structure deviations than the filtered complexes. Previous studies reported that RMSD is important to determine structural stability of compounds at the receptor active pocket [39,40,41,42]. This was further validated by RMSF assay, which determines residue-based stability. RMSF analysis is important to underline the key regions of enzymes that play a vital role in compound recognition and long-term binding stability. RMSF, as can be seen in Figure 2B, evaluated that most of the receptor residues lie within the stable range in the presence of compounds. However, the N-terminal and C-terminal of the enzyme are a bit more flexible due to loops that are naturally flexible. The flexibility provides structure adaption to the enzyme during the catalytic process, thereby creating enough room for compounds to attain stable energy. Next, the compact and relaxed 3D conformation of TgAPN2 was investigated using RoG [32,43,44]. This analysis is important to understand the compactness and relaxed nature of the enzyme in the presence of compounds during simulation. Higher RoG determines loose packing of the enzyme’s secondary structure elements which may results in ligand detachment. The RoG depicted that the systems are compact with no major variations seen. This analysis complements the RMSD in gaining confidence in the system’s intermolecular docked conformation. The RoG value of the system ranged between 43 Å and 45 Å. These values are within the stable range and depict the stable compound binding with the enzyme and higher complex intermolecular stability.

2.4. Binding Free Energy Calculations

The binding free energy calculation by MM-GBSA/MM-PBSA was carried out in order to reassess the docked affinity of compounds for Toxoplasma gondii ME49 TgAPN2. The details of binding free energies are given in Table 2. All the complexes revealed stable energies, with an average value <−36 kcal/mol. The overall energy was dominated by van der Waals forces, which illustrated the complexes were stabilized mainly by weak hydrophobic energy. The TgAPN2- LAS_52160953 complex, LAS_51177972 complex and LAS_52506311 complex had net MM-GBSA binding energy scores of −39.45 kcal/mol, −36.19 kcal/mol, and −39.4 kcal/mol, respectively. High contributions were seen from van der Waals energy: −36.89 kcal/mol for LAS_52160953 complex, −33.60 kcal/mol for LAS_51177972 complex, and −36.02 kcal/mol for LAS_52506311 complex. The electrostatic energy contribution, though, was less contributing than the van der Waals energy, but still played a favorable role in complex stabilization. There were high contributions for LAS_52160953 and LAS_52506311: −14.23 kcal/mol and −14.85 kcal/mol, respectively. The net gas phase energy involved insignificant energy contribution from internal energy. The solvation energy contribution for each complex in MM-GBSA was seen as unfavorable. Similarly, on MM-PBSA analysis, the complexes also revealed stable energies. The net energy MM-PBSA values for TgAPN2-LAS_52160953 complex, LAS_51177972 complex and LAS_52506311 complex were −38.57 kcal/mol, −35.23 kcal/mol and −37.12 kcal/mol, respectively. These energy values were more stable for complexes than reported by MM-GBSA. Binding entropy values of LAS_52160953, LAS_51177972 and LAS_52506311 were 16.14 kcal/mol, 17.31 kcal/mol and 16.78 kcal/mol, respectively.

2.5. WaterSwap Analysis

WaterSwap analysis was done to gain more confidence in the intermolecular affinity of complexes. The WaterSwap method is more sophisticated than the MM-GBSA/MM-PBSA method and considers the water molecules’ role in bridging ligand to receptor residues. This method has been successfully used in different studies and provided promising results in term of predicting compound binding affinity for the targeted biological macromolecule [27,45]. Three algorithms are used in WaterSwap: thermodynamic integration (TI), free energy perturbation (FEP), and Bennett’s. The three algorithms’ energy value for each complex is illustrated in Figure 3. All the three systems were confirmed to get maximum stable energy, depicting a strong interacting network.

2.6. SwissADME Analysis

In silico prediction of compounds ADMET properties is vital in present drug-discovery pipelines, as it saves cost and time associated with drug failure in clinical trials. The SwissADME details of LAS_52160953, LAS_51177972 and LAS_52506311 are given in Table 3. All the three compounds showed favorable adsorption, distribution, metabolism and excretion [46]. The compounds also revealed good topological polar surface area (TPSA) to allow efficient drug adsorption across the cell membrane [47]. Similarly, the compounds have good gastrointestinal absorption, and thus a high concentration of drugs can reach the site of action. Additionally, the compounds are drug-like, as classified by most rules, specifically Lipinski’s rule of five [48], Veber [49], and Egan [50]. Lipinski’s rule of five is a famous drug rule and describes that a biologically active drug molecule should possess several physical and chemical properties to be orally active. These properties include suitable molecular weight (<500 Da), hydrogen bond donors (<5), hydrogen bond acceptors (<10), and partition coefficient less than 5. Additionally, the compounds have no pan-assay interference compounds (PAINS) alerts and thus cannot bind to multiple targets and will specifically bind to a single biomolecule [51]. This prevents cross-reaction of the compounds with different proteins/enzymes when administered to the host and may also reduce unwanted and adverse effects. The compounds are synthetically feasible to synthesize. The easy synthesis of the compounds allows their use in experimental testing to disclose their real biological potency by blocking the biological function of the TgAPN enzyme.

3. Materials and Methods

3.1. Preparation of Protein Structure and Asinex Library

The crystallographic structure of T. gondii ME49 TgAPN2 was retrieved from the Protein Data Bank using PDB ID 6OIU [2,52]. The PDB structure was selected due to its recent submission to the database and good resolution quality compared to other available structures. The experimental data snapshot of the structure is as follows: method X-ray diffraction, resolution 2.20 Å, R-value free 0.236, R-value work 0.190, and R-value observed 0.193. The global stoichiometry of the structure is monomer. The native ligands were removed from the crystal structures along with water molecules. The missing amino acids were added using the CHARMM-GUI PDB manipulator [53]. The protonation of the structure was carried out at pH 7.4 using the H++ web server [54]. The structure was energy-minimized using steepest descent and conjugate gradient algorithms in UCSF Chimera v1.16 to eliminate steric clashes [55]. After energy optimization, the structure was saved in the PDB for additional analysis. The drug library used in this study was Asinex (https://www.asinex.com/ accessed on 15 November 2022). This library contains 575,302 compounds (updated on February 2023). It provides a cost-effective, diverse, drug-like chemical space. Most of the compounds in the library have a high degree of drug-likeness. The library was imported to PyRx 0.8 software [56]. The library was then energy-minimized using MM2 force field and converted to pdbqt format to make the compounds ready for virtual screening [57].

3.2. Virtual Screening

Virtual screening is a regularly used technique to screen drug libraries against any given enzyme/protein active pocket with the goal of identifying molecules that show best fitting with the receptor enzyme [58]. This technique has been used in different studies and provided lead molecules that produced good biological potency against the receptor [59]. The virtual screening process was initiated by docking all the Asinex library using the PyRx 0.8 AutoDock Vina program. The screening was performed against catalytic domain II by specifying the residues His835, His839 and Glu858 [2].The grid box dimensions were 25 Å on the XYZ axes. The grid box was centered on the X-axis = 90.81 Å, Y-axis = 25.40 Å, and Z-axis = 126.93 Å. The grid box in virtual screening was set as such to provide flexibility to the library compounds for binding the cavity across the T. gondii ME49 TgAPN2, where they show high conformational stability. The number of iterations set for each compound in docking was 100. The docking process was validated by redocking of cocrystalized ligand to the same site reported in the crystal structure. The superimposition of crystalized ligand and docked conformation revealed a root mean square deviation (RMSD) value of 0.27 Å, illustrating docking procedure accuracy [38]. For a positive control, compound 1 from Marijanovic et al., 2020 was used [2].

3.3. Molecular Dynamic Simulation

Molecular dynamic simulation is a computer-based approach to study physical movements of atoms/molecules in a fixed period of time [20]. The molecules’ movements are solved by Newton’s law of motion, and behavior is plotted in the form of a graph [60]. The top complexes (LAS_52160953, LAS_51177972 and LAS_52506311) were subjected to molecular dynamic simulations in explicit water. The starting atomic coordinates used were obtained from molecular docking and processed through the AMBER20 simulation package [61]. Amber ff14sb was employed for M1 aminopeptidase TgAPN2, while generalized Amber force field (GAFF) was applied for nonstandard residues as well as the TIP3P water model [62,63,64,65]. Amber ff14sb was used as it has improved potential for side chain torsion potentials. The compounds’ partial charges were added using the AM1-BCC method via antechamber program [66]. Each complex was placed at the center of a dodecahedral box with padding distance of 12 Å, followed by charge neutralization by adding an appropriate number of counterions. The number of sodium ions added were 12, 11 and 13 for LAS_52160953 complex, LAS_51177972 complex and LAS_52506311 complex, respectively. The complexes were next subjected to energy minimization for 5000 steps of the steepest descent algorithm in order to remove irregular geometry and steric clashes. During energy minimization, hydrogen atoms were minimized first for 500 cycles. This was followed by entire systems’ atom minimization with restraint of 10 kcal/mol Å2. Then, nonheavy atoms were lastly energy-minimized for 500 rounds in the presence of 200 kcal/mol Å2 on the remaining system. Afterward, equilibration was done for 100 ps under canonical ensemble. Then, another equilibration was performed for each system under isothermal isobaric ensemble in the presence of 1 bar pressure. Heating of each system was done gradually and achieved a constant temperature of 310 K. The production run was performed for 200 ns in periodic boundary conditions. The production run was facilitated in NVT (constant temperature, constant volume) ensemble and in the presence of Berendsen algorithm. The SHAKE algorithm and Langevin were considered for maintain restraint on hydrogen bonds and constant temperature, respectively [67,68]. The particle-mesh Ewald method was employed for treating long-range interactions. The cutoff value used for long-range interactions was 1.0 nm. The CPPTRAJ module of AMBER was used for structure stability analysis of systems [69]. Plots were generated using XMGRACE v 5.1 [70] and snapshots of simulation were visualized using VMD v 1.93 [71].

3.4. MM/PBSA Binding Free Energy Calculations

The docking programs offer simplified scoring functions for predicting ligand binding affinity for receptors. Therefore, they suffer from several limitations in accurate prediction of binding energy. Validation of the binding affinity of studied compounds was accomplished using molecular mechanics Poisson–Boltzmann surface area (MM/PBSA) and molecular mechanics generalized Born surface area (MM/GBSA) methods as implemented in AMBER MMPBSA.py [29,30,72,73,74]. In the process, 1000 snapshots were extracted from the simulation trajectories. The MM/PBSA net binding energy was determined using the following equation:
ΔGbind = Gcomplex − (Greceptor + Gligand)
During the analysis, the nonpolar energies were estimated using linear combination of pairwise overlaps with water probe radius of 1.4 Å. The β and γ values of nonpolar solvation energy employed were 0 kcal/mol in MMGBSA and 0.092 kcal/mol in MMPBSA. The dielectric constant value in MM-GBSA was set to 1.0. The entropy energy contribution in complex formation was determined using a bash script described by Duan et al., 2016 [75]. The MM-GBSA and MM-PBSA binding energies were reconfirmed by WaterSwap, which is considered more sophisticated in terms of determining water molecules’ contribution in ligand binding to the receptor [76,77].

3.5. Computational Pharmacokinetic Analysis

Computational studies on compounds’ physicochemical properties, lipophilicity, water solubility, pharmacokinetics, drug-likeness, and medicinal chemistry were done using SwissADME webserver [78].

4. Conclusions

In this work, three high-affinity binders were identified—LAS_52160953, LAS_51177972, and LAS_52506311—against TgAPN2. The molecules achieved stable binding at the enzyme active pocket and produced a stable interacting network of hydrophilic and hydrophobic contacts. The docked complexes also demonstrated stable dynamics with no major structure variations. The compound binding of compounds was confirmed by different binding free energy methods. The binding energy of intermolecular complexes was dominated by van der Waals energy. Additionally, the electrostatic energy had a vital contribution in making the compounds’ stable binding mode. Though results of the study are promising, revalidation of biological potency can be ensured by subjecting the compounds to in vitro and in vivo studies. Further, the lead compounds’ structures might be used in designing novel promising derivatives that have fewer structure and physicochemical limitations.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/molecules28073186/s1, Figure S1: Docked conformation and interactions of control compound 1 with TgAPN2 (shown by secondary structure cartoon representation); Figure S2: RMSD of control compound 1.

Funding

This study was supported via funding from Prince Sattam Bin Abdulaziz University project number (PSAU/2023/R/1444).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available within the article.

Acknowledgments

This research was financially supported by Prince Sattam bin Abdulaziz University (PSAU) under the umbrella of the Deanship of Scientific Research. I am extremely grateful to PSAU for their kind and generous support.

Conflicts of Interest

The author declares no conflict of interest.

Sample Availability

Samples of the compounds are not available from the author.

References

  1. Chen, L.; Lin, Y.-L.; Peng, G.; Li, F. Structural basis for multifunctional roles of mammalian aminopeptidase N. Proc. Natl. Acad. Sci. USA 2012, 109, 17966–17971. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Marijanovic, E.M.; Swiderska, K.W.; Andersen, J.; Aschenbrenner, J.C.; Webb, C.T.; Drag, M.; Drinkwater, N.; McGowan, S. X-ray crystal structure and specificity of the Toxoplasma gondii ME49 TgAPN2. Biochem. J. 2020, 477, 3819–3832. [Google Scholar] [CrossRef] [PubMed]
  3. Drinkwater, N.; Lee, J.; Yang, W.; Malcolm, T.R.; McGowan, S. M1 Aminopeptidases as Drug Targets: Broad Applications or Therapeutic Niche? FEBS J. 2017, 284, 1473–1488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Mathew, R.; Wunderlich, J.; Thivierge, K.; Cwiklinski, K.; Dumont, C.; Tilley, L.; Rohrbach, P.; Dalton, J.P. Biochemical and cellular characterisation of the Plasmodium falciparum M1 alanyl aminopeptidase (PfM1AAP) and M17 leucyl aminopeptidase (PfM17LAP). Sci. Rep. 2021, 11, 2854. [Google Scholar] [CrossRef] [PubMed]
  5. Bounaadja, L.; Schmitt, M.; Albrecht, S.; Mouray, E.; Tarnus, C.; Florent, I. Selective inhibition of PfA-M1, over PfA-M17, by an amino-benzosuberone derivative blocks malaria parasites development in vitro and in vivo. Malar. J. 2017, 16, 382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Hill, D.; Dubey, J. Toxoplasma gondii: Transmission, diagnosis and prevention. Clin. Microbiol. Infect. 2002, 8, 634–640. [Google Scholar] [CrossRef] [Green Version]
  7. Liu, Q.; Wang, Z.-D.; Huang, S.-Y.; Zhu, X.-Q. Diagnosis of toxoplasmosis and typing of Toxoplasma gondii. Parasites Vectors 2015, 8, 292. [Google Scholar] [CrossRef] [Green Version]
  8. Ali, A.; Naz, A.; Soares, S.C.; Bakhtiar, M.; Tiwari, S.; Hassan, S.S.; Hanan, F.; Ramos, R.; Pereira, U.; Barh, D.; et al. Pan-Genome Analysis of Human Gastric Pathogen, H. Pylori: Comparative Genomics and Pathogenomics Approaches to Identify Regions Associated with Pathogenicity and Prediction of Potential Core Therapeutic Targets. Biomed. Res. Int. 2015, 2015. [Google Scholar]
  9. McGovern, O.L.; Rivera-Cuevas, Y.; Carruthers, V.B. Emerging Mechanisms of Endocytosis in Toxoplasma gondii. Life 2021, 11, 84. [Google Scholar] [CrossRef]
  10. Milne, G.; Webster, J.P.; Walker, M. Toxoplasma Gondii: An Underestimated Threat? Trends Parasitol. 2020, 36, 959–969. [Google Scholar] [CrossRef]
  11. Montazeri, M.; Mehrzadi, S.; Sharif, M.; Sarvi, S.; Tanzifi, A.; Aghayan, S.A.; Daryani, A. Drug Resistance in Toxoplasma gondii. Front. Microbiol. 2018, 9, 2587. [Google Scholar] [CrossRef] [Green Version]
  12. Macalino, S.J.Y.; Gosu, V.; Hong, S.; Choi, S. Role of computer-aided drug design in modern drug discovery. Arch. Pharm. Res. 2015, 38, 1686–1701. [Google Scholar] [CrossRef] [PubMed]
  13. Van Drie, J.H. Computer-aided drug design: The next 20 years. J. Comput. Aided Mol. Des. 2007, 21, 591–601. [Google Scholar] [CrossRef] [PubMed]
  14. Yu, W.; MacKerell, A.D. Computer-Aided Drug Design Methods. In Antibiotics; Springer: Berlin/Heidelberg, Germany, 2017; pp. 85–106. [Google Scholar]
  15. Alamri, M.A.; Tariq, M.H.; Tahir ul Qamar, M.; Alabbas, A.B.; Alqahtani, S.M.; Ahmad, S. Discovery of Potential Phytochemicals as Inhibitors of TcdB, a Major Virulence Factors of Clostridioides Difficile. J. Biomol. Struct. Dyn. 2023, 1–9. [Google Scholar] [CrossRef]
  16. Altharawi, A.; Riadi, Y.; Qamar, M.T.U. An in Silico Quest for Next-Generation Antimalarial Drugs by Targeting Plasmodium Falciparum Hexose Transporter Protein: A Multi-Pronged Approach. J. Biomol. Struct. Dyn. 2023, 1–10. [Google Scholar] [CrossRef] [PubMed]
  17. Bechelane-Maia, E.H.; Assis, L.C.; Alves de Oliveira, T.; Marques da Silva, A.; Gutterres Taranto, A. Structure-based virtual screening: From classical to artificial intelligence. Front. Chem. 2020, 8, 343. [Google Scholar] [CrossRef]
  18. Meng, X.-Y.; Zhang, H.-X.; Mezei, M.; Cui, M. Molecular Docking: A Powerful Approach for Structure-Based Drug Discovery. Curr. Comput. Aided-Drug Des. 2011, 7, 146–157. [Google Scholar] [CrossRef]
  19. Huey, R.; Morris, G.M. Using AutoDock 4 with AutoDocktools: A Tutorial. Scripps Res. Inst. USA 2008, 54, 56. [Google Scholar]
  20. Hansson, T.; Oostenbrink, C.; van Gunsteren, W. Molecular Dynamics Simulations. Curr. Opin. Struct. Biol. 2002, 12, 190–196. [Google Scholar] [CrossRef]
  21. Qamar, M.T.U.; Ahmad, S.; Fatima, I.; Ahmad, F.; Shahid, F.; Naz, A.; Abbasi, S.W.; Khan, A.; Mirza, M.U.; Ashfaq, U.A.; et al. Designing multi-epitope vaccine against Staphylococcus aureus by employing subtractive proteomics, reverse vaccinology and immuno-informatics approaches. Comput. Biol. Med. 2021, 132, 104389. [Google Scholar] [CrossRef]
  22. Ali, S.; Ali, S.; Javed, S.O.; Shoukat, S.; Ahmad, S.; Ali, S.S.; Hussain, Z.; Waseem, M.; Rizwan, M.; Suleman, M. Proteome wide vaccine targets prioritization and designing of antigenic vaccine candidate to trigger the host immune response against the Mycoplasma genitalium infection. Microb. Pathog. 2021, 152, 104771. [Google Scholar] [CrossRef] [PubMed]
  23. Durrant, J.D.; McCammon, J.A. Molecular dynamics simulations and drug discovery. BMC Biol. 2011, 9, 71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Ahmad, S.; Raza, S.; Uddin, R.; Azam, S.S. Comparative subtractive proteomics based ranking for antibiotic targets against the dirtiest superbug: Acinetobacter baumannii. J. Mol. Graph. Model. 2018, 82, 74–92. [Google Scholar] [CrossRef]
  25. Sanober, G.; Ahmad, S.; Azam, S.S. Identification of plausible drug targets by investigating the druggable genome of MDR Staphylococcus epidermidis. Gene Rep. 2017, 7, 147–153. [Google Scholar] [CrossRef]
  26. Muneer, I.; Ahmad, S.; Naz, A.; Abbasi, S.W.; Alblihy, A.; Aloliqi, A.A.; Alkhayl, F.F.; Alrumaihi, F.; Ahmad, S.; El Bakri, Y. Discovery of Novel Inhibitors from Medicinal Plants for V-Domain Ig Suppressor of T-Cell Activation (VISTA). Front. Mol. Biosci. 2021, 8, 951. [Google Scholar] [CrossRef]
  27. Ahmad, S.; Waheed, Y.; Ismail, S.; Abbasi, S.W.; Najmi, M.H. A computational study to disclose potential drugs and vaccine ensemble for COVID-19 conundrum. J. Mol. Liq. 2020, 324, 114734. [Google Scholar] [CrossRef]
  28. Ehsan, N.; Ahmad, S.; Navid, A.; Azam, S.S. Identification of potential antibiotic targets in the proteome of multi-drug resistant Proteus mirabilis. Meta Gene 2018, 18, 167–173. [Google Scholar] [CrossRef]
  29. Genheden, S.; Ryde, U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 2015, 10, 449–461. [Google Scholar] [CrossRef]
  30. Wang, E.; Sun, H.; Wang, J.; Wang, Z.; Liu, H.; Zhang, J.Z.H.; Hou, T. End-Point Binding Free Energy Calculation with MM/PBSA and MM/GBSA: Strategies and Applications in Drug Design. Chem. Rev. 2019, 119, 9478–9508. [Google Scholar] [CrossRef]
  31. El Bakri, Y.; Anouar, E.H.; Ahmad, S.; Nassar, A.A.; Taha, M.L.; Mague, J.T.; El Ghayati, L.; Essassi, E.M. Synthesis and Identification of Novel Potential Molecules Against COVID-19 Main Protease Through Structure-Guided Virtual Screening Approach. Appl. Biochem. Biotechnol. 2021, 193, 3602–3623. [Google Scholar] [CrossRef]
  32. Javed, N.; Ahmad, S.; Raza, S.; Azam, S.S. Subtractive Proteomics Supported with Rational Drug Design Approach Revealed ZINC23121280 as a Potent Lead Inhibitory Molecule for Multi-Drug Resistant Francisella Tularensis: Drug Designing for Multidrug-Resistant Francisella Tularensis. Proc. Pak. Acad. Sci. B Life Environ. Sci. 2021, 58, 1–42. [Google Scholar] [CrossRef]
  33. Alamri, M.A.; Ahmad, S.; Alqahtani, S.M.; Irfan, M.; Alabbas, A.B.; Qamar, M.T.U. Screening of marine natural products for potential inhibitors targeting biotin biosynthesis pathway in Mycobacterium tuberculosis. J. Biomol. Struct. Dyn. 2022, 1–9. [Google Scholar] [CrossRef] [PubMed]
  34. Alamri, M.A.; Tahir ul Qamar, M.; Alabbas, A.B.; Alqahtani, S.M.; Alossaimi, M.A.; Azam, S.; Hashmi, M.H.; Rajoka, M.S.R. Molecular and Structural Analysis of Specific Mutations from Saudi Isolates of SARS-CoV-2 RNA-Dependent RNA Polymerase and Their Implications on Protein Structure and Drug–Protein Binding. Molecules 2022, 27, 6475. [Google Scholar] [CrossRef] [PubMed]
  35. Tahir ul Qamar, M.; Zhu, X.-T.; Chen, L.-L.; Alhussain, L.; Alshiekheid, M.A.; Theyab, A.; Algahtani, M. Target-Specific Machine Learning Scoring Function Improved Structure-Based Virtual Screening Performance for SARS-CoV-2 Drugs Development. Int. J. Mol. Sci. 2022, 23, 11003. [Google Scholar] [CrossRef]
  36. Alamri, M.A.; Mirza, M.U.; Adeel, M.M.; Ashfaq, U.A.; Qamar, M.T.U.; Shahid, F.; Ahmad, S.; Alatawi, E.A.; Albalawi, G.M.; Allemailem, K.S. Structural Elucidation of Rift Valley Fever Virus L Protein towards the Discovery of Its Potential Inhibitors. Pharmaceuticals 2022, 15, 659. [Google Scholar] [CrossRef]
  37. Ahmad, F.; Albutti, A.; Tariq, M.H.; Din, G.; Qamar, M.T.U.; Ahmad, S. Discovery of Potential Antiviral Compounds against Hendra Virus by Targeting Its Receptor-Binding Protein (G) Using Computational Approaches. Molecules 2022, 27, 554. [Google Scholar] [CrossRef]
  38. Maiorov, V.N.; Crippen, G.M. Significance of Root-Mean-Square Deviation in Comparing Three-dimensional Structures of Globular Proteins. J. Mol. Biol. 1994, 235, 625–634. [Google Scholar] [CrossRef] [Green Version]
  39. Ahmad, S.; Raza, S.; Uddin, R.; Azam, S.S. Binding mode analysis, dynamic simulation and binding free energy calculations of the MurF ligase from Acinetobacter baumannii. J. Mol. Graph. Model. 2017, 77, 72–85. [Google Scholar] [CrossRef] [PubMed]
  40. Lobanov, M.Y.; Bogatyreva, N.S.; Galzitskaya, O.V. Radius of Gyration as an Indicator of Protein Structure Compactness. Mol. Biol. 2008, 42, 623–628. [Google Scholar] [CrossRef]
  41. Wahedi, H.M.; Ahmad, S.; Abbasi, S.W. Stilbene-based natural compounds as promising drug candidates against COVID-19. J. Biomol. Struct. Dyn. 2021, 39, 3225–3234. [Google Scholar] [CrossRef]
  42. Humayun, F.; Khan, A.; Ahmad, S.; Yuchen, W.; Wei, G.; Nizam-Uddin, N.; Hussain, Z.; Khan, W.; Zaman, N.; Rizwan, M.; et al. Abrogation of SARS-CoV-2 Interaction with Host (NRP1) Neuropilin-1 Receptor through High-Affinity Marine Natural Compounds to Curtail the Infectivity: A Structural-Dynamics Data. Comput. Biol. Med. 2021, 141, 104714. [Google Scholar] [CrossRef] [PubMed]
  43. Khan, S.; Irfan, M.; Hameed, A.R.; Ullah, A.; Abideen, S.A.; Ahmad, S.; Haq, M.U.; El Bakri, Y.; Al-Harbi, A.I.; Ali, M. Vaccinomics to design a multi-epitope-based vaccine against monkeypox virus using surface-associated proteins. J. Biomol. Struct. Dyn. 2022, 1–10. [Google Scholar] [CrossRef] [PubMed]
  44. Ahmad, S.; Ranaghan, K.E.; Azam, S.S. Combating tigecycline resistant Acinetobacter baumannii: A leap forward towards multi-epitope based vaccine discovery. Eur. J. Pharm. Sci. 2019, 132, 1–17. [Google Scholar] [CrossRef] [Green Version]
  45. Ain, Q.U.; Ahmad, S.; Azam, S.S. Subtractive proteomics and immunoinformatics revealed novel B-cell derived T-cell epitopes against Yersinia enterocolitica: An etiological agent of Yersiniosis. Microb. Pathog. 2018, 125, 336–348. [Google Scholar] [CrossRef]
  46. Van De Waterbeemd, H.; Gifford, E. ADMET in Silico Modelling: Towards Prediction Paradise? Nat. Rev. Drug Discov. 2003, 2, 192–204. [Google Scholar] [CrossRef]
  47. Assenov, Y.; Ramírez, F.; Schelhorn, S.-E.; Lengauer, T.; Albrecht, M. Computing topological parameters of biological networks. Bioinformatics 2007, 24, 282–284. [Google Scholar] [CrossRef] [Green Version]
  48. Lipinski, C.A. Lead- and drug-like compounds: The rule-of-five revolution. Drug Discov. Today Technol. 2004, 1, 337–341. [Google Scholar] [CrossRef] [PubMed]
  49. Veber, D.F.; Johnson, S.R.; Cheng, H.-Y.; Smith, B.R.; Ward, K.W.; Kopple, K.D. Molecular Properties That Influence the Oral Bioavailability of Drug Candidates. J. Med. Chem. 2002, 45, 2615–2623. [Google Scholar] [CrossRef]
  50. Egan, W.J.; Merz, K.M.; Baldwin, J.J. Prediction of Drug Absorption Using Multivariate Statistics. J. Med. Chem. 2000, 43, 3867–3877. [Google Scholar] [CrossRef]
  51. Whitty, A. Growing PAINS in academic drug discovery. Futur. Med. Chem. 2011, 3, 797–801. [Google Scholar] [CrossRef] [Green Version]
  52. Sussman, J.L.; Lin, D.; Jiang, J.; Manning, N.O.; Prilusky, J.; Ritter, O.; Abola, E.E. Protein Data Bank (PDB): Database of Three-Dimensional Structural Information of Biological Macromolecules. Acta Crystallogr. Sect. D Biol. Crystallogr. 1998, 54, 1078–1084. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Jo, S.; Cheng, X.; Islam, M.S.; Huang, L.; Rui, H.; Zhu, A.; Lee, H.S.; Qi, Y.; Han, W.; Vanommeslaeghe, K.; et al. CHARMM-GUI PDB Manipulator for Advanced Modeling and Simulations of Proteins Containing Nonstandard Residues. Adv. Protein Chem. Struct. Biol. 2014, 96, 235–265. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Anandakrishnan, R.; Aguilar, B.; Onufriev, A.V. H++ 3.0: Automating pK prediction and the preparation of biomolecular structures for atomistic molecular modeling and simulations. Nucleic Acids Res. 2012, 40, W537–W541. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Kaliappan, S.; Bombay, I.I.T. UCSF Chimera-Overview. 2018. Available online: http://doer.col.org/handle/123456789/9120 (accessed on 15 November 2022).
  56. Dallakyan, S.; Olson, A.J. Small-Molecule Library Screening by Docking with PyRx. In Chemical Biology; Springer: Berlin/Heidelberg, Germany, 2015; pp. 243–250. [Google Scholar]
  57. Halgren, T. A Merck Molecular Force Field. J. Comput. Chem. 1996, 17, 490–519. [Google Scholar] [CrossRef]
  58. Horoiwa, S.; Yokoi, T.; Masumoto, S.; Minami, S.; Ishizuka, C.; Kishikawa, H.; Ozaki, S.; Kitsuda, S.; Nakagawa, Y.; Miyagawa, H. Structure-based virtual screening for insect ecdysone receptor ligands using MM/PBSA. Bioorg. Med. Chem. 2019, 27, 1065–1075. [Google Scholar] [CrossRef]
  59. Lionta, E.; Spyrou, G.; Vassilatis, D.K.; Cournia, Z. Structure-Based Virtual Screening for Drug Discovery: Principles, Applications and Recent Advances. Curr. Top. Med. Chem. 2014, 14, 1923–1938. [Google Scholar] [CrossRef]
  60. Shaker, B.; Ahmad, S.; Lee, J.; Jung, C.; Na, D. In silico methods and tools for drug discovery. Comput. Biol. Med. 2021, 137, 104851. [Google Scholar] [CrossRef]
  61. Case, D.A.; Belfon, K.; Ben-Shalom, I.; Brozell, S.R.; Cerutti, D.; Cheatham, T.; Cruzeiro, V.W.D.; Darden, T.; Duke, R.E.; Giambasu, G.; et al. Amber 2020, University of California, San Fransisco. J. Amer. Chem. Soc. 2020, 142, 3823–3835. [Google Scholar]
  62. Case, D.A.; Babin, V.; Berryman, J.T.; Betz, R.M.; Cai, Q.; Cerutti, D.S.; Cheatham, T.E., III; Darden, T.A.; Duke, R.E.; Gohlke, H.; et al. The FF14SB Force Field. Amber 2014, 14, 29–31. [Google Scholar]
  63. Dickson, C.J.; Rosso, L.; Betz, R.M.; Walker, R.C.; Gould, I.R. GAFFlipid: A General Amber Force Field for the accurate molecular dynamics simulation of phospholipid. Soft Matter 2012, 8, 9617–9627. [Google Scholar] [CrossRef]
  64. He, X.; Liu, S.; Lee, T.-S.; Ji, B.; Man, V.H.; York, D.M.; Wang, J. Fast, Accurate, and Reliable Protocols for Routine Calculations of Protein–Ligand Binding Affinities in Drug Design Projects Using AMBER GPU-TI with ff14SB/GAFF. ACS Omega 2020, 5, 4611–4619. [Google Scholar] [CrossRef] [PubMed]
  65. Maier, J.A.; Martinez, C.; Kasavajhala, K.; Wickstrom, L.; Hauser, K.E.; Simmerling, C. ff14SB: Improving the accuracy of protein side chain and backbone parameters from ff99SB. J. Chem. Theory Comput. 2015, 11, 3696–3713. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Wang, J.; Wang, W.; Kollman, P.A.; Case, D.A. Antechamber: An Accessory Software Package for Molecular Mechanical Calculations. J. Am. Chem. Soc. 2001, 222, U403. [Google Scholar]
  67. Kräutler, V.; Van Gunsteren, W.F.; Hünenberger, P.H. A Fast SHAKE Algorithm to Solve Distance Constraint Equations for Small Molecules in Molecular Dynamics Simulations. J. Comput. Chem. 2001, 22, 501–508. [Google Scholar] [CrossRef]
  68. Izaguirre, J.A.; Catarello, D.P.; Wozniak, J.M.; Skeel, R.D. Langevin stabilization of molecular dynamics. J. Chem. Phys. 2001, 114, 2090–2098. [Google Scholar] [CrossRef] [Green Version]
  69. Roe, D.R.; Cheatham III, T.E. PTRAJ and CPPTRAJ: Software for Processing and Analysis of Molecular Dynamics Trajectory Data. J. Chem. Theory Comput. 2013, 9, 3084–3095. [Google Scholar] [CrossRef]
  70. Turner, P.J. XMGRACE, Version 5.1. 19; Center Coastal Land-Margin Research, Oregon Graduate Institute of Science and Technology: Beaverton, OR, USA, 2005. [Google Scholar]
  71. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  72. Genheden, S.; Kuhn, O.; Mikulskis, P.; Hoffmann, D.; Ryde, U. The Normal-Mode Entropy in the MM/GBSA Method: Effect of System Truncation, Buffer Region, and Dielectric Constant. J. Chem. Inf. Model. 2012, 52, 2079–2088. [Google Scholar] [CrossRef] [Green Version]
  73. Miller, B.R.; McGee, T.D.; Swails, J.M.; Homeyer, N.; Gohlke, H.; Roitberg, A.E. MMPBSA.py: An Efficient Program for End-State Free Energy Calculations. J. Chem. Theory Comput. 2012, 8, 3314–3321. [Google Scholar] [CrossRef]
  74. Sahakyan, H. Improving virtual screening results with MM/GBSA and MM/PBSA rescoring. J. Comput. Mol. Des. 2021, 35, 731–736. [Google Scholar] [CrossRef]
  75. Duan, L.; Liu, X.; Zhang, J.Z. Interaction Entropy: A New Paradigm for Highly Efficient and Reliable Computation of Protein–Ligand Binding Free Energy. J. Am. Chem. Soc. 2016, 138, 5722–5728. [Google Scholar] [CrossRef] [PubMed]
  76. Woods, C.J.; Malaisree, M.; Michel, J.; Long, B.; McIntosh-Smith, S.; Mulholland, A.J. Rapid decomposition and visualisation of protein–ligand binding free energies by residue and by water. Faraday Discuss. 2014, 169, 477–499. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Woods, C.J.; Malaisree, M.; Hannongbua, S.; Mulholland, A.J. A water-swap reaction coordinate for the calculation of absolute protein–ligand binding free energies. J. Chem. Phys. 2011, 134, 054114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [Green Version]
Molecules 28 03186 g001 550
Figure 1. Docking conformation of compounds (shown by color ball and sticks) with TgAPN2 (shown by solid cartoon ribbon). (A) LAS_52160953, (B) LAS_51177972, and (C) LAS_52506311.
Figure 1. Docking conformation of compounds (shown by color ball and sticks) with TgAPN2 (shown by solid cartoon ribbon). (A) LAS_52160953, (B) LAS_51177972, and (C) LAS_52506311.
Molecules 28 03186 g001
Molecules 28 03186 g002 550
Figure 2. Molecular dynamic simulation analysis done using the CPPTRAJ module. (A) RMSD, (B) RMSF, and (C) RoG. All the analyses are plotted versus ångström.
Figure 2. Molecular dynamic simulation analysis done using the CPPTRAJ module. (A) RMSD, (B) RMSF, and (C) RoG. All the analyses are plotted versus ångström.
Molecules 28 03186 g002
Molecules 28 03186 g003 550
Figure 3. WaterSwap analysis energy values of complexes.
Figure 3. WaterSwap analysis energy values of complexes.
Molecules 28 03186 g003
Table
Table 1. Compounds identified via virtual screening that showed stable binding affinity for TgAPN2 binding pocket.
Table 1. Compounds identified via virtual screening that showed stable binding affinity for TgAPN2 binding pocket.
LigandStructureIUPAC NameBinding Affinity
LAS_52160953Molecules 28 03186 i0011-(3-(3-(2-methyl-2,3-dihydro-1H-imidazol-1-yl)propoxy)benzyl)-4-((3,4,5-trimethylphenoxy)methyl)piperidin-4-ol−8.6
LAS_51177972Molecules 28 03186 i0025-carboxy-2-(2-(2,5-dimethoxyphenyl)-1H-pyrrol-1-yl)-4-((3,4,5-trimethylphenoxy)methyl)pyrimidin-1-ium−8.5
LAS_52506311Molecules 28 03186 i0033-(phenoxymethyl)-1-(4-(2-(piperidin-1-yl)ethoxy)benzyl)piperidin-3-ol−8.3
LAS_52160943Molecules 28 03186 i0044-methoxy-1-(3-(3-(2-methyl-2,3-dihydro-1H-imidazol-1-yl)propoxy)benzyl)-4-((m-tolyloxy)methyl)piperidine−7.8
LAS_52157607Molecules 28 03186 i0053-methyl-1-(4-(2-(2-methyl-2,3-dihydro-1H-imidazol-1-yl)ethoxy)benzyl)-4-(m-tolyloxy)piperidin-3-ol−7.7
LAS_52160863Molecules 28 03186 i0061-(3-(3-(2-methyl-2,3-dihydro-1H-imidazol-1-yl)propoxy)benzyl)-3-(phenoxymethyl)piperidin-3-ol−7.7
LAS_52171211Molecules 28 03186 i0073-(phenoxymethyl)-1-(4-(2-(piperidin-1-yl)ethoxy)benzyl)piperidin-3-ol−7.6
LAS_52506188Molecules 28 03186 i0081-(3-(3-(2-methyl-2,3-dihydro-1H-imidazol-1-yl)propoxy)benzyl)-3-((2,3,5-trimethylphenoxy)methyl)piperidine-3,4-diol−7.4
LAS_52157615Molecules 28 03186 i0091-(4-(2-aminoethoxy)benzyl)-3-methyl-4-(m-tolyloxy)piperidin-3-ol−7.1
LAS_32135590Molecules 28 03186 i0104-((4-hydroxytetrahydro-2H-pyran-4-yl)methyl)-N,N,2-trimethylmorpholine-2-carboxamide−6.1
Table
Table 2. MM-GBSA/MM-PBSA binding free energies of docked complexes, given in kcal/mol.
Table 2. MM-GBSA/MM-PBSA binding free energies of docked complexes, given in kcal/mol.
Energy ParameterLAS_52160953 ComplexLAS_51177972 ComplexLAS_52506311 Complex
MM-GBSA
van der Waals Energy−36.89−33.60−36.02
Electrostatic Energy−14.23−12.78−14.85
Delta Gas Phase Energy−51.12−46.38−50.87
Delta Solvation Energy11.6710.1911.47
Net Energy−39.45−36.19−39.4
MM-PBSA
van der Waals Energy−36.89−33.60−36.02
Electrostatic Energy−14.23−12.78−14.85
Delta Gas Phase Energy−51.12−46.38−50.87
Delta Solvation Energy12.5511.1513.75
Net Energy−38.57−35.23−37.12
Table
Table 3. SwissADME analysis of compounds. # stands for number.
Table 3. SwissADME analysis of compounds. # stands for number.
Physicochemical PropertiesLAS_52160953LAS_51177972LAS_52506311
FormulaC29H41N3O3C27H29N3O5+C26H36N2O3
Molecular weight479.65 g/mol475.54 g/mol424.58 g/mol
Num. heavy atoms353531
Num. arom. heavy atoms122312
Fraction Csp30.520.220.54
Num. rotatable bonds1089
Num. H-bond acceptors455
Num. H-bond donors231
Molar Refractivity153.33134.40132.19
TPSA57.20 Å298.20 Å245.17 Å2
Lipophilicity
Log Po/w (iLOGP)4.713.604.50
Log Po/w (XLOGP3)5.224.663.97
Log Po/w (WLOGP)3.213.843.04
Log Po/w (MLOGP)2.832.472.72
Log Po/w (SILICOS-IT)4.795.034.29
Consensus Log Po/w4.153.923.70
Water Solubility
Log S (ESOL)−5.70−5.68−4.67
Solubility9.66 × 10−4 mg/mL; 2.01 × 10−6 mol/L9.88 × 10−4 mg/mL; 2.08 × 10−6 mol/L9.16 × 10−3 mg/mL; 2.16 × 10−5 mol/L
ClassModerately solubleModerately solubleModerately soluble
Log S (Ali)−6.17−6.45−4.62
Solubility3.25 × 10−4 mg/mL; 6.77 × 10−7 mol/L1.69 × 10−4 mg/mL; 3.56 × 10−7 mol/L1.02 × 10−2 mg/mL; 2.40 × 10−5 mol/L
ClassPoorly solublePoorly solubleModerately soluble
Log S (SILICOS-IT)−7.47−8.38−6.78
Solubility1.63 × 10−5 mg/mL; 3.39 × 10−8 mol/L2.00 × 10−6 mg/mL; 4.21 × 10−9 mol/L7.11 × 10−5 mg/mL; 1.67 × 10−7 mol/L
ClassPoorly solublePoorly solublePoorly soluble
Pharmacokinetics
GI absorptionHighHighHigh
BBB permeantYesNoYes
P-gp substrateYesNoYes
CYP1A2 inhibitorNoNoNo
CYP2C19 inhibitorNoYesNo
CYP2C9 inhibitorNoYesNo
CYP2D6 inhibitorYesYesYes
CYP3A4 inhibitorYesYesYes
Log Kp (skin permeation)−5.52 cm/s−5.89 cm/s−6.07 cm/s
Druglikeness
LipinskiYes; 0 violationYes; 0 violationYes; 0 violation
GhoseNo; 2 violations: MR > 130, #atoms > 70No; 1 violation: MR > 130No; 1 violation: MR > 130
VeberYesYesYes
EganYesYesYes
MueggeNo; 1 violation: XLOGP3 > 5YesYes
Bioavailability Score0.550.550.55
Medicinal Chemistry
PAINS0 alert0 alert0 alert
Brenk0 alertNo; 3 violations: MW > 350, Rotors > 7, XLOGP3 > 3.5No; 3 violations: MW > 350, Rotors > 7, XLOGP3 > 3.5
Lead-likenessNo; 3 violations: MW > 350, Rotors > 7, XLOGP3 > 3.5No; 1 violations: MW > 350No; 1 violations: MW > 350
Synthetic accessibility4.874.764.81
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Altharawi, A. Targeting Toxoplasma gondii ME49 TgAPN2: A Bioinformatics Approach for Antiparasitic Drug Discovery. Molecules 2023, 28, 3186. https://doi.org/10.3390/molecules28073186

AMA Style

Altharawi A. Targeting Toxoplasma gondii ME49 TgAPN2: A Bioinformatics Approach for Antiparasitic Drug Discovery. Molecules. 2023; 28(7):3186. https://doi.org/10.3390/molecules28073186

Chicago/Turabian Style

Altharawi, Ali. 2023. "Targeting Toxoplasma gondii ME49 TgAPN2: A Bioinformatics Approach for Antiparasitic Drug Discovery" Molecules 28, no. 7: 3186. https://doi.org/10.3390/molecules28073186

APA Style

Altharawi, A. (2023). Targeting Toxoplasma gondii ME49 TgAPN2: A Bioinformatics Approach for Antiparasitic Drug Discovery. Molecules, 28(7), 3186. https://doi.org/10.3390/molecules28073186

Article Metrics

Back to TopTop