Assessing Differential Binding of Aggregation-Induced Emission-Based Luminogens to Host Interacting Surface Proteins of SARS-CoV-2 and Influenza Virus–An in silico Approach

Graphical AbstractProposed design for AIE-based differential detection of SARS-CoV-2 and influenza viruses.

Replication Study: RAF inhibitors prime wild-type RAF to activate the MAPK pathway and enhance growth

As part of the Reproducibility Project: Cancer Biology, we published a Registered Report (Bhargava et al., 2016) that described how we intended to replicate selected experiments from the paper "RAF inhibitors prime wild-type RAF to activate the MAPK pathway and enhance growth" (Hatzivassiliou et al., 2010). Here we report the results. We found two unrelated RAF inhibitors, PLX4720 or GDC-0879, selectively inhibited BRAF(V600E) cell proliferation, while the MEK inhibitor, PD0325901, inhibited BRAF(V600E), wild-type RAF/RAS, and mutant RAS cancer cell proliferation, similar to the original study (Figure 1A; Hatzivassiliou et al., 2010). We found knockdown of CRAF, but not BRAF, in mutant RAS cells attenuated the phospho-MEK induction observed after PLX4720 treatment, similar to the original study (Figure 2B; Hatzivassiliou et al., 2010). The original study reported analogous results with GDC-0879, which was not observed in this replication, although unexpected control results confound the interpretation. We also attempted a replication of an assay with recombinant proteins to test the differential effect of RAF inhibitors on BRAF-CRAF heterodimerization (Figure 4A; Hatzivassiliou et al., 2010). Although we were unable to conduct the experiment as planned, we observed differential binding of BRAF by RAF inhibitors; however, it was between BRAF and beads, independent of CRAF. While these data were unable to address whether, under the conditions of the original study, the same observations could be observed, we discuss key differences between the original study and this replication that are important to consider for further experiments. Finally, where possible, we report meta-analyses for each result.

Identification of ligand binding sites in intrinsically disordered proteins with a differential binding score

Screening ligands directly binding to an ensemble of intrinsically disordered proteins (IDP) to discover potential hits or leads for new drugs is an emerging but challenging area as IDPs lack well-defined and ordered 3D-protein structures. To explore a new IDP-based rational drug discovery strategy, a differential binding score (DIBS) is defined. The basis of DIBS is to quantitatively determine the binding preference of a ligand to an ensemble of conformations specified by IDP versus such preferences to an ensemble of random coil conformations of the same protein. Ensemble docking procedures performed on repeated sampling of conformations, and the results tested for statistical significance determine the preferential ligand binding sites of the IDP. The results of this approach closely reproduce the experimental data from recent literature on the binding of the ligand epigallocatechin gallate (EGCG) to the intrinsically disordered N-terminal domain of the tumor suppressor p53. Combining established approaches in developing a new method to screen ligands against IDPs could be valuable as a screening tool for IDP-based drug discovery.

Nothing Regular about the Regulins: Distinct Functional Properties of SERCA Transmembrane Peptide Regulatory Subunits

The sarco-endoplasmic reticulum calcium ATPase (SERCA) is responsible for maintaining calcium homeostasis in all eukaryotic cells by actively transporting calcium from the cytosol into the sarco-endoplasmic reticulum (SR/ER) lumen. Calcium is an important signaling ion, and the activity of SERCA is critical for a variety of cellular processes such as muscle contraction, neuronal activity, and energy metabolism. SERCA is regulated by several small transmembrane peptide subunits that are collectively known as the “regulins”. Phospholamban (PLN) and sarcolipin (SLN) are the original and most extensively studied members of the regulin family. PLN and SLN inhibit the calcium transport properties of SERCA and they are required for the proper functioning of cardiac and skeletal muscles, respectively. Myoregulin (MLN), dwarf open reading frame (DWORF), endoregulin (ELN), and another-regulin (ALN) are newly discovered tissue-specific regulators of SERCA. Herein, we compare the functional properties of the regulin family of SERCA transmembrane peptide subunits and consider their regulatory mechanisms in the context of the physiological and pathophysiological roles of these peptides. We present new functional data for human MLN, ELN, and ALN, demonstrating that they are inhibitors of SERCA with distinct functional consequences. Molecular modeling and molecular dynamics simulations of SERCA in complex with the transmembrane domains of MLN and ALN provide insights into how differential binding to the so-called inhibitory groove of SERCA—formed by transmembrane helices M2, M6, and M9—can result in distinct functional outcomes.

Evolution of binding preferences among whole-genome duplicated transcription factors

Throughout evolution, new transcription factors (TFs) emerge by gene duplication, promoting growth and rewiring of transcriptional networks. How TF duplicates diverge is known for only a few studied cases. To provide a genome-scale view, we considered the 35% of budding yeast TFs, classified as whole-genome duplication (WGD)-retained paralogs. Using high-resolution profiling, we find that ~60% of paralogs evolved differential binding preferences. We show that this divergence results primarily from variations outside the DNA binding domains (DBDs), while DBD preferences remain largely conserved. Analysis of non-WGD orthologs revealed that ancestral preferences are unevenly split between duplicates, while new targets are acquired preferentially by the least conserved paralog (biased sub/neo-functionalization). Dimer-forming paralogs evolved mostly one-sided dependency, while other paralogs interacted through low-magnitude DNA-binding competition that minimized paralog interference. We discuss the implications of our findings for the evolutionary design of transcriptional networks.

Differential Binding of Autoantibodies to MOG Isoforms in Inflammatory Demyelinating Diseases

ObjectiveTo analyze serum immunoglobulin G (IgG) antibodies to major isoforms of myelin oligodendrocyte glycoprotein (MOG-alpha 1-3 and beta 1-3) in patients with inflammatory demyelinating diseases.MethodsRetrospective case-control study using 378 serum samples from patients with multiple sclerosis (MS), patients with non-MS demyelinating disease, and healthy controls with MOG alpha-1-IgG positive (n = 202) or negative serostatus (n = 176). Samples were analyzed for their reactivity to human, mouse, and rat MOG isoforms with and without mutations in the extracellular MOG Ig domain (MOG-ecIgD), soluble MOG-ecIgD, and myelin from multiple species using live cell-based, tissue immunofluorescence assays and ELISA.ResultsThe strongest IgG reactivities were directed against the longest MOG isoforms alpha-1 (the currently used standard test for MOG-IgG) and beta-1, whereas the other isoforms were less frequently recognized. Using principal component analysis, we identified 3 different binding patterns associated with non-MS disease: (1) isolated reactivity to MOG-alpha-1/beta-1 (n = 73), (2) binding to MOG-alpha-1/beta-1 and at least one other alpha, but no beta isoform (n = 64), and (3) reactivity to all 6 MOG isoforms (n = 65). The remaining samples were negative (n = 176) for MOG-IgG. These MOG isoform binding patterns were associated with a non-MS demyelinating disease, but there were no differences in clinical phenotypes or disease course. The 3 MOG isoform patterns had distinct immunologic characteristics such as differential binding to soluble MOG-ecIgD, sensitivity to MOG mutations, and binding to human MOG in ELISA.ConclusionsThe novel finding of differential MOG isoform binding patterns could inform future studies on the refinement of MOG-IgG assays and the pathophysiologic role of MOG-IgG.

A flexible ChIP-sequencing simulation toolkit

Background A major challenge in evaluating quantitative ChIP-seq analyses, such as peak calling and differential binding, is a lack of reliable ground truth data. Accurate simulation of ChIP-seq data can mitigate this challenge, but existing frameworks are either too cumbersome to apply genome-wide or unable to model a number of important experimental conditions in ChIP-seq. Results We present ChIPs, a toolkit for rapidly simulating ChIP-seq data using statistical models of key experimental steps. We demonstrate how ChIPs can be used for a range of applications, including benchmarking analysis tools and evaluating the impact of various experimental parameters. ChIPs is implemented as a standalone command-line program written in C++ and is available from https://github.com/gymreklab/chips. Conclusions ChIPs is an efficient ChIP-seq simulation framework that generates realistic datasets over a flexible range of experimental conditions. It can serve as an important component in various ChIP-seq analyses where ground truth data are needed.

Differential binding of SARS-CoV-2 Spike protein variants to its cognate receptor hACE2 using molecular modeling based binding analysis

The current emergence of novel coronavirus, SARS-CoV-2 and its ceaseless expansion worldwide has posed a global health emergency that has adversely affected the humans. With the entire world striving to understand the newly emerged virus, differences in morbidity and infection rate of SARS-CoV-2 have been observed across varied geographic areas, which have been ascribed to viral mutation and evolution over time. The homotrimeric Spike (S) glycoprotein on the viral envelope surface is responsible for binding, priming, and initiating infection in the host. Our phylogeny analysis of 1947 sequences of S proteins indicated there is a change in amino acid (aa) from aspartate (Group-A) to glycine (Group-B) at position 614, near the receptor- binding domain (RBD; aa positions 331-524). The two variants are reported to be in circulation, disproportionately across the world, with Group-A dominant in Asia and Group-B in North America. The trimeric, monomeric, and RBD of S protein of both the variant groups (A & B) were modeled using the Swiss-Model server and were docked with the human receptor angiotensin-converting enzyme 2 (hACE2) employing the PatchDock server and visualized in PyMol. Group-A S protein’s RBD bound imperceptibly to the two binding clefts of the hACE2 protein, on the other hand, Group-B S protein’s RBD perfectly interacted inside the binding clefts of hACE2, with higher number of hydrogen and hydrophobic interactions. This implies that the S protein’s amino acid at position 614 near the core RBD influences its interaction with the cognate hACE2 receptor, which may induce its infectivity that should be explored further with molecular and biochemical studies.