Common biochemical properties of metabolic genes recurrently dysregulated in tumors

BACKGROUND

Tumor initiation and progression are associated with numerous metabolic alterations. However, the biochemical drivers and constraints that contribute to metabolic gene dysregulation are unclear.

METHODS

Here, we present MetOncoFit, a computational model that integrates 142 metabolic features that can impact tumor fitness, including enzyme catalytic activity, pathway association, network topology, and reaction flux. MetOncoFit uses genome-scale metabolic modeling and machine-learning to quantify the relative importance of various metabolic features in predicting cancer metabolic gene expression, copy number variation, and survival data.

RESULTS

Using MetOncoFit, we performed a meta-analysis of 9 cancer types and over 4500 samples from TCGA, Prognoscan, and COSMIC tumor databases. MetOncoFit accurately predicted enzyme differential expression and its impact on patient survival using the 142 attributes of metabolic enzymes. Our analysis revealed that enzymes with high catalytic activity were frequently upregulated in many tumors and associated with poor survival. Topological analysis also identified specific metabolites that were hot spots of dysregulation.

CONCLUSIONS

MetOncoFit integrates a broad range of datasets to understand how biochemical and topological features influence metabolic gene dysregulation across various cancer types. MetOncoFit was able to achieve significantly higher accuracy in predicting differential expression, copy number variation, and patient survival than traditional modeling approaches. Overall, MetOncoFit illuminates how enzyme activity and metabolic network architecture influences tumorigenesis.

The entropic force generated by intrinsically disordered segments tunes protein function

Protein structures are dynamic and can explore a large conformational landscape^1,2. Only some of these structural substates are important for protein function (such as ligand binding, catalysis and regulation)^3-5. How evolution shapes the structural ensemble to optimize a specific function is poorly understood^3,4. One of the constraints on the evolution of proteins is the stability of the folded 'native' state. Despite this, 44% of the human proteome contains intrinsically disordered peptide segments greater than 30 residues in length⁶, the majority of which have no known function^7-9. Here we show that the entropic force produced by an intrinsically disordered carboxy terminus (ID-tail) shifts the conformational ensemble of human UDP-α-D-glucose-6-dehydrogenase (UGDH) towards a substate with a high affinity for an allosteric inhibitor. The function of the ID-tail does not depend on its sequence or chemical composition. Instead, the affinity enhancement can be accurately predicted based on the length of the intrinsically disordered segment, and is consistent with the entropic force generated by an unstructured peptide attached to the protein surface^10-13. Our data show that the unfolded state of the ID-tail rectifies the dynamics and structure of UGDH to favour inhibitor binding. Because this entropic rectifier does not have any sequence or structural constraints, it is an easily acquired adaptation. This model implies that evolution selects for disordered segments to tune the energy landscape of proteins, which may explain the persistence of intrinsic disorder in the proteome.

Identification and classification of small molecule kinases: insights into substrate recognition and specificity

BACKGROUND

Many prokaryotic kinases that phosphorylate small molecule substrates, such as antibiotics, lipids and sugars, are evolutionarily related to Eukaryotic Protein Kinases (EPKs). These Eukaryotic-Like Kinases (ELKs) share the same overall structural fold as EPKs, but differ in their modes of regulation, substrate recognition and specificity-the sequence and structural determinants of which are poorly understood.

RESULTS

To better understand the basis for ELK specificity, we applied a Bayesian classification procedure designed to identify sequence determinants responsible for functional divergence. This reveals that a large and diverse family of aminoglycoside kinases, characterized members of which are involved in antibiotic resistance, fall into major sub-groups based on differences in putative substrate recognition motifs. Aminoglycoside kinase substrate specificity follows simple rules of alternating hydroxyl and amino groups that is strongly correlated with variations at the DFG + 1 position.

CONCLUSIONS

Substrate specificity determining features in small molecule kinases are mostly confined to the catalytic core and can be identified based on quantitative sequence and crystal structure comparisons.

Inhibiting EGFR dimerization using triazolyl-bridged dimerization arm mimics

The epidermal growth factor receptor (EGFR) is overexpressed in multiple carcinomas and is the focus of a variety of targeted therapies. Here we report the design of peptide-based compounds that mimic the EGFR dimerization arm and inhibit allosteric activation of EGFR. These peptides are modified to contain a triazolyl bridge between the peptide strands to constrain the EGFR dimerization arm β-loop. In this study, we demonstrate that these peptides have significantly improved proteolytic stability over the non-modified peptide sequence, and their inhibitory effects are dependent on the number of the methylene units and orientation of the introduced triazolyl bridge. We identified a peptide, EDA2, which downregulates receptor phosphorylation and dimerization and reduces cell viability. This is the first example of a biologically active triazolyl-bridged peptide targeting the EGFR dimerization interface that effectively downregulates EGFR activation.

Mitochondrial ADCK3 employs an atypical protein kinase-like fold to enable coenzyme Q biosynthesis

The ancient UbiB protein kinase-like family is involved in isoprenoid lipid biosynthesis and is implicated in human diseases, but demonstration of UbiB kinase activity has remained elusive for unknown reasons. Here, we quantitatively define UbiB-specific sequence motifs and reveal their positions within the crystal structure of a UbiB protein, ADCK3. We find that multiple UbiB-specific features are poised to inhibit protein kinase activity, including an N-terminal domain that occupies the typical substrate binding pocket and a unique A-rich loop that limits ATP binding by establishing an unusual selectivity for ADP. A single alanine-to-glycine mutation of this loop flips this coenzyme selectivity and enables autophosphorylation but inhibits coenzyme Q biosynthesis in vivo, demonstrating functional relevance for this unique feature. Our work provides mechanistic insight into UbiB enzyme activity and establishes a molecular foundation for further investigation of how UbiB family proteins affect diseases and diverse biological pathways.

Deciphering the structural basis of eukaryotic protein kinase regulation

Biochemical and structural analysis of two features of kinase structure, the “R-spine” and “Shell,” afford a detailed insight into the regulation of eukaryotic protein kinases.

Eukaryotic protein kinases (EPKs) regulate numerous signaling processes by phosphorylating targeted substrates through the highly conserved catalytic domain. Our previous computational studies proposed a model stating that a properly assembled nonlinear motif termed the Regulatory (R) spine is essential for catalytic activity of EPKs. Here we define the required intramolecular interactions and biochemical properties of the R-spine and the newly identified “Shell” that surrounds the R-spine using site-directed mutagenesis and various in vitro phosphoryl transfer assays using cyclic AMP-dependent protein kinase as a representative of the entire kinome. Analysis of the 172 available Apo EPK structures in the protein data bank (PDB) revealed four unique structural conformations of the R-spine that correspond with catalytic inactivation of various EPKs. Elucidating the molecular entities required for the catalytic activation of EPKs and the identification of these inactive conformations opens new avenues for the design of efficient therapeutic EPK inhibitors.

Eukaryotic protein kinases (EPKs) have a highly conserved enzymatic kinase core that is involved in the regulation of numerous cell signaling processes through the transfer of a phosphate group from adenosine triphosphate (ATP) to more than 30% of human proteins. EPKs have been implicated in numerous human diseases, including cancer, cardiovascular diseases, and diabetes, making them one of the most sought-after therapeutic drug targets. The lack of structural diversity of the active kinase core has created a bottle-neck for designing successful therapeutic inhibitors. Here we describe the intramolecular interactions required for differentiating between the active and inactive states of EPKs. Kinases contain a hydrophobic regulatory spine (“R-spine”) that is disassembled in inactive kinases, and here we define an additional hydrophobic “Shell” that surrounds one end of the R-spine. Biochemical analysis of the five nonconsecutive R-spine residues and three nonconsecutive Shell residues shows that proper assembly of the R-spine and Shell is essential for maintaining kinase activity. Structural analysis of the 172 known structures of EPKs without bound ligands led to the identification of four inactive conformations that correlate with the disassembly of the R-spine. Understanding the molecular elements involved in the regulation of kinase activity and the identification of these diverse groups of inactive conformations should aid the design of more specific therapeutic EPK inhibitors.

Design principles underpinning the regulatory diversity of protein kinases

Protein phosphorylation in eukaryotes is carried out by a large and diverse family of protein kinases, which display remarkable diversity and complexity in their modes of regulation. The complex modes of regulation have evolved as a consequence of natural selection operating on protein kinase sequences for billions of years. Here we describe how quantitative comparisons of protein kinase sequences from diverse organisms, in particular prokaryotes, have contributed to our understanding of the structural organization and evolution of allosteric regulation in the protein kinase domain. An emerging view from these studies is that regulatory diversity and complexity in the protein kinase domain evolved in a ‘modular’ fashion through elaboration of an ancient core component, which existed before the emergence of eukaryotes. The core component provided the conformational flexibility required for ATP binding and phosphoryl transfer in prokaryotic kinases, but evolved into a highly regulatable domain in eukaryotes through the addition of exaggerated structural features that facilitated tight allosteric control. Family and group-specific features are built upon the core component in eukaryotes to provide additional layers of control. We propose that ‘modularity’ and ‘conformational flexibility’ are key evolvable traits of the protein kinase domain that contributed to its extensive regulatory diversity and complexity.