MM-PBSA Studies Research Articles

Menthol carbonates as potent antiparasitic agents: synthesis and in vitro studies along with computer-aided approaches

IntroductionDespite the number of deaths and the significant economic and social costs associated with Chagas, Leishmaniasis and Malaria diseases worldwide, available drugs are limited and have serious side effects and high toxicity for the patient. Therefore, there is an urgent need for safe, low-cost, and effective treatments. Natural products are an important source of bioactive compounds and there is current interest in finding natural bioactive molecules that can be used for treating these parasitic diseases. In the present study we proposed to evaluate the in vitro antiparasitic activity of new menthol derivatives against Trypanosoma cruzi, Leishmania braziliensis and Plasmodium falciparum; moreover, we propose to explore their mode of action through in silico approaches.Material and methodsA series of carbonate prodrugs (1–9) were synthesized from menthol with different aliphatic alcohols. Spectroscopic techniques were used to confirm the structures of the synthesized compounds. The cytotoxicity of the compounds was assessed using U-937 cells. In vitro trypanocidal, leishmanicidal and antiplasmodial activity were evaluated using a T. cruzi, L. braziliensis and P. falciparum organism, respectively. In addition, in silico studies were also performed through molecular dynamics simulations and MM-PBSA analysis.ResultsThe assay revealed that most of the compounds were highly active against intracellular amastigotes of T. cruzi and L. braziliensis, and had moderate activity against the total forms of P. falciparum. Compound 2 was one of the drugs that showed a high selectivity index (SI) for the three organisms evaluated. The prediction of the ADME properties suggests that all the compounds have drug-like molecular properties and the probability to be lead candidates. Finally, molecular dynamics simulations, and MM-PBSA studies indicate that menthol at the substrate binding site of TcDHODH, LbDHODH and PfDHODH is structurally stable in the same order as the natural substrate; also, interactions of menthol with residues involved in the inhibition of TcDHODH and PfDHODH proteins were predicted.ConclusionsThe present study demonstrates that menthol prodrugs are promising antiparasitic agents; however, the mechanisms of action proposed in this study need to be experimentally verified by future enzymatic assays.

Are the current commercially available oximes capable of reactivating acetylcholinesterase inhibited by the nerve agents of the A-series?

The misuse of novichok agents in assassination attempts has been reported in the international media since 2018. These relatively new class of neurotoxic agents is claimed to be more toxic than the agents of the G and V series and so far, there is no report yet in literature about potential antidotes against them. To shed some light into this issue, we report here the design and synthesis of NTMGMP, a surrogate of A-242 and also the first surrogate of a novichok agent useful for experimental evaluation of antidotes. Furthermore, the efficiency of the current commercial oximes to reactivate NTMGMP-inhibited acetylcholinesterase (AChE) was evaluated. The Ellman test was used to confirm the complete inhibition of AChE, and to compare the subsequent rates of reactivation in vitro as well as to evaluate aging. In parallel, molecular docking, molecular dynamics and MM-PBSA studies were performed on a computational model of the human AChE (HssAChE)/NTMGMP complex to assess the reactivation performances of the commercial oximes in silico. Experimental and theoretical studies matched the exact hierarchy of efficiency and pointed to trimedoxime as the most promising commercial oxime for reactivation of AChE inhibited by A-242.

Applying polypharmacology approach for drug repurposing for SARS-CoV2.

Exploring the new therapeutic indications of known drugs for treating COVID-19, popularly known as drug repurposing, is emerging as a pragmatic approach especially owing to the mounting pressure to control the pandemic. Targeting multiple targets with a single drug by employing drug repurposing known as the polypharmacology approach may be an optimised strategy for the development of effective therapeutics. In this study, virtual screening has been carried out on seven popular SARS-CoV-2 targets (3CLpro, PLpro, RdRp (NSP12), NSP13, NSP14, NSP15, and NSP16). A total of 4015 approved drugs were screened against these targets. Four drugs namely venetoclax, tirilazad, acetyldigitoxin, and ledipasvir have been selected based on the docking score, ability to interact with four or more targets and having a reasonably good number of interactions with key residues in the targets. The MD simulations and MM-PBSA studies showed reasonable stability of protein-drug complexes and sustainability of key interactions between the drugs with their respective targets throughout the course of MD simulations. The identified four drug molecules were also compared with the known drugs namely elbasvir and nafamostat. While the study has provided a detailed account of the chosen protein-drug complexes, it has explored the nature of seven important targets of SARS-CoV-2 by evaluating the protein-drug complexation process in great detail.Graphical abstractDrug repurposing strategy against SARS-CoV2 drug targets. Computational analysis was performed to identify repurposable approved drug candidates against SARS-CoV2 using approaches such as virtual screening, molecular dynamics simulation and MM-PBSA calculations. Four drugs namely venetoclax, tirilazad, acetyldigitoxin, and ledipasvir have been selected as potential candidates. Supplementary InformationThe online version contains supplementary material available at 10.1007/s12039-022-02046-0.

Ligand and Structure-Based In Silico Determination of the Most Promising SARS-CoV-2 nsp16-nsp10 2'-o-Methyltransferase Complex Inhibitors among 3009 FDA Approved Drugs.

As a continuation of our earlier work against SARS-CoV-2, seven FDA-approved drugs were designated as the best SARS-CoV-2 nsp16-nsp10 2′-o-methyltransferase (2′OMTase) inhibitors through 3009 compounds. The in silico inhibitory potential of the examined compounds against SARS-CoV-2 nsp16-nsp10 2′-o-methyltransferase (PDB ID: (6W4H) was conducted through a multi-step screening approach. At the beginning, molecular fingerprints experiment with SAM (S-Adenosylmethionine), the co-crystallized ligand of the targeted enzyme, unveiled the resemblance of 147 drugs. Then, a structural similarity experiment recommended 26 compounds. Therefore, the 26 compounds were docked against 2′OMTase to reveal the potential inhibitory effect of seven promising compounds (Protirelin, (1187), Calcium folinate (1913), Raltegravir (1995), Regadenoson (2176), Ertapenem (2396), Methylergometrine (2532), and Thiamine pyrophosphate hydrochloride (2612)). Out of the docked ligands, Ertapenem (2396) showed an ideal binding mode like that of the co-crystallized ligand (SAM). It occupied all sub-pockets of the active site and bound the crucial amino acids. Accordingly, some MD simulation experiments (RMSD, RMSF, Rg, SASA, and H-bonding) have been conducted for the 2′OMTase—Ertapenem complex over 100 ns. The performed MD experiments verified the correct binding mode of Ertapenem against 2′OMTase exhibiting low energy and optimal dynamics. Finally, MM-PBSA studies indicated that Ertapenem bonded advantageously to the targeted protein with a free energy value of −43 KJ/mol. Furthermore, the binding free energy analysis revealed the essential amino acids of 2′OMTase that served positively to the binding. The achieved results bring hope to find a treatment for COVID-19 via in vitro and in vivo studies for the pointed compounds.

Structural insights into inhibition of PRRSV Nsp4 revealed by structure-based virtual screening, molecular dynamics, and MM-PBSA studies

BackgroundPorcine reproductive and respiratory syndrome respiratory sickness in weaned and growing pigs, as well as sow reproductive failure, and its infection is regarded as one of the most serious swine illnesses worldwide. Given the current lack of an effective treatment, in this study, we identified natural compounds capable of inhibiting non-structural protein 4 (Nsp4) of the virus, which is involved in their replication and pathogenesis.ResultsWe screened natural compounds (n = 97,999) obtained from the ZINC database against Nsp4 and selected the top 10 compounds for analysing protein–ligand interactions and physicochemical properties. The five compounds demonstrating strong binding affinity were then subjected to molecular dynamics simulations (100 ns) and binding free energy calculations. Based on analysis, we identified four possible lead compounds that represent potentially effective drug-like inhibitors.ConclusionsThese methods identified that these natural compounds are capable of inhibiting Nsp4 and possibly effective as antiviral therapeutics against PRRSV.

Computational investigation of the alkaloids of Pilocarpus microphyllus species as phytopharmaceuticals for the inhibition of sterol 14α-demethylase protease of Trypanosoma cruzi

Trypanosoma cruzi is a protozoan transmitted by the insect Triatoma infestans, popularly known as kissing bug. This protozoan causes the Chagas disease, a Neglected Tropical Disease. This study aimed to investigate, through DFT method and B3LYP hybrid functional, the physicochemical, pharmacokinetic, and pharmacodynamic properties of the alkaloids present in the leaves of the species Pilocarpus microphyllus (jaborandi) as a potential inhibitory activity on the protease sterol 14α-demethylase of T. cruzi associated with the techniques of molecular docking, molecular dynamics, MM-PBSA and ADMET predictions. The molecules of isopilosine, epiisopiloturine, epiisopilosine, and pilosine showed up the lowest binding energies by molecular docking, good human intestinal absorption, low penetration in the blood-brain barrier, antiprotozoal and anticarcinogenic activities in ADMET studies. It has been observed a better binding affinity of the sterol 14α-demethylase protease with isopilosine in molecular dynamics and MM-PBSA studies, which indicates it as a potential drug candidate for Chagas disease. Communicated by Ramaswamy H. Sarma

Repurposing potential of Ayurvedic medicinal plants derived active principles against SARS-CoV-2 associated target proteins revealed by molecular docking, molecular dynamics and MM-PBSA studies

All the plants and their secondary metabolites used in the present study were obtained from Ayurveda, with historical roots in the Indian subcontinent. The selected secondary metabolites have been experimentally validated and reported as potent antiviral agents against genetically-close human viruses. The plants have also been used as a folk medicine to treat cold, cough, asthma, bronchitis, and severe acute respiratory syndrome in India and across the globe since time immemorial. The present study aimed to assess the repurposing possibility of potent antiviral compounds with SARS-CoV-2 target proteins and also with host-specific receptor and activator protease that facilitates the viral entry into the host body. Molecular docking (MDc) was performed to study molecular affinities of antiviral compounds with aforesaid target proteins. The top-scoring conformations identified through docking analysis were further validated by 100 ns molecular dynamic (MD) simulation run. The stability of the conformation was studied in detail by investigating the binding free energy using MM-PBSA method. Finally, the binding affinities of all the compounds were also compared with a reference ligand, remdesivir, against the target protein RdRp. Additionally, pharmacophore features, 3D structure alignment of potent compounds and Bayesian machine learning model were also used to support the MDc and MD simulation. Overall, the study emphasized that curcumin possesses a strong binding ability with host-specific receptors, furin and ACE2. In contrast, gingerol has shown strong interactions with spike protein, and RdRp and quercetin with main protease (Mpro) of SARS-CoV-2. In fact, all these target proteins play an essential role in mediating viral replication, and therefore, compounds targeting aforesaid target proteins are expected to block the viral replication and transcription. Overall, gingerol, curcumin and quercetin own multitarget binding ability that can be used alone or in combination to enhance therapeutic efficacy against COVID-19. The obtained results encourage further in vitro and in vivo investigations and also support the traditional use of antiviral plants preventively.

In silico screening of antiviral compounds from Moringa oleifera for inhibition of SARS-CoV-2 main protease

COVID-19 as per early April 2021, has globally exceeded 129 million cases and is the cause of about 2.8 million deaths. Since the illness is known to have a wide array of symptoms, in addition to its ability to spread rapidly through contact and the complexity involved in developing drugs, there is an immediate need to look into alternative treatment regimes. This study was intended to identify potential antiviral compounds from Moringa oleifera against the selected target main protease (Mpro), which is vital for the survival of the virus. In silico molecular docking and dynamics was performed to determine a potential inhibitor against Mpro. Phytochemicals of Moringa olifera reported in literature were retrieved from public databases and employed for molecular docking studies against Mpro. PyRx software was used to perform docking analysis. Visualization of amino acid interactions between the ligand and target was performed using Maestro and Discovery studio visualizer to analyze the type of interactions. Compounds displaying high binding affinities were subjected for analysis of pharmacokinetic studies, later molecular dynamics (MD) and MM-PBSA studies was conducted over selected compounds using GROMACS. Rutin and Isorhamnetin-3-O-rutinoside, both flavonoids thoroughly studied for their medicinal properties showcased strong interactions and the highest binding affinity of −8.9 ​kcal/mol with the Mpro. The binding energy calculated employing MM-PBSA for Rutin and Isorhamnetin-3-O-rutinoside were −86.832 ​kJ/mol and −72.984 ​kJ/mol, respectively. The overall studies revealed that Rutin and Isorhamnetin-3-O-rutinoside are portenial in inhibiting the SARS-CoV-2 Mpro and can be validated through in-vitro and in-vivo studies.

Computational and network pharmacology analysis of bioflavonoids as possible natural antiviral compounds in COVID-19.

Bioflavonoids are the largest group of plant-derived polyphenolic compounds with diverse biological potential and have also been proven efficacious in the treatment of Severe Acute Respiratory Syndrome (SARS) and Middle East Respiratory Syndrome (MERS). The present investigation validates molecular docking, simulation, and MM-PBSA studies of fifteen bioactive bioflavonoids derived from plants as a plausible potential antiviral in the treatment of COVID-19. Molecular docking studies for 15 flavonoids on the three SARS CoV-2 proteins, non-structural protein-15 Endoribonuclease (NSP15), the receptor-binding domain of spike protein (RBD of S protein), and main protease (Mpro/3CLpro) were performed and selected protein-ligand complexes were subjected to Molecular Dynamics simulations. The molecular dynamics trajectories were subjected to free energy calculation by the MM-PBSA method. All flavonoids were further assessed for their effectiveness as adjuvant therapy by network pharmacology analysis on the target proteins. The network pharmacology analysis suggests the involvement of selected bioflavonoids in the modulation of multiple signaling pathways like p53, FoxO, MAPK, Wnt, Rap1, TNF, adipocytokine, and leukocyte transendothelial migration which plays a significant role in immunomodulation, minimizing the oxidative stress and inflammation. Molecular docking and molecular dynamics simulation studies illustrated the potential of glycyrrhizic acid, amentoflavone, and mulberroside in inhibiting key SARS-CoV-2 proteins and these results could be exploited further in designing future ligands from natural sources.

Structure based peptide design, molecular dynamics and MM-PBSA studies for targeting C terminal dimerization of NFAT5 DNA binding domain

NFAT5 as a transcription factor with an established role in osmotic stress response, has also been revealed to be active under numerous settings, including pathological conditions such as diabetic microvascular complications, chronic arthritis and cancer. Despite these links, current strategies for downregulating NFAT5 activity only relies on indirect modulators, not directly targeting NFAT5, itself. With this study, through using a computational approach, an original peptide was explored to directly target C terminal dimerization of NFAT5 RHR, located in its DNA binding domain. At first, homodimeric NFAT5 RHR bound to its consensus DNA was used for prediction of a preliminary peptide sequence. Possible amino acid replacements for this preliminary peptide were predicted for optimization, which was followed by addition of a cell penetrating peptide sequence. These attempts yielded a small peptide library, which was further investigated for peptide affinities towards C terminal of NFAT5 RHR through molecular docking, 50 ns and 250 ns molecular dynamics simulations, followed by estimation of MM-PBSA based relative binding free energies. Results indicated that after receiving mutations on the preliminary peptide sequence for optimization, a unique peptide could target C terminal dimerization region of NFAT5 RHR through using its cell penetrating peptide sequence. In conclusion, this is the first study presenting computational evidence on identification of a novel peptide capable of directly targeting NFAT5 dimerization. Besides, future implications of these observations were also discussed in terms of methodology and possible applications.

Characterization of dye-decolorizing peroxidase from Bacillus subtilis

The dye-decolorizing peroxidases (DyPs) belong to a unique heme peroxidase family for their biotechnological potential to detoxify synthetic dyes. In this work, we have biochemically and structurally characterized the dye-decolorizing peroxidase from Bacillus subtilis (BsDyP). The biochemical studies of BsDyP demonstrate that pH 4.0 is optimum for the oxidation of malachite green (MG) and methyl violet (MV). However, it oxidizes the MG with higher catalytic efficiency (kcat/Km = 6.3 × 102 M−1s−1), than MV (kcat/Km = 5.0 × 102 M−1s−1). While reactive black 5 (RB5) is oxidized at pH 3.0 with the catalytic efficiency of kcat/Km = 3.6 × 102 M−1s−1. The calculated thermodynamic parameters by isothermal titration calorimetry (ITC) reveal the feasibility and spontaneity of dyes binding with BsDyP. Further, the crystal structures of a HEPES bound and unbound of BsDyP provide insight into the probable binding sites of the substrates. In BsDyP-HEPES bound structure, the HEPES-1 molecule is found in the heme cavity at the γ-edge, and another HEPES-2 molecule is bound ~16 Å away from the heme that is fenced by Ile231, Arg234, Ser235, Asp239, Glu334, and surface-exposed Tyr335 residues. Furthermore, the molecular docking, simulation, and MMPBSA studies support the binding of dyes at both the sites of BsDyP and produce lower-energy stable BsDyP-dyes complexes. Here, the BsDyP study allows the identification of its two potential binding sites and shows the oxidation of a variety of dyes. Structural and functional insight of BsDyP will facilitate its engineering for the improved decolorization of dyes.

Molecular docking, simulation and MM-PBSA studies of nigella sativa compounds: a computational quest to identify potential natural antiviral for COVID-19 treatment

Nigella sativa or black seed is used as a medicinal plant around the globe. Oil and seeds have a long tradition of folklore use in various medicinal and food systems. The conventional therapeutic use of Nigella sativa, in different ways, has been reported in several studies to treat different diseases including influenza, headache, hypertension, diabetes, inflammation, eczema, fever, cough, asthma, bronchitis, and fever. Based on previously reported potential therapeutic uses of N. sativa compounds, and keeping in mind the dire need of time for the development of potent antiviral, a combined docking, ADMET properties calculation, molecular dynamics, and MM-PBSA approaches were applied in the current study to check the therapeutic potentials of N. sativa chief constituents against COVID-19. Among the studied compounds, we found that dithymoquinone (DTQ), with binding affinity of −8.6 kcal/mol compared to a positive control (chloroquine, −7.2 kcal/mol) , has the high potential of binding at SARS-CoV-2:ACE2 interface and thus could be predicted as a plausible inhibitor to disrupt viral-host interactions. Molecular dynamics simulation of 100 ns well complemented binding affinity of the compound and revealed strong stability of DTQ at the docked site. Additionally, MM-PBSA also affirms the docking results. Compound DTQ of the present study, if validated in wet lab experiments, could be used to treat COVID-19 and could serve as a lead in the future for development of more effective natural antivirals against COVID-19. Communicated by Ramaswamy H. Sarma

Mechanistic explanation of structural and functional changes induced by methionine mutation in G-CSF protein

Abstract Potency of a biotherapeutic is influenced by its binding affinity to the designated receptor., This study examines how substitution of methionine residues with Leucine/Alanine in granulocyte colony stimulating factor (G-CSF) brings about changes in its affinity to G-CSF receptor. To this effect, two mutants M1 (M122A, M127A and M138A) and M2 (M122L, M127L and M138L) were designed computationally and their stability and functional attributes were compared with the native GCSF. The simulation results show that although all atom root mean square deviations during 100 ns window were within 2–3 A for the mutants as well as the native G-CSF, yet some local conformational changes were evident. Amino acid residues in the stretch 125–145 exhibited shift in their secondary structure composition from predominantly turn form to coiled and gamma turn forms. In addition, increase in the helical content in the amino acid stretch 60–70 was also observed when compared to native G-CSF. Docking and MMPBSA studies of the GCSF variants with G-CSFR show evidence of increased number of interactions with G-CSFR owing to the observed conformational perturbations in the variants. Experimental data on the stability and the affinity of these variants are in agreement with the in-silico predictions. Thus, we show that careful substitution of key residues in G-CSF potentially presents avenues for improving the stability as well therapeutic prospects of the molecule via enhancement in their affinities to the pertinent receptors.

Novel fluorescent antifolates that target folate receptors α and β: Molecular dynamics and density functional theory study

Nine novel fluorescent antifolates, 1–9, were designed and docked with FRα and FRβ. The binding energies of the bound complexes were determined by molecular docking and MM-PBSA studies. The structural properties of the complexes FR-FOL, FR-7, FR-8 and FR-9 were analyzed in detail via molecular docking and molecular dynamics studies. We further calculated the root mean square displacement and root mean square fluctuation of the bound complexes using molecular dynamics simulations. Since compounds 7, 8 and 9 are promising candidate in distinguishing FRα from FRβ, the hydrogen bond properties of complexes FRα-7, FRα-8 and FRα-9 were studied by a dispersion complemented density functional tight-binding method. The purpose of this study is to provide a rationale for the design of novel fluorescent antifolates targeted with FRα and FRβ.

Assessing the ligand selectivity of sphingosine kinases using molecular dynamics and MM-PBSA binding free energy calculations.

The dynamic balance of sphingolipids plays a crucial role in diverse biological processes such as mitogenesis, cell migration and angiogenesis. Sphingosine kinases (SKs) including SK1 and SK2 phosphorylate sphingosine to sphingosine 1-phosphate (S1P), and control the critical balance. SK1 overexpression was reported to increase cell survival and proliferation. Although several SK1 selective inhibitors have been reported, detailed analysis toward their selectivity to understand the molecular mechanism has not been performed to our knowledge. Herein, the crystal structure of SK1 and a homology model of SK2 were used to dock five inhibitors (1, 2, 3, 4 and 5). Protein-ligand complexes were then subjected to a molecular dynamics study and MM-PBSA binding free energy calculations. By analyzing the binding model of these inhibitors, we found that residues ILE170, PHE188 and THR192 in SK1 significantly contribute a favorable binding energy to the selectivity.