Crystal structure of Escherichia coli HdeA

Do All Paths Lead to Rome? How Reliable is Umbrella Sampling Along a Single Path?

Molecular dynamics (MD) simulations are widely applied to estimate absolute binding free energies of protein-ligand and protein-protein complexes. A routinely used method for binding free energy calculations with MD is umbrella sampling (US), which calculates the potential of mean force (PMF) along a single reaction coordinate. Surprisingly, in spite of its widespread use, few validation studies have focused on the convergence of the free energy computed along a single path for specific cases, not addressing the reproducibility of such calculations in general. In this work, we therefore investigate the reproducibility and convergence of US along a standard distance-based reaction coordinate for various protein-protein and protein-ligand complexes, following commonly used guidelines for the setup. We show that repeating the complete US workflow can lead to differences of 2-20 kcal/mol in computed binding free energies. We attribute those discrepancies to small differences in the binding pathways. While these differences are unavoidable in the established US protocol, the popularity of the latter could hint at a lack of awareness of such reproducibility problems. To test if the convergence of PMF profiles can be improved if multiple pathways are sampled simultaneously, we performed additional simulations with an adaptive-biasing method, here the accelerated weight histogram (AWH) approach. Indeed, the PMFs obtained from AHW simulations are consistent and reproducible for the systems tested. To the best of our knowledge, our work is the first to attempt a systematic assessment of the pitfalls in one the most widely used protocols for computing binding affinities. We anticipate therefore that our results will provide an incentive for a critical reassessment of the validity of PMFs computed with US, and make a strong case to further benchmark the performance of adaptive-biasing methods for computing binding affinities.

The Chaperone-Active State of HdeB at pH 4 Arises from Its Conformational Rearrangement and Enhanced Stability Instead of Surface Hydrophobicity

HdeA and HdeB are dimeric ATP-independent acid-stress chaperones, which protect the periplasmic proteins of enteric bacteria at pH 2.0 and 4.0, respectively, during their passage through the acidic environment of the mammalian stomach. Despite being structurally similar, they exhibit distinct functional pH optima and conformational prerequisite for their chaperone action. HdeA undergoes a dimer-to-monomer transition at pH 2.0, whereas HdeB remains dimeric at pH 4.0. The monomerization of HdeA exposes its hydrophobic motifs, which facilitates its interaction with the partially folded substrates. How HdeB functions despite maintaining its dimeric conformation has been poorly elucidated in the literature. Herein, we characterized the conformational states and stability of HdeB at its physiologically relevant pH and compared the data with those of HdeA. At pH 4.0, HdeB exhibited distinct spectroscopic signatures and higher stability against heat and guanidine-HCl-induced denaturation than at pH 7.5. We affirm that the pH 4.0 conformer of HdeB was distinct from that at pH 7.5 and that these two conformational states were hierarchically unrelated. Salt-bridge mutations that perturbed HdeB's intersubunit interactions resulted in the loss of its stability and function at pH 4.0. In contrast, mutations affecting intrasubunit interactions enhanced its function, albeit with a reduction in stability. These findings suggest that, unlike HdeA, HdeB acts as a noncanonical chaperone, where pH-dependent stability and conformational rearrangement at pH 4.0 play a core role in its chaperone function rather than its surface hydrophobicity. Such rearrangement establishes a stability-function trade-off that allows HdeB to function while maintaining its stable dimeric state.

Formation of Fibrils by the Periplasmic Molecular Chaperone HdeB from Escherichia coli

The molecular chaperones HdeA and HdeB of the Escherichia coli (E. coli) periplasm protect client proteins from acid denaturation through a unique mechanism that utilizes their acid denatured states to bind clients. We previously demonstrated that the active, acid-denatured form of HdeA is also prone to forming inactive, amyloid fibril-like aggregates in a pH-dependent, reversible manner. In this study, we report that HdeB also displays a similar tendency to form fibrils at low pH. HdeB fibrils were observed at pH < 3 in the presence of NaCl. Similar to HdeA, HdeB fibrils could be resolubilized by a simple shift to neutral pH. In the case of HdeB, however, we found that after extended incubation at low pH, HdeB fibrils were converted into a form that could not resolubilize at pH 7. Fresh fibrils seeded from these “transformed” fibrils were also incapable of resolubilizing at pH 7, suggesting that the transition from reversible to irreversible fibrils involved a specific conformational change that was transmissible through fibril seeds. Analyses of fibril secondary structure indicated that HdeB fibrils retained significant alpha helical content regardless of the conditions under which fibrils were formed. Fibrils that were formed from HdeB that had been treated to remove its intrinsic disulfide bond also were incapable of resolubilizing at pH 7, suggesting that certain residual structures that are retained in acid-denatured HdeB are important for this protein to recover its soluble state from the fibril form.

Stress-Responsive Periplasmic Chaperones in Bacteria

Periplasmic proteins are involved in a wide range of bacterial functions, including motility, biofilm formation, sensing environmental cues, and small-molecule transport. In addition, a wide range of outer membrane proteins and proteins that are secreted into the media must travel through the periplasm to reach their final destinations. Since the porous outer membrane allows for the free diffusion of small molecules, periplasmic proteins and those that travel through this compartment are more vulnerable to external environmental changes, including those that result in protein unfolding, than cytoplasmic proteins are. To enable bacterial survival under various stress conditions, a robust protein quality control system is required in the periplasm. In this review, we focus on several periplasmic chaperones that are stress responsive, including Spy, which responds to envelope-stress, DegP, which responds to temperature to modulate chaperone/protease activity, HdeA and HdeB, which respond to acid stress, and UgpB, which functions as a bile-responsive chaperone.

The complex role of the N-terminus and acidic residues of HdeA as pH-dependent switches in its chaperone function

HdeA is a small acid-stress chaperone protein found in the periplasm of several pathogenic gram-negative bacteria. In neutral pH environments HdeA is an inactive folded homodimer but when exposed to strong acidic environments it partially unfolds and, once activated, binds to other periplasmic proteins, protecting them from irreversible aggregation. Here we use a combination of hydrogen/deuterium exchange NMR experiments and constant pH molecular dynamics simulations to elucidate the role of HdeA's N-terminus in its activation mechanism. Previous work indicates that the N-terminus is flexible and unprotected at high pH while exhibiting interactions with some HdeA client binding site residues. It, however, becomes partially solvent-protected at pH 2.6 - 2.8 and then loses protection again at pH 2.0. This protection is not due to the appearance of new secondary structure, but rather increased contacts between N-terminal residues and the C-terminus of the other protomer in the dimer, as well as concurrent loosening of its hold on the client binding site residues, priming HdeA for interactions with periplasmic client proteins. This work also uncovers unusual protonation profiles of some titratable residues and suggests their complex role in chaperone function.

Biogenesis, quality control, and structural dynamics of proteins as explored in living cells via site-directed photocrosslinking

Protein biogenesis and quality control are essential to maintaining a functional pool of proteins and involve numerous protein factors that dynamically and transiently interact with each other and with the substrate proteins in living cells. Conventional methods are hardly effective for studying dynamic, transient, and weak protein-protein interactions that occur in cells. Herein, we review how the site-directed photocrosslinking approach, which relies on the genetic incorporation of a photoreactive unnatural amino acid into a protein of interest at selected individual amino acid residue positions and the covalent trapping of the interacting proteins upon ultraviolent irradiation, has become a highly efficient way to explore the aspects of protein contacts in living cells. For example, in the past decade, this approach has allowed the profiling of the in vivo substrate proteins of chaperones or proteases under both physiologically optimal and stressful (e.g., acidic) conditions, mapping residues located at protein interfaces, identifying new protein factors involved in the biogenesis of membrane proteins, trapping transiently formed protein complexes, and snapshotting different structural states of a protein. We anticipate that the site-directed photocrosslinking approach will play a fundamental role in dissecting the detailed mechanisms of protein biogenesis, quality control, and dynamics in the future.

Periplasmic Chaperones and Prolyl Isomerases

The biogenesis of periplasmic and outer membrane proteins (OMPs) in Escherichia coli is assisted by a variety of processes that help with their folding and transport to their final destination in the cellular envelope. Chaperones are macromolecules, usually proteins, that facilitate the folding of proteins or prevent their aggregation without becoming part of the protein's final structure. Because chaperones often bind to folding intermediates, they often (but not always) act to slow protein folding. Protein folding catalysts, on the other hand, act to accelerate specific steps in the protein folding pathway, including disulfide bond formation and peptidyl prolyl isomerization. This review is primarily concerned with E. coli and Salmonella periplasmic and cellular envelope chaperones; it also discusses periplasmic proline isomerization.

Structural basis and mechanism of the unfolding-induced activation of HdeA, a bacterial acid response chaperone

The role of protein structural disorder in biological functions has gained increasing attention in the past decade. The bacterial acid-resistant chaperone HdeA belongs to a group of "conditionally disordered" proteins, because it is inactive in its well-structured state and becomes activated via an order-to-disorder transition under acid stress. However, the mechanism for unfolding-induced activation remains unclear because of a lack of experimental information on the unfolded state conformation and the chaperone-client interactions. Herein, we used advanced solution NMR methods to characterize the activated-state conformation of HdeA under acidic conditions and identify its client-binding sites. We observed that the structure of activated HdeA becomes largely disordered and exposes two hydrophobic patches essential for client interactions. Furthermore, using the pH-dependent chemical exchange saturation transfer (CEST) NMR method, we identified three acid-sensitive regions that act as structural locks in regulating the exposure of the two client-binding sites during the activation process, revealing a multistep activation mechanism of HdeA's chaperone function at the atomic level. Our results highlight the role of intrinsic protein disorder in chaperone function and the self-inhibitory role of ordered structures under nonstress conditions, offering new insights for improving our understanding of protein structure-function paradigms.

Acid-denatured small heat shock protein HdeA from Escherichia coli forms reversible fibrils with an atypical secondary structure

The periplasmic small heat shock protein HdeA from Escherichia coli is inactive under normal growth conditions (at pH 7) and activated only when E. coli cells are subjected to a sudden decrease in pH, converting HdeA into an acid-denatured active state. Here, using in vitro fibrillation assays, transmission EM, atomic-force microscopy, and CD analyses, we found that when HdeA is active as a molecular chaperone, it is also capable of forming inactive aggregates that, at first glance, resemble amyloid fibrils. We noted that the molecular chaperone activity of HdeA takes precedence over fibrillogenesis under acidic conditions, as the presence of denatured substrate protein was sufficient to suppress HdeA fibril formation. Further experiments suggested that the secondary structure of HdeA fibrils deviates somewhat from typical amyloid fibrils and contains α-helices. Strikingly, HdeA fibrils that formed at pH 2 were immediately resolubilized by a simple shift to pH 7 and from there could regain molecular chaperone activity upon a return to pH 1. HdeA, therefore, provides an unusual example of a "reversible" form of protein fibrillation with an atypical secondary structure composition. The competition between active assistance of denatured polypeptides (its "molecular chaperone" activity) and the formation of inactive fibrillary deposits (its "fibrillogenic" activity) provides a unique opportunity to probe the relationship among protein function, structure, and aggregation in detail.

Electrostatic Interactions in Protein Structure, Folding, Binding, and Condensation

Charged and polar groups, through forming ion pairs, hydrogen bonds, and other less specific electrostatic interactions, impart important properties to proteins. Modulation of the charges on the amino acids, e.g., by pH and by phosphorylation and dephosphorylation, have significant effects such as protein denaturation and switch-like response of signal transduction networks. This review aims to present a unifying theme among the various effects of protein charges and polar groups. Simple models will be used to illustrate basic ideas about electrostatic interactions in proteins, and these ideas in turn will be used to elucidate the roles of electrostatic interactions in protein structure, folding, binding, condensation, and related biological functions. In particular, we will examine how charged side chains are spatially distributed in various types of proteins and how electrostatic interactions affect thermodynamic and kinetic properties of proteins. Our hope is to capture both important historical developments and recent experimental and theoretical advances in quantifying electrostatic contributions of proteins.

The Mechanism of HdeA Unfolding and Chaperone Activation

HdeA is a periplasmic chaperone that is rapidly activated upon shifting the pH to acidic conditions. This activation is thought to involve monomerization of HdeA. There is evidence that monomerization and partial unfolding allow the chaperone to bind to proteins denatured by low pH, thereby protecting them from aggregation. We analyzed the acid-induced unfolding of HdeA using NMR spectroscopy and fluorescence measurements, and obtained experimental evidence suggesting a complex mechanism in HdeA's acid-induced unfolding pathway, as previously postulated from molecular dynamics simulations. Counterintuitively, dissociation constant measurements show a stabilization of the HdeA dimer upon exposure to mildly acidic conditions. We provide experimental evidence that protonation of Glu37, a glutamate residue embedded in a hydrophobic pocket of HdeA, is important in controlling HdeA stabilization and thus the acid activation of this chaperone. Our data also reveal a sharp transition from folded dimer to unfolded monomer between pH3 and pH 2, and suggest the existence of a low-populated, partially folded intermediate that could assist in chaperone activation or function. Overall, this study provides a detailed experimental investigation into the mechanism by which HdeA unfolds and activates.

Characterizations of the Interactions between Escherichia coli Periplasmic Chaperone HdeA and Its Native Substrates during Acid Stress

The bacterial acid-resistant chaperone HdeA is a "conditionally disordered" protein that functions at low pH when it undergoes a transition from a well-folded dimer to an unfolded monomer. The dimer dissociation and unfolding processes result in exposure of hydrophobic surfaces that allows binding to a broad range of client proteins. To fully elucidate the chaperone mechanism of HdeA, it is crucial to understand how the activated HdeA interacts with its native substrates during acid stress. Herein, we present a nuclear magnetic resonance study of the pH-dependent HdeA-substrate interactions. Our results show that the activation of HdeA is not only induced by acidification but also regulated by the presence of unfolded substrates. The variable extent of unfolding of substrates differentially regulates the HdeA-substrate interaction, and the binding further affects the HdeA conformation. Finally, we show that HdeA binds its substrates heterogeneously, and the "amphiphilic" model for HdeA-substrate interaction is discussed.

Probing the Structure of the Escherichia coli Periplasmic Proteins HdeA and YmgD by Molecular Dynamics Simulations

HdeA and YmgD are structurally homologous proteins in the periplasm of Escherichia coli. HdeA has been shown to represent an acid-activated chaperone, whereas the function of YmgD has not yet been characterized. We performed pH-titrating molecular dynamics simulations (pHtMD) to investigate the structural changes of both proteins and to assess whether YmgD may also exhibit an unfolding behavior similar to that of HdeA. The unfolding pathway of HdeA includes partially unfolded dimer structures, which represent a prerequisite for subsequent dissociation. In contrast to the coupled unfolding and dissociation of HdeA, YmgD displays dissociation of the folded subunits, and the subunits do not undergo significant unfolding even at low pH values. The differences in subunit stability between HdeA and YmgD may be explained by the structural features of helix D, which represents the starting point of unfolding in HdeA. In summary, the present study suggests that YmgD either is not an acid-activated chaperone or, at least, does not require unfolding for activation.

Comparative proteomics reveal distinct chaperone-client interactions in supporting bacterial acid resistance

HdeA and HdeB constitute the essential chaperone system that functions in the unique periplasmic space of Gram-negative enteric bacteria to confer acid resistance. How this two-chaperone machinery cooperates to protect a broad range of client proteins from acid denaturation while avoiding nonspecific binding during bacterial passage through the highly acidic human stomach remains unclear. We have developed a comparative proteomic strategy that combines the genetically encoded releasable protein photocross-linker with 2D difference gel electrophoresis, which allows an unbiased side-by-side comparison of the entire client pools from these two acid-activated chaperones in Escherichia coli Our results reveal distinct client specificities between HdeA and HdeB in vivo that are determined mainly by their different responses to pH stimulus. The intracellular acidity serves as an environmental cue to determine the folding status of both chaperones and their clients, enabling specific chaperone-client binding and release under defined pH conditions. This cooperative and synergistic mode of action provides an efficient, economical, flexible, and finely tuned protein quality control strategy for coping with acid stress.

Coupled folding and binding with 2D Window-Exchange Umbrella Sampling

Intrinsically disordered regions of proteins can gain structure by binding to a partner. This process, of coupled folding and binding (CFaB), is a fundamental part of many important biological processes. Structure-based models have proven themselves capable of revealing fundamental aspects of how CFaB occurs, however, typical methods to enhance the sampling of these transitions, such as replica exchange, do not adequately sample the transition state region of this extremely rare process. Here, we use a variant of Umbrella Sampling to enforce sampling of the transition states of CFaB of HdeA monomers at neutral pH, an extremely rare process that occurs over timescales ranging from seconds to hours. Using high resolution sampling in the transition state region, we cluster states along the principal transition path to obtain a detailed description of coupled binding and folding for the HdeA dimer, revealing new insight into the ensemble of states that are accessible to client recognition. We then demonstrate that exchanges between umbrella sampling windows, as done in previous work, significantly improve relaxation in variables orthogonal to the restraints used. Altogether, these results show that Window-Exchange Umbrella Sampling is a promising approach for systems that exhibit flexible binding, and can reveal transition state ensembles of these systems in high detail. © 2015 Wiley Periodicals, Inc.

Mimicking titration experiments with MD simulations: A protocol for the investigation of pH-dependent effects on proteins

Protein structure and function are highly dependent on the environmental pH. However, the temporal or spatial resolution of experimental approaches hampers direct observation of pH-induced conformational changes at the atomic level. Molecular dynamics (MD) simulation strategies (e.g. constant pH MD) have been developed to bridge this gap. However, one frequent problem is the sampling of unrealistic conformations, which may also lead to poor pKa predictions. To address this problem, we have developed and benchmarked the pH-titration MD (pHtMD) approach, which is inspired by wet-lab titration experiments. We give several examples how the pHtMD protocol can be applied for pKa calculation including peptide systems, Staphylococcus nuclease (SNase), and the chaperone HdeA. For HdeA, pHtMD is also capable of monitoring pH-dependent dimer dissociation in accordance with experiments. We conclude that pHtMD represents a versatile tool for pKa value calculation and simulation of pH-dependent effects in proteins.

Roles of structural plasticity in chaperone HdeA activity are revealed by 19F NMR

Multiple conformations of acid chaperone HdeA and their roles in activity.

HdeA, a minimal ATP-independent acid chaperone, is crucial for the survival of enteric pathogens as they transit the acidic (pH 1–3) environment of the stomach. Although protein disorder (unfolding) and structural plasticity have been elegantly linked to HdeA function, the details of the linkage are lacking. Here, we apply ¹⁹F NMR to reveal the structural transition associated with activation. We find that unfolding is necessary but not sufficient for activation. Multiple conformations are present in the functional state at low pH, but the partially folded conformation is essential for HdeA chaperone activity, and HdeA's intrinsic disulfide bond is required to maintain the partially folded conformation. The results show that both disorder and order are key to function. The ability of ¹⁹F NMR to reveal and quantify multiple conformational states makes it a powerful tool for studying other chaperones.

HdeB chaperone activity is coupled to its intrinsic dynamic properties

Enteric bacteria encounter extreme acidity when passing through hosts' stomach. Since the bacterial periplasmic space quickly equilibrates with outer environment, an efficient acid resistance mechanism is essential in preventing irreversible protein denaturation/aggregation and maintaining bacteria viability. HdeB, along with its homolog HdeA, was identified as a periplasmic acid-resistant chaperone. Both proteins exist as homodimers and share similar monomeric structures under neutral pH, while showing different dimeric packing interfaces. Previous investigations show that HdeA functions through an acid-induced dimer-to-monomer transition and partial unfolding at low pH (pH 2-3), resulting in exposure of hydrophobic surfaces that bind substrate proteins. In contrast, HdeB appears to have a much higher optimal activation pH (pH 4-5), under which condition the protein maintains a well-folded dimer and the mechanism for its chaperone activity remains elusive. Herein, we present an NMR study of HdeB to investigate its dynamic properties. Our results reveal that HdeB undergoes significant micro- to milli-second timescale conformational exchanges at neutral to near-neutral pH, under the later condition it exhibits optimal activity. The current study indicates that HdeB activation is coupled to its intrinsic dynamics instead of structural changes, and therefore its functional mechanism is apparently different from HdeA.

Multiscale modeling of a conditionally disordered pH-sensing chaperone

The pH-sensing chaperone HdeA promotes the survival of enteropathogenic bacteria during transit through the harshly acidic environment of the mammalian stomach. At low pH, HdeA transitions from an inactive, folded, dimer to chaperone-active, disordered, monomers to protect against the acid-induced aggregation of periplasmic proteins. Toward achieving a detailed mechanistic understanding of the pH response of HdeA, we develop a multiscale modeling approach to capture its pH-dependent thermodynamics. Our approach combines pK(a) (logarithmic acid dissociation constant) calculations from all-atom constant pH molecular dynamics simulations with coarse-grained modeling and yields new, atomic-level, insights into HdeA chaperone function that can be directly tested by experiment. "pH triggers" that significantly destabilize the dimer are each located near the N-terminus of a helix, suggesting that their neutralization at low pH destabilizes the helix macrodipole as a mechanism of monomer disordering. Moreover, we observe a non-monotonic change in the pH-dependent stability of HdeA, with maximal stability of the dimer near pH5. This affect is attributed to the protonation Glu37, which exhibits an anomalously high pK(a) value and is located within the hydrophobic dimer interface. Finally, the pH-dependent binding pathway of HdeA comprises a partially unfolded, dimeric intermediate that becomes increasingly stable relative to the native dimer at lower pH values and displays key structural features for chaperone-substrate interaction. We anticipate that the insights from our model will help inform ongoing NMR and biochemical investigations.

The Escherichia coli Acid Stress Response and Its Significance for Pathogenesis

Escherichia coli has a remarkable ability to survive low pH and possesses a number of different genetic systems that enable it to do this. These may be expressed constitutively, typically in stationary phase, or induced by growth under a variety of conditions. The activities of these systems have been implicated in the ability of E. coli to pass the acidic barrier of the stomach and to become established in the gastrointestinal tract, something causing serious infections. However, much of the work characterizing these systems has been done on standard laboratory strains of E. coli and under conditions which do not closely resemble those found in the human gut. Here we review what is known about acid resistance in E. coli as a model laboratory organism and in the context of its lifestyle as an inhabitant-sometimes an unwelcome one-of the human gut.

Folding mechanisms of periplasmic proteins

More than one fifth of the proteins encoded by the genome of Escherichia coli are destined to the bacterial cell envelope. Over the past 20years, the mechanisms by which envelope proteins reach their three-dimensional structure have been intensively studied, leading to the discovery of an intricate network of periplasmic folding helpers whose members have distinct but complementary roles. For instance, the correct assembly of ß-barrel proteins containing disulfide bonds depends both on chaperones like SurA and Skp for transport across the periplasm and on protein folding catalysts like DsbA and DsbC for disulfide bond formation. In this review, we provide an overview of the current knowledge about the complex network of protein folding helpers present in the periplasm of E. coli and highlight the questions that remain unsolved. This article is part of a Special Issue entitled: Protein trafficking and secretion in bacteria. Guest Editors: Anastassios Economou and Ross Dalbey.

Binding and folding of the small bacterial chaperone HdeA

The small pH stress-sensing chaperone HdeA helps pathogenic enteric E. coli survive passage through the severely acidic environment of the mammalian stomach. Under stress conditions, HdeA transitions from an inactive folded dimer to a chaperone-active unfolded monomer to prevent the acid-induced aggregation of periplasmic proteins. Here we use a topology-based Gō-like model to delineate the relationship between dimer interface formation and monomer folding and to better understand the structural details of the chaperone activation mechanism. Free energy surfaces show that dimer interface formation and monomer folding proceed concurrently through an on-pathway dimeric intermediate in which one monomer is partially unfolded. The absence of a preexisting fully folded monomer in the proposed binding mechanism is in agreement with HdeA's rapid chaperone response. Binding between unfolded monomers exhibits an enhancement of molecular recognition reminiscent of the fly-casting mechanism. Overall, our simulations further highlight the efficient nature of HdeA's chaperone response and we anticipate that knowledge of a dimeric intermediate will facilitate the interpretation of experimental studies.

Chaperone activation by unfolding

Conditionally disordered proteins can alternate between highly ordered and less ordered configurations under physiological conditions. Whereas protein function is often associated with the ordered conformation, for some of these conditionally unstructured proteins, the opposite applies: Their activation is associated with their unfolding. An example is the small periplasmic chaperone HdeA, which is critical for the ability of enteric bacterial pathogens like Escherichia coli to survive passage through extremely acidic environments, such as the human stomach. At neutral pH, HdeA is a chaperone-inactive dimer. On a shift to low pH, however, HdeA monomerizes, partially unfolds, and becomes rapidly active in preventing the aggregation of substrate proteins. By mutating two aspartic acid residues predicted to be responsible for the pH-dependent monomerization of HdeA, we have succeeded in isolating an HdeA mutant that is active at neutral pH. We find this HdeA mutant to be substantially destabilized, partially unfolded, and mainly monomeric at near-neutral pH at a concentration at which it prevents aggregation of a substrate protein. These results provide convincing evidence for direct activation of a protein by partial unfolding.

Probing pH-dependent dissociation of HdeA dimers

HdeA protein is a small, ATP-independent, acid stress chaperone that undergoes a dimer-to-monomer transition in acidic environments. The HdeA monomer binds a broad range of proteins to prevent their acid-induced aggregation. To understand better HdeA's function and mechanism, we perform constant-pH molecular dynamics simulations (CPHMD) to elucidate the details of the HdeA dimer dissociation process. First the pK(a) values of all the acidic titratable groups in HdeA are obtained and reveal a large pK(a) shift only for Glu(37). However, the pH-dependent monomer charge exhibits a large shift from -4 at pH > 6 to +6 at pH = 2.5, suggesting that the dramatic change in charge on each monomer may drive dissociation. By combining the CPHMD approach with umbrella sampling, we demonstrate a significant stability decrease of the HdeA dimer when the environmental pH changes from 4.0 to 3.5 and identify the key acidic residue-lysine interactions responsible for the observed pH sensing in HdeA chaperon activity function.

Salt bridges regulate both dimer formation and monomeric flexibility in HdeB and may have a role in periplasmic chaperone function

Escherichia coli and Gram-negative bacteria that live in the human gut must be able to tolerate rapid and large changes in environmental pH. Low pH irreversibly denatures and precipitates many bacterial proteins. While cytoplasmic proteins are well buffered against such swings, periplasmic proteins are not. Instead, it appears that some bacteria utilize chaperone proteins that stabilize periplasmic proteins, preventing their precipitation. Two highly expressed and related proteins, HdeA and HdeB, have been identified as acid-activated chaperones. The structure of HdeA is known and a mechanism for activation has been proposed. In this model, dimeric HdeA dissociates at low pH, and the exposed dimeric interface binds exposed hydrophobic surfaces of acid-denatured proteins, preventing their irreversible aggregation. We now report the structure and biophysical characterization of the HdeB protein. The monomer of HdeB shares a similar structure with HdeA, but its dimeric interface is different in composition and spatial location. We have used fluorescence to study the behavior of HdeB as pH is lowered, and like HdeA, it dissociates to monomers. We have identified one of the key intersubunit interactions that controls pH-induced monomerization. Our analysis identifies a structural interaction within the HdeB monomer that is disrupted as pH is lowered, leading to enhanced structural flexibility.

Structural biology of periplasmic chaperones

Proteins often require specific helper proteins, chaperones, to assist with their correct folding and to protect them from denaturation and aggregation. The cell envelope of Gram-negative bacteria provides a particularly challenging environment for chaperones to function in as it lacks readily available energy sources such as adenosine 5' triphosphate (ATP) to power reaction cycles. Periplasmic chaperones have therefore evolved specialized mechanisms to carry out their functions without the input of external energy and in many cases to transduce energy provided by protein folding or ATP hydrolysis at the inner membrane. Structural and biochemical studies have in recent years begun to elucidate the specific functions of many important periplasmic chaperones and how these functions are carried out. This includes not only specific carrier chaperones, such as those involved in the biosynthesis of adhesive fimbriae in pathogenic bacteria, but also more general pathways including the periplasmic transport of outer membrane proteins and the extracytoplasmic stress responses. This chapter aims to provide an overview of protein chaperones so far identified in the periplasm and how structural biology has assisted with the elucidation of their functions.

Role of the multidrug resistance regulator MarA in global regulation of the hdeAB acid resistance operon in Escherichia coli

MarA, a transcriptional regulator in Escherichia coli, affects functions such as multiple-antibiotic resistance (Mar) and virulence. Usually an activator, MarA is a repressor of hdeAB and other acid resistance genes. We found that, in wild-type cells grown in LB medium at pH 7.0 or pH 5.5, repression of hdeAB by MarA occurred only in stationary phase and was reduced in the absence of H-NS and GadE, the main regulators of hdeAB. Moreover, repression of hdeAB by MarA was greater in the absence of GadX or Lrp in exponential phase at pH 7.0 and in the absence of GadW or RpoS in stationary phase at pH 5.5. In turn, MarA enhanced repression of hdeAB by H-NS and hindered activation by GadE in stationary phase and also reduced the activity of GadX, GadW, RpoS, and Lrp on hdeAB under some conditions. As a result of its direct and indirect effects, overexpression of MarA prevented most of the induction of hdeAB expression as cells entered stationary phase and made the cells sevenfold more sensitive to acid challenge at pH 2.5. These findings show that repression of hdeAB by MarA depends on pH, growth phase, and other regulators of hdeAB and is associated with reduced resistance to acid conditions.

Periplasmic Protein HdeA Exhibits Chaperone-like Activity Exclusively within Stomach pH Range by Transforming into Disordered Conformation

The extremely acidic environment of the mammalian stomach, with a pH range usually between 1 and 3, represents a stressful challenge for enteric pathogenic bacteria such as Escherichia coli before they enter into the intestine. The hdeA gene of E. coli was found to be acid inducible and was revealed by genetic studies to be important for the acid survival of the strain. This study was performed in an attempt to characterize the mechanism of the activity of the HdeA protein. Our data provided in this report strongly suggest that HdeA employs a novel strategy to modulate its chaperone activity: it possesses an ordered conformation that is unable to bind denatured substrate proteins under normal physiological conditions (i.e. at neutral pH) and transforms into a globally disordered conformation that is able to bind substrate proteins under stress conditions (i.e. at a pH below 3). Furthermore, our data indicate that HdeA exposes hydrophobic surfaces that appear to be involved in the binding of denatured substrate proteins at extremely low pH values. In light of our observations, models are proposed to explain the action of HdeA in both a physiological and a molecular context.

The first crystal structure of an archaeal helical repeat protein

The crystal structure of ST1625p, a protein encoded by a hypothetical open reading frame ST1625 in the genome of the hyperthermophilic archaeon Sulfolobus tokodaii, was determined at 2.2 A resolution. The only sequence similarity exhibited by the amino-acid sequence of ST1625p was a 33% identity with the sequence of SSO0983p from S. solfataricus. The 19 kDa monomeric protein was observed to consist of a right-handed superhelix assembled from a tandem repeat of ten alpha-helices. A structural homology search using the DALI and MATRAS algorithms indicates that this protein can be classified as a helical repeat protein.

Water in protein structure prediction

Proteins have evolved to use water to help guide folding. A physically motivated, nonpairwise-additive model of water-mediated interactions added to a protein structure prediction Hamiltonian yields marked improvement in the quality of structure prediction for larger proteins. Free energy profile analysis suggests that long-range water-mediated potentials guide folding and smooth the underlying folding funnel. Analyzing simulation trajectories gives direct evidence that water-mediated interactions facilitate native-like packing of supersecondary structural elements. Long-range pairing of hydrophilic groups is an integral part of protein architecture. Specific water-mediated interactions are a universal feature of biomolecular recognition landscapes in both folding and binding.

Crystal structure of a conserved hypothetical protein from Escherichia coli

The crystal structure of a conserved hypothetical protein from Escherichia coli has been determined using X-ray crystallography. The protein belongs to the Cluster of Orthologous Group COG1553 (National Center for Biotechnology Information database, NLM, NIH), for which there was no structural information available until now. Structural homology search with DALI algorism indicated that this protein has a new fold with no obvious similarity to those of other proteins with known three-dimensional structures. The protein quaternary structure consists of a dimer of trimers, which makes a characteristic cylinder shape. There is a large closed cavity with approximate dimensions of 16 A x 16 A x 20 A in the center of the hexameric structure. Six putative active sites are positioned along the equatorial surface of the hexamer. There are several highly conserved residues including two possible functional cysteines in the putative active site. The possible molecular function of the protein is discussed.

Annotation in three dimensions. PINTS: Patterns in Non-homologous Tertiary Structures

The detection of local structural patterns in proteins (e.g. active sites) can provide insights into protein function in the absence of sequence or fold similarity. Methods to detect such similarities are key during structural annotation, for example with results from Structural Genomics initiatives. PINTS (Patterns in Non-homologous Tertiary Structures, http://pints.embl.de) performs database searches for such patterns and most importantly provides a measure of statistical significance for any similarity uncovered. To aid functional annotation of proteins, we allow comparisons of pre-defined patterns against databases of complete structures and of entire structures to databases of particular residues likely to be functionally important.

Structural proteomics: lessons learnt from the early case studies

The genomics efforts have identified a large number of novel genes and thus provided a pool of interesting but not functionally characterized target proteins. It has been suggested that structural proteomics will significantly impact the success rate of functional characterization of such identified genes and proteins by providing structure-function hypotheses by fold and feature recognition and analysis. Structural proteomics initiatives, both in academic and industrial settings, are today generating protein structures at an unprecedented rate although relatively few large-scale efforts have been displayed in the public domain. However, a number of individual studies have provided a 'road-map' for selected approaches that hold the promise to significantly impact the process of deriving function from structure.

Crystal structure of conserved hypothetical protein Aq1575 from Aquifex aeolicus

The crystal structure of a conserved hypothetical protein, Aq1575, from Aquifex aeolicus has been determined by using x-ray crystallography. The protein belongs to the domain of unknown function DUF28 in the Pfam and PALI databases for which there was no structural information available until now. A structural homology search with the DALI algorithm indicates that this protein has a new fold with no obvious similarity to those of other proteins of known three-dimensional structure. The protein reveals a monomer consisting of three domains arranged along a pseudo threefold symmetry axis. There is a large cleft with approximate dimensions of 10 A x 10 A x 20 A in the center of the three domains along the symmetry axis. Two possible active sites are suggested based on the structure and multiple sequence alignment. There are several highly conserved residues in these putative active sites. The structure based molecular properties and thermostability of the protein are discussed.

Structural similarity to link sequence space: new potential superfamilies and implications for structural genomics

The current pace of structural biology now means that protein three-dimensional structure can be known before protein function, making methods for assigning homology via structure comparison of growing importance. Previous research has suggested that sequence similarity after structure-based alignment is one of the best discriminators of homology and often functional similarity. Here, we exploit this observation, together with a merger of protein structure and sequence databases, to predict distant homologous relationships. We use the Structural Classification of Proteins (SCOP) database to link sequence alignments from the SMART and Pfam databases. We thus provide new alignments that could not be constructed easily in the absence of known three-dimensional structures. We then extend the method of Murzin (1993b) to assign statistical significance to sequence identities found after structural alignment and thus suggest the best link between diverse sequence families. We find that several distantly related protein sequence families can be linked with confidence, showing the approach to be a means for inferring homologous relationships and thus possible functions when proteins are of known structure but of unknown function. The analysis also finds several new potential superfamilies, where inspection of the associated alignments and superimpositions reveals conservation of unusual structural features or co-location of conserved amino acids and bound substrates. We discuss implications for Structural Genomics initiatives and for improvements to sequence comparison methods.

Structural proteomics: developments in structure-to-function predictions

The major challenge for post-genomic research is to functionally assign and validate a large number of novel target genes and their corresponding proteins. Functional genomics approaches have, therefore, gained considerable attention in the quest to convert this massive data set into useful information. One of the crucial components for the functional understanding of unassigned proteins is the analysis of their experimental or modeled 3D structures. Structural proteomics initiatives are generating protein structures at an unprecedented rate but our current knowledge of 3D-structural space is still limited. Estimates on the completeness of the 3D-structural coverage of proteins vary but it is generally accepted that only a minority of the structural proteome has a template structure from which reliable conclusions can be drawn. Thus, structural proteomics has set out to build a map of protein structures that will represent all protein folds included in the 'global proteome'.

Predicting function from structure: examples of the serine protease inhibitor canonical loop conformation found in extracellular proteins

The prediction of protein function from structure is becoming of growing importance in the age of structural genomics. We have focused on the problem of identifying sites of potential serine protease inhibitor interactions on the surface of proteins of known structure. Given that there is no sequence conservation within canonical loops from different inhibitor families we first compare representative loops to all fragments of equal length among proteins of known structure by calculating main-chain RMS deviation. Fragments with RMS deviation below a certain threshold (hits) are removed if residues have solvent accessibilities appreciably lower than those observed in the search structure. These remaining hits are further filtered to remove those occurring largely within secondary structure elements. Likely functional significance is restricted further by considering only extracellular protein domains. Also a test is performed to see if the loop can dock into the binding site of the serine protease trypsin without unacceptable steric clashes. By comparing different canonical loop structures to the protein structure database we show that the method was able to detect previously known inhibitors. In addition, we discuss potentially new canonical loop structures found in secreted hydrolases, toxins, viral proteins, cytokines and other proteins. We discuss the possible functional significance of several of the examples found.

Solution structure and function of a conserved protein SP14.3 encoded by an essential Streptococcus pneumoniae gene

Streptococcus pneumoniae is a major human pathogen that causes high mortality and morbidity rates and has developed resistance to many antibiotics. The genome of S. pneumoniae has recently been completely sequenced revealing many genes encoding hypothetical proteins of unknown function. We have found that the gene encoding one such conserved protein, SP14.3, is essential for growth of S. pneumonia. Since it is essential, SP14.3 represents a potential target for drug discovery. Here, we describe the three-dimensional solution structure of SP14.3 as determined by NMR spectroscopy. The structure consists of two domains each with an alpha/beta-fold. The N-terminal domain contains two alpha-helices and a three-stranded beta-sheet, while the C-terminal domain is composed of one alpha-helix and a five-stranded beta-sheet. The N-terminal domain of the protein contains a highly negatively charged surface and resembles the fold of the N-terminal domain of Thermus thermophilus ribosomal protein S3. The C-terminal domain has a protein fold similar to human small nuclear ribonucleoprotein Sm D3 and Haloarcula marismortui ribosomal protein L21E. The two domains of the protein tumble in solution overall as a whole with an overall molecular rotational correlation time (tau(m)) of 12.9 ns at 25 degrees C. The relative orientation of the two domains is not defined by the nuclear Overhauser effect data. Indeed, residual dipolar couplings and the structure calculations indicate that the relative orientation of the two domains is not rigidly oriented with respect to one another in solution.

Determination of protein function, evolution and interactions by structural genomics

The genome sequencing projects and knowledge of the entire protein repertoires of many organisms have prompted new procedures and techniques for the large-scale determination of protein structure, function and interactions. Recently, new work has been carried out on the determination of the function and evolutionary relationships of proteins by experimental structural genomics, and the discovery of protein-protein interactions by computational structural genomics.

Fold predictions for bacterial genomes

Fold assignments for newly sequenced genomes belong to the most important and interesting applications of the booming field of protein structure prediction. We present a brief survey and a discussion of such assignments completed to date, using as an example several fold assignment projects for proteins from the Escherichia coli genome. This review focuses on steps that are necessary to go beyond the simple assignment projects and into the development of tools extending our understanding of functions of proteins in newly sequenced genomes. This paper also discusses several problems seldom addressed in the literature, such as the problem of domain prediction and complementary predictions (e.g., transmembrane regions and flexible regions) and cross-correlation of predictions from different servers. The influence of sequence and structure database growth on prediction success is also addressed. Finally, we discuss the perspectives of the field in the context of massive sequence and structure determination projects, as well as the development of novel prediction methods.

Evolution of function in protein superfamilies, from a structural perspective

The recent growth in protein databases has revealed the functional diversity of many protein superfamilies. We have assessed the functional variation of homologous enzyme superfamilies containing two or more enzymes, as defined by the CATH protein structure classification, by way of the Enzyme Commission (EC) scheme. Combining sequence and structure information to identify relatives, the majority of superfamilies display variation in enzyme function, with 25 % of superfamilies in the PDB having members of different enzyme types. We determined the extent of functional similarity at different levels of sequence identity for 486,000 homologous pairs (enzyme/enzyme and enzyme/non-enzyme), with structural and sequence relatives included. For single and multi-domain proteins, variation in EC number is rare above 40 % sequence identity, and above 30 %, the first three digits may be predicted with an accuracy of at least 90 %. For more distantly related proteins sharing less than 30 % sequence identity, functional variation is significant, and below this threshold, structural data are essential for understanding the molecular basis of observed functional differences. To explore the mechanisms for generating functional diversity during evolution, we have studied in detail 31 diverse structural enzyme superfamilies for which structural data are available. A large number of variations and peculiarities are observed, at the atomic level through to gross structural rearrangements. Almost all superfamilies exhibit functional diversity generated by local sequence variation and domain shuffling. Commonly, substrate specificity is diverse across a superfamily, whilst the reaction chemistry is maintained. In many superfamilies, the position of catalytic residues may vary despite playing equivalent functional roles in related proteins. The implications of functional diversity within supefamilies for the structural genomics projects are discussed. More detailed information on these superfamilies is available at http://www.biochem.ucl.ac.uk/bsm/FAM-EC/.

The structure of the yrdC gene product from Escherichia coli reveals a new fold and suggests a role in RNA binding

The yrdC family of genes codes for proteins that occur both independently and as a domain in proteins that have been implicated in regulation. An example for the latter case is the sua5 gene from yeast. SuaS was identified as a suppressor of a translation initiation defect in cytochrome c and is required for normal growth in yeast (Na JG, Pinto I, Hampsey M, 1992, Genetics 11:791-801). However, the function of the Sua5 protein remains unknown; Sua5 could act either at the transcriptional or the posttranscriptional levels to compensate for an aberrant translation start codon in the cyc gene. To potentially learn more about the function of YrdC and proteins featuring this domain, the crystal structure of the YrdC protein from Escherichia coli was determined at a resolution of 2.0 A. YrdC adopts a new fold with no obvious similarity to those of other proteins with known three-dimensional (3D) structure. The protein features a large concave surface on one side that exhibits a positive electrostatic potential. The dimensions of this depression, its curvature, and the fact that conserved basic amino acids are located at its floor suggest that YrdC may be a nucleic acid binding protein. An investigation of YrdC's binding affinities for single- and double-stranded RNA and DNA fragments as well as tRNAs demonstrates that YrdC binds preferentially to double-stranded RNA. Our work provides evidence that 3D structures of functionally uncharacterized gene products with unique sequences can yield novel folds and functional insights.

Structure-based functional classification of hypothetical protein MTH538 from Methanobacterium thermoautotrophicum

The structure of MTH538, a previously uncharacterized hypothetical protein from Methanobacterium thermoautotrophicum, has been determined by NMR spectroscopy. MTH538 is one of numerous structural genomics targets selected in a genome-wide survey of uncharacterized sequences from this organism. MTH538 is a so-called singleton, a sequence not closely related to any other (known) sequences. The structure of MTH538 closely resembles the known structures of receiver domains from two component response regulator systems, such as CheY, and is similar to the structures of flavodoxins and GTP-binding proteins. Tests on MTH538 for characteristic activities of CheY and flavodoxin were negative. MTH538 did not become phosphorylated in the presence of acetyl phosphate and Mg(2+), although it appeared to bind Mg(2+). MTH538 also did not bind flavin mononucleotide (FMN) or coenzyme F(420). Nevertheless, sequence and structure parallels between MTH538/CheY and two families of ATPase/phosphatase proteins suggest that MTH538 may have a role in a phosphorylation-independent two-component response regulator system.