MHC2AffyPred: A machine-learning approach to estimate affinity of MHC class II peptides based on structural interaction fingerprints.

Understanding how MHC class II (MHC-II) binding peptides with differing lengths exhibit specific interaction at the core and extended sites within the large MHC-II pocket is a very important aspect of immunological research for designing peptides. Certain efforts were made to generate peptide conformations amenable for MHC-II binding and calculate the binding energy of such complex formation but not directed toward developing a relationship between the peptide conformation in MHC-II structures and the binding affinity (BA) (IC50 ). We present here a machine-learning approach to calculate the BA of the peptides within the MHC-II pocket for HLA-DRA1, HLA-DRB1, HLA-DP, and HLA-DQ allotypes. Instead of generating ensembles of peptide conformations conventionally, the biased mode of conformations was created by considering the peptides in the crystal structures of pMHC-II complexes as the templates, followed by site-directed peptide docking. The structural interaction fingerprints generated from such docked pMHC-II structures along with the Moran autocorrelation descriptors were trained using a random forest regressor specific to each MHC-II peptide lengths (9-19). The entire workflow is automated using Linux shell and Perl scripts to promote the utilization of MHC2AffyPred program to any characterized MHC-II allotypes and is made for free access at https://github.com/SiddhiJani/MHC2AffyPred. The MHC2AffyPred attained better performance (correlation coefficient [CC] of .612-.898) than MHCII3D (.03-.594) and NetMHCIIpan-3.2 (.289-.692) programs in the HLA-DRA1, HLA-DRB1 types. Similarly, the MHC2AffyPred program achieved CC between .91 and .98 for HLA-DP and HLA-DQ peptides (13-mer to 17-mer). Further, a case study on MHC-II binding 15-mer peptides of severe acute respiratory syndrome coronavirus-2 showed very close competency in computing the IC50 values compared to the sequence-based NetMHCIIpan v3.2 and v4.0 programs with a correlation of .998 and .570, respectively.

Pancancer network analysis reveals key master regulators for cancer invasiveness.

BACKGROUND: Tumor invasiveness reflects numerous biological changes, including tumorigenesis, progression, and metastasis. To decipher the role of transcriptional regulators (TR) involved in tumor invasiveness, we performed a systematic network-based pan-cancer assessment of master regulators of cancer invasiveness. MATERIALS AND METHODS: We stratified patients in The Cancer Genome Atlas (TCGA) into invasiveness high (INV-H) and low (INV-L) groups using consensus clustering based on an established robust 24-gene signature to determine the prognostic association of invasiveness with overall survival (OS) across 32 different cancers. We devise a network-based protocol to identify TRs as master regulators (MRs) unique to INV-H and INV-L phenotypes. We validated the activity of MRs coherently associated with INV-H phenotype and worse OS across cancers in TCGA on a series of additional datasets in the Prediction of Clinical Outcomes from the Genomic Profiles (PRECOG) repository. RESULTS: Based on the 24-gene signature, we defined the invasiveness score for each patient sample and stratified patients into INV-H and INV-L clusters. We observed that invasiveness was associated with worse survival outcomes in almost all cancers and had a significant association with OS in ten out of 32 cancers. Our network-based framework identified common invasiveness-associated MRs specific to INV-H and INV-L groups across the ten prognostic cancers, including COL1A1, which is also part of the 24-gene signature, thus acting as a positive control. Downstream pathway analysis of MRs specific to INV-H phenotype resulted in the identification of several enriched pathways, including Epithelial into Mesenchymal Transition, TGF-ß signaling pathway, regulation of Toll-like receptors, cytokines, and inflammatory response, and selective expression of chemokine receptors during T-cell polarization. Most of these pathways have connotations of inflammatory immune response and feasibility for metastasis. CONCLUSION: Our pan-cancer study provides a comprehensive master regulator analysis of tumor invasiveness and can suggest more precise therapeutic strategies by targeting the identified MRs and downstream enriched pathways for patients across multiple cancers.

Excavating phytochemicals from plants possessing antiviral activities for identifying SARS-CoV hemagglutinin-esterase inhibitors by diligent computational workflow.

Coronaviruses (CoVs) belong to a group of RNA viruses that cause diseases in vertebrates including. Newer and deadlier than SARS CoV-2 are sought to appear in future for which the scientific community must be prepared with the strategies for their control. Spike protein (S-protein) of all the CoVs require angiotensin-converting enzyme2 (ACE2), while CoVs also require hemagglutinin-acetylesterase (HE) glycoprotein receptor to simultaneously interact with O-acetylated sialic acids on host cells, both these interactions enable viral particle to enter host cell leading to its infection. Target inhibition of viral S-protein and HE glycoprotein receptor can lead to a development of therapy against the SARS CoV-2. The proposition is to recognize molecules from the bundle of phytochemicals of medicinal plants known to possess antiviral potentials as a lead that could interact and mask the active site of, HE glycoprotein which would ideally bind to O-acetylated sialic acids on human host cells. Such molecules can be addressed as 'HE glycoprotein blockers'. A library of 110 phytochemicals from Withania somnifera, Asparagus racemosus, Zinziber officinalis, Allium sativum, Curcuma longa and Adhatoda vasica was constructed and was used under present study. In silico analysis was employed with plant-derived phytochemicals. The molecular docking, molecular dynamics simulations over the scale of 1000 ns (1 µs) and ADMET prediction revealed that the Withania somnifera (ashwagandha) and Asparagus racemosus (shatavari) plants possessed various steroidal saponins and alkaloids which could potentially inhibit the COVID-19 virus and even other CoVs targeted HE glycoprotein receptor.Communicated by Ramaswamy H. Sarma.

Computational investigation of natural compounds as potential main protease (Mpro) inhibitors for SARS-CoV-2 virus.

The coronavirus disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is significantly impacting human lives, overburdening the healthcare system and weakening global economies. Plant-derived natural compounds are being largely tested for their efficacy against COVID-19 targets to combat SARS-CoV-2 infection. The SARS-CoV-2 Main protease (Mpro) is considered an appealing target because of its role in replication in host cells. We curated a set of 7809 natural compounds by combining the collections of five databases viz Dr Duke's Phytochemical and Ethnobotanical database, IMPPAT, PhytoHub, AromaDb and Zinc. We applied a rigorous computational approach to identify lead molecules from our curated compound set using docking, dynamic simulations, the free energy of binding and DFT calculations. Theaflavin and ginkgetin have emerged as better molecules with a similar inhibition profile in both SARS-CoV-2 and Omicron variants.

Identification of antiviral phytochemicals as a potential SARS-CoV-2 main protease (Mpro) inhibitor using docking and molecular dynamics simulations.

Novel SARS-CoV-2, an etiological factor of Coronavirus disease 2019 (COVID-19), poses a great challenge to the public health care system. Among other druggable targets of SARS-Cov-2, the main protease (Mpro) is regarded as a prominent enzyme target for drug developments owing to its crucial role in virus replication and transcription. We pursued a computational investigation to identify Mpro inhibitors from a compiled library of natural compounds with proven antiviral activities using a hierarchical workflow of molecular docking, ADMET assessment, dynamic simulations and binding free-energy calculations. Five natural compounds, Withanosides V and VI, Racemosides A and B, and Shatavarin IX, obtained better binding affinity and attained stable interactions with Mpro key pocket residues. These intermolecular key interactions were also retained profoundly in the simulation trajectory of 100 ns time scale indicating tight receptor binding. Free energy calculations prioritized Withanosides V and VI as the top candidates that can act as effective SARS-CoV-2 Mpro inhibitors.