Exploring Ligand Stability in Protein Crystal Structures Using Binding Pose Metadynamics

The potential of Chlorella spp. as antiviral source against African swine fever virus through a virtual screening pipeline

African swine fever (ASF) causes high mortality in pigs and threatens global swine production. There is still a lack of therapeutics available, with two vaccines under scrutiny and no approved small-molecule drugs. Eleven (11) viral proteins were used to identify potential antivirals in in silico screening of secondary metabolites (127) from Chlorella spp. The metabolites were screened for affinity and binding selectivity. High-scoring compounds were assessed through in silico ADMET (Absorption, Distribution, Metabolism, Excretion, Toxicity) predictions, compared to structurally similar drugs, and checked for off-target docking with prepared swine receptors. Molecular dynamics (MD) simulations determined binding stability while binding energy was measured in Molecular Mechanics - Generalized Born Surface Area (MMGBSA) or Poisson-Boltzmann Surface Area (MMPBSA). Only six (6) compounds passed until MD analyses, of which five (5) were stable after 100 ns of MD runs. Of these five compounds, only three had binding affinities that were comparable to or stronger than controls. Specifically, phytosterols 24,25-dihydrolanosterol and CID 4206521 that interact with the RNA capping enzyme (pNP868R), and ergosterol which bound to the Erv-like thioreductase (pB119L). The compounds identified in this study can be used as a theoretical basis for in vitro screening to develop potent antiviral drugs against ASFV.

Edaravone N-benzyl pyridinium derivatives: BACE-1 inhibition, kinetics and in silico binding pose determination

BACE-1 plays a pivotal role in the production of β-amyloid (Aβ) peptides, implicated in Alzheimer's Disease (AD) pathology. We previously described edaravone N-benzyl pyridinium derivatives (EBPDs) that exhibited multifunctional activity against multiple AD targets. In this study we explored the EBPDs BACE-1 inhibitory activity to potentially enhance the compounds therapeutic profile. The EBPDs exhibited moderate BACE-1 inhibitory activity (IC50 = 44.10 µM - 123.70 µM) and obtained IC50 values between 2.0 and 5.8-fold greater than resveratrol, a known BACE-1 inhibitor (IC50 = 253.20 µM), in this assay. Compound 3 was the most potent inhibitor with an IC50 of 44.10 µM and a Ki of 19.96 µM and a mixed-type mode of inhibition that favored binding in a competitive manner. Molecular docking identified crucial interactions with BACE-1 active site residues, supported by 100 ns MD simulations. The study highlighted the EBPDs therapeutic potential as BACE-1 inhibitors and multifunctional anti-AD therapeutic agents.

Collaborative GSK-University of Strathclyde doctoral research and training programmes: Transforming approaches to industry-academia engagement

A global biopharma company, GSK, and the University of Strathclyde have developed an expansive and transformative research and training partnership originating in chemistry-aligned disciplines, with subsequent extensive expansion across further areas of the company. This has opened unique approaches for the delivery of collaborative research innovations while also enhancing the professional development and learning of GSK personnel, in addition to other embedded researchers and collaborating scientists, on a pathway towards more rapid and efficient discovery of new medicines.

Discovery of an 8-oxoguanine regulator PCBP1 inhibitor by virtual screening and its synergistic effects with ROS-modulating agents in pancreatic cancer

Introduction: Drugs that target reactive oxygen species (ROS) metabolism have progressed the treatment of pancreatic cancer treatment, yet their efficacy remains poor because of the adaptation of cancer cells to high concentration of ROS. Cells cope with ROS by recognizing 8-oxoguanine residues and processing severely oxidized RNA, which make it feasible to improve the efficacy of ROS-modulating drugs in pancreatic cancer by targeting 8-oxoguanine regulators. Methods: Poly(rC)-binding protein 1 (PCBP1) was identified as a potential oncogene in pancreatic cancer through datasets of The Cancer Genome Atlas (TCGA) project and Gene Expression Omnibus (GEO). High-throughput virtual screening was used to screen out potential inhibitors for PCBP1. Computational molecular dynamics simulations was used to verify the stable interaction between the two compounds and PCBP1 and their structure-activity relationships. In vitro experiments were performed for functional validation of silychristin. Results: In this study, we identified PCBP1 as a potential oncogene in pancreatic cancer. By applying high-throughput virtual screening, we identified Compound 102 and Compound 934 (silychristin) as potential PCBP1 inhibitors. Computational molecular dynamics simulations and virtual alanine mutagenesis verified the structure-activity correlation between PCBP1 and the two identified compounds. These two compounds interfere with the PCBP1-RNA interaction and impair the ability of PCBP1 to process RNA, leading to intracellular R loop accumulation. Compound 934 synergized with ROS agent hydrogen peroxide to strongly improve induced cell death in pancreatic cancer cells. Discussion: Our results provide valuable insights into the development of drugs that target PCBP1 and identified promising synergistic agents for ROS-modulating drugs in pancreatic cancer.

Leveraging binding pose metadynamics to optimise target fishing predictions for three diverse ligands and their true targets

Computational target fishing plays an important role in target identification, particularly in drug discovery campaigns utilizing phenotypic screening. Numerous approaches exist to predict potential targets for a given ligand, but true targets may be inconsistently ranked. More advanced simulation methods may provide benefit in such cases by reranking these initial predictions. We evaluated the ability of binding pose metadynamics to improve the predicted rankings for three diverse ligands and their six true targets. Initial predictions using pharmacophore mapping showed no true targets ranked in the top 50 and two targets each ranked within the 50-100, 100-150, and 250-300 ranges respectively. Following binding pose metadynamics, ranking of true targets improved for four out of the six targets and included the highest ranked predictions overall, while rankings deteriorated for two targets. The revised rankings predicted two true targets ranked within the top 50, and one target each within the 50-100, 100-150, 150-200, and 200-250 ranges respectively. The findings of this study demonstrate that binding pose metadynamics may be of benefit in refining initial predictions from structure-based target fishing algorithms, thereby improving the efficiency of the target identification process in drug discovery efforts.

Structural basis for urate recognition and apigenin inhibition of human GLUT9

Urate, the physiological form of uric acid and a potent antioxidant in serum, plays a pivotal role in scavenging reactive oxygen species. Yet excessive accumulation of urate, known as hyperuricemia, is the primary risk factor for the development of gout. The high-capacity urate transporter GLUT9 represents a promising target for gout treatment. Here, we present cryo-electron microscopy structures of human GLUT9 in complex with urate or its inhibitor apigenin at overall resolutions of 3.5 Å and 3.3 Å, respectively. In both structures, GLUT9 exhibits an inward open conformation, wherein the substrate binding pocket faces the intracellular side. These structures unveil the molecular basis for GLUT9's substrate preference of urate over glucose, and show that apigenin acts as a competitive inhibitor by occupying the substrate binding site. Our findings provide critical information for the development of specific inhibitors targeting GLUT9 as potential therapeutics for gout and hyperuricemia.

Active and machine learning-enhanced discovery of new FGFR3 inhibitor, Rhapontin, through virtual screening of receptor structures and anti-cancer activity assessment

Introduction: This study bridges traditional remedies and modern pharmacology by exploring the synergy between natural compounds and Ceritinib in treating Non-Small Cell Lung Cancer (NSCLC), aiming to enhance efficacy and reduce toxicities. Methods: Using a combined approach of computational analysis, machine learning, and experimental procedures, we identified and analyzed PD173074, Isoquercitrin, and Rhapontin as potential inhibitors of fibroblast growth factor receptor 3 (FGFR3). Machine learning algorithms guided the initial selection, followed by Quantitative Structure-Activity Relationship (QSAR) modeling and molecular dynamics simulations to evaluate the interaction dynamics and stability of Rhapontin. Physicochemical assessments further verified its drug-like properties and specificity. Results: Our experiments demonstrate that Rhapontin, when combined with Ceritinib, significantly suppresses tumor activity in NSCLC while sparing healthy cells. The molecular simulations and physicochemical evaluations confirm Rhapontin's stability and favorable interaction with FGFR3, highlighting its potential as an effective adjunct in NSCLC therapy. Discussion: The integration of natural compounds with established cancer therapies offers a promising avenue for enhancing treatment outcomes in NSCLC. By combining the ancient wisdom of natural remedies with the precision of modern science, this study contributes to evolving cancer treatment paradigms, potentially mitigating the side effects associated with current therapies.

JD-02, a novel Hsp90 inhibitor, induces ROS/SRC axis-dependent cytoprotective autophagy in colorectal cancer cells

Heat shock protein 90 (Hsp90) is a tumor marker that accelerates cancer growth by disrupting protein homeostasis. However, concerns such as low clinical efficacy and drug resistance continue to be obstacles to the successful marketing of Hsp90 inhibitors. The cytoprotective function of autophagy has been identified as one of the mechanisms by which tumor cells gain resistance to chemotherapy. JD-02 was identified as a new Hsp90 inhibitor that suppressed colorectal cancer (CRC) growth by lowering client protein levels in vivo and in vitro. We found that JD-02 increased cellular autophagy, which inhibited apoptosis. JD-02 enhanced cytoprotective autophagy and regulated apoptotic suppression by increasing intracellular reactive oxygen species and inhibiting SRC protein levels, as demonstrated by quantitative proteomics, bioinformatic analysis, western blotting, and flow cytometry. This effect was reversed by autophagy inhibition. Therefore, due to the synergistic effects of Hsp90 and autophagy inhibitors in efficiently activating apoptotic pathways, they could potentially serve as promising therapeutic options for CRC.

Chaetoxylariones A-G: undescribed chromone-derived polyketides from co-culture of Chaetomium virescens and Xylaria grammica enabled via the molecular networking strategy

By co-culturing two endophytic fungi (Chaetomium virescens and Xylaria grammica) collected from the medicinal and edible plant Smilax glabra Roxb. and analyzing them with MolNetEnhancer module on GNPS platform, seven undescribed chromone-derived polyketides (chaetoxylariones A-G), including three pairs of enantiomer ones (2a/2b, 4a/4b and 6a/6b) and four optical pure ones (1, 3, 5 and 7), as well as five known structural analogues (8-12), were obtained. The structures of these new compounds were characterized by NMR spectroscopy, single-crystal X-ray diffraction, 13C NMR calculation and DP4+ probability analyses, as well as the comparison of the experimental electronic circular dichroism (ECD) data. Structurally, compound 1 featured an unprecedented chromone-derived sulfonamide tailored by two isoleucine-derived δ-hydroxy-3-methylpentenoic acids via the acylamide and NO bonds, respectively; compound 2 represented the first example of enantiomeric chromone derivative bearing a unique spiro-[3.3]alkane ring system; compound 3 featured a decane alkyl side chain that formed an undescribed five-membered lactone ring between C-7' and C-10'; compound 4 contained an unexpected highly oxidized five-membered carbocyclic system featuring rare adjacent keto groups; compound 7 featured a rare methylsulfonyl moiety. In addition, compound 10 showed a significant inhibition towards SW620/AD300 cells with an IC50 value of PTX significantly decreased from 4.09 μM to 120 nM, and a further study uncovered that compound 10 could obviously reverse the MDR of SW620/AD300 cells.

Structural investigation of Trypanosoma cruzi Akt-like kinase as drug target against Chagas disease

According to the World Health Organization, Chagas disease (CD) is the most prevalent poverty-promoting neglected tropical disease. Alarmingly, climate change is accelerating the geographical spreading of CD causative parasite, Trypanosoma cruzi, which additionally increases infection rates. Still, CD treatment remains challenging due to a lack of safe and efficient drugs. In this work, we analyze the viability of T. cruzi Akt-like kinase (TcAkt) as drug target against CD including primary structural and functional information about a parasitic Akt protein. Nuclear Magnetic Resonance derived information in combination with Molecular Dynamics simulations offer detailed insights into structural properties of the pleckstrin homology (PH) domain of TcAkt and its binding to phosphatidylinositol phosphate ligands (PIP). Experimental data combined with Alpha Fold proposes a model for the mechanism of action of TcAkt involving a PIP-induced disruption of the intramolecular interface between the kinase and the PH domain resulting in an open conformation enabling TcAkt kinase activity. Further docking experiments reveal that TcAkt is recognized by human inhibitors PIT-1 and capivasertib, and TcAkt inhibition by UBMC-4 and UBMC-6 is achieved via binding to TcAkt kinase domain. Our in-depth structural analysis of TcAkt reveals potential sites for drug development against CD, located at activity essential regions.

Identification of Riluzole derivatives as novel calmodulin inhibitors with neuroprotective activity by a joint synthesis, biosensor, and computational guided strategy

The development of new molecules for the treatment of calmodulin related cardiovascular or neurodegenerative diseases is an interesting goal. In this work, we introduce a novel strategy with four main steps: (1) chemical synthesis of target molecules, (2) Förster Resonance Energy Transfer (FRET) biosensor development and in vitro biological assay of new derivatives, (3) Cheminformatics models development and in vivo activity prediction, and (4) Docking studies. This strategy is illustrated with a case study. Firstly, a series of 4-substituted Riluzole derivatives 1-3 were synthetized through a strategy that involves the construction of the 4-bromoriluzole framework and its further functionalization via palladium catalysis or organolithium chemistry. Next, a FRET biosensor for monitoring Ca2+-dependent CaM-ligands interactions has been developed and used for the in vitro assay of Riluzole derivatives. In particular, the best inhibition (80%) was observed for 4-methoxyphenylriluzole 2b. Besides, we trained and validated a new Networks Invariant, Information Fusion, Perturbation Theory, and Machine Learning (NIFPTML) model for predicting probability profiles of in vivo biological activity parameters in different regions of the brain. Next, we used this model to predict the in vivo activity of the compounds experimentally studied in vitro. Last, docking study conducted on Riluzole and its derivatives has provided valuable insights into their binding conformations with the target protein, involving calmodulin and the SK4 channel. This new combined strategy may be useful to reduce assay costs (animals, materials, time, and human resources) in the drug discovery process of calmodulin inhibitors.

Acetylcholinesterase inhibition studies of alkaloid components from Crinum asiaticum var. sinicum: in vitro assessments by molecular docking and molecular dynamics simulations

Alkaloids are among the most important and best-known secondary metabolites as sources of new drugs from medicinal plants and marine organisms. A phytochemical investigation of the whole plant of Crinum asiaticum var. sinicum resulted in the isolation of seven alkaloids (1-7), including one new dimeric compound, bis-(-)-8-demethylmaritidine (1). Their structures were elucidated using NMR and HR-ESI-MS. The absolute configuration of new compound 1 was established by circular dichroism spectroscopy. All isolated compounds were evaluated for their inhibitory effects on acetylcholinesterase (AChE) activity in vitro. Among them, compound 1 exhibited the most potent AChE inhibition. Moreover, molecular docking and molecular dynamics simulations were carried out for the most active compound to investigate their binding interactions and dynamics behavior of the AChE protein-ligand complex. Therefore, compound 1 may be a potential candidate for effectively treating Alzheimer's disease.

Investigations of Ni(II)Cysteine-Tyrosine Dithiocarbamate Complex: Synthesis, Characterization, Molecular Docking, Molecular Dynamic, and Anticancer Activity on MCF-7 Breast Cancer Cell Line

OBJECTIVE

Breast cancer ranks second in terms of the highest number of cancer deaths for women worldwide and is one of the leading causes of death from cancer in women. The drug that is often used for chemotherapy is cisplatin. However, cisplatin drugs have a number of problems, including lack of selectivity, unwanted side effects, resistance, and toxicity in the body. In this work, we investigated Ni(II) cysteine-tyrosine dithiocarbamate complex against breast cancer.

METHODS

Research on the new complex compound Ni(II) cysteine-tyrosine dithiocarbamate have several stages including synthesis, characterization, in-silico and in-vitro testing of MCF-7 cells for anticancer drugs. The synthesis involved reacting cysteine, CS2, KOH and tyrosine with Mn metal. The new complex compound Ni(II) cysteine-tyrosine dithiocarbamate has been synthesized, characterized, and tested in vitro MCF-7 cells for anticancer drugs. Characterization tests such as melting point, conductivity, SEM-EDS, UV Vis, XRD, and FT-IR spectroscopy have been carried out.

RESULT

The synthesis yielded a 60,16%, conversion with a melting point of 216-218 oC and a conductivity value of 0.4 mS/cm. In vitro test results showed morphological changes (apoptosis) in MCF-7 cancer cells starting at a sample concentration of 250 µg/mL and an IC50 value of 618.40 µg/mL. Molecular docking study of Ni(II) cysteine-tyrosine dithiocarbamate complex identified with 4,4',4''-[(2R)-butane-1,1,2-triyl]triphenol - Estrogen α showing active site with acidic residue amino E323, M388, L387, G390 and I389. Hydrophobic and hydrophobic bonds are seen in Ni(II) cysteine-tyrosine dithiocarbamate - Estrogen α has a binding energy of -80.9429 kJ /mol.

CONCLUSION

there were 5 residues responsible for maintaining the ligand binding stable. The compound had significant Hbond contact intensity, however, it was not strong enough to make a significant anticancer effect. Though the synthesized compound shows low bioactivity, this research is expected to give valuable insight into the effect of molecular structure on anticancer activity.

Discovery of novel inhibitors for Pseudomonas aeruginosa lipase enzyme from in silico and in vitro studies

Pseudomonas aeruginosa is an opportunistic pathogen prone to developing drug-resistance and is a major cause of infection for burn patients and patients suffering from cystic fibrosis or are hospitalized in intensive care units. One of the virulence factors of this bacterium is the lipase enzyme that degrades the extracellular matrix of the host tissue and promotes invasion. Bromhexine is a mucolytic drug and has recently been reported to function as a competitive inhibitor of lipase with an IC50 value of 49 µM. In the present study, an attempt was made to identify stronger inhibitors from the ChEMBL database of bioactive compounds, as compared to the reference compound Bromhexine. Following docking and MD simulations, four hit compounds (N1-N4) were selected that showed promising binding modes and low RMSD values indicative of stable protein-ligand complexes. From subsequent binding pose metadynamics (BPMD) simulations, two of these (N2 and N4) stood out as more potent than Bromhexine, displaying stable interactions with residues in the catalytic site of the enzyme. Biological investigations were performed for all four compounds. Among them, the same two hit compounds were found to be the most effective binders with IC50 values of 22.1 and 27.5 µM, respectively; i.e. roughly twice as efficient as the reference Bromhexine. Taken together, our results show that these hits can be promising new candidates to use as leads for the development of drugs targeting the P. aeruginosa lipase enzyme.Communicated by Ramaswamy H. Sarma.

Comparative Structure-Based Virtual Screening Utilizing Optimized AlphaFold Model Identifies Selective HDAC11 Inhibitor

HDAC11 is a class IV histone deacylase with no crystal structure reported so far. The catalytic domain of HDAC11 shares low sequence identity with other HDAC isoforms, which makes conventional homology modeling less reliable. AlphaFold is a machine learning approach that can predict the 3D structure of proteins with high accuracy even in absence of similar structures. However, the fact that AlphaFold models are predicted in the absence of small molecules and ions/cofactors complicates their utilization for drug design. Previously, we optimized an HDAC11 AlphaFold model by adding the catalytic zinc ion and minimization in the presence of reported HDAC11 inhibitors. In the current study, we implement a comparative structure-based virtual screening approach utilizing the previously optimized HDAC11 AlphaFold model to identify novel and selective HDAC11 inhibitors. The stepwise virtual screening approach was successful in identifying a hit that was subsequently tested using an in vitro enzymatic assay. The hit compound showed an IC50 value of 3.5 µM for HDAC11 and could selectively inhibit HDAC11 over other HDAC subtypes at 10 µM concentration. In addition, we carried out molecular dynamics simulations to further confirm the binding hypothesis obtained by the docking study. These results reinforce the previously presented AlphaFold optimization approach and confirm the applicability of AlphaFold models in the search for novel inhibitors for drug discovery.

Deciphering Campylobacter jejuni DsbA1 protein dynamics in the presence of anti-virulent compounds: a multi-pronged computer-aided approach

The current study aims to evaluate Asinex library compounds against Campylobacter jejuni DsbA1 protein, a thiol disulfide oxidoreductase enzyme that plays a major role in the oxidative folding of bacterial virulence proteins, making it a promising anti-viral drug target. By employing several techniques of computer-aided drug design, BDC25697459, BDD33601083, and BDC30129064 were identified with binding energy scores of -8.8 kcal/mol, -8.8 kcal/mol, and -8.3 kcal/mol, respectively. However, the control molecule, tetraethylene glycol, exhibited a binding energy score of -7.0 kcal/mol. The control, BDD33601083, and BDC30129064 were unveiled to bind the same co-crystallized binding site (pocket 1), while BDC25697459 interacted with a new binding pocket (pocket 2) adjacent to the control binding region. The molecular dynamics simulation showed that complexes exhibit stable dynamics without significant global or residue-level fluctuations. The average RMSD values were in the range of 2.07 Å-2.45 Å. Similarly, mean RMSF was recorded between 1.30 and 1.42 Å. The C. jejuni DsbA1 was also observed as compact in the presence of the compounds, showing a mean RoG value in the range of 16.42 Å-16.55 Å. In terms of MM/PBSA binding energy, the BDC30129064 complex was ranked top with -44.88 ± 4.14 kcal/mol, whereas the positive control molecule exhibited -22.22 ± 3.33 kcal/mol. From a pharmacokinetic perspective, the compounds are suitable candidates for clinical trial investigation. Preliminary computational analysis of these virtual hits indicates that these compounds have a low potential for ADME and toxicity-associated liabilities. In summary, the compounds displayed a high affinity for the C. jejuni DsbA1 protein, indicating potential efficacy that requires further investigation.Communicated by Ramaswamy H. Sarma.

In Silico Identification of a Potential TNF-Alpha Binder Using a Structural Similarity: A Potential Drug Repurposing Approach to the Management of Alzheimer's Disease

Introduction

Alzheimer's disease (AD) is a neurodegenerative disorder with no conclusive remedy. Yohimbine, found in Rauwolfia vomitoria, may reduce brain inflammation by targeting tumour necrosis factor-alpha (TNFα), implicated in AD pathogenesis. Metoserpate, a synthetic compound, may inhibit TNFα. The study is aimed at assessing the potential utility of repurposing metoserpate for TNFα inhibition to reduce neuronal damage and inflammation in AD. The development of safe and effective treatments for AD is crucial to address the growing burden of the disease, which is projected to double over the next two decades.

Methods

Our study repurposed an FDA-approved drug as TNFα inhibitor for AD management using structural similarity studies, molecular docking, and molecular dynamics simulations. Yohimbine was used as a reference compound. Molecular docking used SeeSAR, and molecular dynamics simulation used GROMACS.

Results

Metoserpate was selected from 10 compounds similar to yohimbine based on pharmacokinetic properties and FDA approval status. Molecular docking and simulation studies showed a stable interaction between metoserpate and TNFα over 100 ns (100000 ps). This suggests a reliable and robust interaction between the protein and ligand, supporting the potential utility of repurposing metoserpate for TNFα inhibition in AD treatment.

Conclusion

Our study has identified metoserpate, a previously FDA-approved antihypertensive agent, as a promising candidate for inhibiting TNFα in the management of AD.

Different binding modalities of quercetin to inositol-requiring enzyme 1 of S. cerevisiae and human lead to opposite regulation

The flavonoid Quercetin (Qe) was identified as an activator of Inositol-requiring enzyme 1 (IRE1) in S. cerevisiae (scIre1p), but its impact on human IRE1 (hIRE1) remains controversial due to the absence of a conserved Qe binding site. We have explored the binding modes and effect of Qe on both scIre1p and hIRE1 dimers using in silico and in vitro approaches. The activation site in scIre1p stably accommodates both Qe and its derivative Quercitrin (Qi), thus enhancing the stability of the RNase pocket. However, the corresponding region in hIRE1 does not bind any of the two molecules. Instead, we show that both Qe and Qi block the RNase activity of hIRE1 in vitro, with sub-micromolar IC50 values. Our results provide a rationale for why Qe is an activator in scIre1p but a potent inhibitor in hIRE1. The identification of a new allosteric site in hIRE1 opens a promising window for drug development and UPR modulation.

Optimization of Novel α7 Nicotinic Acetylcholine Receptor Positive Allosteric Modulators and the Discovery of a Preclinical Development Candidate Molecule (RGH-560)

During optimization of a previously identified lead compound, attempts were made to optimize the reactive indole structural element, the suboptimal metabolic stability, as well as the low kinetic solubility. It was concluded that the indole was important for in vitro activity. With the aim of further improvements, more thorough modifications were also carried out. As a result, a new chemotype (the azetidinespirochromone family) was identified, which proved to be 1 order of magnitude less lipophilic retaining the same high level of in vitro potency as the lead series itself, however, with improved metabolic stability and kinetic solubility. Compound 53 showed the most balanced physicochemical and pharmacological profile with significant in vivo efficacy in the scopolamine-induced amnesia test. Based on these promising results, cognitive enhancement through the positive modulation of α7 nAChRs appears to be a viable approach. Compound 53 was selected to be a preclinical development candidate (as RGH-560).

Harnessing the druggability at orthosteric and allosteric sites of PD-1 for small molecule discovery by an integrated in silico pipeline

The PD-1/PD-L1 interaction is a promising target for small molecule inhibitors in cancer immunotherapy, but targeting this interface has been challenging. While efforts have been made to identify compounds that target the orthosteric sites, no reports have explored the potential of small molecules to target the allosteric region of PD-1. Therefore, our study aims to establish a pipeline to identify small molecules that can effectively bind to either the orthosteric or allosteric pockets of PD-1. We categorized the PD-1 interface into two hot-spot zones (P-and N-zones) based on extensive analysis of its structural, dynamical, and energetic properties. These zones correspond to the orthosteric and allosteric PPI sites, respectively, targeted by monoclonal antibodies. We used a guided virtual screening workflow to identify hits from ∼7 million compounds library, which were then clustered based on structural similarity and assessed by interaction fingerprinting. The selective and diverse chemical representatives were subjected to MD simulations and binding energetics calculations to filter out false positives and identify actual binders. Binding poses metadynamics calculations confirmed the stability of the final hits in the pocket. This study emphasizes the need for an integrated pipeline that uses molecular dynamics simulations and binding energetics to identify potential binders for the dynamic PD-1/PD-L1 interface, due to the lack of small molecule co-crystals. Only a few potential binders were discovered from a large pool of molecules targeting both the allosteric and orthosteric zones. Our results suggest that the allosteric site has more potential than the orthosteric site for inhibitor design. The identified "computational hits" hold potential as starting points for in vitro evaluations followed by hit-to-lead optimization. Overall, this study represents an effort to establish a computational pipeline for exploring and enriching both the allosteric and orthosteric sites of PPI interfaces, "a tough but indispensable nut to crack".

Cheminformatics identification of modulators of key carbohydrate-metabolizing enzymes from C. cujete for type-2 diabetes mellitus intervention

Purpose

The therapeutic use of oral hypoglycaemic agents in the management of type-2 diabetes mellitus (T2DM) is without adverse effects; thus, calls for alternative and novel candidates from natural products in medicinal plants.

Method

The study explored molecular docking and molecular dynamics (MD) simulation approaches to identify key antidiabetic metabolites from Crescentia cujete.

Results

Molecular docking results identified four and/or five best compounds against each target enzyme (alpha-glucosidase, dipeptidyl peptidase-IV, aldose reductase, and protein tyrosine phosphatase-1B (PTP-1B)) implicated in diabetes. The resulting complexes (except against PTP-1B) had higher docking scores above respective standards (acarbose, Diprotin A, ranirestat). The MD simulation results revealed compounds such as benzoic acid (-48.414 kcal/mol) and phytol (-45.112 kcal/mol) as well as chlorogenic acid (-42.978 kcal/mol) and naringenin (-31.292 kcal/mol) had higher binding affinities than the standards [acarbose (-28.248 kcal/mol), ranirestat (-21.042 kcal/mol)] against alpha-glucosidase and aldose reductase, respectively while Diprotin A (-45.112 kcal/mol) and ursolic acid (-18.740 kcal/mol) presented superior binding affinities than the compounds [luteolin (-41.957 kcal/mol and naringenin (-16.518 kcal/mol)] against DPP-IV and PTP-1B respectively.

Conclusion

While isoflavone (alpha-glucosidase), xylocaine (DPP-IV), luteolin (aldose reductase,) and chlorogenic acid (PTP-1B) were affirmed as the best inhibitors of respective enzyme targets, luteolin, and chlorogenic acid may be suggested and proposed as probable candidates against T2DM and related retinopathy complication based on their structural stability, compactness and affinity for three (DPP-IV, aldose reductase, and PTP-1B) of the four targets investigated. Further studies are warranted in vitro and in vivo on the antihyperglycaemic effects of these drug candidates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40200-023-01249-7.

Assessing Metadynamics and Docking for Absolute Binding Free Energy Calculations Using Severe Acute Respiratory Syndrome Coronavirus 2 Main Protease Inhibitors

Absolute binding free energy (ABFE) calculations can be an important part of the drug discovery process by identifying molecules that have the potential to be strong binders for a biomolecular target. Recent work has used free energy perturbation (FEP) theory for these calculations, focusing on a set of 16 inhibitors of the severe acute respiratory syndrome coronavirus 2 main protease (Mpro). Herein, the same data set is evaluated by metadynamics (MetaD), four different docking programs, and molecular mechanics with generalized Born and surface area solvation. MetaD yields a Kendall τ distance of 0.28 and Pearson r2 of 0.49, which reflect somewhat less accuracy than that from the ABFE FEP results. Notably, it is demonstrated that an ensemble docking protocol by which each ligand is docked into the 13 crystal structures in this data set provides improved performance, particularly when docking is carried out with Glide XP (Kendall τ distance = 0.20, Pearson r2 = 0.71), Glide SP (Kendall τ distance = 0.19, Pearson r2 = 0.66), or AutoDock 4 (Kendall τ distance = 0.21, Pearson r2 = 0.55). The best results are obtained with "superconsensus" docking by averaging the 52 results for each compound using the 4 docking protocols and all 13 crystal structures (Kendall τ distance = 0.18, Pearson r2 = 0.73).

Computational Investigations on Reaction Mechanisms of the Covalent Inhibitors Ponatinib and Analogs Targeting the Extracellular Signal-Regulated Kinases

As an important cancer therapeutic target, extracellular signal-regulated kinases (ERK) are involved in triggering various cellular responses in tumors. Regulation of the ERK signaling pathway by the small molecular inhibitors is highly desired for the sake of cancer therapy. In contrast to the routine inhibitors targeting ERKs through long-range non-bonding interactions, Ponatinib, a covalent inhibitor to ERK2 with a macrocyclic structure characterized by the α,β-C=C unsaturated ketone, can form the stable -C(S)-C(H)-type complex via the four-center barrier due to the nucleophilic addition reaction of the thiol group of the Cys166 residue of ERK2 with the C=C double bond of Ponatinib with reaction free-energy barrier of 47.2 kcal/mol. Reaction mechanisms for the covalent binding were calculated using QM/MM methods and molecular dynamics simulations. The interaction modes and the corresponding binding free energies were obtained for the non-covalent and covalent complexation. The binding free energies of the non-covalent and covalent inhibitions are 14.8 kcal/mol and 33.4 kcal/mol, respectively. The mechanistic study stimulated a rational design on the modified Ponatinib structure by substituting the C=C bond with the C=N bond. It was demonstrated that the new compound exhibits better inhibition activity toward ERK2 in term of both thermodynamic and kinetic aspects through the covalent binding with a lower reaction free-energy barrier of 23.1 kcal/mol. The present theoretical work sheds new light on the development of the covalent inhibitors for the regulation of ERKs.

A systematic pipeline of protein structure selection for computer-aided drug discovery: A case study on T790M/L858R mutant EGFR structures

Virtual screening (VS) is a routine method to evaluate chemical libraries for lead identification. Therefore, the selection of appropriate protein structures for VS is an essential prerequisite to identify true actives during docking. But the presence of several crystal structures of the same protein makes it difficult to select one or few structures rationally for screening. Therefore, a computational prioritization protocol has been developed for shortlisting crystal structures that identify true active molecules with better efficiency. As identification of small-molecule inhibitors is an important clinical requirement for the T790M/L858R (TMLR) EGFR mutant, it has been selected as a case study. The approach involves cross-docking of 21 co-crystal ligands with all the structures of the same protein to select structures that dock non-native ligands with lower RMSD. The cross docking performance was then correlated with ligand similarity and binding-site conformational similarity. Eventually, structures were shortlisted by integrating cross-docking performance, and ligand and binding-site similarity. Thereafter, binding pose metadynamics was employed to identify structures having stable co-crystal ligands in their respective binding pockets. Finally, different enrichment metrics like BEDROC, RIE, AUAC, and EF1% were evaluated leading to the identification of five TMLR structures (5HCX, 5CAN, 5CAP, 5CAS, and 5CAO). These structures docked a number of non-native ligands with low RMSD, contain structurally dissimilar ligands, have conformationally dissimilar binding sites, harbor stable co-crystal ligands, and also identify true actives early. The present approach can be implemented for shortlisting protein targets of any other important therapeutic kinases.

Preclinical Characterization of 22-(4'-Pyridinecarbonyl) Jorunnamycin A against Lung Cancer Cell Invasion and Angiogenesis via AKT/mTOR Signaling

Non-small-cell lung cancer (NSCLC), the most prevalent form of lung cancer, is associated with an unfavorable prognosis owing to its high rate of metastasis. Thus, the identification of new drugs with potent anticancer activities is essential to improve the clinical outcome of this disease. Marine organisms exhibit a diverse source of biologically active compounds with anticancer effects. The anticancer effects of jorunnamycin A (JA) derived from the Thai blue sponge (Xestospongia sp.) and 22-(4'-pyridinecarbonyl) jorunnamycin A (22-(4'-py)-JA), the semisynthetic derivative of JA, have been reported. The present study aimed to investigate the impact of 22-(4'-py)-JA on NSCLC metastasis using in vitro, in vivo, and in silico approaches. The JA derivative inhibited tumor cell invasion and tube formation in human umbilical vein endothelial cells (HUVECs). The computational analysis demonstrated strong and stable interactions between 22-(4'-py)-JA and the AKT protein. Further examinations into the molecular mechanisms revealed the suppression of AKT/mTOR/p70S6K signaling by 22-(4'-py)-JA, leading to the downregulation of matrix metalloproteinases (MMP-2 and MMP-9), hypoxia-inducible factor-1α (HIF-1α), and vascular endothelial growth factor (VEGF). Furthermore, 22-(4'-py)-JA suppressed in vivo metastasis by decreasing the number of colonies in the lung. These findings indicated the antimetastasis activity of 22-(4'-py)-JA, which might prove useful for further clinical applications.

Identification and characterization of 7-azaindole derivatives as inhibitors of the SARS-CoV-2 spike-hACE2 protein interaction

The COVID-19 pandemic, caused by SARS-CoV-2, has become a global public health crisis. The entry of SARS-CoV-2 into host cells is facilitated by the binding of its spike protein (S1-RBD) to the host receptor hACE2. Small molecule compounds targeting S1-RBD-hACE2 interaction could provide an alternative therapeutic strategy sensitive to viral mutations. In this study, we identified G7a as a hit compound that targets the S1-RBD-hACE2 interaction, using high-throughput screening in the SARS2-S pseudovirus model. To enhance the antiviral activity of G7a, we designed and synthesized a series of novel 7-azaindole derivatives that bind to the S1-RBD-hACE2 interface. Surprisingly, ASM-7 showed excellent antiviral activity and low cytotoxicity, as confirmed by pseudovirus and native virus assays. Molecular docking and molecular dynamics simulations revealed that ASM-7 could stably bind to the binding interface of S1-RBD-hACE2, forming strong non-covalent interactions with key residues. Furthermore, the binding of ASM-7 caused alterations in the structural dynamics of both S1-RBD and hACE2, resulting in a decrease in their binding affinity and ultimately impeding the viral invasion of host cells. Our findings demonstrate that ASM-7 is a promising lead compound for developing novel therapeutics against SARS-CoV-2.

Computational Design, Synthesis, and Biological Evaluation of Diimidazole Analogues Endowed with Dual PCSK9/HMG-CoAR-Inhibiting Activity

Proprotein convertase subtilisin/kexin 9 (PCSK9) is responsible for the degradation of the hepatic low-density lipoprotein receptor (LDLR), which regulates circulating cholesterol levels. Consequently, the PCSK9 inhibition is a valuable therapeutic approach for the treatment of hypercholesterolemia and cardiovascular diseases. In our studies, we discovered Rim13, a polyimidazole derivative reducing the protein-protein interaction between PCSK9 and LDLR with an IC₅₀ of 1.6 μM. The computational design led to the optimization of the shape of the PCSK9/ligand complementarity, enabling the discovery of potent diimidazole derivatives. In fact, carrying out biological assays to fully characterize the cholesterol-lowering activity of the new analogues and using both biochemical and cellular techniques, compound Dim16 displayed improved PCSK9 inhibitory activity (IC₅₀ 0.9 nM). Interestingly, similar to other lupin-derived peptides and their synthetic analogues, some compounds in this series showed dual hypocholesterolemic activity since some of them complementarily inhibited the 3-hydroxy-3-methylglutaryl coenzyme A reductase.

Study of SQ109 analogs binding to mycobacterium MmpL3 transporter using MD simulations and alchemical relative binding free energy calculations

N-geranyl-N΄-(2-adamantyl)ethane-1,2-diamine (SQ109) is a tuberculosis drug that has high potency against Mycobacterium tuberculosis (Mtb) and may function by blocking cell wall biosynthesis. After the crystal structure of MmpL3 from Mycobacterium smegmatis in complex with SQ109 became available, it was suggested that SQ109 inhibits Mmpl3 mycolic acid transporter. Here, we showed using molecular dynamics (MD) simulations that the binding profile of nine SQ109 analogs with inhibitory potency against Mtb and alkyl or aryl adducts at C-2 or C-1 adamantyl carbon to MmpL3 was consistent with the X-ray structure of MmpL3 - SQ109 complex. We showed that rotation of SQ109 around carbon-carbon bond in the monoprotonated ethylenediamine unit favors two gauche conformations as minima in water and lipophilic solvent using DFT calculations as well as inside the transporter's binding area using MD simulations. The binding assays in micelles suggested that the binding affinity of the SQ109 analogs was increased for the larger, more hydrophobic adducts, which was consistent with our results from MD simulations of the SQ109 analogues suggesting that sizeable C-2 adamantyl adducts of SQ109 can fill a lipophilic region between Y257, Y646, F260 and F649 in MmpL3. This was confirmed quantitatively by our calculations of the relative binding free energies using the thermodynamic integration coupled with MD simulations method with a mean assigned error of 0.74 kcal mol^-1 compared to the experimental values.

Computational characteristics of the structure-activity relationship of inhibitors targeting Pks13-TE domain

Multiple studies have established the Pks13-TE domain as a promising target for anti-tuberculosis drug development. However, recent findings have revealed that the lead compound currently in the pipeline for Pks13-TE has significant cardiotoxicity issues. Given the pressing need for new chemical structures for Pks13-TE inhibitors, this study aims to provide a detailed understanding of the Pks13-TE domain binding site through the application of computational chemical biology techniques. Our results highlight the size and shape of the Pks13-TE domain binding pocket, key residues including Asp1644, Asn1640, Phe1670, and Tyr1674 within the pocket, and inhibitor pharmacophore characteristics such as aromatic ring sites, positively charged sites, and hydrogen bond donors. To our knowledge, these simulation results are novel and contribute to the discovery of next-generation Pks13-TE inhibitors without similar prior studies.

AI-driven drug repurposing and binding pose meta dynamics identifies novel targets for Monkeypox virus

Monkeypox virus (MPXV) was confirmed in May 2022 and designated a global health emergency by WHO in July 2022. MPX virions are big, enclosed, brick-shaped, and contain a linear, double-stranded DNA genome as well as enzymes. MPXV particles bind to the host cell membrane via a variety of viral-host protein interactions. As a result, the wrapped structure is a potential therapeutic target. DeepRepurpose, an artificial intelligence-based compound-viral proteins interaction framework, was used via a transfer learning setting to prioritize a set of FDA approved and investigational drugs which can potentially inhibit MPXV viral proteins. To filter and narrow down the lead compounds from curated collections of pharmaceutical compounds, we used a rigorous computational framework that included homology modeling, molecular docking, dynamic simulations, binding free energy calculations, and binding pose metadynamics. We identified Elvitegravir as a potential inhibitor of MPXV virus using our comprehensive pipeline.

Identification of potential ATP-competitive cyclin-dependent kinase 1 inhibitors: De novo drug generation, molecular docking, and molecular dynamics simulation

Cyclin-dependent kinases 1 (CDK1) has been identified as a potential target for the search for new antitumor drugs. However, no clinically effective CDK1 inhibitors are now available for cancer treatment. Therefore, this study aimed to offer potential CDK1 inhibitors using de novo drug generation, molecular docking, and molecular dynamics (MD) simulation studies. We first utilized the BREED algorithm (a de novo drug generation approach) to produce a novel library of small molecules targeting CDK1. To initially obtain novel potential CDK1 inhibitors with favorable physicochemical properties and excellent druggability, we performed a virtual rule-based rational drug screening on our generated library and found ten initial hits. Then, the molecular interactions and dynamic stability of these ten initial hits and CDK1 complexes during their all-atom MD simulations (total 18 μs) and binding pose metadynamics simulations were investigated, resulting in five final hits. Furthermore, another MD simulation (total 2.1 μs) with different force fields demonstrated the binding ability of the five hits to CDK1. It was found that these five hits, CBMA001 (ΔG = -29.88 kcal/mol), CBMA002 (ΔG = -34.89 kcal/mol), CBMA004 (ΔG = -32.47 kcal/mol), CBMA007 (ΔG = -31.16 kcal/mol), and CBMA008 (ΔG = -34.78 kcal/mol) possessed much greater binding affinity to CDK1 than positive compound Flavopiridol (FLP, ΔG = -25.38 kcal/mol). Finally, CBMA002 and CBMA004 were identified as excellent selective CDK1 inhibitors in silico. Together, this study provides a workflow for rational drug design and two promising selective CDK1 inhibitors that deserve further investigation.

Application of Molecular Simulation Methods in Food Science: Status and Prospects

Molecular simulation methods, such as molecular docking, molecular dynamic (MD) simulation, and quantum chemical (QC) calculation, have become popular as characterization and/or virtual screening tools because they can visually display interaction details that in vitro experiments can not capture and quickly screen bioactive compounds from large databases with millions of molecules. Currently, interdisciplinary research has expanded molecular simulation technology from computer aided drug design (CADD) to food science. More food scientists are supporting their hypotheses/results with this technology. To understand better the use of molecular simulation methods, it is necessary to systematically summarize the latest applications and usage trends of molecular simulation methods in the research field of food science. However, this type of review article is rare. To bridge this gap, we have comprehensively summarized the principle, combination usage, and application of molecular simulation methods in food science. We also analyzed the limitations and future trends and offered valuable strategies with the latest technologies to help food scientists use molecular simulation methods.

Computational Approaches to the Rational Design of Tubulin-Targeting Agents

Microtubules are highly dynamic polymers of α,β-tubulin dimers which play an essential role in numerous cellular processes such as cell proliferation and intracellular transport, making them an attractive target for cancer and neurodegeneration research. To date, a large number of known tubulin binders were derived from natural products, while only one was developed by rational structure-based drug design. Several of these tubulin binders show promising in vitro profiles while presenting unacceptable off-target effects when tested in patients. Therefore, there is a continuing demand for the discovery of safer and more efficient tubulin-targeting agents. Since tubulin structural data is readily available, the employment of computer-aided design techniques can be a key element to focus on the relevant chemical space and guide the design process. Due to the high diversity and quantity of structural data available, we compiled here a guide to the accessible tubulin-ligand structures. Furthermore, we review different ligand and structure-based methods recently used for the successful selection and design of new tubulin-targeting agents.

Structural Elucidation of Relevant Gibberellic Acid Impurities and In Silico Investigation of Their Interaction with Soluble Gibberellin Receptor GID1

Gibberellin derivatives are a family of tetracyclic diterpenoid plant hormones used in agriculture as plant growth regulators included in the European Directive 91/414. In the pesticide peer review process and to assess their toxicological relevance and product chemical equivalence, the European Food Safety Authority (EFSA) highlighted data gaps such as the identification of hydrolysis products and unknown impurities. The aspect of impurity characterization and quantitation is challenging and requires the use of hyphenated analytical techniques. In this regard, we used an LC-QTOF/MS and NMR analysis for the characterization of gibberellic acid impurities found in technical products. Gibberellic acid impurities such as gibberellin A1 (GA ₁ ), 3-isolactone gibberellic acid (iso-GA ₃ ), gibberellenic acid, 1α,2α-epoxygibberellin A3 (2-epoxy- GA ₃ ), and (1α,2β,3α,4bβ,10β)-2,3,7-trihydroxy-1-methyl-8-methylenegibb-4-ene-1,10-dicarboxylic acid were identified and successfully characterized. Moreover, an in silico investigation on selected gibberellic acid impurities and derivatives and their interactions with a gibberellin insensitive dwarf1 (GID1) receptor has been carried out by means of induced fit docking (IFD), generalized-Born surface area (MM-GBSA), and metadynamics (MTD) experiments. A direct HPLC method with DAD and MS for the detection of gibberellic acid and its impurities in a technical sample has been developed. Moreover, by means of the in silico characterization of the GID1 receptor-binding pocket, we investigated the receptor affinity of the selected gibberellins, identifying compounds (2) and (4) as the most promising hit to lead compounds.

In silico and in vitro studies confirm Ondansetron as a novel acetylcholinesterase and butyrylcholinesterase inhibitor

Alzheimer's disease (AD) is a progressive neurodegenerative disorder that is growing rapidly among the elderly population around the world. Studies show that a lack of acetylcholine and butyrylcholine due to the overexpression of enzymes Acetylcholinesterase (AChE) and Butyrylcholinesterase (BChE) may lead to reduced communication between neuron cells. As a result, seeking novel inhibitors targeting these enzymes might be vital for the future treatment of AD. Ondansetron is used to prevent nausea and vomiting caused by chemotherapy or radiation treatments and is herein shown to be a potent inhibitor of cholinesterase. Comparison is made between Ondansetron and FDA-approved cholinesterase inhibitors Rivastigmine and Tacrine. Molecular docking demonstrates that interactions between the studied ligand and aromatic residues in the peripheral region of the active site are important in binding. Molecular dynamics simulations and binding pose metadynamics show that Ondansetron is highly potent against both enzymes and far better than Rivastigmine. Inhibitor activities evaluated by in vitro studies confirm that the drug inhibits AChE and BChE by non-competitive and mixed inhibition, respectively, with IC50 values 33 µM (AChE) and 2.5 µM (BChE). Based on the findings, we propose that Ondansetron may have therapeutic applications in inhibiting cholinesterase, especially for BChE.

Drug repurposing combined with MM/PBSA based validation strategies towards MEK inhibitors screening

Emergence of oncogenic mutations in the MAPK pathway gaining more impact in the recent years. Importantly, MEK is a core element of this pathway as it is easy to inhibit and is a gatekeeper of multiple malignancies. Therefore, we performed in-silico strategy to screen repurposed candidate for MEK protein using a library of 11,808 compounds from different clusters in the DrugBank database. Glide docking, Prime-MM/GBSA and QikProp analysis were implemented to retrieve the hits with high precision. The stability of the binding mode and binding affinity of the resultant hit were explored using molecular dynamic simulations and MM/PBSA approach. The results highlight that Nebivolol (DB04861) not only achieved a stable conformation in the MEK binding pocket but also displayed highest binding affinity than the other molecules investigated in our study. Taken together, we hypothesized that Nebivolol is an excellent candidate for the inhibition of MEK in NSCLC patients in future.Communicated by Ramaswamy H. Sarma.

Open Binding Pose Metadynamics: An Effective Approach for the Ranking of Protein-Ligand Binding Poses

Predicting the correct pose of a ligand binding to a protein and its associated binding affinity is of great importance in computer-aided drug discovery. A number of approaches have been developed to these ends, ranging from the widely used fast molecular docking to the computationally expensive enhanced sampling molecular simulations. In this context, methods such as coarse-grained metadynamics and binding pose metadynamics (BPMD) use simulations with metadynamics biasing to probe the binding affinity without trying to fully converge the binding free energy landscape in order to decrease the computational cost. In BPMD, the metadynamics bias perturbs the ligand away from the initial pose. The resistance of the ligand to this bias is used to calculate a stability score. The method has been shown to be useful in reranking predicted binding poses from docking. Here, we present OpenBPMD, an open-source Python reimplementation and reinterpretation of BPMD. OpenBPMD is powered by the OpenMM simulation engine and uses a revised scoring function. The algorithm was validated by testing it on a wide range of targets and showing that it matches or exceeds the performance of the original BPMD. We also investigated the role of accurate water positioning on the performance of the algorithm and showed how the combination with a grand-canonical Monte Carlo algorithm improves the accuracy of the predictions.

Exploring new catechin derivatives as SARS-CoV-2 Mpro inhibitors from tea by molecular networking, surface plasma resonance, enzyme inhibition, induced fit docking, and metadynamics simulations

SARS-CoV-2 M^pro (Mpro) is the critical cysteine protease in coronavirus viral replication. Tea polyphenols are effective M^pro inhibitors. Therefore, we aim to isolate and synthesize more novel tea polyphenols from Zhenghedabai (ZHDB) white tea methanol-water (MW) extracts that might inhibit COVID-19. Through molecular networking, 33 compounds were identified and divided into 5 clusters. Further, natural products molecular network (MN) analysis showed that MN1 has new phenylpropanoid-substituted ester-catechin (PSEC), and MN5 has the important basic compound type hydroxycinnamoylcatechins (HCCs). Thus, a new PSEC (1, PSEC636) was isolated, which can be further detected in 14 green tea samples. A series of HCCs were synthesized (2-6), including three new acetylated HCCs (3-5). Then we used surface plasmon resonance (SPR) to analyze the equilibrium dissociation constants (K_D) for the interaction of 12 catechins and M^pro. The K_D values of PSEC636 (1), EGC-C (2), and EC-CDA (3) were 2.25, 2.81, and 2.44 μM, respectively. Moreover, compounds 1, 2, and 3 showed the potential M^pro inhibition with IC₅₀ 5.95 ± 0.17, 9.09 ± 0.22, and 23.10 ± 0.69 μM, respectively. Further, we used induced fit docking (IFD), binding pose metadynamics (BPMD), and molecular dynamics (MD) to explore the stable binding pose of M^pro-1, showing that 1 could tightly bond with the amino acid residues THR²⁶, HIS⁴¹, CYS⁴⁴, TYR⁵⁴, GLU¹⁶⁶, and ASP¹⁸⁷. The computer modeling studies reveal that the ester, acetyl, and pyrogallol groups could improve inhibitory activity. Our research suggests that these catechins are effective M^pro inhibitors, and might be developed as therapeutics against COVID-19.

Structure-Based Discovery of a Novel Class of Small-Molecule Pure Antagonists of Integrin αVβ3

Inhibitors of integrin αVβ3 have therapeutic promise for a variety of diseases. Most αVβ3-targeting small molecules patterned after the RGD motif are partial agonists because they induce a high-affinity, ligand-binding conformation and prime the receptor to bind the ligand without an activating stimulus, in part via a charge-charge interaction between their aspartic acid carboxyl group and the metal ion in the metal-ion-dependent adhesion site (MIDAS). Building upon our previous studies on the related integrin αIIbβ3, we searched for pure αVβ3 antagonists that lack this typical aspartic acid carboxyl group and instead engage through direct binding to one of the coordinating residues of the MIDAS metal ion, specifically β3 E220. By in silico screening of two large chemical libraries for compounds interacting with β3 E220, we indeed discovered a novel molecule that does not contain an acidic carboxyl group and does not induce the high-affinity, ligand-binding state of the receptor. Functional and structural characterization of a chemically optimized version of this compound led to the discovery of a novel small-molecule pure αVβ3 antagonist that (i) does not prime the receptor to bind the ligand and does not induce hybrid domain swing-out or receptor extension as judged by antibody binding and negative-stain electron microscopy, (ii) binds at the RGD-binding site as predicted by metadynamics rescoring of induced-fit docking poses and confirmed by a cryo-electron microscopy structure of the compound-bound integrin, and (iii) coordinates the MIDAS metal ion via a quinoline moiety instead of an acidic carboxyl group.

An in silico comparative study of curcumin and 2-deoxyuridine nucleoside derivatives: Reveals the role of angiogenin in ER stress-induced apoptosis signaling

Angiogenin (ANG) protein plays a crucial role in angiogenesis, neovascularization, and cancer metastasis in NSCLC (non-small cell lung cancer) via non-coding tiRNA. It protects the cell under ER (endoplasmic reticulum) stress-induced apoptosis through the translational reprogramming process. Although B82 (Curcumin derivatives) induces ER stress-induced apoptosis, its mechanism of action was not studied. Therefore, it was hypothesized that the ribonucleolytic activity of ANG may be regulated by B82, resulting in modulated ER stress signaling for apoptosis. Hence, we designed and proposed a synthesis scheme for RNA-based anti-angiogenic derivatives of 2-deoxyuridine nucleoside forming peptide bond with amino acids like serine (Ser-3) and para-hydroxy-phenyl glycine (Normtyr-1) and compared B82 with them to know the binding affinity with ANG, anti-angiogenic potential, and its probable mechanism of anti-RNase activity through MD simulation study. Therefore, using Gromos96 43a1 and 43a2 force fields, MD simulation was performed to investigate binding affinity, ligand-induced molecular surface area change, conformational change, and dynamics of catalytic site residues to predict ligand binding to ANG in this study. The obtained binding free energy (∆G_bind ) result showed the total average ∆G_bind as -113.480 ± 1.682 (Normtyr-1) > -53.038 ± 33.069 (B82) > -27.909 ± 16.438 (Ser-3) kJ/mole specify role of B82 in regulating ER stress signaling induced apoptosis through ANG ribonucleolytic activity inhibition, suitability of 43a2 force fields and methodology in ligand screening. It shows the crucial role of Leu115 and His13 residue involvement in total ∆G_bind contribution. Hence, based on the MD result, novel conformation of catalytic residues, and ∆G_bind , a promising combination candidate could be proposed for metastatic NSCLC therapy.

Computational molecular interaction between SARS-CoV-2 main protease and theaflavin digallate using free energy perturbation and molecular dynamics

Our objective was to identify the molecule which can inhibit SARS-CoV-2 main protease and can be easily procured. Natural products may provide such molecules and can supplement the current custom chemical synthesis-based drug discovery for this objective. A combination of docking approaches, scoring functions, classical molecular dynamic simulation, binding pose metadynamics, and free energy perturbation calculations have been employed in this study. Theaflavin digallate has been observed in top-scoring compounds after the three independent virtual screening simulations of 598435 compounds (unique 27256 chemical entities). The main protease-theaflavin digallate complex interacts with critical active site residues of the main protease in molecular dynamics simulation independent of the explored computational framework, simulation time, initial structure, and force field used. Theaflavin digallate forms approximately three hydrogen bonds with Glutamate166 of main protease, primarily through hydroxyl groups in the benzene ring of benzo(7)annulen-6-one, along with other critical residues. Glu166 is the most critical amino acid for main protease dimerization, which is necessary for catalytic activity. The estimated binding free energy, calculated by Amber and Schrodinger MMGBSA module, reflects a high binding free energy between theaflavin digallate and main protease. Binding pose metadynamics simulation shows the highly persistent H-bond and a stable pose for the theaflavin digallate-main protease complex. Using method control, experimental controls, and test set, alchemical transformation studies confirm high relative binding free energy of theaflavin digallate with the main protease. Computational molecular interaction suggests that theaflavin digallate can inhibit the main protease of SARS-CoV-2.

Pharmacological characterization of novel synthetic opioids: Isotonitazene, metonitazene, and piperidylthiambutene as potent MU opioid receptor agonists

Recent trends of opioid abuse and related fatalities have highlighted the critical role of Novel Synthetic Opioids (NSOs). We studied the μ-opioid-like properties of isotonitazene (ITZ), metonitazene (MTZ), and piperidylthiambutene (PTB) using different approaches. In vitro studies showed that ITZ and MTZ displayed a higher potency in both rat membrane homogenates (EC₅₀: 0.99 and 19.1 nM, respectively) and CHO-MOR (EC₅₀: 0.71 and 10.0 nM, respectively) than [D-Ala², NMe-Phe⁴, Gly-ol⁵]-enkephalin (DAMGO), with no difference in maximal efficacy (Emax) between DAMGO and NSOs. ITZ also has higher affinity (Ki: 0.06 and 0.05 nM) at the MOR than DAMGO in both systems, whilst MTZ has higher affinity in CHO-MOR (Ki = 0.23 nM) and similar affinity in rat cerebral cortex (Ki = 0.22 nM). PTB showed lower affinity and potency than DAMGO. In vivo, ITZ displayed higher analgesic potency than fentanyl and morphine (ED₅₀: 0.00156, 0.00578, 2.35 mg/kg iv, respectively); ITZ (0.01 mg/kg iv) and MTZ (0.03 mg/kg iv) reduced behavioral activity and increased dialysate dopamine (DA) in the NAc shell (max. about 200% and 170% over basal value, respectively. Notably, ITZ elicited an increase in DA comparable to that of higher dose of morphine (1 mg/kg iv), but higher than the same dose of fentanyl (0.01 mg/kg iv). In silico, induced fit docking (IFD) and metadynamic simulations (MTD) showed that binding modes and structural changes at the receptor, ligand stability, and the overall energy score of NSOs were consistent with the results of the biological assays.

A combination strategy of structure-based virtual screening, MM-GBSA, cross docking, molecular dynamics and metadynamics simulations used to investigate natural compounds as potent and specific inhibitors of tumor linked human carbonic anhydrase IX

Cancer remains a serious health concern representing one of the leading causes of deaths worldwide. The enzyme human carbonic anhydrase IX (hCA IX) is found to be over-expressed in many cancer types and its selective inhibition over its cytosolic off-target isoform, human carbonic anhydrase II (hCA II), represents a potential area of research in the development of novel anticancer compounds. This work is concerned with the use of various in silico tools for the identification of natural product based molecules that can selectively inhibit hCA IX over hCA II. MM-GBSA assisted structure-based virtual screening against hCA IX was performed for nearly 225,000 natural products imported from the ZINC15 database. The obtained hits were checked for their potency by considering acetazolamide, the bound inhibitor of hCA IX, as the reference molecule, and 121 molecules were identified as potent hCA IX inhibitors. After ensuring their potency, cross-docking, followed by MM-GBSA calculations of the hits with hCA II, was performed, and their selectivity was assessed by considering the hCA IX selective compound SLC-0111 as the reference molecule, and 50 natural products were identified as potent as well as selective hCA IX inhibitors. Molecules with the quinoline scaffold showed the highest selectivity, and their selectivity was attributed to the strong electrostatic interactions of the zinc binding group (ZBG) with the active site Zn(II) ion. Furthermore, the stability of the binding modes of the top hCA IX selective hits was ensured by performing molecular dynamics (MD) simulations, which clearly proved that one of the short-listed molecules is truly selective, as it does not interact with the active site Zn(II) ion of hCA II, but interacts strongly with this ion in hCA IX. Bonding pose metadynamics studies revealed that the ligand moves to a more stable binding site from the one predicted by the docking studies and shows stronger interaction with the protein and Zn(II) at this binding site. The ligand is not likely to have issues with bioavailability. As a result, this ligand can be taken for bioassay testing and subsequently used as a feasible therapeutic treatment for a variety of cancer types. Communicated by Ramaswamy H. Sarma.

The Current Perspectives in Clinical Research: Computer-Assisted Drug Designing, Ethics, and Good Clinical Practice

In the era of emerging microbial and non-communicable diseases and re-emerging microbial infections, the medical fraternity and the public are plagued by under-preparedness. It is evident by the severity of the Coronavirus disease (COVID-19) pandemic that novel microbial diseases are a challenge and are challenging to control. This is mainly attributed to the lack of complete knowledge of the novel microbe’s biology and pathogenesis and the unavailability of therapeutic drugs and vaccines to treat and control the disease. Clinical research is the only answer utilizing which can handle most of these circumstances. In this review, we highlight the importance of computer-assisted drug designing (CADD) and the aspects of molecular docking, molecular superimposition, 3D-pharmacophore technology, ethics, and good clinical practice (GCP) for the development of therapeutic drugs, devices, and vaccines.

Identification of Potential WSB1 Inhibitors by AlphaFold Modeling, Virtual Screening, and Molecular Dynamics Simulation Studies

WD40 repeat and SOCS box containing 1 (WSB1) consists of seven WD40 repeat structural domains at the N-terminal end and one SOCS box structural domain at the C-terminal end. WSB1 promotes cancer progression by affecting the Von Hippel–Lindau tumor suppressor protein (pVHL) and upregulating hypoxia inducible factor-1α (HIF-1α) target gene expression. However, the crystal structure of WSB1 has not been reported, which is not beneficial to the research on WSB1 inhibitors. Therefore, we focused on specific small molecule inhibitors of WSB1. This study applied virtual screening and molecular dynamics simulations; finally, 20 compounds were obtained. Among them, compound G490-0341 showed the best stable structure and was a promising composite for further development of WSB1 inhibitors.

Computational-Based Discovery of the Anti-Cancer Activities of Pyrrole-Based Compounds Targeting the Colchicine-Binding Site of Tubulin

Despite the tubulin-binding agents (TBAs) that are widely used in the clinic for cancer therapy, tumor resistance to TBAs (both inherited and acquired) significantly impairs their effectiveness, thereby decreasing overall survival (OS) and progression-free survival (PFS) rates, especially for the patients with metastatic, recurrent, and unresectable forms of the disease. Therefore, the development of novel effective drugs interfering with the microtubules' dynamic state remains a big challenge in current oncology. We report here about the novel ethyl 2-amino-1-(furan-2-carboxamido)-5-(2-aryl/tert-butyl-2-oxoethylidene)-4-oxo-4,5-dihydro-1H-pyrrole-3-carboxylates (EAPCs) exhibiting potent anti-cancer activities against the breast and lung cancer cell lines in vitro. This was due to their ability to inhibit tubulin polymerization and induce cell cycle arrest in M-phase. As an outcome, the EAPC-treated cancer cells exhibited a significant increase in apoptosis, which was evidenced by the expression of cleaved forms of PARP, caspase-3, and increased numbers of Annexin-V-positive cells. By using the in silico molecular modeling methods (e.g., induced-fit docking, binding metadynamics, and unbiased molecular dynamics), we found that EAPC-67 and -70 preferentially bind to the colchicine-binding site of tubulin. Lastly, we have shown that the EAPCs indicated above and colchicine utilizes a similar molecular mechanism to inhibit tubulin polymerization via targeting the T7 loop in the β-chain of tubulin, thereby preventing the conformational changes in the tubulin dimers required for their polymerization. Collectively, we identified the novel and potent TBAs that bind to the colchicine-binding site and disrupt the microtubule network. As a result of these events, the compounds induced a robust cell cycle arrest in M-phase and exhibited potent pro-apoptotic activities against the epithelial cancer cell lines in vitro.

Computational-Based Discovery of the Anti-Cancer Activities of Pyrrole-Based Compounds Targeting the Colchicine-Binding Site of Tubulin

Despite the tubulin-binding agents (TBAs) that are widely used in the clinic for cancer therapy, tumor resistance to TBAs (both inherited and acquired) significantly impairs their effectiveness, thereby decreasing overall survival (OS) and progression-free survival (PFS) rates, especially for the patients with metastatic, recurrent, and unresectable forms of the disease. Therefore, the development of novel effective drugs interfering with the microtubules' dynamic state remains a big challenge in current oncology. We report here about the novel ethyl 2-amino-1-(furan-2-carboxamido)-5-(2-aryl/tert-butyl-2-oxoethylidene)-4-oxo-4,5-dihydro-1H-pyrrole-3-carboxylates (EAPCs) exhibiting potent anti-cancer activities against the breast and lung cancer cell lines in vitro. This was due to their ability to inhibit tubulin polymerization and induce cell cycle arrest in M-phase. As an outcome, the EAPC-treated cancer cells exhibited a significant increase in apoptosis, which was evidenced by the expression of cleaved forms of PARP, caspase-3, and increased numbers of Annexin-V-positive cells. By using the in silico molecular modeling methods (e.g., induced-fit docking, binding metadynamics, and unbiased molecular dynamics), we found that EAPC-67 and -70 preferentially bind to the colchicine-binding site of tubulin. Lastly, we have shown that the EAPCs indicated above and colchicine utilizes a similar molecular mechanism to inhibit tubulin polymerization via targeting the T7 loop in the β-chain of tubulin, thereby preventing the conformational changes in the tubulin dimers required for their polymerization. Collectively, we identified the novel and potent TBAs that bind to the colchicine-binding site and disrupt the microtubule network. As a result of these events, the compounds induced a robust cell cycle arrest in M-phase and exhibited potent pro-apoptotic activities against the epithelial cancer cell lines in vitro.

Dual-Inhibition of Human N-Myristoyltransferase Subtypes Halts Common Cold Pathogenesis: Atomistic Perspectives from the Case of IMP-1088

The pharmacological inhibition of human N-myristoyltransferase (HsNMT) has emerged as an efficient strategy to completely prevent the replication process of rhinoviruses, a potential treatment for the common cold. This was corroborated by the recent discovery of compound IMP-1088, a novel inhibitor that demonstrated a dual-inhibitory activity against the two HsNMT subtypes 1 and 2 without inducing cytotoxicity. However, the molecular and structural basis for the dual-inhibitory potential of IMP-1088 has not been investigated. As such, we employ molecular modelling techniques to resolve the structural mechanisms that account for the dual-inhibitory prowess of IMP-1088. Sequence and nanosecond-based analyses identified Tyr296, Phe190, Tyr420, Leu453, Gln496, Val181, Leu474, Glu182, and Asn246 as residues common within the binding pockets of both HsNMT1 and HsNMT2 subtypes whose consistent interactions with IMP-1088 underpin the basis for its dual inhibitory potency. Nano-second-based assessment of interaction dynamics revealed that Tyr296 consistently elicited high-affinity π-π stacked interaction with IMP-1088, thus further highlighting its cruciality corroborating previous report. An exploration of resulting structural changes upon IMP-1088 binding further revealed a characteristic impeding of residue fluctuations, structural compactness, and a consequential burial of crucial hydrophobic residues, features required for HsNMT1/2 functionality. Findings present essential structural perspectives that augment previous experimental efforts and could also advance drug development for treating respiratory tract infections, especially those mediated by rhinoviruses.

Positioning of an unprecedented 1,5-oxaza spiroquinone scaffold into SMYD2 inhibitors in epigenetic space

Lysine methyltransferases are important regulators of epigenetic signaling and are emerging as a novel drug target for drug discovery. This work demonstrates the positioning of novel 1,5-oxaza spiroquinone scaffold into selective SET and MYND domain-containing proteins 2 methyltransferases inhibitors. Selectivity of the scaffold was identified by epigenetic target screening followed by SAR study for the scaffold. The optimization was performed iteratively by two-step optimization consisting of iterative synthesis and computational studies (docking, metadynamics simulations). Computational binding studies guided the important interactions of the spiro[5.5]undeca scaffold in pocket 1 and Lysine channel and suggested extension of tail length for the improvement of potency (IC₅₀: up to 399 nM). The effective performance of cell proliferation assay for chosen compounds (IC₅₀: up to 11.9 nM) led to further evaluation in xenograft assay. The potent compound 24 demonstrated desirable in vivo efficacy with growth inhibition rate of 77.7% (4 fold decrease of tumor weight and 3 fold decrease of tumor volume). Moreover, mirosomal assay and pharmacokinetic profile suggested further developability of this scaffold through the identification of major metabolites (dealkylation at silyl group, reversible hydration product, the absence of toxic quinone fragments) and enough exposure of the testing compound 24 in plasma. Such spiro[5.5]undeca framework or ring system was neither been reported nor suggested as a modulator of methyltransferases. The chemo-centric target positioning and structural novelty can lead to potential pharmacological benefit.