Identification of RdRp inhibitors against SARS-CoV-2 through E-pharmacophore-based virtual screening, molecular docking and MD simulations approaches.

The outbreak of novel Coronavirus, an enduring pandemic declared by WHO, has consequences to an alarming ongoing public health menace which has already claimed several million human lives. In addition to numerous vaccinations and medications for mild to moderate COVID-19 infection, lack of promising medication or therapeutic pharmaceuticals remains a serious concern to counter the ongoing coronavirus infections and to hinder its dreadful spread. Global health emergencies have called for urgency for potential drug discovery and time is the biggest constraint apart from the financial and human resources required for the high throughput drug screening. However, computational screening or in-silico approaches appeared to be an effective and faster approach to discover potential molecules without sacrificing the model animals. Accumulated shreds of evidence on computational studies against viral diseases have revealed significance of in-silico drug discovery approaches especially in the time of urgency. The central role of RdRp in SARS-CoV-2 replication makes it promising drug target to curtain on going infection and its spread. The present study aimed to employ E-pharmacophore-based virtual screening to reveal potent inhibitors of RdRp as potential leads to block the viral replication. An energy-optimised pharmacophore model was generated to screen the Enamine REAL DataBase (RDB). Then, ADME/T profiles were determined to validate the pharmacokinetics and pharmacodynamics properties of the hit compounds. Moreover, High Throughput Virtual Screening (HTVS) and molecular docking (SP & XP) were employed to screen the top hits from pharmacophore-based virtual screening and ADME/T screen. The binding free energies of the top hits were calculated by conducting MM-GBSA analysis followed by MD simulations to determine the stability of molecular interactions between top hits and RdRp protein. These virtual investigations revealed six compounds having binding free energies of -57.498, -45.776, -46.248, -35.67, -25.15 and -24.90 kcal/mol respectively as calculated by the MM-GBSA method. The MD simulation studies confirmed the stability of protein ligand complexes, hence, indicating as potent RdRp inhibitors and are promising candidate drugs to be further validated and translated into clinics in future.

Identification of NS2B-NS3 Protease Inhibitors for Therapeutic Application in ZIKV Infection: A Pharmacophore-Based High-Throughput Virtual Screening and MD Simulations Approaches.

Zika virus (ZIKV) pandemic and its implication in congenital malformations and severe neurological disorders had created serious threats to global health. ZIKV is a mosquito-borne flavivirus which spread rapidly and infect a large number of people in a shorter time-span. Due to the lack of effective therapeutics, this had become paramount urgency to discover effective drug molecules to encounter the viral infection. Various anti-ZIKV drug discovery efforts during the past several years had been unsuccessful to develop an effective cure. The NS2B-NS3 protein was reported as an attractive therapeutic target for inhibiting viral proliferation, due to its central role in viral replication and maturation of non-structural viral proteins. Therefore, the current in silico drug exploration aimed to identify the novel inhibitors of Zika NS2B-NS3 protease by implementing an e-pharmacophore-based high-throughput virtual screening. A 3D e-pharmacophore model was generated based on the five-featured (ADPRR) pharmacophore hypothesis. Subsequently, the predicted model is further subjected to the high-throughput virtual screening to reveal top hit molecules from the various small molecule databases. Initial hits were examined in terms of binding free energies and ADME properties to identify the candidate hit exhibiting a favourable pharmacokinetic profile. Eventually, molecular dynamic (MD) simulations studies were conducted to evaluate the binding stability of the hit molecule inside the receptor cavity. The findings of the in silico analysis manifested affirmative evidence for three hit molecules with -64.28, -55.15 and -50.16 kcal/mol binding free energies, as potent inhibitors of Zika NS2B-NS3 protease. Hence, these molecules holds the promising potential to serve as a prospective candidates to design effective drugs against ZIKV and related viral infections.

Correction: Computational investigation of the increased virulence and pathogenesis of SARS-CoV-2 lineage B.1.1.7.

Correction for 'Computational investigation of the increased virulence and pathogenesis of SARS-CoV-2 lineage B.1.1.7' by N. Arul Murugan et al., Phys. Chem. Chem. Phys., 2022, 24, 20371-20380, https://doi.org/10.1039/D2CP00469K.

Computational investigation of the increased virulence and pathogenesis of SARS-CoV-2 lineage B.1.1.7.

New variants of SARS-CoV-2 are being reported worldwide. The World Health Organization has reported Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2) and Omicron (B.1.1.529) as the variants of concern. There are speculations that the variants might evade the host immune responses induced by currently available vaccines and develop resistance to drugs under consideration. The first step of viral infection in COVID-19 occurs through the interaction of the spike protein's receptor-binding domain (RBD) with the peptidase domain of the human ACE-2 (hACE-2) receptor. This study aims to get a molecular-level understanding of the mechanism behind the increased infection rate in the alpha variant. We have computationally studied the spike protein interaction in both the wild-type and B.1.1.7 variant with the hACE-2 receptor using molecular dynamics and MM-GBSA based binding free energy calculations. The binding free energy difference shows that the mutant variant of the spike protein has increased binding affinity for the hACE-2 receptor (i.e. ΔG(N501Y,A570D) is in the range -7.2 to -7.6 kcal mol-1) and the results were validated using Density functional theory. We demonstrate that with the use of state-of-the-art computational approaches, we can, in advance, predict the virulent nature of variants of SARS-CoV-2 and alert the world healthcare system.

Common and mutation specific phenotypes of KRAS and BRAF mutations in colorectal cancer cells revealed by integrative -omics analysis.

BACKGROUND: Genes in the Ras pathway have somatic mutations in at least 60 % of colorectal cancers. Despite activating the same pathway, the BRAF V600E mutation and the prevalent mutations in codon 12 and 13 of KRAS have all been linked to different clinical outcomes, but the molecular mechanisms behind these differences largely remain to be clarified. METHODS: To characterize the similarities and differences between common activating KRAS mutations and between KRAS and BRAF mutations, we used genome editing to engineer KRAS G12C/D/V and G13D mutations in colorectal cancer cells that had their mutant BRAF V600E allele removed and subjected them to transcriptome sequencing, global proteomics and metabolomics analyses. RESULTS: By intersecting differentially expressed genes, proteins and metabolites, we uncovered (i) two-fold more regulated genes and proteins when comparing KRAS to BRAF mutant cells to those lacking Ras pathway mutation, (ii) five differentially expressed proteins in KRAS mutants compared to cells lacking Ras pathway mutation (IFI16, S100A10, CD44, GLRX and AHNAK2) and 6 (CRABP2, FLNA, NXN, LCP1, S100A10 and S100A2) compared to BRAF mutant cells, (iii) 19 proteins expressed differentially in a KRAS mutation specific manner versus BRAF V600E cells, (iv) regulation of the Integrin Linked Kinase pathway by KRAS but not BRAF mutation, (v) regulation of amino acid metabolism, particularly of the tyrosine, histidine, arginine and proline pathways, the urea cycle and purine metabolism by Ras pathway mutations, (vi) increased free carnitine in KRAS and BRAF mutant RKO cells. CONCLUSIONS: This comprehensive integrative -omics analysis confirms known and adds novel genes, proteins and metabolic pathways regulated by mutant KRAS and BRAF signaling in colorectal cancer. The results from the new model systems presented here can inform future development of diagnostic and therapeutic approaches targeting tumors with KRAS and BRAF mutations.

Monoallelic and bi-allelic variants in NCDN cause neurodevelopmental delay, intellectual disability, and epilepsy.

Neurochondrin (NCDN) is a cytoplasmatic neural protein of importance for neural growth, glutamate receptor (mGluR) signaling, and synaptic plasticity. Conditional loss of Ncdn in mice neural tissue causes depressive-like behaviors, impaired spatial learning, and epileptic seizures. We report on NCDN missense variants in six affected individuals with variable degrees of developmental delay, intellectual disability (ID), and seizures. Three siblings were found homozygous for a NCDN missense variant, whereas another three unrelated individuals carried different de novo missense variants in NCDN. We assayed the missense variants for their capability to rescue impaired neurite formation in human neuroblastoma (SH-SY5Y) cells depleted of NCDN. Overexpression of wild-type NCDN rescued the neurite-phenotype in contrast to expression of NCDN containing the variants of affected individuals. Two missense variants, associated with severe neurodevelopmental features and epilepsy, were unable to restore mGluR5-induced ERK phosphorylation. Electrophysiological analysis of SH-SY5Y cells depleted of NCDN exhibited altered membrane potential and impaired action potentials at repolarization, suggesting NCDN to be required for normal biophysical properties. Using available transcriptome data from human fetal cortex, we show that NCDN is highly expressed in maturing excitatory neurons. In combination, our data provide evidence that bi-allelic and de novo variants in NCDN cause a clinically variable form of neurodevelopmental delay and epilepsy, highlighting a critical role for NCDN in human brain development.

The Therapeutic Prospects of Naturally Occurring and Synthetic Indole Alkaloids for Depression and Anxiety Disorders.

Depression and anxiety are the most common disorders among all age groups. Several antidepressant drugs including benzodiazepine, antidepressant tricyclics, azapirone, noradrenaline reuptake inhibitors, serotonin selective reuptake inhibitors, serotonin, noradrenaline reuptake inhibitors, and monoamine oxidase inhibitors have been used to treat these psychiatric disorders. However, these antidepressants are generally synthetic agents and can cause a wide range of side effects. The potential efficacy of plant-derived alkaloids has been reviewed against various neurodegenerative diseases including Alzheimer's disease, Huntington disease, Parkinson's disease, schizophrenia, and epilepsy. However, data correlating the indole alkaloids and antidepressant activity are limited. Natural products, especially plants and the marine environment, are rich sources of potential new drugs. Plants possess a variety of indole alkaloids, and compounds that have an indole moiety are related to serotonin, which is a neurotransmitter that regulates brain function and cognition, which in turn alleviates anxiety, and ensures a good mood and happiness. The present review is a summary of the bioactive compounds from plants and marine sources that contain the indole moiety, which can serve as potent antidepressants. The prospects of naturally occurring as well as synthetic indole alkaloids for the amelioration of anxiety and depression-related disorders, structure-activity relationship, and their therapeutic prospects have been discussed.

Acinetobacter baumannii Sequence Types Harboring Genes Encoding Aminoglycoside Modifying Enzymes and 16SrRNA Methylase; a Multicenter Study from Pakistan.

INTRODUCTION: The aminoglycosides are widely used for the therapeutic management of infections caused by gram-negative bacteria, including the Acinetobacter baumannii strains. However, the resistance to the members of the aminoglycoside family, such as amikacin, gentamicin, and tobramycin, is increasingly being common among the clinical isolates. PURPOSE: This study aimed to investigate the presence of 16SrRNA methylases and aminoglycoside modifying enzymes (AMEs) genes among aminoglycoside resistant A. baumannii isolates and to study the genetic diversity of the clinical population of A. baumannii in local hospitals. MATERIAL AND METHODS: The 143 A. baumannii clinical strains were analyzed for antimicrobial susceptibility, genetic screening for enzymes conferring aminoglycosides resistance followed by the multilocus sequence typing. RESULTS: The 133/143 (93%) isolates were non-susceptible to at least one of the tested aminoglycosides, including amikacin, gentamicin, and tobramycin. The MIC distribution has shown that 87.486.7% strains were resistant to amikacin and gentamicin, respectively. The aphA6, aadB, aacC1, and aphA1 were found in 74.1%, 59.4%, 16.1%, and 11.2% isolates, respectively, whereas the armA was found in 28% of the strains having a higher MIC value (MIC; ≥256µg/mL). The MLST data have shown that the ST589 and ST2 were the most common STs and corresponded to 51 (35.7%) and 38 (26.6%) isolates, respectively, and few of the isolates corresponding to these STs were found to harbor the armA gene with a variable genotypic profile for AMEs. DISCUSSION: The study has reported the incidence of various enzymes conferring aminoglycoside resistance among the A. baumannii clones for the first time from Pakistan. The findings suggest the possibility of transmission of aminoglycoside resistance determinants through the lateral gene transfer as well as clonal dissemination.

Restoration of KMT2C/MLL3 in human colorectal cancer cells reinforces genome-wide H3K4me1 profiles and influences cell growth and gene expression.

BACKGROUND: The histone 3 lysine 4 (H3K4) monomethylase KMT2C is mutated across several cancer types; however, the effects of mutations on epigenome organization, gene expression, and cell growth are not clear. A frequently recurring mutation in colorectal cancer (CRC) with microsatellite instability is a single nucleotide deletion within the exon 38 poly-A(9) repeat (c.8390delA) which results in frameshift preceding the functional carboxy-terminal SET domain. To study effects of KMT2C expression in CRC cells, we restored one allele to wild type KMT2C in the two CRC cell lines RKO and HCT116, which both are homozygous c.8390delA mutant. RESULTS: Gene editing resulted in increased KMT2C expression, increased H3K4me1 levels, altered gene expression profiles, and subtle negative effects on cell growth, where higher dependence and stronger effects of KMT2C expression were observed in RKO compared to HCT116 cells. Surprisingly, we found that the two RKO and HCT116 CRC cell lines have distinct baseline H3K4me1 epigenomic profiles. In RKO cells, a flatter genome-wide H3K4me1 profile was associated with more increased H3K4me1 deposition at enhancers, reduced cell growth, and more differential gene expression relative to HCT116 cells when KMT2C was restored. Profiling of H3K4me1 did not indicate a highly specific regulation of gene expression as KMT2C-induced H3K4me1 deposition was found globally and not at a specific enhancer sub-set in the engineered cells. Although we observed variation in differentially regulated gene sets between cell lines and individual clones, differentially expressed genes in both cell lines included genes linked to known cancer signaling pathways, estrogen response, hypoxia response, and aspects of immune system regulation. CONCLUSIONS: Here, KMT2C restoration reduced CRC cell growth and reinforced genome-wide H3K4me1 deposition at enhancers; however, the effects varied depending upon the H3K4me1 status of KMT2C deficient cells. Results indicate that KMT2C inactivation may promote colorectal cancer development through transcriptional dysregulation in several pathways with known cancer relevance.

Linking FOXO3, NCOA3, and TCF7L2 to Ras pathway phenotypes through a genome-wide forward genetic screen in human colorectal cancer cells.

BACKGROUND: The Ras pathway genes KRAS, BRAF, or ERBBs have somatic mutations in ~ 60% of human colorectal carcinomas. At present, it is unknown whether the remaining cases lack mutations activating the Ras pathway or whether they have acquired mutations in genes hitherto unknown to belong to the pathway. METHODS: To address the second possibility and extend the compendium of Ras pathway genes, we used genome-wide transposon mutagenesis of two human colorectal cancer cell systems deprived of their activating KRAS or BRAF allele to identify genes enabling growth in low glucose, a Ras pathway phenotype, when targeted. RESULTS: Of the 163 recurrently targeted genes in the two different genetic backgrounds, one-third were known cancer genes and one-fifth had links to the EGFR/Ras/MAPK pathway. When compared to cancer genome sequencing datasets, nine genes also mutated in human colorectal cancers were identified. Among these, stable knockdown of FOXO3, NCOA3, and TCF7L2 restored growth in low glucose but reduced MEK/MAPK phosphorylation, reduced anchorage-independent growth, and modulated expressions of GLUT1 and Ras pathway related proteins. Knockdown of NCOA3 and FOXO3 significantly decreased the sensitivity to cetuximab of KRAS mutant but not wild-type cells. CONCLUSIONS: This work establishes a proof-of-concept that human cell-based genome-wide forward genetic screens can assign genes to pathways with clinical importance in human colorectal cancer.

Somatic PRDM2 c.4467delA mutations in colorectal cancers control histone methylation and tumor growth.

The chromatin modifier PRDM2/RIZ1 is inactivated by mutation in several forms of cancer and is a putative tumor suppressor gene. Frameshift mutations in the C-terminal region of PRDM2, affecting (A)8 or (A)9 repeats within exon 8, are found in one third of colorectal cancers with microsatellite instability, but the contribution of these mutations to colorectal tumorigenesis is unknown. To model somatic mutations in microsatellite unstable tumors, we devised a general approach to perform genome editing while stabilizing the mutated nucleotide repeat. We then engineered isogenic cell systems where the PRDM2 c.4467delA mutation in human HCT116 colorectal cancer cells was corrected to wild-type by genome editing. Restored PRDM2 increased global histone 3 lysine 9 dimethylation and reduced migration, anchorage-independent growth and tumor growth in vivo. Gene set enrichment analysis revealed regulation of several hallmark cancer pathways, particularly of epithelial-to-mesenchymal transition (EMT), with VIM being the most significantly regulated gene. These observations provide direct evidence that PRDM2 c.4467delA is a driver mutation in colorectal cancer and confirms PRDM2 as a cancer gene, pointing to regulation of EMT as a central aspect of its tumor suppressive action.

Loss of DIP2C in RKO cells stimulates changes in DNA methylation and epithelial-mesenchymal transition.

BACKGROUND: The disco-interacting protein 2 homolog C (DIP2C) gene is an uncharacterized gene found mutated in a subset of breast and lung cancers. To understand the role of DIP2C in tumour development we studied the gene in human cancer cells. METHODS: We engineered human DIP2C knockout cells by genome editing in cancer cells. The growth properties of the engineered cells were characterised and transcriptome and methylation analyses were carried out to identify pathways deregulated by inactivation of DIP2C. Effects on cell death pathways and epithelial-mesenchymal transition traits were studied based on the results from expression profiling. RESULTS: Knockout of DIP2C in RKO cells resulted in cell enlargement and growth retardation. Expression profiling revealed 780 genes for which the expression level was affected by the loss of DIP2C, including the tumour-suppressor encoding CDKN2A gene, the epithelial-mesenchymal transition (EMT) regulator-encoding ZEB1, and CD44 and CD24 that encode breast cancer stem cell markers. Analysis of DNA methylation showed more than 30,000 sites affected by differential methylation, the majority of which were hypomethylated following loss of DIP2C. Changes in DNA methylation at promoter regions were strongly correlated to changes in gene expression, and genes involved with EMT and cell death were enriched among the differentially regulated genes. The DIP2C knockout cells had higher wound closing capacity and showed an increase in the proportion of cells positive for cellular senescence markers. CONCLUSIONS: Loss of DIP2C triggers substantial DNA methylation and gene expression changes, cellular senescence and epithelial-mesenchymal transition in cancer cells.

Computational and molecular tools for scalable rAAV-mediated genome editing.

The rapid discovery of potential driver mutations through large-scale mutational analyses of human cancers generates a need to characterize their cellular phenotypes. Among the techniques for genome editing, recombinant adeno-associated virus (rAAV)-mediated gene targeting is suited for knock-in of single nucleotide substitutions and to a lesser degree for gene knock-outs. However, the generation of gene targeting constructs and the targeting process is time-consuming and labor-intense. To facilitate rAAV-mediated gene targeting, we developed the first software and complementary automation-friendly vector tools to generate optimized targeting constructs for editing human protein encoding genes. By computational approaches, rAAV constructs for editing ~71% of bases in protein-coding exons were designed. Similarly, ~81% of genes were predicted to be targetable by rAAV-mediated knock-out. A Gateway-based cloning system for facile generation of rAAV constructs suitable for robotic automation was developed and used in successful generation of targeting constructs. Together, these tools enable automated rAAV targeting construct design, generation as well as enrichment and expansion of targeted cells with desired integrations.

In situ sequencing identifies TMPRSS2-ERG fusion transcripts, somatic point mutations and gene expression levels in prostate cancers.

Translocations contribute to the genesis and progression of epithelial tumours and in particular to prostate cancer development. To better understand the contribution of fusion transcripts and visualize the clonal composition of multifocal tumours, we have developed a technology for multiplex in situ detection and identification of expressed fusion transcripts. When compared to immunohistochemistry, TMPRSS2-ERG fusion-negative and fusion-positive prostate tumours were correctly classified. The most prevalent TMPRSS2-ERG fusion variants were visualized, identified, and quantitated in human prostate cancer tissues, and the ratio of the variant fusion transcripts could for the first time be directly determined by in situ sequencing. Further, we demonstrate concurrent in situ detection of gene expression, point mutations, and gene fusions of the prostate cancer relevant targets AMACR, AR, TP53, and TMPRSS2-ERG. This unified approach to in situ analyses of somatic mutations can empower studies of intra-tumoural heterogeneity and future tissue-based diagnostics of mutations and translocations.

Molecular pathways in tumor progression: from discovery to functional understanding.

The advent of large scale sequencing methods has enabled analyses of the protein-coding parts of cancer genomes to find the mutated genes that cause common human cancers. Unbiased mutation analyses of human tumors originating in the breast, colon, brain, and pancreas have revealed genomic landscapes composed of a few frequently mutated genes alongside a multitude of infrequently mutated genes. These analyses have revealed a stark heterogeneity in the compendium of mutated genes even among tumors of the same tissue origin, and provide evidence for a larger number of driver mutations during tumorigenesis than hitherto presumed. From the multitude of mutated genes, a limited number of central molecular pathways are emerging. Systems biology approaches will be increasingly important to identify and better define these core pathways. Downstream of genetic analyses, scalable methods for prediction and experimental determination of the phenotypes of mutant alleles and pathways will be instrumental for improved mechanistic understanding of cancer as well as future drug discovery efforts.