Virtual Screening of Chemical Compounds for Discovery of Complement C3 Ligands

In silico analyses of acetylcholinesterase (AChE) and its genetic variants in interaction with the anti-Alzheimer drug Rivastigmine

Alzheimer's disease (AD) is the leading cause of dementia worldwide. Despite causing great social and economic impact, there is currently no cure for AD. The most effective therapy to manage AD symptoms is based on acetylcholinesterase inhibitors (AChEi), from which rivastigmine presented numerous benefits. However, mutations in AChE, which affect approximately 5% of the population, can modify protein structure and function, changing the individual response to Alzheimer's treatment. In this study, we performed computer simulations of AChE wild type and variants R34Q, P135A, V333E, and H353N, identified by one or more genome-wide association studies, to evaluate their effects on protein structure and interaction with rivastigmine. The functional effects of AChE variants were predicted using eight machine learning algorithms, while the evolutionary conservation of AChE residues was analyzed using the ConSurf server. Autodock4.2.6 was used to predict the binding modes for the hAChE-rivastigmine complex, which is still unknown. Molecular dynamics (MD) simulations were performed in triplicates for the AChE wild type and mutants using the GROMACS packages. Among the analyzed variants, P135A was classified as deleterious by all the functional prediction algorithms, in addition to occurring at highly conserved positions, which may have harmful consequences on protein function. The molecular docking results suggested that rivastigmine interacts with hAChE at the upper active-site gorge, which was further confirmed by MD simulations. Our MD findings also suggested that the complex hAChE-rivastigmine remains stable over time. The essential dynamics revealed flexibility alterations at the active-site gorge upon mutations P135A, V333E, and H353N, which may lead to strong and nonintuitive consequences to hAChE binding. Nonetheless, similar binding affinities were registered in the MMPBSA analysis for the hAChE wild type and variants when complexed to rivastigmine. Finally, our findings indicated that the rivastigmine binding to hAChE is an energetically favorable process mainly driven by negatively charged amino acids.

Computational Study on Potential Novel Anti-Ebola Virus Protein VP35 Natural Compounds

Ebola virus (EBOV) is one of the most lethal pathogens that can infect humans. The Ebola viral protein VP35 (EBOV VP35) inhibits host IFN-α/β production by interfering with host immune responses to viral invasion and is thus considered as a plausible drug target. The aim of this study was to identify potential novel lead compounds against EBOV VP35 using computational techniques in drug discovery. The 3D structure of the EBOV VP35 with PDB ID: 3FKE was used for molecular docking studies. An integrated library of 7675 African natural product was pre-filtered using ADMET risk, with a threshold of 7 and, as a result, 1470 ligands were obtained for the downstream molecular docking using AutoDock Vina, after an energy minimization of the protein via GROMACS. Five known inhibitors, namely, amodiaquine, chloroquine, gossypetin, taxifolin and EGCG were used as standard control compounds for this study. The area under the curve (AUC) value, evaluating the docking protocol obtained from the receiver operating characteristic (ROC) curve, generated was 0.72, which was considered to be acceptable. The four identified potential lead compounds of NANPDB4048, NANPDB2412, ZINC000095486250 and NANPDB2476 had binding affinities of −8.2, −8.2, −8.1 and −8.0 kcal/mol, respectively, and were predicted to possess desirable antiviral activity including the inhibition of RNA synthesis and membrane permeability, with the probable activity (Pa) being greater than the probable inactivity (Pi) values. The predicted anti-EBOV inhibition efficiency values (IC₅₀), found using a random forest classifier, ranged from 3.35 to 11.99 μM, while the Ki values ranged from 0.97 to 1.37 μM. The compounds NANPDB4048 and NANPDB2412 had the lowest binding energy of −8.2 kcal/mol, implying a higher binding affinity to EBOV VP35 which was greater than those of the known inhibitors. The compounds were predicted to possess a low toxicity risk and to possess reasonably good pharmacological profiles. Molecular dynamics (MD) simulations of the protein–ligand complexes, lasting 50 ns, and molecular mechanisms Poisson-Boltzmann surface area (MM-PBSA) calculations corroborated the binding affinities of the identified compounds and identified novel critical interacting residues. The antiviral potential of the molecules could be confirmed experimentally, while the scaffolds could be optimized for the design of future novel anti-EBOV chemotherapeutics.

Network Pharmacology-Based Study to Uncover Potential Pharmacological Mechanisms of Korean Thistle (Cirsium japonicum var. maackii (Maxim.) Matsum.) Flower against Cancer

Cirsium japonicum var. maackii (Maxim.) Matsum. or Korean thistle flower is a herbal plant used to treat tumors in Korean folk remedies, but its essential bioactives and pharmacological mechanisms against cancer have remained unexplored. This study identified the main compounds(s) and mechanism(s) of the C. maackii flower against cancer via network pharmacology. The bioactives from the C. maackii flower were revealed by gas chromatography-mass spectrum (GC-MS), and SwissADME evaluated their physicochemical properties. Next, target(s) associated with the obtained bioactives or cancer-related targets were retrieved by public databases, and the Venn diagram selected the overlapping targets. The networks between overlapping targets and bioactives were visualized, constructed, and analyzed by RPackage. Finally, we implemented a molecular docking test (MDT) to explore key target(s) and compound(s) on AutoDockVina and LigPlot+. GC-MS detected a total of 34 bioactives and all were accepted by Lipinski's rules and therefore classified as drug-like compounds (DLCs). A total of 597 bioactive-related targets and 4245 cancer-related targets were identified from public databases. The final 51 overlapping targets were selected between the bioactive targets network and cancer-related targets. With Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment, a total of 20 signaling pathways were manifested, and a hub signaling pathway (PI3K-Akt signaling pathway), a key target (Akt1), and a key compound (Urs-12-en-24-oic acid, 3-oxo, methyl ester) were selected among the 20 signaling pathways via MDT. Overall, Urs-12-en-24-oic acid, 3-oxo, methyl ester from the C. maackii flower has potent anti-cancer efficacy by inactivating Akt1 on the PI3K-Akt signaling pathway.

Phytochemical Moieties From Indian Traditional Medicine for Targeting Dual Hotspots on SARS-CoV-2 Spike Protein: An Integrative in-silico Approach

SARS-CoV-2 infection across the world has led to immense turbulence in the treatment modality, thus demanding a swift drug discovery process. Spike protein of SARS-CoV-2 binds to ACE2 receptor of human to initiate host invasion. Plethora of studies demonstrate the inhibition of Spike-ACE2 interactions to impair infection. The ancient Indian traditional medicine has been of great interest of Virologists worldwide to decipher potential antivirals. Hence, in this study, phytochemicals (1,952 compounds) from eight potential medicinal plants used in Indian traditional medicine were meticulously collated, based on their usage in respiratory disorders, along with immunomodulatory and anti-viral potential from contemporary literature. Further, these compounds were virtually screened against Receptor Binding Domain (RBD) of Spike protein. The potential compounds from each plant were prioritized based on the binding affinity, key hotspot interactions at ACE2 binding region and glycosylation sites. Finally, the potential hits in complex with spike protein were subjected to Molecular Dynamics simulation (450 ns), to infer the stability of complex formation. Among the compounds screened, Tellimagrandin-II (binding energy of −8.2 kcal/mol and binding free energy of −32.08 kcal/mol) from Syzygium aromaticum L. and O-Demethyl-demethoxy-curcumin (binding energy of −8.0 kcal/mol and binding free energy of −12.48 kcal/mol) from Curcuma longa L. were found to be highly potential due to their higher binding affinity and significant binding free energy (MM-PBSA), along with favorable ADMET properties and stable intermolecular interactions with hotspots (including the ASN343 glycosylation site). The proposed hits are highly promising, as these are resultant of stringent in silico checkpoints, traditionally used, and are documented through contemporary literature. Hence, could serve as promising leads for subsequent experimental validations.

Computational Screening of Potential Inhibitors of Desulfobacter postgatei for Pyrite Scale Prevention in Oil and Gas Wells

Sulfate-reducing bacteria (SRB), such as Desulfobacter postgatei are found in oil wells. However, they lead to the release of hydrogen sulfide. This in turn leads to the iron sulfide scale formation (pyrite). ATP sulfurylase is an enzyme present in SRB, which catalyzes the formation of adenylyl sulfate (APS) and inorganic pyrophosphatase (PPi) from ATP and sulfate. This reaction is the first among many in hydrogen sulfide production by D. postgatei . Consensus scoring using molecular docking and machine learning was used to identify three potential inhibitors of ATP sulfurylase from a database of about 40 million compounds. These selected hits ((S,E)-1-(4-methoxyphenyl)-3-(9-((m-tolylimino)methyl)-9,10-dihydroanthracen-9-yl)pyrrolidine-2,5-dione; methyl 2-[[(1S)-5-cyano-2-imino-1-(4-phenylthiazol-2-yl)-3-azaspiro[5.5]undec-4-en-4-yl]sulfanyl]acetate; and (4S)-4-(3-chloro-4-hydroxy-phenyl)-1-(6-hydroxypyridazin-3-yl)-3-methyl-4,5-dihydropyrazolo[3,4-b]pyridin-6-ol), known as A, B, and C, respectively) all had good binding affinities with ATP sulfurylase and were further analyzed for their toxicological properties. Compound A had the highest docking score. However, based on the physicochemical and toxicological properties, only compound C was predicted to be both safe and effective as a potential inhibitor of ATP sulfurylase, hence the preferred choice. The molecular interactions of compound C revealed favorable interactions with the following residues: LEU213, ASP308, ARG307, TRP347, LEU224, GLN212, MET211, and HIS309.

A network pharmacology study on main chemical compounds from Hibiscus cannabinus L. leaves

Hibiscus cannabinus L. leaves (HCLLs) are considered a favorable source of natural antiobesity substances. However, actual bioactive compound(s) in it and their mechanism(s) against obesity have not been confirmed. Hence, network pharmacology was conducted to identify its key compounds and mechanism(s) against obesity. Compounds in HCLLs were identified through GC-MS analysis and screened by Lipinski's rule. Genes related to the selected compounds and obesity were obtained from public databases, and overlapping genes between HCLL compound-related genes and obesity target genes were selected using a Venn diagram. The networking between selected compounds and overlapping genes was then constructed, visualized, and analyzed by RStudio. Finally, the binding affinity between compounds and genes was evaluated via molecular docking (MD). A total of 30 compounds in HCLLs were detected via GC-MS, and Lipinski's rule accepted all compounds. The compound-related genes (570 genes) and obesity targeted genes (3028 genes) were identified, and between them, 64 overlapping genes were selected. Gene Set Enrichment Analysis (GSEA) displayed that the mechanisms of HCLLs against obesity were associated with 13 signaling pathways on 22 compounds in HCLLs. Superficially, AKT1, vitamin E, and RAS signaling pathways were noted as a hub gene, an uppermost bioactive compound, and a hub signaling pathway, respectively. However, the binding affinity of ligands and proteins on the RAS signaling pathway was very low; instead, the PPAR signalling pathway was evaluated with potent efficacy against obesity through MD. On the PPAR signaling pathway, α-amyrin was found as the most significant compound for the amelioration of obesity. α-Amyrin manifested the strongest binding affinity on six target proteins associated with the PPAR signaling pathway. Our study suggests that an auxiliary (PPAR) signaling pathway of HCLLs might intervene efficiently against obesity over the hub (RAS) signaling pathway.

Computational analysis of complement inhibitor compstatin using molecular dynamics

The complement system plays a major role in human immunity, but its abnormal activation can have severe pathological impacts. By mimicking a natural mechanism of complement regulation, the small peptide compstatin has proven to be a very promising complement inhibitor. Over the years, several compstatin analogs have been created, with improved inhibitory potency. A recent analog is being developed as a candidate drug against several pathological conditions, including COVID-19. However, the reasons behind its higher potency and increased binding affinity to complement proteins are not fully clear. This computational study highlights the mechanistic properties of several compstatin analogs, thus complementing previous experimental studies. We perform molecular dynamics simulations involving six analogs alone in solution and two complexes with compstatin bound to complement component 3. These simulations reveal that all the analogs we consider, except the original compstatin, naturally adopt a pre-bound conformation in solution. Interestingly, this set of analogs adopting a pre-bound conformation includes analogs that were not known to benefit from this behavior. We also show that the most recent compstatin analog (among those we consider) forms a stronger hydrogen bond network with its complement receptor than an earlier analog.

Interactions between Curcumin Derivatives and Amyloid-β Fibrils: Insights from Molecular Dynamics Simulations

The aggregation of amyloid-β (Aβ) peptides into senile plaques is a hallmark of Alzheimer's disease (AD) and is hypothesized to be the primary cause of AD related neurodegeneration. Previous studies have shown the ability of curcumin to both inhibit the aggregation of Aβ peptides into oligomers or fibrils and reduce amyloids in vivo. Despite the promise of curcumin and its derivatives to serve as diagnostic, preventative, and potentially therapeutic AD molecules, the mechanism by which curcumin and its derivatives bind to and inhibit Aβ fibrils' formation remains elusive. Here, we investigated curcumin and a set of curcumin derivatives in complex with a hexamer peptide model of the Aβ1-42 fibril using nearly exhaustive docking, followed by multi-ns molecular dynamics simulations, to provide atomistic-detail insights into the molecules' binding and inhibitory properties. In the vast majority of the simulations, curcumin and its derivatives remain firmly bound in complex with the fibril through primarily three different principle binding modes, in which the molecules interact with residue domain 17LVFFA21, in line with previous experiments. In a small subset of these simulations, the molecules partly dissociate the outermost peptide of the Aβ1-42 fibril by disrupting β-sheets within the residue domain 12VHHQKLVFF20. A comparison between binding modes leading or not leading to partial dissociation of the outermost peptide suggests that the latter is attributed to a few subtle key structural and energetic interaction-based differences. Interestingly, partial dissociation appears to be either an outcome of high affinity interactions or a cause leading to high affinity interactions between the molecules and the fibril, which could partly serve as a compensation for the energy loss in the fibril due to partial dissociation. In conjunction with this, we suggest a potential inhibition mechanism of Αβ1-42 aggregation by the molecules, where the partially dissociated 16KLVFF20 domain of the outermost peptide could either remain unstructured or wrap around to form intramolecular interactions with the same peptide's 29GAIIG33 domain, while the molecules could additionally act as a patch against the external edge of the second outermost peptide's 16KLVFF20 domain. Thereby, individually or concurrently, these could prohibit fibril elongation.

Molecular Mechanism for Attractant Signaling to DHMA by E. coli Tsr

The attractant chemotaxis response of Escherichia coli to norepinephrine requires that it be converted to 3,4-dihydroxymandelic acid (DHMA) by the monoamine oxidase TynA and the aromatic aldehyde dehydrogenase FeaB. DHMA is sensed by the serine chemoreceptor Tsr, and the attractant response requires that at least one subunit of the periplasmic domain of the Tsr homodimer (pTsr) has an intact serine-binding site. DHMA that is generated in vivo by E. coli is expected to be a racemic mixture of the (R) and (S) enantiomers, so it has been unclear whether one or both chiral forms are active. Here, we used a combination of state-of-the-art tools in molecular docking and simulations, including an in-house simulation-based docking protocol, to investigate the binding properties of (R)-DHMA and (S)-DHMA to E. coli pTsr. Our studies computationally predicted that (R)-DHMA should promote a stronger attractant response than (S)-DHMA because of a consistently greater-magnitude piston-like pushdown of the pTsr α-helix 4 toward the membrane upon binding of (R)-DHMA than upon binding of (S)-DHMA. This displacement is caused primarily by interaction of DHMA with Tsr residue Thr156, which has been shown by genetic studies to be critical for the attractant response to L-serine and DHMA. These findings led us to separate the two chiral species and test their effectiveness as chemoattractants. Both the tethered cell and motility migration coefficient assays validated the prediction that (R)-DHMA is a stronger attractant than (S)-DHMA. Our study demonstrates that refined computational docking and simulation studies combined with experiments can be used to investigate situations in which subtle differences between ligands may lead to diverse chemotactic responses.

Isoflavones as Ah Receptor Agonists in Colon-Derived Cell Lines: Structure-Activity Relationships

Many of the protective responses observed for flavonoids in the gastrointestinal track resemble aryl hydrocarbon receptor (AhR)-mediated effects. Therefore, we examined the structure-activity relationships of isoflavones and isomeric flavone and flavanones as AhR ligands on the basis of their induction of CYP1A1, CYP1B1, and UGT1A1 gene expression in colon cancer Caco2 cells and young adult mouse colonocyte (YAMC) cells. Caco2 cells were significantly more Ah-responsive than YAMC cells, and this was due, in part, to flavonoid-induced cytotoxicity in the latter cell lines. The structure-activity relationships for the flavonoids were complex and both response and cell context specific; however, there was significant variability in the AhR activities of the isomeric substituted isoflavones and flavones. For example, 4',5,7-trihydroxyisoflavone (genistein) was AhR-inactive whereas 4',5,7-trihydroxyflavone (apigenin) induced CYP1A1, CYP1B1, and UGT1A1 in Caco2 cells. In contrast, both 5,7-dihydroxy-4-methoxy substituted isoflavone (biochanin A) and flavone (acacetin) induced all three AhR-responsive genes; 4',5,7-trimethoxyisoflavone was a potent AhR agonist, and the isomeric flavone was AhR-inactive. These results coupled with simulation studies modeling flavonoid interaction within the AhR binding pocket demonstrate that the orientation of the substituted phenyl ring at C-2 (flavones) or C-3 (isoflavones) on the common 4-H-chromen-4-one ring strongly influences the activities of isoflavones and flavones as AhR agonists.

Activation of COUP-TFI by a Novel Diindolylmethane Derivative

Chicken ovalbumin upstream promoter-transcription factor I (COUP-TFI) is an orphan receptor and member of the nuclear receptor superfamily. Among a series of methylene substituted diindolylmethanes (C-DIMs) containing substituted phenyl and heteroaromatic groups, we identified 1,1-bis(3'-indolyl)-1-(4-pyridyl)-methane (DIM-C-Pyr-4) as an activator of COUP-TFI. Structure activity studies with structurally diverse heteroaromatic C-DIMs showed that the pyridyl substituted compound was active and the 4-pyridyl substituent was more potent than the 2- or 3-pyridyl analogs in transactivation assays in breast cancer cells. The DIM-C-Pyr-4 activated chimeric GAL4-COUP-TFI constructs containing full length, C- or N-terminal deletions, and transactivation was inhibited by phosphatidylinositol-3-kinase and protein kinase A inhibitors. However, DIM-C-Pyr-4 also induced transactivation and interactions of COUP-TFI and steroid receptor coactivators-1 and -2 in mammalian two-hybrid assays, and ligand-induced interactions of the C-terminal region of COUP-TFI were not affected by kinase inhibitors. We also showed that DIM-C-Pyr-4 activated COUP-TFI-dependent early growth response 1 (Egr-1) expression and this response primarily involved COUP-TFI interactions with Sp3 and to a lesser extent Sp1 bound to the proximal region of the Egr-1 promoter. Modeling studies showed interactions of DIM-C-Pyr-4 within the ligand binding domain of COUP-TFI. This report is the first to identify a COUP-TFI agonist and demonstrate activation of COUP-TFI-dependent Egr-1 expression.