Ligand-Based Virtual Screening, Molecular Docking, Molecular Dynamics, and MM-PBSA Calculations towards the Identification of Potential Novel Ricin Inhibitors

From antibiotic to antiviral: computational screening reveals a multi-targeting antibiotic from Streptomyces spp. against Nipah virus fusion proteins

Nipah Virus is a re-emerging zoonotic paramyxovirus that poses a significant threat to both swine industry and human health. The pursuit of potential antiviral agents with both preventive and therapeutic properties holds promise for targeting such viruses. To expedite this search, leveraging computational biology is essential. Streptomyces is renowned for its capacity to produce large and diverse metabolites with promising bioactivities. In the current study, we conducted a comprehensive structure-based virtual screening of 6524 Streptomyces spp. metabolites sourced from the StreptomeDB database to evaluate their potential inhibitory effects on three Nipah virus fusion (NiVF) protein conformations: NiVF pre-fusion 1-mer (NiVF-1mer), pre-fusion 3-mer (NiVF-3mer), and NiVF post-fusion (NiVF-PoF). Prior to virtual screening, the drug-likeness of Streptomyces spp. compounds was profiled using ADMET properties. From the 913 ADMET-filtered compounds, the subsequent targeted and confirmatory blind docking analysis revealed that S896 or virginiamycin M1, a known macrolide antibiotic, showed a maximum binding affinity with the NiVF proteins, suggesting a multi-targeting inhibitory property. In addition, the 200-ns molecular dynamics simulation and MM/PBSA analyses revealed stable and strong binding affinity between the NiVF-S896 complexes, indicating favorable interactions between S896 and the target proteins. These findings suggest the potential of virginiamycin M1, an antibiotic, as a promising multi-targeting antiviral drug. However, in vitro and in vivo experimental validations are necessary to assess their safety and efficacy.

The Search for Antidotes Against Ricin

The castor plant (Ricinus communis) is primarily known for its seeds, which contain a unique fatty acid called ricinoleic acid with several industrial and commercial applications. Castor seeds also contain ricin, a toxin considered a chemical and biological warfare agent. Despite years of investigation, there is still no effective antidote or vaccine available. However, some progress has been made, and the development of an effective treatment may be on the horizon. To provide an updated overview of this issue, we have conducted a comprehensive review of the literature on the current state of research in the fight against ricin. This review is based on the reported research and aims to address the challenges faced by researchers, as well as highlight the most successful cases achieved thus far. Our goal is to encourage the scientific community to continue their efforts in this critical search.

Açaí (Euterpe oleracea Mart.) Seed Oil Exerts a Cytotoxic Role over Colorectal Cancer Cells: Insights of Annexin A2 Regulation and Molecular Modeling

Açaí, Euterpe oleracea Mart., is a native plant from the Amazonian and is rich in several phytochemicals with anti-tumor activities. The aim was to analyze the effects of açaí seed oil on colorectal adenocarcinoma (ADC) cells. In vitro analyses were performed on CACO-2, HCT-116, and HT-29 cell lines. The strains were treated with açaí seed oil for 24, 48, and 72 h, and cell viability, death, and morphology were analyzed. Molecular docking was performed to evaluate the interaction between the major compounds in açaí seed oil and Annexin A2. The viability assay showed the cytotoxic effect of the oil in colorectal adenocarcinoma cells. Acai seed oil induced increased apoptosis in CACO-2 and HCT-116 cells and interfered with the cell cycle. Western blotting showed an increased expression of LC3-B, suggestive of autophagy, and Annexin A2, an apoptosis regulatory protein. Molecular docking confirmed the interaction of major fatty acids with Annexin A2, suggesting a role of açaí seed oil in modulating Annexin A2 expression in these cancer cell lines. Our results suggest the anti-tumor potential of açaí seed oil in colorectal adenocarcinoma cells and contribute to the development of an active drug from a known natural product.

In Vitro and In Silico Studies of Kinase Inhibitor of MAPK3 Protein to Determine Leishmania martiniquensis Treatment

PURPOSE

Leishmaniasis is a parasitic disease transmitted by the bite of the phlebotomine female sand fly. Currently, no reported effective vaccines are available for the treatment of leishmaniasis; consequently, restricting this disease completely depends on controlling its transmission. Mitogen-activated protein (MAP) kinases have been reported to be involved in the regulation of the flagellum length and hence play an important role in disease transmission, especially the MAPK3 protein. Therefore, the current work focused on identifying approved drugs that can inhibit the MAPK3 protein.

METHODS

First, the recombinant plasmid (pET28b( +) MAPK3) was cloned into E. coli strain BL21 using the heatshock method. Afterward, E. coli was induced using IPTG, and cells were harvested for protein purification in the next step. After that, the MAPK3 protein was purified using Ni-NTA column. Then, the inhibition kinase activity of the purified MAPK3 protein was performed using an ADP-Glo^™ Kinase Assay kit. Furthermore, the cytotoxicity of Leishmania cells were detected by alamarBlue^™ Cell Viability Reagent. Finally, the binding affinity within the binding site of MAPK3 protein was performed by computational methods.

RESULTS

Purification of the MAPK3 protein was done using an Ni-NTA column and a protein band was identified at the expected 44 kDa molecular weight. Afterward, the ability of commercial drugs (afatinib and lapatinib) to inhibit the purified MAPK3 kinase activity was performed using an ADP-Glo^™ Kinase Assay kit. The half-maximal inhibitory concentrations (IC₅₀) of two drugs inhibited the MAPK3 protein within the same range of IC₅₀ values (3.27 and 2.22 µM for afatinib and lapatinib, respectively). Furthermore, the cytotoxicity assay of compounds toward the extracellular promastigote and intracellular amastigote stages was investigated using alamarBlue^™ Cell Viability Reagent. The results showed that both drugs were more efficient against extracellular promastigotes and intracellular amastigotes of both Leishmania donovani and Leishmania martiniquensis. Finally, the molecular dynamics simulation (MD) was performed to study the intermolecular interactions of both drugs with MAPK3 protein. From 100 ns molecular dynamics simulation, the structural stability of both drugs in a complex with MAPK3 was quite stable.

CONCLUSION

This work was suggesting that afatinib and lapatinib act as MAPK3 inhibitors and might be developed for leishmaniasis treatment.

Intrinsic Dynamics of the ClpXP Proteolytic Machine Using Elastic Network Models

ClpXP complex is an ATP-dependent mitochondrial matrix protease that binds, unfolds, translocates, and subsequently degrades specific protein substrates. Its mechanisms of operation are still being debated, and several have been proposed, including the sequential translocation of two residues (SC/2R), six residues (SC/6R), and even long-pass probabilistic models. Therefore, it has been suggested to employ biophysical-computational approaches that can determine the kinetics and thermodynamics of the translocation. In this sense, and based on the apparent inconsistency between structural and functional studies, we propose to apply biophysical approaches based on elastic network models (ENM) to study the intrinsic dynamics of the theoretically most probable hydrolysis mechanism. The proposed models ENM suggest that the ClpP region is decisive for the stabilization of the ClpXP complex, contributing to the flexibility of the residues adjacent to the pore, favoring the increase in pore size and, therefore, with the energy of interaction of its residues with a larger portion of the substrate. It is predicted that the complex may undergo a stable configurational change once assembled and that the deformability of the system once assembled is oriented, to increase the rigidity of the domains of each region (ClpP and ClpX) and to gain flexibility of the pore. Our predictions could suggest under the conditions of this study the mechanism of the interaction of the system, of which the substrate passes through the unfolding of the pore in parallel with a folding of the bottleneck. The variations in the distance calculated by molecular dynamics could allow the passage of a substrate with a size equivalent to ∼3 residues. The theoretical behavior of the pore and the stability and energy of binding to the substrate based on ENM models suggest that in this system, there are thermodynamic, structural, and configurational conditions that allow a possible translocation mechanism that is not strictly sequential.

Flavanones: A potential natural inhibitor of the ATP binding site of PknG of Mycobacterium tuberculosis

Over the years, Mycobacterium tuberculosis has been one of the major causes of death worldwide. As several clinical isolates of the bacteria have developed drug resistance against the target sites of the current therapeutic agents, the development of a novel drug is the pressing priority. According to recent studies on Mycobacterium tuberculosis, ATP binding sites of Mycobacterium tuberculosis serine/threonine protein kinases (MTPKs) have been identified as the new promising drug target. Among the several other protein kinases (PKs), Protein kinase G (PknG) was selected for the study because of its crucial role in modulating bacterium's metabolism to survive in host macrophages. In this work, we have focused on the H37Rv strain of Mycobacterium tuberculosis. A list of 477 flavanones obtained from the PubChem database was docked one by one against the crystallized and refined structure of PknG by in-silico techniques. Initially, potential inhibitors were narrowed down by preliminary docking. Flavanones were then selected using binding energies ranging from -7.9 kcal.mol^-1 to -10.8 kcal.mol^-1. This was followed by drug-likeness prediction, redocking analysis, and molecular dynamics simulations. Here, we have used experimentally confirmed drug AX20017 as a reference to determine candidate compounds that can act as potential inhibitors for PknG. PubChem165506, PubChem242065, PubChem688859, PubChem101367767, PubChem3534982, and PubChem42607933 were identified as possible target site inhibitors for PknG with a desirable negative binding energy of -8.1, -8.3, -8.4, -8.8, -8.6 and -7.9 kcal.mol^-1 respectively. Communicated by Ramaswamy H. Sarma.

Modeling studies on the role of vitamins B1 (thiamin), B3 (nicotinamide), B6 (pyridoxamine), and caffeine as potential leads for the drug design against COVID-19

In response to the COVID-19 pandemic, and the lack of effective and safe antivirals against it, we adopted a new approach in which food supplements with vital antiviral characteristics, low toxicity, and fast excretion have been targeted. The structures and chemical properties of the food supplements were compared to the promising antivirals against SARS-COV-2. Our goal was to exploit the food supplements to mimic the topical antivirals' functions but circumventing their severe side effects, which has limited the necessary dosage needed to exhibit the desired antiviral activity. On this line, after a comparative structural analysis of the chemicals mentioned above, and investigation of their potential mechanisms of action, we selected caffeine and some compounds of the vitamin B family and further applied molecular modeling techniques to evaluate their interactions with the RDB domain of the Spike protein of SARS-CoV-2 (SC2Spike) and its corresponding binding site on human ACE-2 (HssACE2). Our results pointed to vitamins B1 and B6 in the neutral form as potential binders to the HssACE2 RDB binding pocket that might be able to impair the SARS-CoV-2 mechanism of cell invasion, qualifying as potential leads for experimental investigation against COVID-19.

Theoretical Investigation of Repurposed Drugs Potentially Capable of Binding to the Catalytic Site and the Secondary Binding Pocket of Subunit A of Ricin

Recently, we reported a library of 82 compounds, selected from different databanks through virtual screening and docking studies, and pointed to 6 among them as potential repurposed dual binders to both the catalytic site and the secondary binding pockets of subunit A of ricin (RTA). Here, we report additional molecular modeling studies of an extended list of compounds from the original library. Rounds of flexible docking followed by molecular dynamics simulations and further rounds of MM-PBSA calculations using a more robust protocol, enabled a better investigation of the interactions of these compounds inside RTA, the elucidation of their dynamical behaviors, and updating the list of the most important residues for the ligand binding. Four compounds were pointed as potential repurposed ricin inhibitors that are worth being experimentally investigated.

Theoretical assessment of the performances of commercial oximes on the reactivation of acetylcholinesterase inhibited by the nerve agent A-242 (novichok)

The nerve agents of the A-series are relatively recent chemical weapons with no antidote available yet. Once inside the human body, those chemicals act similarly to the classic nerve agents, by binding to the catalytic residue Serine 203 (Ser203) of human acetylcholinesterase (HssAChE) and thus preventing the proper function of this enzyme. However, there is no experimental evidence yet if the current antidotes for intoxication by nerve agents are also capable of restoring AChE inhibited by the nerve agents of the A-series. In order to launch some light on this issue, we used computational techniques (molecular docking, molecular dynamics and MM-PBSA interaction energy calculations) to assess the performances of the four currently available commercial oximes (2-PAM, HI-6, obidoxime and trimedoxime) when in contact with HssAChE inhibited by the agent A-242. Based on the near-attack conformation (NAC) criterion, our results suggest that the commercial oximes would have limited efficacy to reactivate the enzyme since they are not able to properly approach the adduct Ser203-A-242. Among those oximes, trimedoxime seems to be the most promising, since it showed lower values of energy in the MM-PBSA calculations, a higher stability inside the catalytic anionic center (CAS) of HssAChE, and was able to adopt a position closer to the NAC that could enable the reactivation mechanism.

Immunotoxins and nanobody-based immunotoxins: review and update

Immunotoxins (ITs) are protein-based drugs that compose of targeting and cytotoxic moieties. After binding the IT to the specific cell-surface antigen, the IT internalises into the target cell and kills it. Targeting and cytotoxic moieties usually include monoclonal antibodies and protein toxins with bacterial or plant origin, respectively. ITs have been successful in haematologic malignancies treatment. However, ITs penetrate poorly into solid tumours because of their large size. Use of camelid antibody fragments known as nanobodies (Nbs) as a targeting moiety may overcome this problem. Nbs are the smallest fragment of antibodies with excellent tumour tissue penetration. The ability to recognise cryptic (immuno-evasive) target antigens, low immunogenicity, and high-affinity are other fundamental characteristics of Nbs that make them suitable candidates in targeted therapy. Here, we reviewed and discussed the structure and function of ITs, Nbs, and nanobody-based ITs. To gain sound insight into the issue at hand, we focussed on nanobody-based ITs.

Recent Developments in Free Energy Calculations for Drug Discovery

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.

A Multi-target Drug Designing for BTK, MMP9, Proteasome and TAK1 for the Clinical Treatment of Mantle Cell Lymphoma

BACKGROUND

Mantle cell lymphoma (MCL) is a type of non-Hodgkin lymphoma characterized by the mutation and overexpression of the cyclin D1 protein by the reciprocal chromosomal translocation t(11;14)(q13:q32).

AIM

The present study aims to identify potential inhibition of MMP9, Proteasome, BTK, and TAK1 and determine the most suitable and effective protein target for the MCL.

METHODOLOGY

Nine known inhibitors for MMP9, 24 for proteasome, 15 for BTK and 14 for TAK1 were screened. SB-3CT (PubChem ID: 9883002), oprozomib (PubChem ID: 25067547), zanubrutinib (PubChem ID: 135565884) and TAK1 inhibitor (PubChem ID: 66760355) were recognized as drugs with high binding capacity with their respective protein receptors. 41, 72, 102 and 3 virtual screened compounds were obtained after the similarity search with compound (PubChem ID:102173753), PubChem compound SCHEMBL15569297 (PubChem ID:72374403), PubChem compound SCHEMBL17075298 (PubChem ID:136970120) and compound CID: 71814473 with best virtual screened compounds.

RESULT

MMP9 inhibitors show commendable affinity and good interaction profile of compound holding PubChem ID:102173753 over the most effective established inhibitor SB-3CT. The pharmacophore study of the best virtual screened compound reveals its high efficacy based on various interactions. The virtual screened compound's better affinity with the target MMP9 protein was deduced using toxicity and integration profile studies.

CONCLUSION

Based on the ADMET profile, the compound (PubChem ID: 102173753) could be a potent drug for MCL treatment. Similar to the established SB-3CT, the compound was non-toxic with LD50 values for both the compounds lying in the same range.

Discovery of RTA ricin subunit inhibitors: a computational study using PM7 quantum chemical method and steered molecular dynamics

Ricin is a potent toxin derived from the castor bean plant and comprises two subunits, RTA and RTB. Because of its cytotoxicity, ricin has alarmed world authorities for its potential use as a chemical weapon. Ricin also affects castor bean agribusiness, given the risk of animal and human poisoning. Over the years, many groups attempted to propose small-molecules that bind to the RTA active site, the catalytic chain. Despite such efforts, there is still no effective countermeasure against ricin poisoning. The computational study carried out in the present work renews the discussion about small-molecules that may inhibit this toxin. Here, a structure-based virtual screening protocol capable of discerning active RTA inhibitors from inactive ones was performed to screen over 2 million compounds from the ZINC database to find novel scaffolds that strongly bind into the active site of the RTA. Besides, a novel score method based on ligand undocking force profiles and semi-empirical quantum chemical calculations provided insights into the rescore of docking poses. Summing up, the filtering steps pointed out seven main compounds, with the SCF00-451 as a promising candidate to inhibit the killing activity of such potent phytotoxin.

Searching for potential drugs against SARS-CoV-2 through virtual screening on several molecular targets

The acute respiratory syndrome caused by the SARS-CoV-2, known as COVID-19, has been ruthlessly tormenting the world population for more than six months. However, so far no effective drug or vaccine against this plague have emerged yet, despite the huge effort in course by researchers and pharmaceutical companies worldwide. Willing to contribute with this fight to defeat COVID-19, we performed a virtual screening study on a library containing Food and Drug Administration (FDA) approved drugs, in a search for molecules capable of hitting three main molecular targets of SARS-CoV-2 currently available in the Protein Data Bank (PDB). Our results were refined with further molecular dynamics (MD) simulations and MM-PBSA calculations and pointed to 7 multi-target hits which we propose here for experimental evaluation and repurposing as potential drugs against COVID-19. Additional rounds of docking, MD simulations and MM-PBSA calculations with remdesivir suggested that this compound can also work as a multi-target drug against SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

Identification of novel potential ricin inhibitors by virtual screening, molecular docking, molecular dynamics and MM-PBSA calculations: a drug repurposing approach

Ricin is a potent cytotoxin with no available antidote. Its catalytic subunit, RTA, damages the ribosomal RNA (rRNA) of eukaryotic cells, preventing protein synthesis and eventually leading to cell death. The combination between easiness of obtention and high toxicity turns ricin into a potential weapon for terrorist attacks, urging the need of discovering effective antidotes. On this context, we used computational techniques, in order to identify potential ricin inhibitors among approved drugs. Two libraries were screened by two different docking algorithms, followed by molecular dynamics simulations and MM-PBSA calculations in order to corroborate the docking results. Three drugs were identified as potential ricin inhibitors: deferoxamine, leucovorin and plazomicin. Our calculations showed that these compounds were able to, simultaneously, form hydrogen bonds with residues of the catalytic site and the secondary binding site of RTA, qualifying as potential antidotes against intoxication by ricin.Communicated by Ramaswamy H. Sarma.