In silico discovery of potent inhibitors against monkeypox's major structural proteins

Exploring the chemical space for potential inhibitors against cell surface binding protein of Mpox virus using molecular fingerprint based screening approach

Mpox virus is the latest member of the Poxviridae family of which small pox virus is a member. Monekypox virus has led to thousands of infections across the globe. Poxvirus gains entry into the cell making use of glycosaminoglycans like chondroitin sulphate and heparan sulphate. The interaction of the Mpox virus protein E8L also called cell surface binding protein is crucial for host cell attachment, membrane fusion and viral entry into the host cell leading to establishment of infection thus making this protein a very attractive therapeutic target. In this study we have tried to utilize the chondroitin sulphate binding groove present in the protein and identify molecules which are structurally similar to chondroitin sulphate. These molecules can thus occupy the same pocket but with a better binding affinity than chondroitin sulphate in order to outcompete the latter molecule from binding to the E8L protein and thus prevent it from performing its function. This study may pave the way for development of highly efficient therapeutics against the Mpox virus and further curb its infective potential.Communicated by Ramaswamy H. Sarma.

Repositioning of anti-infective compounds against monkeypox virus core cysteine proteinase: a molecular dynamics study

Monkeypox virus (MPXV) core cysteine proteinase (CCP) is one of the major drug targets used to examine the inhibitory action of chemical moieties. In this study, an in silico technique was applied to screen 1395 anti-infective compounds to find out the potential molecules against the MPXV-CCP. The top five hits were selected after screening and processed for exhaustive docking based on the docked score of ≤ -9.5 kcal/mol. Later, the top three hits based on the exhaustive-docking score and interaction profile were selected to perform MD simulations. The overall RMSD suggested that two compounds, SC75741 and ammonium glycyrrhizinate, showed a highly stable complex with a standard deviation of 0.18 and 0.23 nm, respectively. Later, the MM/GBSA binding free energies of complexes showed significant binding strength with ΔGTOTAL from -21.59 to -15 kcal/mol. This report reported the potential inhibitory activity of SC75741 and ammonium glycyrrhizinate against MPXV-CCP by competitively inhibiting the binding of the native substrate.

Cracking a cancer code DNA methylation in epigenetic modification: an in-silico approach on efficacy assessment of Sri Lanka-oriented nutraceuticals

DNA methyltransferase (DNMTs) are essential epigenetic modifiers that play a critical role in gene regulation. These enzymes add a methyl group to cytosine's 5'-carbon, specifically within CpG dinucleotides, using S-adenosyl-L-methionine. Abnormal overexpression of DNMTs can alter the gene expression patterns and contribute to cancer development in the human body. Therefore, the inhibition of DNMT is a promising therapeutic approach to cancer treatment. This study was aimed to identify potential nutraceutical inhibitors from the Sri Lanka Flora database using computational methods, which provided an atomic-level description of the drug binding site and examined the interactions between nutraceuticals and amino acids of the DNMT enzyme. A series of nutraceuticals from Sri Lanka-oriented plants were selected and evaluated to assess their inhibitory effects on DNMT using absorption, distribution, metabolism, excretion and toxicity analysis, virtual screening, molecular docking, molecular dynamics simulation and trajectory analysis. Azacitidine, a DNMT inhibitor approved by the US Food and Drug Administration, was selected as a reference inhibitor. The complexes with more negative binding energies were selected and further assessed for their potency. Seven molecules were identified from 200 nutraceuticals, demonstrating significantly negative binding energies against the DNMT enzyme. Various trajectory analyses were conducted to investigate the stability of the DNMT enzyme. The results indicated that petchicine (NP#0003), ouregidione (NP#0011) and azacitidine increased the stability of the DNMT enzyme. Consequently, these two nutraceuticals showed inhibitory efficacies similar to azacitidine, making them potential candidates for therapeutic interventions targeting DNMT enzyme-related cancers. Additional bioassay testing is recommended to confirm the efficacies of these nutraceuticals and explore their applicability in clinical treatments.Communicated by Ramaswamy H. Sarma.

Bioprospecting of Meliaceae family phytomolecules for the treatment of monkeypox virus infection: a QSAR modeling and MD simulation approach

Recent monkeypox virus (MPXV) infections show the risk of MPXV transmission that persists today and the significance of surveillance and quick response methods to stop the virus's spread. Currently, the monkeypox virus infection is not specifically treated. In this study, QSAR models were designed using known inhibitors of cysteine proteinase from the vaccinia virus, where the Random Forest model and Ridge model had showed the best correlation between predicted and observed EC50. These models were used to screen Maliaceae family phytochemicals against MPXV cysteine proteinase. The compound, IMPHY010637 was detected in top 5 from both the QSAR screening models and showed best docked score (-8.6 kcal/mol) and thus selected for further investigation. Further, the IMPHY010637 showed interaction with the catalytic residue His241 of the protein as reported in earlier studies. The ADMET analysis of the compound showed the acceptable drug-like properties of IMPHY010637. However, these properties could be improved after experimental validation of protein-ligand binding. Both docked complex and poses created in 100 ns MD simulation of the protein-ligand complex showed the presence of multiple hydrogen bonds. RMSD and conformation analysis showed stable binding of IMPHY010637 with the cysteine proteinase of MPXV at its active site. Compared to the known inhibitor, IMPHY010637 showed better binding with the protein as observed by the PCA and MM/GBSA analysis. This study concluded IMPHY010637 as a potential inhibitor for the cysteine proteinase of MPXV using computational methods that could be tested in in-vitro experiments.Communicated by Ramaswamy H. Sarma.

Discovery of non-peptide GLP-1r natural agonists for enhancing coronary safety in type 2 diabetes patients

This study explores the computational discovery of non-peptide agonists targeting the Glucagon-Like Peptide-1 Receptor (GLP-1R) to enhance the safety of major coronary outcomes in individuals affected by Type 2 Diabetes. The objective is to identify novel compounds that can activate the GLP-1R pathway without the limitations associated with peptide agonists. Type 2 diabetes mellitus (T2DM) is associated with an increased risk of cardiovascular disease (CVD) and mortality, which is attributed to the accumulation of fat in organs, including the heart. Glucagon-like peptide-1 receptor agonists (GLP-1RAs) are frequently used to manage T2DM and could potentially offer cardiovascular benefits. Therefore, this study examines non-peptide agonists of GLP-1R to improve coronary safety in type 2 diabetes patients. After rigorous assessments, two standout candidates were identified, with natural compound 12 emerging as the most promising. This study represents a notable advancement in enhancing the management of coronary outcomes among individuals with type 2 diabetes. The computational methodology employed successfully pinpointed potential GLP-1R natural agonists, providing optimism for the development of safer and more effective therapeutic interventions. Although computational methodologies have provided crucial insights, realizing the full potential of these compounds requires extensive experimental investigations, crucial in advancing therapeutic strategies for this critical patient population.Communicated by Ramaswamy H. Sarma.

Role of Conformational Dynamics of Sulfotransferases SULT1A1 and SULT1A3 in Substrate Specificity

Sulfotransferases (SULTs) are phase II metabolizing enzymes catalyzing the sulfoconjugation from the co-factor 3'-Phosphoadenosine 5'-Phosphosulfate (PAPS) to a wide variety of endogenous compounds, drugs and natural products. Although SULT1A1 and SULT1A3 share 93% identity, SULT1A1, the most abundant SULT isoform in humans, exhibits a broad substrate range with specificity for small phenolic compounds, while SULT1A3 displays a high affinity toward monoamine neurotransmitters like dopamine. To elucidate the factors determining the substrate specificity of the SULT1 isoenzymes, we studied the dynamic behavior and structural specificities of SULT1A1 and SULT1A3 by using molecular dynamics (MD) simulations and ensemble docking of common and specific substrates of the two isoforms. Our results demonstrated that while SULT1A1 exhibits a relatively rigid structure by showing lower conformational flexibility except for the lip (loop L1), the loop L2 and the cap (L3) of SULT1A3 are extremely flexible. We identified protein residues strongly involved in the recognition of different substrates for the two isoforms. Our analyses indicated that being more specific and highly flexible, the structure of SULT1A3 has particularities in the binding site, which are crucial for its substrate selectivity.

A high-throughput structural dynamics approach for identification of potential agonists of FFAR4 for type 2 diabetes mellitus therapy

Diabetes mellitus is a metabolic disorder that persists as a global threat to the world. A G-protein coupled receptor (GPCR), free fatty acid receptor 4 (FFAR4), has emerged as a potential target for type 2 diabetes mellitus (T2DM) and obesity-related disorders. The current study has investigated the FFAR4, deploying 3-dimensional structure modeling, molecular docking, machine learning, and high-throughput virtual screening methods to unravel the receptor's crucial and non-crucial binding site residues. We screened four lakh compounds and shortlisted them based on binding energy, stereochemical considerations, non-bonded interactions, and pharmacokinetic profiling. Out of the screened compounds, four compounds were selected for ligand-bound simulations. The molecular dynamic simulations were carried out for 1µs for native FFAR4 and 500 ns each for complexes of FFAR4 with compound 1, compound 2, compound 3, and compound 4. Our findings showed that in addition to reported binding site residues ARG99, ARG183, and VAL98 in known agonists like TUG-891, the amino acids ARG22, ARG24, THR23, TRP305, and GLU43 were also critical binding site residues. These amino acids impart stability to the FFAR4 complexes and contribute to the stronger binding affinity of the compounds. The study also indicated that aromatic residues like PHE211 are crucial for recognizing the active site's pi-pi and C-C double bonds. Since FFAR4 is a membrane protein, the simulation studies give an insight into the mechanisms of the crucial protein-lipid and lipid-water interactions. The analysis of the molecular dynamics trajectories showed all four compounds as potential hit molecules that can be developed further into potential agonists for T2DM therapy. Amongst the four compounds, compound 4 showed relatively better binding affinity, stronger non-bonded interactions, and a stable complex.Communicated by Ramaswamy H. Sarma.

Exploration of phytochemical compounds against Marburg virus using QSAR, molecular dynamics, and free energy landscape

Marburg virus disease (MVD) is caused by the Marburg virus, a one-of-a-kind zoonotic RNA virus from the genus Filovirus. Thus, this current study employed AI-based QSAR and molecular docking-based virtual screening for identifying potential binders against the target protein (nucleoprotein (NP)) of the Marburg virus. A total of 2727 phytochemicals were used for screening, out of which the top three compounds (74977521, 90470472, and 11953909) were identified based on their predicted bioactivity (pIC50) and binding score (< - 7.4 kcal/mol). Later, MD simulation in triplicates and trajectory analysis were performed which showed that 11953909 and 74977521 had the most stable and consistent complex formations and had the most significant interactions with the highest number of hydrogen bonds. PCA (principal component analysis) and FEL (free energy landscape) analysis indicated that these compounds had favourable energy states for most of the conformations. The total binding free energy of the compounds using the MM/GBSA technique showed that 11953909 (ΔGTOTAL = - 30.78 kcal/mol) and 74977521 (ΔGTOTAL = - 30 kcal/mol) had the highest binding affinity with the protein. Overall, this in silico pipeline proposed that the phytochemicals 11953909 and 74977521 could be the possible binders of NP. This study aimed to find phytochemicals inhibiting the protein's function and potentially treating MVD.

In-vitro and computational analysis of Urolithin-A for anti-inflammatory activity on Cyclooxygenase 2 (COX-2)

Cyclooxygenase 2 (COX-2) participates in the inflammation process by converting arachidonic acid into prostaglandin G2 which increases inflammation, pain and fever. COX-2 has an active site and a heme pocket and blocking these sites stops the inflammation. Urolithin A is metabolite of ellagitannin produced from humans and animals gut microbes. In the current study, Urolithin A showed good pharmacokinetic properties. Molecular docking of the complex of Urolithin A and COX-2 revealed the ligand affinity of -7.97 kcal/mol with the ligand binding sites at TYR355, PHE518, ILE517 and GLN192 with the 4-H bonds at a distance of 2.8 Å, 2.3 Å, 2.5 Å and 1.9 Å. The RMSD plot for Urolithin A and COX-2 complex was observed to be constant throughout the duration of dynamics. A total of 3 pair of hydrogen bonds was largely observed on average of 3 simulation positions for dynamics duration of 500 ns. The MMPBSA analysis showed that active site amino acids had a binding energy of -22.0368 kJ/mol indicating that throughout the simulation the protein of target was bounded by Urolithin A. In-silico results were validated by biological assays. Urolithin A strongly revealed to exhibit anti-inflammatory effect on COX-2 with an IC50 value of 44.04 µg/mL. The anti-inflammatory capability was also depicted through reduction of protein denaturation that showed 37.6 ± 0.1 % and 43.2 ± 0.07 % reduction of protein denaturation for BSA and egg albumin respectively at 500 µg/mL. The present study, suggests Urolithin A to be an effective anti-inflammatory compound for therapeutic use.