A profound computational study to prioritize the natural compound inhibitors against the P. falciparum orotidine-5-monophosphate decarboxylase enzyme

Abstract

In the current contribution, a multicomplex-based pharmacophore modeling approach was employed on the structural proteome of Plasmodium falciparum orotidine-5-monophosphate decarboxylase enzyme (PfOMPDC). Among the constructed pharmacophore models, the representative hypotheses were selected as the primary filter to screen the molecules with the complementary features responsible for showing inhibition. Thereafter, auxiliary evaluations were performed on the screened candidates via drug-likeness and molecular docking studies. Subsequently, the stability of the docked protein-ligand complexes was scrutinized by employing molecular dynamics simulations and molecular mechanics-Poisson Boltzmann surface area based free binding energy calculations. The stability the docked candidates was compared with the highly active crystallized inhibitor (3S9Y-FNU) to seek more potential candidates. All the docked molecules displayed stable dynamic behavior and high binding free energy in comparison to 3S9Y-FNU. The employed workflow resulted in the retrieval of five drug-like candidates with diverse scaffolds that may show inhibitory activity against PfOMPDC and could be further used as the novel scaffold to develop novel antimalarials.

Communicated by Ramaswamy H. Sarma     

Insight into the structural requirements of thiophene-3-carbonitriles-based MurF inhibitors by 3D-QSAR,molecular docking and molecular dynamics study

The discovery of clinically relevant inhibitors against MurF enzyme has proven to be a challenging task. In order to get further insight into the structural features required for the MurF inhibitory activity, we performed pharmacophore and atom-based three-dimensional quantitative structure–activity relationship studies for novel thiophene-3-carbonitriles based MurF inhibitors. The five-feature pharmacophore model was generated using 48 inhibitors having IC₅₀ values ranging from 0.18 to 663?μm. The best-fitted model showed a higher coefficient of determination (R²?=?0.978), cross-validation coefficient (Q²?=?0.8835) and Pearson coefficient (0.9406) at four component partial least-squares factor. The model was validated with external data set and enrichment study. The effectiveness of the docking protocol was validated by docking the co-crystallized ligand into the catalytic pocket of MurF enzyme. Further, binding free energy calculated by the molecular mechanics generalized Born surface area approach showed that van der Waals and non-polar solvation energy terms are the main contributors to ligand binding in the active site of MurF enzyme. A 10-ns molecular dynamic simulation was performed to confirm the stability of the 3ZM6-ligand complex. Four new molecules are also designed as potent MurF inhibitors. These results provide insights regarding the development of novel MurF inhibitors with better binding affinity.     

Discover potential inhibitors for PFKFB3 using 3D-QSAR,virtual screening,molecular docking and molecular dynamics simulation

Abstract

The 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase-3 (PFKFB3) is a master regulator of glycolysis in cancer cells by synthesizing fructose-2,6-bisphosphate (F-2,6-BP), a potent allosteric activator of phosphofructokinase-1 (PFK-1), which is a rate-limiting enzyme of glycolysis. PFKFB3 is an attractive target for cancer treatment. It is valuable to discover promising inhibitors by using 3D-QSAR pharmacophore modeling, virtual screening, molecular docking and molecular dynamics simulation. Twenty molecules with known activity were used to build 3D-QSAR pharmacophore models. The best pharmacophore model was ADHR called Hypo1, which had the highest correlation value of 0.98 and the lowest RMSD of 0.82. Then, the Hypo1 was validated by cost value method, test set method and decoy set validation method. Next, the Hypo1 combined with Lipinski's rule of five and ADMET properties were employed to screen databases including Asinex and Specs, total of 1,048,159 molecules. The hits retrieved from screening were docked into protein by different procedures including HTVS, SP and XP. Finally, nine molecules were picked out as potential PFKFB3 inhibitors. The stability of PFKFB3-lead complexes was verified by 40?ns molecular dynamics simulation. The binding free energy and the energy contribution of per residue to the binding energy were calculated by MM-PBSA based on molecular dynamics simulation.     

Computational investigations of gram-negative bacteria phosphopantetheine adenylyltransferase inhibitors using 3D-QSAR,molecular docking and molecular dynamic simulations

Abstract

Phosphopantetheine adenylyltransferase (PPAT) has been recognized as a promising target to develop novel antimicrobial agents, which is a hexameric enzyme that catalyzes the penultimate step in coenzyme A biosynthesis. In this work, molecular modeling study was performed with a series of PPAT inhibitors using molecular docking, three-dimensional qualitative structure-activity relationship (3D-QSAR) and molecular dynamic (MD) simulations to reveal the structural determinants for their bioactivities. Molecular docking study was applied to understand the binding mode of PPAT with its inhibitors. Subsequently, 3D-QSAR model was constructed to find the features required for different substituents on the scaffolds. For the best comparative molecular field analysis (CoMFA) model, the Q² and R² values of which were calculated as 0.702 and 0.989, while they were calculated as 0.767 and 0.983 for the best comparative molecular similarity index analysis model. The statistical data verified the significance and accuracy of our 3D-QSAR models. Furthermore, MD simulations were carried out to evaluate the stability of the receptor–ligand contacts in physiological conditions, and the results were consistent with molecular docking studies and 3D-QSAR contour map analysis. Binding free energy was calculated with molecular mechanics generalized born surface area approach, the result of which coincided well with bioactivities and demonstrated that van der Waals accounted for the largest portion. Overall, our study provided a valuable insight for further research work on the recognition of potent PPAT inhibitors.

Communicated by Ramaswamy H. Sarma     

Pharmacophore generation,atom-based 3D-QSAR,molecular docking and molecular dynamics simulation studies on benzamide analogues as FtsZ inhibitors

FtsZ is an appealing target for the design of antimicrobial agent that can be used to defeat the multidrug-resistant bacterial pathogens. Pharmacophore modelling, molecular docking and molecular dynamics (MD) simulation studies were performed on a series of three-substituted benzamide derivatives. In the present study a five-featured pharmacophore model with one hydrogen bond acceptors, one hydrogen bond donors, one hydrophobic and two aromatic rings was developed using 97 molecules having MIC values ranging from .07 to 957 μM. A statistically significant 3D-QSAR model was obtained using this pharmacophore hypothesis with a good correlation coefficient (R² = .8319), cross validated coefficient (Q² = .6213) and a high Fisher ratio (F = 103.9) with three component PLS factor. A good correlation between experimental and predicted activity of the training (R² = .83) and test set (R² = .67) molecules were displayed by ADHRR.1682 model. The generated model was further validated by enrichment studies using the decoy test and MAE-based criteria to measure the efficiency of the model. The docking studies of all selected inhibitors in the active site of FtsZ protein showed crucial hydrogen bond interactions with Val 207, Asn 263, Leu 209, Gly 205 and Asn-299 residues. The binding free energies of these inhibitors were calculated by the molecular mechanics/generalized born surface area VSGB 2.0 method. Finally, a 15 ns MD simulation was done to confirm the stability of the 4DXD–ligand complex. On a wider scope, the prospect of present work provides insight in designing molecules with better selective FtsZ inhibitory potential.     

Virtual screening of small molecules databases for discovery of novel PARP-1 inhibitors: combination of in silico and in vitro studies

Poly(ADP-ribose) polymerase-1 (PARP-1) enzyme has critical roles in DNA replication repair and recombination. Thus, PARP-1 inhibitors play an important role in the cancer therapy. In the current study, we have performed combination of in silico and in vitro studies in order to discover novel inhibitors against PARP-1 target. Structure-based virtual screening was carried out for an available small molecules database. A total of 257,951 ligands from Otava database were screened at the binding pocket of PARP-1 using high-throughput virtual screening techniques. Filtered structures based on predicted binding energy results were then used in more sophisticated molecular docking simulations (i.e. Glide/standard precision, Glide/XP, induced fit docking – IFD, and quantum mechanics polarized ligand docking – QPLD). Potential high binding affinity compounds that are predicted by molecular simulations were then tested by in vitro methods. Computationally proposed compounds as PARP-1 inhibitors (Otava Compound Codes: 7111620047 and 7119980926) were confirmed by in vitro studies. In vitro results showed that compounds 7111620047 and 7119980926 have IC₅₀ values of 0.56 and 63 μM against PARP-1 target, respectively. The molecular mechanism analysis, free energy perturbation calculations using long multiple molecular dynamics simulations for the discovered compounds which showed high binding affinity against PARP-1 enzyme, as well as structure-based pharmacophore development (E-pharmacophore) studies were also studied.     

Structure-based screening and molecular dynamics simulations offer novel natural compounds as potential inhibitors of Mycobacterium tuberculosis isocitrate lyase

Mycobacterium tuberculosis is the etiological agent of tuberculosis in humans and is responsible for more than two million deaths annually. M. tuberculosis isocitrate lyase (MtbICL) catalyzes the first step in the glyoxylate cycle, plays a pivotal role in the persistence of M. tuberculosis, which acts as a potential target for an anti-tubercular drug. To identify the potential anti-tuberculosis compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,748) against the MtbICL structure. The ligands were docked against MtbICL in three sequential docking modes that resulted in 340 ligands having better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 27 compounds were found to fit well with re-docking studies. After refinement by molecular docking and drug-likeness analyses, three potential inhibitors (ZINC1306071, ZINC2111081, and ZINC2134917) were identified. These three ligands and the reference compounds were further subjected to molecular dynamics simulation and binding energy analyses to compare the dynamic structure of protein after ligand binding and the stability of the MtbICL and bound complexes. The binding free energy analyses were calculated to validate and capture the intermolecular interactions. The results suggested that the three compounds had a negative binding energy with ?96.462, ?143.549, and ?122.526 kJ mol^?1 for compounds with IDs ZINC1306071, ZINC2111081, and ZINC2134917, respectively. These lead compounds displayed substantial pharmacological and structural properties to be drug candidates. We concluded that ZINC2111081 has a great potential to inhibit MtbICL and would add to the drug discovery process against tuberculosis.     

Probing the binding mechanism of mercaptoguanine derivatives as inhibitors of HPPK by docking and molecular dynamics simulations

6-Hydroxymethyl-7,8-dihydropterin pyrophosphokinase (HPPK) is a promising antimicrobial target involved in the folate biosynthesis pathway. Although, the results from crystallographic studies of HPPK have attracted a great interest in the design of novel HPPK inhibitors, the mechanism of action of HPPK due to inhibitor binding remains questionable. Recently, mercaptoguanine derivatives were reported to inhibit the pyrophosphoryl transfer mechanism of Staphylococcus aureus HPPK (SaHPPK). The present study is an attempt to understand the SaHPPK-inhibitors binding mechanism and to highlight the key residues that possibly involve in the complex formation. To decipher these questions, we used the state-of-the-art advanced insilico approach such as molecular docking, molecular dynamics (MD), molecular mechanics-generalized Born surface area approach. Domain cross correlation and principle component analysis were applied to the snapshots obtained from MD revealed that the compounds with high binding affinity stabilize the conformational dynamics of SaHPPK. The binding free energy estimation showed that the van der Waals and electrostatic interactions played a vital role for the binding mechanism. Additionally, the predicted binding free energy was in good agreement with the experimental values (R² = .78). Moreover, the free energy decomposition on per-residue confirms the key residues that significantly contribute to the complex formation. These results are expected to be useful for rational design of novel SaHPPK inhibitors.     

Structure-based virtual screening,molecular docking,ADMET and molecular simulations to develop benzoxaborole analogs as potential inhibitor against Leishmania donovani trypanothione reductase

Visceral leishmaniasis (VL) is the most fatal form of leishmaniasis and it affects 70 countries worldwide. Increasing drug resistant for antileishmanial drugs such as miltefosine, sodium stibogluconate and pentamidine has been reported in the VL endemic region. Amphotericin B has shown potential antileishmanial activity in different formulations but its cost of treatment and associated nephrotoxicity have limited its use by affected people living in the endemic zone. To control the VL infection in the affected countries, it is necessary to develop new antileishmanial compounds with high efficacy and negligible toxicity. Computer aided programs such as binding free energy estimation; ADMET prediction and molecular dynamics simulation can be used to investigate novel antileishmanial molecules in shorter duration. To develop antileishmanial lead molecule, we performed standard precision (SP) docking for 1160 benzoxaborole analogs along with reference inhibitors against trypanothione reductase of Leishmania parasite. Furthermore, extra precision (XP) docking, ADMET prediction, prime MM-GBSA was conducted over 115 ligands, showing better docking score than reference inhibitors to get potential antileishmanial compounds. Simultaneously, area under the curve (AUC) was estimated using ROC plot to validate the SP and XP docking protocol. Later on, two benzoxaborole analogs with best MM-GBSA ΔG-bind were subjected to molecular simulation and docking confirmation to ensure the ligand interaction with TR. The presented drug discovery based on computational study confirms that BOB27 can be used as a potential drug candidate and warrants further experimental investigation to fight against VL in endemic areas.     

Identification of anti-filarial leads against aspartate semialdehyde dehydrogenase of Wolbachia endosymbiont of Brugia malayi: combined molecular docking and molecular dynamics approaches

Lymphatic filariasis is a debilitating vector borne parasitic disease that infects human lymphatic system by nematode Brugia malayi. Currently available anti-filarial drugs are effective only on the larval stages of parasite. So far, no effective drugs are available for humans to treat filarial infections. In this regard, aspartate semialdehyde dehydrogenase (ASDase) in lysine biosynthetic pathway from Wolbachia endosymbiont Brugia malayi represents an attractive therapeutic target for the development of novel anti-filarial agents. In this present study, molecular modeling combined with molecular dynamics simulations and structure-based virtual screening were performed to identify potent lead molecules against ASDase. Based on Glide score, toxicity profile, binding affinity and mode of interactions with the ASDase, five potent lead molecules were selected. The molecular docking and dynamics results revealed that the amino acid residues Arg103, Asn133, Cys134, Gln161, Ser164, Lys218, Arg239, His246, and Asn321 plays a crucial role in effective binding of Top leads into the active site of ASDase. The stability of the ASDase-lead complexes was confirmed by running the 30 ns molecular dynamics simulations. The pharmacokinetic properties of the identified lead molecules are in the acceptable range. Furthermore, density functional theory and binding free energy calculations were performed to rank the lead molecules. Thus, the identified lead molecules can be used for the development of anti-filarial agents to combat the pathogenecity of Brugia malayi.     

2D-QSAR study,molecular docking,and molecular dynamics simulation studies of interaction mechanism between inhibitors and transforming growth factor-beta receptor I (ALK5)

Transforming growth factor type 1 receptor (ALK5) is kinase associated with a wide variety of pathological processes, and inhibition of ALK5 is a good strategy to treat many kinds of cancer and fibrotic diseases. Recently, a series of compounds have been synthesized as ALK5 inhibitors. However, the study of their selectivity against other potential targets remains elusive. In this research, a data-set of ALK5 inhibitors were collected and studied based on the combination of 2D-QSAR, molecular docking and molecular dynamics simulation. The quality of QSAR models were assessed statistically by F, R², and R²_ADJ, proved to be credible. The cross-validations for the models (q²_LOO = 0.571 and 0.629, respectively) showed their robustness, while the external validations (r²_test = 0.703 and 0.764, respectively) showed their predictive power. Besides, the predicted binding free energy results calculated by MM/GBSA method were in accordance with the experimental data, and the van der Waals energy term was the factor that had the most significant impact on ligand binding. What is more, several important residues were found to significantly affect the binding affinity. Finally, based on our analyses above, a proposed series of molecules were designed.     

The influence of glycerol on the binding of methotrexate to Mycobacterium tuberculosis dihydrofolate reductase: a molecular modelling study

The enzyme, dihydrofolate reductase (DHFR), from Mycobacterium tuberculosis (mt-DHFR) is believed to be a potential drug target for the treatment of tuberculosis. The co-crystal structure of mt-DHFR bound with glycerol (GOL), NAPDH and methotrexate (MTX) reveals a GOL binding site on the enzyme. This GOL binding site could be very important for the design of novel, selective mt-DHFR inhibitors, because this binding site is absent on human DHFR (h-DHFR). We have performed molecular dynamic simulations and free energy calculations to evaluate the binding affinity of GOL and its free energy contribution to the binding of MTX to mt-DHFR. The results showed that GOL does not bind tightly to mt-DHFR. Although GOL itself contributed free energy on MTX binding to mt-DHFR, GOL also increased the flexibilities of MTX, so that MTX cannot maintain strong electronic interactions with ARG32 and ARG60, which caused the total binding free energy to decrease. These data suggest that GOL binding is weak and it could be expelled from the binding site, to allow inhibitors containing appropriate side chains to bind. This observation can be used to inform future drug design studies, especially those aimed at improving drug selectivity against h-DHFR.     

Molecular modeling study for the design of novel acetyl-CoA carboxylase inhibitors using 3D QSAR,molecular docking and dynamic simulations

Acetyl-CoA carboxylase (ACC) enzyme plays an important role in the regulation of biosynthesis and oxidation of fatty acids. ACC is a recognized drug target for the treatment of obesity and diabetes. Combination of ligand and structure-based in silico methods along with activity and toxicity prediction provides best lead compounds in the drug discovery process. In this study, a data-set of 100 ACC inhibitors were used for the development of comparative molecular field analysis (CoMFA) and comparative molecular similarity index matrix analysis (CoMSIA) models. The generated contour maps were used for the design of novel ACC inhibitors. CoMFA and CoMSIA models were used for the predication of activity of designed compounds. In silico toxicity risk prediction study was carried out for the designed compounds. Molecular docking and dynamic simulations studies were performed to know the binding mode of designed compounds with the ACC enzyme. The designed compounds showed interactions with key amino acid residues important for catalysis, and good correlation was observed between binding free energy and inhibition of ACC.     

Exploration of potential RSK2 inhibitors by pharmacophore modelling,structure-based 3D-QSAR,molecular docking study and molecular dynamics simulation

Abstract

The p90 ribosomal s6 kinase 2 (RSK2) is a promising target because of its over expression and activation in human cancer cells and tissues. Over the last few years, significant efforts have been made in order to develop RSK2 inhibitors to treat myeloma, prostatic cancer, skin cancer and etc., but with limited success so far. In this paper, pharmacophore modelling, molecular docking study and molecular dynamics (MD) simulation have been performed to explore the novel inhibitors of RSK2. Pharmacophore models were developed by 95 molecules having pIC₅₀ ranging from 4.577 to 9.000. The pharmacophore model includes one hydrogen bond acceptor (A), one hydrogen bond donor (D), one hydrophobic feature (H) and one aromatic ring (R). It is the best pharmacophore hypothesis that has the highest correlation coefficient (R² = 0.91) and cross validation coefficient (Q² = 0.71) at 5 component PLS factor. It was evaluated using enrichment analysis and the best model was used for virtual screening. The constraints used in this study were docking score, ADME properties, binding free energy estimates and IFD Score to screen the database. Ultimately, 12 hits were identified as potent and novel RSK2 inhibitors. A 15 ns molecular dynamics (MD) simulation was further employed to validate the reliability of the docking results.     

Structural insight into Mycobacterium tuberculosis maltosyl transferase inhibitors: pharmacophore-based virtual screening,docking, and molecular dynamics simulations

Pharmacophore-based virtual screening, subsequent docking, and molecular dynamics (MD) simulations have been done to identify potential inhibitors of maltosyl transferase of Mycobacterium tuberculosis (mtb GlgE). Ligand and structure-based pharmacophore models representing its primary binding site (pbs) and unique secondary binding site 2 (sbs2), respectively, were constructed based on the three dimensional structure of mtb GlgE. These pharmacophore models were further used for screening of ZINC and antituberculosis compounds database (ATD). Virtually screened molecules satisfying Lipinski’s rule of five were then analyzed using docking studies and have identified 23 molecules with better binding affinity than its natural substrate, maltose. Four top scoring ligands from ZINC and ATD that either binds to pbs or sbs2 have been subjected to 10 ns each MD simulations and binding free energy calculations. Results of these studies have confirmed stable protein ligand binding. Results reported in the article are likely to be helpful in antitubercular therapeutic development research.     

Deciphering the crucial molecular properties of a series of Benzothiazole Hydrazone inhibitors that targets anti-apoptotic Bcl-xL protein

The Bcl-2 family proteins are the central regulators of apoptosis. Due to its predominant role in cancer progression, the Bcl-2 family proteins act as attractive therapeutic targets. Recently, molecular series of Benzothiazole Hydrazone (BH) inhibitors that exhibits drug-likeness characteristics, which selectively targets Bcl-xL have been reported. In the present study, docking was used to explore the plausible binding mode of the highly active BH inhibitor with Bcl-xL; and Molecular Dynamics (MD) simulation was applied to investigate the stability of predicted conformation over time. Furthermore, the molecular properties of the series of BH inhibitors were extensively investigated by pharmacophore based 3D-QSAR model. The docking correctly predicted the binding mode of the inhibitor inside the Bcl-xL hydrophobic groove, whereas the MD-based free energy calculation exhibited the binding strength of the complex over the time period. Furthermore, the residue decomposition analysis revealed the major energy contributing residues – F105, L108, L130, N136, and R139 – involved in complex stability. Additionally, a six-featured pharmacophore model – AAADHR.89 – was developed using the series of BH inhibitors that exhibited high survival score. The statistically significant 3D-QSAR model exhibited high correlation co-efficient (R² = .9666) and cross validation co-efficient (Q² = .9015) values obtained from PLS regression analysis. The results obtained from the current investigation might provide valuable insights for rational drug design of Bcl-xL inhibitor synthesis.     

Thermodynamics of site-specific small molecular ion interactions with DNA duplex: a molecular dynamics study

The stability and dynamics of a double-stranded DNA (dsDNA) is affected by the preferential occupancy of small monovalent molecular ions. Small metal and molecular ions such as sodium and alkyl ammonium have crucial biological functions in human body, affect the thermodynamic stability of the duplex DNA and exhibit preferential binding. Here, using atomistic molecular dynamics simulations, we investigate the preferential binding of metal ion such as Na⁺ and molecular ions such as tetramethyl ammonium (TMA⁺) and 2-hydroxy-N,N,N-trimethylethanaminium (CHO⁺) to double-stranded DNA. The thermodynamic driving force for a particular molecular ion-DNA interaction is determined by decomposing the free energy of binding into its entropic and enthalpic contributions. Our simulations show that each of these molecular ions preferentially binds to the minor groove of the DNA and the extent of binding is highest for CHO⁺. The ion binding processes are found to be entropically favourable. In addition, the contribution of hydrophobic effects towards the entropic stabilisation (in case of TMA⁺) and the effect of hydrogen bonding contributing to enthalpic stabilisation (in case of CHO⁺) have also been investigated.     

Insight into the binding modes and inhibition mechanisms of adamantyl-based 1,3-disubstituted urea inhibitors in the active site of the human soluble epoxide hydrolase

Soluble epoxide hydrolase (sEH) is a promising new target for treating hypertension and inflammation. Considerable efforts have been devoted to develop novel inhibitors. In this study, the binding modes and interaction mechanisms of a series of adamantyl-based 1,3-disubstituted urea inhibitors were investigated by molecular docking, molecular dynamics simulations, binding free energy calculations, and binding energy decomposition analysis. Based on binding affinity, the most favorable binding mode was determined for each inhibitor. The calculation results indicate that the total binding free energy (ΔG_TOT, the sum of enthalpy ΔG_MM-GB/SA, and entropy ?TΔS) presents a good correlation with the experimental inhibitory activity (IC₅₀, r²?=?.99). The van der Waals energy contributes most to the total binding free energy (ΔG_TOT). A detailed discussion on the interactions between inhibitors and those residues located in the active pocket is made based on hydrogen bond and binding modes analysis. According to binding energy decomposition, the residues Asp333 and Trp334 contribute the most to binding free energy in all systems. Furthermore, Hip523 plays a major role in determining this class of inhibitor-binding orientations. Combined with the results of hydrogen bond analysis and binding free energy, we believe that the conserved hydrogen bonds play a role only in anchoring the inhibitors to the exact site for binding and the number of hydrogen bonds may not directly relate to the binding free energy. The results we obtained will provide valuable information for the design of high potency sEH inhibitors.     

Identification of potent inhibitors against snake venom metalloproteinase (SVMP) using molecular docking and molecular dynamics studies

Snake venom metalloproteinase (SVMP) (Echis coloratus (Carpet viper) is a multifunctional enzyme that is involved in producing several symptoms that follow a snakebite, such as severe local hemorrhage, nervous system effects and tissue necrosis. Because the three-dimensional (3D) structure of SVMP is not known, models were constructed, and the best model was selected based on its stereo-chemical quality. The stability of the modeled protein was analyzed through molecular dynamics (MD) simulation studies. Structure-based virtual screening was performed, and 15 potential molecules with the highest binding energies were selected. Further analysis was carried out with induced fit docking, Prime/MM–GBSA (ΔG_Bind calculations), quantum-polarized ligand docking, and density functional theory calculations. Further, the stability of the lead molecules in the SVMP-active site was examined using MD simulation. The results showed that the selected lead molecules were highly stable in the active site of SVMP. Hence, these molecules could potentially be selective inhibitors of SVMP. These lead molecules can be experimentally validated, and their backbone structural scaffold could serve as building blocks in designing drug-like molecules for snake antivenom.     

Identification of potential inhibitors of Fasciola gigantica thioredoxin1: computational screening,molecular dynamics simulation,and binding free energy studies

Fasciola gigantica is the causative organism of fascioliasis and is responsible for major economic losses in livestock production globally. F. gigantica thioredoxin1 (FgTrx1) is an important redox-active enzyme involved in maintaining the redox homeostasis in the cell. To identify a potential anti-fasciolid compound, we conducted a structure-based virtual screening of natural compounds from the ZINC database (n = 1,67,740) against the FgTrx1 structure. The ligands were docked against FgTrx1 and 309 ligands were found to have better docking score. These compounds were evaluated for Lipinski and ADMET prediction, and 30 compounds were found to fit well for re-docking studies. After refinement by molecular docking and drug-likeness analysis, three potential inhibitors (ZINC15970091, ZINC9312362, and ZINC9312661) were identified. These three ligands were further subjected to molecular dynamics simulation (MDS) to compare the dynamics and stability of the protein structure after binding of the ligands. The binding free energy analyses were calculated to determine the intermolecular interactions. The results suggested that the two compounds had a binding free energy of –82.237, and –109.52 kJ.mol^?1 for compounds with IDs ZINC9312362 and ZINC9312661, respectively. These predicted compounds displayed considerable pharmacological and structural properties to be drug candidates. We concluded that these two compounds could be potential drug candidates to fight against F. gigantica parasites.