Association and dissociation kinetics for CheY interacting with the P2 domain of CheA

Strategies adopted by gastric pathogen Helicobacter pylori for a mature biofilm formation: Antimicrobial peptides as a visionary treatment

Enormous efforts in recent past two decades to eradicate the pathogen that has been prevalent in half of the world's population have been problematic. The biofilm formed by Helicobacter pylori provides resistance towards innate immune cells, various combinatorial antibiotics, and human antimicrobial peptides, despite the fact that these all are potent enough to eradicate it in vitro. Biofilm provides the opportunity to secrete various virulence factors that strengthen the interaction between host and pathogen helping in evading the innate immune system and ultimately leading to persistence. To our knowledge, this review is the first of its kind to explain briefly the journey of H. pylori starting with the chemotaxis, the mechanism for selecting the site for colonization, the stress faced by the pathogen, and various adaptations to evade these stress conditions by forming biofilm and the morphological changes acquired by the pathogen in mature biofilm. Furthermore, we have explained the human GI tract antimicrobial peptides and the reason behind the failure of these AMPs, and how encapsulation of Pexiganan-A(MSI-78A) in a chitosan microsphere increases the efficiency of eradication.

Cooperatively enhanced reactivity and "stabilitaxis" of dissociating oligomeric proteins

Many functional units in biology, such as enzymes or molecular motors, are composed of several subunits that can reversibly assemble and disassemble. This includes oligomeric proteins composed of several smaller monomers, as well as protein complexes assembled from a few proteins. By studying the generic spatial transport properties of such proteins, we investigate here whether their ability to reversibly associate and dissociate may confer on them a functional advantage with respect to nondissociating proteins. In uniform environments with position-independent association-dissociation, we find that enhanced diffusion in the monomeric state coupled to reassociation into the functional oligomeric form leads to enhanced reactivity with localized targets. In nonuniform environments with position-dependent association-dissociation, caused by, for example, spatial gradients of an inhibiting chemical, we find that dissociating proteins generically tend to accumulate in regions where they are most stable, a process that we term "stabilitaxis."

Disordered linkers in multidomain allosteric proteins: Entropic effect to favor the open state or enhanced local concentration to favor the closed state?

There are many multidomain allosteric proteins where an allosteric signal at the allosteric domain modifies the activity of the functional domain. Intrinsically disordered regions (linkers) are widely involved in this kind of regulation process, but the essential role they play therein is not well understood. Here, we investigated the effect of linkers in stabilizing the open or the closed states of multidomain proteins using combined thermodynamic deduction and coarse-grained molecular dynamics simulations. We revealed that the influence of linker can be fully characterized by an effective local concentration [B]₀ . When K_d is smaller than [B]₀ , the closed state would be favored; while the open state would be preferred when K_d is larger than [B]₀ . We used four protein systems with markedly different domain-domain binding affinity and structural order/disorder as model systems to understand the relationship between [B]₀ and the linker length as well as its flexibility. The linker length is the main practical determinant of [B]₀ . [B]₀ of a flexible linker with 40-60 residues was determined to be in a narrow range of 0.2-0.6 mM, while a too short or too long length would dramatically decrease [B]₀ . With the revealed [B]₀ range, the introduction of a flexible linker makes the regulation of weakly interacting partners possible.

Dynamic domain arrangement of CheA-CheY complex regulates bacterial thermotaxis, as revealed by NMR

Bacteria utilize thermotaxis signal transduction proteins, including CheA, and CheY, to switch the direction of the cell movement. However, the thermally responsive machinery enabling warm-seeking behavior has not been identified. Here we examined the effects of temperature on the structure and dynamics of the full-length CheA and CheY complex, by NMR. Our studies revealed that the CheA-CheY complex exists in equilibrium between multiple states, including one state that is preferable for the autophosphorylation of CheA, and another state that is preferable for the phosphotransfer from CheA to CheY. With increasing temperature, the equilibrium shifts toward the latter state. The temperature-dependent population shift of the dynamic domain arrangement of the CheA-CheY complex induced changes in the concentrations of phosphorylated CheY that are comparable to those induced by chemical attractants or repellents. Therefore, the dynamic domain arrangement of the CheA-CheY complex functions as the primary thermally responsive machinery in warm-seeking behavior.

Interplay between binding affinity and kinetics in protein-protein interactions

To clarify the interplay between the binding affinity and kinetics of protein-protein interactions, and the possible role of intrinsically disordered proteins in such interactions, molecular simulations were carried out on 20 protein complexes. With bias potential and reweighting techniques, the free energy profiles were obtained under physiological affinities, which showed that the bound-state valley is deep with a barrier height of 12 - 33 RT. From the dependence of the affinity on interface interactions, the entropic contribution to the binding affinity is approximated to be proportional to the interface area. The extracted dissociation rates based on the Arrhenius law correlate reasonably well with the experimental values (Pearson correlation coefficient R = 0.79). For each protein complex, a linear free energy relationship between binding affinity and the dissociation rate was confirmed, but the distribution of the slopes for intrinsically disordered proteins showed no essential difference with that observed for ordered proteins. A comparison with protein folding was also performed. Proteins 2016; 84:920-933. © 2016 Wiley Periodicals, Inc.

Preformed Soluble Chemoreceptor Trimers That Mimic Cellular Assembly States and Activate CheA Autophosphorylation

Bacterial chemoreceptors associate with the histidine kinase CheA and coupling protein CheW to form extended membrane arrays that receive and transduce environmental signals. A receptor trimers-of-dimers resides at each vertex of the hexagonal protein lattice. CheA is fully activated and regulated when it is integrated into the receptor assembly. To mimic these states in solution, we have engineered chemoreceptor cytoplasmic kinase-control modules (KCMs) based on the Escherichia coli aspartate receptor Tar that are covalently fused and trimerized by a foldon domain (Tar(FO)). Small-angle X-ray scattering, multi-angle light scattering, and pulsed-dipolar electron spin resonance spectroscopy of spin-labeled proteins indicate that the Tar(FO) modules assemble into homogeneous trimers wherein the protein interaction regions closely associate at the end opposite to the foldon domains. The Tar(FO) variants greatly increase the saturation levels of phosphorylated CheA (CheA-P), indicating that the association with a trimer of receptor dimers changes the fraction of active kinase. However, the rate constants for CheA-P formation with the Tar variants are low compared to those for autophosphorylation by free CheA, and net phosphotransfer from CheA to CheY does not increase commensurately with CheA autophosphorylation. Thus, the Tar variants facilitate slow conversion to an active form of CheA that then undergoes stable autophosphorylation and is capable of subsequent phosphotransfer to CheY. Free CheA is largely incapable of phosphorylation but contains a small active fraction. Addition of Tar(FO) to CheA promotes a planar conformation of the regulatory domains consistent with array models for the assembly state of the ternary complex and different from that observed with a single inhibitory receptor. Introduction of Tar(FO) into E. coli cells activates endogenous CheA to produce increased clockwise flagellar rotation, with the effects increasing in the presence of the chemotaxis methylation system (CheB/CheR). Overall, the Tar(FO) modules demonstrate that trimerized signaling tips self-associate, bind CheA and CheW, and facilitate conversion of CheA to an active conformation.

Remarkably fast coupled folding and binding of the intrinsically disordered transactivation domain of cMyb to CBP KIX

Association rates for interactions between folded proteins have been investigated extensively, allowing the development of computational and theoretical prediction methods. Less is known about association rates for complexes where one or more partner is initially disordered, despite much speculation about how they may compare to those for folded proteins. We have attached a fluorophore to the N-terminus of the 25 amino acid cMyb peptide used previously in NMR and equilibrium studies (termed FITC-cMyb), and used this to monitor the kinetics of its interaction with the KIX protein. We have investigated the ionic strength and temperature dependence of the kinetics, and conclude that the association process is extremely fast, apparently exceeding the rates predicted by formulations applicable to interactions between pairs of folded proteins. This is despite the fact that not all collisions result in complex formation (there is an observable activation energy for the association process). We propose that this is partially a result of the disordered nature of the FITC-cMyb peptide itself.

Folding and binding of an intrinsically disordered protein: fast, but not 'diffusion-limited'

Coupled folding and binding of intrinsically disordered proteins (IDPs) is prevalent in biology. As the first step toward understanding the mechanism of binding, it is important to know if a reaction is 'diffusion-limited' as, if this speed limit is reached, the association must proceed through an induced fit mechanism. Here, we use a model system where the 'BH3 region' of PUMA, an IDP, forms a single, contiguous α-helix upon binding the folded protein Mcl-1. Using stopped-flow techniques, we systematically compare the rate constant for association (k(+)) under a number of solvent conditions and temperatures. We show that our system is not 'diffusion-limited', despite having a k(+) in the often-quoted 'diffusion-limited' regime (10(5)-10(6) M(-1) s(-1) at high ionic strength) and displaying an inverse dependence on solvent viscosity. These standard tests, developed for folded protein-protein interactions, are not appropriate for reactions where one protein is disordered.

Slow, reversible, coupled folding and binding of the spectrin tetramerization domain

Many intrinsically disordered proteins (IDPs) are significantly unstructured under physiological conditions. A number of these IDPs have been shown to undergo coupled folding and binding reactions whereby they can gain structure upon association with an appropriate partner protein. In general, these systems display weaker binding affinities than do systems with association between completely structured domains, with micromolar K(d) values appearing typical. One such system is the association between α- and β-spectrin, where two partially structured, incomplete domains associate to form a fully structured, three-helix bundle, the spectrin tetramerization domain. Here, we use this model system to demonstrate a method for fitting association and dissociation kinetic traces where, using typical biophysical concentrations, the association reactions are expected to be highly reversible. We elucidate the unusually slow, two-state kinetics of spectrin assembly in solution. The advantages of studying kinetics in this regime include the potential for gaining equilibrium constants as well as rate constants, and for performing experiments with low protein concentrations. We suggest that this approach would be particularly appropriate for high-throughput mutational analysis of two-state reversible binding processes.

Automated prediction of protein association rate constants

The association rate constants (k(a)) of proteins with other proteins or other macromolecular targets are a fundamental biophysical property. Observed rate constants span over ten orders of magnitude, from 1 to 10(10) M(-1)s(-1). Protein association can be rate limited either by the diffusional approach of the subunits to form a transient complex, with near-native separation and orientation but without short-range native interactions, or by the subsequent conformational rearrangement to form the native complex. Our transient-complex theory showed promise in predicting k(a) in the diffusion-limited regime. Here, we develop it into a web server called TransComp (http://pipe.sc.fsu.edu/transcomp/) and report on the server's accuracy and robustness based on applications to over 100 protein complexes. We expect this server to be a valuable tool for systems biology applications and for kinetic characterization of protein-protein and protein-nucleic acid association in general.

Two component systems: physiological effect of a third component

Signal transduction systems mediate the response and adaptation of organisms to environmental changes. In prokaryotes, this signal transduction is often done through Two Component Systems (TCS). These TCS are phosphotransfer protein cascades, and in their prototypical form they are composed by a kinase that senses the environmental signals (SK) and by a response regulator (RR) that regulates the cellular response. This basic motif can be modified by the addition of a third protein that interacts either with the SK or the RR in a way that could change the dynamic response of the TCS module. In this work we aim at understanding the effect of such an additional protein (which we call "third component") on the functional properties of a prototypical TCS. To do so we build mathematical models of TCS with alternative designs for their interaction with that third component. These mathematical models are analyzed in order to identify the differences in dynamic behavior inherent to each design, with respect to functionally relevant properties such as sensitivity to changes in either the parameter values or the molecular concentrations, temporal responsiveness, possibility of multiple steady states, or stochastic fluctuations in the system. The differences are then correlated to the physiological requirements that impinge on the functioning of the TCS. This analysis sheds light on both, the dynamic behavior of synthetically designed TCS, and the conditions under which natural selection might favor each of the designs. We find that a third component that modulates SK activity increases the parameter space where a bistable response of the TCS module to signals is possible, if SK is monofunctional, but decreases it when the SK is bifunctional. The presence of a third component that modulates RR activity decreases the parameter space where a bistable response of the TCS module to signals is possible.

Sinorhizobium meliloti CheA complexed with CheS exhibits enhanced binding to CheY1, resulting in accelerated CheY1 dephosphorylation

Retrophosphorylation of the histidine kinase CheA in the chemosensory transduction chain is a widespread mechanism for efficient dephosphorylation of the activated response regulator. First discovered in Sinorhizobium meliloti, the main response regulator CheY2-P shuttles its phosphoryl group back to CheA, while a second response regulator, CheY1, serves as a sink for surplus phosphoryl groups from CheA-P. We have identified a new component in this phospho-relay system, a small 97-amino-acid protein named CheS. CheS has no counterpart in enteric bacteria but revealed distinct similarities to proteins of unknown function in other members of the α subgroup of proteobacteria. Deletion of cheS causes a phenotype similar to that of a cheY1 deletion strain. Fluorescence microscopy revealed that CheS is part of the polar chemosensory cluster and that its cellular localization is dependent on the presence of CheA. In vitro binding, as well as coexpression and copurification studies, gave evidence of CheA/CheS complex formation. Using limited proteolysis coupled with mass spectrometric analyses, we defined CheA(163-256) to be the CheS binding domain, which overlaps with the N-terminal part of the CheY2 binding domain (CheA(174-316)). Phosphotransfer experiments using isolated CheA-P showed that dephosphorylation of CheY1-P but not CheY2-P is increased in the presence of CheS. As determined by surface plasmon resonance spectroscopy, CheY1 binds ∼100-fold more strongly to CheA/CheS than to CheA. We propose that CheS facilitates signal termination by enhancing the interaction of CheY1 and CheA, thereby promoting CheY1-P dephosphorylation, which results in a more efficient drainage of the phosphate sink.

Electrostatics in proteins and protein-ligand complexes

Accurate computational methods for predicting electrostatic energies are of major importance for our understanding of protein energetics in general for computer-aided drug design as well as for the design of novel biocatalysts and protein therapeutics. Electrostatic energies are of particular importance in such applications as virtual screening, drug design and protein-protein docking due to the high charge density of protein ligands and small-molecule drugs, and the frequent protonation state changes observed when drugs bind to their protein targets. Therefore, the development of a reliable and fast algorithm for the evaluation of electrostatic free energies, as an important contributor to the overall protein energy function, has been the focus for many scientists over the past three decades. In this review we describe the current state-of-the-art in modeling electrostatic effects in proteins and protein-ligand complexes. We focus mainly on the merits and drawbacks of the continuum methodology, and speculate on future directions in refining algorithms for calculating electrostatic energies in proteins using experimental data.

Identification of an anchor residue for CheA-CheY interactions in the chemotaxis system of Escherichia coli

Transfer of a phosphoryl group from autophosphorylated CheA (P-CheA) to CheY is an important step in the bacterial chemotaxis signal transduction pathway. This reaction involves CheY (i) binding to the P2 domain of P-CheA and then (ii) acquiring the phosphoryl group from the P1 domain. Crystal structures indicated numerous side chain interactions at the CheY-P2 binding interface. To investigate the individual contributions of the P2 side chains involved in these contacts, we analyzed the effects of eight alanine substitution mutations on CheA-CheY binding interactions. An F214A substitution in P2 caused ∼1,000-fold reduction in CheA-CheY binding affinity, while Ala substitutions at other P2 positions had small effects (E171A, E178A, and I216A) or no detectable effects (H181A, D202A, D207A, and C213A) on binding affinity. These results are discussed in relation to previous in silico predictions of hot-spot and anchor positions at the CheA-CheY interface. We also investigated the consequences of these mutations for chemotaxis signal transduction in living cells. CheA(F214A) was defective in mediating localization of CheY-YFP to the large clusters of signaling proteins that form at the poles of Escherichia coli cells, while the other CheA variants did not differ from wild-type (wt) CheA (CheA(wt)) in this regard. In our set of mutants, only CheA(F214A) exhibited a markedly diminished ability to support chemotaxis in motility agar assays. Surprisingly, however, in FRET assays that monitored receptor-regulated production of phospho-CheY, CheA(F214A) (and each of the other Ala substitution mutants) performed just as well as CheA(wt). Overall, our findings indicate that F214 serves as an anchor residue at the CheA-CheY interface and makes an important contribution to the binding energy in vitro and in vivo; however, loss of this contribution does not have a large negative effect on the overall ability of the signaling pathway to modulate P-CheY levels in response to chemoattractants.

Structure of the ternary complex formed by a chemotaxis receptor signaling domain, the CheA histidine kinase, and the coupling protein CheW as determined by pulsed dipolar ESR spectroscopy

The signaling apparatus that controls bacterial chemotaxis is composed of a core complex containing chemoreceptors, the histidine autokinase CheA, and the coupling protein CheW. Site-specific spin labeling and pulsed dipolar ESR spectroscopy (PDS) have been applied to investigate the structure of a soluble ternary complex formed by Thermotoga maritima CheA (TmCheA), CheW, and receptor signaling domains. Thirty-five symmetric spin-label sites (SLSs) were engineered into the five domains of the CheA dimer and CheW to provide distance restraints within the CheA:CheW complex in the absence and presence of a soluble receptor that inhibits kinase activity (Tm14). Additional PDS restraints among spin-labeled CheA, CheW, and an engineered single-chain receptor labeled at six different sites allow docking of the receptor structure relative to the CheA:CheW complex. Disulfide cross-linking between selectively incorporated Cys residues finds two pairs of positions that provide further constraints within the ternary complex: one involving Tm14 and CheW and another involving Tm14 and CheA. The derived structure of the ternary complex indicates a primary site of interaction between CheW and Tm14 that agrees well with previous biochemical and genetic data for transmembrane chemoreceptors. The PDS distance distributions are most consistent with only one CheW directly engaging one dimeric Tm14. The CheA dimerization domain (P3) aligns roughly antiparallel to the receptor-conserved signaling tip but does not interact strongly with it. The angle of the receptor axis with respect to P3 and the CheW-binding P5 domains is bound by two limits differing by approximately 20 degrees . In one limit, Tm14 aligns roughly along P3 and may interact to some extent with the hinge region near the P3 hairpin loop. In the other limit, Tm14 tilts to interact with the P5 domain of the opposite subunit in an interface that mimics that observed with the P5 homologue CheW. The time domain ESR data can be simulated from the model only if orientational variability is introduced for the P5 and, especially, P3 domains. The Tm14 tip also binds beside one of the CheA kinase domains (P4); however, in both bound and unbound states, P4 samples a broad range of distributions that are only minimally affected by Tm14 binding. The CheA P1 domains that contain the substrate histidine are also broadly distributed in space under all conditions. In the context of the hexagonal lattice formed by trimeric transmembrane chemoreceptors, the PDS structure is best accommodated with the P3 domain in the center of a honeycomb edge.

Rate theories for biologists

Some of the rate theories that are most useful for modeling biological processes are reviewed. By delving into some of the details and subtleties in the development of the theories, the review will hopefully help the reader gain a more than superficial perspective. Examples are presented to illustrate how rate theories can be used to generate insight at the microscopic level into biomolecular behaviors. An attempt is made to clear up a number of misconceptions in the literature regarding popular rate theories, including the appearance of Planck's constant in the transition-state theory and the Smoluchowski result as an upper limit for protein-protein and protein-DNA association rate constants. Future work in combining the implementation of rate theories through computer simulations with experimental probes of rate processes, and in modeling effects of intracellular environments so that theories can be used for generating rate constants for systems biology studies is particularly exciting.

Stochastic kinetic model of two component system signalling reveals all-or-none, graded and mixed mode stochastic switching responses

Two-component systems (TCSs) are prevalent signal transduction systems in bacteria that control innumerable adaptive responses to environmental cues and host-pathogen interactions. We constructed a detailed stochastic kinetic model of two component signalling based on published data. Our model has been validated with flow cytometry data and used to examine reporter gene expression in response to extracellular signal strength. The model shows that, depending on the actual kinetic parameters, TCSs exhibit all-or-none, graded or mixed mode responses. In accordance with other studies, positively autoregulated TCSs exhibit all-or-none responses. Unexpectedly, our model revealed that TCSs lacking a positive feedback loop exhibit not only graded but also mixed mode responses, in which variation of the signal strength alters the level of gene expression in induced cells while the regulated gene continues to be expressed at the basal level in a substantial fraction of cells. The graded response of the TCS changes to mixed mode response by an increase of the translation initiation rate of the histidine kinase. Thus, a TCS is an evolvable design pattern capable of implementing deterministic regulation and stochastic switches associated with both graded and threshold responses. This has implications for understanding the emergence of population diversity in pathogenic bacteria and the design of genetic circuits in synthetic biology applications. The model is available in systems biology markup language (SBML) and systems biology graphical notation (SBGN) formats and can be used as a component of large-scale biochemical reaction network models.

Prediction of salt and mutational effects on the association rate of U1A protein and U1 small nuclear RNA stem/loop II

We have developed a computational approach for predicting protein-protein association rates (Alsallaq and Zhou, Structure 2007, 15, 215). Here we expand the range of applicability of this approach to protein-RNA binding and report the first results for protein-RNA binding rates predicted from atomistic modeling. The system studied is the U1A protein and stem/loop II of the U1 small nuclear RNA. Experimentally it was observed that the binding rate is significantly reduced by increasing salt concentration while the dissociation changes little with salt concentration, and charges distant from the binding site make marginal contribution to the binding rate. These observations are rationalized. Moreover, predicted effects of salt and charge mutations are found to be in quantitative agreement with experimental results.

Distribution, structure and diversity of "bacterial" genes encoding two-component proteins in the Euryarchaeota

The publicly available annotated archaeal genome sequences (23 complete and three partial annotations, October 2005) were searched for the presence of potential two-component open reading frames (ORFs) using gene category lists and BLASTP. A total of 489 potential two-component genes were identified from the gene category lists and BLASTP. Two-component genes were found in 14 of the 21 Euryarchaeal sequences (October 2005) and in neither the Crenarchaeota nor the Nanoarchaeota. A total of 20 predicted protein domains were identified in the putative two-component ORFs that, in addition to the histidine kinase and receiver domains, also includes sensor and signalling domains. The detailed structure of these putative proteins is shown, as is the distribution of each class of two-component genes in each species. Potential members of orthologous groups have been identified, as have any potential operons containing two or more two-component genes. The number of two-component genes in those Euryarchaeal species which have them seems to be linked more to lifestyle and habitat than to genome complexity, with most examples being found in Methanospirillum hungatei, Haloarcula marismortui, Methanococcoides burtonii and the mesophilic Methanosarcinales group. The large numbers of two-component genes in these species may reflect a greater requirement for internal regulation. Phylogenetic analysis of orthologous groups of five different protein classes, three probably involved in regulating taxis, suggests that most of these ORFs have been inherited vertically from an ancestral Euryarchaeal species and point to a limited number of key horizontal gene transfer events.

Diversity in times of adversity: probabilistic strategies in microbial survival games

Population diversification strategies are ubiquitous among microbes, encompassing random phase-variation (RPV) of pathogenic bacteria, viral latency as observed in some bacteriophage and HIV, and the non-genetic diversity of bacterial stress responses. Precise conditions under which these diversification strategies confer an advantage have not been well defined. We develop a model of population growth conditioned on dynamical environmental and cellular states. Transitions among cellular states, in turn, may be biased by possibly noisy readings of the environment from cellular sensors. For various types of environmental dynamics and cellular sensor capability, we apply game-theoretic analysis to derive the evolutionarily stable strategy (ESS) for an organism and determine when that strategy is diversification. We find that: (1) RPV, effecting a sort of Parrondo paradox wherein random alternations between losing strategies produce a winning strategy, is selected when transitions between different selective environments cannot be sensed, (2) optimal RPV cell switching rates are a function of environmental lifecycle asymmetries and environmental autocorrelation, (3) probabilistic diversification upon entering a new environment is selected when sensors can detect environmental transitions but have poor precision in identifying new environments, and (4) in the presence of excess additive noise, low-pass filtering is required for evolutionary stability. We show that even when RPV is not the ESS, it may minimize growth rate variance and the risk of extinction due to 'unlucky' environmental dynamics.

Kinetic analysis of YPD1-dependent phosphotransfer reactions in the yeast osmoregulatory phosphorelay system

In Saccharomyces cerevisiae, the histidine-containing phosphotransfer (HPt) protein YPD1 transfers phosphoryl groups between the three different response regulator domains of SLN1, SSK1, and SKN7 (designated R1, R2, and R3, respectively). Together these proteins form a branched histidine-aspartic acid phosphorelay system through which cells can respond to hyperosmotic and other environmental stresses. The in vivo order of phosphotransfer reactions is believed to proceed from SLN1-R1 to YPD1 and then subsequently to SSK1-R2 or SKN7-R3. The individual phosphoryl transfer reactions between YPD1 and the response regulator domains have been examined kinetically. A maximum forward rate constant of 29 s(-)(1) was determined for the reaction between SLN1-R1 approximately P and YPD1 with a K(d) of 1.4 microM for the SLN1-R1 approximately P.YPD1 complex. In the subsequent reactions, phosphotransfer from YPD1 to SSK1-R2 is very rapid (160 s(-)(1)) and is strongly favored over phosphotransfer to SKN7-R3. Phosphotransfer reactions between YPD1 and SLN1-R1 or SKN7-R3 were reversible. In contrast, no reverse transfer from SSK1-R2 approximately P to YPD1 was observed. These findings are consistent with the notion that SSK1 is constitutively phosphorylated under normal osmotic conditions. In addition, we have examined the roles of several conserved amino acid residues surrounding the phosphorylatable histidine (H64) of YPD1 using phosphoryl transfer reactions involving YPD1 mutants. With respect to phosphoryl transfer from SLN1-R1 approximately P, only one YPD1 mutant (K67A) exhibited an increase in K(d) and thus affects binding of YPD1 to SLN1-R1 approximately P, whereas other mutants (R90A, Q86A, and G68Q) showed a decrease in phosphoryl transfer rate. Only the G68Q-YPD1 mutant was significantly affected in phosphotransfer to SSK1-R2 ( approximately 680-fold decrease in rate in comparison to wild-type). This is the first report of a kinetic analysis of a eukaryotic "two-component" histidine-aspartic acid phosphotransfer system, enabling a comparison of the transfer rates and binding constants to the few bacterial systems that have been studied this way.

Phosphorylation and Binding Interactions of CheY Studied by Use of Badan-Labeled Protein†

In the chemotaxis signal transduction pathway of Escherichia coli, the response regulator protein CheY is phosphorylated by the receptor-coupled protein kinase CheA. Previous studies of CheY phosphorylation and CheY interactions with other proteins in the chemotaxis pathway have exploited the fluorescence properties of Trp(58), located immediately adjacent to the phosphorylation site of CheY (Asp(57)). Such studies can be complicated by the intrinsic fluorescence and absorbance properties of CheA and other proteins of interest. To circumvent these difficulties, we generated a derivative of CheY carrying a covalently attached fluorescent label that serves as a sensitive reporter of phosphorylation and binding events and that absorbs and emits light at wavelengths well removed from potential interference by other proteins. This labeled version of CheY has the (dimethylamino)naphthalene fluorophore from Badan [6-bromoacetyl-2-(dimethylamino)naphthalene] attached to the thiol group of a cysteine introduced at position 17 of CheY by site-directed mutagenesis. Under phosphorylating conditions (or in the presence of beryllofluoride), the fluorescence emission of Badan-labeled CheY(M17C) exhibited an approximately 10 nm blue shift and an approximately 30% increase in signal intensity at 490 nm. The fluorescence of Badan-labeled CheY(M17C) also served as a sensitive reporter of CheY-CheA binding interactions, exhibiting an approximately 50% increase in emission intensity in the presence of saturating levels of CheA. Compared to wild-type CheY, Badan-labeled CheY exhibited reduced ability to autodephosphorylate and could not interact productively with the phosphatase CheZ. However, with respect to autophosphorylation and interactions with CheA, Badan-CheY performed identically to wild-type CheY, allowing us to explore CheA-CheY phosphotransfer kinetics and binding kinetics without interference from the fluorescence/absorbance properties of CheA and ATP. These results provide insights into CheY interactions with CheA, CheZ, and other components of the chemotaxis signaling pathway.