Crystal structure and functional analysis of the SurE protein identify a novel phosphatase family

Predicting pKa values with continuous constant pH molecular dynamics

Knowledge of pK(a) values is important for understanding structure and function relationships in proteins. Over the past two decades, theoretical methods for pK(a) calculations have been mainly based on macroscopic models, in which the protein is considered as a low-dielectric cavity embedded in a high-dielectric continuum. In recent years, constant pH molecular dynamics methods have been developed based on a microscopic description of the protein. We describe here the methodology of continuous constant pH molecular dynamics (CPHMD), which has emerged as one of the most robust and accurate tools for predicting protein pK(a)s and for the study of pH-modulated conformational dynamics. We illustrate the utility of CPHMD by the calculation of pK(a)s for surface residues in ribonuclease A, buried residues in staphylococcal nuclease, and titratable groups in the intrinsically flexible protein α-lactalbumin. We will compare the CPHMD results with experimental data as well as calculations from PB-based and empirical methods. These examples demonstrate the accuracy and robustness of the CPHMD method and its ability to capture the correlation between ionization equilibria and conformational dynamics as well as the local dielectric response to structural rearrangement. Finally, we discuss future improvement of the CPHMD method.

Integration of cell-free protein coexpression with an enzyme-linked immunosorbent assay enables rapid analysis of protein-protein interactions directly from DNA

Assays that integrate detection of binding with cell-free protein expression directly from DNA can dramatically increase the pace at which protein-protein interactions (PPIs) can be analyzed by mutagenesis. In this study, we present a method that combines in vitro protein production with an enzyme-linked immunosorbent assay (ELISA) to measure PPIs. This method uses readily available commodity instrumentation and generic antibody-affinity tag interactions. It is straightforward and rapid to execute, enabling many interactions to be assessed in parallel. In traditional ELISAs, reporter complexes are assembled stepwise with one layer at a time. In the method presented here, all the members of the reporter complex are present and assembled together. The signal strength is dependent on all the intercomponent interaction affinities and concentrations. Although this assay is straightforward to execute, establishing proper conditions and analysis of the results require a thorough understanding of the processes that determine the signal strength. The formation of the fully assembled reporter sandwich can be modeled as a competition between Langmuir adsorption isotherms for the immobilized components and binding equilibria of the solution components. We have shown that modeling this process provides semiquantitative understanding of the effects of affinity and concentration and can guide strategies for the development of experimental protocols. We tested the method experimentally using the interaction between a synthetic ankyrin repeat protein (Off7) and maltose-binding protein. Measurements obtained for a collection of alanine mutations in the interface between these two proteins demonstrate that a range of affinities can be analyzed.

Architecture of the flagellar rotor

Rotation and switching of the bacterial flagellum depends on a large rotor-mounted protein assembly composed of the proteins FliG, FliM and FliN, with FliG most directly involved in rotation. The crystal structure of a complex between the central domains of FliG and FliM, in conjunction with several biochemical and molecular-genetic experiments, reveals the arrangement of the FliG and FliM proteins in the rotor. A stoichiometric mismatch between FliG (26 subunits) and FliM (34 subunits) is explained in terms of two distinct positions for FliM: one where it binds the FliG central domain and another where it binds the FliG C-terminal domain. This architecture provides a structural framework for addressing the mechanisms of motor rotation and direction switching and for unifying the large body of data on motor performance. Recently proposed alternative models of rotor assembly, based on a subunit contact observed in crystals, are not supported by experiment.

Functional dissection of adenylate cyclase R, an inducer of spore encapsulation

Cyclic AMP acting on protein kinase A controls sporulation and encystation in social and solitary amoebas. In Dictyostelium discoideum, adenylate cyclase R (ACR), is essential for spore encapsulation. In addition to its cyclase (AC) domain, ACR harbors seven transmembrane helices, a histidine kinase domain, and two receiver domains. We investigated the role of these domains in the regulation of AC activity. Expression of an ACR-YFP fusion protein in acr(-) cells rescued their sporulation defective phenotype and revealed that ACR is associated with the nuclear envelope and endoplasmic reticulum. Loss of the transmembrane helices (ΔTM) caused a 60% reduction of AC activity, but ΔTM-ACR still rescued the acr(-) phenotype. The isolated AC domain was properly expressed but inactive. Mutation of three essential ATP-binding residues in the histidine kinase domain did not affect the AC activity or phenotypic rescue. Mutation of the essential phosphoryl-accepting aspartate in receivers 1, 2, or both had only modest effects on AC activity and did not affect phenotypic rescue, indicating that AC activity is not critically regulated by phosphorelay. Remarkably, the dimerizing histidine phosphoacceptor subdomain, which in ACR lacks the canonical histidine for autophosphorylation, was essential for AC activity. Transformation of wild-type cells with an ACR allele (ΔCRA) that is truncated after this domain inhibited AC activity of endogenous ACR and replicated the acr(-) phenotype. Combined with the observation that the isolated AC domain was inactive, the dominant-negative effect of ΔCRA strongly suggests that the defunct phosphoacceptor domain acquired a novel role in enforcing dimerization of the AC domain.

Prediction of protein-protein interaction sites using electrostatic desolvation profiles

Protein-protein complex formation involves removal of water from the interface region. Surface regions with a small free energy penalty for water removal or desolvation may correspond to preferred interaction sites. A method to calculate the electrostatic free energy of placing a neutral low-dielectric probe at various protein surface positions has been designed and applied to characterize putative interaction sites. Based on solutions of the finite-difference Poisson equation, this method also includes long-range electrostatic contributions and the protein solvent boundary shape in contrast to accessible-surface-area-based solvation energies. Calculations on a large set of proteins indicate that in many cases (>90%), the known binding site overlaps with one of the six regions of lowest electrostatic desolvation penalty (overlap with the lowest desolvation region for 48% of proteins). Since the onset of electrostatic desolvation occurs even before direct protein-protein contact formation, it may help guide proteins toward the binding region in the final stage of complex formation. It is interesting that the probe desolvation properties associated with residue types were found to depend to some degree on whether the residue was outside of or part of a binding site. The probe desolvation penalty was on average smaller if the residue was part of a binding site compared to other surface locations. Applications to several antigen-antibody complexes demonstrated that the approach might be useful not only to predict protein interaction sites in general but to map potential antigenic epitopes on protein surfaces.

3D-QSAR studies on triclosan derivatives as Plasmodium falciparum enoyl acyl carrier reductase inhibitors

3D-QSAR studies were carried out on a training set of 53 structurally highly diverse analogues of triclosan to investigate the correlation of the structural properties of triclosan derivatives with the inhibition of the activity of enoyl acyl carrier protein reductase in Plasmodium falciparum (PfENR) by employing Comparative Molecular Field Analysis (CoMFA) and Comparative Molecular Similarity Indices Analysis (CoMSIA). The crystal structure bound conformation of triclosan, was used as a template for aligning molecules. The probable binding mode conformations of other inhibitors were explored according to molecular docking and molecular mechanics poisson-boltzmann surface area (MM/PBSA) solvation free energy estimation methods using grid based linear Poisson-Boltzmann calculations. Predictive 3D-QSAR models, established using routine database alignment rule based on crystallographic-bound conformation of template molecule, produced statistically significant results with cross-validated r2 cv values of 0.64 and 0.54 and non-cross-validated r2 ncv values of 0.96 and 0.97 for CoMFA and CoMSIA models, respectively. The statistically significant models were validated by a test set of nine compounds with predictive r(2) values of 0.534 and 0.765 for CoMFA and CoMSIA respectively. Our QSAR model is able to successfully explain the geometric and electrostatic complementarities between ligands and receptor and provides useful guidelines to design novel triclosan derivatives as Plasmodium falciparum enoyl acyl carrier reductase inhibitors.

Chemotaxis signaling protein CheY binds to the rotor protein FliN to control the direction of flagellar rotation in Escherichia coli

The direction of rotation of the Escherichia coli flagellum is controlled by an assembly called the switch complex formed from multiple subunits of the proteins FliG, FliM, and FliN. Structurally, the switch complex corresponds to a drum-shaped feature at the bottom of the basal body, termed the C-ring. Stimulus-regulated reversals in flagellar motor rotation are the basis for directed movement such as chemotaxis. In E. coli, the motors turn counterclockwise (CCW) in their default state, allowing the several filaments on a cell to join together in a bundle and propel the cell smoothly forward. In response to the chemotaxis signaling molecule phospho-CheY (CheY(P)), the motors can switch to clockwise (CW) rotation, causing dissociation of the filament bundle and reorientation of the cell. CheY(P) has previously been shown to bind to a conserved segment near the N terminus of FliM. Here, we show that this interaction serves to capture CheY(P) and that the switch to CW rotation involves the subsequent interaction of CheY(P) with FliN. FliN is located at the bottom of the C-ring, in close association with the C-terminal domain of FliM (FliM(C)), and the switch to CW rotation has been shown to involve relative movement of FliN and FliM(C). Using a recently developed structural model for the FliN/FliM(C) array, and the CheY(P)-binding site here identified on FliN, we propose a mechanism by which CheY(P) binding could induce the conformational switch to CW rotation.

Acetylation represses the binding of CheY to its target proteins

The ability of CheY, the response regulator of bacterial chemotaxis, to generate clockwise rotation is regulated by two covalent modifications - phosphorylation and acetylation. While the function and signal propagation of the former are widely understood, the mechanism and role of the latter are still obscure. To obtain information on the function of this acetylation, we non-enzymatically acetylated CheY to a level similar to that found in vivo, and examined its binding to its kinase CheA, its phosphatase CheZ and the switch protein FliM - its target at the flagellar switch complex. Acetylation repressed the binding to all three proteins. These results suggest that both phosphorylation and acetylation determine CheY's ability to bind to its target proteins, thus providing two levels of regulation, fast and slow respectively. The fast level is modulated by environmental signals (e.g. chemotactic and thermotactic stimuli). The slow one is regulated by the metabolic state of the cell and it determines, at each metabolic state, the fraction of CheY molecules that can participate in signalling.

CobB regulates Escherichia coli chemotaxis by deacetylating the response regulator CheY

The silent information regulator (Sir2) family proteins are NAD+-dependent deacetylases. Although a few substrates have been identified, functions of the bacteria Sir2-like protein (CobB) still remain unclear. Here the role of CobB on Escherichia coli chemotaxis was investigated. We used Western blotting and mass spectrometry to show that the response regulator CheY is a substrate of CobB. Surface plasmon resonance (SPR) indicated that acetylation affects the interaction between CheY and the flagellar switch protein FliM. The presence of intact flagella in knockout strains DeltacobB, Deltaacs, Delta(cobB) Delta(acs), Delta(cheA) Delta(cheZ), Delta(cheA) Delta(cheZ) Delta(cobB) and Delta(cheA) Delta(cheZ) Delta(acs) was confirmed by electron microscopy. Genetic analysis of these knockout strains showed that: (i) the DeltacobB mutant exhibited reduced responses to chemotactic stimuli in chemotactic assays, whereas the Deltaacs mutant was indistinguishable from the parental strain, (ii) CheY from the DeltacobB mutant showed a higher level of acetylation, indicating that CobB can mediate the deacetylation of CheY in vivo, and (iii) deletion of cobB reversed the phenotype of Delta(cheA) Delta(cheZ). Our findings suggest that CobB regulates E. coli chemotaxis by deacetylating CheY. Thus a new function of bacterial cobB was identified and also new insights of regulation of bacterial chemotaxis were provided.

Topological frustration in beta alpha-repeat proteins: sequence diversity modulates the conserved folding mechanisms of alpha/beta/alpha sandwich proteins

The thermodynamic hypothesis of Anfinsen postulates that structures and stabilities of globular proteins are determined by their amino acid sequences. Chain topology, however, is known to influence the folding reaction, in that motifs with a preponderance of local interactions typically fold more rapidly than those with a larger fraction of nonlocal interactions. Together, the topology and sequence can modulate the energy landscape and influence the rate at which the protein folds to the native conformation. To explore the relationship of sequence and topology in the folding of beta alpha-repeat proteins, which are dominated by local interactions, we performed a combined experimental and simulation analysis on two members of the flavodoxin-like, alpha/beta/alpha sandwich fold. Spo0F and the N-terminal receiver domain of NtrC (NT-NtrC) have similar topologies but low sequence identity, enabling a test of the effects of sequence on folding. Experimental results demonstrated that both response-regulator proteins fold via parallel channels through highly structured submillisecond intermediates before accessing their cis prolyl peptide bond-containing native conformations. Global analysis of the experimental results preferentially places these intermediates off the productive folding pathway. Sequence-sensitive Gō-model simulations conclude that frustration in the folding in Spo0F, corresponding to the appearance of the off-pathway intermediate, reflects competition for intra-subdomain van der Waals contacts between its N- and C-terminal subdomains. The extent of transient, premature structure appears to correlate with the number of isoleucine, leucine, and valine (ILV) side chains that form a large sequence-local cluster involving the central beta-sheet and helices alpha2, alpha 3, and alpha 4. The failure to detect the off-pathway species in the simulations of NT-NtrC may reflect the reduced number of ILV side chains in its corresponding hydrophobic cluster. The location of the hydrophobic clusters in the structure may also be related to the differing functional properties of these response regulators. Comparison with the results of previous experimental and simulation analyses on the homologous CheY argues that prematurely folded unproductive intermediates are a common property of the beta alpha-repeat motif.

Quorum sensing-controlled Evr regulates a conserved cryptic pigment biosynthetic cluster and a novel phenomycin-like locus in the plant pathogen, Pectobacterium carotovorum

Pectobacterium carotovorum SCRI193 is a phytopathogenic Gram-negative bacterium. In this study, we have identified a novel cryptic pigment biosynthetic locus in P. carotovorum SCRI193 which we have called the Pectobacterium orange pigment (pop) cluster. The pop cluster is flanked by two tRNA genes and contains genes that encode non-ribosomal peptide synthases and polyketide synthase and produces a negatively charged polar orange pigment. Orange pigment production is activated when an adjacent transcriptional activator sharing sequence similarity with the Erwinia virulence regulator (Evr) is overexpressed. Evr was shown to positively activate its own transcription and that of the pigment biosynthetic genes and an unlinked locus encoding a phenomycin homologue. In addition, the expression of Evr and orange pigment production was shown to be regulated by N-(3-oxohexanoyl)-HSL (OHHL) quorum sensing and have a virulence phenotype in potato. Finally, by comparative genomics and Southern blotting we demonstrate that this pigment biosynthetic cluster is present in multiple P. carotovorum spp., Pectobacterium brasiliensis 1692 and a truncated version of the cluster is present in Pectobacterium atrosepticum. The conserved nature of this cluster in P. carotovorum and P. brasiliensis suggests that the pop cluster has an important function in these broad-host-range soft rotting bacteria, which is no longer required in the narrow-host-range P. atrosepticum SCRI1043.

Auxiliary phosphatases in two-component signal transduction

Signal termination in two-component systems occurs by loss of the phosphoryl group from the response regulator protein. This review explores our current understanding of the structures, catalytic mechanisms and means of regulation of the known families of phosphatases that catalyze response regulator dephosphorylation. The CheZ and CheC/CheX/FliY families, despite different overall structures, employ identical catalytic strategies using an amide side chain to orient a water molecule for in-line attack of the aspartyl phosphate. Spo0E phosphatases contain sequence and structural features that suggest a strategy similar to the chemotaxis phosphatases but the mechanism used by the Rap phosphatases is not yet elucidated. Identification of features shared by phosphatase families may aid in the identification of currently unrecognized classes of response regulator phosphatases.

Molecular strategies for phosphorylation-mediated regulation of response regulator activity

Response regulator (RR) proteins exploit different molecular surfaces in their inactive and active conformations for a variety of regulatory intramolecular and/or intermolecular protein-protein interactions that either inhibit or activate effector domain activities. This versatile strategy enables numerous regulatory mechanisms among RRs. The recent accumulation of structures of inactive and active forms of multidomain RRs and RR complexes has revealed many different domain arrangements that have provided insight into regulatory mechanisms. Although diversity is the rule, even among subfamily members containing homologous domains, several structural modes of interaction and mechanisms of regulation recur frequently. These themes involve interactions at the alpha4-beta5-alpha5 face of the receiver domain, modes of dimerization of receiver domains, and inhibitory or activating heterodomain interactions.

Structure of SurE protein from Aquifex aeolicus VF5 at 1.5 A resolution

SurE is a stationary-phase survival protein found in bacteria, eukaryotes and archaea that exhibits a divalent-metal-ion-dependent phosphatase activity and acts as a nucleotidase and polyphosphate phosphohydrolase. The structure of the SurE protein from the hyperthermophile Aquifex aeolicus has been solved at 1.5 A resolution using molecular replacement with one dimer in the asymmetric unit and refined to an R factor of 15.6%. The crystal packing reveals that two dimers assemble to form a tetramer, although gel-filtration chromatography showed the presence of only a dimer in solution. The phosphatase active-site pocket was occupied by sulfate ions from the crystallization medium.

Three (and more) component regulatory systems - auxiliary regulators of bacterial histidine kinases

Two-component signal transduction (TCST) is the most prevalent mechanism employed by microbes to sense and respond to environmental changes. It is characterized by the signal-induced transfer of phosphate from a sensor histidine kinase (HK) to a response regulator (RR), resulting in a cellular response. An emerging theme in the field of TCST signalling is the discovery of auxiliary factors, distinct from the HK and RR, which are capable of influencing phosphotransfer. One group of TCST auxiliary proteins accomplishes this task by acting on HKs. Auxiliary regulators of HKs are widespread and have been identified in all cellular compartments, where they can influence HK activity through interactions with the sensing, transmembrane or enzymatic domains of the HK. The effects of an auxiliary regulator are controlled by its regulated expression, modification and/or through ligand binding. Ultimately, auxiliary regulators can connect a given TCST system to other regulatory networks in the cell or result in regulation of the TCST system in response to an expanded range of stimuli. The studies highlighted in this review draw attention to an emerging view of bacterial TCST systems as core signalling units upon which auxiliary factors act.

Subunit organization and reversal-associated movements in the flagellar switch of Escherichia coli

Bacterial flagella contain a rotor-mounted protein complex termed the switch complex that functions in flagellar assembly, rotation, and clockwise/counterclockwise direction control. In Escherichia coli and Salmonella, the switch complex contains the proteins FliG, FliM, and FliN and corresponds structurally with the C-ring in the flagellar basal body. Certain features of subunit organization in the switch complex have been deduced previously, but details of subunit organization in the lower part of the C-ring and the molecular movements responsible for motor switching remain unclear. In this study, we use cross-linking, binding, and mutational experiments to examine subunit organization in the bottom of the C-ring and to probe movements that occur upon switching. The results show that FliN tetramers alternate with FliM C-terminal domains to form the bottom of the C-ring in an arrangement that closely reproduces the major features observed in electron microscopic reconstructions. When motors were switched to clockwise rotation by a repellent stimulus, cross-link yields were altered in a pattern indicating relative movement of FliN and FliM(C). These results are discussed in the framework of a structurally grounded hypothesis for the switching mechanism.

Functional and small-angle X-ray scattering studies of a new stationary phase survival protein E (SurE) from Xylella fastidiosa--evidence of allosteric behaviour

The genome data of bacterium Xylella fastidiosa strain 9a5c has identified several orfs related to its phytopathogenic adaptation and survival. Among these genes, the surE codifies a survival protein E (XfSurE) whose function is not so well understood, but functional assays in Escherichia coli revealed nucleotidase and exopolyphosphate activity. In the present study, we report the XfSurE protein overexpression in E. coli and its purification. The overall secondary structure was analyzed by CD. Small-angle X-ray scattering and gel filtration techniques demonstrated that the oligomeric state of the protein in solution is a tetramer. In addition, functional kinetics experiments were carried out with several monophosphate nucleoside substrates and revealed a highly positive cooperativity. An allosteric mechanism involving torsion movements in solution is proposed to explain the cooperative behaviour of XfSurE. This is the first characterization of a SurE enzyme from a phytopathogen organism and, to our knowledge, the first solution structure of a SurE protein to be described.

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(6). 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(3)(-)) 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.

A molecular mechanism of bacterial flagellar motor switching

The high-resolution structures of nearly all the proteins that comprise the bacterial flagellar motor switch complex have been solved; yet a clear picture of the switching mechanism has not emerged. Here, we used NMR to characterize the interaction modes and solution properties of a number of these proteins, including several soluble fragments of the flagellar motor proteins FliM and FliG, and the response-regulator CheY. We find that activated CheY, the switch signal, binds to a previously unidentified region of FliM, adjacent to the FliM-FliM interface. We also find that activated CheY and FliG bind with mutual exclusivity to this site on FliM, because their respective binding surfaces partially overlap. These data support a model of CheY-driven motor switching wherein the binding of activated CheY to FliM displaces the carboxy-terminal domain of FliG (FliGC) from FliM, modulating the FliGC-MotA interaction, and causing the motor to switch rotational sense as required for chemotaxis.

Structure-based design of a T-cell receptor leads to nearly 100-fold improvement in binding affinity for pepMHC

T-cell receptors (TCRs) are proteins that recognize peptides from foreign proteins bound to the major histocompatibility complex (MHC) on the surface of an antigen-presenting cell. This interaction enables the T cells to initiate a cell-mediated immune response to terminate cells displaying the foreign peptide on their MHC. Naturally occurring TCRs have high specificity but low affinity toward the peptide-MHC (pepMHC) complex. This prevents the usage of solubilized TCRs for diagnosis and treatment of viral infections or cancers. Efforts to enhance the binding affinity of several TCRs have been reported in recent years, through randomized libraries and in vitro selection. However, there have been no reported efforts to enhance the affinity via structure-based design, which allows more control and understanding of the mechanism of improvement. Here, we have applied structure-based design to a human TCR to improve its pepMHC binding. Our design method evolved based on iterative steps of prediction, testing, and generating more predictions based on the new data. The final design function, named ZAFFI, has a correlation of 0.77 and average error of 0.35 kcal/mol with the binding free energies of 26 point mutations for this system that we measured by surface plasmon resonance (SPR). Applying the filter that we developed to remove nonbinding predictions, this correlation increases to 0.85, and the average error decreases to 0.3 kcal/mol. Using this algorithm, we predicted and tested several point mutations that improved binding, with one giving over sixfold binding improvement. Four of the point mutations that improved binding were then combined to give a mutant TCR that binds the pepMHC 99 times more strongly than the wild-type TCR.

Single domain response regulators: molecular switches with emerging roles in cell organization and dynamics

Single domain response regulators (SD-RRs) are signaling components of two-component phosphorylation pathways that harbor a phosphoryl receiver domain but lack a dedicated output domain. The Escherichia coli protein CheY, the paradigm member of this family, regulates chemotaxis by relaying information between chemoreceptors and the flagellar motor switch. New data provide a more complex picture of CheY-mediated motility control in several bacteria and suggest diverging mechanisms in control of cellular motors. Moreover, advances have been made in understanding cellular functions of SD-RRs beyond chemotaxis. We review recent reports indicating that SD-RRs constitute a family of versatile molecular switches that contribute to cellular organization and dynamics as spatial organizer and/or as allosteric regulators of histidine protein kinases.

Structural and functional studies on a mesophilic stationary phase survival protein (Sur E) from Salmonella typhimurium

SurE, the stationary-phase survival protein of Salmonella typhimurium, forms part of a stress survival operon regulated by the stationary-phase RNA polymerase alternative sigma factor. SurE is known to improve bacterial viability during stress conditions. It functions as a phosphatase specific to nucleoside monophosphates. In the present study we reported the X-ray crystal structure of SurE from Salmonella typhimurium. The protein crystallized in two forms: orthorhombic F222; and monoclinic C2. The two structures were determined to resolutions of 1.7 and 2.7 A, respectively. The protein exists as a domain-swapped dimer. The residue D230 is involved in several interactions that are probably crucial for domain swapping. A divalent metal ion is found at the active site of the enzyme, which is consistent with the divalent metal ion-dependent activity of the enzyme. Interactions of the conserved DD motif present at the N-terminus with the phosphate and the Mg(2+) present in the active site suggest that these residues play an important role in enzyme activity. The divalent metal ion specificity and the kinetic constants of SurE were determined using the generic phosphatase substrate para-nitrophenyl phosphate. The enzyme was inactive in the absence of divalent cations and was most active in the presence of Mg(2+). Thermal denaturation studies showed that S. typhimurium SurE is much less stable than its homologues and an attempt was made to understand the molecular basis of the lower thermal stability based on solvation free-energy. This is the first detailed crystal structure analysis of SurE from a mesophilic organism.

Interpreting the Coulomb-field approximation for generalized-Born electrostatics using boundary-integral equation theory

The importance of molecular electrostatic interactions in aqueous solution has motivated extensive research into physical models and numerical methods for their estimation. The computational costs associated with simulations that include many explicit water molecules have driven the development of implicit-solvent models, with generalized-Born (GB) models among the most popular of these. In this paper, we analyze a boundary-integral equation interpretation for the Coulomb-field approximation (CFA), which plays a central role in most GB models. This interpretation offers new insights into the nature of the CFA, which traditionally has been assessed using only a single point charge in the solute. The boundary-integral interpretation of the CFA allows the use of multiple point charges, or even continuous charge distributions, leading naturally to methods that eliminate the interpolation inaccuracies associated with the Still equation. This approach, which we call boundary-integral-based electrostatic estimation by the CFA (BIBEE/CFA), is most accurate when the molecular charge distribution generates a smooth normal displacement field at the solute-solvent boundary, and CFA-based GB methods perform similarly. Conversely, both methods are least accurate for charge distributions that give rise to rapidly varying or highly localized normal displacement fields. Supporting this analysis are comparisons of the reaction-potential matrices calculated using GB methods and boundary-element-method (BEM) simulations. An approximation similar to BIBEE/CFA exhibits complementary behavior, with superior accuracy for charge distributions that generate rapidly varying normal fields and poorer accuracy for distributions that produce smooth fields. This approximation, BIBEE by preconditioning (BIBEE/P), essentially generates initial guesses for preconditioned Krylov-subspace iterative BEMs. Thus, iterative refinement of the BIBEE/P results recovers the BEM solution; excellent agreement is obtained in only a few iterations. The boundary-integral-equation framework may also provide a means to derive rigorous results explaining how the empirical correction terms in many modern GB models significantly improve accuracy despite their simple analytical forms.

Accurate solution of multi-region continuum biomolecule electrostatic problems using the linearized Poisson-Boltzmann equation with curved boundary elements

We present a boundary-element method (BEM) implementation for accurately solving problems in biomolecular electrostatics using the linearized Poisson-Boltzmann equation. Motivating this implementation is the desire to create a solver capable of precisely describing the geometries and topologies prevalent in continuum models of biological molecules. This implementation is enabled by the synthesis of four technologies developed or implemented specifically for this work. First, molecular and accessible surfaces used to describe dielectric and ion-exclusion boundaries were discretized with curved boundary elements that faithfully reproduce molecular geometries. Second, we avoided explicitly forming the dense BEM matrices and instead solved the linear systems with a preconditioned iterative method (GMRES), using a matrix compression algorithm (FFTSVD) to accelerate matrix-vector multiplication. Third, robust numerical integration methods were employed to accurately evaluate singular and near-singular integrals over the curved boundary elements. Fourth, we present a general boundary-integral approach capable of modeling an arbitrary number of embedded homogeneous dielectric regions with differing dielectric constants, possible salt treatment, and point charges. A comparison of the presented BEM implementation and standard finite-difference techniques demonstrates that for certain classes of electrostatic calculations, such as determining absolute electrostatic solvation and rigid-binding free energies, the improved convergence properties of the BEM approach can have a significant impact on computed energetics. We also demonstrate that the improved accuracy offered by the curved-element BEM is important when more sophisticated techniques, such as nonrigid-binding models, are used to compute the relative electrostatic effects of molecular modifications. In addition, we show that electrostatic calculations requiring multiple solves using the same molecular geometry, such as charge optimization or component analysis, can be computed to high accuracy using the presented BEM approach, in compute times comparable to traditional finite-difference methods.

Poisson-Boltzmann model analysis of binding mRNA cap analogues to the translation initiation factor eIF4E

The electrostatic free energy of binding of two analogues of the 5'-mRNA cap, differing in size and electric charge, to the wild type and mutated eukaryotic initiation factor eIF4E was computed using the finite difference solutions to the Poisson-Boltzmann equation. Two definitions of the solute-solvent dielectric boundary were used: van der Waals model, solvent exclusion (SE) model. The computed electrostatic energies were supplemented by estimations of the non polar and entropic contributions. A comparison with experimental data for the investigated systems was done. It appears that the SE model with additional contribution fits experimental findings better than the van der Waals model does.

The structures of T87I phosphono-CheY and T87I/Y106W phosphono-CheY help to explain their binding affinities to the FliM and CheZ peptides

CheY is a response regulator in bacterial chemotaxis. Escherichia coli CheY mutants T87I and T87I/Y106W CheY are phosphorylatable on Asp57 but unable to generate clockwise rotation of the flagella. To understand this phenotype in terms of structure, stable analogs of the two CheY-P mutants were synthesized: T87I phosphono-CheY and T87I phosphono-CheY. Dissociation constants for peptides derived from flagellar motor protein FliM and phosphatase CheZ were determined for phosphono-CheY and the two mutants. The peptides bind phosphono-CheY almost as strongly as CheY-P; however, they do not bind T87I phosphono-CheY or T87I/Y106W phosphono-CheY, implying that the mutant proteins cannot bind FliM or CheZ tightly in vivo. The structures of T87I phosphono-CheY and T87I/Y106W phosphono-CheY were solved to resolutions of 1.8 and 2.4A, respectively. The increased bulk of I87 forces the side-chain of Y106 or W106, into a more solvent-accessible conformation, which occludes the peptide-binding site.

Subdomain competition, cooperativity, and topological frustration in the folding of CheY

The folding of multidomain proteins often proceeds in a hierarchical fashion with individual domains folding independent of one another. A large single-domain protein, however, can consist of multiple modules whose folding may be autonomous or interdependent in ways that are unclear. We used coarse-grained simulations to explore the folding landscape of the two-subdomain bacterial response regulator CheY. Thermodynamic and kinetic characterization shows the landscape to be highly analogous to the four-state landscape reported for another two-subdomain protein, T4 lysozyme. An on-pathway intermediate structured in the more stable nucleating subdomain was observed, as were transient states frustrated in off-pathway contacts prematurely structured in the weaker subdomain. Local unfolding, or backtracking, was observed in the frustrated state before the native conformation could be reached. Nonproductive frustration was attributable to competition for van der Waals contacts between the two subdomains. In an accompanying article, stopped-flow kinetic measurements support an off-pathway burst-phase intermediate, seemingly consistent with our prediction of early frustration in the folding landscape of CheY. Comparison of the folding mechanisms for CheY, T4 lysozyme, and interleukin-1 beta leads us to postulate that subdomain competition is a general feature of large single-domain proteins with multiple folding modules.

Allostery and cooperativity revisited

Although phenomenlogical models that account for cooperativity in allosteric systems date back to the early and mid-60's (e.g., the KNF and MWC models), there is resurgent interest in the topic due to the recent experimental and computational studies that attempted to reveal, at an atomistic level, how allostery actually works. In this review, using systems for which atomistic simulations have been carried out in our groups as examples, we describe the current understanding of allostery, how the mechanisms go beyond the classical MWC/Pauling-KNF descriptions, and point out that the "new view" of allostery, emphasizing "population shifts," is, in fact, an "old view." The presentation offers not only an up-to-date description of allostery from a theoretical/computational perspective, but also helps to resolve several outstanding issues concerning allostery.

Evaluation of electrostatic interactions

Described here are several computational procedures for the analysis of electrostatic interactions in molecular complexes, all based on a continuum model of solvation. The first section describes how to compute the residual potential, a measure of how electrostatically complementary a ligand is for its receptor. The second procedure describes electrostatic component analysis, a method by which the electrostatic contribution to the binding free energy can be broken up into terms directly attributable to individual chemical groups. Finally, electrostatic affinity optimization is described. This procedure is particularly useful in determining what portions of a ligand are the most suboptimal, and thus provide the greatest opportunity for the design of improvements.

Crystal structural analysis of protein-protein interactions drastically destabilized by a single mutation

The complex of barnase (bn) and barstar (bs), which has been widely studied as a model for quantitative analysis of protein-protein interactions, is significantly destabilized by a single mutation, namely, bs Asp39 --> Ala, which corresponds to a change of 7.7 kcal x mol(-1) in the free energy of binding. However, there has been no structural information available to explain such a drastic destabilization. In the present study, we determined the structure of the mutant complex at 1.58 A resolution by X-ray crystallography. The complex was similar to the wild-type complex in terms of overall and interface structures; however, the hydrogen bond network mediated by water molecules at the interface was significantly different. Several water molecules filled the cavity created by the mutation and consequently caused rearrangement of the hydrated water molecules at the interface. The water molecules were redistributed into a channel-like structure that penetrated into the complex. Furthermore, molecular dynamics simulations showed that the mutation increased the mobility of water molecules at the interface. Since such a drastic change in hydration was not observed in other mutant complexes of bn and bs, the significant destabilization of the interaction may be due to this channel-like structure of hydrated water molecules.

Computational design and experimental study of tighter binding peptides to an inactivated mutant of HIV-1 protease

Drug resistance in HIV-1 protease, a barrier to effective treatment, is generally caused by mutations in the enzyme that disrupt inhibitor binding but still allow for substrate processing. Structural studies with mutant, inactive enzyme, have provided detailed information regarding how the substrates bind to the protease yet avoid resistance mutations; insights obtained inform the development of next generation therapeutics. Although structures have been obtained of complexes between substrate peptide and inactivated (D25N) protease, thermodynamic studies of peptide binding have been challenging due to low affinity. Peptides that bind tighter to the inactivated protease than the natural substrates would be valuable for thermodynamic studies as well as to explore whether the structural envelope observed for substrate peptides is a function of weak binding. Here, two computational methods-namely, charge optimization and protein design-were applied to identify peptide sequences predicted to have higher binding affinity to the inactivated protease, starting from an RT-RH derived substrate peptide. Of the candidate designed peptides, three were tested for binding with isothermal titration calorimetry, with one, containing a single threonine to valine substitution, measured to have more than a 10-fold improvement over the tightest binding natural substrate. Crystal structures were also obtained for the same three designed peptide complexes; they show good agreement with computational prediction. Thermodynamic studies show that binding is entropically driven, more so for designed affinity enhanced variants than for the starting substrate. Structural studies show strong similarities between natural and tighter-binding designed peptide complexes, which may have implications in understanding the molecular mechanisms of drug resistance in HIV-1 protease.

Novel method for probing the specificity binding profile of ligands: applications to HIV protease

A detailed understanding of factors influencing the binding specificity of a ligand to a set of desirable targets and undesirable decoys is a key step in the design of potent and selective therapeutics. We have developed a general method for optimizing binding specificity in ligand-receptor complexes based on the theory of electrostatic charge optimization. This methodology can be used to tune the binding of a ligand to a panel of potential targets and decoys, along the continuum from narrow binding to only one partner to broad binding to the entire panel. Using HIV-1 protease as a model system, we probe specificity in three distinct ways. First, we probe interactions that could make the promiscuous protease inhibitor pepstatin more selective toward HIV-1 protease. Next, we study clinically approved HIV-1 protease inhibitors and probe ways to broaden the binding profiles toward both wild-type HIV-1 protease and drug-resistant mutants. Finally, we study a conformational ensemble of wild-type HIV-1 protease to 'design in' broad specificity to known drugs before resistance mutations arise. The results from this conformational ensemble were similar to those from the drug-resistant ensemble, suggesting the use of a conformational wild-type ensemble as a tool to develop escape-mutant-resistant inhibitors.

Expression, purification, crystallization and preliminary X-ray characterization of two crystal forms of stationary-phase survival E protein from Campylobacter jejuni

Survival E (SurE) protein from Campylobacter jejuni, a Gram-negative mesophile, has been overexpressed in Escherichia coli as a soluble protein, successfully purified and crystallized in two distinct crystal forms. The first form belongs to space group P2(1)2(1)2(1), with a tetramer in the asymmetric unit and unit-cell parameters a = 80.5, b = 119.0, c = 135.3 A. The second form belongs to space group C2, with unit-cell parameters a = 121.4, b = 47.1, c = 97.8 A, and contains a dimer in the asymmetric unit. Diffraction data have been collected from these crystal forms to 2.5 and 2.95 A resolution, respectively.

In vivo acetylation of CheY, a response regulator in chemotaxis of Escherichia coli

CheY, the excitatory response regulator in the chemotaxis system of Escherichia coli, can be modulated by two covalent modifications: phosphorylation and acetylation. Both modifications have been detected in vitro only. The role of CheY acetylation is still obscure, although it is known to be involved in chemotaxis and to occur in vitro by two mechanisms--acetyl-CoA synthetase-catalyzed transfer of acetyl groups from acetate to CheY and autocatalyzed transfer from AcCoA. Here, we succeeded in detecting CheY acetylation in vivo by three means--Western blotting with a specific anti-acetyl-lysine antibody, mass spectrometry, and radiolabeling with [(14)C]acetate in the presence of protein-synthesis inhibitor. Unexpectedly, the level and rate of CheY acetylation in vivo were much higher than that in vitro. Thus, before any treatment, 9-13% of the lysine residues were found acetylated, depending on the growth phase, meaning that, on average, essentially every CheY molecule was acetylated in vivo. This high level was mainly the outcome of autoacetylation. Addition of acetate caused an incremental increase in the acetylation level, in which acetyl-CoA synthetase was involved too. These findings may have far-reaching implications for the structure-function relationship of CheY.

Flagellar motility in bacteria structure and function of flagellar motor

Bacterial flagella are filamentous organelles that drive cell locomotion. They thrust cells in liquids (swimming) or on surfaces (swarming) so that cells can move toward favorable environments. At the base of each flagellum, a reversible rotary motor, which is powered by the proton- or the sodium-motive force, is embedded in the cell envelope. The motor consists of two parts: the rotating part, or rotor, that is connected to the hook and the filament, and the nonrotating part, or stator, that conducts coupling ion and is responsible for energy conversion. Intensive genetic and biochemical studies of the flagellum have been conducted in Salmonella typhimurium and Escherichia coli, and more than 50 gene products are known to be involved in flagellar assembly and function. The energy-coupling mechanism, however, is still not known. In this chapter, we survey our current knowledge of the flagellar system, based mostly on studies from Salmonella, E. coli, and marine species Vibrio alginolyticus, supplemented with distinct aspects of other bacterial species revealed by recent studies.

Crystal structure of a complex between the phosphorelay protein YPD1 and the response regulator domain of SLN1 bound to a phosphoryl analog

The crystal structure of the yeast SLN1 response regulator (RR) domain bound to both a phosphoryl analog [beryllium fluoride (BeF(3)(-))] and Mg(2+), in complex with its downstream phosphorelay signaling partner YPD1, has been determined at a resolution of 1.70 A. Comparisons between the BeF(3)(-)-activated complex and the unliganded (or apo) complex determined previously reveal modest but important differences. The SLN1-R1 x Mg(2+) x BeF(3)(-) structure from the complex provides evidence for the first time that the mechanism of phosphorylation-induced activation is highly conserved between bacterial RR domains and this example from a eukaryotic organism. Residues in and around the active site undergo slight rearrangements in order to form bonds with the essential divalent cation and fluorine atoms of BeF(3)(-). Two conserved switch-like residues (Thr1173 and Phe1192) occupy distinctly different positions in the apo versus BeF(3)(-)-bound structures, consistent with the "Y-T" coupling mechanism proposed for the activation of CheY and other bacterial RRs. Several loop regions and the alpha 4-beta 5-alpha 5 surface of the SLN1-R1 domain undergo subtle conformational changes ( approximately 1-3 A displacements relative to the apo structure) that lead to significant changes in terms of contacts that are formed with YPD1. Detailed structural comparisons of protein-protein interactions in the apo and BeF(3)(-)-bound complexes suggest at least a two-state equilibrium model for the formation of a transient encounter complex, in which phosphorylation of the RR promotes the formation of a phosphotransfer-competent complex. In the BeF(3)(-)-activated complex, the position of His64 from YPD1 needs to be within ideal distance of and in near-linear geometry with Asp1144 from the SLN1-R1 domain for phosphotransfer to occur. The ground-state structure presented here suggests that phosphoryl transfer will likely proceed through an associative mechanism involving the formation of a pentacoordinate phosphorus intermediate.

Integration of DIGE and bioinformatics analyses reveals a role of the antiobesity agent tungstate in redox and energy homeostasis pathways in brown adipose tissue

Our previous results demonstrated that tungstate decreased weight gain and adiposity in obese rats through increased thermogenesis and lipid oxidation, suggesting that brown adipose tissue was one of the targets of its antiobesity effect. To identify potential targets of tungstate, we used DIGE to compare brown adipose tissue protein extracts from the following experimental groups: untreated lean, tungstate-treated lean, untreated obese, and tungstate-treated obese rats. To distinguish direct targets of tungstate action from those that are secondary to body weight loss, we also included in the analysis an additional group consisting of obese rats that lose weight by caloric restriction. Hierarchical clustering of analysis of variance and t test contrasts clearly separated the different experimental groups. DIGE analysis identified 20 proteins as tungstate obesity direct targets involved in Krebs cycle, glycolysis, lipolysis and fatty acid oxidation, electron transport, and redox. Protein oxidation was decreased by tungstate treatment, confirming a role in redox processes; however, palmitate oxidation, as a measure of fatty acid beta-oxidation, was not altered by tungstate, thus questioning its putative function in fatty acid oxidation. Protein network analyses using Ingenuity Pathways Analysis highlighted peroxisome proliferator-activated receptor gamma coactivator 1alpha (PGC-1alpha) as a potential target. We confirmed by real time PCR that indeed tungstate up-regulates PGC-1alpha, and its major target, uncoupling protein 1, was also increased as shown by Western blot. These results illustrate the utility of proteomics and bioinformatics approaches to identify targets of obesity therapies and suggest that in brown adipose tissue tungstate modulates redox processes and increases energy dissipation through uncoupling and PGC-1alpha up-regulation, thus contributing to its overall antiobesity effect.

Robust protein protein interactions in crowded cellular environments

The capacity of proteins to interact specifically with one another underlies our conceptual understanding of how living systems function. Systems-level study of specificity in protein-protein interactions is complicated by the fact that the cellular environment is crowded and heterogeneous; interaction pairs may exist at low relative concentrations and thus be presented with many more opportunities for promiscuous interactions compared with specific interaction possibilities. Here we address these questions by using a simple computational model that includes specifically designed interacting model proteins immersed in a mixture containing hundreds of different unrelated ones; all of them undergo simulated diffusion and interaction. We find that specific complexes are quite robust to interference from promiscuous interaction partners only in the range of temperatures T(design) > T > T(rand). At T > T(design), specific complexes become unstable, whereas at T < T(rand), formation of specific complexes is suppressed by promiscuous interactions. Specific interactions can form only if T(design) > T(rand). This condition requires an energy gap between binding energy in a specific complex and set of binding energies between randomly associating proteins, providing a general physical constraint on evolutionary selection or design of specific interacting protein interfaces. This work has implications for our understanding of how the protein repertoire functions and evolves within the context of cellular systems.

Optimization of electrostatic interactions in protein-protein complexes

In this article, we present a statistical analysis of the electrostatic properties of 298 protein-protein complexes and 356 domain-domain structures extracted from the previously developed database of protein complexes (ProtCom, http://www.ces.clemson.edu/compbio/protcom). For each structure in the dataset we calculated the total electrostatic energy of the binding and its two components, Coulombic and reaction field energy. It was found that in a vast majority of the cases (>90%), the total electrostatic component of the binding energy was unfavorable. At the same time, the Coulombic component of the binding energy was found to favor the complex formation while the reaction field component of the binding energy opposed the binding. It was also demonstrated that the components in a wild-type (WT) structure are optimized/anti-optimized with respect to the corresponding distributions, arising from random shuffling of the charged side chains. The degree of this optimization was assessed through the Z-score of WT energy in respect to the random distribution. It was found that the Z-scores of Coulombic interactions peak at a considerably negative value for all 654 cases considered while the Z-score of the reaction field energy varied among different types of complexes. All these findings indicate that the Coulombic interactions within WT protein-protein complexes are optimized to favor the complex formation while the total electrostatic energy predominantly opposes the binding. This observation was used to discriminate WT structures among sets of structural decoys and showed that the electrostatic component of the binding energy is not a good discriminator of the WT; while, Coulombic or reaction field energies perform better depending upon the decoy set used.

Why are proteins charged? Networks of charge-charge interactions in proteins measured by charge ladders and capillary electrophoresis

Almost all proteins contain charged amino acids. While the function in catalysis or binding of individual charges in the active site can often be identified, it is less clear how to assign function to charges beyond this region. Are they necessary for solubility? For reasons other than solubility? Can manipulating these charges change the properties of proteins? A combination of capillary electrophoresis (CE) and protein charge ladders makes it possible to study the roles of charged residues on the surface of proteins outside the active site. This method involves chemical modification of those residues to generate a large number of derivatives of the protein that differ in charge. CE separates those derivatives into groups with the same number of modified charged groups. By studying the influence of charge on the properties of proteins using charge ladders, it is possible to estimate the net charge and hydrodynamic radius and to infer the role of charged residues in ligand binding and protein folding.

Crystal structure of the stationary phase survival protein SurE with metal ion and AMP

The stationary phase survival protein SurE is a metal ion-dependent phosphatase distributed among eubacteria, archaea, and eukaryotes. In Escherichia coli, SurE has activities as nucleotidase and exopolyphosphatase, and is thought to be involved in stress response. However, its physiological role and reaction mechanism are unclear. We report here the crystal structures of the tetramer of SurE from Thermus thermophilus HB8 (TtSurE) both alone and crystallized with Mn(2+) and substrate AMP. In the presence of Mn(2+) and AMP, differences between the protomers were observed in the active site and in the loop located near the active site; AMP-bound active sites with the loops in a novel open conformation were found in the two protomers, and AMP-free active sites with the loops in a conventional closed conformation were found in the other two protomers. The two loops in the open conformation are entwined with each other, and this entwining is suggested to be required for enzymatic activity by site-directed mutagenesis. TtSurE exists as an equilibrium mixture of dimer and tetramer in solution. The loop-entwined structure indicates that SurE acts as a tetramer. The structural features and the absence of negative cooperativity imply the half-of-the-sites reactivity mechanism resulting from a pre-existing tendency toward structural asymmetry.

Phosphorylation events implicating p38 and PI3K mediate tungstate-effects in MIN6 beta cells

Oral administration of sodium tungstate is an effective treatment for diabetes in animal models. Several lines of evidence indicate the pancreatic beta cell as one of the targets of tungstate action. Here, we examined the molecular mechanism by which this compound exerts its effects on the beta cell line MIN6. Tungstate treatment induced phosphorylation and subsequent activation of p38 and PI3K which in turn are implicated in tungstate PDX-1 nuclear localization and activation. Although no effect was observed in glucose-induced insulin secretion we found that tungstate activates basal insulin release, a process driven, at least in part, by activation of p38. These results show a direct involvement of p38 and PI3K phosphorylation in the mechanism of action of tungstate in the beta cell.

The use of a generalized born model for the analysis of protein conformational transitions: a comparative study with explicit solvent simulations for chemotaxis Y protein (CheY)

To investigate whether implicit solvent models are appropriate for mechanistic studies of conformational transition in proteins, a recently developed generalized Born model (GBSW) was applied to a small signaling protein, chemotaxis protein Y (CheY), with different combinations of the phosphorylation state and conformation of the system; the results were compared to explicit solvent simulations using a stochastic boundary condition. The subtle but distinct conformational transitions involved in CheY activation makes the system ideally suited for comparing implicit and explicit solvent models because these conformational transitions are potentially accessible in both types of simulations. The structural and dynamical properties analyzed include not only those localized to the active site region but also throughout the protein, such as sidechain methyl group order parameters, backbone hydrogen bonding lifetime and occupancy as well as principal components of the trajectories. Overall, many properties were well reproduced by the GBSW simulations when compared with the explicit solvent calculations, although a number of observations consistently point to the suggestion that the current parameterization of the GBSW model tends to overestimate hydrogen-bonding interactions involving both charged groups and (charge-neutral) backbone atoms. This deficiency led to overstabilization of certain secondary structural motifs and more importantly, qualitatively different behaviors for the active site groups (Thr 87, Ala 88, the beta4-alpha4 loop) in response to phosphorylation, when compared with explicit solvent simulations. The current study highlights the value of carrying out both explicit and implicit solvent simulations for complementary mechanistic insights in the analysis of conformational transition in biomolecules.

Coupled atomic charge selectivity for optimal ligand-charge distributions at protein binding sites

Charge optimization as a tool for both analyzing and enhancing binding electrostatics has become an attractive approach over the past few years. An interesting feature of this method for molecular design is that it provides not only the optimal charge magnitudes, but also the selectivity of a particular atomic center for its optimal charge. The current approach to compute the charge selectivity at a given atomic center of a ligand in a particular binding process is based on the binding-energy cost incurred upon the perturbation of the optimal charge distribution by a unit charge at the given atomic center, while keeping the other atomic partial charges at their optimal values. A limitation of this method is that it does not take into account the possible concerted changes in the other atomic charges that may incur a lower energetic cost than perturbing a single charge. Here, we describe a novel approach for characterizing charge selectivity in a concerted manner, taking into account the coupling between the ligand charge centers in the binding process. We apply this novel charge selectivity measure to the celecoxib molecule, a nonsteroidal anti-inflammatory agent binding to cyclooxygenase 2 (COX2), which has been recently shown to also exhibit cross-reactivity toward carbonic anhydrase II (CAII), to which it binds with nanomolar affinity. The uncoupled and coupled charge selectivity profiles over the atomic centers of the celecoxib ligand, binding independently to COX2 and CAII, are analyzed comparatively and rationalized with respect to available experimental data. Very different charge selectivity profiles are obtained for the uncoupled versus coupled selectivity calculations.

The bidirectional polar and unidirectional lateral flagellar motors of Vibrio alginolyticus are controlled by a single CheY species

The bacterial flagellar motor is an elaborate molecular machine that converts ion-motive force into mechanical force (rotation). One of its remarkable features is its swift switching of the rotational direction or speed upon binding of the response regulator phospho-CheY, which causes the changes in swimming that achieve chemotaxis. Vibrio alginolyticus has dual flagellar systems: the Na(+)-driven polar flagellum (Pof) and the H(+)-driven lateral flagella (Laf), which are used for swimming in liquid and swarming over surfaces respectively. Here we show that both swimming and surface-swarming of V. alginolyticus involve chemotaxis and are regulated by a single CheY species. Some of the substitutions of CheY residues conserved in various bacteria have different effects on the Pof and Laf motors, implying that CheY interacts with the two motors differently. Furthermore, analyses of tethered cells revealed that their switching modes are different: the Laf motor rotates exclusively counterclockwise and is slowed down by CheY, whereas the Pof motor turns both counterclockwise and clockwise, and CheY controls its rotational direction.