基于弹性网络模型的蛋白质变构路径与关键残基识别研究

目的 变构效应在蛋白质生物学功能执行过程中发挥着重要的调控作用，如何基于蛋白质空间结构，有效识别变构信号的传播路径和关键的残基位点是蛋白质结构-功能关系研究领域的热点科学问题。方法 本研究利用基于弹性网络模型（elastic network model，ENM）的力分布计算方法，通过分析蛋白质对外力的响应过程，来识别体系的变构路径以及变构过程中的关键残基。在该方法中，对蛋白质的关键变构位点施加外力，通过对体系形变以及内力分布情况的分析，有效识别与外力承载区域形变相耦合的关键残基，从而得到力信号在蛋白质结构内的传播路径。结果 利用该方法研究了人类磷酸甘油酸激酶（human phosphoglycerate kinase，hPGK）和蛋白质酪氨酸磷酸酶（protein tyrosine phosphatase，PTP）PDZ2结构域的变构调控路径和关键残基。对于hPGK，识别出从底物结合位点到铰链区的两条变构信号传导路径。对于PTP PDZ2，也成功识别出从配体结合位点传递到蛋白质远端的两条长程变构调控路径。计算结果与实验和分子动力学（molecular dynamics，MD）模拟得到的结果一致。结论 本研究为蛋白质体系关键残基识别及变构路径研究提供了有效的分析方法。     

Chemosensory Regulation of a HEAT-Repeat Protein Couples Aggregation and Sporulation in Myxococcus xanthus

Chemosensory systems are complex, highly modified two-component systems (TCS) used by bacteria to control various biological functions ranging from motility to sporulation. Chemosensory systems and TCS both modulate phosphorelays comprised of histidine kinases and response regulators, some of which are single-domain response regulators (SD-RRs) such as CheY. In this study, we have identified and characterized the Che7 chemosensory system of Myxococcus xanthus, a common soil bacterium which displays multicellular development in response to stress. Both genetic and biochemical analyses indicate that the Che7 system regulates development via a direct interaction between the SD-RR CheY7 and a HEAT repeat domain-containing protein, Cpc7. Phosphorylation of the SD-RR affects the interaction with its target, and residues within the α4-β5-α5 fold of the REC domain govern this interaction. The identification of the Cpc7 interaction with CheY7 extends the diversity of known targets for SD-RRs in biological systems.     

Matching Biochemical Reaction Kinetics to the Timescales of Life: Structural Determinants That Influence the Autodephosphorylation Rate of Response Regulator Proteins

In two-component regulatory systems, covalent phosphorylation typically activates the response regulator signaling protein, and hydrolysis of the phosphoryl group reestablishes the inactive state. Despite highly conserved three-dimensional structures and active-site features, the rates of catalytic autodephosphorylation for different response regulators vary by a factor of almost 10⁶. Previous studies identified two variable active-site residues, corresponding to Escherichia coli CheY residues 59 and 89, that modulate response regulator autodephosphorylation rates about 100-fold. Here, a set of five CheY mutants, which match other “model” response regulators (ArcA, CusR, DctD, FixJ, PhoB, or Spo0F) at variable active-site positions corresponding to CheY residues 14, 59, and 89, were characterized functionally and structurally in an attempt to identify mechanisms that modulate autodephosphorylation rate. As expected, the autodephosphorylation rates of the CheY mutants were reduced 6- to 40-fold relative to wild-type CheY, but all still autodephosphorylated 12- to 80-fold faster than their respective model response regulators. Comparison of X-ray crystal structures of the five CheY mutants (complexed with the phosphoryl group analogue BeF₃⁻) to wild-type CheY or corresponding model response regulator structures gave strong evidence for steric obstruction of the phosphoryl group from the attacking water molecule as one mechanism to enhance phosphoryl group stability. Structural data also suggested that impeding the change of a response regulator from the active to the inactive conformation might retard the autodephosphorylation reaction if the two processes are coupled, and that the residue at position ‘58’ may contribute to rate modulation. A given combination of amino acids at positions ‘14’, ‘59’, and ‘89’ adopted similar conformations regardless of protein context (CheY or model response regulator), suggesting that knowledge of residue identity may be sufficient to predict autodephosphorylation rate, and hence the kinetics of the signaling response, in the response regulator family of proteins.     

CheY's acetylation sites responsible for generating clockwise flagellar rotation in Escherichia coli

Stimulation of Escherichia coli with acetate elevates the acetylation level of the chemotaxis response regulator CheY. This elevation, in an unknown mechanism, activates CheY to generate clockwise rotation. Here, using quantitative selective reaction monitoring mass spectrometry and high‐resolution targeted mass spectrometry, we identified K91 and K109 as the major sites whose acetylation level in vivo increases in response to acetate. Employing single and multiple lysine replacements in CheY, we found that K91 and K109 are also the sites mainly responsible for acetate‐dependent clockwise generation. Furthermore, we showed that clockwise rotation is repressed when residue K91 is nonmodified, as evidenced by an increased ability of CheY to generate clockwise rotation when K91 was acetylated or replaced by specific amino acids. Using molecular dynamics simulations, we show that K91 repression is manifested in the conformational dynamics of the β4α4 loop, shifted toward an active state upon mutation. Removal of β4α4 loop repression may represent a general activation mechanism in CheY, pertaining also to the canonical phosphorylation activation pathway as suggested by crystal structures of active and inactive CheY from Thermotoga maritima. By way of elimination, we further suggest that K109 acetylation is actively involved in generating clockwise rotation.     

Determinants of Endogenous Ligand Specificity Divergence among Metabotropic Glutamate Receptors

To determine the structural origins of diverse ligand response specificities among metabotropic glutamate receptors (mGluRs), we combined computational approaches with mutagenesis and ligand response assays to identify specificity-determining residues in the group I receptor, mGluR1, and the group III receptors, mGluR4 and mGluR7. Among these, mGluR1 responds to l-glutamate effectively, whereas it binds weakly to another endogenous ligand, l-serine-O-phosphate (l-SOP), which antagonizes the effects of l-glutamate. In contrast, mGluR4 has in common with other group III mGluR that it is activated with higher potency and efficacy by l-SOP. mGluR7 differs from mGluR4 and other group III mGluR in that l-glutamate and l-SOP activate it with low potency and efficacy. Enhanced versions of the evolutionary trace (ET) algorithm were used to identify residues that when swapped between mGluR1 and mGluR4 increased the potency of l-SOP inhibition relative to the potency of l-glutamate activation in mGluR1 mutants and others that diminished the potency/efficacy of l-SOP for mGluR4 mutants. In addition, combining ET identified swaps from mGluR4 with one identified by computational docking produced mGluR7 mutants that respond with dramatically enhanced potency/efficacy to l-SOP. These results reveal that an early functional divergence between group I/II and group III involved variation at positions primarily at allosteric sites located outside of binding pockets, whereas a later divergence within group III occurred through sequence variation both at the ligand-binding pocket and at loops near the dimerization interface and interlobe hinge region. They also demonstrate the power of ET for identifying allosteric determinants of evolutionary importance.     

The switching mechanism of the bacterial rotary motor combines tight regulation with inherent flexibility

Regulatory switches are wide spread in many biological systems. Uniquely among them, the switch of the bacterial flagellar motor is not an on/off switch but rather controls the motor’s direction of rotation in response to binding of the signaling protein CheY. Despite its extensive study, the molecular mechanism underlying this switch has remained largely unclear. Here, we resolved the functions of each of the three CheY‐binding sites at the switch in E. coli, as well as their different dependencies on phosphorylation and acetylation of CheY. Based on this, we propose that CheY motor switching activity is potentiated upon binding to the first site. Binding of potentiated CheY to the second site produces unstable switching and at the same time enables CheY binding to the third site, an event that stabilizes the switched state. Thereby, this mechanism exemplifies a unique combination of tight motor regulation with inherent switching flexibility.     

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.     

Segmental Motions, Not a Two-State Concerted Switch, Underlie Allostery in CheY

The switch between an inactive and active conformation is an important transition for signaling proteins, yet the mechanisms underlying such switches are not clearly understood. Escherichia coli CheY, a response regulator protein from the two-component signal transduction system that regulates bacterial chemotaxis, is an ideal protein for the study of allosteric mechanisms. By using (15)N CPMG relaxation dispersion experiments, we monitored the inherent dynamic switching of unphosphorylated CheY. We show that CheY does not undergo a two-state concerted switch between the inactive and active conformations. Interestingly, partial saturation of Mg(2+) enhances the intrinsic allosteric motions. Taken together with chemical shift perturbations, these data indicate that the μs-ms timescale motions underlying CheY allostery are segmental in nature. We propose an expanded allosteric network of residues, including W58, that undergo asynchronous, local switching between inactive and active-like conformations as the primary basis for the allosteric mechanism.     

CheZ Has No Effect on Flagellar Motors Activated by CheY13DK106YW

The behaviors of both cheZ-deleted and wild-type cells of Escherichia coli were found to be very sensitive to the level of expression of CheZ, a protein known to accelerate the dephosphorylation of the response regulator CheY-phosphate (CheY-P). However, cells induced to run and tumble by the unphosphorylated mutant protein CheY^13DK106YW (CheY**) failed to respond to CheZ, even when CheZ was expressed at high levels. Therefore, CheZ neither affects the flagellar motors directly nor sequesters CheY**. In in vitro cross-linking studies, CheY** promoted trimerization of CheZ to the same extent as wild-type CheY but failed to induce the formation of complexes of higher molecular weight observed with CheY-P. Also, CheY** could be cross-linked to FliM, the motor receptor protein, nearly as well as CheY-P. Thus, to CheZ, CheY** looks like CheY, but to FliM, it looks like CheY-P.     

Understanding the functional difference between growth arrest-specific protein 6 and protein S: an evolutionary approach

Although protein S (PROS1) and growth arrest-specific protein 6 (GAS6) proteins are homologous with a high degree of structural similarity, they are functionally different. The objectives of this study were to identify the evolutionary origins from which these functional differences arose. Bioinformatics methods were used to estimate the evolutionary divergence time and to detect the amino acid residues under functional divergence between GAS6 and PROS1. The properties of these residues were analysed in the light of their three-dimensional structures, such as their stability effects, the identification of electrostatic patches and the identification potential protein–protein interaction. The divergence between GAS6 and PROS1 probably occurred during the whole-genome duplications in vertebrates. A total of 78 amino acid sites were identified to be under functional divergence. One of these sites, Asn463, is involved in N-glycosylation in GAS6, but is mutated in PROS1, preventing this post-translational modification. Sites experiencing functional divergence tend to express a greater diversity of stabilizing/destabilizing effects than sites that do not experience such functional divergence. Three electrostatic patches in the LG1/LG2 domains were found to differ between GAS6 and PROS1. Finally, a surface responsible for protein–protein interactions was identified. These results may help researchers to analyse disease-causing mutations in the light of evolutionary and structural constraints, and link genetic pathology to clinical phenotypes.     

Cooperative Changes in Solvent Exposure Identify Cryptic Pockets,Switches, and Allosteric Coupling

Proteins are dynamic molecules that undergo conformational changes to a broad spectrum of different excited states. Unfortunately, the small populations of these states make it difficult to determine their structures or functional implications. Computer simulations are an increasingly powerful means to identify and characterize functionally relevant excited states. However, this advance has uncovered a further challenge: it can be extremely difficult to identify the most salient features of large simulation data sets. We reasoned that many functionally relevant conformational changes are likely to involve large, cooperative changes to the surfaces that are available to interact with potential binding partners. To examine this hypothesis, we introduce a method that returns a prioritized list of potentially functional conformational changes by segmenting protein structures into clusters of residues that undergo cooperative changes in their solvent exposure, along with the hierarchy of interactions between these groups. We term these groups exposons to distinguish them from other types of clusters that arise in this analysis and others. We demonstrate, using three different model systems, that this method identifies experimentally validated and functionally relevant conformational changes, including conformational switches, allosteric coupling, and cryptic pockets. Our results suggest that key functional sites are hubs in the network of exposons. As a further test of the predictive power of this approach, we apply it to discover cryptic allosteric sites in two different β-lactamase enzymes that are widespread sources of antibiotic resistance. Experimental tests confirm our predictions for both systems. Importantly, we provide the first evidence, to our knowledge, for a cryptic allosteric site in CTX-M-9 β-lactamase. Experimentally testing this prediction did not require any mutations and revealed that this site exerts the most potent allosteric control over activity of any pockets found in β-lactamases to date. Discovery of a similar pocket that was previously overlooked in the well-studied TEM-1 β-lactamase demonstrates the utility of exposons.     

Chemotactic control of the two flagellar systems of Rhodobacter sphaeroides is mediated by different sets of CheY and FliM proteins

Rhodobacter sphaeroides expresses two different flagellar systems, a subpolar flagellum (fla1) and multiple polar flagella (fla2). These structures are encoded by different sets of flagellar genes. The chemotactic control of the subpolar flagellum (fla1) is mediated by three of the six different CheY proteins (CheY6, CheY4, or CheY3). We show evidence that CheY1, CheY2, and CheY5 control the chemotactic behavior mediated by fla2 flagella and that RSP6099 encodes the fla2 FliM protein.     

Vitis Phylogenomics: Hybridization Intensities from a SNP Array Outperform Genotype Calls

Understanding relationships among species is a fundamental goal of evolutionary biology. Single nucleotide polymorphisms (SNPs) identified through next generation sequencing and related technologies enable phylogeny reconstruction by providing unprecedented numbers of characters for analysis. One approach to SNP-based phylogeny reconstruction is to identify SNPs in a subset of individuals, and then to compile SNPs on an array that can be used to genotype additional samples at hundreds or thousands of sites simultaneously. Although powerful and efficient, this method is subject to ascertainment bias because applying variation discovered in a representative subset to a larger sample favors identification of SNPs with high minor allele frequencies and introduces bias against rare alleles. Here, we demonstrate that the use of hybridization intensity data, rather than genotype calls, reduces the effects of ascertainment bias. Whereas traditional SNP calls assess known variants based on diversity housed in the discovery panel, hybridization intensity data survey variation in the broader sample pool, regardless of whether those variants are present in the initial SNP discovery process. We apply SNP genotype and hybridization intensity data derived from the Vitis9kSNP array developed for grape to show the effects of ascertainment bias and to reconstruct evolutionary relationships among Vitis species. We demonstrate that phylogenies constructed using hybridization intensities suffer less from the distorting effects of ascertainment bias, and are thus more accurate than phylogenies based on genotype calls. Moreover, we reconstruct the phylogeny of the genus Vitis using hybridization data, show that North American subgenus Vitis species are monophyletic, and resolve several previously poorly known relationships among North American species. This study builds on earlier work that applied the Vitis9kSNP array to evolutionary questions within Vitis vinifera and has general implications for addressing ascertainment bias in array-enabled phylogeny reconstruction.     

Structural insight into the low affinity between Thermotoga maritima CheA and CheB compared to their Escherichia coli/Salmonella typhimurium counterparts

CheA-mediated CheB phosphorylation and the subsequent CheB-mediated demethylation of the chemoreceptors are important steps required for the bacterial chemotactic adaptation response. Although Escherichia coli CheB has been reported to interact with CheA competitively against CheY, we have observed that Thermotoga maritima CheB has no detectable CheA-binding. By determining the CheY-like domain crystal structure of T. maritima CheB, and comparing against the T. maritima CheY and Salmonella typhimurium CheB structures, we propose that the two consecutive glutamates in the β4/α4 loop of T. maritima CheB that is absent in T. maritima CheY and in E. coli/S. typhimurium CheB may be one factor contributing to the low CheA affinity.     

Comparative genomic analysis of the <Emphasis Type="Italic">Lipase3</Emphasis> gene family in five plant species reveals distinct evolutionary origins

Lipases are physiologically important and ubiquitous enzymes that share a conserved domain and are classified into eight different families based on their amino acid sequences and fundamental biological properties. The Lipase3 family of lipases was reported to possess a canonical fold typical of α/β hydrolases and a typical catalytic triad, suggesting a distinct evolutionary origin for this family. Genes in the Lipase3 family do not have the same functions, but maintain the conserved Lipase3 domain. There have been extensive studies of Lipase3 structures and functions, but little is known about their evolutionary histories. In this study, all lipases within five plant species were identified, and their phylogenetic relationships and genetic properties were analyzed and used to group them into distinct evolutionary families. Each identified lipase family contained at least one dicot and monocot Lipase3 protein, indicating that the gene family was established before the split of dicots and monocots. Similar intron/exon numbers and predicted protein sequence lengths were found within individual groups. Twenty-four tandem Lipase3 gene duplications were identified, implying that the distinctive function of Lipase3 genes appears to be a consequence of translocation and neofunctionalization after gene duplication. The functional genes EDS1, PAD4, and SAG101 that are reportedly involved in pathogen response were all located in the same group. The nucleotide diversity (Dxy) and the ratio of nonsynonymous to synonymous nucleotide substitutions rates (Ka/Ks) of the three genes were significantly greater than the average across the genomes. We further observed evidence for selection maintaining diversity on three genes in the Toll-Interleukin-1 receptor type of nucleotide binding/leucine-rich repeat immune receptor (TIR-NBS LRR) immunity-response signaling pathway, indicating that they could be vulnerable to pathogen effectors.     

Allosteric Coupling in the Bacterial Adhesive Protein FimH

The protein FimH is expressed by the majority of commensal and uropathogenic strains of Escherichia coli on the tips of type 1 fimbriae and mediates adhesion via a catch bond to its ligand mannose. Crystal structures of FimH show an allosteric conformational change, but it remains unclear whether all of the observed structural differences are part of the allosteric mechanism. Here we use the protein structural analysis tool RosettaDesign combined with human insight to identify and synthesize 10 mutations in four regions that we predicted would stabilize one of the conformations of that region. The function of each variant was characterized by measuring binding to the ligand mannose, whereas the allosteric state was determined using a conformation-specific monoclonal antibody. These studies demonstrated that each region investigated was indeed part of the FimH allosteric mechanism. However, the studies strongly suggested that some regions were more tightly coupled to mannose binding and others to antibody binding. In addition, we identified many FimH variants that appear locked in the low affinity state. Knowledge of regulatory sites outside the active and effector sites as well as the ability to make FimH variants locked in the low affinity state may be crucial to the future development of novel antiadhesive and antimicrobial therapies using allosteric regulation to inhibit FimH.     

Multiple CheY Homologs Control Swimming Reversals and Transient Pauses in Azospirillum brasilense

Chemotaxis, together with motility, helps bacteria foraging in their habitat. Motile bacteria exhibit a variety of motility patterns, often controlled by chemotaxis, to promote dispersal. Motility in many bacteria is powered by a bidirectional flagellar motor. The flagellar motor has been known to briefly pause during rotation because of incomplete reversals or stator detachment. Transient pauses were previously observed in bacterial strains lacking CheY, and these events could not be explained by incomplete motor reversals or stator detachment. Here, we systematically analyzed swimming trajectories of various chemotaxis mutants of the monotrichous soil bacterium, Azospirillum brasilense. Like other polar flagellated bacterium, the main swimming pattern in A. brasilense is run and reverse. A. brasilense also uses run-pauses and putative run-reverse-flick-like swimming patterns, although these are rare events. A. brasilense mutant derivatives lacking the chemotaxis master histidine kinase, CheA4, or the central response regulator, CheY7, also showed transient pauses. Strikingly, the frequency of transient pauses increased dramatically in the absence of CheY4. Our findings collectively suggest that reversals and pauses are controlled through signaling by distinct CheY homologs, and thus are likely to be functionally important in the lifestyle of this soil organism.     

Proposed Signal Transduction Role for Conserved CheY Residue Thr87, a Member of the Response Regulator Active-Site Quintet

CheY serves as a structural prototype for the response regulator proteins of two-component regulatory systems. Functional roles have previously been defined for four of the five highly conserved residues that form the response regulator active site, the exception being the hydroxy amino acid which corresponds to Thr87 in CheY. To investigate the contribution of Thr87 to signaling, we characterized, genetically and biochemically, several cheY mutants with amino acid substitutions at this position. The hydroxyl group appears to be necessary for effective chemotaxis, as a Thr→Ser substitution was the only one of six tested which retained a Che⁺ swarm phenotype. Although nonchemotactic, cheY mutants with amino acid substitutions T87A and T87C could generate clockwise flagellar rotation either in the absence of CheZ, a protein that stimulates dephosphorylation of CheY, or when paired with a second site-activating mutation, Asp13→Lys, demonstrating that a hydroxy amino acid at position 87 is not essential for activation of the flagellar switch. All purified mutant proteins examined phosphorylated efficiently from the CheA kinase in vitro but were impaired in autodephosphorylation. Thus, the mutant CheY proteins are phosphorylated to a greater degree than wild-type CheY yet support less clockwise flagellar rotation. The data imply that Thr87 is important for generating and/or stabilizing the phosphorylation-induced conformational change in CheY. Furthermore, the various position 87 substitutions differentially affected several properties of the mutant proteins. The chemotaxis and autodephosphorylation defects were tightly linked, suggesting common structural elements, whereas the effects on self-catalyzed and CheZ-mediated dephosphorylation of CheY were uncorrelated, suggesting different structural requirements for the two dephosphorylation reactions.