Gene Promoter Evolution Targets the Center of the Human Protein Interaction Network

Assessing the contribution of promoters and coding sequences to gene evolution is an important step toward discovering the major genetic determinants of human evolution. Many specific examples have revealed the evolutionary importance of cis-regulatory regions. However, the relative contribution of regulatory and coding regions to the evolutionary process and whether systemic factors differentially influence their evolution remains unclear. To address these questions, we carried out an analysis at the genome scale to identify signatures of positive selection in human proximal promoters. Next, we examined whether genes with positively selected promoters (Prom+ genes) show systemic differences with respect to a set of genes with positively selected protein-coding regions (Cod+ genes). We found that the number of genes in each set was not significantly different (8.1% and 8.5%, respectively). Furthermore, a functional analysis showed that, in both cases, positive selection affects almost all biological processes and only a few genes of each group are located in enriched categories, indicating that promoters and coding regions are not evolutionarily specialized with respect to gene function. On the other hand, we show that the topology of the human protein network has a different influence on the molecular evolution of proximal promoters and coding regions. Notably, Prom+ genes have an unexpectedly high centrality when compared with a reference distribution (P = 0.008, for Eigenvalue centrality). Moreover, the frequency of Prom+ genes increases from the periphery to the center of the protein network (P = 0.02, for the logistic regression coefficient). This means that gene centrality does not constrain the evolution of proximal promoters, unlike the case with coding regions, and further indicates that the evolution of proximal promoters is more efficient in the center of the protein network than in the periphery. These results show that proximal promoters have had a systemic contribution to human evolution by increasing the participation of central genes in the evolutionary process.

Vacuolar SNARE Protein Transmembrane Domains Serve as Nonspecific Membrane Anchors with Unequal Roles in Lipid Mixing.

Membrane fusion is induced by SNARE complexes that are anchored in both fusion partners. SNAREs zipper up from the N to C terminus bringing the two membranes into close apposition. Their transmembrane domains (TMDs) might be mere anchoring devices, deforming bilayers by mechanical force. Structural studies suggested that TMDs might also perturb lipid structure by undergoing conformational transitions or by zipping up into the bilayer. Here, we tested this latter hypothesis, which predicts that the activity of SNAREs should depend on the primary sequence of their TMDs. We replaced the TMDs of all vacuolar SNAREs (Nyv1, Vam3, and Vti1) by a lipid anchor, by a TMD from a protein unrelated to the membrane fusion machinery, or by artificial leucine-valine sequences. Individual exchange of the native SNARE TMDs against an unrelated transmembrane anchor or an artificial leucine-valine sequence yielded normal fusion activities. Fusion activity was also preserved upon pairwise exchange of the TMDs against unrelated peptides, which eliminates the possibility for specific TMD-TMD interactions. Thus, a specific primary sequence or zippering beyond the SNARE domains is not a prerequisite for fusion. Lipid-anchored Vti1 was fully active, and lipid-anchored Nyv1 permitted the reaction to proceed up to hemifusion, and lipid-anchored Vam3 interfered already before hemifusion. The unequal contribution of proteinaceous TMDs on Vam3 and Nyv1 suggests that Q- and R-SNAREs might make different contributions to the hemifusion intermediate and the opening of the fusion pore. Furthermore, our data support the view that SNARE TMDs serve as nonspecific membrane anchors in vacuole fusion.

Protein interaction studies point to new functions for Escherichia coli glyceraldehyde-3-phosphate dehydrogenase

Glyceraldehyde-3-phosphate dehydrogenase (GAPDH) is considered a multifunctional protein with defined functions in numerous mammalian cellular processes. GAPDH functional diversity depends on various factors such as covalent modifications, subcellular localization, oligomeric state and intracellular concentration of substrates or ligands, as well as protein-protein interactions. In bacteria, alternative GAPDH functions have been associated with its extracellular location in pathogens or probiotics. In this study, new intracellular functions of E. coli GAPDH were investigated following a proteomic approach aimed at identifying interacting partners using in vivo formaldehyde cross-linking followed by mass spectrometry. The identified proteins were involved in metabolic processes, protein synthesis and folding or DNA repair. Some interacting proteins were also identified in immunopurification experiments in the absence of cross-linking. Pull-down experiments and overlay immunoblotting were performed to further characterize the interaction with phosphoglycolate phosphatase (Gph). This enzyme is involved in the metabolism of 2-phosphoglycolate formed in the DNA repair of 3"-phosphoglycolate ends generated by bleomycin damage. We show that interaction between Gph and GAPDH increases in cells challenged with bleomycin, suggesting involvement of GAPDH in cellular processes linked to DNA repair mechanisms.

Protein loop compaction and the origin of the effect of arginine and glutamic acid mixtures on solubility, stability and transient oligomerization of proteins

Addition of a 50 mM mixture of l-arginine and l-glutamic acid (RE) is extensively used to improve protein solubility and stability, although the origin of the effect is not well understood. We present Small Angle X-ray Scattering (SAXS) and Nuclear Magnetic Resonance (NMR) results showing that RE induces protein compaction by collapsing flexible loops on the protein core. This is suggested to be a general mechanism preventing aggregation and improving resistance to proteases and to originate from the polyelectrolyte nature of RE. Molecular polyelectrolyte mixtures are expected to display long range correlation effects according to dressed interaction site theory. We hypothesize that perturbation of the RE solution by dissolved proteins is proportional to the volume occupied by the protein. As a consequence, loop collapse, minimizing the effective protein volume, is favored in the presence of RE.

Mechanisms of G protein activation via the D-2 dopamine receptor: evidence for persistent receptor/G protein interaction after agonist stimulation

Background and purpose: The aim of this report is to study mechanisms of G protein activation by agonists. Experimental approach: The association and dissociation of guanosine 5'-O-(3-[S-35] thio) triphosphate ([S-35] GTP gamma S) binding at G proteins in membranes of CHO cells stably transfected with the human dopamine D-2short receptor was studied in the presence of a range of agonists. Key results: Binding of [S-35] GTPgS was dissociable in the absence of agonist and dissociation was accelerated both in rate and extent by dopamine, an effect which was blocked by the dopamine D-2 receptor antagonist raclopride and by suramin, which inhibits receptor/G protein interaction. A range of agonists of varying efficacy increased the rate of dissociation of [S-35] GTPgS binding, with the more efficacious agonists resulting in faster dissociation. Agonists were able to dissociate about 70% of the pre-bound [S-35] GTPgS, leaving a component which may not be accessible to the agonist-bound receptor. The dissociable component of the [S-35] GTPgS binding was reduced with longer association times and increased [S-35] GTPgS concentrations. Conclusions and implications: These data are consistent with [S-35] GTPgS binding being initially to receptor-linked G proteins and then to G proteins which have separated from the agonist bound receptor. Under the conditions used typically for [S-35] GTPgS binding assays, therefore, much of the agonist-receptor complex remains in proximity to G proteins after they have been activated by agonist.

Role of alpha 7 Nicotinic Acetylcholine Receptor in Calcium Signaling Induced by Prion Protein Interaction with Stress-inducible Protein 1

The prion protein (PrP(C)) is a conserved glycosylphosphatidyl-inositol-anchored cell surface protein expressed by neurons and other cells. Stress-inducible protein 1 (STI1) binds PrP(C) extracellularly, and this activated signaling complex promotes neuronal differentiation and neuroprotection via the extracellular signal-regulated kinase 1 and 2 (ERK1/2) and cAMP-dependent protein kinase 1 (PKA) pathways. However, the mechanism by which the PrPC-STI1 interaction transduces extracellular signals to the intracellular environment is unknown. We found that in hippocampal neurons, STI1-PrP(C) engagement induces an increase in intracellular Ca(2+) levels. This effect was not detected in PrP(C)-null neurons or wild-type neurons treated with an STI1 mutant unable to bind PrP(C). Using a best candidate approach to test for potential channels involved in Ca(2+) influx evoked by STI1-PrP(C), we found that alpha-bungarotoxin, a specific inhibitor for alpha 7 nicotinic acetylcholine receptor (alpha 7nAChR), was able to block PrP(C)-STI1-mediated signaling, neuroprotection, and neuritogenesis. Importantly, when alpha 7nAChR was transfected into HEK 293 cells, it formed a functional complex with PrP(C) and allowed reconstitution of signaling by PrP(C)-STI1 interaction. These results indicate that STI1 can interact with the PrP(C).alpha 7nAChR complex to promote signaling and provide a novel potential target for modulation of the effects of prion protein in neurodegenerative diseases.

Resilience of protein-protein interaction networks as determined by their large-scale topological features

The relationship between the structure and function of biological networks constitutes a fundamental issue in systems biology. Particularly, the structure of protein-protein interaction networks is related to important biological functions. In this work, we investigated how such a resilience is determined by the large scale features of the respective networks. Four species are taken into account, namely yeast Saccharomyces cerevisiae, worm Caenorhabditis elegans, fly Drosophila melanogaster and Homo sapiens. We adopted two entropy-related measurements (degree entropy and dynamic entropy) in order to quantify the overall degree of robustness of these networks. We verified that while they exhibit similar structural variations under random node removal, they differ significantly when subjected to intentional attacks (hub removal). As a matter of fact, more complex species tended to exhibit more robust networks. More specifically, we quantified how six important measurements of the networks topology (namely clustering coefficient, average degree of neighbors, average shortest path length, diameter, assortativity coefficient, and slope of the power law degree distribution) correlated with the two entropy measurements. Our results revealed that the fraction of hubs and the average neighbor degree contribute significantly for the resilience of networks. In addition, the topological analysis of the removed hubs indicated that the presence of alternative paths between the proteins connected to hubs tend to reinforce resilience. The performed analysis helps to understand how resilience is underlain in networks and can be applied to the development of protein network models.

Sdo1p, the yeast orthologue of Shwachman-Bodian-Diamond syndrome protein, binds RNA and interacts with nuclear rRNA-processing factors

The Shwachman-Bodian-Diamond syndrome protein (SBDS) is a member of a highly conserved protein family of not well understood function, with putative orthologues found in different organisms ranging from Archaea, yeast and plants to vertebrate animals. The yeast orthologue of SBDS, Sdo1p, has been previously identified in association with the 60S ribosomal subunit and is proposed to participate in ribosomal recycling. Here we show that Sdo1p interacts with nucleolar rRNA processing factors and ribosomal proteins, indicating that it might bind the pre-60S complex and remain associated with it during processing and transport to the cytoplasm. Corroborating the protein interaction data, Sdo1p localizes to the nucleus and cytoplasm and co-immunoprecipitates precursors of 60S and 40S subunits, as well as the mature rRNAs. Sdo1p binds RNA directly, suggesting that it may associate with the ribosomal subunits also through RNA interaction. Copyright (C) 2009 John Wiley & Sons, Ltd.

Protein, methionine+cystine and lysine levels for Japanese quails during the production phase

An experiment was conducted at Faculdade de Medicina Veterinária e Zootecnia/Unesp - Botucatu for 168 days to evaluate the effects of protein, Met + Cys and lysine diet levels on egg production and egg quality of laying Japanese quails. Quails with 42 days of age were reared in a completely randomized design. There were 1,944 quails distributed in four replicates of 27 birds per pen, according to a factorial 3x3x2 with three crude protein levels (16, 18 and 20% CP), three Met + Cys levels (0.700; 0.875 and 1.050%) and two lysine levels (1.100 and 1.375%). Birds fed diets with 18 and 20% CP had higher feed intake and egg production than those fed diets with 16% CP. There was significant interaction (p<0.05) between protein and Met + Cys levels on egg weight. There was no effects (p>0.05) of the protein level on feed conversion per dozen eggs; however, improved feed conversion per egg mass was seen for birds fed diets with 20% CP compared to those fed diets with 16% and 18% CP. Protein and lipid percentage in the yolk increased when dietary protein level increased from 16 to 18%. Increasing Met + Cys from 0.700% to 0.875% reduced yolk protein percentage. Higher lipid percentage in the yolk was seen in eggs from quails fed diets with 1.050% Met + Cys, whereas 1.375% lysine in the diet of resulted in decreased egg production and egg mass, besides poorer feed conversion per dozen eggs and per egg mass.

Aspects of poriferan (Suberites domuncula) apoptosis and innate immune system and evolutionary implications

Survivin, a unique member of the family of inhibitors of apoptosis (IAP) proteins, orchestrates intracellular pathways during cell division and apoptosis. Its central regulatory function in vertebrate molecular pathways as mitotic regulator and inhibitor of apoptotic cell death has major implications for tumor cell proliferation and viability, and has inspired several approaches that target survivin for cancer therapy. Analyses in early-branching Metazoa so far propose an exclusive role of survivin as a chromosomal passenger protein, whereas only later during evolution the second, complementary antiapoptotic function might have arisen, concurrent with increased organismal complexity. To lift the veil on the ancestral function(s) of this key regulatory molecule, a survivin homologue of the phylogenetically oldest extant metazoan taxon (phylum Porifera) was identified and functionally characterized. SURVL of the demosponge Suberites domuncula shares significant similarities with its metazoan homologues, ranging from conserved exon/intron structures to the presence of localization signal and protein-interaction domains, characteristic of IAP proteins. Whereas sponge tissue displayed a very low steady-state level, SURVL expression was significantly up-regulated in rapidly proliferating primmorph cells. In addition, challenge of sponge tissue and primmorphs with cadmium and the lipopeptide Pam3Cys-Ser-(Lys)4 stimulated SURVL expression, concurrent with the expression of newly discovered poriferan caspases (CASL and CASL2). Complementary functional analyses in transfected HEK-293 revealed that heterologous expression of poriferan survivin in human cells not only promotes cell proliferation but also augments resistance to cadmium-induced cell death. Taken together, these results demonstrate both a deep evolutionary conserved and fundamental dual role of survivin, and an equally conserved central position of this key regulatory molecule in interconnected pathways of cell cycle and apoptosis. Additionally, SDCASL, SDCASL2, and SDTILRc (TIR-LRR containing protein) may represent new components of the innate defense sentinel in sponges. SDCASL and SDCASL2 are two new caspase-homolog proteins with a singular structure. In addition to their CASc domains, SDCASL and SDCASL2 feature a small prodomain NH2-terminal (effector caspases) and a remarkably long COOH-terminal domain containing one or several functional double stranded RNA binding domains (dsrm). This new caspase prototype can characterize a caspase specialization coupling pathogen sensing and apoptosis, and could represent a very efficient defense mechanism. SDTILRc encompasses also a unique combination of domains: several leucine rich repeats (LRR) and a Toll/IL-1 receptor (TIR) domain. This unusual domain association may correspond to a new family of intracellular sensing protein, forming a subclass of pattern recognition receptors (PRR).

Topological aspects of the human protein interaction network

It is currently widely accepted that the understanding of complex cell functions depends on an integrated network theoretical approach and not on an isolated view of the different molecular agents. Aim of this thesis was the examination of topological properties that mirror known biological aspects by depicting the human protein network with methods from graph- and network theory. The presented network is a partial human interactome of 9222 proteins and 36324 interactions, consisting of single interactions reliably extracted from peer-reviewed scientific publications. In general, one can focus on intra- or intermodular characteristics, where a functional module is defined as "a discrete entity whose function is separable from those of other modules". It is found that the presented human network is also scale-free and hierarchically organised, as shown for yeast networks before. The interactome also exhibits proteins with high betweenness and low connectivity which are biologically analyzed and interpreted here as shuttling proteins between organelles (e.g. ER to Golgi, internal ER protein translocation, peroxisomal import, nuclear pores import/export) for the first time. As an optimisation for finding proteins that connect modules, a new method is developed here based on proteins located between highly clustered regions, rather than regarding highly connected regions. As a proof of principle, the Mediator complex is found in first place, the prime example for a connector complex. Focusing on intramodular aspects, the measurement of k-clique communities discriminates overlapping modules very well. Twenty of the largest identified modules are analysed in detail and annotated to known biological structures (e.g. proteasome, the NFκB-, TGF-β complex). Additionally, two large and highly interconnected modules for signal transducer and transcription factor proteins are revealed, separated by known shuttling proteins. These proteins yield also the highest number of redundant shortcuts (by calculating the skeleton), exhibit the highest numbers of interactions and might constitute highly interconnected but spatially separated rich-clubs either for signal transduction or for transcription factors. This design principle allows manifold regulatory events for signal transduction and enables a high diversity of transcription events in the nucleus by a limited set of proteins. Altogether, biological aspects are mirrored by pure topological features, leading to a new view and to new methods that assist the annotation of proteins to biological functions, structures and subcellular localisations. As the human protein network is one of the most complex networks at all, these results will be fruitful for other fields of network theory and will help understanding complex network functions in general.

NETWORK TOPOLOGY IN HUMAN PROTEIN INTERACTION DATA PREDICTS FUNCTIONAL ASSOCIATION

High-throughput assays, such as yeast two-hybrid system, have generated a huge amount of protein-protein interaction (PPI) data in the past decade. This tremendously increases the need for developing reliable methods to systematically and automatically suggest protein functions and relationships between them. With the available PPI data, it is now possible to study the functions and relationships in the context of a large-scale network. To data, several network-based schemes have been provided to effectively annotate protein functions on a large scale. However, due to those inherent noises in high-throughput data generation, new methods and algorithms should be developed to increase the reliability of functional annotations. Previous work in a yeast PPI network (Samanta and Liang, 2003) has shown that the local connection topology, particularly for two proteins sharing an unusually large number of neighbors, can predict functional associations between proteins, and hence suggest their functions. One advantage of the work is that their algorithm is not sensitive to noises (false positives) in high-throughput PPI data. In this study, we improved their prediction scheme by developing a new algorithm and new methods which we applied on a human PPI network to make a genome-wide functional inference. We used the new algorithm to measure and reduce the influence of hub proteins on detecting functionally associated proteins. We used the annotations of the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) as independent and unbiased benchmarks to evaluate our algorithms and methods within the human PPI network. We showed that, compared with the previous work from Samanta and Liang, our algorithm and methods developed in this study improved the overall quality of functional inferences for human proteins. By applying the algorithms to the human PPI network, we obtained 4,233 significant functional associations among 1,754 proteins. Further comparisons of their KEGG and GO annotations allowed us to assign 466 KEGG pathway annotations to 274 proteins and 123 GO annotations to 114 proteins with estimated false discovery rates of <21% for KEGG and <30% for GO. We clustered 1,729 proteins by their functional associations and made pathway analysis to identify several subclusters that are highly enriched in certain signaling pathways. Particularly, we performed a detailed analysis on a subcluster enriched in the transforming growth factor β signaling pathway (P<10-50) which is important in cell proliferation and tumorigenesis. Analysis of another four subclusters also suggested potential new players in six signaling pathways worthy of further experimental investigations. Our study gives clear insight into the common neighbor-based prediction scheme and provides a reliable method for large-scale functional annotations in this post-genomic era.

NADPH-CYTOCHROME P-450 REDUCTASE: EVIDENCE FOR TWO SEPARATE STRUCTURAL DOMAINS AND ISOLATION AND CHARACTERIZATION OF THE MEMBRANE ASSOCIATED SEGMENT

Homogenous detergent-solubilized NADPH-Cytochrome P-450 reductase was incorporated into microsomes and liposomes. This binding occurred spontaneously at temperatures between 4(DEGREES) and 37(DEGREES) and appeared to involve hydrophobic forces as the binding was not disrupted by 0.5 M sodium chloride. This exogenously-added reductase was active catalytically towards native cytochrome P-450, suggesting an association with the microsomal membrane similar to endogenous reductase. Homogeneous detergent-solubilized reductase was disaggregated by Renex-690 micelles, confirming the presence of a hydrophobic combining region on the enzyme. In contrast to these results, steapsin protease-solubilized reductase was incapable of microsomal attachment and did not interact with Renex-690 micelles. Detergent-solubilized reductase (76,500 daltons) was converted into a form with the electrophoretic mobility of steapsin protease-solubilized reductase (68,000 daltons) and a 12,500 dalton peptide (as determined by polyacrylamide-SDS gel electrophoresis) when the liposomal-incorporated enzyme was incubated with steapsin protease. The 68,000 dalton fragment thus obtained had properties identical with steapsin protease-solubilized reductase, i.e. it was catalytically active towards cytochrome c but inactive towards cytochrome P-450 and did not bind liposomes. The 12,500 dalton fragment remained associated with the liposomes when the digest was fractionated by gel filtration, suggesting that this is the segment of the enzyme which is embedded in the phospholipid bilayer. Thus, detergent-solubilized reductase appears to contain a soluble catalytic domain and a separate and separable membrane-binding domain. This latter domain is required for attaching the enzyme to the membrane and also to facilitate the catalytic interaction between the reductase and its native electron acceptor, cytochrome P-450. The membrane-binding segment of the reductase was isolated by preparative gel electrophoresis in SDS following its generation by proteolytic treatment of liposome-incorporated reductase. The peptide has a molecular weight of 6,400 as determined by gel filtration in 8 M guanidine hydrochloride and has an amino acid composition which is not especially hydrophobic. Following removal of SDS and dialysis out of 6 M urea, the membrane-binding peptide was unable to inhibit the activity of a reconstituted system containing purified reductase and cytochrome P-450. Moreover, when reductase and cytochrome P-450 were added to liposomes which contained the membrane-binding peptide, it was determined that mixed function oxidase activity was reconstituted as effectively as when vesicles without the membrane-binding peptide were used. Thus, the membrane-binding peptide was ineffective as an inhibitor of mixed function oxidase activity, suggesting perhaps that it facilitates catalysis by anchoring the catalytic domain of the reductase proximal to cytochrome P-450 (i.e. in the same mixed micelle) rather than through a specific interaction with cytochrome P-450. ^

Inhibitor of cytochrome c messenger RNA -protein interaction in stimulated rat skeletal muscle

The purpose of this work was to examine the possible mechanisms for the regulation of cytochrome c gene expression in response to increased contractile activity in rat skeletal muscle. The working hypothesis was that increased contractile activity enhances cytochrome c gene expression through a cis-element. A 110% increase in cytochrome c mRNA concentration was observed in tibialis anterior (TA) muscle after 9 days of chronic stimulation. Similar difference (120%) exists between soleus (SO) muscle of higher contractile activity and white vastus lateralis (WV) muscle of lower contractile activity. These results suggest that the endogenous cytochrome c gene expression is regulated by contractile activity. Cytochrome c-reporter genes were injected into skeletal muscles to identify the cis-element that is responsible for the regulation. Although the data was inconclusive, part of it suggested the importance of the 3$\sp\prime$-untranslated region (3$\sp\prime$-UTR) in mediating the response to increased contractile activity.^ RNA gel mobility shift (GMSA) and ultraviolet (UV) cross-linking assays revealed specific RNA-protein interaction in a 50-nucleotide region of the 3$\sp\prime$-UTR in unstimulated TA muscle. Computer analysis predicted a stem-loop structure of 17 nucleotides, which provides a structural basis for RNA-protein interaction. These 17 nucleotides are 100% conserved among rat, mouse and human cytochrome c genes and their 13 pseudogenes, suggesting a functional role for this region. The RNA-protein interaction was significantly less in highly active SO muscle than in inactive WV muscle and was dramatically decreased in stimulated TA muscle due to a protein inhibitor(s) associated with ribosome. It is possible that cytochrome c mRNAs undergoing translation are subject to a compartmentalized regulatory influence.^ The conclusion from these results is that increases in contractile activity induce or activate a protein inhibitor(s) associated with ribosome in rat skeletal muscle. The inhibitor decreases RNA-protein interaction in the 3$\sp\prime$-UTR of cytochrome c mRNA, which may result in increased mRNA stability and/or translation. ^

Inactivation and self -association of CaM kinase: Effects of ATP, substrate binding domains, and isozymic forms

Calcium/calmodulin-dependent protein kinase II (CaM kinase) is a multifunctional Ser/Thr protein kinase, that is highly enriched in brain and is involved in regulating many aspects of neuronal function. We observed that forebrain CaM kinase from crude homogenates, cytosolic fractions and purified preparations inactivates and translocates into the particulate fraction following autophosphorylation. Using purified forebrain CaM kinase as well as recombinant $\alpha$ isozyme, we determined that the formation of particulate enzyme was due to enzyme self-association. The conditions of autophosphorylation determine whether enzyme self-association and/or inactivation will occur. Self-association of CaM kinase is sensitive to pH, ATP concentration, and enzyme autophosphorylation. This process is prevented by saturating concentrations of ATP. However, in limiting ATP, pH is the dominant factor, and enzyme self-association occurs at pH values $\rm{<}7.0.$ Site-specific mutants were produced by substituting Ala for Thr286, Thr253, or Thr305,306 to determine whether these sites of autophosphorylation affect enzyme inactivation and self-association. The only mutation that influenced these processes was Ala286, which removed the protective effect afforded by autophosphorylation in saturating ATP. Enzyme inactivation occurs in the presence and absence of self-association and appears predominantly sensitive to nucleotide concentration, because saturating concentrations of $\rm Mg\sp{2+}/ADP$ or $\rm Mg\sp{2+}/ATP$ prevent this process. These data implicate the ATP binding pocket in both inactivation and self-association. We also observed that select peptide substrates and peptide inhibitors modeled after the autoregulatory domain of CaM kinase prevented these processes. The $\alpha$ and $\beta$ isozymes of CaM kinase were characterized independently, and were observed to exhibit differences in both enzyme inactivation and self-association. The $\beta$ isozyme was less sensitive to inactivation, and was never observed to self-associate. Biophysical characterization, and transmission electron microscopy coupled with image analysis indicated both isozymes were multimeric, however, the $\alpha$ and $\beta$ isozymes appeared structurally different. We hypothesize that the $\alpha$ subunit of CaM kinase plays both a structural and enzymatic role, and the $\beta$ subunit plays an enzymatic role. The ramifications for the functional differences observed for inactivation and self-association are discussed based on potential structural differences and autoregulation of the $\alpha$ and $\beta$ isozymes in both calcium-induced physiological and pathological processes. ^