ProteinsPlus: a web portal for structure analysis of macromolecules

New pyridine-based chalcones and pyrazolines with anticancer, antibacterial, and antiplasmodial activities

New pyridine-based chalcones 4a-h and pyrazolines 5a-h (N-acetyl), 6a-h (N-phenyl), and 7a-h (N-4-chlorophenyl) were synthesized and evaluated by the National Cancer Institute (NCI) against 60 different human cancer cell lines. Pyrazolines 6a, 6c-h, and 7a-h satisfied the pre-determined threshold inhibition criteria, obtaining that compounds 6c and 6f exhibited high antiproliferative activity, reaching submicromolar GI50 values from 0.38 to 0.45 μM, respectively. Moreover, compound 7g (4-CH3) exhibited the highest cytostatic activity of these series against different cancer cell lines from leukemia, nonsmall cell lung, colon, ovarian, renal, and prostate cancer, with LC50 values ranging from 5.41 to 8.35 μM, showing better cytotoxic activity than doxorubicin. Furthermore, the compounds were tested for antibacterial and antiplasmodial activities. Chalcone 4c was the most active with minimal inhibitory concentration (MIC) = 2 μg/mL against methicillin-resistant Staphylococcus aureus (MRSA), while the pyrazoline 6h showed a MIC = 8 μg/mL against Neisseria gonorrhoeae. For anti-Plasmodium falciparum activity, the chalcones display higher activity with EC50 values ranging from 10.26 to 10.94 μg/mL. Docking studies were conducted against relevant proteins from P. falciparum, exhibiting the minimum binding energy with plasmepsin II. In vivo toxicity assay in Galleria mellonella suggests that most compounds are low or nontoxic.

Qijiao Shengbai Capsule alleviated leukopenia by interfering leukotriene pathway: Integrated network study of multi-omics

BACKGROUND

Leukopenia could be induced by chemotherapy, which leads to bone marrow suppression and even affects the therapeutic progression of cancer. Qijiao Shengbai Capsule (QSC) has been used for the treatment of leukopenia in clinic, but its bioactive components and mechanisms have not yet been elucidated clearly.

PURPOSE

This study aimed to elucidate the molecular mechanisms of QSC in treating leukopenia.

STUDY DESIGN

Serum pharmacochemistry, multi-omics, network pharmacology, and validation experiment were combined to study the effect of QSC in murine leukopenia model.

METHODS

First, UPLC-QTOF-MS was used to clarify the absorbed components of QSC. Then, cyclophosphamide (CTX) was used to induce mice model with leukopenia, and the therapeutic efficacy of QSC was assessed by an integrative approach of multi-omics and network pharmacology strategy. Finally, molecular mechanisms and potential therapeutic targets were identified by validated experiments.

RESULTS

121 compounds absorbed in vivo were identified. QSC significantly increase the count of white blood cells (WBCs) in peripheral blood of leukopenia mice with 15 days treatment. Multi-omics and network pharmacology revealed that leukotriene pathway and MAPK signaling pathway played crucial roles during the treatment of leukopenia with QSC. Six targets (ALOX5, LTB4R, CYSLTR1, FOS, JUN, IL-1β) and 13 prototype compounds were supposed to be the key targets and potential active components, respectively. The validation experiment further confirmed that QSC could effectively inhibit the inflammatory response induced by leukopenia. The inhibitors of ALOX5 activity can significantly increase the number of WBCs in leukopenia mice. Molecular docking of ALOX5 suggested that calycosin, daidzein, and medicarpin were the potentially active compounds of QSC.

CONCLUSION

Leukotriene pathway was found for the first time to be a key role in the development of leukopenia, and ALOX5 was conformed as the potential target. QSC may inhibit the inflammatory response and interfere the leukotriene pathway, it is able to improve hematopoiesis and achieve therapeutic effects in the mice with leukopenia.

User-centric design of a 3D search interface for protein-ligand complexes

In this work, we present the frontend of GeoMine and showcase its application, focusing on the new features of its latest version. GeoMine is a search engine for ligand-bound and predicted empty binding sites in the Protein Data Bank. In addition to its basic text-based search functionalities, GeoMine offers a geometric query type for searching binding sites with a specific relative spatial arrangement of chemical features such as heavy atoms and intermolecular interactions. In contrast to a text search that requires simple and easy-to-formulate user input, a 3D input is more complex, and its specification can be challenging for users. GeoMine's new version aims to address this issue from the graphical user interface perspective by introducing an additional visualization concept and a new query template type. In its latest version, GeoMine extends its query-building capabilities primarily through input formulation in 2D. The 2D editor is fully synchronized with GeoMine's 3D editor and provides the same functionality. It enables template-free query generation and template-based query selection directly in 2D pose diagrams. In addition, the query generation with the 3D editor now supports predicted empty binding sites for AlphaFold structures as query templates. GeoMine is freely accessible on the ProteinsPlus web server ( https://proteins.plus ).

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

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

Prenylated Flavanone Isolated from Dalea Species as a Potential Multitarget-Neuroprotector in an In Vitro Alzheimer's Disease Mice Model

Alzheimer's disease (AD) involves a neurodegenerative process that has not yet been prevented, reversed, or stopped. Continuing with the search for natural pharmacological treatments, flavonoids are a family of compounds with proven neuroprotective effects and multi-targeting behavior. The American genus Dalea L. (Fabaceae) is an important source of bioactive flavonoids. In this opportunity, we tested the neuroprotective potential of three prenylated flavanones isolated from Dalea species in a new in vitro pre-clinical AD model previously developed by us. Our approach consisted in exposing neural cells to conditioned media (3xTg-AD ACM) from neurotoxic astrocytes derived from hippocampi and cortices of old 3xTg-AD mice, mimicking a local neurodegenerative microenvironment. Flavanone 1 and 3 showed a neuroprotective effect against 3xTg-AD ACM, being 1 more active than 3. The structural requirements to afford neuroprotective activity in this model are a 5'-dimethylallyl and 4'-hydroxy at the B ring. In order to search the mechanistic performance of the most active flavanone, we focus on the flavonoid-mediated regulation of GSK-3β-mediated tau phosphorylation previously reported. Flavanone 1 treatment decreased the rise of hyperphosphorylated tau protein neuronal levels induced after 3xTg-AD ACM exposure and inhibited the activity of GSK-3β. Finally, direct exposure of these neurotoxic 3xTg-AD astrocytes to flavanone 1 resulted in toxicity to these cells and reduced the neurotoxicity of 3xTg-AD ACM as well. Our results allow us to present compound 1 as a natural prenylated flavanone that could be used as a precursor to development and design of future drug therapies for AD.

Identifying the Key Role of Plutella xylostella General Odorant Binding Protein 2 in Perceiving a Larval Attractant, (E,E)-2,6-Farnesol

There is currently a lack of effective olfaction-based techniques to control diamondback moth (DBM) larvae. Identifying behaviorally active odorants for DBM larvae and exploring their recognition mechanisms can provide insights into olfaction-based larval control strategies. Through the two-choice assay, (E,E)-2,6-farnesol (farnesol) was identified as a compound exhibiting significant attractant activity toward DBM larvae, achieving an attraction index of 0.48 ± 0.13. PxylGOBP1 and PxylGOBP2, highly expressed in the antennae of DBM larvae, both showed high affinity toward farnesol. RNAi technology was used to knock down PxylGOBP1 and PxylGOBP2, revealing that the attraction of DBM larvae to farnesol nearly vanished following the knockdown of PxylGOBP2, indicating its critical role in recognizing farnesol. Further investigation into the PxylGOBP2-farnesol interaction revealed the importance of residues like Thr9, Trp37, and Phe118 in PxylGOBP2's binding to farnesol. This research is significant for unveiling the olfactory mechanisms of DBM larvae and developing larval behavior regulation techniques.

Highly Sensitive and Rapid Screening Technique for the Detection of Organophosphate Pesticides and Copper Compounds Using Bifunctional Recombinant TrxA-PvCarE1

Enabling the detection of organophosphate pesticide (OP) residues through enzyme inhibition-based technology is crucial for ensuring food safety and human health. However, the use of acetylcholinesterase, the primary target enzyme for OPs, isolated from animals in practical production poses challenges in terms of sensitivity and batch stability. To address this issue, we identified a highly sensitive and reproducible biorecognition element, TrxA-PvCarE1, derived from red kidney beans and successfully overexpressed it in Escherichia coli. The resulting recombinant TrxA-PvCarE1 exhibited remarkable sensitivity toward 10 OPs, surpassing that of commercial acetylcholinesterase. Additionally, this approach demonstrated the capability to simultaneously detect copper compounds with high sensitivity, expanding the range of pesticides detectable using the traditional enzyme inhibition method. Spiking recovery tests conducted on cowpea and carrot samples verified the suitability of the TrxA-PvCarE1-based technique for real-life sample analysis. In summary, this study highlights a promising comprehensive candidate for the rapid detection of pesticide residues.

Investigating the association between CYP2J2 inhibitors and QT prolongation: a literature review

Drug withdrawal post-marketing due to cardiotoxicity is a major concern for drug developers, regulatory agencies, and patients. One common mechanism of cardiotoxicity is through inhibition of cardiac ion channels, leading to prolongation of the QT interval and sometimes fatal arrythmias. Recently, oxylipin signaling compounds have been shown to bind to and alter ion channel function, and disruption in their cardiac levels may contribute to QT prolongation. Cytochrome P450 2J2 (CYP2J2) is the predominant CYP isoform expressed in cardiomyocytes, where it oxidizes arachidonic acid to cardioprotective epoxyeicosatrienoic acids (EETs). In addition to roles in vasodilation and angiogenesis, EETs bind to and activate various ion channels. CYP2J2 inhibition can lower EET levels and decrease their ability to preserve cardiac rhythm. In this review, we investigated the ability of known CYP inhibitors to cause QT prolongation using Certara's Drug Interaction Database. We discovered that among the multiple CYP isozymes, CYP2J2 inhibitors were more likely to also be QT-prolonging drugs (by approximately 2-fold). We explored potential binding interactions between these inhibitors and CYP2J2 using molecular docking and identified four amino acid residues (Phe61, Ala223, Asn231, and Leu402) predicted to interact with QT-prolonging drugs. The four residues are located near the opening of egress channel 2, highlighting the potential importance of this channel in CYP2J2 binding and inhibition. These findings suggest that if a drug inhibits CYP2J2 and interacts with one of these four residues, then it may have a higher risk of QT prolongation and more preclinical studies are warranted to assess cardiovascular safety.

Investigation of Zika virus methyl transferase inhibitors using steered molecular dynamics

Zika virus (ZIKV) spread is considered a major public health threat by the World Health Organization (WHO). There are no vaccines or drugs available to control the infection of the Zika virus, therefore a highly effective medicinal molecule is urgently required. In this study, a computationally intensive investigation was performed to identify a potent natural compound that could inhibit the ZIKV NS5 methyltransferase. This research approach is based on target-based drug identification principles where the native inhibitor SAH (S-adenosylhomocysteine) of ZIKV NS5 methyltransferase was selected as a reference. High-throughput virtual screening and tanimoto similarity coefficient were applied to the natural compound library for ranking the potential candidates. The top five compounds were selected for interaction analysis, MD simulation, total binding free energy through MM/GBSA, and steered MD simulation. Among these compounds, Adenosine 5'-monophosphate monohydrate, Tubercidin, and 5-Iodotubercidin showed stable binding to the protein compared to the native compound, SAH. These three compounds also showed less fluctuations in RMSF in contrast to native compound. Additionally, the same interacting residues observed in SAH also made strong interactions with these three compounds. Adenosine 5'-monophosphate monohydrate and 5-Iodotubercidin had greater total binding free energies than the reference ligand. Moreover, the dissociation resistance of all three compounds was equivalent to that of the reference ligand. This study suggested binding properties of three-hit compounds that could be used to develop drugs against Zika virus infections.Communicated by Ramaswamy H. Sarma.

Machine Learning Application for Medicinal Chemistry: Colchicine Case, New Structures, and Anticancer Activity Prediction

In the contemporary era, the exploration of machine learning (ML) has gained widespread attention and is being leveraged to augment traditional methodologies in quantitative structure-activity relationship (QSAR) investigations. The principal objective of this research was to assess the anticancer potential of colchicine-based compounds across five distinct cell lines. This research endeavor ultimately sought to construct ML models proficient in forecasting anticancer activity as quantified by the IC50 value, while concurrently generating innovative colchicine-derived compounds. The resistance index (RI) is computed to evaluate the drug resistance exhibited by LoVo/DX cells relative to LoVo cancer cell lines. Meanwhile, the selectivity index (SI) is computed to determine the potential of a compound to demonstrate superior efficacy against tumor cells compared to its toxicity against normal cells, such as BALB/3T3. We introduce a novel ML system adept at recommending novel chemical structures predicated on known anticancer activity. Our investigation entailed the assessment of inhibitory capabilities across five cell lines, employing predictive models utilizing various algorithms, including random forest, decision tree, support vector machines, k-nearest neighbors, and multiple linear regression. The most proficient model, as determined by quality metrics, was employed to predict the anticancer activity of novel colchicine-based compounds. This methodological approach yielded the establishment of a library encompassing new colchicine-based compounds, each assigned an IC50 value. Additionally, this study resulted in the development of a validated predictive model, capable of reasonably estimating IC50 values based on molecular structure input.

α-L-Fucosidases from an Alpaca Faeces Metagenome: Characterisation of Hydrolytic and Transfucosylation Potential

In various life forms, fucose-containing glycans play vital roles in immune recognition, developmental processes, plant immunity, and host-microbe interactions. Together with glucose, galactose, N-acetylglucosamine, and sialic acid, fucose is a significant component of human milk oligosaccharides (HMOs). Fucosylated HMOs benefit infants by acting as prebiotics, preventing pathogen attachment, and potentially protecting against infections, including HIV. Although the need for fucosylated derivatives is clear, their availability is limited. Therefore, synthesis methods for various fucosylated oligosaccharides are explored, employing enzymatic approaches and α-L-fucosidases. This work aimed to characterise α-L-fucosidases identified in an alpaca faeces metagenome. Based on bioinformatic analyses, they were confirmed as members of the GH29A subfamily. The recombinant α-L-fucosidases were expressed in Escherichia coli and showed hydrolytic activity towards p-nitrophenyl-α-L-fucopyranoside and 2'-fucosyllactose. Furthermore, the enzymes' biochemical properties and kinetic characteristics were also determined. All four α-L-fucosidases could catalyse transfucosylation using a broad diversity of fucosyl acceptor substrates, including lactose, maltotriose, L-serine, and L-threonine. The results contribute insights into the potential use of α-L-fucosidases for synthesising fucosylated amino acids.

Redocking the PDB

Molecular docking is a standard technique in structure-based drug design (SBDD). It aims to predict the 3D structure of a small molecule in the binding site of a receptor (often a protein). Despite being a common technique, it often necessitates multiple tools and involves manual steps. Here, we present the JAMDA preprocessing and docking workflow that is easy to use and allows fully automated docking. We evaluate the JAMDA docking workflow on binding sites extracted from the complete PDB and derive key factors determining JAMDA's docking performance. With that, we try to remove most of the bias due to manual intervention and provide a realistic estimate of the redocking performance of our JAMDA preprocessing and docking workflow for any PDB structure. On this large PDBScan22 data set, our JAMDA workflow finds a pose with an RMSD of at most 2 Å to the crystal ligand on the top rank for 30.1% of the structures. When applying objective structure quality filters to the PDBScan22 data set, the success rate increases to 61.8%. Given the prepared structures from the JAMDA preprocessing pipeline, both JAMDA and the widely used AutoDock Vina perform comparably on this filtered data set (the PDBScan22-HQ data set).

Ligand 1 H NMR Chemical Shifts as Accurate Reporters for Protein-Ligand Binding Interfaces in Solution

The availability of high-resolution 3D structural information is crucial for investigating guest-host systems across a wide range of fields. In the context of drug discovery, the information is routinely used to establish and validate structure-activity relationships, grow initial hits from screening campaigns, and to guide molecular docking. For the generation of protein-ligand complex structural information, X-ray crystallography is the experimental method of choice, however, with limited information on protein flexibility. An experimentally verified structural model of the binding interface in the native solution-state would support medicinal chemists in their molecular design decisions. Here we demonstrate that protein-bound ligand 1 H NMR chemical shifts are highly sensitive and accurate probes for the immediate chemical environment of protein-ligand interfaces. By comparing the experimental ligand 1 H chemical shift values with those computed from the X-ray structure using quantum mechanics methodology, we identify significant disagreements for parts of the ligand between the two experimental techniques. We show that quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) ensembles can be used to refine initial X-ray co-crystal structures resulting in a better agreement with experimental 1 H ligand chemical shift values. Overall, our findings highlight the usefulness of ligand 1 H NMR chemical shift information in combination with a QM/MM MD workflow for generating protein-ligand ensembles that accurately reproduce solution structural data.

In silico prediction and in vitro assessment of novel heterocyclics with antimalarial activity

The development of new antimalarials is paramount to keep the goals on reduction of malaria cases in endemic regions. The search for quality hits has been challenging as many inhibitory molecules may not progress to the next development stage. The aim of this work was to screen an in-house library of heterocyclic compounds (HCUV) for antimalarial activity combining computational predictions and phenotypic techniques to find quality hits. The physicochemical determinants, pharmacokinetic properties (ADME), and drug-likeness of HCUV were evaluated in silico, and compounds were selected for structure-based virtual screening and in vitro analysis. Seven Plasmodium target proteins were selected from the DrugBank Database, and ligands and receptors were processed using UCSF Chimera and Open Babel before being subjected to docking using Autodock Vina and Autodock 4. Growth inhibition of P. falciparum (3D7) cultures was tested by SYBR Green assays, and toxicity was assessed using hemolytic activity tests and the Galleria mellonella in vivo model. From a total of 792 compounds, 341 with good ADME properties, drug-likeness, and no interference structures were subjected to in vitro analysis. Eight compounds showed IC50 ranging from 0.175 to 0.990 µM, and active compounds included pyridyl-diaminopyrimido-diazepines, pyridyl-N-acetyl- and pyridyl-N-phenyl-pyrazoline derivatives. The most potent compound (UV802, IC50 0.178 µM) showed no toxicophoric and was predicted to interact with P. falciparum 1-cysperoxidredoxin (PfPrx1). For the remaining 7 hits (IC50 < 1 μM), 3 showed in silico binding to PfPrx1, one was predicted to bind the haloacid dehalogenase-like hydrolase and plasmepsin II, and one interacted with the plasmodial heat shock protein 90.

Design, synthesis, and antitumor screening of new thiazole, thiazolopyrimidine, and thiazolotriazine derivatives as potent inhibitors of VEGFR-2

New thiazole, thiazolopyrimidine, and thiazolotriazine derivatives 3-12 and 14a-f were synthesized. The newly synthesized analogs were tested for in vitro antitumor activity against HepG2, HCT-116, MCF-7, HeP-2, and Hela cancer cells. Results indicated that compound 5 displayed the highest potency toward the tested cancer cells. Compound 11b possessed enhanced effectiveness over MCF-7, HepG2, HCT-116, and Hela cancer cells. In addition, compounds 4 and 6 showed promising activity toward HCT-116, MCF-7, and Hela cancer cells and eminent activity against HepG2 and HeP-2 cells. Moreover, compounds 3-6 and 11b were tested for their capability to inhibit vascular endothelial growth factor receptor-2 (VEGFR-2) activity. The obtained results showed that compound 5 displayed significant inhibitory activity against VEGFR-2 (half-maximal inhibitory concentration [IC50 ] = 0.044 μM) comparable to sunitinib (IC50  = 0.100 μM). Also, the synthesized compounds 3-6 and 11b were subjected to in vitro cytotoxicity tests over WI38 and WISH normal cells. It was found that the five tested compounds displayed significantly lower cytotoxicity than doxorubicin toward normal cell lines. Cell cycle analysis proved that compound 5 induces cell cycle arrest in the S phase for HCT-116 and Hela cancer cell lines and in the G2/M phase for the MCF-7 cancer cell line. Moreover, compound 5 induced cancer cell death through apoptosis accompanied by a high ratio of BAX/BCL-2 in the screened cancer cells. Furthermore, docking results revealed that compound 5 showed the essential interaction bonds with VEGFR-2, which agreed with in vitro enzyme assay results. In silico studies showed that most of the analyzed compounds complied with the requirements of good oral bioavailability with minimal toxicity threats in humans.

Unlocking therapeutic potential of trigonelline through molecular docking as a promising approach for treating diverse neurological disorders

Neurological disorders pose significant challenges in terms of treatment options, necessitating the exploration of novel therapeutic approaches. Trigonelline, a naturally occurring alkaloid found in various plants, has emerged as a potential treatment option. It has also been reported that trigonelline is involved in several pathways like; Oxidative Stress and Antioxidant, Inflammatory, Neuroprotection and Neurotrophic, Mitochondrial Function and Energy Metabolism. This study aims to investigate the therapeutic potential of trigonelline for diverse neurological disorders using a molecular docking approach. Molecular docking simulations were performed to predict the binding affinity and interaction between trigonelline and target proteins implicated in neurological disorders. The structural requirements for effective binding were also explored. The molecular docking results revealed strong binding interactions and favorable binding affinities between trigonelline and the target proteins involved in diverse neurological disorders like Alzheimer's disease, Parkinson's disease, epilepsy, and depression etc. The predicted binding modes provided insights into the key molecular interactions governing the ligand-protein complexes. The findings suggest that trigonelline holds promise as a therapeutic approach for several neurological disorders. The molecular docking approach employed in this study provides a valuable tool for rational drug design and optimization of trigonelline-based compounds. Further experimental validation and preclinical studies are warranted to confirm the efficacy and safety of trigonelline as a potential treatment option, paving the way for the development of more effective and targeted therapies for neurological disorders.

Rutin from Begonia roxburghii modulates iNOS and Sep A activity in treatment of Shigella flexneri induced diarrhoea in rats: An in vitro, in vivo and computational analysis

In developing countries, diarrhoea is a major issue of concern, where consistent use of antibiotics has resulted in several side effects along with development of resistance among pathogens against these antibiotics. Since natural products are becoming the treatment of choice, therefore present investigation involves mechanistic evaluation of antidiarrhoeal potential of Begonia roxburghii and its marker rutin against Shigella flexneri (SF) induced diarrhoea in rats following in vitro, in vivo and in silico protocols. The roots of the plant are used as vegetable in the North East India and are also used traditionally in treating diarrhoea. Phytochemically standardized ethanolic extract of B. roxburghii (EBR) roots and its marker rutin were first subjected to in vitro antibacterial evaluation against SF. Diarrhoea was induced in rats using suspension of SF and various diarrhoeagenic parameters were examined after first, third and fifth day of treatment at 100, 200 and 300 mg/kg, p.o. with EBR and 50 mg/kg, p.o. with rutin respectively. Additionally, density of SF in stools, stool water content, haematological and biochemical parameters, cytokine profiling, ion concentration, histopathology and Na+/K+-ATPase activity were also performed. Molecular docking and dynamics simulation studies of ligand rutin was studied against secreted extracellular protein A (Sep A, PDB: 5J44) from SF and Inducible nitric oxide synthase (iNOS, PDB: 1DD7) followed by network pharmacology. EBR and rutin demonstrated a potent antibacterial activity against SF and also showed significant recovery from diarrhoea (EBR: 81.29 ± 0.91% and rutin: 75.27 ± 0.89%) in rats after five days of treatment. EBR and rutin also showed significant decline in SF density in stools, decreased cytokine expression, potential antioxidant activity, cellular proliferative nature and recovered ion loss due to enhanced Na+/K+-ATPase activity, which was also supported by histopathology. Rutin showed a very high docking score of -11.61 and -9.98 kcal/mol against iNOS and Sep A respectively and their stable complex was also confirmed through dynamics, while network pharmacology suggested that, rutin is quite capable of modulating the pathways of iNOS and Sep A. Thus, we may presume that rutin played a key role in the observed antidiarrhoeal activity of B. roxburghii against SF induced diarrhoea.

An arginine-to-histidine mutation in flavanone-3-hydroxylase results in pink strawberry fruits

Fruit color is a very important external commodity factor for consumers. Compared to the most typical red octoploid strawberry (Fragaria × ananassa), the pink strawberry often sells for a more expensive price and has a higher economic benefit due to its outstanding color. However, few studies have examined the molecular basis of pink-colored strawberry fruit. Through an EMS mutagenesis of woodland strawberry (Fragaria vesca), we identified a mutant with pink fruits and green petioles. Bulked-segregant analysis sequencing analysis and gene function verification confirmed that the responsible mutation resides in a gene encoding flavanone-3-hydroxylase (F3H) in the anthocyanin synthesis pathway. This nonsynonymous mutation results in an arginine-to-histidine change at position 130 of F3H. Molecular docking experiments showed that the arginine-to-histidine mutation results in a reduction of intermolecular force-hydrogen bonding between the F3H protein and its substrates. Enzymatic experiments showed a greatly reduced ability of the mutated F3H protein to catalyze the conversion of the substrates and hence a blockage of the anthocyanin synthesis pathway. The discovery of a key residue in the F3H gene controlling anthocyanin synthesis provides a clear target of modification for the molecular breeding of strawberry varieties with pink-colored fruits, which may be of great commercial value.

Antcin-B, a phytosterol-like compound from Taiwanofungus camphoratus inhibits SARS-CoV-2 3-chymotrypsin-like protease (3CLPro) activity in silico and in vitro

Despite the remarkable development of highly effective vaccines, including mRNA-based vaccines, within a limited timeframe, coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is not been entirely eradicated. Thus, it is crucial to identify new effective anti-3CLPro compounds, pivotal for the replication of SARS-CoV-2. Here, we identified an antcin-B phytosterol-like compound from Taiwanofungus camphoratus that targets 3CLPro activity. MTT assay and ADMET prediction are employed for assessing potential cytotoxicity. Computational molecular modeling was used to screen various antcins and non-antcins for binding affinity and interaction type with 3CLPro. Further, these compounds were subjected to study their inhibitory effects on 3CLPro activity in vitro. Our results indicate that antcin-B has the best binding affinity by contacting residues like Leu141, Asn142, Glu166, and His163 via hydrogen bond and salt bridge and significantly inhibits 3CLPro activity, surpassing the positive control compound (GC376). The 100 ns molecular dynamics simulation studies showed that antcin-B formed consistent, long-lasting water bridges with Glu166 for their inhibitory activity. In summary, antcin-B could be useful to develop therapeutically viable drugs to inhibit SARS-CoV-2 replication alone or in combination with medications specific to other SARS-CoV-2 viral targets.

Assessing the inhibitory potential of anti-dengue compounds against Japanese encephalitis virus RNA dependent RNA polymerase: an in silico study

Japanese encephalitis (JE), a neurological infection of severe nature, is caused by the Japanese encephalitis virus (JEV) and is transmitted by the mosquito vector. The polymerase domain of Non-structural 5 (NS5), which is also referred to as RdRp (RNA-dependent RNA polymerase), is considered a potential therapeutic target for JEV. The present study employed molecular dynamics modelling and high-throughput virtual screening to evaluate the possible antiviral activity of anti-dengue drugs against JEV RdRp. Furthermore, a ranking was performed utilising the MM/GBSA analysis to identify the three most promising compounds. Compound ID 57409246 exhibited the highest binding affinity with the protein, as evidenced by its minimum binding free energy of -72.96 kcal/mole. In contrast, the other two compounds had minimum binding free energies of -67.57 and -59.19 kcal/mole, respectively. Upon conducting a 100 nanosecond molecular dynamics simulation to confirm the binding of the chemical complexes, it was observed that the three hits, namely 57409246, 70683874, and 44577154, exhibited a consistent and stable RMSD. Subsequently, the binding strength of the trajectory was confirmed through MM/GBSA analysis. The compounds 70683874 and 57409246 exhibited the lowest binding free energies, which were -97.58 kcal/mol and -96.38 kcal/mol, respectively. The binding free energy (ΔG Bind) values for the native ligand ATP and molecule 44577154 were -65.64 kcal/mol and -69.44 kcal/mol, respectively. Overall, compared to the native ligand ATP, all three compounds exhibited higher binding affinity. The study proposes three anti-dengue molecules as a potential remedy for JE, which can be confirmed through in vitro and in vivo investigations.Communicated by Ramaswamy H. Sarma.

PoseEdit: enhanced ligand binding mode communication by interactive 2D diagrams

In this article, we present PoseEdit, a new, interactive frontend of the popular pose visualization tool PoseView. PoseEdit automatically produces high-quality 2D diagrams of intermolecular interactions in 3D binding sites calculated from ligands in complex with protein, DNA, and RNA. The PoseView diagrams have been improved in several aspects, most notably in their interactivity. Thanks to the easy-to-use 2D editor of PoseEdit, the diagrams are extensively editable and extendible by the user, can be merged with other diagrams, and even be created from scratch. A large variety of graphical objects in the diagram can be moved, rotated, selected and highlighted, mirrored, removed, or even newly added. Furthermore, PoseEdit enables a synchronized 2D-3D view of macromolecule-ligand complexes simplifying the analysis of structural features and interactions. The representation of individual diagram objects regarding their visualized chemical properties, like stereochemistry, and general graphical styles, like the color of interactions, can additionally be edited. The primary objective of PoseEdit is to support scientists with an enhanced way to communicate ligand binding mode information through graphical 2D representations optimized with the scientist's input in accordance with objective criteria and individual needs. PoseEdit is freely available on the ProteinsPlus web server ( https://proteins.plus ).

Functional and structural insights into the multi-step activation and catalytic mechanism of bacterial ExoY nucleotidyl cyclase toxins bound to actin-profilin

ExoY virulence factors are members of a family of bacterial nucleotidyl cyclases (NCs) that are activated by specific eukaryotic cofactors and overproduce cyclic purine and pyrimidine nucleotides in host cells. ExoYs act as actin-activated NC toxins. Here, we explore the Vibrio nigripulchritudo Multifunctional-Autoprocessing Repeats-in-ToXin (MARTX) ExoY effector domain (Vn-ExoY) as a model for ExoY-type members that interact with monomeric (G-actin) instead of filamentous (F-actin) actin. Vn-ExoY exhibits moderate binding affinity to free or profilin-bound G-actin but can capture the G-actin:profilin complex, preventing its spontaneous or VASP- or formin-mediated assembly at F-actin barbed ends in vitro. This mechanism may prolong the activated cofactor-bound state of Vn-ExoY at sites of active actin cytoskeleton remodelling. We present a series of high-resolution crystal structures of nucleotide-free, 3'-deoxy-ATP- or 3'-deoxy-CTP-bound Vn-ExoY, activated by free or profilin-bound G-actin-ATP/-ADP, revealing that the cofactor only partially stabilises the nucleotide-binding pocket (NBP) of NC toxins. Substrate binding induces a large, previously-unidentified, closure of their NBP, confining catalytically important residues and metal cofactors around the substrate, and facilitating the recruitment of two metal ions to tightly coordinate the triphosphate moiety of purine or pyrimidine nucleotide substrates. We validate critical residues for both the purinyl and pyrimidinyl cyclase activity of NC toxins in Vn-ExoY and its distantly-related ExoY from Pseudomonas aeruginosa, which specifically interacts with F-actin. The data conclusively demonstrate that NC toxins employ a similar two-metal-ion mechanism for catalysing the cyclisation of nucleotides of different sizes. These structural insights into the dynamics of the actin-binding interface of actin-activated ExoYs and the multi-step activation of all NC toxins offer new perspectives for the specific inhibition of class II bacterial NC enzymes.

LifeSoaks: a tool for analyzing solvent channels in protein crystals and obstacles for soaking experiments

Due to the structural complexity of proteins, their corresponding crystal arrangements generally contain a significant amount of solvent-occupied space. These areas allow a certain degree of intracrystalline protein flexibility and mobility of solutes. Therefore, knowledge of the geometry of solvent-filled channels and cavities is essential whenever the dynamics inside a crystal are of interest. Especially in soaking experiments for structure-based drug design, ligands must be able to traverse the crystal solvent channels and reach the corresponding binding pockets. Unsuccessful screenings are sometimes attributed to the geometry of the crystal packing, but the underlying causes are often difficult to understand. This work presents LifeSoaks, a novel tool for analyzing and visualizing solvent channels in protein crystals. LifeSoaks uses a Voronoi diagram-based periodic channel representation which can be efficiently computed. The size and location of channel bottlenecks, which might hinder molecular diffusion, can be directly derived from this representation. This work presents the calculated bottleneck radii for all crystal structures in the PDB and the analysis of a new, hand-curated data set of structures obtained by soaking experiments. The results indicate that the consideration of bottleneck radii and the visual inspection of channels are beneficial for planning soaking experiments.

Fluoroquinolone Analogs, SAR Analysis, and the Antimicrobial Evaluation of 7-Benzimidazol-1-yl-fluoroquinolone in In Vitro, In Silico, and In Vivo Models

Structure-activity relationship (SAR) studies allow the evaluation of the relationship between structural chemical changes and biological activity. Fluoroquinolones have chemical characteristics that allow their structure to be modified and new analogs with different therapeutic properties to be generated. The objective of this research is to identify and select the C-7 heterocycle fluoroquinolone analog (FQH 1-5) with antibacterial activity similar to the reference fluoroquinolone through in vitro, in silico, and in vivo evaluations. First, SAR analysis was conducted on the FQH 1-5, using an in vitro antimicrobial sensibility model in order to select the best compound. Then, an in silico model mechanism of action analysis was carried out by molecular docking. The non-bacterial cell cytotoxicity was evaluated, and finally, the antimicrobial potential was determined by an in vivo model of topical infection in mice. The results showed antimicrobial differences between the FQH 1-5 and Gram-positive and Gram-negative bacteria, identifying the 7-benzimidazol-1-yl-fluoroquinolone (FQH-2) as the most active against S. aureus. Suggesting the same mechanism of action as the other fluoroquinolones; no cytotoxic effects on non-bacterial cells were found. FQH-2 was demonstrated to decrease the amount of bacteria in infected wound tissue.

P. Discovery of novel inhibitor of 11 beta-hydroxysteroid dehydrogenase type 1 using in silico structure-based screening approach for the treatment of type 2 diabetes

Purpose

The current study is aimed to perform structure-based screening of FDA-approved drugs that can act as novel inhibitor of the 11beta- hydroxysteroid dehydrogenase type 1 (11β-HSD1) enzyme.

Methods

Structural analogs of carbenoxolone (CBX) were selected from DrugBank database and their Absorption, Distribution, Metabolism, Excretion, and Toxicity (ADMET) parameters were investigated by SwissADME. Molecular docking of CBX analogs against 11β-HSD1 was performed by AutoDock tool, their binding patterns were visualized using PyMOL and the interacting amino acids were determined by ProteinPlus tool. Molecular dynamics simulation was performed on the docked structure of 11β-HSD1 (Protein Data Bank (PDB) code: 2ILT) using GROMACS 2018.1.

Results

The binding energies of hydrocortisone succinate, medroxyprogesterone acetate, testolactone, hydrocortisone cypionate, deoxycorticosterone acetate, and hydrocortisone probutate were lower than that of substrate corticosterone. The molecular dynamics simulation of 11β-HSD1 and hydrocortisone cypionate docked structure showed that it formed a stable complex with the inhibitor. The Root mean square deviation (RMSD) of the protein (0.37 ± 0.05 nm) and ligand (0.41 ± 0.06 nm) shows the stability of the ligand-protein interaction.

Conclusion

The docking study revealed that hydrocortisone cypionate has a higher binding affinity than carbenoxolone and its other analogs. The molecular dynamics simulation indicated the stability of the docked complex of 11β-HSD1 and hydrocortisone cypionate. These findings indicate the potential use of this FDA approved drug in the treatment of type 2 diabetes. However, validation by in vitro inhibitory studies and clinical trials on type 2 diabetes patients is essential to confirm the current findings.

Binding Site Detection Remastered: Enabling Fast, Robust, and Reliable Binding Site Detection and Descriptor Calculation with DoGSite3

Binding site prediction on protein structures is a crucial step in early phase drug discovery whenever experimental or predicted structure models are involved. DoGSite belongs to the widely used tools for this task. It is a grid-based method that uses a Difference-of-Gaussian filter to detect cavities on the protein surface. We recently reimplemented the first version of this method, released in 2010, focusing on improved binding site detection in the presence of ligands and optimized parameters for more robust, reliable, and fast predictions and binding site descriptor calculations. Here, we introduce the new version, DoGSite3, compare it to its predecessor, and re-evaluate DoGSite on published data sets for a large-scale comparative performance evaluation.

Triptolide alleviates collagen-induced arthritis in mice by modulating Treg/Th17 imbalance through the JAK/PTEN-STAT3 pathway

BACKGROUND

This study aimed to investigate the effects of triptolide (TP) on collagen-induced arthritis (CIA) mice and the related mechanisms.

METHODS

CIA mice were administered TP for 35 days. Mouse ankle joints and serum antibodies and cytokines were examined to assess the therapeutic effects of TP. The ratios of Treg, Th1, and Th17 cells were measured by flow cytometry and RT-qPCR. Reverse docking was used to characterize the binding modes of TP against target proteins. The expression of the STAT3 pathway in CIA mice was evaluated by western blotting and immunofluorescence staining. Mouse spleen lymphocytes were extracted and the expression of the STAT3 pathway after IL-6 stimulation was analyzed.

RESULTS

TP could significantly alleviate joint swelling, reduce bone destruction, and downregulate serum inflammation levels. TP improved the imbalance of Treg/Th17 cells in CIA mice. TP could form stable complexes with target proteins. TP significantly inhibited the activation of the JAK/PTEN-STAT3 pathway in mice. Moreover, TP regulated the activation of the JAK1/2-STAT3 signaling pathway in mouse spleen lymphocytes under inflammatory stimulation.

CONCLUSION

TP can inhibit inflammation and alleviate bone destruction in CIA mice. The underlying mechanism is related to the regulation of the imbalance of Treg/Th17 cells through the JAK/PTEN-STAT3 pathway.

Structural analysis of factors related to FAM3C/ILEI dimerization and identification of inhibitor candidates targeting cancer treatment

FAM3 is a superfamily of four cytokines that maintain a single globular structure β -β -α of three classes: FAM3A, B, C and D. FAM3C was the first member of this family related to cancer and is functionally characterized as an essential factor for the epithelial-mesenchymal transition (EMT), leading to late delays in tumor progression. Due to its crucial role in EMT and metastasis, FAM3C has been termed an interleukin-like EMT (ILEI) inducer. There are several studies on the part of FAM3C in the progression of cancer and other diseases. However, little is known about its cellular receptors and possible inhibitors. In this study, based on in silico approaches, we performed structural analyses of factors related to FAM3C/ILEI dimerization. We also identified four possible inhibitor candidates, expected to be exciting prototypes and could be submitted to future biological tests targeting cancer treatment.

Lymphocyte antigen 6K signaling to aurora kinase promotes advancement of the cell cycle and the growth of cancer cells, which is inhibited by LY6K-NSC243928 interaction

Lymphocyte antigen 6K (LY6K) is a small GPI-linked protein that is normally expressed in testes. Increased expression of LY6K is significantly associated with poor survival outcomes in many solid cancers, including cancers of the breast, ovary, gastrointestinal tract, head and neck, brain, bladder, and lung. LY6K is required for ERK-AKT and TGF-β pathways in cancer cells and is required for in vivo tumor growth. In this report, we describe a novel role for LY6K in mitosis and cytokinesis through aurora B kinase and its substrate histone H3 signaling axis. Further, we describe the structural basis of the molecular interaction of small molecule NSC243928 with LY6K protein and the disruption of LY6K-aurora B signaling in cell cycle progression due to LY6K-NSC243928 interaction. Overall, disruption of LY6K function via NSC243928 led to failed cytokinesis, multinucleated cells, DNA damage, senescence, and apoptosis of cancer cells. LY6K is not required for vital organ function, thus inhibition of LY6K signaling is an ideal therapeutic approach for hard-to-treat cancers that lack targeted therapy such as triple-negative breast cancer.

Identification of a new class of proteasome inhibitors based on a naphthyl-azotricyclic-urea-phenyl scaffold

Proteasomes play an important role in protein degradation and regulation of many cellular pathways by maintaining protein balance. Inhibitors of proteasomes disrupt this balance affecting proteins that are key in malignancies and as such have found applications in the treatment of multiple myeloma and mantle cell lymphoma. However, resistance mechanisms have been reported for these proteasome inhibitors including mutations at the β5 site which necessitates the constant development of new inhibitors. In this work, we report the identification of a new class of proteasome inhibitors, polycyclic molecules bearing a naphthyl-azotricyclic-urea-phenyl scaffold, from screening of the ZINC library of natural products. The most potent of these compounds showed evidence of dose dependency through proteasome assays with IC50 values in the low micromolar range, and kinetic analysis revealed competitive binding at the β5c site with an estimated inhibition constant, K i, of 1.15 μM. Inhibition was also shown for the β5i site of the immunoproteasome at levels similar to those of the constitutive proteasome. Structure-activity relationship studies identified the naphthyl substituent to be crucial for activity and this was attributed to enhanced hydrophobic interactions within β5c. Further to this, halogen substitution within the naphthyl ring enhanced the activity and allowed for π-π interactions with Y169 in β5c and Y130 and F124 in β5i. The combined data highlight the importance of hydrophobic and halogen interactions in β5 binding and assist in the design of next generation inhibitors of proteasomes.

Structure-Based Profiling of Potential Phytomolecules with AKT1 a Key Cancer Drug Target

Identifying cancer biomarkers is imperative, as upregulated genes offer a better microenvironment for the tumor; hence, targeted inhibition is preferred. The theme of our study is to predict molecular interactions between cancer biomarker proteins and selected natural compounds. We identified an overexpressed potential molecular target (AKT1) and computationally evaluated its inhibition by four dietary ligands (isoliquiritigenin, shogaol, tehranolide, and theophylline). The three-dimensional structures of protein and phytochemicals were retrieved from the RCSB PDB database (4EKL) and NCBI’s PubChem, respectively. Rational structure-based docking studies were performed using AutoDock. Results were analyzed based primarily on the estimated free binding energy (kcal/mol), hydrogen bonds, and inhibition constant, Ki, to identify the most effective anti-cancer phytomolecule. Toxicity and drug-likeliness prediction were performed using OSIRIS and SwissADME. Amongst the four phytocompounds, tehranolide has better potential to suppress the expression of AKT1 and could be used for anti-cancer drug development, as inhibition of AKT1 is directly associated with the inhibition of growth, progression, and metastasis of the tumor. Docking analyses reveal that tehranolide has the most efficiency in inhibiting AKT1 and has the potential to be used for the therapeutic management of cancer. Natural compounds targeting cancer biomarkers offer less rejection, minimal toxicity, and fewer side effects.

Exploration of potent antiviral phytomedicines from Lauraceae family plants against SARS-CoV-2 RNA-dependent RNA polymerase

RNA-dependent RNA polymerase, also known as RdRp, is a possible therapeutic target that could be used to suppress the proliferation of RNA viruses such as SARS-CoV-2. This protein has two major functional sites (a) catalytic and (b) substrate entry, which regulate the natural substrate entry and its corresponding interaction with the protein. In this study, a computational drug design pipeline was applied to investigate potential inhibitors against SARS-CoV-2 RdRp from Lauraceae plants, and five top hits were selected based on the docked score (< -7 kcal/mol). The docking study suggested that the Glochidioboside had a minimum binding score of -7.8 kcal/mol. This compound showed total five hydrogen bonds while two of them were with catalytic residues Asp618 and Asp760. However, another compound, Sitogluside showed a binding score of -7.3 kcal/mol with four hydrogen bonds targeting three functional residues (Arg555, Ser759, and Asp760). Later, 100 ns explicit solvent molecular dynamics (MD) simulation was performed to evaluate the stability of the protein-ligand docked system. These compounds translocated their positions from the catalytic site to the substrate entry site, as observed in the MD simulation trajectory. However, translocation did not affect the binding strength of these compounds, and they retained the strong binding affinity (ΔG < -11.5 kcal/mol), estimated using the MM/GBSA method. In general, the findings of this study indicated the potential therapeutic compounds that may be used targeting SARS-CoV-2 RdRp. However, these compounds still need to be validated by experimentation in order to determine their inhibitory function.Communicated by Ramaswamy H. Sarma.

Virtual Screening of Benzimidazole Derivatives as Potential Triose Phosphate Isomerase Inhibitors with Biological Activity against Leishmania mexicana

Leishmania mexicana (L. mexicana) is a causal agent of cutaneous leishmaniasis (CL), a “Neglected disease”, for which the search for new drugs is a priority. Benzimidazole is a scaffold used to develop antiparasitic drugs; therefore, it is interesting molecule against L. mexicana. In this work, a ligand-based virtual screening (LBVS) of the ZINC15 database was performed. Subsequently, molecular docking was used to predict the compounds with potential binding at the dimer interface of triosephosphate isomerase (TIM) of L. mexicana (LmTIM). Compounds were selected on binding patterns, cost, and commercial availability for in vitro assays against L. mexicana blood promastigotes. The compounds were analyzed by molecular dynamics simulation on LmTIM and its homologous human TIM. Finally, the physicochemical and pharmacokinetic properties were determined in silico. A total of 175 molecules with docking scores between −10.8 and −9.0 Kcal/mol were obtained. Compound E2 showed the best leishmanicidal activity (IC50 = 4.04 µM) with a value similar to the reference drug pentamidine (IC50 = 2.23 µM). Molecular dynamics analysis predicted low affinity for human TIM. Furthermore, the pharmacokinetic and toxicological properties of the compounds were suitable for developing new leishmanicidal agents.

Small molecule inhibiting microglial nitric oxide release could become a potential treatment for neuroinflammation

Microglia are the immune effector cells of the central nervous system (CNS) and react to pathologic events with a complex process including the release of nitric oxide (NO). NO is a free radical, which is toxic for all cells at high concentrations. To target an exaggerated NO release, we tested a library of 16 544 chemical compounds for their effect on lipopolysaccharide (LPS)-induced NO release in cell line and primary neonatal microglia. We identified a compound (C1) which significantly reduced NO release in a dose-dependent manner, with a low IC50 (252 nM) and no toxic side effects in vitro or in vivo. Target finding strategies such as in silico modelling and mass spectroscopy hint towards a direct interaction between C1 and the nitric oxide synthase making C1 a great candidate for specific intra-cellular interaction with the NO producing machinery.

Toxicological impacts and likely protein targets of bisphenol a in Paramecium caudatum

Bisphenol A (BPA) is a widely used plasticizer agent and a well-known ubiquitous endocrine disruptor, which is frequently associated with a series of reproductive, developmental, and transgenerational effects over wildlife, livestocks, and humans. Although extensive toxicological data is available for metazoans, the impact of BPA over unicellular eukaryotes, which represents a considerable proportion of eukaryotic diversity, remains largely overlooked. Here, we used acute end-point toxicological assay and an inverted virtual-screening (IVS) approach to evaluate cellular impairments infringed by BPA over the cosmopolitan ciliated protist, Paramecium caudatum. Our data indicate a clear time-dependent effect over P. caudatum survival, which seems to be a consequence of disruptions to multiple core cellular functions, such as DNA and cell replication, transcription, translation and signaling pathways. Finally, the use of this ciliate as a biosensor to monitor BPA within environments and the relevance of bioinformatic methods to leverage our current knowledge on the impacts of emerging contaminants to biological systems are discussed.

Comparative Genomics of Histoplasma capsulatum and Prediction of New Vaccines and Drug Targets

Histoplasma capsulatum is a thermodymorphic fungus that causes histoplasmosis, a systemic mycosis that presents different clinical manifestations, ranging from self-limiting to acute lung infection, chronic lung infection and disseminated infection. Usually, it affects severely immunocompromised patients although immunocompetent patients can also be infected. Currently, there are no vaccines to prevent histoplasmosis and the available antifungal treatment presents moderate to high toxicity. Additionally, there are few options of antifungal drugs. Thus, the aim of this study was to predict possible protein targets for the construction of potential vaccine candidates and predict potential drug targets against H. capsulatum. Whole genome sequences from four previously published H. capsulatum strains were analyzed and submitted to different bioinformatic approaches such as reverse vaccinology and subtractive genomics. A total of four proteins were characterized as good protein candidates (vaccine antigens) for vaccine development, three of which are membrane-bound and one is secreted. In addition, it was possible to predict four cytoplasmic proteins which were classified as good protein candidates and, through molecular docking performed for each identified target, we found four natural compounds that showed favorable interactions with our target proteins. Our study can help in the development of potential vaccines and new drugs that can change the current scenario of the treatment and prevention of histoplasmosis.

SARS-CoV-2 proteases Mpro and PLpro: Design of inhibitors with predicted high potency and low mammalian toxicity using artificial neural networks, ligand-protein docking, molecular dynamics simulations, and ADMET calculations

The main (Mpro) and papain-like (PLpro) proteases are highly conserved viral proteins essential for replication of the COVID-19 virus, SARS-COV-2. Therefore, a logical plan for producing new drugs against this pathogen is to discover inhibitors of these enzymes. Accordingly, the goal of the present work was to devise a computational approach to design, characterize, and select compounds predicted to be potent dual inhibitors - effective against both Mpro and PLpro. The first step employed LigDream, an artificial neural network, to create a virtual ligand library. Ligands with computed ADMET profiles indicating drug-like properties and low mammalian toxicity were selected for further study. Initial docking of these ligands into the active sites of Mpro and PLpro was done with GOLD, and the highest-scoring ligands were redocked with AutoDock Vina to determine binding free energies (ΔG). Compounds 89-00, 89-07, 89-32, and 89-38 exhibited favorable ΔG values for Mpro (-7.6 to -8.7 kcal/mol) and PLpro (-9.1 to -9.7 kcal/mol). Global docking of selected compounds with the Mpro dimer identified prospective allosteric inhibitors 89-00, 89-27, and 89-40 (ΔG -8.2 to -8.9 kcal/mol). Molecular dynamics simulations performed on Mpro and PLpro active site complexes with the four top-scoring ligands from Vina demonstrated that the most stable complexes were formed with compounds 89-32 and 89-38. Overall, the present computational strategy generated new compounds with predicted drug-like characteristics, low mammalian toxicity, and high inhibitory potencies against both target proteases to form stable complexes. Further preclinical studies will be required to validate the in silico findings before the lead compounds could be considered for clinical trials.

Water Networks in Complexes between Proteins and FDA-Approved Drugs

Water molecules at protein-ligand interfaces are often of significant pharmaceutical interest, owing in part to the entropy which can be released upon the displacement of an ordered water by a therapeutic compound. Protein structures may not, however, completely resolve all critical bound water molecules, or there may be no experimental data available. As such, predicting the location of water molecules in the absence of a crystal structure is important in the context of rational drug design. Grand canonical Monte Carlo (GCMC) is a computational technique that is gaining popularity for the simulation of buried water sites. In this work, we assess the ability of GCMC to accurately predict water binding locations, using a dataset that we have curated, containing 108 unique structures of complexes between proteins and Food and Drug Administration (FDA)-approved small-molecule drugs. We show that GCMC correctly predicts 81.4% of nonbulk crystallographic water sites to within 1.4 Å. However, our analysis demonstrates that the reported performance of water prediction methods is highly sensitive to the way in which the performance is measured. We also find that crystallographic water sites with more protein/ligand hydrogen bonds and stronger electron density are more reliably predicted by GCMC. An analysis of water networks revealed that more than half of the structures contain at least one ligand-contacting water network. In these cases, displacement of a water site by a ligand modification might yield unexpected results if the larger network is destabilized. Cooperative effects between waters should therefore be explicitly considered in structure-based drug design.

Anti-Entry Activity of Natural Flavonoids against SARS-CoV-2 by Targeting Spike RBD

COVID-19 is still a global public health concern, and the SARS-CoV-2 mutations require more effective antiviral agents. In this study, the antiviral entry activity of thirty-one flavonoids was systematically evaluated by a SARS-CoV-2 pseudovirus model. Twenty-four flavonoids exhibited antiviral entry activity with IC50 values ranging from 10.27 to 172.63 µM and SI values ranging from 2.33 to 48.69. The structure-activity relationship of these flavonoids as SARS-CoV-2 entry inhibitors was comprehensively summarized. A subsequent biolayer interferometry assay indicated that flavonoids bind to viral spike RBD to block viral interaction with ACE2 receptor, and a molecular docking study also revealed that flavonols could bind to Pocket 3, the non-mutant regions of SARS-CoV-2 variants, suggesting that flavonols might be also active against virus variants. These natural flavonoids showed very low cytotoxic effects on human normal cell lines. Our findings suggested that natural flavonoids might be potential antiviral entry agents against SARS-CoV-2 via inactivating the viral spike. It is hoped that our study will provide some encouraging evidence for the use of natural flavonoids as disinfectants to prevent viral infections.

Sperm activate TLR2/TLR1 heterodimerization to induce a weak proinflammatory response in the bovine uterus

Toll-like receptor 2 (TLR2) signaling pathway is involved in the sperm-triggered uterine inflammatory response at insemination, but its precise mechanism at molecular-level remains unknown. According to the ligand specificity, TLR2 forms a heterodimer with TLR1 or TLR6 as an initial step to mediate intracellular signaling, leading to a specific type of immune response. Hence, the present study aimed to identify the active TLR2 heterodimer (TLR2/1 or TLR2/6) that is involved in sperm-uterine immune crosstalk in bovine using various models. First, in-vitro (bovine endometrial epithelial cells, BEECs) and ex-vivo (bovine uterine explant) models were employed to test different TLR2 dimerization pathways in endometrial epithelia after exposure to sperm or TLR2 agonists; PAM3 (TLR2/1 agonist), and PAM2 (TLR2/6 agonist). Additionally, in-silico approaches were performed to confirm the dimer stability using de novo protein structure prediction model for bovine TLRs. The in-vitro approach revealed that sperm triggered the mRNA and protein expression of TLR1 and TLR2 but not TLR6 in BEECs. Moreover, this model disclosed that activation of TLR2/6 heterodimer, triggers a much stronger inflammatory response than TLR2/1 and sperm in bovine uterine epithelia. In the ex-vivo model that mimics the intact uterine tissue at insemination, sperm also induced the protein expression of both TLR1 and TLR2, but not TLR6, in bovine endometrium, particularly in uterine glands. Importantly, PAM3 and sperm induced similar and low mRNA expression of pro-inflammatory cytokines and TNFA protein to a lesser extent than PAM2 in endometrial epithelia. This implied that sperm might trigger a weak inflammatory response via TLR2/TLR1 activation which is similar to that of PAM3. Additionally, the in-silico analyses showed that the existence of bridging ligands is essential for heterodimer stability in bovine TLR2 with either TLR1 or TLR6. Altogether, the present findings revealed that sperm utilize TLR2/1, but not TLR2/6, heterodimerization to trigger a weak physiological inflammatory response in the bovine uterus. This might be the way to remove excess dead sperm remaining in the uterine lumen without tissue damage for providing an ideal uterine environment for early embryo reception and implantation.

Identifying NFKB1, STAT3, and CDKN1A as Baicalein's Potential Hub Targets in Parkinson's Disease-related α-synuclein-mediated Pathways by Integrated Bioinformatics Strategies

BACKGROUND

The overexpression, accumulation, and cell-to-cell transmission of α-synuclein leads to the deterioration of Parkinson's disease (PD). Previous studies suggest that Baicalein (BAI) can bind to α-synuclein and inhibit α-synuclein aggregation and secretion. However, it is still unclear whether BAI can intervene with the pathogenic molecules in α-synuclein-mediated PD pathways beyond directly targeting α-synuclein per se.

METHODS

This study aimed to systematically investigate BAI's potential targets in PD-related A53T mutant α-synuclein-mediated pathways by integrating data mining, network pharmacological analysis, and molecular docking simulation techniques.

RESULTS

The results suggest that BAI may target genes that are dysregulated in synaptic transmission, vesicle trafficking, gene transcription, protein binding, extracellular matrix formation, and kinase activity in α-synucleinmediated pathways. NFKB1, STAT3, and CDKN1A are BAI's potential hub targets in these pathways.

CONCLUSION

Our findings highlight BAI's potentiality to modulate α-synuclein-mediated pathways beyond directly targeting α-synuclein per se.

Exploring the Gallic and Cinnamic Acids Chimeric Derivatives as Anticancer Agents over HeLa Cell Line: An in silico and in vitro Study

Cervical cancer is one of the most aggressive and important cancer types in the female population, due to its low survival rate. Actually, the search for new bioactive compounds, like gallic and cinnamic acid, is one of the most employed options to finding a treatment. In the present study, 134 phenolic compounds with cytotoxic activity over HeLa cell line were used to generate a descriptive ( R 2 ${{R}^{2}}$ =0.76) and predictive ( Q 2 ${{Q}^{2}}$ =0.69 and Q e x t 2 ${{Q}_{{\rm e}{\rm x}{\rm t}}^{2}}$ =0.62) QSAR model. Structural, electronic, steric, and hydrophobic features are represented as different molecular descriptors in our QSAR model. From this model, nine gallate-cinnamate ester derivatives (N1-N9) were designed and synthesized. Furthermore, in vitro cytotoxic activity was evaluated against HeLa and non-tumorigenic cells. Derivatives N6, N5, N1, and N9 were the most active molecules with IC_50ExpHeLa values from 7.26 to 11.95 μM. Finally, the binding of the synthesized compounds to the colchicine binding site on tubulin was evaluated by molecular docking as a possible action mechanism. N1, N5 and N6 can be considered as templates for the design of new cervical anticancer compounds.

Binding of RPR260243 at the intracellular side of the hERG1 channel pore domain slows closure of the helix bundle crossing gate

The opening and closing of voltage-dependent potassium channels is dependent on a tight coupling between movement of the voltage sensing S4 segments and the activation gate. A specific interaction between intracellular amino- and carboxyl-termini is required for the characteristically slow rate of channel closure (deactivation) of hERG1 channels. Compounds that increase hERG1 channel currents represent a novel approach for prevention of arrhythmia associated with prolonged ventricular repolarization. RPR260243 (RPR), a quinoline oxo-propyl piperidine derivative, inhibits inactivation and dramatically slows the rate of hERG1 channel deactivation. Here we report that similar to its effect on wild-type channels, RPR greatly slows the deactivation rate of hERG1 channels missing their amino-termini, or of split channels lacking a covalent link between the voltage sensor domain and the pore domain. By contrast, RPR did not slow deactivation of C-terminal truncated hERG1 channels or D540K hERG1 mutant channels activated by hyperpolarization. Together, these findings indicate that ability of RPR to slow deactivation requires an intact C-terminus, does not slow deactivation by stabilizing an interaction involving the amino-terminus or require a covalent link between the voltage sensor and pore domains. All-atom molecular dynamics simulations using the cryo-EM structure of the hERG1 channel revealed that RPR binds to a pocket located at the intracellular ends of helices S5 and S6 of a single subunit. The slowing of channel deactivation by RPR may be mediated by disruption of normal S5-S6 interactions.

Ligand-Based Drug Design of Novel Antimicrobials against Staphylococcus aureus by Targeting Bacterial Transcription

Staphylococcus aureus is a common human commensal pathogen that causes a wide range of infectious diseases. Due to the generation of antimicrobial resistance, the pathogen becomes resistant to more and more antibiotics, resulting in methicillin-resistant S. aureus (MRSA) and even multidrug-resistant S. aureus (MDRSA), namely 'superbugs'. This situation highlights the urgent need for novel antimicrobials. Bacterial transcription, which is responsible for bacterial RNA synthesis, is a valid but underutilized target for developing antimicrobials. Previously, we reported a novel class of antimicrobials, coined nusbiarylins, that inhibited bacterial transcription by interrupting the protein-protein interaction (PPI) between two transcription factors NusB and NusE. In this work, we developed a ligand-based workflow based on the chemical structures of nusbiarylins and their activity against S. aureus. The ligand-based models-including the pharmacophore model, 3D QSAR, AutoQSAR, and ADME/T calculation-were integrated and used in the following virtual screening of the ChemDiv PPI database. As a result, four compounds, including J098-0498, 1067-0401, M013-0558, and F186-026, were identified as potential antimicrobials against S. aureus, with predicted pMIC values ranging from 3.8 to 4.2. The docking study showed that these molecules bound to NusB tightly with the binding free energy ranging from -58 to -66 kcal/mol.

New putative therapeutic targets against Serratia marcescens using reverse vaccinology and subtractive genomics

The Gram-negative bacillus Serratia marcescens, a member of Enterobacteriaceae family, is an opportunistic nosocomial pathogen commonly found in hospital outbreaks that can cause infections in the urinary tract, bloodstream, central nervous system and pneumonia. Because S. marcescens strains are resistant to several antibiotics, it is critical the need for effective treatments, including new drugs and vaccines. Here, we applied reverse vaccinology and subtractive genomic approaches for the in silico prediction of potential vaccine and drug targets against 59 strains of S. marcescens. We found 759 core non-host homologous proteins, of which 87 are putative surface-exposed proteins, 183 secreted proteins, and 80 membrane proteins. From these proteins, we predicted seven candidates vaccine targets: a sn-glycerol-3-phosphate-binding periplasmic protein UgpB, a vitamin B12 TonB-dependent receptor, a ferrichrome porin FhuA, a divisome-associated lipoprotein YraP, a membrane-bound lytic murein transglycosylase A, a peptidoglycan lytic exotransglycosylase, and a DUF481 domain-containing protein. We also predicted two drug targets: a N(4)-acetylcytidine amidohydrolase, and a DUF1428 family protein. Using the molecular docking approach for each drug target, we identified and selected ZINC04259491 and ZINC04235390 molecules as the most favorable interactions with the target active site residues. Our findings may contribute to the development of vaccines and new drug targets against S. marcescens. Communicated by Ramaswamy H. Sarma.

Comparative genomics of Bordetella pertussis and prediction of new vaccines and drug targets

Pertussis is a highly contagious respiratory disease caused by Bordetella pertussis, a Gram-negative bacterium described over a century ago. Despite broad vaccine coverage and treatment options, the disease is remerging as a public health problem especially in infants and older children. Recent data indicate re-emergence of the disease is related to bacterial resistance to immune defences and decreased vaccine effectiveness, which obviously suggests the need of new effective vaccines and drugs. In an attempt to contribute with solutions to this great challenge, bioinformatics tools were used to genetically comprehend the species of these bacteria and predict new vaccines and drug targets. In fact, approaches were used to analysis genomic plasticity, gene synteny and species similarities between the 20 genomes of Bordetella pertussis already available. Furthermore, it was conducted reverse vaccinology and docking analysis to identify proteins with potential to become vaccine and drug targets, respectively. The analyses showed the 20 genomes belongs to a homogeneous group that has preserved most of the genes over time. Besides that, were found genomics islands and good proteins to be candidates for vaccine and drugs. Taken together, these results suggests new possibilities that may be useful to develop new vaccines and drugs that will help the prevention and treatment strategies of pertussis disease caused by these Bordetella strains. Communicated by Ramaswamy H. Sarma.

JM-20, a Benzodiazepine-Dihydropyridine Hybrid Molecule, Inhibits the Formation of Alpha-Synuclein-Aggregated Species

Studies showed that JM-20, a benzodiazepine-dihydropyridine hybrid molecule, protects against rotenone and 6-hydroxydopamine neurotoxicity. However, its protective effects against cytotoxicity induced by endogenous neurotoxins involved in Parkinson's disease (PD) pathogenesis have never been investigated. In this study, we evaluated the ability of JM-20 to inhibit alpha-synuclein (aSyn) aggregation. We also evaluated the interactions of JM-20 with aSyn by molecular docking and molecular dynamics and assessed the protective effect of JM-20 against aminochrome cytotoxicity. We demonstrated that JM-20 induced the formation of heterogeneous amyloid fibrils, which were innocuous to primary cultures of mesencephalic cells. Moreover, JM-20 reduced the average size of aSyn positive inclusions in H4 cells transfected with SynT wild-type and synphilin-1-V5, but not in HEK cells transfected with synphilin-1-GFP. In silico studies showed the interaction between JM-20 and the aSyn-binding site. Additionally, we showed that JM-20 protects SH-SY5Y cells against aminochrome cytotoxicity. These results reinforce the potential of JM-20 as a neuroprotective compound for PD and suggest aSyn as a molecular target for JM-20.

A Scaffold-based Deep Generative Model Considering Molecular Stereochemical Information

Designing molecules with specific scaffolds can facilitate the discovery and optimization of lead compounds. Some scaffold-based molecular generation models have been developed using deep-learning methods based on specific scaffolds, although incorporating scaffold generalization is expected to achieve scaffold hopping. Moreover, most of the existing models focus on the 2D shape of the scaffold and overlook the stereochemical properties of the compound, especially for natural products. In this study, we optimized the scaffold-based molecular generation model designed by Lim et al. (Chemical Science 2020, 11, 1153-1164). Real-time ultrafast shape recognition with pharmacophore constraints (USRCAT) was introduced into the model to search for molecules similar to the 3D conformation and pharmacophore of the input scaffold sourced from the training set; the searched molecules were then used as new scaffolds to execute scaffold hopping. The optimized model could generate new molecules with the same chirality as the input scaffold. Furthermore, the probability distribution of the molecular structure and various physicochemical properties were analyzed to evaluate the model's generation capability. We thus believe that the optimized model can provide a basis for medicinal chemists to explore a wider chemical space toward optimization of the lead compounds and to screen the virtual compound library.