Molecular dynamics of the FixJ receiver domain: movement of the beta4-alpha4 loop correlates with the in and out flip of Phe101

Thirty years of molecular dynamics simulations on posttranslational modifications of proteins

Posttranslational modifications (PTMs) are an integral component to how cells respond to perturbation. While experimental advances have enabled improved PTM identification capabilities, the same throughput for characterizing how structural changes caused by PTMs equate to altered physiological function has not been maintained. In this Perspective, we cover the history of computational modeling and molecular dynamics simulations which have characterized the structural implications of PTMs. We distinguish results from different molecular dynamics studies based upon the timescales simulated and analysis approaches used for PTM characterization. Lastly, we offer insights into how opportunities for modern research efforts on in silico PTM characterization may proceed given current state-of-the-art computing capabilities and methodological advancements.

Oxygen regulatory mechanisms of nitrogen fixation in rhizobia

Rhizobia are α- and β-proteobacteria that form a symbiotic partnership with legumes, fixing atmospheric dinitrogen to ammonia and providing it to the plant. Oxygen regulation is key in this symbiosis. Fixation is performed by an oxygen-intolerant nitrogenase enzyme but requires respiration to meet its high energy demands. To satisfy these opposing constraints the symbiotic partners cooperate intimately, employing a variety of mechanisms to regulate and respond to oxygen concentration. During symbiosis rhizobia undergo significant changes in gene expression to differentiate into nitrogen-fixing bacteroids. Legumes host these bacteroids in specialized root organs called nodules. These generate a near-anoxic environment using an oxygen diffusion barrier, oxygen-binding leghemoglobin and control of mitochondria localization. Rhizobia sense oxygen using multiple interconnected systems which enable a finely-tuned response to the wide range of oxygen concentrations they experience when transitioning from soil to nodules. The oxygen-sensing FixL-FixJ and hybrid FixL-FxkR two-component systems activate at relatively high oxygen concentration and regulate fixK transcription. FixK activates the fixNOQP and fixGHIS operons producing a high-affinity terminal oxidase required for bacterial respiration in the microaerobic nodule. Additionally or alternatively, some rhizobia regulate expression of these operons by FnrN, an FNR-like oxygen-sensing protein. The final stage of symbiotic establishment is activated by the NifA protein, regulated by oxygen at both the transcriptional and protein level. A cross-species comparison of these systems highlights differences in their roles and interconnections but reveals common regulatory patterns and themes. Future work is needed to establish the complete regulon of these systems and identify other regulatory signals.

Small phospho-donors phosphorylate MorR without inducing protein conformational changes

Phosphorylation is an essential mechanism of protein control and plays an important role in biology. The two-component system (TCS) is a bacterial regulation mechanism mediated by a response regulator (RR) protein and a kinase protein, which synchronize the regulatory circuit according to the environment. Phosphorylation is a key element in TCS function as it controls RR activity. In the present study, we characterize the behavior of MorR, an RR associated with Mo homeostasis, upon acetylphosphate and phosphoramidate treatment in vitro. Our results show that MorR was phosphorylated by both phospho-donors. Fluorescence experiments showed that MorR tryptophan emission is quenched by phosphoramidate. Furthermore, theoretical and computational results demonstrate that phosphorylation by phosphoramidate is more favorable than that by acetylphosphate. In conclusion, phosphorylated MorR is a monomeric protein and phosphorylation does not appear to induce observable conformational changes in the protein structure.

Production, characterization, and assessment of a stable analog of the response regulator CheY-phosphate from Thermotoga maritima

Phosphorylation of CheY promotes association with the flagellar motor and ultimately controls the directional bias of the motor. However, biochemical studies of activated CheY-phosphate have been challenging due to the rapid hydrolysis of the aspartyl-phosphate in vitro. An inert analog of Tm CheY-phosphate, phosphono-CheY, was synthesized by chemical modification and purified by cation-exchange chromatography. Changes in HPLC retention times, chemical assays for phosphate and free thiol, and mass spectrometry experiments demonstrate modification of Cys54 with a phosphonomethyl group. Additionally, a crystal structure showed electron density for the phosphonomethyl group at Cys54, consistent with a modification at that position. Subsequent biochemical experiments confirmed that protein crystals were phosphono-CheY. Isothermal titration calorimetry and fluorescence polarization binding assays demonstrated that phosphono-CheY bound a peptide derived from FliM, a native partner of CheY-phosphate, with a dissociation constant of ∼29 µM, at least sixfold more tightly than unmodified CheY. Taken together these results suggest that Tm phosphono-CheY is a useful and unique analog of Tm CheY-phosphate.

Evidence against the "Y-T coupling" mechanism of activation in the response regulator NtrC

The dominant theory on the mechanism of response regulators activation in two-component bacterial signaling systems is the "Y-T coupling" mechanism, wherein the χ1 rotameric state of a highly conserved aromatic residue correlates with the activation of the protein via structural rearrangements coupled to a conserved tyrosine. In this paper, we present evidence that, in the receiver domain of the response regulator nitrogen regulatory protein C (NtrC(R)), the interconversion of this tyrosine (Y101) between its rotameric states is actually faster than the rate of inactive/active conversion and is not correlated to the activation process. Data gathered from NMR relaxation dispersion experiments show that a subset of residues surrounding the conserved tyrosine sense a process that is occurring at a faster rate than the inactive/active conformational transition. We show that this process is related to χ1 rotamer exchange of Y101 and that mutation of this aromatic residue to a leucine eliminated this second faster process without affecting activation. Computational simulations of NtrC(R) in its active conformation further demonstrate that the rotameric state of Y101 is uncorrelated with the global conformational transition during activation. Moreover, the tyrosine does not appear to be involved in the stabilization of the active form upon phosphorylation and is not essential in propagating the signal downstream for ATPase activity of the central domain. Our data provide experimental evidence against the generally accepted "Y-T coupling" mechanism of activation in NtrC(R).

Novel redox-sensing modules: accessory protein- and nucleic acid-mediated signaling

SIGNIFICANCE

Organisms have evolved both enzymatic and nonenzymatic pathways to prevent oxidative damage to essential macromolecules, including proteins and nucleic acids. Pathways modulated by different protein-based sensory and regulatory modules ensure a rapid and appropriate response.

RECENT ADVANCES

In contrast to classical two-component systems that possess internal sensory and regulatory modules, an accessory protein-dependent redox-signaling system has been recently characterized in bacteria. This system senses extracellular iron-mediated oxidative stress signals via an extracellularly located protein (HbpS). In vivo and in vitro studies allowed the elucidation of molecular mechanisms governing this system. Moreover, recent studies show that nucleic acids may also participate in redox-signaling during antioxidative stress response.

CRITICAL ISSUES

Research for novel redox-signaling systems is often focused on known types of sensory and regulatory modules. It is also often considered that the oxidative attack of macromolecules, leading to modification and degradation processes, is the final step during oxidative stress. However, recent studies have demonstrated that oxidatively modified macromolecules can be intermediary states in the process of redox-signaling.

FUTURE DIRECTIONS

Analyses of adjacent regions of genes encoding for known sensory and regulatory modules can identify potential accessory modules that may increase the complexity of sensing systems. Despite the fact that the involvement of DNA-mediated signaling in the modulation of one bacterial regulator protein has been analyzed in detail, further studies are necessary to identify additional regulators. Given the role of DNA in oxidative-stress response, it is tempting to hypothesize that RNA modules may also mediate redox-signaling.

Detection of allosteric signal transmission by information-theoretic analysis of protein dynamics

Allostery offers a highly specific way to modulate protein function. Therefore, understanding this mechanism is of increasing interest for protein science and drug discovery. However, allosteric signal transmission is difficult to detect experimentally and to model because it is often mediated by local structural changes propagating along multiple pathways. To address this, we developed a method to identify communication pathways by an information-theoretical analysis of molecular dynamics simulations. Signal propagation was described as information exchange through a network of correlated local motions, modeled as transitions between canonical states of protein fragments. The method was used to describe allostery in two-component regulatory systems. In particular, the transmission from the allosteric site to the signaling surface of the receiver domain NtrC was shown to be mediated by a layer of hub residues. The location of hubs preferentially connected to the allosteric site was found in close agreement with key residues experimentally identified as involved in the signal transmission. The comparison with the networks of the homologues CheY and FixJ highlighted similarities in their dynamics. In particular, we showed that a preorganized network of fragment connections between the allosteric and functional sites exists already in the inactive state of all three proteins.

EVOLUTION OF A CATALYTICALLY EFFECTIVE MODEL ENZYME: THE IMPORTANCE OF TUNED CONFORMATIONAL FLUCTUATIONS

Possible causal connections between the dynamics of a thermally fluctuating model enzyme molecule and catalysis are explored. The model is motivated by observations from experiment and simulation that amino acid residues residing in different enzymatic domains may show markedly different degrees of conformational freedom. Consequently, we are interested in the catalytic efficacy of an enzyme as a function of long-range many-atom cooperative effects resulting from strong, moderate, and weak interactions between enzymatic residues. Here we show and quantify through molecular dynamics simulations how the number and distribution of these interactions affects an enzyme's conformational fluctuation dynamics and its effectiveness as a catalyst. For any given distribution of "stiff" and "loose" enzymatic domains, catalytic fitness is defined as the number of chemical events — specifically the number of times a catalytic residue and substrate surmount a chemical reaction barrier — during molecular dynamics simulation. Through mutation, recombination, and a selection procedure following the ideas of Darwinian evolution, a genetic algorithm drives a population of enzyme molecules to greater catalytic fitness by modifying the mix of stiff and loose interactions. Approximately 30,000 different enzyme molecules are generated by the genetic algorithm — each with a unique number and distribution of strong, moderate, and weak inter-residue interactions. While the catalytically least fit enzyme exhibits 16 chemical events, the fittest boasts 253. That point mutations far from the active-site chemistry in the fittest enzyme have a strong effect on the number of chemical events suggests that catalysis depends, in part, on long-range many-atom globally correlated dynamical fluctuations.

Structural basis of response regulator dephosphorylation by Rap phosphatases

Bacterial Rap family proteins have been most extensively studied in Bacillus subtilis, where they regulate activities including sporulation, genetic competence, antibiotic expression, and the movement of the ICEBs1 transposon. One subset of Rap proteins consists of phosphatases that control B. subtilis and B. anthracis sporulation by dephosphorylating the response regulator Spo0F. The mechanistic basis of Rap phosphatase activity was unknown. Here we present the RapH-Spo0F X-ray crystal structure, which shows that Rap proteins consist of a 3-helix bundle and a tetratricopeptide repeat domain. Extensive biochemical and genetic functional studies reveal the importance of the observed RapH-Spo0F interactions, including the catalytic role of a glutamine in the RapH 3-helix bundle that inserts into the Spo0F active site. We show that in addition to dephosphorylating Spo0F, RapH can antagonize sporulation by sterically blocking phosphoryl transfer to and from Spo0F. Our structure-function analysis of the RapH-Spo0F interaction identified Rap protein residues critical for Spo0F phosphatase activity. This information enabled us to assign Spo0F phosphatase activity to a Rap protein based on sequence alone, which was not previously possible. Finally, as the ultimate test of our newfound understanding of the structural requirements for Rap phosphatase function, a non-phosphatase Rap protein that inhibits the binding of the response regulator ComA to DNA was rationally engineered to dephosphorylate Spo0F. In addition to revealing the mechanistic basis of response regulator dephosphorylation by Rap proteins, our studies support the previously proposed T-loop-Y allostery model of receiver domain regulation that restricts the aromatic "switch" residue to an internal position when the β4-α4 loop adopts an active-site proximal conformation.

Insights into correlated motions and long-range interactions in CheY derived from molecular dynamics simulations

CheY is a response regulator protein involved in bacterial chemotaxis. Much is known about its active and inactive conformations, but little is known about the mechanisms underlying long-range interactions or correlated motions. To investigate these events, molecular dynamics simulations were performed on the unphosphorylated, inactive structure from Salmonella typhimurium and the CheY-BeF(3)(-) active mimic structure (with BeF(3)(-) removed) from Escherichia coli. Simulations utilized both sequences in each conformation to discriminate sequence- and structure-specific behavior. The previously identified conformational differences between the inactive and active conformations of the strand-4-helix-4 loop, which are present in these simulations, arise from the structural, and not the sequence, differences. The simulations identify previously unreported structure-specific flexibility features in this loop and sequence-specific flexibility features in other regions of the protein. Both structure- and sequence-specific long-range interactions are observed in the active and inactive ensembles. In the inactive ensemble, two distinct mechanisms based on Thr-87 or Ile-95 rotameric forms, are observed for the previously identified g+ and g- rotamer sampling by Tyr-106. These molecular dynamics simulations have thus identified both sequence- and structure-specific differences in flexibility, long-range interactions, and rotameric form of key residues. Potential biological consequences of differential flexibility and long-range correlated motion are discussed.

Reconciling the "old" and "new" views of protein allostery: a molecular simulation study of chemotaxis Y protein (CheY)

A combination of thirty-two 10-ns-scale molecular dynamics simulations were used to explore the coupling between conformational transition and phosphorylation in the bacteria chemotaxis Y protein (CheY), as a simple but representative example of protein allostery. Results from these simulations support an activation mechanism in which the beta4-alpha4 loop, at least partially, gates the isomerization of Tyr106. The roles of phosphorylation and the conserved Thr87 are deemed indirect in that they stabilize the active configuration of the beta4-alpha4 loop. The indirect role of the activation event (phosphorylation) and/or conserved residues in stabilizing, rather than causing, specific conformational transition is likely a feature in many signaling systems. The current analysis of CheY also helps to make clear that neither the "old" (induced fit) nor the "new" (population shift) views for protein allostery are complete, because they emphasize the kinetic (mechanistic) and thermodynamic aspects of allosteric transitions, respectively. In this regard, an issue that warrants further analysis concerns the interplay of concerted collective motion and sequential local structural changes in modulating cooperativity between distant sites in biomolecules.

Molecular dynamic simulations of the N-terminal receiver domain of NtrC reveal intrinsic conformational flexibility in the inactive state

The N-terminal receiver domain of NtrC is the molecular switch in the two-component signal transduction. It is the first protein where structures of both the active (phosphyroylated) and inactive (unphosphyroylated) states are determined experimentally. Phosphorylation of the NtrC at the active site induces large structural change. NMR experiments suggested that the wild type unphosphorylated NtrC adopts both the active and the inactive conformations and the phosphorylation stabilizes the active conformations. We applied free (unconstrained) molecular dynamic (MD) simulation to examine the intrinsic flexibilities and stabilities of the NtrC receiver domain in both the active and inactive conformations. Molecular dynamic simulations showed that the inactive state of NtrC receiver domain is more flexible than the active state. There were large movements in helix 4 and loop beta3-alpha3 which coincide with major structural differences between the inactive and active states. We observed large root-mean-square deviations from the initial starting structure and the large root-mean-square fluctuations during MD simulation for the inactive state. We then investigated the activation pathway with Targeted MD simulation. We show that the intrinsic flexibility in the loop beta3-alpha3 plays an important role in triggering the conformational change. Phosphorylation at the active site may serve to stabilize the conformational change. These results together suggest that the unphosphorylated NtrC receiver domain could be involved in a conformational equilibrium between two different states.

Conformational changes in protein loops and helices induced by post-translational phosphorylation

Post-translational phosphorylation is a ubiquitous mechanism for modulating protein activity and protein-protein interactions. In this work, we examine how phosphorylation can modulate the conformation of a protein by changing the energy landscape. We present a molecular mechanics method in which we phosphorylate proteins in silico and then predict how the conformation of the protein will change in response to phosphorylation. We apply this method to a test set comprised of proteins with both phosphorylated and non-phosphorylated crystal structures, and demonstrate that it is possible to predict localized phosphorylation-induced conformational changes, or the absence of conformational changes, with near-atomic accuracy in most cases. Examples of proteins used for testing our methods include kinases and prokaryotic response regulators. Through a detailed case study of cyclin-dependent kinase 2, we also illustrate how the computational methods can be used to provide new understanding of how phosphorylation drives conformational change, why substituting Glu or Asp for a phosphorylated amino acid does not always mimic the effects of phosphorylation, and how a phosphatase can “capture” a phosphorylated amino acid. This work illustrates how computational methods can be used to elucidate principles and mechanisms of post-translational phosphorylation, which can ultimately help to bridge the gap between the number of known sites of phosphorylation and the number of structures of phosphorylated proteins.

Many proteins are chemically modified after they are synthesized in the cell. These post-translational modifications can modulate the ability of a protein to perform chemical reactions and to interact with other proteins. At the cellular level, for example, these chemical modifications are critical for allowing the cell to respond to its environment and control its division. One of the most common mechanisms by which proteins can be modified is by phosphorylation—the addition of a phosphate group to an amino acid side chain of the protein. Thousands of proteins are known to be modified by phosphorylation, but only for a small minority of these do we have any detailed understanding of how the chemical modification regulates the function of the protein. The authors describe a computational method that can make testable predictions about the structural changes that occur in a protein induced by post-translational phosphorylation. Their results show that the method can produce structural models of the phosphorylated proteins with near-atomic accuracy, and provide insight into the energetics of conformational switches driven by phosphorylation. As such, the computational method complements experiments aimed at understanding the mechanisms of protein regulation by phosphorylation.

Crystal structures of beryllium fluoride-free and beryllium fluoride-bound CheY in complex with the conserved C-terminal peptide of CheZ reveal dual binding modes specific to CheY conformation

Chemotaxis, the environment-specific swimming behavior of a bacterial cell is controlled by flagellar rotation. The steady-state level of the phosphorylated or activated form of the response regulator CheY dictates the direction of flagellar rotation. CheY phosphorylation is regulated by a fine equilibrium of three phosphotransfer activities: phosphorylation by the kinase CheA, its auto-dephosphorylation and dephosphorylation by its phosphatase CheZ. Efficient dephosphorylation of CheY by CheZ requires two spatially distinct protein-protein contacts: tethering of the two proteins to each other and formation of an active site for dephosphorylation. The former involves interaction of phosphorylated CheY with the small highly conserved C-terminal helix of CheZ (CheZ(C)), an indispensable structural component of the functional CheZ protein. To understand how the CheZ(C) helix, representing less than 10% of the full-length protein, ascertains molecular specificity of binding to CheY, we have determined crystal structures of CheY in complex with a synthetic peptide corresponding to 15 C-terminal residues of CheZ (CheZ(200-214)) at resolutions ranging from 2.0 A to 2.3A. These structures provide a detailed view of the CheZ(C) peptide interaction both in the presence and absence of the phosphoryl analog, BeF3-. Our studies reveal that two different modes of binding the CheZ(200-214) peptide are dictated by the conformational state of CheY in the complex. Our structures suggest that the CheZ(C) helix binds to a "meta-active" conformation of inactive CheY and it does so in an orientation that is distinct from the one in which it binds activated CheY. Our dual binding mode hypothesis provides implications for reverse information flow in CheY and extends previous observations on inherent resilience in CheY-like signaling domains.

A Circadian Timing Mechanism in the Cyanobacteria

Cyanobacteria such as Synechococcus elongatus PCC 7942, Thermosynechococcus elongatus BP-1, and Synechocystis species strain PCC 6803 have an endogenous timing mechanism that can generate and maintain a 24 h (circadian) periodicity to global (whole genome) gene expression patterns. This rhythmicity extends to many other physiological functions, including chromosome compaction. These rhythmic patterns seem to reflect the periodicity of availability of the primary energy source for these photoautotrophic organisms, the Sun. Presumably, eons of environmentally derived rhythmicity--light/dark cycles--have simply been mechanistically incorporated into the regulatory networks of these cyanobacteria. Genetic and biochemical experimentation over the last 15 years has identified many key components of the primary timing mechanism that generates rhythmicity, the input pathways that synchronize endogenous rhythms to exogenous rhythms, and the output pathways that transduce temporal information from the timekeeper to the regulators of gene expression and function. Amazingly, the primary timing mechanism has evidently been extracted from S. elongatus PCC 7942 and can also keep time in vitro. Mixing the circadian clock proteins KaiA, KaiB, and KaiC from S. elongatus PCC 7942 in vitro and adding ATP results in a circadian rhythm in the KaiC protein phosphorylation state. Nonetheless, many questions still loom regarding how this circadian clock mechanism works, how it communicates with the environment and how it regulates temporal patterns of gene expression. Many details regarding structure and function of the individual clock-related proteins are provided here as a basis to discuss these questions. A strong, data-intensive foundation has been developed to support the working model for the cyanobacterial circadian regulatory system. The eventual addition to that model of the metabolic parameters participating in the command and control of this circadian global regulatory system will ultimately allow a fascinating look into whole-cell physiology and metabolism and the consequential organization of global gene expression patterns.

Direct determination of reaction paths and stationary points on potential of mean force surfaces

A simulation approach is introduced for directly determining reaction paths and stationary points on potential of mean force (PMF) surfaces associated with molecular events that occur in complex environments. The nudged elastic band approach was employed to search for steepest descent paths on the PMF surface using the relevant PMF derivatives from a series of local simulations. The steepest descent path on the PMF surface corresponds to the minimum PMF path (i.e. the path with the lowest local PMF barrier), which contains important information about stationary points (e.g. saddle points) on the PMF surface, which in turn can provide useful insights into the thermodynamics and kinetics for the process of interest. By working with the PMF defined in a low-dimensional sub-space rather than a potential energy function of full molecular dimensionality, the main features of the process under study are concisely represented and the orthogonal degrees of freedom are adequately sampled with the appropriate canonical distribution at the desired temperature (e.g. 300 K). Therefore, minimum PMF paths carry statistically meaningful mechanistic information and are complementary to reaction paths of full molecular dimensionality proposed in previous studies. The NEB based path optimization method is direct in the sense that no information regarding the global PMF surface is necessary for the determination of the local reaction path and stationary points along this path. Since only low-dimensional quantities (paths) are searched for, the PMF-path method is expected to scale better in terms of dimension of the PMF sub-space than those aims to fully explore multi-dimensional PMF surfaces. Test applications on simple molecular systems, the alanine di-peptide in vacuum and in solution and a micro-solvated proton-wire, indicate that reliable PMF paths can be determined for both conformational isomerization and chemical reaction processes. However, highly accurate PMF derivatives are required for determining more quantitative observables, such as the free energy profile along the minimum PMF path. Therefore, effective numerical algorithms for calculating local PMF derivatives and systematic protocols for defining the relevant sub-space are the main focus in the near future. Finally, we emphasize that the minimum PMF path defined here includes thermal (e.g. entropic) effects associated with the orthogonal degrees of freedom, but finite kinetic energies associated with the PMF degrees of freedom are not included; this can be improved by adopting a different definition of the reaction path, such as the maximum flux path, on the PMF surface, or thermally sampling all degrees of freedom orthogonal to the one-dimensional path.

Competition between Intramolecular Hydrogen Bonds and Solvation in Phosphorylated Peptides: Simulations with Explicit and Implicit Solvent

The atomic-level mechanisms of protein regulation by post-translational phosphorylation remain poorly understood, except in a few well-studied systems. Molecular mechanics simulations can in principle be used to help understand and predict the effects of protein phosphorylation, but the accuracy of the results will of course depend on the quality of the force field parameters for the phosphorylated residues as well as the quality of the solvent model. The phosphorylated residues typically carry a -2 charge at physiological pH; however, the effects of phosphorylation can sometimes be mimicked by substituting Asp or Glu for the phosphorylated residue. Here we examine the suitability of explicit and implicit solvent models for simulating phospho-serine in both the -1 and -2 charge states. Specifically, we simulate a capped phosphorylated peptide, Ace-Gly-Ser-pSer-Ser-Nme, and compare the results to each other and to experimental observables from an NMR experiment. The first major conclusion is that explicit water models (TIP3P, TIP4P and SPC/E) and a Generalized Born implicit solvent model provide reasonable agreement with the experimental observables, given appropriate partial charges for the phosphate group. The Generalized Born results, however, show greater hydrogen bonding propensity than the explicit solvent results. Distance dependent dielectric treatments perform poorly. The second major conclusion is that many ensemble-averaged properties obtained for the phosphopeptide in the -1 and -2 charge states are strikingly similar; the -1 species has a slightly higher propensity to form internal hydrogen bonds. All of the results can be rationalized by quantifying the strength of the P-O/H-N hydrogen bond, which depends on a sensitive balance between strongly favorable charge/dipole and dipole/dipole interactions and strongly unfavorable desolvation.

A Cyanobacterial Circadian Timing Mechanism

Cyanobacteria such as Synechococcus elongatus PCC 7942 exhibit 24-h rhythms of gene expression that are controlled by an endogenous circadian clock that is mechanistically distinct from those described for diverse eukaryotes. Genetic and biochemical experiments over the past decade have identified key components of the circadian oscillator, input pathways that synchronize the clock with the daily environment, and output pathways that relay temporal information to downstream genes. The mechanism of the cyanobacterial circadian clock that is emerging is based principally on the assembly and disassembly of a large complex at whose heart are the proteins KaiA, KaiB, and KaiC. Signal transduction pathways that feed into and out of the clock employ protein domains that are similar to those in two-component regulatory systems of bacteria.