Systems-wide Identification of cis-Regulatory Elements in Proteins

Widespread alteration of protein autoinhibition in human cancers

Autoinhibition is a prevalent allosteric regulatory mechanism in signaling proteins. Reduced autoinhibition underlies the tumorigenic effect of some known cancer drivers, but whether autoinhibition is altered generally in cancer remains elusive. Here, we demonstrate that cancer-associated missense mutations, in-frame insertions/deletions, and fusion breakpoints are enriched within inhibitory allosteric switches (IASs) across all cancer types. Selection for IASs that are recurrently mutated in cancers identifies established and unknown cancer drivers. Recurrent missense mutations in IASs of these drivers are associated with distinct, cancer-specific changes in molecular signaling. For the specific case of PPP3CA, the catalytic subunit of calcineurin, we provide insights into the molecular mechanisms of altered autoinhibition by cancer mutations using biomolecular simulations, and demonstrate that such mutations are associated with transcriptome changes consistent with increased calcineurin signaling. Our integrative study shows that autoinhibition-modulating genetic alterations are positively selected for by cancer cells.

Mechanistic and evolutionary insights into isoform-specific 'supercharging' in DCLK family kinases

Catalytic signaling outputs of protein kinases are dynamically regulated by an array of structural mechanisms, including allosteric interactions mediated by intrinsically disordered segments flanking the conserved catalytic domain. The doublecortin-like kinases (DCLKs) are a family of microtubule-associated proteins characterized by a flexible C-terminal autoregulatory 'tail' segment that varies in length across the various human DCLK isoforms. However, the mechanism whereby these isoform-specific variations contribute to unique modes of autoregulation is not well understood. Here, we employ a combination of statistical sequence analysis, molecular dynamics simulations, and in vitro mutational analysis to define hallmarks of DCLK family evolutionary divergence, including analysis of splice variants within the DCLK1 sub-family, which arise through alternative codon usage and serve to 'supercharge' the inhibitory potential of the DCLK1 C-tail. We identify co-conserved motifs that readily distinguish DCLKs from all other calcium calmodulin kinases (CAMKs), and a 'Swiss Army' assembly of distinct motifs that tether the C-terminal tail to conserved ATP and substrate-binding regions of the catalytic domain to generate a scaffold for autoregulation through C-tail dynamics. Consistently, deletions and mutations that alter C-terminal tail length or interfere with co-conserved interactions within the catalytic domain alter intrinsic protein stability, nucleotide/inhibitor binding, and catalytic activity, suggesting isoform-specific regulation of activity through alternative splicing. Our studies provide a detailed framework for investigating kinome-wide regulation of catalytic output through cis-regulatory events mediated by intrinsically disordered segments, opening new avenues for the design of mechanistically divergent DCLK1 modulators, stabilizers, or degraders.

Mechanistic and evolutionary insights into isoform-specific 'supercharging' in DCLK family kinases

Catalytic signaling outputs of protein kinases are dynamically regulated by an array of structural mechanisms, including allosteric interactions mediated by intrinsically disordered segments flanking the conserved catalytic domain. The Doublecortin Like Kinases (DCLKs) are a family of microtubule-associated proteins characterized by a flexible C-terminal autoregulatory 'tail' segment that varies in length across the various human DCLK isoforms. However, the mechanism whereby these isoform-specific variations contribute to unique modes of autoregulation is not well understood. Here, we employ a combination of statistical sequence analysis, molecular dynamics simulations and in vitro mutational analysis to define hallmarks of DCLK family evolutionary divergence, including analysis of splice variants within the DCLK1 sub-family, which arise through alternative codon usage and serve to 'supercharge' the inhibitory potential of the DCLK1 C-tail. We identify co-conserved motifs that readily distinguish DCLKs from all other Calcium Calmodulin Kinases (CAMKs), and a 'Swiss-army' assembly of distinct motifs that tether the C-terminal tail to conserved ATP and substrate-binding regions of the catalytic domain to generate a scaffold for auto-regulation through C-tail dynamics. Consistently, deletions and mutations that alter C-terminal tail length or interfere with co-conserved interactions within the catalytic domain alter intrinsic protein stability, nucleotide/inhibitor-binding, and catalytic activity, suggesting isoform-specific regulation of activity through alternative splicing. Our studies provide a detailed framework for investigating kinome-wide regulation of catalytic output through cis-regulatory events mediated by intrinsically disordered segments, opening new avenues for the design of mechanistically-divergent DCLK1 modulators, stabilizers or degraders.

Where protein structure and cell diversity meet

Protein-protein interaction networks - interactomes - are charted with the hope to understand how phenotypes emerge and how they are altered in disease states. Early efforts to map interactomes have focused on the assembly of context agnostic, reference networks. However, recent studies have mapped interactomes across different cell lines and tissues, finding highly variable interactomes due to the rewiring of protein-protein interactions in different contexts. Increasing evidence points to significant links between protein structure and interactome diversity seen across cell types and tissues. We discuss how recent findings support the key role of alternative splicing and phosphorylation, two well-established regulators of protein structural and functional diversity, in defining cell type- and tissue-specific interactomes. Moreover, we show that intrinsically disordered protein regions are most favorably equipped to support interactome rewiring by acting as hubs of protein structure and function regulation.

Connecting sequence features within the disordered C-terminal linker of Bacillus subtilis FtsZ to functions and bacterial cell division

Intrinsically disordered regions (IDRs) can function as autoregulators of folded enzymes to which they are tethered. One example is the bacterial cell division protein FtsZ. This includes a folded core and a C-terminal tail (CTT) that encompasses a poorly conserved, disordered C-terminal linker (CTL) and a well-conserved 17-residue C-terminal peptide (CT17). Sites for GTPase activity of FtsZs are formed at the interface between GTP binding sites and T7 loops on cores of adjacent subunits within dimers. Here, we explore the basis of autoregulatory functions of the CTT in Bacillus subtilis FtsZ (Bs-FtsZ). Molecular simulations show that the CT17 of Bs-FtsZ makes statistically significant CTL-mediated contacts with the T7 loop. Statistical coupling analysis of more than 1,000 sequences from FtsZ orthologs reveals clear covariation of the T7 loop and the CT17 with most of the core domain, whereas the CTL is under independent selection. Despite this, we discover the conservation of nonrandom sequence patterns within CTLs across orthologs. To test how the nonrandom patterns of CTLs mediate CTT-core interactions and modulate FtsZ functionalities, we designed Bs-FtsZ variants by altering the patterning of oppositely charged residues within the CTL. Such alterations disrupt the core-CTT interactions, lead to anomalous assembly and inefficient GTP hydrolysis in vitro and protein degradation, aberrant assembly, and disruption of cell division in vivo. Our findings suggest that viable CTLs in FtsZs are likely to be IDRs that encompass nonrandom, functionally relevant sequence patterns that also preserve three-way covariation of the CT17, the T7 loop, and core domain.

Tracking mutation and drug-driven alterations of oncokinase conformations

Numerous kinases act as central nodes of cellular signaling networks. As such, many of these enzymes function as molecular switches for coordinating spatiotemporal signal transmission. Typically, it is the compartmentalized phosphorylation of protein substrates which relays the transient input signal to determine decisive physiological cell responses. Genomic alterations affect kinase abundance and/or their activities which contribute to the malignant transformation, progression, and metastasis of human cancers. Thus, major drug discovery efforts have been made to identify lead molecules targeting clinically relevant oncokinases. The concept of personalized medicine aims to apply the therapeutic agent with the highest efficacy towards a patient-specific mutation. Here, we discuss the implementation of a cell-based reporter system which may foster the decision-making process to identify the most promising lead-molecules. We present a modular kinase conformation (KinCon) biosensor platform for live-cell analyses of kinase activity states. This biosensor facilitates the recording of kinase activity conformations of the wild-type and the respective mutated kinase upon lead molecule exposure. We reflect proof-of-principle studies demonstrating how this technology has been extended to profile drug properties of the full-length kinases BRAF and MEK1 in intact cells. Further, we pinpoint how this technology may open new avenues for systematic and patient-tailored drug discovery efforts. Overall, this precision-medicineoriented biosensor concept aims to determine kinase inhibitor specificity and anticipate their drug efficacies.

Transient Tertiary Structures of Disordered Dynein Intermediate Chain Regulate its Interactions with Multiple Partners

The N-terminal domain of dynein intermediate chain (N-IC) is central to the cytoplasmic dynein 'cargo attachment subcomplex' and regulation of motor activity. It is a prototypical intrinsically disordered protein (IDP), serving as a primarily disordered polybivalent molecular scaffold for numerous binding partners, including three dimeric dynein light chains and coiled coil domains of dynein partners dynactin p150^Glued and NudE. At the very N-terminus, a 40 amino acid single alpha helix (SAH) forms the major binding site for both p150^Glued and NudE, while a shorter nascent helix (H2) separated from SAH by a disordered linker, is necessary for tight binding to dynactin p150^Glued but not to NudE. Here we demonstrate that transient tertiary interactions in this highly dynamic protein underlie the differences in its interactions with p150^Glued and NudE. NMR paramagnetic relaxation enhancement experiments and restrained molecular dynamics simulations identify interactions between the two non-contiguous SAH and H2 helical regions, the extent of which correlates with the length and stability of H2, showing clearly that tertiary and secondary structure formation are coupled in IDPs. These interactions are significantly attenuated when N-IC is bound to NudE, suggesting that NudE binding shifts the conformational ensemble to one that is more extended and with less structure in H2. While the intrinsic disorder and flexibility in N-IC modulate its ability to serve as a binding platform for numerous partners, deviations of this protein from random-coil behavior provide a process for regulating these binding interactions and potentially the dynein motor.

Mutation-oriented profiling of autoinhibitory kinase conformations predicts RAF inhibitor efficacies

Kinase-targeted therapies have the potential to improve the survival of patients with cancer. However, the cancer-specific spectrum of kinase alterations exhibits distinct functional properties and requires mutation-oriented drug treatments. Besides post-translational modifications and diverse intermolecular interactions of kinases, it is the distinct disease mutation which reshapes full-length kinase conformations, affecting their activity. Oncokinase mutation profiles differ between cancer types, as it was shown for BRAF in melanoma and non-small-cell lung cancers. Here, we present the target-oriented application of a kinase conformation (KinCon) reporter platform for live-cell measurements of autoinhibitory kinase activity states. The bioluminescence-based KinCon biosensor allows the tracking of conformation dynamics of full-length kinases in intact cells and real time. We show that the most frequent BRAF cancer mutations affect kinase conformations and thus the engagement and efficacy of V600E-specific BRAF inhibitors (BRAFi). We illustrate that the patient mutation harboring KinCon reporters display differences in the effectiveness of the three clinically approved BRAFi vemurafenib, encorafenib, and dabrafenib and the preclinical paradox breaker PLX8394. We confirmed KinCon-based drug efficacy predictions for BRAF mutations other than V600E in proliferation assays using patient-derived lung cancer cell lines and by analyzing downstream kinase signaling. The systematic implementation of such conformation reporters will allow to accelerate the decision process for the mutation-oriented RAF-kinase cancer therapy. Moreover, we illustrate that the presented kinase reporter concept can be extended to other kinases which harbor patient mutations. Overall, KinCon profiling provides additional mechanistic insights into full-length kinase functions by reporting protein-protein interaction (PPI)-dependent, mutation-specific, and drug-driven changes of kinase activity conformations.

KinCon: Cell-based recording of full-length kinase conformations

The spectrum of kinase alterations displays distinct functional characteristics and requires kinase mutation-oriented strategies for therapeutic interference. Besides phosphotransferase activity, protein abundance, and intermolecular interactions, particular patient-mutations promote pathological kinase conformations. Despite major advances in identifying lead molecules targeting clinically relevant oncokinase functions, still many kinases are neglected and not part of drug discovery efforts. One explanation is attributed to challenges in tracking kinase activities. Chemical probes are needed to functionally annotate kinase functions, whose activities may not always depend on catalyzing phospho-transfer. Such non-catalytic kinase functions are related to transitions of full-length kinase conformations. Recent findings underline that cell-based reporter systems can be adapted to record conformation changes of kinases. Here, we discuss the possible applications of an extendable kinase conformation (KinCon) reporter toolbox for live-cell recording of kinase states. KinCon is a genetically encoded bioluminescence-based biosensor platform, which can be subjected for measurements of conformation dynamics of mutated kinases upon small molecule inhibitor exposure. We hypothesize that such biosensors can be utilized to delineate the molecular modus operandi for kinase and pseudokinase regulation. This should pave the path for full-length kinase-targeted drug discovery efforts aiming to identify single and combinatory kinase inhibitor therapies with increased specificity and efficacy.

Intrinsically disordered domains: Sequence ➔ disorder ➔ function relationships

Disordered domains are long regions of intrinsic disorder that ideally have conserved sequences, conserved disorder, and conserved functions. These domains were first noticed in protein-protein interactions that are distinct from the interactions between two structured domains and the interactions between structured domains and linear motifs or molecular recognition features (MoRFs). So far, disordered domains have not been systematically characterized. Here, we present a bioinformatics investigation of the sequence-disorder-function relationships for a set of probable disordered domains (PDDs) identified from the Pfam database. All the Pfam seed proteins from those domains with at least one PDD sequence were collected. Most often, if a set contains one PDD sequence, then all members of the set are PDDs or nearly so. However, many seed sets have sequence collections that exhibit diverse proportions of predicted disorder and structure, thus giving the completely unexpected result that conserved sequences can vary substantially in predicted disorder and structure. In addition to the induction of structure by binding to protein partners, disordered domains are also induced to form structure by disulfide bond formation, by ion binding, and by complex formation with RNA or DNA. The two new findings, (a) that conserved sequences can vary substantially in their predicted disorder content and (b) that homologues from a single domain can evolve from structure to disorder (or vice versa), enrich our understanding of the sequence ➔ disorder ensemble ➔ function paradigm.

BRAF inhibitors promote intermediate BRAF(V600E) conformations and binary interactions with activated RAS

Oncogenic BRAF mutations initiate tumor formation by unleashing the autoinhibited kinase conformation and promoting RAS-decoupled proliferative RAF-MEK-ERK signaling. We have engineered luciferase-based biosensors to systematically track full-length BRAF conformations and interactions affected by tumorigenic kinase mutations and GTP loading of RAS. Binding of structurally diverse αC-helix-OUT BRAF inhibitors (BRAFi) showed differences in specificity and efficacy by shifting patient mutation-containing BRAF reporters from the definitive opened to more closed conformations. Unexpectedly, BRAFi engagement with the catalytic pocket of V600E-mutated BRAF stabilized an intermediate and inactive kinase conformation that enhanced binary RAS:RAF interactions, also independently of RAF dimerization in melanoma cells. We present evidence that the interference with RAS interactions and nanoclustering antagonizes the sequential formation of drug-induced RAS:RAF tetramers. This suggests a previously unappreciated allosteric effect of anticancer drug-driven intramolecular communication between the kinase and RAS-binding domains of mutated BRAF, which may further promote paradoxical kinase activation and drug resistance mechanisms.

In silico prediction of potential chemical reactions mediated by human enzymes

Background

Administered drugs are often converted into an ineffective or activated form by enzymes in our body. Conventional in silico prediction approaches focused on therapeutically important enzymes such as CYP450. However, there are more than thousands of different cellular enzymes that potentially convert administered drug into other forms.

Result

We developed an in silico model to predict which of human enzymes including metabolic enzymes as well as CYP450 family can catalyze a given chemical compound. The prediction is based on the chemical and physical similarity between known enzyme substrates and a query chemical compound. Our in silico model was developed using multiple linear regression and the model showed high performance (AUC = 0.896) despite of the large number of enzymes. When evaluated on a test dataset, it also showed significantly high performance (AUC = 0.746). Interestingly, evaluation with literature data showed that our model can be used to predict not only enzymatic reactions but also drug conversion and enzyme inhibition.

Conclusion

Our model was able to predict enzymatic reactions of a query molecule with a high accuracy. This may foster to discover new metabolic routes and to accelerate the computational development of drug candidates by enabling the prediction of the potential conversion of administered drugs into active or inactive forms.

Biophysical Characterization of the Tandem FHA Domain Regulatory Module from the Mycobacterium tuberculosis ABC Transporter Rv1747

The Mycobacterium tuberculosis ATP-binding cassette transporter Rv1747 is a putative exporter of cell wall biosynthesis intermediates. Rv1747 has a cytoplasmic regulatory module consisting of two pThr-interacting Forkhead-associated (FHA) domains connected by a conformationally disordered linker with two phospho-acceptor threonines (pThr). The structures of FHA-1 and FHA-2 were determined by X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy, respectively. Relative to the canonical 11-strand β-sandwich FHA domain fold of FHA-1, FHA-2 is circularly permuted and lacking one β-strand. Nevertheless, the two share a conserved pThr-binding cleft. FHA-2 is less stable and more dynamic than FHA-1, yet binds model pThr peptides with moderately higher affinity (∼50 μM versus 500 μM equilibrium dissociation constants). Based on NMR relaxation and chemical shift perturbation measurements, when joined within a polypeptide chain, either FHA domain can bind either linker pThr to form intra- and intermolecular complexes. We hypothesize that this enables tunable phosphorylation-dependent multimerization to regulate Rv1747 transporter activity.

Regulation of apoptosis by an intrinsically disordered region of Bcl-xL

Intrinsically disordered regions (IDRs) of proteins often regulate function upon posttranslational modifications (PTMs) through interactions with folded domains. An IDR linking two α-helices (α1–α2) of the anti-apoptotic protein, Bcl-xL, experiences several PTMs, which reduce anti-apoptotic activity. Here, we report that PTMs within the α1–α2 IDR promote its interaction with the folded core of Bcl-xL that inhibits the pro-apoptotic activity of two types of regulatory targets, BH3-only proteins and p53. This autoregulation utilizes an allosteric pathway where, in one direction, the IDR induces a direct displacement of p53 from Bcl-xL coupled to allosteric displacement of simultaneously bound BH3-only partners. This pathway operates in the opposite direction when the BH3-only protein PUMA binds to the BH3 binding groove of Bcl-xL, directly displacing other bound BH3-only proteins, and allosterically remodeling the distal site, displacing p53. Our findings show how an IDR enhances functional versatility through PTM-dependent, allosteric regulation of a folded protein domain.

Sequence fingerprints distinguish erroneous from correct predictions of intrinsically disordered protein regions

More than 60 prediction methods for intrinsically disordered proteins (IDPs) have been developed over the years, many of which are accessible on the World Wide Web. Nearly, all of these predictors give balanced accuracies in the ~65%-~80% range. Since predictors are not perfect, further studies are required to uncover the role of amino acid residues in native IDP as compared to predicted IDP regions. In the present work, we make use of sequences of 100% predicted IDP regions, false positive disorder predictions, and experimentally determined IDP regions to distinguish the characteristics of native versus predicted IDP regions. A higher occurrence of asparagine is observed in sequences of native IDP regions but not in sequences of false positive predictions of IDP regions. The occurrences of certain combinations of amino acids at the pentapeptide level provide a distinguishing feature in the IDPs with respect to globular proteins. The distinguishing features presented in this paper provide insights into the sequence fingerprints of amino acid residues in experimentally determined as compared to predicted IDP regions. These observations and additional work along these lines should enable the development of improvements in the accuracy of disorder prediction algorithm.

Hunting for Cis-Regulatory Elements in Proteins

A new bioinformatics tool predicts natively disordered protein control elements that function in cis, opening the door to more systematic studies of this biomedically important class of protein modules.