Search

Abstract The innate immune stimulator of interferon genes (STING) pathway is known to activate type I interferons (IFN-I) and participate in generating antitumor immunity. We previously produced hDT806, a recombinant diphtheria immunotoxin, and demonstrated its efficacy against head and neck squamous cell carcinoma (HNSCC). However, it’s unknown whether the tumor-intrinsic STING plays a role in the anti-HNSCC effects of hDT806. In this study, we investigated the innate immune modulation of hDT806 on HNSCC. hDT806 significantly upregulated the level of STING and the ratio of p-TBK1/TBK1 in the HNSCC cells. Moreover, intratumoral hDT806 treatment increased the expression of STING-IFN-I signaling proteins including IFNA1, IFNB, CXCL10 and MX1, a marker of IFN-I receptor activity, in the HNSCC xenografts. Overexpression of STING mimicked the hDT806-induced upregulation of the STING-IFN-I signaling and induced apoptosis in the HNSCC cells. In the mouse xenograft models of HNSCC with STING overexpression, we observed a significant suppression of tumor growth and reduced tumor weight with increased apoptosis compared to their control xenograft counterparts without STING overexpression. Collectively, our data revealed that hDT806 may act as a stimulator of tumor-intrinsic STING-IFN-I signaling to inhibit tumor growth in HNSCC.

The global effort to combat the COVID-19 pandemic faces ongoing uncertainty with the emergence of Variants of Concern featuring numerous mutations on the Spike (S) protein. In particular, the Omicron Variant is distinguished by 32 mutations, including 10 within its receptor-binding domain (RBD). These mutations significantly impact viral infectivity and the efficacy of vaccines and antibodies currently in use for therapeutic purposes. In our study, we employed structure-based computational saturation mutagenesis approaches to predict the effects of Omicron missense mutations on RBD stability and binding affinity, comparing them to the original Wuhan-Hu-1 strain. Our results predict that mutations such as G431W and P507W induce the most substantial destabilizations in the Wuhan-Hu-1-S/Omicron-S RBD. Notably, we postulate that mutations in the Omicron-S exhibit a higher percentage of enhancing binding affinity compared to Wuhan-S. We found that the mutations at residue positions G447, Y449, F456, F486, and S496 led to significant changes in binding affinity. In summary, our findings may shed light on the widespread prevalence of Omicron mutations in human populations. The Omicron mutations that potentially enhance their affinity for human receptors may facilitate increased viral binding and internalization in infected cells, thereby enhancing infectivity. This informs the development of new neutralizing antibodies capable of targeting Omicron’s immune-evading mutations, potentially aiding in the ongoing battle against the COVID-19 pandemic.

The ongoing SARS-CoV-2/COVID-19 pandemic caused a global public health crisis. Yet, everyone’s response to SARS-CoV-2 infection varies, and different viral variants confer diverse pathogenicity. Thus, it is imperative to understand how viral determinants contribute to COVID-19. Viral ORF3a protein is one of those viral determinants, as its functions are linked to induction of cell and tissues damages, disease severity and cytokine storm that is a major cause of COVID-19-related death. ORF3a is a membrane-associated protein. Upon synthesis, it is transported from endoplasmic reticulum, Golgi apparatus to plasma membrane and subcellular endomembranes including endosomes and lysosomes. However, how ORF3a is transported intracellularly remains elusive. The goal of this study was to carry out a systematic mutagenesis study to determine the structural relationship of ORF3a protein with its subcellular locations. Single amino acid (aa) and deletion mutations were generated in the putative function-relevant motifs and other regions of interest. Immunofluorescence and ImageJ analyses were used to determine and quantitate subcellular locations of ORF3a mutants in comparison with wildtype ORF3a. The wildtype ORF3a localizes predominantly (Pearson’s coefficients about 0.8) on the membranes of endosomes and lysosomes. Consistent with earlier findings, deletion of the YXXΦ motif, which is required for protein export, retained ORF3a in the Golgi apparatus. Interestingly, mutations in a double glycine (diG) region (aa 187–188) displayed a similar phenotype to the YXXΦ deletion, implicating a similar role of the diG motif in intracellular transport. Indeed, interrupting any one of the two glycine residues such as deletion of a single (dG188), both (dG187/dG188) or substitution (G188Y) of these residues led to ORF3a retention in the Golgi apparatus (Pearson’s coefficients ≥0.8). Structural analyses further suggest that the diG motif supports a type-II β-turn between the anti-parallel β4 and β5 sheets and connects to the YXXΦ motif via hydrogen bonds between two monomers. The diG- YXXΦ interaction forms a hand-in-hand configuration that could facilitate dimerization. Together, these observations suggest a functional role of the diG motif in intracellular transport of ORF3a.

SARS-CoV-2 is a novel zoonotic positive-sense RNA virus (ssRNA+) belonging to the genus beta coronaviruses (CoVs) in the Coronaviridae family. It is the causative agent for the outbreak of the disease, COVID-19. It is the third CoV causing pneumonia around the world in the past 2 decades. To date, it has caused significant deaths worldwide. Notably, the emergence of new genetic variants conferring efficient transmission and immune evasion remained a challenge, despite the reduction in the number of death cases, owing to effective vaccination regimen (boosting) and safety protocols. Thus, information harnessed from SARS-CoV-2 genomic organization is indispensable for seeking laboratory diagnosis and treatment options. Here in, we review previously circulating variants of SARS-CoV-2 designated variant of concern (VOC) including the Alpha (United Kingdom), Beta (South Africa), Gamma (Brazil), Delta (India), and recently circulating VOC, Omicron (South Africa) and its divergent subvariants (BA.1, BA.2, BA.3, BA.2.12.1, BA.4 and BA.5) with BA.5 currently becoming dominant and prolonging the COVID pandemic. In addition, we address the role of computational models for mutagenesis analysis which can predict important residues that contribute to transmissibility, virulence, immune evasion, and molecular detections of SARS-CoV-2. Concomitantly, the importance of harnessing the immunobiology of SARS‐CoV‐2 and host interaction for therapeutic purpose; and use of an in slilico based biocomputational approaches to achieve this purpose via predicting novel therapeutic agents targeting PRR such as toll like receptor, design of universal vaccine and chimeric antibodies tailored to the emergent variant have been highlighted.

While worldwide efforts for improving COVID-19 vaccines are currently considered a top priority, the role of the genetic variants responsible for virus receptor protein stability is less studied. Angiotensin-converting enzyme-2 is the primary target of the SARS-CoV-1/SARS-CoV-2 spike (S) glycoprotein, enabling entry into the human body. Here, we applied computational saturation mutagenesis approaches to determine the folding energy caused by all possible mutations in ACE2 proteins within ACE2 - SARS-CoV-1-S/ACE2 - SARS-CoV-2-S complexes. We observed ACE2 mutations at residue D350 causing the most stabilizing effects on the protein. In addition, we identified ACE2 genetic variations in African Americans (rs73635825, rs766996587, and rs780574871), Latino Americans (rs924799658), and both groups (rs4646116 and rs138390800) affecting stability in the ACE2 - SARS-CoV-2-S complex. The findings in this study may aid in targeting the design of stable neutralizing peptides for treating minority patients.

Middle East Respiratory Syndrome Coronavirus (MERS-CoV) causes severe pneumonia-like symptoms and is still pose a significant threat to global public health. A key component in the virulence of MERS-CoV is the Spike (S) protein, which binds with the host membrane receptor dipeptidyl peptidase 4 (DPP4). The goal of the present investigation is to examine the effects of missense mutations in the MERS-CoV S protein on protein stability and binding affinity with DPP4 to provide insight that is useful in developing vaccines to prevent coronavirus infection. We utilized a saturation mutagenesis approach to simulate all possible mutations in the MERS-CoV full-length S, S Receptor Binding Domain (RBD) and DPP4. We found the mutations in MERS-CoV S protein residues, G552, C503, C526, N468, G570, S532, S451, S419, S465, and S435, affect protein stability. We identified key residues, G538, E513, V555, S557, L506, L507, R511, M452, D537, and S454 in the S protein RBD region are important in the binding of MERS-CoV S protein to the DPP4 receptor. We investigated the effects of MERS-CoV S protein viral mutations on protein stability and binding affinity. In addition, we studied all DPP4 mutations and found the functional substitution R336T weakens both DPP4 protein stability and S-DPP4 binding affinity. We compared the S protein structures of MERS-CoV, SARS-CoV, and SARS-CoV-2 viruses and identified the residues like C526, C383, and N468 located in equivalent positions of these viruses have effects on S protein structure. These findings provide further information on how mutations in coronavirus S proteins effect protein function.

COVID-19, caused by SARS-CoV-2, is a systemic illness due to its multiorgan effects in patients. The disease has a detrimental impact on respiratory and cardiovascular systems. One early symptom of infection is anosmia or lack of smell; this implicates the involvement of the olfactory bulb in COVID-19 disease and provides a route into the central nervous system. However, little is known about how SARS-CoV-2 affects neurological or psychological symptoms. SARS-CoV-2 exploits host receptors that converge on pathways that impact psychological symptoms. This systemic review discusses the ways involved by coronavirus infection and their impact on mental health disorders. We begin by briefly introducing the history of coronaviruses, followed by an overview of the essential proteins to viral entry. Then, we discuss the downstream effects of viral entry on host proteins. Finally, we review the literature on host factors that are known to play critical roles in neuropsychiatric symptoms and mental diseases and discuss how COVID-19 could impact mental health globally. Our review details the host factors and pathways involved in the cellular mechanisms, such as systemic inflammation, that play a significant role in the development of neuropsychological symptoms stemming from COVID-19 infection.

ACE2 plays a critical role in SARS-CoV-2 infection to cause COVID-19 and SARS-CoV-2 spike protein binds to ACE2 and probably functionally inhibits ACE2 to aggravate the underlying diseases of COVID-19. The important factors that affect the severity and fatality of COVID-19 include patients' underlying diseases and ages. Therefore, particular care to the patients with underlying diseases is needed during the treatment of COVID-19 patients.

Severe Acute respiratory syndrome coronavirus (SARS-CoV-1) attaches to the host cell surface to initiate the interaction between the receptor-binding domain (RBD) of its spike glycoprotein (S) and the human Angiotensin-converting enzyme (hACE2) receptor. SARS-CoV-1 mutates frequently because of its RNA genome, which challenges the antiviral development. Here, we per-formed computational saturation mutagenesis of the S protein of SARS-CoV-1 to identify the residues crucial for its functions. We used the structure-based energy calculations to analyze the effects of the missense mutations on the SARS-CoV-1 S stability and the binding affinity with hACE2. The sequence and structure alignment showed similarities between the S proteins of SARS-CoV-1 and SARS-CoV-2. Interestingly, we found that target mutations of S protein amino acids generate similar effects on their stabilities between SARS-CoV-1 and SARS-CoV-2. For example, G839W of SARS-CoV-1 corresponds to G857W of SARS-CoV-2, which decrease the stability of their S glycoproteins. The viral mutation analysis of the two different SARS-CoV-1 isolates showed that mutations, T487S and L472P, weakened the S-hACE2 binding of the 2003–2004 SARS-CoV-1 isolate. In addition, the mutations of L472P and F360S destabilized the 2003–2004 viral isolate. We further predicted that many mutations on N-linked glycosylation sites would increase the stability of the S glycoprotein. Our results can be of therapeutic importance in the design of antivirals or vaccines against SARS-CoV-1 and SARS-CoV-2.

Abstract Background Lymphotoxin-α (LTα), located in the Major Histocompatibility Complex (MHC) class III region on chromosome 6, encodes a cytotoxic protein that mediates a variety of antiviral responses among other biological functions. Furthermore, several genotypes at this gene have been implicated in the onset of a number of complex diseases, including myocardial infarction, autoimmunity, and various types of cancer. However, little is known about levels of nucleotide variation and linkage disequilibrium (LD) in and near LTα, which could also influence phenotypic variance. To address this gap in knowledge, we examined sequence variation across ~ 10 kilobases (kbs), encompassing LTα and the upstream region, in 2039 individuals from the 1000 Genomes Project originating from 21 global populations. Results Here, we observed striking patterns of diversity, including an excess of intermediate-frequency alleles, the maintenance of multiple common haplotypes and a deep coalescence time for variation (dating > 1.0 million years ago), in global populations. While these results are generally consistent with a model of balancing selection, we also uncovered a signature of positive selection in the form of long-range LD on chromosomes with derived alleles primarily in Eurasian populations. To reconcile these findings, which appear to support different models of selection, we argue that selective sweeps (particularly, soft sweeps) of multiple derived alleles in and/or near LTα occurred in non-Africans after their ancestors left Africa. Furthermore, these targets of selection were predicted to alter transcription factor binding site affinity and protein stability, suggesting they play a role in gene function. Additionally, our data also showed that a subset of these functional adaptive variants are present in archaic hominin genomes. Conclusions Overall, this study identified candidate functional alleles in a biologically-relevant genomic region, and offers new insights into the evolutionary origins of these loci in modern human populations.

The ongoing outbreak of COVID-19 has been a serious threat to human health worldwide. The virus SARS-CoV-2 initiates its infection to the human body via the interaction of its spike (S) protein with the human Angiotensin-Converting Enzyme 2 (ACE2) of the host cells. Therefore, understanding the fundamental mechanisms of how SARS-CoV-2 S protein receptor binding domain (RBD) binds to ACE2 is highly demanded for developing treatments for COVID-19. Here we implemented multi-scale computational approaches to study the binding mechanisms of human ACE2 and S proteins of both SARS-CoV and SARS-CoV-2. Electrostatic features, including electrostatic potential, electric field lines, and electrostatic forces of SARS-CoV and SARS-CoV-2 were calculated and compared in detail. The results demonstrate that SARS-CoV and SARS-CoV-2 S proteins are both attractive to ACE2 by electrostatic forces even at different distances. However, the residues contributing to the electrostatic features are quite different due to the mutations between SARS-CoV S protein and SARS-CoV-2 S protein. Such differences are analyzed comprehensively. Compared to SARS-CoV, the SARS-CoV-2 binds with ACE2 using a more robust strategy: The electric field line related residues are distributed quite differently, which results in a more robust binding strategy of SARS-CoV-2. Also, SARS-CoV-2 has a higher electric field line density than that of SARS-CoV, which indicates stronger interaction between SARS-CoV-2 and ACE2, compared to that of SARS-CoV. Key residues involved in salt bridges and hydrogen bonds are identified in this study, which may help the future drug design against COVID-19.

COVID-19 has been threatening human health since the late 2019, and has a significant impact on human health and economy. Understanding SARS-CoV-2 and other coronaviruses is important to develop effective treatments for COVID-19 and other coronavirus-caused diseases. In this work, we applied multi-scale computational approaches to study the electrostatic features of spike (S) proteins for SARS-CoV and SARS-CoV-2. From our results, we found that SARS-CoV and SARS-CoV-2 have similar charge distributions and electrostatic features when binding with the human angiotensin-converting enzyme 2 (hACE2). Energy pH-dependence calculations revealed that the complex structures of hACE2 and the S proteins of SARS-CoV/SARS-CoV-2 are stable at pH values ranging from 7.5 to 9. Three independent 100 ns molecular dynamics (MD) simulations were performed using NAMD to investigate the hydrogen bonds between S proteins RBD and hACE2 RBD. From MD simulations, we found that SARS-CoV-2 forms 19 pairs (average of three simulations) of hydrogen bonds with high occupancy (>50%) to hACE2, compared to 16 pairs between SARS-CoV and hACE2. Additionally, SARS-CoV viruses prefer sticking to the same hydrogen bond pairs, while SARS-CoV-2 tends to have a larger range of selections on hydrogen bonds acceptors. We also labelled key residues involved in forming the top five hydrogen bonds that were found in all three independent 100 ns simulations. This identification is important to potential drug designs for COVID-19 treatments. Our work will shed the light on current and future coronavirus-caused diseases.

Abstract The ability to predict if a given mutation is disease-causing or not has enormous potential to impact human health. Typically, these predictions are made by assessing the effects of mutation on macromolecular stability and amino acid conservation. Here we report a novel feature: the electrostatic component of the force acting between a kinesin motor domain and tubulin. We demonstrate that changes in the electrostatic component of the binding force are able to discriminate between disease-causing and non-disease-causing mutations found in human kinesin motor domains using the receiver operating characteristic (ROC). Because diseases may originate from multiple effects not related to kinesin-microtubule binding, the prediction rate of 0.843 area under the ROC plot due to the change in magnitude of the electrostatic force alone is remarkable. These results reflect the dependence of kinesin’s function on motility along the microtubule, which suggests a precise balance of microtubule binding forces is required.

After publication of our article [1] we were notified that a few duplicate sentences were included on Figure 3 and Figure 4 legends.

Schizophrenia (SCZ) is a serious psychiatric disorder that affects 1% of general population and places a heavy burden worldwide. The underlying genetic mechanism of SCZ remains unknown, but studies indicate that the disease is associated with a global gene expression disturbance across many genes. Next-generation sequencing, particularly of RNA sequencing (RNA-Seq), provides a powerful genome-scale technology to investigate the pathological processes of SCZ. RNA-Seq has been used to analyze the gene expressions and identify the novel splice isoforms and rare transcripts associated with SCZ. This paper provides an overview on the genetics of SCZ, the advantages of RNA-Seq for transcriptome analysis, the accomplishments of RNA-Seq in SCZ cohorts, and the applications of induced pluripotent stem cells and RNA-Seq in SCZ research.

Many studies have shown that missense mutations might play an important role in carcinogenesis. However, the extent to which cancer mutations might affect biomolecular interactions remains unclear. Here, we map glioblastoma missense mutations on the human protein interactome, model the structures of affected protein complexes and decipher the effect of mutations on protein-protein, protein-nucleic acid and protein-ion binding interfaces. Although some missense mutations over-stabilize protein complexes, we found that the overall effect of mutations is destabilizing, mostly affecting the electrostatic component of binding energy. We also showed that mutations on interfaces resulted in more drastic changes of amino acid physico-chemical properties than mutations occurring outside the interfaces. Analysis of glioblastoma mutations on interfaces allowed us to stratify cancer-related interactions, identify potential driver genes, and propose two dozen additional cancer biomarkers, including those specific to functions of the nervous system. Such an analysis also offered insight into the molecular mechanism of the phenotypic outcomes of mutations, including effects on complex stability, activity, binding and turnover rate. As a result of mutated protein and gene network analysis, we observed that interactions of proteins with mutations mapped on interfaces had higher bottleneck properties compared to interactions with mutations elsewhere on the protein or unaffected interactions. Such observations suggest that genes with mutations directly affecting protein binding properties are preferably located in central network positions and may influence critical nodes and edges in signal transduction networks.