Profiling Subcellular Protein Phosphatase Responses to Coxsackievirus B3 Infection of Cardiomyocytes

Proteome-wide copy-number estimation from transcriptomics

Protein copy numbers constrain systems-level properties of regulatory networks, but proportional proteomic data remain scarce compared to RNA-seq. We related mRNA to protein statistically using best-available data from quantitative proteomics and transcriptomics for 4366 genes in 369 cell lines. The approach starts with a protein's median copy number and hierarchically appends mRNA-protein and mRNA-mRNA dependencies to define an optimal gene-specific model linking mRNAs to protein. For dozens of cell lines and primary samples, these protein inferences from mRNA outmatch stringent null models, a count-based protein-abundance repository, empirical mRNA-to-protein ratios, and a proteogenomic DREAM challenge winner. The optimal mRNA-to-protein relationships capture biological processes along with hundreds of known protein-protein complexes, suggesting mechanistic relationships. We use the method to identify a viral-receptor abundance threshold for coxsackievirus B3 susceptibility from 1489 systems-biology infection models parameterized by protein inference. When applied to 796 RNA-seq profiles of breast cancer, inferred copy-number estimates collectively re-classify 26-29% of luminal tumors. By adopting a gene-centered perspective of mRNA-protein covariation across different biological contexts, we achieve accuracies comparable to the technical reproducibility of contemporary proteomics.

Three Modes of Viral Adaption by the Heart

Viruses elicit long-term adaptive responses in the tissues they infect. Understanding viral adaptions in humans is difficult in organs such as the heart, where primary infected material is not routinely collected. In search of asymptomatic infections with accompanying host adaptions, we mined for cardio-pathogenic viruses in the unaligned reads of nearly one thousand human hearts profiled by RNA sequencing. Among virus-positive cases (~20%), we identified three robust adaptions in the host transcriptome related to inflammatory NFκB signaling and post-transcriptional regulation by the p38-MK2 pathway. The adaptions are not determined by the infecting virus, and they recur in infections of human or animal hearts and cultured cardiomyocytes. Adaptions switch states when NFκB or p38-MK2 are perturbed in cells engineered for chronic infection by the cardio-pathogenic virus, coxsackievirus B3. Stratifying viral responses into reversible adaptions adds a targetable systems-level simplification for infections of the heart and perhaps other organs.

A novel site on dual-specificity phosphatase MKP7/DUSP16 is required for catalysis and MAPK binding

The dual-specificity phosphatases responsible for the inactivation of the mitogen-activated protein kinases (MAPKs) are designated as the MAPK phosphatases (MKPs). We demonstrated previously that MKP5 is regulated through a novel allosteric site suggesting additional regulatory mechanisms of catalysis exist amongst the MKPs. Here, we sought to determine whether the equivalent site within the phosphatase domain of a highly similar MKP family member, MKP7, is also important for phosphatase function. We found that mutation of tyrosine 271 (Y271) in MKP7, which represents the comparable Y435 within the MKP5 allosteric pocket, inhibited MKP7 catalytic activity. Consistent with this, when MKP7 Y271 mutants were overexpressed in cells, the substrates of MKP7, p38 MAPK or JNK, failed to undergo dephosphorylation. The binding efficiency of MKP7 to p38 MAPK and JNK1/2 was also reduced when MKP7 Y271 is mutated. Consistent with reduced MAPK binding, we observed a greater accumulation of nuclear p38 MAPK and JNK when the MKP7 Y271 mutants are expressed in cells as compared with WT MKP7, which sequesters p38 MAPK/JNK in the cytoplasm. Therefore, we propose that Y271 is critical for effective MAPK dephosphorylation through a mechanism whereby binding to this residue precedes engagement of the catalytic site and upon overexpression, MKP7 allosteric site mutants potentiate MAPK signaling. These results provide insight into the regulatory mechanisms of MKP7 catalysis and interactions with the MAPKs. Furthermore, these data support the generality of the MKP allosteric site and provide a basis for small molecule targeting of MKP7.

Rta is the principal activator of Epstein-Barr virus epithelial lytic transcription

The transition from latent Epstein-Barr virus (EBV) infection to lytic viral replication is mediated by the viral transcription factors Rta and Zta. Although both are required for virion production, dissecting the specific roles played by Rta and Zta is challenging because they induce each other's expression. To circumvent this, we constructed an EBV mutant deleted for the genes encoding Rta and Zta (BRLF1 and BZLF1, respectively) in the Akata strain BACmid. This mutant, termed EBVΔRZ, was used to infect several epithelial cell lines, including telomerase-immortalized normal oral keratinocytes, a highly physiologic model of EBV epithelial cell infection. Using RNA-seq, we determined the gene expression induced by each viral transactivator. Surprisingly, Zta alone only induced expression of the lytic origin transcripts BHLF1 and LF3. In contrast, Rta activated the majority of EBV early gene transcripts. As expected, Zta and Rta were both required for expression of late gene transcripts. Zta also cooperated with Rta to enhance a subset of early gene transcripts (Rtasynergy transcripts) that Zta was unable to activate when expressed alone. Interestingly, Rta and Zta each cooperatively enhanced the other's binding to EBV early gene promoters, but this effect was not restricted to promoters where synergy was observed. We demonstrate that Zta did not affect Rtasynergy transcript stability, but increased Rtasynergy gene transcription despite having no effect on their transcription when expressed alone. Our results suggest that, at least in epithelial cells, Rta is the dominant transactivator and that Zta functions primarily to support DNA replication and co-activate a subset of early promoters with Rta. This closely parallels the arrangement in KSHV where ORF50 (Rta homolog) is the principal activator of lytic transcription and K8 (Zta homolog) is required for DNA replication at oriLyt.

Modeling the complete kinetics of coxsackievirus B3 reveals human determinants of host-cell feedback

Complete kinetic models are pervasive in chemistry but lacking in biological systems. We encoded the complete kinetics of infection for coxsackievirus B3 (CVB3), a compact and fast-acting RNA virus. The model consists of separable, detailed modules describing viral binding-delivery, translation-replication, and encapsidation. Specific module activities are dampened by the type I interferon response to viral double-stranded RNAs (dsRNAs), which is itself disrupted by viral proteinases. The experimentally validated kinetics uncovered that cleavability of the dsRNA transducer mitochondrial antiviral signaling protein (MAVS) becomes a stronger determinant of viral outcomes when cells receive supplemental interferon after infection. Cleavability is naturally altered in humans by a common MAVS polymorphism, which removes a proteinase-targeted site but paradoxically elevates CVB3 infectivity. These observations are reconciled with a simple nonlinear model of MAVS regulation. Modeling complete kinetics is an attainable goal for small, rapidly infecting viruses and perhaps viral pathogens more broadly. A record of this paper's transparent peer review process is included in the Supplemental information.

Fragile epitopes-Antibody's guess is as good as yours

Monoclonal antibodies recognize epitopes so specifically that altering a single residue can disrupt binding. In this issue of Science Signaling, Schüchner et al and Frohner et al report that flanking amino acids and an underappreciated posttranslational modification perturb epitope affinity for two groups of widely used monoclonal antibodies.

An ultrasensitive fiveplex activity assay for cellular kinases

Protein kinases are enzymes whose abundance, protein-protein interactions, and posttranslational modifications together determine net signaling activity in cells. Large-scale data on cellular kinase activity are limited, because existing assays are cumbersome, poorly sensitive, low throughput, and restricted to measuring one kinase at a time. Here, we surmount the conventional hurdles of activity measurement with a multiplexing approach that leverages the selectivity of individual kinase-substrate pairs. We demonstrate proof of concept by designing an assay that jointly measures activity of five pleiotropic signaling kinases: Akt, IκB kinase (IKK), c-jun N-terminal kinase (JNK), mitogen-activated protein kinase (MAPK)-extracellular regulated kinase kinase (MEK), and MAPK-activated protein kinase-2 (MK2). The assay operates in a 96-well format and specifically measures endogenous kinase activation with coefficients of variation less than 20%. Multiplex tracking of kinase-substrate pairs reduces input requirements by 25-fold, with ~75 µg of cellular extract sufficient for fiveplex activity profiling. We applied the assay to monitor kinase signaling during coxsackievirus B3 infection of two different host-cell types and identified multiple differences in pathway dynamics and coordination that warrant future study. Because the Akt–IKK–JNK–MEK–MK2 pathways regulate many important cellular functions, the fiveplex assay should find applications in inflammation, environmental-stress, and cancer research.

Microenvironmental Signals and Biochemical Information Processing: Cooperative Determinants of Intratumoral Plasticity and Heterogeneity

Intra-tumor cellular heterogeneity is a major challenge in cancer therapy. Tumors are composed of multiple phenotypic subpopulations that vary in their ability to initiate metastatic tumors and in their sensitivity to chemotherapy. In many cases, cells can transition between these subpopulations, not by genetic mutation, but instead through reversible changes in signal transduction or gene expression programs. This plasticity begins at the level of the microenvironment where local autocrine and paracrine signals, exosomes, tumor–stroma interactions, and extracellular matrix (ECM) composition create a signaling landscape that varies over space and time. The integration of this complex array of signals engages signaling pathways that control gene expression. The resulting modulation of gene expression programs causes individual cells to sample a wide array of phenotypic states that support tumor growth, dissemination, and therapeutic resistance. In this review, we discuss how information flows dynamically within the microenvironmental landscape to inform cell state decisions and to create intra-tumoral heterogeneity. We address the role of plasticity in the acquisition of transient and prolonged drug resistant states and discuss how targeted pharmacological modification of the signaling landscape may be able to constrain phenotypic plasticity, leading to improved treatment responses.

Pin1 facilitates isoproterenol‑induced cardiac fibrosis and collagen deposition by promoting oxidative stress and activating the MEK1/2‑ERK1/2 signal transduction pathway in rats

Peptidyl-prolyl cis/trans isomerase, NIMA-interacting 1 (Pin1) is a member of a large superfamily of phosphorylation-dependent peptidyl-prolyl cis/trans isomerases, which not only regulates multiple targets at various stages of cellular processes, but is also involved in the pathogenesis of several diseases, including microbial infection, cancer, asthma and Alzheimer's disease. However, the role of Pin1 in cardiac fibrosis remains to be fully elucidated. The present study investigated the potential mechanism of Pin1 in isoprenaline (ISO)-induced myocardial fibrosis in rats. The rats were randomly divided into three groups. Echocardiography was used to evaluate changes in the size, shape and function of the heart, and histological staining was performed to visualize inflammatory cell infiltration and fibrosis. Reverse transcription-quantitative polymerase chain reaction analysis, immunohistochemistry and Picrosirius red staining were used to differentiate collagen subtypes. Additionally, cardiac-specific phosphorylation of mitogen-activated protein kinase kinase 1/2 (MEK1/2) and extracellular-signal regulated protein kinase 1/2 (ERK1/2), and the activities of Pin1 and α-smooth muscle actin (α-SMA) and other oxidative stress parameters were estimated in the heart. The administration of ISO resulted in an increase in cardiac parameters and elevated the heart-to-body weight ratio. Histopathological examination of heart tissues revealed interstitial inflammatory cellular infiltrate and disorganized collagen fiber deposition. In addition, lipid peroxidation products and oxidative stress marker activity in plasma and tissues were significantly increased in the ISO-treated rats. Western blot analysis showed significantly elevated protein levels of phosphorylated Pin1, MEK1/2, ERK1/2 and α-SMA in remodeling hearts. Treatment with juglone following intraperitoneal injection of ISO significantly prevented inflammatory cell infiltration, improved cardiac function, and suppressed oxidative stresses and fibrotic alterations. In conclusion, the results of the present study suggested that the activation of Pin1 promoted cardiac extracellular matrix deposition and oxidative stress damage by regulating the phosphorylation of the MEK1/2-ERK1/2 signaling pathway and the expression of α-SMA. By contrast, the inhibition of Pin1 alleviated cardiac damage and fibrosis in the experimental models, suggesting that Pin1 contributed to the development of cardiac remodeling in ISO-administered rats, and that the inactivation of Pin1 may be a novel therapeutic candidate for the treatment of cardiovascular disease and heart failure.

Linear Integration of ERK Activity Predominates over Persistence Detection in Fra-1 Regulation

ERK signaling regulates the expression of target genes, but it is unclear how ERK activity dynamics are interpreted. Here, we investigate this question using simultaneous, live, single-cell imaging of two ERK activity reporters and expression of Fra-1, a target gene controlling epithelial cell identity. We find that Fra-1 is expressed in proportion to the amplitude and duration of ERK activity. In contrast to previous "persistence detector" and "selective filter" models in which Fra-1 expression only occurs when ERK activity persists beyond a threshold duration, our observations demonstrate that the network regulating Fra-1 expression integrates total ERK activity and responds to it linearly. However, exploration of a generalized mathematical model of the Fra-1 coherent feedforward loop demonstrates that it can perform either linear integration or persistence detection, depending on the basal mRNA production rate and protein production delays. Our data indicate that significant basal expression and short delays cause Fra-1 to respond linearly to integrated ERK activity.