Speed vs Accuracy: Effect on Ligand Pose Accuracy of Varying Box Size and Exhaustiveness in AutoDock Vina

Integrated network pharmacology and cellular assay reveal the biological mechanisms of Limonium sinense (Girard) Kuntze against Breast cancer

BACKGROUND

Limonium Sinense (Girard) Kuntze (L. sinense) has been widely used for the treatment of anaemia, bleeding, cancer, and other disorders in Chinese folk medicine. The aim of this study is to predict the therapeutic effects of L. sinense and investigate the potential mechanisms using integrated network pharmacology methods and in vitro cellular experiments.

METHODS

The active ingredients of L. sinense were collected from published literature, and the potential targets related to L. sinense were obtained from public databases. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and DisGeNET enrichment analyses were performed to explore the underlying mechanisms. Molecular docking, cellular experiments, RNA-sequencing (RNA-seq) and Gene Expression Omnibus (GEO) datasets were employed to further evaluate the findings.

RESULTS

A total of 15 active ingredients of L. sinense and their corresponding 389 targets were obtained. KEGG enrichment analysis revealed that the biological effects of L. sinense were primarily associated with "Pathways in cancer". DisGeNET enrichment analysis highlighted the potential role of L. sinense in the treatment of breast cancer. Apigenin within L. sinense showed promising potential against cancer. Cellular experiments demonstrated that the L. sinense ethanol extract (LSE) exhibited a significant growth inhibitory effect on multiple breast cancer cell lines in both 2D and 3D cultures. RNA-seq analysis revealed a potential impact of LSE on breast cancer. Additionally, analysis of GEO datasets verified the significant enrichment of breast cancer and several cancer-related pathways upon treatment with Apigenin in human breast cancer cells.

CONCLUSION

This study predicts the biological activities of L. sinense and demonstrates the inhibitory effect of LSE on breast cancer cells, highlighting the potential application of L. sinense in cancer treatment.

Synthesis of Pyrrolo[3,4-b]pyridin-5-ones via Ugi-Zhu Reaction and In Vitro-In Silico Studies against Breast Carcinoma

An Ugi-Zhu three-component reaction (UZ-3CR) coupled in a one-pot manner to a cascade process (N-acylation/aza Diels-Alder cycloaddition/decarboxylation/dehydration) was performed to synthesize a series of pyrrolo[3,4-b]pyridin-5-ones in 20% to 92% overall yields using ytterbium triflate as a catalyst, toluene as a solvent, and microwaves as a heat source. The synthesized molecules were evaluated in vitro against breast cancer cell lines MDA-MB-231 and MCF-7, finding that compound 1f, at a concentration of 6.25 μM, exhibited a potential cytotoxic effect. Then, to understand the interactions between synthesized compounds and the main proteins related to the cancer cell lines, docking studies were performed on the serine/threonine kinase 1 (AKT1) and Orexetine type 2 receptor (Ox2R), finding moderate to strong binding energies, which matched accurately with the in vitro results. Additionally, molecular dynamics were performed between proteins related to the studied cell lines and the three best ligands.

ReverseDock: a web server for blind docking of a single ligand to multiple protein targets using AutoDock Vina

Several platforms exist to perform molecular docking to computationally predict binders to a specific protein target from a library of ligands. The reverse, that is, docking a single ligand to various protein targets, can currently be done by very few web servers, which limits the search to a small set of pre-selected human proteins. However, the possibility to in silico predict which targets a compound identified in a high-throughput drug screen bind would help optimize and reduce the costs of the experimental workflow needed to reveal the molecular mechanism of action of a ligand. Here, we present ReverseDock, a blind docking web server based on AutoDock Vina specifically designed to allow users with no computational expertise to dock a ligand to 100 protein structures of their choice. ReverseDock increases the number and type of proteins a ligand can be docked to, making the task of in silico docking of a ligand to entire families of proteins straightforward. We envision ReverseDock will support researchers by providing the possibility to apply inverse docking computations using web browser. ReverseDock is available at: https://reversedock.biologie.uni-freiburg.de/.

New Benzamides as Multi-Targeted Compounds: A Study on Synthesis, AChE and BACE1 Inhibitory Activity and Molecular Docking

The synthesis of eleven new and previously undescribed benzamides was designed. These compounds were specifically projected as potential inhibitors of the enzymes acetylcholinesterase (AChE) and β-secretase (BACE1). N,N'-(1,4-phenylene)bis(3-methoxybenzamide) was most active against AChE, with an inhibitory concentration of AChE IC50 = 0.056 µM, while the IC50 for donepezil was 0.046 µM. This compound was also the most active against the BACE1 enzyme. The IC50 value was 9.01 µM compared to that for quercetin, with IC50 = 4.89 µM. Quantitative results identified this derivative to be the most promising. Molecular modeling was performed to elucidate the potential mechanism of action of this compound. Dynamic simulations showed that new ligands only had a limited stabilizing effect on AChE, but all clearly reduced the flexibility of the enzyme. It can, therefore, be concluded that a possible mechanism of inhibition increases the stiffness and decreases the flexibility of the enzyme, which obviously impedes its proper function. An analysis of the H-bonding patterns suggests a different mechanism (from other ligands) when interacting the most active derivative with the enzyme.

Neuroprotective Properties of Oleanolic Acid-Computational-Driven Molecular Research Combined with In Vitro and In Vivo Experiments

Oleanolic acid (OA), as a ubiquitous compound in the plant kingdom, is studied for both its neuroprotective and neurotoxic properties. The mechanism of acetylcholinesterase (AChE) inhibitory potential of OA is investigated using molecular dynamic simulations (MD) and docking as well as biomimetic tests. Moreover, the in vitro SH-SY5Y human neuroblastoma cells and the in vivo zebrafish model were used. The inhibitory potential towards the AChE enzyme is examined using the TLC-bioautography assay (the IC50 value is 9.22 μM). The CH-π interactions between the central fragment of the ligand molecule and the aromatic cluster created by the His440, Phe288, Phe290, Phe330, Phe331, Tyr121, Tyr334, Trp84, and Trp279 side chains are observed. The results of the in vitro tests using the SH-SY5Y cells indicate that the viability rate is reduced to 71.5%, 61%, and 43% at the concentrations of 100 µg/mL, 300 µg/mL, and 1000 µg/mL, respectively, after 48 h of incubation, whereas cytotoxicity against the tested cell line with the IC50 value is 714.32 ± 32.40 µg/mL. The in vivo tests on the zebrafish prove that there is no difference between the control and experimental groups regarding the mortality rate and morphology (p > 0.05).

Deciphering Molecular Aspects of Potential α-Glucosidase Inhibitors within Aspergillus terreus: A Computational Odyssey of Molecular Docking-Coupled Dynamics Simulations and Pharmacokinetic Profiling

Hyperglycemia, as a hallmark of the metabolic malady diabetes mellitus, has been an overwhelming healthcare burden owing to its high rates of comorbidity and mortality, as well as prospective complications affecting different body organs. Available therapeutic agents, with α-glucosidase inhibitors as one of their cornerstone arsenal, control stages of broad glycemia while showing definitive characteristics related to their low clinical efficiency and off-target complications. This has propelled the academia and industrial section into discovering novel and safer candidates. Herein, we provided a thorough computational exploration of identifying candidates from the marine-derived Aspergillus terreus isolates. Combined structural- and ligand-based approaches using a chemical library of 275 metabolites were adopted for pinpointing promising α-glucosidase inhibitors, as well as providing guiding insights for further lead optimization and development. Structure-based virtual screening through escalating precision molecular docking protocol at the α-glucosidase canonical pocket identified 11 promising top-docked hits, with several being superior to the market drug reference, acarbose. Comprehensive ligand-based investigations of these hits' pharmacokinetics ADME profiles, physiochemical characterizations, and obedience to the gold standard Lipinski's rule of five, as well as toxicity and mutagenicity profiling, proceeded. Under explicit conditions, a molecular dynamics simulation identified the top-stable metabolites: butyrolactone VI (SK-44), aspulvinone E (SK-55), butyrolactone I 4''''-sulfate (SK-72), and terrelumamide B (SK-173). They depicted the highest free binding energies and steadiest thermodynamic behavior. Moreover, great structural insights have been revealed, including the advent of an aromatic scaffold-based interaction for ligand-target complex stability. The significance of introducing balanced hydrophobic/polar moieties, like triazole and other bioisosteres of carboxylic acid, has been highlighted across docking, ADME/Tox profiling, and molecular dynamics studies for maximizing binding interactions while assuring safety and optimal pharmacokinetics for targeting the intestinal-localized α-glucosidase enzyme. Overall, this study provided valuable starting points for developing new α-glucosidase inhibitors based on nature-derived unique scaffolds, as well as guidance for prospective lead optimization and development within future pre-clinical and clinical investigations.

In Silico Identification and In Vitro Validation of Repurposed Compounds Targeting the RSV Polymerase

Respiratory Syncytial Virus (RSV) is the top cause of infant hospitalization globally, with no effective treatments available. Researchers have sought small molecules to target the RNA-dependent RNA Polymerase (RdRP) of RSV, which is essential for replication and transcription. Based on the cryo-EM structure of the RSV polymerase, in silico computational analysis including molecular docking and the protein-ligand simulation of a database, including 6554 molecules, is currently undergoing phases 1-4 of clinical trials and has resulted in the top ten repurposed compound candidates against the RSV polymerase, including Micafungin, Totrombopag, and Verubecestat. We performed the same procedure to evaluate 18 small molecules from previous studies and chose the top four compounds for comparison. Among the top identified repurposed compounds, Micafungin, an antifungal medication, showed significant inhibition and binding affinity improvements over current inhibitors such as ALS-8112 and Ribavirin. We also validated Micafungin's inhibition of the RSV RdRP using an in vitro transcription assay. These findings contribute to RSV drug development and hold promise for broad-spectrum antivirals targeting the non-segmented negative-sense (NNS) RNA viral polymerases, including those of rabies (RABV) and Ebola (EBOV).

A Study Based on Network Pharmacology Decoding the Multi-Target Mechanism of Duhuo Jisheng Decoction for the Treatment of Intervertebral Disc Degeneration

Intervertebral disc degeneration (IDD) poses a grim public health impact. Duhuo Jisheng Decoction (DJD), a traditional Chinese medicine formula, has recently received significant attention for its efficacy and safety in treating IDD. However, the pathological processes of IDD in which DJD interferes and molecular mechanism involved are poorly understood, which brings difficulties to the clinical practice of DJD for the treatment of IDD. This study systematically investigated the underlying mechanism of DJD treatment of IDD. Network pharmacology approaches were employed, integrating molecular docking and random walk with restart (RWR) algorithm, to identify key compounds and targets for DJD in the treatment of IDD. Bioinformatics approaches were used to further explore the biological insights in DJD treatment of IDD. The analysis identifies AKT1, PIK3R1, CHUK, ALB, TP53, MYC, NR3C1, IL1B, ERBB2, CAV1, CTNNB1, AR, IGF2, and ESR1 as key targets. Responses to mechanical stress, oxidative stress, cellular inflammatory responses, autophagy, and apoptosis are identified as the critical biological processes involved in DJD treatment of IDD. The regulation of DJD targets in extracellular matrix components, ion channel regulation, transcriptional regulation, synthesis and metabolic regulation of reactive oxygen products in the respiratory chain and mitochondria, fatty acid oxidation, the metabolism of Arachidonic acid, and regulation of Rho and Ras protein activation are found to be potential mechanisms in disc tissue response to mechanical stress and oxidative stress. MAPK, PI3K/AKT, and NF-κB signaling pathways are identified as vital signaling pathways for DJD to treat IDD. Quercetin and Kaempferol are assigned a central position in the treatment of IDD. This study contributes to a more comprehensive understanding of the mechanism of DJD in treating IDD. It provides a reference for applying natural products to delay the pathological process of IDD.

Structural analysis of human G-protein-coupled receptor 17 ligand binding sites

The human G protein coupled membrane receptor (GPR17), the sensor of brain damage, is identified as a biomarker for many neurological diseases. In human brain tissue, GPR17 exist in two isoforms, long and short. While cryo-electron microscopy technology has provided the structure of the long isoform of GPR17 with Gi complex, the structure of the short isoform and its activation mechanism remains unclear. Recently, we theoretically modeled the structure of the short isoform of GPR17 with Gi signaling protein and identified novel ligands. In the present work, we demonstrated the presence of two distinct ligand binding sites in the short isoform of GPR17. The molecular docking of GPR17 with endogenous (UDP) and synthetic ligands (T0510.3657, MDL29950) found the presence of two distinct binding pockets. Our observations revealed that endogenous ligand UDP can bind stronger in two different binding pockets as evidenced by glide and autodock vina scores, whereas the other two ligand's binding with GPR17 has less docking score. The analysis of receptor-UDP interactions shows complexes' stability in the lipid environment by 100 ns atomic molecular dynamics simulations. The amino acid residues VAL83, ARG87, and PHE111 constitute ligand binding site 1, whereas site 2 constitutes ASN67, ARG129, and LYS232. Root mean square fluctuation analysis showed the residues 83, 87, and 232 with higher fluctuations during molecular dynamics simulation in both binding pockets. Our findings imply that the residues of GPR17's two binding sites are crucial, and their interaction with UDP reveals the protein's hidden signaling and communication properties. Furthermore, this finding may assist in the development of targeted therapies for the treatment of neurological diseases.

Approaching the Dimerization Mechanism of Small Molecule Inhibitors Targeting PD-L1 with Molecular Simulation

Inhibitors blocking the PD-1/PD-L1 immune checkpoint demonstrate impressive anti-tumor immunity, and small molecule inhibitors disclosed by the Bristol-Myers Squibb (BMS) company have become a hot topic. In this work, by modifying the carbonyl group of BMS-202 into a hydroxyl group to achieve two enantiomers (MS and MR) with a chiral center, we found that this is an effective way to regulate its hydrophobicity and thus to reduce the negative effect of polar solvation free energy, which enhances the stability of PD-L1 dimer/inhibitor complexes. Moreover, we studied the binding modes of BMS-200 and BMS-202-related small molecule inhibitors by molecular dynamics simulation to explore their inhibitory mechanism targeting PD-L1 dimerization. The results showed that the size exclusion effect of the inhibitors triggered the rearrangement of the residue _ATyr56, leading to the formation of an axisymmetric tunnel-shaped pocket, which is an important structural basis for improving the binding affinity of symmetric inhibitors with PD-L1. Furthermore, after inhibitor dissociation, the conformation of _ATyr123 and _BMet115 rearranged, which blocked the entrance of the binding pocket, while the reverse rearrangements of the same residues occurred when the PD-L1 monomer was complexed with the inhibitors, preparing PD-L1 for dimerization. Overall, this study casts a new light on the inhibitory mechanism of BMS inhibitors targeting PD-L1 dimerization and provides an idea for designing novel small molecule inhibitors for future cancer immunotherapy.