Evolution of allostery in the cyclic nucleotide binding module
© Kannan et al.; licensee BioMed Central Ltd. 2007
Received: 29 August 2007
Accepted: 12 December 2007
Published: 12 December 2007
The cyclic nucleotide binding (CNB) domain regulates signaling pathways in both eukaryotes and prokaryotes. In this study, we analyze the evolutionary information embedded in genomic sequences to explore the diversity of signaling through the CNB domain and also how the CNB domain elicits a cellular response upon binding to cAMP.
Identification and classification of CNB domains in Global Ocean Sampling and other protein sequences reveals that they typically are fused to a wide variety of functional domains. CNB domains have undergone major sequence variation during evolution. In particular, the sequence motif that anchors the cAMP phosphate (termed the PBC motif) is strikingly different in some families. This variation may contribute to ligand specificity inasmuch as members of the prokaryotic cooA family, for example, harbor a CNB domain that contains a non-canonical PBC motif and that binds a heme ligand in the cAMP binding pocket. Statistical comparison of the functional constraints imposed on the canonical and non-canonical PBC containing sequences reveals that a key arginine, which coordinates with the cAMP phosphate, has co-evolved with a glycine in a distal β2-β3 loop that allosterically couples cAMP binding to distal regulatory sites.
Our analysis suggests that CNB domains have evolved as a scaffold to sense a wide variety of second messenger signals. Based on sequence, structural and biochemical data, we propose a mechanism for allosteric regulation by CNB domains.
The cyclic nucleotide binding (CNB) domain is a conserved signaling module that has evolved to respond to second messenger signals such as cAMP and cGMP [1, 2]. The CNB domain is ubiquitous in eukaryotes and controls a variety of cellular functions in a cAMP/cGMP dependent manner. Some of the well characterized CNB domain containing families in eukaryotes include: the protein kinase A (PKA) regulatory subunit that regulates the activity of PKA [3, 4]; the guanine nucleotide exchange factor that regulates nucleotide exchange in small GTPases ; and the ion channels that regulate metal ion gating (reviewed in ).
CNB domains also occur in prokaryotes. The first characterized family containing a CNB domain in prokaryotes is the CAP (catabolite gene activator protein) family of transcriptional regulators  that contain a DNA binding helix-turn-helix (HTH) domain covalently linked to the CNB domain . This domain organization is important for CAP function as it couples cAMP binding functions of the CNB domain with DNA binding functions of the HTH domain . The CAP family is functionally diverse and, in addition to cAMP, responds to other exogenous signals, such as carbon monoxide (CO) and nitric oxide (NO) (reviewed in ). The cooA subfamily, for instance, responds to CO signals and binds a heme ligand in the cAMP binding pocket . Likewise, the CprK subfamily of transcriptional regulators binds to ortho-chlorophenolic compounds in the cAMP binding pocket .
Crystal structures of CNB domains from both eukaryotes and prokaryotes have been determined and their structural comparison reveals a conserved mode of cAMP recognition  and regulation (reviewed in ). CNB domains are characterized by an eight stranded beta barrel domain (beta subdomain)  that is conserved among all CNB domain containing proteins . A key structural region within the beta subdomain is the phosphate binding cassette (PBC) that anchors the phosphate group of cAMP . CNB domains also contain a helical subdomain (henceforth called alpha subdomain), which, unlike the beta subdomain, is more variable in sequence and structure. The helical subdomain is also a docking site for the catalytic subunit of PKA .
An emerging theme in CNB domain signaling is the allosteric control of CNB domain functions. In the PKA regulatory subunit, for instance, cAMP binding to the beta subdomain causes conformational changes in the distal alpha subdomain, thereby releasing its inhibitory interactions with the catalytic subunit . This propagation of the cAMP signal to distal regulatory sites was suggested to involve specific regions in the beta subdomain . Specifically, a loop connecting the β2 and β3 strands (β2-β3 loop) was shown to undergo large chemical shift changes upon binding to cAMP . While these and other studies have provided important insights into PKA allostery, it is not known whether this mode of regulation is unique to the PKA regulatory subunit or is conserved among other members of the CNB domain superfamily. Here, we address this question by extracting and analyzing the evolutionary information encoded within CNB domain containing sequences. Towards this end, we have identified nearly 7,700 CNB domain containing proteins, and classified them into 30 distinct families. A systematic comparison of these families reveals that the CNB domains recombine with a wide variety of functional domains to respond to diverse cellular signals. Statistical comparison of the evolutionary constraints imposed on CNB domain sequences reveals that the residues that anchor the phosphate group of cAMP (within the beta subdomain) have co-evolved with residues in the β2-β3 loop. Analyzing these residues in light of existing structural and biochemical data provides a model of allostery that is conserved through evolution.
In the following sections, we first describe the identification and classification of CNB domains to illustrate the diversity of this protein family, and later show how a comparative analysis of CNB domain sequences has provided insights into the evolution of allostery.
Results and discussion
Identification and classification of CNB domains in the public and Global Ocean Sampling data
Cyclic nucleotide binding domains in the National Center for Biotechnology Information's non-redundant amino acid database (NR) and Global Ocean Sampling (GOS) [19, 20] data were identified using a combination of psi-blast profiles and motif models (see Materials and methods). This resulted in nearly 5,241 significant hits in NR and 2,455 hits in the GOS data. Most of the identified sequences were multi-domain proteins in that they contained other functional domains covalently linked to the CNB domain. Because these functional domains play an important role in CNB domain functions, they were used as markers for annotation and classification (see below).
Classification of CNB domains in the public and GOS data
PBC consensus motif
cAMP dependent regulatory subunit that activates PKA
cGMP activated proteins that are typically attached to a kinase domain
A distinct group of PKGs in parasites that are also attached to kinase domains
CNB domains from metazoans and plants. These are attached to various functional domains such as PKs, PAS domains, PP2C like phosphatases and phospholipases
cAMP-dependent guanine nucleotide exchange factors. Typically attached to an amino-terminal DEP domain and a carboxy-terminal RasGEF domain
A distinct class of Epac's, also called Epac6, which contains a PDZ domain in between the CNB and RasGEF domain. Epac's of this class contain a non-canonical PBC
Potassium channels specific to plants. Most of them contain an Ankryin repeat carboxy-terminal to the CNB domain
CNB domains found in metazoans and fungi, usually occur in tandem like the PKA regulatory subunit and contain a carboxy-terminal F-box domain and leucine rich domain
cGMP-gated cation channels. Mostly present in metazoans
Potassium channels that contain a PAC motif (motif carboxy-terminal of PAS) amino-terminal of the trans-membrane segment. This subfamily also contains a non-canonical PBC
Likely HCN channels from the single celled eukaryote Tetrahymena thermophila. This subfamily is quite distinct from the HCN channels in higher eukaryotes
Other HCN channels in protozoans
Tandem CNB domains that are attached to an amino-terminal pyridine nucleotide-disulphide oxidoreductase domain
Bacterial CNBs that are attached to mechanosensitive ion channels
Bacterial CNBs that contain a HisK like ATPase, carboxy-terminal of the CNB domain
A distinct sub-group containing AAA-ATPase domains attached to the CNB domain. Several members of this group contain an ABC-transporter like transmembrane region. The PBC arginine (Arg209) is quite variable within this family
Nitrogen responsive regulatory protein that contains a DNA binding domain (HTH) carboxy-terminal of the CNB domain
Involved in nitrogen fixation and contains a HTH motif
Transcriptional regulators that are implicated in oxygen sensing
Transcriptional regulator that is implicated in the aerobic arginase reaction. Arginine is used as a source of energy in bacteria
Transcriptional regulators that act on the nir and nor operons to achieve expression under aerobic conditions
This group contains tandem CBS domain located carboxy-terminal of the CNB domain
Bacterial CNB domains that are attached to various functional domains such as CheY response regulators, Rhodanese homology domain, kinases and DNA binding domains
Transcriptional regulator that is implicated in the repression of the acetate operon (also known as glyoxylate bypass operon) in Escherichia coli and Salmonella typhimurium
Transcriptional regulator containing a HTH domain and implicated in the repression of the gluconate operon
Involved in the bacterial oxidative stress response
Functions as a transcriptional repressor of an arsenic resistance operon. Dissociates from DNA in the presence of the metal
Transcriptional regulation of the crp operon
Repressor of genes that activate the multiple antibiotic resistance and oxidative stress regulons
An autogenously regulated activator of asparagine synthetase A transcription in Escherichia coli
CNB domains are also prevalent in prokaryotes and some of the major groups include: the CRP family members (Marr, Arsr, AsnC, ICLR, GNTR) that contain a DNA binding domain covalently linked to the CNB domain; and a distinct class of DNA binding domain containing proteins (NnR, ArcR, Fnr and FixK) that are activated by second messenger signals such as NO, oxygen and heme . In addition, our analysis reveals several novel families (CBS, HisK and AAA ATPases) in prokaryotes that lack the DNA binding domain, but conserve other functional domains (Table 1) such as histidine kinases (HisKs), cystathionine beta synthase (CBS) domains and AAA ATPases (AAA_Atpases in Table 1).
Expansion of transcriptional regulators in the Global Ocean Sampling data
Most of the GOS sequences, as expected, are prokaryotic in origin since they belong to families that are exclusively prokaryotic (Table 1). In particular, the CAP/CRP family, which contains a DNA binding domain covalently linked to the CNB domain and is implicated in the transcriptional regulation of genes, is greatly expanded in the GOS data (Table 1). The expansion of this family in the GOS data suggests that transcriptional regulation of many genes in oceanic microorganisms may be controlled in a cAMP or cGMP dependent manner. Also, the diversity displayed by the GOS sequences in the CAP family suggests that this family may regulate a wide variety of operons, in addition to the well studied lac operon . In addition to the CAP family, the NtcA family (Table 1), which is involved in nitrogen fixing in cyanobacteria , is also expanded in the GOS data. More than half the GOS sequences fall into the 'Other_Bacterial' family (table 1), which is poorly characterized. This family is highly diverse and contains several distinct sub-families that are associated with functional domains such as Rhodanases, Chey response regulators and DUF domains (Table 1). Thus, GOS data greatly contribute to the diversity of the CNB superfamily and enable the use of statistical methods to understand how sequence divergence contributes to functional divergence (see below).
Diversity in prokaryotes
Until now, the primary function of CNB domains in prokaryotes was believed to be in the transcriptional regulation of genes. However, our analysis suggests that other cellular processes, such as ATP production, protein phosphorylation and NADH production, may also involve CNB domain functions (Table 1). Of particular interest is the CBS domain associated CNB domains. CBS domains are known to function as sensors of cellular energy levels in eukaryotes as they are activated by AMP and inhibited by ATP. They are also implicated in various hereditary diseases in humans . The function of CBS domains in prokaryotes, however, is poorly understood, although the crystal structure of a CBS domain from Thermotoga maritime has been determined as part of the structural genomics initiative . The occurrence of both a CBS domain and a CNB domain in the same open reading frame suggests that, in some bacteria, ATP levels may be regulated in a cAMP-dependent manner. Structurally characterizing the full-length protein (CBS + CNB domain) may shed light on this regulatory mechanism in prokaryotes.
Other novel domains in prokaryotes that are fused to CNB domains include the HisKs that are involved in bacterial two component signaling, and the AAA class of ATPases (AAA_Atpases in Table 1) that control a wide variety of cellular functions in both eukaryotes and prokaryotes .
A conserved core shared by the entire superfamily
Functional diversity of the CNB module: a common scaffold to sense diverse ligands
Sequence variation within the PBC region contributes to ligand specificity
Several families in prokaryotes conserve a non-canonical PBC motif. Some of these include the transcriptional regulators FixK, FnR, ArcR, NnR and ARSR (Figure 5b). Within the FixK, or cooA family, for instance, the observed sequence variation within the PBC region appears to contribute to ligand specificity inasmuch as the cooA family binds to a heme ligand in the cAMP binding pocket (Figure 5e). In the crystal structure of cooA, a conserved histidine, which occupies a position that is structurally analogous to Arg209 in PKA, coordinates with the heme and plays a key role in cooA activation . Likewise, in the crystal structure of the transcriptional regulator CrpK bound to chlorophenolacetic acid , a structurally analogous asparagine (Asn92) residue hydrogen bonds to chlorophenolacetic acid (Figure 5f).
Evolution of allostery in the CNB module
The ability of the CNB domain to bind to diverse ligands raises an important question: what features distinguish the cAMP binding families (ones that conserve a canonical PBC motif) from those that bind to other ligands? In order to address this question we used the CHAIN (Contrast Hierarchical Alignment and Interaction Network analysis) program, which quantifies the differences between two functionally divergent groups of sequences using statistical methods . Using this program, we identified sequence features that distinguish the canonical PBC motif containing CNB domains from those that lack the canonical PBC motif. Analyzing these features in light of existing structural and biochemical data provides a model for allosteric regulation, which is likely conserved in all cAMP binding modules.
Selective constraints distinguishing the canonical PBC containing sequences
Recent NMR studies on the PKA regulatory subunit had suggested a key role for the β2-β3 loop in coupling cAMP signals to distal regulatory sites . Specifically, the backbone amide of Gly169 was shown to undergo large chemical shift changes upon binding to cAMP. This change was proposed to alter the conformation of an adjacent aspartate (Asp170), the backbone of which forms an N-cap to the B/C-helix (Figure 6b). Because the B/C helix forms a docking site for the catalytic subunit, this coupling between the PBC and the B/C-helix (via the β2-β3 loop) was proposed to play a key role in PKA allostery . The co-conservation of Gly169 with Arg209 suggests that this allosteric coupling may have specifically evolved in CBDs that bind to cAMP. Notably, MARR-bacteria and ASNC-bacteria (Figure 6a) are two families that conserve Arg209 in the PBC, but lack Gly169 in the β2-β3 loop. These two families presumably may have evolved alternative mechanisms of regulation. Future studies will focus on delineating these mechanisms using a combination of computational and experimental techniques.
A global analysis of CNB domain containing sequences in the public and GOS data has provided novel insights into the evolution of CNB domain structure and function. Two evolutionary events appear to have contributed to CNB domain functional divergence, domain recombination and sequence variation. The sequence diversity observed within the PBC suggests that the CNB domain has evolved as a scaffold for not only binding cAMP, but also a wide variety of other ligands, many of which are yet to be characterized. Statistical comparison of the evolutionary constraints acting on the canonical PBC motif containing CNB domains with the non-canonical ones reveals that the residues in the PBC region have co-evolved with residues in the β2-β3 loop. Examining these constraints in light of structural and biochemical data provides a model of allosteric regulation, which is likely conserved in all cAMP binding modules. The results described in this study have implications for protein engineering and for the design of allosteric inhibitors.
Materials and methods
Identification of CNB domains
CNB domains in GOS and NR data were identified using a combination of psi-blast  and Gibbs motif sampling procedures . Psi-blast profiles and motif models were initially built using CNB domains of known structures. These models were then iteratively updated as distant members from NR and GOS data were identified. An e-value cutoff of 0.001 was used for psi-blast searches.
Classification of CNB domains in NR
CNB domains identified from NR (5,241 sequences) were multiply aligned using the CHAIN analysis program . The aligned sequences were clustered into families and sub-families using the clustering option in the CHAIN program and the SECATOR program . Families were annotated by identifying the functional domains linked to the CNB domain. The taxonomic origin of the sequences was also taken into account in the annotation processes. For instance, PKG-like CNB domains from parasitic organisms were annotated as 'PKG_parasites'. Functional domains were identified using rpsblast, which was run against a collection of conserved domains in CDD, Smart and Pfam  with an e-value cutoff of 0.0001.
Classification of Global Ocean Sampling CNB domain containing proteins
Because CNB domains in the GOS data displayed significant sequence similarity to known CNB domains, they were assigned to one of the 30 families by running them against 30 family specific blast profiles. The taxonomic assignment for the GOS sequences was likewise done based on their similarity to known NR sequences . Examination of the domain organization in individual families indicated that while the NR sequence contained both the CNB domain and functional domains, GOS sequences usually contained only the CNB domain. This presumably is due to the fragmentary nature of the GOS data. In any case, nearly all the CNB domain containing GOS sequences could be assigned to one of the 30 families based on the similarity within the CNB domain alone.
Visualization of phylogenetic trees
In order to visually examine the evolutionary relationship between the identified sequences, we first constructed a phylogentic tree of all the 7,696 CNB sequences. The resulting tree, however, was very complex and hard to interpret. Therefore, we decided to take an alternative approach where we depicted each family by a consensus sequence. The 30 consensus sequences, corresponding to each of the 30 families, were generated from multiple alignments of individual families. The neighbor joining algorithm as implemented in the Molecular Evolutionary Genetics and Analysis (MEGA) program  was used for tree construction and visualization. Bootstrap test was done using default settings in MEGA.
Measuring the evolutionary constraints imposed on CNB sequences
The evolutionary constraints imposed on CNB sequences were measured using the CHAIN program . In brief, the CHAIN program identifies co-conserved residues that distinguish two related sets of sequences (foreground and background) by measuring the degree to which aligned residue positions in the foreground set are shifted away from the corresponding position in the background set. Residue positions that are shifted the most (indicated by red histograms above the alignment) contribute to the functional divergence of the foreground set from the background set. In the current study, all the CNB sequences that contain the canonical PBC motif constitute the foreground set, while the ones that lack the canonical motif constitute the background set.
The sequence identifiers for the sequences used in alignments Figures 2, 5b and 6a are: 94370018|PDZ_GEF-mouse; 93138731|K-channel-plant; 9857982|FixK-bacteria; 6759981|Fnr-bacteria; 15675445|ArcR-bacteria; 17989331|NnR-bacteria; 68552962|CBS-bacteria; 15673985|Flp-bacteria; 56419292|ARSR-bacteria; 1942960|PKA-mouse; 37964177|PKG-seahare; 68076807|PKA-parasite; 76609590|Epac-cattle; 68402320|HCN-zebrafish; 89309052|channel_Tetrahymena; 87198326|Bact_Pyrredox; 22298372|channel_Bact; 76259471|HisK-bacteria; 106879720|AAA_Atpase-bacteria; 462748|NtcA-bacteria; 86610079|ICLR-bacteria; 71367866|GNTR-bacteria; 111225891|CRP-bacteria; 115352640|MARR-bacteria; 116183754|ASNC-bacteria; 1FT9|pdb|cooA-bacteria; 2D93|pdb|PDZ_GEF_human; 2H6B|pdb|CprK-human.
catabolite activator protein
cystathionine beta synthase
cyclic nucleotide binding
Global Ocean Sampling
National Center for Biotechnology Information's non-redundant amino acid database
phosphate binding cassette
We thank Doug Rusch at the Venter Institute and Alexander Kornev at the San Diego Supercomputer center for helpful discussions. We thank the Taylor Lab members for useful comments and Sventja in the Taylor Lab for help with the illustrations. This work was supported by funding from the National Institutes of Health grant IP01DK54441 to SST. Grants to AFN from the National Library of Medicine (LM06747) and the Division of General Medicine (GM078541) are also acknowledged. We gratefully acknowledge the US Department of Energy, Office of Science (DE-FG02-02ER63453), the Gordon and Betty Moore Foundation, and the J Craig Venter Science Foundation for funding the GOS expedition.
- Berman HM, Ten Eyck LF, Goodsell DS, Haste NM, Kornev A, Taylor SS: The cAMP binding domain: an ancient signaling module. Proc Natl Acad Sci USA. 2005, 102: 45-50. 10.1073/pnas.0408579102.PubMedPubMed CentralView ArticleGoogle Scholar
- Anantharaman V, Koonin EV, Aravind L: Regulatory potential, phyletic distribution and evolution of ancient, intracellular small-molecule-binding domains. J Mol Biol. 2001, 307: 1271-1292. 10.1006/jmbi.2001.4508.PubMedView ArticleGoogle Scholar
- Gill GN, Garren LD: Role of the receptor in the mechanism of action of adenosine 3':5'-cyclic monophosphate. Proc Natl Acad Sci USA. 1971, 68: 786-790. 10.1073/pnas.68.4.786.PubMedPubMed CentralView ArticleGoogle Scholar
- Taylor SS, Buechler JA, Yonemoto W: cAMP-dependent protein kinase: framework for a diverse family of regulatory enzymes. Annu Rev Biochem. 1990, 59: 971-1005. 10.1146/annurev.bi.59.070190.004543.PubMedView ArticleGoogle Scholar
- de Rooij J, Zwartkruis FJ, Verheijen MH, Cool RH, Nijman SM, Wittinghofer A, Bos JL: Epac is a Rap1 guanine-nucleotide-exchange factor directly activated by cyclic AMP. Nature. 1998, 396: 474-477. 10.1038/24884.PubMedView ArticleGoogle Scholar
- Kaupp UB, Seifert R: Cyclic nucleotide-gated ion channels. Physiol Rev. 2002, 82: 769-824.PubMedView ArticleGoogle Scholar
- Weber IT, Takio K, Titani K, Steitz TA: The cAMP-binding domains of the regulatory subunit of cAMP-dependent protein kinase and the catabolite gene activator protein are homologous. Proc Natl Acad Sci USA. 1982, 79: 7679-7683. 10.1073/pnas.79.24.7679.PubMedPubMed CentralView ArticleGoogle Scholar
- McKay DB, Steitz TA: Structure of catabolite gene activator protein at 2.9 A resolution suggests binding to left-handed B-DNA. Nature. 1981, 290: 744-749. 10.1038/290744a0.PubMedView ArticleGoogle Scholar
- Benoff B, Yang H, Lawson CL, Parkinson G, Liu J, Blatter E, Ebright YW, Berman HM, Ebright RH: Structural basis of transcription activation: the CAP-alpha CTD-DNA complex. Science. 2002, 297: 1562-1566. 10.1126/science.1076376.PubMedView ArticleGoogle Scholar
- Korner H, Sofia HJ, Zumft WG: Phylogeny of the bacterial superfamily of Crp-Fnr transcription regulators: exploiting the metabolic spectrum by controlling alternative gene programs. FEMS Microbiol Rev. 2003, 27: 559-592. 10.1016/S0168-6445(03)00066-4.PubMedView ArticleGoogle Scholar
- Lanzilotta WN, Schuller DJ, Thorsteinsson MV, Kerby RL, Roberts GP, Poulos TL: Structure of the CO sensing transcription activator CooA. Nat Struct Biol. 2000, 7: 876-880. 10.1038/82820.PubMedView ArticleGoogle Scholar
- Joyce MG, Levy C, Gabor K, Pop SM, Biehl BD, Doukov TI, Ryter JM, Mazon H, Smidt H, van den Heuvel RH, et al: CprK crystal structures reveal mechanism for transcriptional control of halorespiration. J Biol Chem. 2006, 281: 28318-28325. 10.1074/jbc.M602654200.PubMedView ArticleGoogle Scholar
- Rehmann H, Wittinghofer A, Bos JL: Capturing cyclic nucleotides in action: snapshots from crystallographic studies. Nat Rev Mol Cell Biol. 2007, 8: 63-73. 10.1038/nrm2082.PubMedView ArticleGoogle Scholar
- Su Y, Dostmann WRG, Herberg FW, Durick K, Xuong NH, Ten Eyck LF, Taylor SS, Varughese KI: Regulatory (RIa) subunit of protein kinase a: structure of deletion mutant with cAMP binding domains. Science. 1995, 269: 807-819. 10.1126/science.7638597.PubMedView ArticleGoogle Scholar
- Diller TC, Madhusudan , Xuong NH, Taylor SS: Molecular basis for regulatory subunit diversity in cAMP-dependent protein kinase: crystal structure of the type II beta regulatory subunit. Structure. 2001, 9: 73-82. 10.1016/S0969-2126(00)00556-6.PubMedView ArticleGoogle Scholar
- Kim C, Xuong NH, Taylor SS: Crystal structure of a complex between the catalytic and regulatory (RIalpha) subunits of PKA. Science. 2005, 307: 690-696. 10.1126/science.1104607.PubMedView ArticleGoogle Scholar
- Kim C, Cheng CY, Saldanha SA, Taylor SS: PKA-I holoenzyme structure reveals a mechanism for cAMP-dependent activation. Cell. 2007, 130: 1032-1043. 10.1016/j.cell.2007.07.018.PubMedView ArticleGoogle Scholar
- Das R, Esposito V, Abu-Abed M, Anand GS, Taylor SS, Melacini G: cAMP activation of PKA defines an ancient signaling mechanism. Proc Natl Acad Sci USA. 2007, 104: 93-98. 10.1073/pnas.0609033103.PubMedPubMed CentralView ArticleGoogle Scholar
- Yooseph S, Sutton G, Rusch DB, Halpern AL, Williamson SJ, Remington K, Eisen JA, Heidelberg KB, Manning G, Li W, et al: The Sorcerer II Global Ocean Sampling expedition: expanding the universe of protein families. PLoS Biol. 2007, 5: e16-10.1371/journal.pbio.0050016.PubMedPubMed CentralView ArticleGoogle Scholar
- Rusch DB, Halpern AL, Sutton G, Heidelberg KB, Williamson S, Yooseph S, Wu D, Eisen JA, Hoffman JM, Remington K, et al: The Sorcerer II Global Ocean Sampling expedition: northwest Atlantic through eastern tropical Pacific. PLoS Biol. 2007, 5: e77-10.1371/journal.pbio.0050077.PubMedPubMed CentralView ArticleGoogle Scholar
- Bridges D, Fraser ME, Moorhead GB: Cyclic nucleotide binding proteins in the Arabidopsis thaliana and Oryza sativa genomes. BMC Bioinformatics. 2005, 6: 6-10.1186/1471-2105-6-6.PubMedPubMed CentralView ArticleGoogle Scholar
- Okamoto K, Hara S, Bhasin R, Freundlich M: Evidence in vivo for autogenous control of the cyclic AMP receptor protein gene (crp) in Escherichia coli by divergent RNA. J Bacteriol. 1988, 170: 5076-5079.PubMedPubMed CentralGoogle Scholar
- Bradley RL, Reddy KJ: Cloning, sequencing, and regulation of the global nitrogen regulator gene ntcA in the unicellular diazotrophic cyanobacterium Cyanothece sp. strain BH68K. J Bacteriol. 1997, 179: 4407-4410.PubMedPubMed CentralGoogle Scholar
- Scott JW, Hawley SA, Green KA, Anis M, Stewart G, Scullion GA, Norman DG, Hardie DG: CBS domains form energy-sensing modules whose binding of adenosine ligands is disrupted by disease mutations. J Clin Invest. 2004, 113: 274-284. 10.1172/JCI200419874.PubMedPubMed CentralView ArticleGoogle Scholar
- Miller MD, Schwarzenbacher R, von Delft F, Abdubek P, Ambing E, Biorac T, Brinen LS, Canaves JM, Cambell J, Chiu HJ, et al: Crystal structure of a tandem cystathionine-beta-synthase (CBS) domain protein (TM0935) from Thermotoga maritima at 1.87 A resolution. Proteins. 2004, 57: 213-217. 10.1002/prot.20024.PubMedView ArticleGoogle Scholar
- Neuwald AF, Aravind L, Spouge JL, Koonin EV: AAA+: A class of chaperone-like ATPases associated with the assembly, operation, and disassembly of protein complexes. Genome Res. 1999, 9: 27-43.PubMedGoogle Scholar
- Bubis J, Neitzel JJ, Saraswat LD, Taylor SS: A point mutation abolishes binding of cAMP to site A in the regulatory subunit of cAMP-dependent protein kinase. J Biol Chem. 1988, 263: 9668-9673.PubMedGoogle Scholar
- Neuwald AF: The CHAIN program: forging evolutionary links to underlying mechanisms. Trends Biochem Sci. 2007, 32: 487-493. 10.1016/j.tibs.2007.08.009.PubMedView ArticleGoogle Scholar
- Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997, 25: 3389-3402. 10.1093/nar/25.17.3389.PubMedPubMed CentralView ArticleGoogle Scholar
- Neuwald AF, Liu JS: Gapped alignment of protein sequence motifs through Monte Carlo optimization of a hidden Markov model. BMC Bioinformatics. 2004, 5: 157-10.1186/1471-2105-5-157.PubMedPubMed CentralView ArticleGoogle Scholar
- Wicker N, Perrin GR, Thierry JC, Poch O: Secator: a program for inferring protein subfamilies from phylogenetic trees. Mol Biol Evol. 2001, 18: 1435-1441.PubMedView ArticleGoogle Scholar
- Marchler-Bauer A, Panchenko AR, Shoemaker BA, Thiessen PA, Geer LY, Bryant SH: CDD: a database of conserved domain alignments with links to domain three-dimensional structure. Nucleic Acids Res. 2002, 30: 281-283. 10.1093/nar/30.1.281.PubMedPubMed CentralView ArticleGoogle Scholar
- Kumar S, Tamura K, Nei M: MEGA: Molecular Evolutionary Genetics Analysis software for microcomputers. Comput Appl Biosci. 1994, 10: 189-191.PubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.