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

Phospholipase A2 (PLA2) is the most abundant protein found in snake venom. PLA2 induces a variety of pharmacological effects such as neurotoxicity, myotoxicity and cardiotoxicity as well as anticoagulant, hemolytic, anti-platelet, hypertensive, hemorrhagic and edema inducing effects. In this study, the three dimensional structure of PLA2 of Naja sputatrix (Malayan spitting cobra) was modeled by I-TASSER, SWISS-MODEL, PRIME and MODELLER programs. The best model was selected based on overall stereo-chemical quality. Further, molecular dynamics simulation was performed to know the stability of the modeled protein using Gromacs software. Average structure was generated during the simulation period of 10 ns. High throughput virtual screening was employed through different databases (Asinex, Hit finder, Maybridge, TOSLab and ZINC databases) against PLA2. The top seven compounds were selected based on the docking score and free energy binding calculations. These compounds were analyzed by quantum polarized ligand docking, induced fit docking and density functional theory calculation. Furthermore, the stability of lead molecules in the active site of PLA2 was employed by MD simulation. The results show that selected lead molecules were highly stable in the active site of PLA2.

Neutralization of the haemorrhagic activities of viperine snake venoms and venom metalloproteinases using synthetic peptide inhibitors and chelators

Non-small cell lung cancer (NSCLC) is an obscure disease whose incidence is increasing worldwide day by day, and PI3Kα is one of the major targets for cell proliferation due to the mutation. Since PI3K is a class of kinase enzyme, and no in silico research has been performed on the inhibition of PI3Kα mutation by small molecules, we have selected the protein kinase inhibitor database and performed the energy minimization process by ligand preparation. The key objective of this research is to identify the potential hits from the protein kinase inhibitor library and further to perform lead optimization by a molecular docking and dynamics approach. And so, the protein was selected (PDB ID: 4JPS), having a unique inhibitor and a specific binding pocket with amino acid residue for the inhibition of kinase activity. After the docking protocol validation, structure-based virtual screening by molecular docking and MMGBSA binding affinity calculations were performed and a total of ten hits were reported. Detailed analysis of the best scoring molecules was performed with ADMET analysis, induced fit docking (IFD) and molecular dynamics (MD) simulation. Two molecules – 6943 and 34100 – were considered lead molecules and showed better results than the PI3K inhibitor Copanlisib in the docking assessment, ADMET analysis, and molecular dynamics simulation. Furthermore, the synthetic accessibility of the two compounds – 6943 and 34100 – was investigated using SwissADME, and the two lead molecules are easier to synthesize than the PI3K inhibitor Copanlisib. Computational drug discovery tools were used for identification of kinase inhibitors as anti-cancer agents for NSCLC in the present research.

Quorum sensing (QS) plays an important role in the biofilm formation, production of virulence factors and stress responses in Vibrio harveyi. Therefore, interrupting QS is a possible approach to modulate bacterial behavior. In the present study, three docking protocols, such as Rigid Receptor Docking (RRD), Induced Fit Docking (IFD), and Quantum Polarized Ligand Docking (QPLD) were used to elucidate the binding mode of boronic acid derivatives into the binding pocket of LuxP protein in V. harveyi. Among the three docking protocols, IFD accurately predicted the correct binding mode of the studied inhibitors. Molecular dynamics (MD) simulations of the protein-ligand complexes indicates that the inter-molecular hydrogen bonds formed between the protein and ligand complex remains stable during the simulation time. Pharmacophore and shape-based virtual screening were performed to find selective and potent compounds from ChemBridge database. Five hit compounds were selected and subjected to IFD and MD simulations to validate the binding mode. In addition, enrichment calculation was performed to discriminate and separate active compounds from the inactive compounds. Based on the computational studies, the potent Bicyclo [2.2.1] hept-5-ene-2,3-dicarboxylic acid-2,6-dimethylpyridine 1-oxide (ChemBridge_5144368) was selected for in vitro assays. The compound exhibited dose dependent inhibition in bioluminescence and also inhibits biofilm formation in V. harveyi to the level of 64.25%. The result from the study suggests that ChemBridge_5144368 could serve as an anti-quorum sensing molecule for V. harveyi.

Collagen binding is a key factor for the hemorrhagic activity of snake venom metalloproteinases

c-Jun-NH2 terminal kinases (JNKs) come under a class of serine/threonine protein kinases and are encoded by three genes, namely JNK1, JNK2 and JNK3. Human JNK1 is a cytosolic kinase belonging to mitogen-activated protein kinase (MAPK) family, which plays a major role in intracrinal signal transduction cascade mechanism. Overexpressed human JNK1, a key kinase interacts with other kinases involved in the etiology of many cancers, such as skin cancer, liver cancer, breast cancer, brain tumors, leukemia, multiple myeloma and lymphoma. Thus, to unveil a novel human JNK1 antagonist, receptor-based pharmacophore modeling was performed with the available eighteen cocrystal structures of JNK1 in the protein data bank. Eighteen e-pharmacophores were generated from the 18 cocrystal structures. Four common e-pharmacophores were developed from the 18 e-pharmacophores, which were used as three-dimensional (3D) query for shape-based similarity screening against more than one million small molecules to generate a JNK1 ligand library. Rigid receptor docking (RRD) performed using GLIDE v6.3 for the 1683 compounds from in-house library and 18 cocrystal ligands with human JNK1 from lower stringency to higher stringency revealed 17 leads. Further to derive the best leads, dock complexes obtained from RRD were studied further with quantum-polarized ligand docking (QPLD), induced fit docking (IFD) and molecular mechanics/generalized Born surface area (MM-GBSA). Four leads have showed lesser binding free energy and better binding affinity towards JNK1 compared to 18 cocrystal ligands. Additionally, JNK1–lead1 complex interaction stability was reasserted using 50 ns MD simulations run and also compared with the best resolute cocrystal structure using Desmond v3.8. Thus, the results obtained from RRD, QPLD, IFD and MD simulations indicated that lead1 might be used as a potent antagonist toward human JNK1 in cancer therapeutics.

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.

Gemini surfactant corrosion inhibitor (CI) is one type of CI mainly used in mitigating corrosion in the complex system of oil/gas production industries. Computer modeling methods such as density functional theory (DFT) calculation and molecular dynamic (MD) simulation are required to develop new CI molecules focusing on their application condition as a prediction or screening process before the physical empirical assessment. In this work, the adsorption inhibition efficiencies of two monomer surfactants (2B and H) and their respective Gemini structures with the addition of different spacers (alkyl, benzene, ester, ether, and ketone) are investigated using DFT calculation and MD simulation method in 3% sodium chloride (NaCl), and 1500 ppm acetic acid solutions. In DFT calculation, 2B-benzene molecules are assumed to have the most promising inhibition efficiency based on their high reactivity and electron-donating ability at their electron-rich benzene ring region based on the lowest bandgap energy (0.765 eV) and highest HOMO energy value (-2.879 eV), respectively. DFT calculation results correlate with the adsorption energy calculated from MD simulation, where 2B-benzene is also assumed to work better as a CI molecule with the most adsorption strength towards Fe (110) metal with the highest negative adsorption energy value (-1837.33 kJ/mol at temperature 323 K). Further, diffusion coefficient and molecular aggregation analysis in different CI concentrations through MD simulation reveals that only a small amount of Gemini surfactant CI is needed in the inhibition application compared to its respective monomer. Computer simulation methods successfully predict and screen the Gemini surfactant CI molecules that can work better as a corrosion inhibitor in acetic acid media. The amount of Gemini surfactant CI that needs to be used is also predicted. The future planning or way forward from this study will be the development of the most promising Gemini surfactant CI based on the results from DFT calculation and MD simulations.

Snake venoms are cocktails of biologically active molecules that have evolved to immobilize prey, but can also induce a severe pathology in humans that are bitten. While animal-derived polyclonal antivenoms are the primary treatment for snakebites, they often have limitations in efficacy and can cause severe adverse side effects. Building on recent efforts to develop improved antivenoms, notably through monoclonal antibodies, requires a comprehensive understanding of venom toxins. Among these toxins, snake venom metalloproteinases (SVMPs) play a pivotal role, particularly in viper envenomation, causing tissue damage, hemorrhage and coagulation disruption. One of the current challenges in the development of neutralizing monoclonal antibodies against SVMPs is the large size of the protein and the lack of existing knowledge of neutralizing epitopes. Here, we screened a synthetic human antibody library to isolate monoclonal antibodies against an SVMP from saw-scaled viper (genus Echis) venom. Upon characterization, several antibodies were identified that effectively blocked SVMP-mediated prothrombin activation. Cryo-electron microscopy revealed the structural basis of antibody-mediated neutralization, pinpointing the non-catalytic cysteine-rich domain of SVMPs as a crucial target. These findings emphasize the importance of understanding the molecular mechanisms of SVMPs to counter their toxic effects, thus advancing the development of more effective antivenoms.

Snakebite envenomation is a major public health concern, particularly in low- and middle-income regions where access to safe and effective antivenoms is limited. Traditional antivenoms, derived from immunization with crude venom, often trigger adverse reactions and lack specificity against key venom components. This study presents a bioinformatics-driven approach to design, construct, and evaluate a novel panel of recombinant multi-epitope toxin-derived immunogens targeting the most clinically significant snake venom toxin families. Five principal toxin families-three-finger toxins (3FT), phospholipase A2 (PLA2), snake venom metalloproteinases (SVMP), snake venom serine proteases (SVSP), and dendrotoxins (DDT)-were computationally analyzed to identify conserved, antigenic, non-toxic, and non-allergenic B-cell and T-cell epitopes. Multi-epitope recombinant constructs were designed using linker systems, a TLR4 agonist adjuvant, and an MITD-signal sequence to enhance immunogenicity. Structural modeling, refinement, and validation were performed using AlphaFold 3 and GalaxyRefine. Protein-TLR interactions were assessed using molecular docking with ClusPro, selecting the top-ranked pose based on the lowest energy score for initial analysis, as these often correspond to the most stable configurations. Normal mode analysis (NMA), molecular dynamics simulation (MDS), and post-simulation analysis were used to evaluate the stability of the complexes. While in silico immune simulations evaluated the immunogenic potential of the constructs. The designed multi-epitope toxin-derived immunogens were predicted to have favorable physicochemical properties, including molecular weights ranging from 33.7 to 56.37 kDa, basic isoelectric points, high thermostability (aliphatic index: 66.72-77.08), and hydrophilicity (negative GRAVY scores). Refined 3D models exhibited more than 93% residues in Ramachandran favored regions, suggesting structural reliability. Molecular docking revealed strong and stable interactions with TLR2 and TLR4, particularly in the SVMP-TLR4 complex (ΔG = -21.0 kcal/mol; Kd = 3.7 × 10-16 M). NMA, MDS, and post-simulation analysis collectively showed that 3FT immunogen complexes with TLR2/4 were the most stable and compact with coordinated motions, PLA2 and DDT displayed moderate flexibility with maintained integrity, SVSP showed intermediate instability, and SVMP, particularly with TLR2, exhibited pronounced conformational instability and dynamic disorder, highlighting clear receptor- and toxin-dependent differences in stability and collective behavior. Immune simulations theoretically predicted robust humoral and cellular immune responses, with early IgM/IgG production, expansion of B-cell and T-cell populations, and balanced cytokine profiles indicative of a safe immunogenic response. This study provides in silico evidence suggesting the potential of recombinant multi-epitope toxin-derived immunogens as a next-generation therapeutic strategy for snakebite management. The designed constructs may offer improvements in specificity, safety, and manufacturability over traditional antivenoms, providing a promising foundation for further experimental validation and clinical translation. Future refinements could incorporate cluster overlap evaluation to mitigate bias from single-pose selection, particularly for poses with overlapping scores.

Ikshusterol3-O-glucoside was isolated from Clematis gouriana Roxb. ex DC. root. A structure of the isolated compound was determined on the basis of various spectroscopic interpretations (UV, NMR, FTIR, and GC-MS-EI). This structure was submitted in the PubChem compound database (SID 249494133). SID 249494133 was carried out by density functional theory calculation to observe the chemical stability and electrostatic potential of this compound. The absorption, distribution, metabolism, and excretion property of this compound was predicted to evaluate the drug likeness and toxicity. In addition, molecular docking, quantum polarized ligand docking, prime MMGBSA calculation, and induced fit docking were performed to predict the binding status of SID 249494133 with the active site of phospholipase A2 (PLA2) (PDB ID: 1A3D). The stability of the compound in the active site of PLA2 was carried out using molecular dynamics simulation. Further, the anti-venom activity of the compound was assessed using the PLA2 assay against Naja naja (Indian cobra) crude venom. The results strongly show that Ikshusterol3-O-glucoside has a potent snake-venom neutralizing capacity and it might be a potential molecule for the therapeutic treatment for snakebites.

Epidermal growth factor receptor (EGFR) is a transmembrane protein belonging to the receptor tyrosine kinase (RTK) superfamily, reported as a promising anticancer target in treating diverse malignancies. Previous studies on microarray gene expression and methylation status in cell and animal models revealed the differential expression of EGFR at the early stages of cellular transformation. Additionally, an unpublished study of methylation analysis of EGFR gene promoters conducted in human cancer-related samples showed a several-fold increase in EGFR gene expression, suggesting epigenetic upregulation in tumors. Considering these findings, in the present study, we selected inactive (DFGout) (D: aspartic acid, F: phenylalanine, G: glycine) and active (DFGin) confirmations of EGFR to identify novel lead molecules against aberrant EGFR activity in cancer. Extra precision (XP) docking, molecular mechanics/generalized born surface area (MM/GBSA), molecular dynamics (MD) simulations, and ADME/T were performed, and the results showed that the lead 1 molecule of each target exhibited a better binding affinity and favorable stability than the existing ligands.

This study aimed to screen putative drug targets associated with biofilm formation of multidrug-resistant Acinetobacter baumannii and Pseudomonas areugenosa and prioritize carbon nano-fullerene as potential lead molecule by structure-based virtual screening. Based on the functional role, 36 and 83 genes that are involved in biofilm formation of A. baumannii and P. areugenosa respectively were selected and metabolic network was computationally constructed. The genes that lack three-dimensional structures were predicted and validated. Carbon nano-fullerene selected as lead molecule and their drug-likeliness and pharmacokinetics properties were computationally predicted. The binding potential of carbon nano-fullerene toward selected drug targets was modeled and compared with the binding of conventional drugs, doripenem, and polymyxin-B with their usual targets. The stabilities of four best-docked complexes were confirmed by molecular dynamic (MD) simulation. This study suggested that selected genes demonstrated relevant interactions in the constructed metabolic pathways. Carbon fullerene exhibited significant binding abilities to most of the prioritized targets in comparison with the binding of last-resort antibiotics and their usual target. The four best ligand–receptor interactions predicted by molecular docking revealed that stability throughout MD simulation. Notably, carbon fullerene exhibited profound binding with outer membrane protein (OmpA) and ribonuclease-HII (rnhB) of A. baumannii and 2-heptyl-4(1H)-quinolone synthase (pqsBC) and chemotaxis protein (wspA) of P. aeruginosa. Thus, the current study suggested that carbon fullerene was probably used as potential lead molecules toward selected targets of A. baumannii and P. aeruginosa and the applied aspects probably scaled up to design promising lead molecules toward these pathogens. Communicated by Ramaswamy H. Sarma.

The Arabian saw-scaled viper (Echis coloratus) is among the snakes of highest medical relevance in the Middle East and North Africa. However, to date, its venom has been investigated in a very limited number of studies, and much remains unknown regarding its compositional and functional properties. By integrating proteotranscriptomics with bioactivity profiling, we present a comprehensive transcriptome-level catalogue of E. coloratus venom components and their associated biological activities. Our analysis identified 183 venom components belonging to 17 distinct protein families. Relative toxin abundances revealed that 92% of the venom proteome is composed of C-type lectin and C-type lectin-related protein (CTL), L-amino acid oxidase (LAAO), phospholipase A2 (PLA2), snake venom serine protease (SVSP), and snake venom metalloproteinase (SVMP), with CTL and PLA2 alone accounting for 73% of the total composition. Bioassays targeting key aspects of viperid envenomation demonstrated potent protease and PLA2 activity in a concentration-dependent manner. In contrast, Factor Xa-like, plasmin-like, and haemolytic activities were negligible. Marked cytotoxicity was observed at the highest concentration tested (i.e., 25 μg/ml) in the mammalian cell lines MDCK II and Calu-3, whereas cytotoxic effects were minimal at lower concentrations. These findings highlight the complexity and potency of E. coloratus venom, and provide a valuable foundation for improving our understanding of envenomation caused by this species.

Identification of potent inhibitors against chorismate synthase of Toxoplasma gondii using molecular dynamics simulations