Towards understanding aromatase inhibitory activity via QSAR modeling

First report on chemometrics-driven multilayered lead prioritization in addressing oxysterol-mediated overexpression of G protein-coupled receptor 183

Contemporary research has convincingly demonstrated that upregulation of G protein-coupled receptor 183 (GPR183), orchestrated by its endogenous agonist, 7α,25-dihydroxyxcholesterol (7α,25-OHC), leads to the development of cancer, diabetes, multiple sclerosis, infectious, and inflammatory diseases. A recent study unveiled the cryo-EM structure of 7α,25-OHC bound GPR183 complex, presenting an untapped opportunity for computational exploration of potential GPR183 inhibitors, which served as our inspiration for the current work. A predictive and validated two-dimensional QSAR model using genetic algorithm (GA) and multiple linear regression (MLR) on experimental GPR183 inhibition data was developed. QSAR study highlighted that structural features like dissimilar electronegative atoms, quaternary carbon atoms, and CH2RX fragment (X: heteroatoms) influence positively, while the existence of oxygen atoms with a topological separation of 3, negatively affects GPR183 inhibitory activity. Post assessment of true external set prediction capability, the MLR model was deployed to screen 12,449 DrugBank compounds, followed by a screening pipeline involving molecular docking, druglikeness, ADMET, protein-ligand stability assessment using deep learning algorithm, molecular dynamics, and molecular mechanics. The current findings strongly evidenced DB05790 as a potential lead for prospective interference of oxysterol-mediated GPR183 overexpression, warranting further in vitro and in vivo validation.

Classification models for predicting the bioactivity of pan-TRK inhibitors and SAR analysis

Tropomyosin receptor kinases (TRKs) are important broad-spectrum anticancer targets. The oncogenic rearrangement of the NTRK gene disrupts the extracellular structural domain and epitopes for therapeutic antibodies, making small-molecule inhibitors essential for treating NTRK fusion-driven tumors. In this work, several algorithms were used to construct descriptor-based and nondescriptor-based models, and the models were evaluated by outer 10-fold cross-validation. To find a model with good generalization ability, the dataset was partitioned by random and cluster-splitting methods to construct in- and cross-domain models, respectively. Among the 48 models built, the model with the combination of the deep neural network (DNN) algorithm and extended connectivity fingerprints 4 (ECFP4) descriptors achieved excellent performance in both dataset divisions. The results indicate that the DNN algorithm has a strong generalization prediction ability, and the richness of features plays a vital role in predicting unknown spatial molecules. Additionally, we combined the clustering results and decision tree models of fingerprint descriptors to perform structure-activity relationship analysis. It was found that nitrogen-containing aromatic heterocyclic and benzo heterocyclic structures play a crucial role in enhancing the activity of TRK inhibitors.

Mpropred: A machine learning (ML) driven Web-App for bioactivity prediction of SARS-CoV-2 main protease (Mpro) antagonists

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic emerged in 2019 and still requiring treatments with fast clinical translatability. Frequent occurrence of mutations in spike glycoprotein of SARS-CoV-2 led the consideration of an alternative therapeutic target to combat the ongoing pandemic. The main protease (Mpro) is such an attractive drug target due to its importance in maturating several polyproteins during the replication process. In the present study, we used a classification structure-activity relationship (CSAR) model to find substructures that leads to to anti-Mpro activities among 758 non-redundant compounds. A set of 12 fingerprints were used to describe Mpro inhibitors, and the random forest approach was used to build prediction models from 100 distinct data splits. The data set's modelability (MODI index) was found to be robust, with a value of 0.79 above the 0.65 threshold. The accuracy (89%), sensitivity (89%), specificity (73%), and Matthews correlation coefficient (79%) used to calculate the prediction performance, was also found to be statistically robust. An extensive analysis of the top significant descriptors unveiled the significance of methyl side chains, aromatic ring and halogen groups for Mpro inhibition. Finally, the predictive model is made publicly accessible as a web-app named Mpropred in order to allow users to predict the bioactivity of compounds against SARS-CoV-2 Mpro. Later, CMNPD, a marine compound database was screened by our app to predict bioactivity of all the compounds and results revealed significant correlation with their binding affinity to Mpro. Molecular dynamics (MD) simulation and molecular mechanics/Poisson Boltzmann surface area (MM/PBSA) analysis showed improved properties of the complexes. Thus, the knowledge and web-app shown herein can be used to develop more effective and specific inhibitors against the SARS-CoV-2 Mpro. The web-app can be accessed from https://share.streamlit.io/nadimfrds/mpropred/Mpropred_app.py.

Towards combating antibiotic resistance by exploring the quantitative structure-activity relationship of NDM-1 inhibitors

The emergence of New Delhi metallo-beta-lactamase-1 (NDM-1) has conferred enteric bacteria resistance to almost all beta-lactam antibiotics. Its capability of horizontal transfer through plasmids, amongst humans, animal reservoirs and the environment, has added up to the totality of antimicrobial resistance control, animal husbandry and food safety. Thus far, there have been no effective drugs for neutralizing NDM-1. This study explores the structure-activity relationship of NDM-1 inhibitors. IC₅₀ values of NDM-1 inhibitors were compiled from both the ChEMBL database and literature. After curation, a final set of 686 inhibitors were used for machine learning model building using the random forest algorithm against 12 sets of molecular fingerprints. Benchmark results indicated that the KlekotaRothCount fingerprint provided the best overall performance with an accuracy of 0.978 and 0.778 for the training and testing set, respectively. Model interpretation revealed that nitrogen-containing features (KRFPC 4080, KRFPC 3882, KRFPC 677, KRFPC 3608, KRFPC 3750, KRFPC 4287 and KRFPC 3943), sulfur-containing substructures (KRFPC 2855 and KRFPC 4843), aromatic features (KRFPC 1566, KRFPC 1564, KRFPC 1642, KRFPC 3608, KRFPC 4287 and KRFPC 3943), carbonyl features (KRFPC 1193 and KRFPC 3025), aliphatic features (KRFPC 2975, KRFPC 297, KRFPC 3224 and KRFPC 669) are features contributing to NDM-1 inhibitory activity. It is anticipated that findings from this study would help facilitate the drug discovery of NDM-1 inhibitors by providing guidelines for further lead optimization.

SCORPION is a stacking-based ensemble learning framework for accurate prediction of phage virion proteins

Fast and accurate identification of phage virion proteins (PVPs) would greatly aid facilitation of antibacterial drug discovery and development. Although, several research efforts based on machine learning (ML) methods have been made for in silico identification of PVPs, these methods have certain limitations. Therefore, in this study, we propose a new computational approach, termed SCORPION, (StaCking-based Predictior fOR Phage VIrion PrOteiNs), to accurately identify PVPs using only protein primary sequences. Specifically, we explored comprehensive 13 different feature descriptors from different aspects (i.e., compositional information, composition-transition-distribution information, position-specific information and physicochemical properties) with 10 popular ML algorithms to construct a pool of optimal baseline models. These optimal baseline models were then used to generate probabilistic features (PFs) and considered as a new feature vector. Finally, we utilized a two-step feature selection strategy to determine the optimal PF feature vector and used this feature vector to develop a stacked model (SCORPION). Both tenfold cross-validation and independent test results indicate that SCORPION achieves superior predictive performance than its constitute baseline models and existing methods. We anticipate SCORPION will serve as a useful tool for the cost-effective and large-scale screening of new PVPs. The source codes and datasets for this work are available for downloading in the GitHub repository (https://github.com/saeed344/SCORPION).

Ensemble machine learning to evaluate the in vivo acute oral toxicity and in vitro human acetylcholinesterase inhibitory activity of organophosphates

Organophosphates (OPs) are hazardous chemicals widely used in industry and agriculture. Distribution of their residues in nature causes serious risks to humans, animals, and plants. To reduce hazards from OPs, quantitative structure-activity relationship (QSAR) models for predicting their acute oral toxicity in rats and mice and inhibition constants concerning human acetylcholinesterase were developed according to the bioactivity data of 456 unique OPs. Based on robust, two-dimensional molecular descriptors and quantum chemical descriptors, which accurately reflect OP electronic structures and reactivities, the influences of eight machine-learning algorithms on the prediction performance of the QSAR models were explored, and consensus QSAR models were constructed. Several strict model validation indices and the results of applicability domain evaluations show that the established consensus QSAR models exhibit good robustness, practical prediction abilities, and wide application scopes. Poor correlation was observed between acute oral toxicity at the mammalian level and the inhibition constants at the molecular level, indicating that the acute toxicity of OPs cannot be evaluated only by the experimental data of enzyme inhibitory activity, their toxicokinetic characteristics must also be considered. The constructed QSAR models described herein provide rapid, theoretical assessment of the bioactivity of unstudied or unknown OPs, as well as guidance for making decisions regarding their regulation.

The Prognostic Value of Combination of Plasma Fibrinogen and CA19-9 in Non-Distant Metastatic Breast Cancer Patients Undergoing Surgery

Purpose

This article aimed to study the prognostic value of preoperative plasma fibrinogen and CA19-9 in non-distant metastatic breast cancer (BC).

Patients and Methods

A total of 343 non-distant metastatic BC patients were included in this study. The optimal cut-off values of plasma fibrinogen and CA19-9 were obtained by receiver operating characteristic (ROC) curve analysis. Univariate and multivariate Cox regression analyses were used to evaluate prognostic factors for overall survival (OS). Survival data were assessed using Kaplan–Meier survival analysis with the Log-rank test. Based on the cut-off values, we classified the fibrinogen-CA19-9 score as follows: 2 (both hyperfibrinogenemia and high CA19-9), 1 (either hyperfibrinogenemia or high CA19-9), and 0 (neither hypefibrinogenemia nor high CA19-9).

Results

Our follow-up time totaled 10 years, the median follow-up time was 77 months (range=2–119 months), and 82 (23.9%) of 343 patients died during the follow-up period. The optimal cut-off values of plasma fibrinogen and CA19-9 were 2.805 g/L and 11.85 U/mL, respectively. The multivariate Cox analysis results suggested that there was a significant association between worse OS and elevated preoperative plasma fibrinogen and CA19-9 levels (HR=2.016, 95% CI=1.216–3.342, P=0.007; and HR=2.042, 95% CI=1.282–3.253, P=0.003). The area under the ROC curve (AUC) increased from 0.589 (for plasma fibrinogen) and 0.594 (for CA19-9) to 0.640 when these two parameters were combined. When we added this combined factor to the multivariate analysis, it was an independent prognostic factor for BC (P<0.001). According to the above results, we chose four prognostic factors to construct our nomogram. The AUC was 0.724, which indicates that the nomogram performs well.

Conclusion

The combination of plasma fibrinogen and CA19-9 could be used as a valid independent prognostic factor for non-distant metastatic BC compared with either parameter alone and could easily be applied in clinical practice.

Synthesis, Docking Studies and Biological Activity of New Benzimidazole- Triazolothiadiazine Derivatives as Aromatase Inhibitor

In the last step of estrogen biosynthesis, aromatase enzyme catalyzes the conversion of androgens to estrogens. Aromatase inhibition is an important way to control estrogen-related diseases and estrogen levels. In this study, sixteen of benzimidazole-triazolothiadiazine derivatives have been synthesized and studied as potent aromatase inhibitors. First, these compounds were tested for their anti-cancer properties against human breast cancer cell line (MCF-7). The most active compounds 5c, 5e, 5k, and 5m on MCF-7 cell line were subject to further in vitro aromatase enzyme inhibition assays to determine the possible mechanisms of action underlying their activity. Compound 5e showed slight less potent aromatase inhibitory activity than that of letrozole with IC₅₀ = 0.032 ± 0.042 µM, compared to IC₅₀ = 0.024 ± 0.001 µM for letrozole. Furthermore, compound 5e and reference drug letrozole were docked into human placental aromatase enzyme to predict their possible binding modes with the enzyme. Finally, ADME parameters (absorption, distribution, metabolism, and excretion) of synthesized compounds (5a–5p) were calculated by QikProp 4.8 software.

Towards an Understanding of the Mode of Action of Human Aromatase Activity for Azoles through Quantum Chemical Descriptors-Based Regression and Structure Activity Relationship Modeling Analysis

Aromatase is an enzyme member of the cytochrome P450 superfamily coded by the CYP19A1 gene. Its main action is the conversion of androgens into estrogens, transforming androstenedione into estrone and testosterone into estradiol. This enzyme is present in several tissues and it has a key role in the maintenance of the balance of androgens and estrogens, and therefore in the regulation of the endocrine system. With regard to chemical safety and human health, azoles, which are used as agrochemicals and pharmaceuticals, are potential endocrine disruptors due to their agonist or antagonist interactions with the human aromatase enzyme. This theoretical study investigated the active agonist and antagonist properties of “chemical classes of azoles” to determine the relationships of azole interaction with CYP19A1, using stereochemical and electronic properties of the molecules through classification and multilinear regression (MLR) modeling. The antagonist activities for the same substituent on diazoles and triazoles vary with its chemical composition and its position and both heterocyclic systems require aromatic substituents. The triazoles require the spherical shape and diazoles have to be in proper proportion of the branching index and the number of ring systems for the inhibition. Considering the electronic aspects, triazole antagonist activity depends on the electrophilicity index that originates from interelectronic exchange interaction (ω_HF) and the LUMO energy (ELUMOPM7), and the diazole antagonist activity originates from the penultimate orbital (EHOMONLPM7) of diazoles. The regression models for agonist activity show that it is opposed by the static charges but favored by the delocalized charges on the diazoles and thiazoles. This study proposes that the electron penetration of azoles toward heme group decides the binding behavior and stereochemistry requirement for antagonist activity against CYP19A1 enzyme.

PVPred-SCM: Improved Prediction and Analysis of Phage Virion Proteins Using a Scoring Card Method

Although, existing methods have been successful in predicting phage (or bacteriophage) virion proteins (PVPs) using various types of protein features and complex classifiers, such as support vector machine and naïve Bayes, these two methods do not allow interpretability. However, the characterization and analysis of PVPs might be of great significance to understanding the molecular mechanisms of bacteriophage genetics and the development of antibacterial drugs. Hence, we herein proposed a novel method (PVPred-SCM) based on the scoring card method (SCM) in conjunction with dipeptide composition to identify and characterize PVPs. In PVPred-SCM, the propensity scores of 400 dipeptides were calculated using the statistical discrimination approach. Rigorous independent validation test showed that PVPred-SCM utilizing only dipeptide composition yielded an accuracy of 77.56%, indicating that PVPred-SCM performed well relative to the state-of-the-art method utilizing a number of protein features. Furthermore, the propensity scores of dipeptides were used to provide insights into the biochemical and biophysical properties of PVPs. Upon comparison, it was found that PVPred-SCM was superior to the existing methods considering its simplicity, interpretability, and implementation. Finally, in an effort to facilitate high-throughput prediction of PVPs, we provided a user-friendly web-server for identifying the likelihood of whether or not these sequences are PVPs. It is anticipated that PVPred-SCM will become a useful tool or at least a complementary existing method for predicting and analyzing PVPs.

iQSP: A Sequence-Based Tool for the Prediction and Analysis of Quorum Sensing Peptides via Chou's 5-Steps Rule and Informative Physicochemical Properties

Understanding of quorum-sensing peptides (QSPs) in their functional mechanism plays an essential role in finding new opportunities to combat bacterial infections by designing drugs. With the avalanche of the newly available peptide sequences in the post-genomic age, it is highly desirable to develop a computational model for efficient, rapid and high-throughput QSP identification purely based on the peptide sequence information alone. Although, few methods have been developed for predicting QSPs, their prediction accuracy and interpretability still requires further improvements. Thus, in this work, we proposed an accurate sequence-based predictor (called iQSP) and a set of interpretable rules (called IR-QSP) for predicting and analyzing QSPs. In iQSP, we utilized a powerful support vector machine (SVM) cooperating with 18 informative features from physicochemical properties (PCPs). Rigorous independent validation test showed that iQSP achieved maximum accuracy and MCC of 93.00% and 0.86, respectively. Furthermore, a set of interpretable rules IR-QSP was extracted by using random forest model and the 18 informative PCPs. Finally, for the convenience of experimental scientists, the iQSP web server was established and made freely available online. It is anticipated that iQSP will become a useful tool or at least as a complementary existing method for predicting and analyzing QSPs.

Meta-iAVP: A Sequence-Based Meta-Predictor for Improving the Prediction of Antiviral Peptides Using Effective Feature Representation

In spite of the large-scale production and widespread distribution of vaccines and antiviral drugs, viruses remain a prominent human disease. Recently, the discovery of antiviral peptides (AVPs) has become an influential antiviral agent due to their extraordinary advantages. With the avalanche of newly-found peptide sequences in the post-genomic era, there is a great demand to develop a sequence-based predictor for timely identifying AVPs as this information is very useful for both basic research and drug development. In this study, we propose a novel sequence-based meta-predictor with an effective feature representation, called Meta-iAVP, for the accurate prediction of AVPs from given peptide sequences. Herein, the effective feature representation was extracted from a set of prediction scores derived from various machine learning algorithms and types of features. To the best of our knowledge, the model proposed herein represents the first meta-based approach for the prediction of AVPs. An overall accuracy and Matthews correlation coefficient of 95.20% and 0.90, respectively, was achieved from the independent test set on an objective benchmark dataset. Comparative analysis suggested that Meta-iAVP was superior to that of existing methods and therefore represents a useful tool for AVP prediction. Finally, in an effort to facilitate high-throughput prediction of AVPs, the model was deployed as the Meta-iAVP web server and is made freely available online at http://codes.bio/meta-iavp/ where users can submit query peptide sequences for determining the likelihood of whether or not these peptides are AVPs.

ACPred: A Computational Tool for the Prediction and Analysis of Anticancer Peptides

Anticancer peptides (ACPs) have emerged as a new class of therapeutic agent for cancer treatment due to their lower toxicity as well as greater efficacy, selectivity and specificity when compared to conventional small molecule drugs. However, the experimental identification of ACPs still remains a time-consuming and expensive endeavor. Therefore, it is desirable to develop and improve upon existing computational models for predicting and characterizing ACPs. In this study, we present a bioinformatics tool called the ACPred, which is an interpretable tool for the prediction and characterization of the anticancer activities of peptides. ACPred was developed by utilizing powerful machine learning models (support vector machine and random forest) and various classes of peptide features. It was observed by a jackknife cross-validation test that ACPred can achieve an overall accuracy of 95.61% in identifying ACPs. In addition, analysis revealed the following distinguishing characteristics that ACPs possess: (i) hydrophobic residue enhances the cationic properties of α-helical ACPs resulting in better cell penetration; (ii) the amphipathic nature of the α-helical structure plays a crucial role in its mechanism of cytotoxicity; and (iii) the formation of disulfide bridges on β-sheets is vital for structural maintenance which correlates with its ability to kill cancer cells. Finally, for the convenience of experimental scientists, the ACPred web server was established and made freely available online.