Further insights into the mechanism of function of the response regulator CheY from crystallographic studies of the CheY--CheA(124--257) complex

A Coarse-Grained Methodology Identifies Intrinsic Mechanisms That Dissociate Interacting Protein Pairs

We address the problem of triggering dissociation events between proteins that have formed a complex. We have collected a set of 25 non-redundant, functionally diverse protein complexes having high-resolution three-dimensional structures in both the unbound and bound forms. We unify elastic network models with perturbation response scanning (PRS) methodology as an efficient approach for predicting residues that have the propensity to trigger dissociation of an interacting protein pair, using the three-dimensional structures of the bound and unbound proteins as input. PRS reveals that while for a group of protein pairs, residues involved in the conformational shifts are confined to regions with large motions, there are others where they originate from parts of the protein unaffected structurally by binding. Strikingly, only a few of the complexes have interface residues responsible for dissociation. We find two main modes of response: In one mode, remote control of dissociation in which disruption of the electrostatic potential distribution along protein surfaces play the major role; in the alternative mode, mechanical control of dissociation by remote residues prevail. In the former, dissociation is triggered by changes in the local environment of the protein, e.g., pH or ionic strength, while in the latter, specific perturbations arriving at the controlling residues, e.g., via binding to a third interacting partner is required for decomplexation. We resolve the observations by relying on an electromechanical coupling model which reduces to the usual elastic network result in the limit of the lack of coupling. We validate the approach by illustrating the biological significance of top residues selected by PRS on select cases where we show that the residues whose perturbation leads to the observed conformational changes correspond to either functionally important or highly conserved residues in the complex.

Properties that rank protein:protein docking poses with high accuracy

The development of docking algorithms to predict near-native structures of protein:protein complexes from the structure of the isolated monomers is of paramount importance for molecular biology and drug discovery.

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.

A Machine Learning Approach for Hot-Spot Detection at Protein-Protein Interfaces

Understanding protein-protein interactions is a key challenge in biochemistry. In this work, we describe a more accurate methodology to predict Hot-Spots (HS) in protein-protein interfaces from their native complex structure compared to previous published Machine Learning (ML) techniques. Our model is trained on a large number of complexes and on a significantly larger number of different structural- and evolutionary sequence-based features. In particular, we added interface size, type of interaction between residues at the interface of the complex, number of different types of residues at the interface and the Position-Specific Scoring Matrix (PSSM), for a total of 79 features. We used twenty-seven algorithms from a simple linear-based function to support-vector machine models with different cost functions. The best model was achieved by the use of the conditional inference random forest (c-forest) algorithm with a dataset pre-processed by the normalization of features and with up-sampling of the minor class. The method has an overall accuracy of 0.80, an F1-score of 0.73, a sensitivity of 0.76 and a specificity of 0.82 for the independent test set.

Solution structure of a complex of the histidine autokinase CheA with its substrate CheY

In the bacterial chemotaxis two-component signaling system, the histidine-containing phosphotransfer domain (the "P1" domain) of CheA receives a phosphoryl group from the catalytic domain (P4) of CheA and transfers it to the cognate response regulator (RR) CheY, which is docked by the P2 domain of CheA. Phosphorylated CheY then diffuses into the cytoplasm and interacts with the FliM moiety of the flagellar motors, thereby modulating the direction of flagellar rotation. Structures of various histidine phosphotransfer domains (HPt) complexed with their cognate RR domains have been reported. Unlike the Escherichia coli chemotaxis system, however, these systems lack the additional domains dedicated to binding to the response regulators, and the interaction of an HPt domain with an RR domain in the presence of such a domain has not been examined on a structural basis. In this study, we used modern nuclear magnetic resonance techniques to construct a model for the interaction of the E. coli CheA P1 domain (HPt) and CheY (RR) in the presence of the CheY-binding domain, P2. Our results indicate that the presence of P2 may lead to a slightly different relative orientation of the HPt and RR domains versus those seen in such complex structures previously reported.

A structural model of the E. coli PhoB dimer in the transcription initiation complex

BACKGROUND

There exist > 78,000 proteins and/or nucleic acids structures that were determined experimentally. Only a small portion of these structures corresponds to those of protein complexes. While homology modeling is able to exploit knowledge-based potentials of side-chain rotomers and backbone motifs to infer structures for new proteins, no such general method exists to extend our understanding of protein interaction motifs to novel protein complexes.

RESULTS

We use a Motif Binding Geometries (MBG) approach, to infer the structure of a protein complex from the database of complexes of homologous proteins taken from other contexts (such as the helix-turn-helix motif binding double stranded DNA), and demonstrate its utility on one of the more important regulatory complexes in biology, that of the RNA polymerase initiating transcription under conditions of phosphate starvation. The modeled PhoB/RNAP/σ-factor/DNA complex is stereo-chemically reasonable, has sufficient interfacial Solvent Excluded Surface Areas (SESAs) to provide adequate binding strength, is physically meaningful for transcription regulation, and is consistent with a variety of known experimental constraints.

CONCLUSIONS

Based on a straightforward and easy to comprehend concept, "proteins and protein domains that fold similarly could interact similarly", a structural model of the PhoB dimer in the transcription initiation complex has been developed. This approach could be extended to enable structural modeling and prediction of other bio-molecular complexes. Just as models of individual proteins provide insight into molecular recognition, catalytic mechanism, and substrate specificity, models of protein complexes will provide understanding into the combinatorial rules of cellular regulation and signaling.

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.

A tale of two ferredoxins: sequence similarity and structural differences

Background

Sequence similarity between proteins is usually considered a reliable indicator of homology. Pyruvate-ferredoxin oxidoreductase and quinol-fumarate reductase contain ferredoxin domains that bind [Fe-S] clusters and are involved in electron transport. Profile-based methods for sequence comparison, such as PSI-BLAST and HMMer, suggest statistically significant similarity between these domains.

Results

The sequence similarity between these ferredoxin domains resides in the area of the [Fe-S] cluster-binding sites. Although overall folds of these ferredoxins bear no obvious similarity, the regions of sequence similarity display a remarkable local structural similarity. These short regions with pronounced sequence motifs are incorporated in completely different structural environments. In pyruvate-ferredoxin oxidoreductase (bacterial ferredoxin), the hydrophobic core of the domain is completed by two β-hairpins, whereas in quinol-fumarate reductase (α-helical ferredoxin), the cluster-binding motifs are part of a larger all-α-helical globin-like fold core.

Conclusion

Functionally meaningful sequence similarity may sometimes be reflected only in local structural similarity, but not in global fold similarity. If detected and used naively, such similarities may lead to incorrect fold predictions.

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.

Searching for potential drug targets in two-component and phosphorelay signal-transduction systems using three-dimensional cluster analysis

Two-component and phosphorelay signal transduction systems are central components in the virulence and antimicrobial resistance responses of a number of bacterial and fungal pathogens; in some cases, these systems are essential for bacterial growth and viability. Herein, we analyze in detail the conserved surface residue clusters in the phosphotransferase domain of histidine kinases and the regulatory domain of response regulators by using complex structure-based three-dimensional cluster analysis. We also investigate the protein-protein interactions that these residue clusters participate in. The Spo0B-Spo0F complex structure was used as the reference structure, and the multiple aligned sequences of phosphotransferases and response regulators were paired correspondingly. The results show that a contiguous conserved residue cluster is formed around the active site, which crosses the interface of histidine kinases and response regulators. The conserved residue clusters of phosphotransferase and the regulatory domains are directly involved in the functional implementation of two-component signal transduction systems and are good targets for the development of novel antimicrobial agents.

In different organisms, the mode of interaction between two signaling proteins is not necessarily conserved

Although interfaces mediating protein-protein interactions are thought to be under strong evolutionary constraints, binding of the chemotaxis histidine kinase CheA to its phosphorylation target CheY suggests otherwise. The structure of Thermotoga maritima CheA domain P2 in complex with CheY reveals a different association than that observed for the same Escherichia coli proteins. Similar regions of CheY bind CheA P2 in the two systems, but the CheA P2 domains differ by an approximately 90 degrees rotation. CheA binds CheY with identical affinity in T. maritima and E. coli at the vastly different temperatures where the respective organisms live. Distinct sets of P2 residues mediate CheY binding in the two complexes; conservation patterns of these residues in CheA and compensations in CheY delineate two families of prokaryotic chemotaxis systems. A protein complex that has the same components and general function in different organisms, but an altered structure, indicates unanticipated complexity in the evolution of protein-protein interactions and cautions against extrapolating structural data from homologs.

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

The chemotaxis system of Escherichia coli makes use of an extended two-component sensory response pathway in which CheA, an autophosphorylating protein histidine kinase (PHK) rapidly passes its phosphoryl group to CheY, a phospho-accepting response regulator protein (RR). The CheA-->CheY phospho-transfer reaction is 100-1000 times faster than the His-->Asp phospho-relays that operate in other (non-chemotaxis) two-component regulatory systems, suggesting that CheA and CheY have unique features that enhance His-->Asp phospho-transfer kinetics. One such feature could be the P2 domain of CheA. P2 encompasses a binding site for CheY, but an analogous RR-binding domain is not found in other PHKs. In previous work, we removed P2 from CheA, and this decreased the catalytic efficiency of CheA-->CheY phospho-transfer by a factor of 50-100. Here we examined the kinetics of the binding interactions between CheY and P2. The rapid association reaction (k(assn) approximately 10(8)M(-1)s(-1) at 25 degrees C and micro=0.03 M) exhibited a simple first-order dependence on P2 concentration and appeared to be largely diffusion-limited. Ionic strength (micro) had a moderate effect on k(assn) in a manner predictable based on the calculated electrostatic interaction energy of the protein binding surfaces and the expected Debye-Hückel shielding. The speed of binding reflects, in part, electrostatic interactions, but there is also an important contribution from the inherent plasticity of the complex and the resulting flexibility that this allows during the process of complex formation. Our results support the idea that the P2 domain of CheA contributes to the overall speed of phospho-transfer by promoting rapid association between CheY and CheA. However, this alone does not account for the ability of the chemotaxis system to operate much more rapidly than other two-component systems: k(cat) differences indicate that CheA and CheY also achieve the chemical events of phospho-transfer more rapidly than do PHK-RR pairs of slower systems.

Nuclear magnetic resonance spectroscopy reveals the functional state of the signalling protein CheY in vivo in Escherichia coli

Two-component signal transduction (TCST) pathways are regulatory systems that are highly homologous throughout the bacterial kingdom. Their established role in virulence and absence in vertebrates has made TCST an attractive target for therapeutic intervention. However, such systems have yet to yield success in the development of novel antibiotics. CheY serves as a prototype for the analysis of response regulator function. The protein structure exhibits several conformations by both X-ray and nuclear magnetic resonance (NMR) analyses. Knowledge of which structures are relevant in vivo would be valuable in a rational drug design project. Our aim was to probe the in vivo conformation and ligand binding of CheY in Escherichia coli under resting conditions by in-cell NMR methods. CheY was selectively labelled with 15N by the control of growth and expression conditions. NMR spectra obtained in vivo demonstrated that the Mg2+ complex was the predominant form even though cells were resuspended in metal-free buffers and the intracellular free Mg2+ was low. In-cell NMR also confirmed the uptake and in vivo binding mode to CheY of small-molecular-weight compounds identified in vitro. This paper reports the first observation of the structure and interactions with a potential drug of a regulator protein in its native host in vivo using NMR spectroscopy.

Different evolutionary constraints on chemotaxis proteins CheW and CheY revealed by heterologous expression studies and protein sequence analysis

CheW and CheY are single-domain proteins from a signal transduction pathway that transmits information from transmembrane receptors to flagellar motors in bacterial chemotaxis. In various bacterial and archaeal species, the cheW and cheY genes are usually encoded within homologous chemotaxis operons. We examined evolutionary changes in these two proteins from distantly related proteobacterial species, Escherichia coli and Azospirillum brasilense. We analyzed the functions of divergent CheW and CheY proteins from A. brasilense by heterologous expression in E. coli wild-type and mutant strains. Both proteins were able to specifically inhibit chemotaxis of a wild-type E. coli strain; however, only CheW from A. brasilense was able to restore signal transduction in a corresponding mutant of E. coli. Detailed protein sequence analysis of CheW and CheY homologs from the two species revealed substantial differences in the types of amino acid substitutions in the two proteins. Multiple, but conservative, substitutions were found in CheW homologs. No severe mismatches were found between the CheW homologs in positions that are known to be structurally or functionally important. Substitutions in CheY homologs were found to be less conservative and occurred in positions that are critical for interactions with other components of the signal transduction pathway. Our findings suggest that proteins from the same cellular pathway encoded by genes from the same operon have different evolutionary constraints on their structures that reflect differences in their functions.

Genetic analysis of response regulator activation in bacterial chemotaxis suggests an intermolecular mechanism

Response regulator proteins of two-component systems are usually activated by phosphorylation. The phosphorylated response regulator protein CheY-P mediates the chemotaxis response in Escherichia coli. We performed random mutagenesis and selected CheY mutants that are constitutively active in the absence of phosphorylation. Although a single amino acid substitution can lead to constitutive activation, no single DNA base change can effect such a transition. Numerous different sets of mutations that activate in synergy were selected in several different combinations. These mutations were all located on the side of CheY defined by alpha4, beta5, alpha5, and alpha1. Our findings argue against the two-state hypothesis for response regulator activation. We propose an alternative intermolecular mechanism that involves a dynamic interplay between response regulators and their effector targets.

A distinct meta-active conformation in the 1.1-A resolution structure of wild-type ApoCheY

CheY is the best characterized member of the response regulator superfamily, and as such it has become the principal model for understanding the initial molecular mechanisms of signaling in two-component systems. Normal signaling by response regulators requires phosphorylation, in combination with an activation mechanism whose conformational effects are not completely understood. CheY activation involves three events, phosphorylation, a conformational change in the beta(4)--alpha(4) loop, and a rotational restriction of the side chain of tyrosine 106. An outstanding question concerns the nature of an active conformation in the apoCheY population. The details of this 1.08-A resolution crystal structure of wild-type apoCheY shows the beta(4)--alpha(4) loop in two distinctly different conformations that sterically correlate with the two rotameric positions of the tyrosine 106 side chain. One of these conformational states of CheY is the inactive form, and we propose that the other is a meta-active form, responsible for the active properties seen in apoCheY.