Accurate proteome-wide missense variant effect prediction with AlphaMissense.

The vast majority of missense variants observed in the human genome are of unknown clinical significance. We present AlphaMissense, an adaptation of AlphaFold fine-tuned on human and primate variant population frequency databases to predict missense variant pathogenicity. By combining structural context and evolutionary conservation, our model achieves state-of-the-art results across a wide range of genetic and experimental benchmarks, all without explicitly training on such data. The average pathogenicity score of genes is also predictive for their cell essentiality, capable of identifying short essential genes that existing statistical approaches are underpowered to detect. As a resource to the community, we provide a database of predictions for all possible human single amino acid substitutions and classify 89% of missense variants as either likely benign or likely pathogenic.

BioPlexR and BioPlexPy: integrated data products for the analysis of human protein interactions.

SUMMARY: The BioPlex project has created two proteome scale, cell-line-specific protein-protein interaction (PPI) networks: the first in 293T cells, including 120k interactions among 15k proteins; and the second in HCT116 cells, including 70k interactions between 10k proteins. Here, we describe programmatic access to the BioPlex PPI networks and integration with related resources from within R and Python. Besides PPI networks for 293T and HCT116 cells, this includes access to CORUM protein complex data, PFAM protein domain data, PDB protein structures, and transcriptome and proteome data for the two cell lines. The implemented functionality serves as a basis for integrative downstream analysis of BioPlex PPI data with domain-specific R and Python packages, including efficient execution of maximum scoring sub-network analysis, protein domain-domain association analysis, mapping of PPIs onto 3D protein structures and analysis of BioPlex PPIs at the interface of transcriptomic and proteomic data. AVAILABILITY AND IMPLEMENTATION: The BioPlex R package is available from Bioconductor (bioconductor.org/packages/BioPlex), and the BioPlex Python package is available from PyPI (pypi.org/project/bioplexpy). Applications and downstream analyses are available from GitHub (github.com/ccb-hms/BioPlexAnalysis).

Structure-function analysis of the SHOC2-MRAS-PP1C holophosphatase complex.

Receptor tyrosine kinase (RTK)-RAS signalling through the downstream mitogen-activated protein kinase (MAPK) cascade regulates cell proliferation and survival. The SHOC2-MRAS-PP1C holophosphatase complex functions as a key regulator of RTK-RAS signalling by removing an inhibitory phosphorylation event on the RAF family of proteins to potentiate MAPK signalling1. SHOC2 forms a ternary complex with MRAS and PP1C, and human germline gain-of-function mutations in this complex result in congenital RASopathy syndromes2-5. However, the structure and assembly of this complex are poorly understood. Here we use cryo-electron microscopy to resolve the structure of the SHOC2-MRAS-PP1C complex. We define the biophysical principles of holoenzyme interactions, elucidate the assembly order of the complex, and systematically interrogate the functional consequence of nearly all of the possible missense variants of SHOC2 through deep mutational scanning. We show that SHOC2 binds PP1C and MRAS through the concave surface of the leucine-rich repeat region and further engages PP1C through the N-terminal disordered region that contains a cryptic RVXF motif. Complex formation is initially mediated by interactions between SHOC2 and PP1C and is stabilized by the binding of GTP-loaded MRAS. These observations explain how mutant versions of SHOC2 in RASopathies and cancer stabilize the interactions of complex members to enhance holophosphatase activity. Together, this integrative structure-function model comprehensively defines key binding interactions within the SHOC2-MRAS-PP1C holophosphatase complex and will inform therapeutic development .

Sparse dictionary learning recovers pleiotropy from human cell fitness screens.

In high-throughput functional genomic screens, each gene product is commonly assumed to exhibit a singular biological function within a defined protein complex or pathway. In practice, a single gene perturbation may induce multiple cascading functional outcomes, a genetic principle known as pleiotropy. Here, we model pleiotropy in fitness screen collections by representing each gene perturbation as the sum of multiple perturbations of biological functions, each harboring independent fitness effects inferred empirically from the data. Our approach (Webster) recovered pleiotropic functions for DNA damage proteins from genotoxic fitness screens, untangled distinct signaling pathways upstream of shared effector proteins from cancer cell fitness screens, and predicted the stoichiometry of an unknown protein complex subunit from fitness data alone. Modeling compound sensitivity profiles in terms of genetic functions recovered compound mechanisms of action. Our approach establishes a sparse approximation mechanism for unraveling complex genetic architectures underlying high-dimensional gene perturbation readouts.

OCT4 cooperates with distinct ATP-dependent chromatin remodelers in naïve and primed pluripotent states in human.

Understanding the molecular underpinnings of pluripotency is a prerequisite for optimal maintenance and application of embryonic stem cells (ESCs). While the protein-protein interactions of core pluripotency factors have been identified in mouse ESCs, their interactome in human ESCs (hESCs) has not to date been explored. Here we mapped the OCT4 interactomes in naïve and primed hESCs, revealing extensive connections to mammalian ATP-dependent nucleosome remodeling complexes. In naïve hESCs, OCT4 is associated with both BRG1 and BRM, the two paralog ATPases of the BAF complex. Genome-wide location analyses and genetic studies reveal that these two enzymes cooperate in a functionally redundant manner in the transcriptional regulation of blastocyst-specific genes. In contrast, in primed hESCs, OCT4 cooperates with BRG1 and SOX2 to promote chromatin accessibility at ectodermal genes. This work reveals how a common transcription factor utilizes differential BAF complexes to control distinct transcriptional programs in naïve and primed hESCs.

Dual proteome-scale networks reveal cell-specific remodeling of the human interactome.

Thousands of interactions assemble proteins into modules that impart spatial and functional organization to the cellular proteome. Through affinity-purification mass spectrometry, we have created two proteome-scale, cell-line-specific interaction networks. The first, BioPlex 3.0, results from affinity purification of 10,128 human proteins-half the proteome-in 293T cells and includes 118,162 interactions among 14,586 proteins. The second results from 5,522 immunoprecipitations in HCT116 cells. These networks model the interactome whose structure encodes protein function, localization, and complex membership. Comparison across cell lines validates thousands of interactions and reveals extensive customization. Whereas shared interactions reside in core complexes and involve essential proteins, cell-specific interactions link these complexes, "rewiring" subnetworks within each cell's interactome. Interactions covary among proteins of shared function as the proteome remodels to produce each cell's phenotype. Viewable interactively online through BioPlexExplorer, these networks define principles of proteome organization and enable unknown protein characterization.

Synthetic Lethal Interaction of SHOC2 Depletion with MEK Inhibition in RAS-Driven Cancers.

The mitogen-activated protein kinase (MAPK) pathway is a critical effector of oncogenic RAS signaling, and MAPK pathway inhibition may be an effective combination treatment strategy. We performed genome-scale loss-of-function CRISPR-Cas9 screens in the presence of a MEK1/2 inhibitor (MEKi) in KRAS-mutant pancreatic and lung cancer cell lines and identified genes that cooperate with MEK inhibition. While we observed heterogeneity in genetic modifiers of MEKi sensitivity across cell lines, several recurrent classes of synthetic lethal vulnerabilities emerged at the pathway level. Multiple members of receptor tyrosine kinase (RTK)-RAS-MAPK pathways scored as sensitizers to MEKi. In particular, we demonstrate that knockout, suppression, or degradation of SHOC2, a positive regulator of MAPK signaling, specifically cooperated with MEK inhibition to impair proliferation in RAS-driven cancer cells. The depletion of SHOC2 disrupted survival pathways triggered by feedback RTK signaling in response to MEK inhibition. Thus, these findings nominate SHOC2 as a potential target for combination therapy.

The ATPase module of mammalian SWI/SNF family complexes mediates subcomplex identity and catalytic activity-independent genomic targeting.

Perturbations to mammalian switch/sucrose non-fermentable (mSWI/SNF) chromatin remodeling complexes have been widely implicated as driving events in cancer1. One such perturbation is the dual loss of the SMARCA4 and SMARCA2 ATPase subunits in small cell carcinoma of the ovary, hypercalcemic type (SCCOHT)2-5, SMARCA4-deficient thoracic sarcomas6 and dedifferentiated endometrial carcinomas7. However, the consequences of dual ATPase subunit loss on mSWI/SNF complex subunit composition, chromatin targeting, DNA accessibility and gene expression remain unknown. Here we identify an ATPase module of subunits that is required for functional specification of the Brahma-related gene-associated factor (BAF) and polybromo-associated BAF (PBAF) mSWI/SNF family subcomplexes. Using SMARCA4/2 ATPase mutant variants, we define the catalytic activity-dependent and catalytic activity-independent contributions of the ATPase module to the targeting of BAF and PBAF complexes on chromatin genome-wide. Finally, by linking distinct mSWI/SNF complex target sites to tumor-suppressive gene expression programs, we clarify the transcriptional consequences of SMARCA4/2 dual loss in SCCOHT.

A non-canonical SWI/SNF complex is a synthetic lethal target in cancers driven by BAF complex perturbation.

Mammalian SWI/SNF chromatin remodelling complexes exist in three distinct, final-form assemblies: canonical BAF (cBAF), PBAF and a newly characterized non-canonical complex (ncBAF). However, their complex-specific targeting on chromatin, functions and roles in disease remain largely undefined. Here, we comprehensively mapped complex assemblies on chromatin and found that ncBAF complexes uniquely localize to CTCF sites and promoters. We identified ncBAF subunits as synthetic lethal targets specific to synovial sarcoma and malignant rhabdoid tumours, which both exhibit cBAF complex (SMARCB1 subunit) perturbation. Chemical and biological depletion of the ncBAF subunit, BRD9, rapidly attenuates synovial sarcoma and malignant rhabdoid tumour cell proliferation. Importantly, in cBAF-perturbed cancers, ncBAF complexes maintain gene expression at retained CTCF-promoter sites and function in a manner distinct from fusion oncoprotein-bound complexes. Together, these findings unmask the unique targeting and functional roles of ncBAF complexes and present new cancer-specific therapeutic targets.

Targeting CDK9 Reactivates Epigenetically Silenced Genes in Cancer.

Cyclin-dependent kinase 9 (CDK9) promotes transcriptional elongation through RNAPII pause release. We now report that CDK9 is also essential for maintaining gene silencing at heterochromatic loci. Through a live cell drug screen with genetic confirmation, we discovered that CDK9 inhibition reactivates epigenetically silenced genes in cancer, leading to restored tumor suppressor gene expression, cell differentiation, and activation of endogenous retrovirus genes. CDK9 inhibition dephosphorylates the SWI/SNF protein BRG1, which contributes to gene reactivation. By optimization through gene expression, we developed a highly selective CDK9 inhibitor (MC180295, IC50 = 5 nM) that has broad anti-cancer activity in vitro and is effective in in vivo cancer models. Additionally, CDK9 inhibition sensitizes to the immune checkpoint inhibitor α-PD-1 in vivo, making it an excellent target for epigenetic therapy of cancer.

Modular Organization and Assembly of SWI/SNF Family Chromatin Remodeling Complexes.

Mammalian SWI/SNF (mSWI/SNF) ATP-dependent chromatin remodeling complexes are multi-subunit molecular machines that play vital roles in regulating genomic architecture and are frequently disrupted in human cancer and developmental disorders. To date, the modular organization and pathways of assembly of these chromatin regulators remain unknown, presenting a major barrier to structural and functional determination. Here, we elucidate the architecture and assembly pathway across three classes of mSWI/SNF complexes-canonical BRG1/BRM-associated factor (BAF), polybromo-associated BAF (PBAF), and newly defined ncBAF complexes-and define the requirement of each subunit for complex formation and stability. Using affinity purification of endogenous complexes from mammalian and Drosophila cells coupled with cross-linking mass spectrometry (CX-MS) and mutagenesis, we uncover three distinct and evolutionarily conserved modules, their organization, and the temporal incorporation of these modules into each complete mSWI/SNF complex class. Finally, we map human disease-associated mutations within subunits and modules, defining specific topological regions that are affected upon subunit perturbation.

Binding of TMPRSS2-ERG to BAF Chromatin Remodeling Complexes Mediates Prostate Oncogenesis.

Chromosomal rearrangements resulting in the fusion of TMPRSS2, an androgen-regulated gene, and the ETS family transcription factor ERG occur in over half of prostate cancers. However, the mechanism by which ERG promotes oncogenic gene expression and proliferation remains incompletely understood. Here, we identify a binding interaction between ERG and the mammalian SWI/SNF (BAF) ATP-dependent chromatin remodeling complex, which is conserved among other oncogenic ETS factors, including ETV1, ETV4, and ETV5. We find that ERG drives genome-wide retargeting of BAF complexes in a manner dependent on binding of ERG to the ETS DNA motif. Moreover, ERG requires intact BAF complexes for chromatin occupancy and BAF complex ATPase activity for target gene regulation. In a prostate organoid model, BAF complexes are required for ERG-mediated basal-to-luminal transition, a hallmark of ERG activity in prostate cancer. These observations suggest a fundamental interdependence between ETS transcription factors and BAF chromatin remodeling complexes in cancer.

Interrogation of Mammalian Protein Complex Structure, Function, and Membership Using Genome-Scale Fitness Screens.

Protein complexes are assemblies of subunits that have co-evolved to execute one or many coordinated functions in the cellular environment. Functional annotation of mammalian protein complexes is critical to understanding biological processes, as well as disease mechanisms. Here, we used genetic co-essentiality derived from genome-scale RNAi- and CRISPR-Cas9-based fitness screens performed across hundreds of human cancer cell lines to assign measures of functional similarity. From these measures, we systematically built and characterized functional similarity networks that recapitulate known structural and functional features of well-studied protein complexes and resolve novel functional modules within complexes lacking structural resolution, such as the mammalian SWI/SNF complex. Finally, by integrating functional networks with large protein-protein interaction networks, we discovered novel protein complexes involving recently evolved genes of unknown function. Taken together, these findings demonstrate the utility of genetic perturbation screens alone, and in combination with large-scale biophysical data, to enhance our understanding of mammalian protein complexes in normal and disease states.

IL-1ß Inhibits Connexin 43 and Disrupts Decidualization of Human Endometrial Stromal Cells Through ERK1/2 and p38 MAP Kinase.

Inflammation can interfere with endometrial receptivity. We examined how interleukin 1ß (IL-1ß) affects expression of the uterine gap junction protein connexin 43 (Cx43), which is known to be critical for embryonic implantation. We used an in vitro model of human endometrial stromal cells (ESCs), Western blotting, and a combination of validated, selective kinase inhibitors to evaluate five canonical IL-1ß signaling pathways. Cx43 and two other markers of ESC differentiation (prolactin and VEGF) were inhibited predominantly via IL-1ß-activated ERK1/2 and p38 MAP kinase cascades. The findings were corroborated using small interfering RNA to silence critical genes in either pathway. By contrast, upregulation of endogenous pro-IL-1α and pro-IL-1ß following recombinant IL-1ß treatment was mediated via the Jun N-terminal kinase pathway. The clinicopharmacological significance of our findings is that multiple signaling cascades may need to be neutralized to reverse deleterious effects of IL-1ß on human endometrial function.

Combined Spinal Epidural Technique for Labor Analgesia Does Not Delay Recognition of Epidural Catheter Failures: A Single-center Retrospective Cohort Survival Analysis.

BACKGROUND: It is unclear whether recognition of epidural catheter failures is delayed with combined spinal epidural technique (CSE) compared to traditional epidural technique (EPID) when used for labor analgesia. The authors hypothesized that recognition of failed catheters is not delayed by CSE. METHODS: Anesthetic, obstetric, and quality assurance records from 2,395 labor neuraxial procedures (1,440 CSE and 955 EPID) performed at Forsyth Medical Center (Winston-Salem, North Carolina) between June 30 and December 31, 2012, were retrospectively analyzed. The primary outcome was catheter survival (failure-free) time during labor analgesia. A proportional hazards model with the counting method was used to assess relationships between the techniques and survival (failure-free) time of catheters, while controlling for subjects' body mass index and providers' level of training in the final best-fit multivariable regression model. RESULTS: Cumulative incidence of epidural catheter failures was 6.6% for CSE and 11.6% for EPID (P = 0.001). In the multivariable regression model, catheters placed with CSE versus epidural were less likely to fail (hazard ratio, 0.58; 95% CI, 0.43 to 0.79; P = 0.0002) for labor analgesia. Among the catheters that failed, there was no overall difference in failure time course between the techniques (hazard ratio, 1.17; 95% CI, 0.89 to 1.54; P = 0.26) even though more failed catheters with CSE (48.4%) than with EPID (30.6%) were recognized within the first 30 min of placement (P = 0.009). CONCLUSIONS: In this cohort, CSE has a significantly lower risk of overall epidural catheter failures than EPID and does not delay recognition of epidural catheter failures. Choice of CSE versus EPID should be based on overall risk of failure, efficacy, and side effects.

CRISPR/Cas9-Directed Genome Editing of Cultured Cells.

Human genome engineering has been transformed by the introduction of the CRISPR (clustered regularly interspaced short palindromic repeats)/Cas (CRISPR-associated) system found in most bacteria and archaea. Type II CRISPR/Cas systems have been engineered to induce RNA-guided genome editing in human cells, where small RNAs function together with Cas9 nucleases for sequence-specific cleavage of target sequences. Here we describe the protocol for Cas9-mediated human genome engineering, including construct building and transfection methods necessary for delivering Cas9 and guide RNA (gRNA) into human-induced pluripotent stem cells (hiPSCs) and HEK293 cells. Following genome editing, we also describe methods to assess genome editing efficiency using next-generation sequencing and isolate monoclonal hiPSCs with the desired modifications for downstream applications.