Long‐range molecular dynamics show that inactive forms of Protein Kinase A are more dynamic than active forms

Many protein kinases are characterized by at least two structural forms corresponding to the highest level of activity (active) and low or no activity, (inactive). Further, protein dynamics is an important consideration in understanding the molecular and mechanistic basis of enzyme function. In this work, we use protein kinase A (PKA) as the model system and perform microsecond range molecular dynamics (MD) simulations on six variants which differ from one another in terms of active and inactive form, with or without bound ligands, C‐terminal tail and phosphorylation at the activation loop. We find that the root mean square fluctuations in the MD simulations are generally higher for the inactive forms than the active forms. This difference is statistically significant. The higher dynamics of inactive states has significant contributions from ATP binding loop, catalytic loop, and αG helix. Simulations with and without C‐terminal tail show this differential dynamics as well, with lower dynamics both in the active and inactive forms if C‐terminal tail is present. Similarly, the dynamics associated with the inactive form is higher irrespective of the phosphorylation status of Thr 197. A relatively stable stature of active kinases may be better suited for binding of substrates and detachment of the product. Also, phosphoryl group transfer from ATP to the phosphosite on the substrate requires precise transient coordination of chemical entities from three different molecules, which may be facilitated by the higher stability of the active state.     

Emerging Methods and Applications to Decrypt Allostery in Proteins and Nucleic Acids

Many large protein-nucleic acid complexes exhibit allosteric regulation. In these systems, the propagation of the allosteric signaling is strongly coupled to conformational dynamics and catalytic function, challenging state-of-the-art analytical methods. Here, we review established and innovative approaches used to elucidate allosteric mechanisms in these complexes. Specifically, we report network models derived from graph theory and centrality analyses in combination with molecular dynamics (MD) simulations, introducing novel schemes that implement the synergistic use of graph theory with enhanced simulations methods and ab-initio MD. Accelerated MD simulations are used to construct “enhanced network models”, describing the allosteric response over long timescales and capturing the relation between allostery and conformational changes. “Ab-initio network models” combine graph theory with ab-initio MD and quantum mechanics/molecular mechanics (QM/MM) simulations to describe the allosteric regulation of catalysis by following the step-by-step dynamics of biochemical reactions. This approach characterizes how the allosteric regulation changes from reactants to products and how it affects the transition state, revealing a tense-to-relaxed allosteric regulation along the chemical step. Allosteric models and applications are showcased for three paradigmatic examples of allostery in protein-nucleic acid complexes: (i) the nucleosome core particle, (ii) the CRISPR-Cas9 genome editing system and (iii) the spliceosome. These methods and applications create innovative protocols to determine allosteric mechanisms in protein-nucleic acid complexes that show tremendous promise for medicine and bioengineering.     

A genomic perspective of protein kinases in Plasmodium falciparum

Protein kinases are central to regulation of cellular signaling in the eukaryotes. Well-conserved and lineage-specific protein kinases have previously been identified from various completely sequenced genomes of eukaryotes. The current work describes a genome-wide analysis for protein kinases encoded in the Plasmodium falciparum genome. Using a few different profile matching methods, we have identified 99 protein kinases or related proteins in the parasite genome. We have classified these kinases into subfamilies and analyzed them in the context of noncatalytic domains that occur in these catalytic kinase domain-containing proteins. Compared to most eukaryotic protein kinases, these sequences vary significantly in terms of their lengths, inserts in catalytic domains, and co-occurring domains. Catalytic and noncatalytic domains contain long stretches of repeats of positively charged and other polar amino acids. Various components of the cell cycle, including 4 cyclin-dependent kinase (CDK) homologues, 2 cyclins, 1 CDK regulatory subunit, and 1 kinase-associated phosphatase, are identified. Identification of putative mitogen-activated protein (MAP) Kinase and MAP Kinase Kinase of P. falciparum suggests a new paradigm in the highly conserved signaling pathway of eukaryotes. The calcium-dependent kinase family, well represented in P. falciparum, shows varying domain combinations with EF-hands and pleckstrin homology domains. The analysis reveals a new subfamily of protein kinases having limited sequence similarity with previously known subfamilies. A new transmembrane kinase with 6 membrane-spanning regions is identified. Putative apicoplast targeting sequences have been detected in some of these protein kinases, suggesting their export to the apicoplast.     

TYK2 Protein-Coding Variants Protect against Rheumatoid Arthritis and Autoimmunity,with No Evidence of Major Pleiotropic Effects on Non-Autoimmune Complex Traits

Despite the success of genome-wide association studies (GWAS) in detecting a large number of loci for complex phenotypes such as rheumatoid arthritis (RA) susceptibility, the lack of information on the causal genes leaves important challenges to interpret GWAS results in the context of the disease biology. Here, we genetically fine-map the RA risk locus at 19p13 to define causal variants, and explore the pleiotropic effects of these same variants in other complex traits. First, we combined Immunochip dense genotyping (n = 23,092 case/control samples), Exomechip genotyping (n = 18,409 case/control samples) and targeted exon-sequencing (n = 2,236 case/controls samples) to demonstrate that three protein-coding variants in TYK2 (tyrosine kinase 2) independently protect against RA: P1104A (rs34536443, OR = 0.66, P = 2.3x10^-21), A928V (rs35018800, OR = 0.53, P = 1.2x10^-9), and I684S (rs12720356, OR = 0.86, P = 4.6x10^-7). Second, we show that the same three TYK2 variants protect against systemic lupus erythematosus (SLE, P_omnibus = 6x10^-18), and provide suggestive evidence that two of the TYK2 variants (P1104A and A928V) may also protect against inflammatory bowel disease (IBD; P_omnibus = 0.005). Finally, in a phenome-wide association study (PheWAS) assessing >500 phenotypes using electronic medical records (EMR) in >29,000 subjects, we found no convincing evidence for association of P1104A and A928V with complex phenotypes other than autoimmune diseases such as RA, SLE and IBD. Together, our results demonstrate the role of TYK2 in the pathogenesis of RA, SLE and IBD, and provide supporting evidence for TYK2 as a promising drug target for the treatment of autoimmune diseases.     

Allosteric activator domain of maintenance human DNA (cytosine-5) methyltransferase and its role in methylation spreading

The human maintenance DNA (cytosine-5) methyltransferase (hDNMT1) consists of a large N-terminal regulatory domain fused to a catalytic C-terminal domain by randomly repeated Gly-Lys dipeptides. Several N-terminal deletion mutants of hDNMT1 were made, purified, and tested for substrate specificity. Deletion mutants lacking 121, 501, 540, or 580 amino acids from the N-terminus still functioned as DNA methyltransferases, methylated CG sequences, and preferred hemimethylated to unmethylated DNA, as did the full-length hDNMT1. Methylated DNA stimulated methylation spreading on unmethylated CpG sequences for the full-length and the 121 amino acid deletion hDNMT1 equally well but not for the mutants lacking 501, 540, or 580 amino acids, indicating the presence of an allosteric activation determinant between amino acids 121 and 501. Peptides from the N- and C-termini bound methylated DNA independently. Point mutation analysis within the allosteric region revealed that amino acids 284-287 (KKHR) were involved in methylated DNA-mediated allosteric activation. Allosteric activation was reduced in the double point mutant enzymes D25 (K284A and K285A) and D12 (H286A and R287A). Retinoblastoma gene product (Rb), a negative regulator of DNA methylation, bound to the allosteric site of hDNMT1 and inhibited methylation, suggesting Rb may regulate methylation spreading.     

An electrostatic network and long-range regulation of Src kinases

The regulatory mechanism of Src tyrosine kinases includes conformational activation by a change in the catalytic domain tertiary structure and in domain-domain contacts between the catalytic domain and the SH2/SH3 regulatory domains. The kinase is activated when tyrosine phosphorylation occurs on the activation loop, but without phosphorylation of the C-terminal tail. Activation also occurs by allostery when contacts between the catalytic domain (CD) and the regulatory SH3 and SH2 domains are released as a result of exogenous protein binding. The aim of this work is to examine the proposed role of an electrostatic network in the conformational transition and to elucidate the molecular mechanism for long-range, allosteric conformational activation by using a combination of experimental enzyme kinetics and nonequilibrium molecular dynamics simulations. Salt dependence of the induction phase is observed in kinetic assays and supports the role of an electrostatic network in the transition. In addition, simulations provide evidence that allosteric activation involves a concerted motion coupling highly conserved residues, and spanning several nanometers from the catalytic site to the regulatory domain interface to communicate between the CD and the regulatory domains.     

Computational modeling of allosteric communication reveals organizing principles of mutation-induced signaling in ABL and EGFR kinases

The emerging structural information about allosteric kinase complexes and the growing number of allosteric inhibitors call for a systematic strategy to delineate and classify mechanisms of allosteric regulation and long-range communication that control kinase activity. In this work, we have investigated mechanistic aspects of long-range communications in ABL and EGFR kinases based on the results of multiscale simulations of regulatory complexes and computational modeling of signal propagation in proteins. These approaches have been systematically employed to elucidate organizing molecular principles of allosteric signaling in the ABL and EGFR multi-domain regulatory complexes and analyze allosteric signatures of the gate-keeper cancer mutations. We have presented evidence that mechanisms of allosteric activation may have universally evolved in the ABL and EGFR regulatory complexes as a product of a functional cross-talk between the organizing αF-helix and conformationally adaptive αI-helix and αC-helix. These structural elements form a dynamic network of efficiently communicated clusters that may control the long-range interdomain coupling and allosteric activation. The results of this study have unveiled a unifying effect of the gate-keeper cancer mutations as catalysts of kinase activation, leading to the enhanced long-range communication among allosterically coupled segments and stabilization of the active kinase form. The results of this study can reconcile recent experimental studies of allosteric inhibition and long-range cooperativity between binding sites in protein kinases. The presented study offers a novel molecular insight into mechanistic aspects of allosteric kinase signaling and provides a quantitative picture of activation mechanisms in protein kinases at the atomic level.     

Structural modulation of factor VIIa by full-length tissue factor (TF1-263): implication of novel interactions between EGF2 domain and TF

Tissue factor (TF)-mediated factor VII (FVII) activation and a subsequent proteolytic TF-FVIIa binary complex formation is the key step initiating the coagulation cascade, with implications in various homeostatic and pathologic scenarios. TF binding allosterically modifies zymogen-like free FVIIa to its highly catalytically active form. As a result of unresolved crystal structure of the full-length TF_1-263-FVIIa binary complex and free FVIIa, allosteric alterations in FVIIa following its binding to full-length TF and the consequences of these on function are not entirely clear. The present study aims to map and identify structural alterations in FVIIa and TF resulting from full-length TF binding to FVIIa and the key events responsible for enhanced FVIIa activity in coagulation. We constructed the full-length TF_1-263-FVIIa membrane bound complex using computational modeling and subjected it to molecular dynamics (MD) simulations. MD simulations showed that TF alters the structure of each domain of FVIIa and these combined alterations contribute to enhanced TF-FVIIa activity. Detailed, domain-wise investigation revealed several new non-covalent interactions between TF and FVIIa that were not found in the truncated soluble TF-FVIIa crystal structure. The structural modulation of each FVIIa domain imparted by TF indicated that both inter and intra-domain communication is crucial for allosteric modulation of FVIIa. Our results suggest that these newly formed interactions can provide additional stability to the protease domain and regulate its activity profile by governing catalytic triad (CT) orientation and localization. The unexplored newly formed interactions between EGF2 and TF provides a possible explanation for TF-induced allosteric activation of FVIIa.     

The slow-onset nature of allosteric inhibition in α-isopropylmalate synthase from Mycobacterium tuberculosis is mediated by a flexible loop

The identification of structure-function relationships in allosteric enzymes is essential to describing a molecular mechanism for allosteric processes. The enzyme α-isopropylmalate synthase from Mycobacterium tuberculosis (MtIPMS) is subject to slow-onset, allosteric inhibition by l-leucine. Here we report that alternate amino acids act as rapid equilibrium noncompetitive inhibitors of MtIPMS failing to display biphasic inhibition kinetics. Amino acid substitutions on a flexible loop covering the regulatory binding pocket generate enzyme variants that have significant affinity for l-leucine but lack biphasic inhibition kinetics. Taken together, these results are consistent with the flexible loop mediating the slow-onset step of allosteric inhibition.     

pheA (Rv3838c) of Mycobacterium tuberculosis encodes an allosterically regulated monofunctional prephenate dehydratase that requires both catalytic and regulatory domains for optimum activity

Prephenate dehydratase (PDT) is a key regulatory enzyme in l-phenylalanine biosynthesis. In Mycobacterium tuberculosis, expression of pheA, the gene encoding PDT, has been earlier reported to be iron-dependent (1, 2). We report that M. tuberculosis pheA is also regulated at the protein level by aromatic amino acids. All of the three aromatic amino acids (phenylalanine, tyrosine, and tryptophan) are potent allosteric activators of M. tuberculosis PDT. We also provide in vitro evidence that M. tuberculosis PDT does not possess any chorismate mutase activity, which suggests that, unlike many other enteric bacteria (where PDT exists as a fusion protein with chorismate mutase), M. tuberculosis PDT is a monofunctional and a non-fusion protein. Finally, the biochemical and biophysical properties of the catalytic and regulatory domains (ACT domain) of M. tuberculosis PDT were studied to observe that, in the absence of the ACT domain, the enzyme not only loses its regulatory activity but also its catalytic activity. These novel results provide evidence for a monofunctional prephenate dehydratase enzyme from a pathogenic bacterium that exhibits extensive allosteric activation by aromatic amino acids and is absolutely dependent upon the presence of catalytic as well as the regulatory domains for optimum enzyme activity.     

Characterization of an autonomously activated plant ADP-glucose pyrophosphorylase

ADP-glucose pyrophosphorylase (AGPase) catalyzes the rate-limiting step in starch biosynthesis in plants and changes in its catalytic and/or allosteric properties can lead to increased starch production. Recently, a maize (Zea mays)/potato (Solanum tuberosum) small subunit mosaic, MP [Mos(1–198)], containing the first 198 amino acids of the small subunit of the maize endosperm enzyme and the last 277 amino acids from the potato tuber enzyme, was expressed with the maize endosperm large subunit and was reported to have favorable kinetic and allosteric properties. Here, we show that this mosaic, in the absence of activator, performs like a wild-type AGPase that is partially activated with 3-phosphoglyceric acid (3-PGA). In the presence of 3-PGA, enzyme properties of Mos(1–198)/SH2 are quite similar to those of the wild-type maize enzyme. In the absence of 3-PGA, however, the mosaic enzyme exhibits greater activity, higher affinity for the substrates, and partial inactivation by inorganic phosphate. The Mos(1–198)/SH2 enzyme is also more stable to heat inactivation. The different properties of this protein were mapped using various mosaics containing smaller portions of the potato small subunit. Enhanced heat stability of Mos(1–198) was shown to originate from five potato-derived amino acids between 322 and 377. These amino acids were shown previously to be important in small subunit/large subunit interactions. These five potato-derived amino acids plus other potato-derived amino acids distributed throughout the carboxyl-terminal portion of the protein are required for the enhanced catalytic and allosteric properties exhibited by Mos(1–198)/SH2.     

Insights into the structural dynamics of Liver kinase B1 (LKB1) by the binding of STe20 Related Adapterα (STRADα) and Mouse protein 25α (MO25α) co-activators

LKB1, the tumour suppressor, is found mutated in Peutz-Jeghers syndrome (PJS). The LKB1 is a serine-threonine kinase protein that is allosterically activated by the binding of STRADα and MO25α without phosphorylating the Thr212 present at activation loop. The present study aims to highlight the structural dynamics and complexation mechanism during the allosteric activation of LKB1 by these co-activators using molecular dynamics simulations. The all atom simulations performed on the complexes of LKB1 with ATP, STRADα, and MO25α for a period of 30 ns reveal that binding of STRADα and MO25α significantly stabilizes the highly flexible regions of LKB1 such as ATP binding region (β1-β2 loop), catalytic & activation loop segments and αG helix. Also, binding of STRADα and MO25α to LKB1 promotes coordinated motion between N- and C-lobes along with the catalytic & activation loops by forming H-bonds between LKB1 and co-activators, which further facilitate to establish the conserved attributes of active LKB1 such as (i) formation of salt bridge between Lys78 and Glu98, (ii) formation of stable hydrophobic R- and C-spines, and (iii) interaction between both catalytic and activation loops. Especially, the residues of LKB1 interacting with STRADα (Arg74, Glu342) and MO25α (Glu165, Pro203 and Phe204) are observed to play a significant role in stabilizing the (LKB1-ATP)-(STRADα-ATP)-MO25α complex. Overall, the present work highlighting the structural dynamics of LKB1 by the binding of allosteric co-activators is expected to provide a basic understanding on drug design specific to PJS syndrome.     

Intermolecular insights into allosteric inhibition of histone lysine-specific demethylase 1

BackgroundHistone lysine-specific demethylase 1 (LSD1) has become a potential anticancer target for the novel drug discovery. Recent reports have shown that SP2509 and its derivatives strongly inhibit LSD1 as allosteric inhibitors. However, the binding mechanism of these allosteric inhibitors in the allosteric site of LSD1 is not known yet.MethodsThe stability and binding mechanism of allosteric inhibitors in the binding site of LSD1 were evaluated by molecular docking, ligand-based pharmacophore, molecular dynamics (MD) simulations, molecular mechanics generalized born surface area (MM/GBSA) analysis, quantum mechanics/molecular mechanics (QM/MM) calculation and Hirshfeld surface analysis.ResultsThe conformational geometry and the intermolecular interactions of allosteric inhibitors showed high binding affinity towards allosteric site of LSD1 with the neighboring amino acids (Gly358, Cys360, Leu362, Asp375 and Glu379). Meanwhile, MD simulations and MM/GBSA analysis were performed on selected allosteric inhibitors in complex with LSD1 protein, which confirmed the high stability and binding affinity of these inhibitors in the allosteric site of LSD1.ConclusionThe simulation results revealed the crucial factors accounting for allosteric inhibitors of LSD1, including different protein–ligand interactions, the positions and conformations of key residues, and the ligands flexibilities. Meanwhile, a halogen bond interaction between chlorine atom of ligand and key residues Trp531 and His532 was recurrent in our analysis confirming its importance.General significanceOverall, our research analyzed in depth the binding modes of allosteric inhibitors with LSD1 and could provide useful information for the design of novel allosteric inhibitors.     

Structural and Functional Characterization of the JH2 Pseudokinase Domain of JAK Family Tyrosine Kinase 2 (TYK2)

JAK (Janus family of cytoplasmic tyrosine kinases) family tyrosine kinase 2 (TYK2) participates in signaling through cytokine receptors involved in immune responses and inflammation. JAKs are characterized by dual kinase domain: a tyrosine kinase domain (JH1) that is preceded by a pseudokinase domain (JH2). The majority of disease-associated mutations in JAKs map to JH2, demonstrating its central regulatory function. JH2s were considered catalytically inactive, but JAK2 JH2 was found to have low autoregulatory catalytic activity. Whether the other JAK JH2s share ATP binding and enzymatic activity has been unclear. Here we report the crystal structure of TYK2 JH2 in complex with adenosine 5′-O-(thiotriphosphate) (ATP-γS) and characterize its nucleotide binding by biochemical and biophysical methods. TYK2 JH2 did not show phosphotransfer activity, but it binds ATP and the nucleotide binding stabilizes the protein without inducing major conformational changes. Mutation of the JH2 ATP-binding pocket increased basal TYK2 phosphorylation and downstream signaling. The overall structural characteristics of TYK2 JH2 resemble JAK2 JH2, but distinct stabilizing molecular interactions around helix αAL in the activation loop provide a structural basis for differences in substrate access and catalytic activities among JAK family JH2s. The structural and biochemical data suggest that ATP binding is functionally important for both TYK2 and JAK2 JH2s, whereas the regulatory phosphorylation appears to be a unique property of JAK2. Finally, the co-crystal structure of TYK2 JH2 complexed with a small molecule inhibitor demonstrates that JH2 is accessible to ATP-competitive compounds, which offers novel approaches for targeting cytokine signaling as well as potential therapeutic applications.     

Structural, kinetic and computational investigation of Vitis vinifera DHDPS reveals new insight into the mechanism of lysine-mediated allosteric inhibition

Lysine is one of the most limiting amino acids in plants and its biosynthesis is carefully regulated through inhibition of the first committed step in the pathway catalyzed by dihydrodipicolinate synthase (DHDPS). This is mediated via a feedback mechanism involving the binding of lysine to the allosteric cleft of DHDPS. However, the precise allosteric mechanism is yet to be defined. We present a thorough enzyme kinetic and thermodynamic analysis of lysine inhibition of DHDPS from the common grapevine, Vitis vinifera (Vv). Our studies demonstrate that lysine binding is both tight (relative to bacterial DHDPS orthologs) and cooperative. The crystal structure of the enzyme bound to lysine (2.4 Å) identifies the allosteric binding site and clearly shows a conformational change of several residues within the allosteric and active sites. Molecular dynamics simulations comparing the lysine-bound (PDB ID 4HNN) and lysine free (PDB ID 3TUU) structures show that Tyr132, a key catalytic site residue, undergoes significant rotational motion upon lysine binding. This suggests proton relay through the catalytic triad is attenuated in the presence of lysine. Our study reveals for the first time the structural mechanism for allosteric inhibition of DHDPS from the common grapevine.     

Conserved internal hydration motifs in protein kinases

Crystal structures of diverse protein kinase catalytic subunits reveal a number of water molecules whose positions within the protein core are better preserved than amino acid types in many functionally important locations. It remains unknown whether they play any particular role, and whether their removal, disturbing local interaction patterns to no smaller degree than amino acid mutations, can affect kinase stability and function. In this study, we apply an array of computational approaches to characterize hydration of kinase catalytic subunits. First, we systematically screen multiple crystal structures with the use of a simplified hydration model in order to determine the distribution of internal solvent and the degree of its conservation. Second, we analyze water structure, dynamics and binding affinity to buried hydration sites in a number of kinases, also taking into account their variable functional state. We find that a large portion of buried solvent is dynamic, possibly contributing to kinase conformational changes related to the activation process. In turn, binding free energies of some of tightly bound conserved water molecules to different kinases tend to shift in a similar manner following the change of their functional state. This finding highlights the likely specific role of internal solvent in fine tuning local protein plasticity.     

An absolute requirement of fructose 1,6-bisphosphate for the Lactobacillus casei L-lactate dehydrogenase activity induced by a single amino acid substitution

Lactobacillus casei allosteric L-lactate dehydrogenase (L-LDH) absolutely requires fructose 1,6-bisphosphate [Fru(1,6)P2] for its catalytic activity under neutral conditions, but exhibits marked catalytic activity in the absence of Fru(1,6)P(2) under acidic conditions through the homotropic activation effect of substrate pyruvate. In this enzyme, a single amino acid replacement, i.e. that of His205 conserved in the Fru(1,6)P(2)-binding site of certain allosteric L-LDHs of lactic acid bacteria with Thr, did not induce a marked loss of the activation effect of Fru(1,6)P(2) or divalent metal ions, which are potent activators that improve the activation function of Fru(1,6)P(2) under neutral conditions. However, this replacement induced a great loss of the Fru(1,6)P(2)-independent activation effect of pyruvate or pyruvate analogs under acidic conditions, consequently indicating an absolute Fru(1,6)P(2) requirement for the enzyme activity. The replacement also induced a significant reduction in the pH-dependent sensitivity of the enzyme to Fru(1,6)P(2), through a slight decrease and increase of the Fru(1,6)P(2) sensitivity under acidic and neutral conditions, respectively, indicating that His205 is also largely involved in the pH-dependent sensitivity of L.casei L-LDH to Fru(1,6)P(2). The role of His205 in the allosteric regulation of the enzyme is discussed on the basis of the known crystal structures of L-LDHs.     

Evolutionary constraints associated with functional specificity of the CMGC protein kinases MAPK, CDK, GSK, SRPK, DYRK, and CK2alpha

Amino acid residues associated with functional specificity of cyclin-dependent kinases (CDKs), mitogen-activated protein kinases (MAPKs), glycogen synthase kinases (GSKs), and CDK-like kinases (CLKs), which are collectively termed the CMGC group, were identified by categorizing and quantifying the selective constraints acting upon these proteins during evolution. Many constraints specific to CMGC kinases correspond to residues between the N-terminal end of the activation segment and a CMGC-conserved insert segment associated with coprotein binding. The strongest such constraint is imposed on a "CMGC-arginine" near the substrate phosphorylation site with a side chain that plays a role both in substrate recognition and in kinase activation. Two nearby buried waters, which are also present in non-CMGC kinases, typically position the main chain of this arginine relative to the catalytic loop. These and other CMGC-specific features suggest a structural linkage between coprotein binding, substrate recognition, and kinase activation. Constraints specific to individual subfamilies point to mechanisms for CMGC kinase specialization. Within casein kinase 2alpha (CK2alpha), for example, the binding of one of the buried waters appears prohibited by the side chain of a leucine that is highly conserved within CK2alpha and that, along with substitution of lysine for the CMGC-arginine, may contribute to the broad substrate specificity of CK2alpha by relaxing characteristically conserved, precise interactions near the active site. This leucine is replaced by a conserved isoleucine or valine in other CMGC kinases, thereby illustrating the potential functional significance of subtle amino acid substitutions. Analysis of other CMGC kinases similarly suggests candidate family-specific residues for experimental follow-up.     

Divergent allosteric patterns verify the regulatory paradigm for aspartate transcarbamylase

The native Escherichia coli aspartate transcarbamoylase (ATCase, E.C. 2.1.3.2) provides a classic allosteric model for the feedback inhibition of a biosynthetic pathway by its end products. Both E. coli and Erwinia herbicola possess ATCase holoenzymes which are dodecameric (2(c3):3(r2)) with 311 amino acid residues per catalytic monomer and 153 and 154 amino acid residues per regulatory (r) monomer, respectively. While the quaternary structures of the two enzymes are identical, the primary amino acid sequences have diverged by 14 % in the catalytic polypeptide and 20 % in the regulatory polypeptide. The amino acids proposed to be directly involved in the active site and nucleotide binding site are strictly conserved between the two enzymes; nonetheless, the two enzymes differ in their catalytic and regulatory characteristics. The E. coli enzyme has sigmoidal substrate binding with activation by ATP, and inhibition by CTP, while the E. herbicola enzyme has apparent first order kinetics at low substrate concentrations in the absence of allosteric ligands, no ATP activation and only slight CTP inhibition. In an apparently important and highly conserved characteristic, CTP and UTP impose strong synergistic inhibition on both enzymes. The co-operative binding of aspartate in the E. coli enzyme is correlated with a T-to-R conformational transition which appears to be greatly reduced in the E. herbicola enzyme, although the addition of inhibitory heterotropic ligands (CTP or CTP+UTP) re-establishes co-operative saturation kinetics. Hybrid holoenzymes assembled in vivo with catalytic subunits from E. herbicola and regulatory subunits from E. coli mimick the allosteric response of the native E. coli holoenzyme and exhibit ATP activation. The reverse hybrid, regulatory subunits from E. herbicola and catalytic subunits from E. coli, exhibited no response to ATP. The conserved structure and diverged functional characteristics of the E. herbicola enzyme provides an opportunity for a new evaluation of the common paradigm involving allosteric control of ATCase.