Probing pH-dependent dissociation of HdeA dimers

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.

Approach to Study pH-Dependent Protein Association Using Constant-pH Molecular Dynamics: Application to the Dimerization of β-Lactoglobulin

Protein-protein association is often mediated by electrostatic interactions and modulated by pH. However, experimental and computational studies have often overlooked the effect of association on the protonation state of the protein. In this work, we present a methodological approach based on constant-pH molecular dynamics (MD), which aims to provide a detailed description of a pH-dependent protein-protein association, and apply it to the dimerization of β-lactoglobulin (BLG). A selection of analyses is performed using the data generated by constant-pH MD simulations of monomeric and dimeric forms of bovine BLG, in the pH range 3-8. First, we estimate free energies of dimerization using a computationally inexpensive approach based on the Wyman-Tanford linkage theory, calculated in a new way through the use of thermodynamically based splines. The individual free energy contribution of each titratable site is also calculated, allowing for identification of relevant residues. Second, the correlations between the proton occupancies of pairs of sites are calculated (using the Pearson coefficient), and extensive networks of correlated sites are observed at acidic pH values, sometimes involving distant pairs. In general, strongly correlated sites are also slow proton exchangers and contribute significantly to the pH-dependency of the dimerization free energy. Third, we use ionic density as a fingerprint of protein charge distribution and observe electrostatic complementarity between the monomer faces that form the dimer interface, more markedly at the isoionic point (where maximum dimerization occurs) than at other pH values, which might contribute to guide the association. Finally, the pH-dependent dimerization modes are inspected using PCA, among other analyses, and two states are identified: a relaxed state at pH 4-8 (with the typical alignment of the crystallographic structure) and a compact state at pH 3-4 (with a tighter association and rotated alignment). This work shows that an approach based on constant-pH MD simulations can produce rich detailed pictures of pH-dependent protein associations, as illustrated for BLG dimerization.

Detection of key sites of dimer dissociation and unfolding initiation during activation of acid-stress chaperone HdeA at low pH

HdeA is a small acid-stress chaperone protein with a unique activity profile. At physiological pH, it forms a folded, but inactive, dimer. Below pH 3.0, HdeA unfolds and dissociates into disordered monomers, utilizing exposed hydrophobic patches to bind other unfolded proteins and prevent their irreversible aggregation. In this way, HdeA has a key role in helping pathogenic bacteria survive our acidic stomach and colonize our intestines, facilitating the spread of dysentery. Despite numerous publications on the topic, there remain questions about the mechanism by which HdeA unfolding and activation are triggered. Previous studies usually assessed HdeA unfolding over pH increments that are too far apart to gain fine detail of the process of unfolding and dimer dissociation, and often employed techniques that prevented thorough evaluation of specific regions of the protein. We used a variety of heteronuclear NMR experiments to investigate changes to backbone and side chain structure and dynamics of HdeA at four pHs between 3.0 and 2.0. We found that the long loop in the dimer interface is an early site of initiation of dimer dissociation, and that a molecular "clasp" near the disulfide bond is broken open at low pH as part, or as a trigger, of unfolding; this process also results in the separation of C-terminal helices and exposure of key hydrophobic client binding sites. Our results highlight important regions of HdeA that may have previously been overlooked because they lie too close to the disulfide bond or are thought to be too dynamic in the folded state to influence unfolding processes.

Exploring pH Dependent Host/Guest Binding Affinities

When the electrostatic environment surrounding binding partners changes between unbound and bound states, the net uptake or release of a proton is possible by either binding partner. This process is pH-dependent in that the free energy required to uptake or release the proton varies with pH. This pH-dependence is typically not considered in conventional free energy methods where the use of fixed protonation states is the norm. In the present paper, we apply a simple two-step approach to calculate the pH-dependent binding free energy of a model cucubit[7]uril host/guest system. By use of λ-dynamics with an enhanced sampling protocol, adaptive landscape flattening, pK_a shifts and reference binding free energies upon complexation were determined. This information enables the construction of pH-dependent binding profiles that accurately capture the pK_a shifts and reproduce binding free energies at the different pH conditions that were observed experimentally. Our calculations illustrate a general framework for computing pH-dependent binding free energies but also point to some issues in modeling the molecular charge distributions within this series of molecules with CGenFF. However, by introducing some minor charge modifications to the CGenFF force field, we saw significant improvement in accuracy of the calculated pK_a shifts.

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.

PCAViz: An Open-Source Python/JavaScript Toolkit for Visualizing Molecular Dynamics Simulations in the Web Browser

Molecular dynamics (MD) simulations reveal molecular motions at atomic resolution. Recent advances in high-performance computing now enable microsecond-long simulations capable of sampling a wide range of biologically relevant events. But the disk space required to store an MD trajectory increases with simulation length and system size, complicating collaborative sharing and visualization. To overcome these limitations, we created PCAViz, an open-source toolkit for sharing and visualizing MD trajectories via the web browser. PCAViz includes two components: the PCAViz Compressor, which compresses and saves simulation data; and the PCAViz Interpreter, which decompresses the data in users’ browsers and feeds it to any of several browser-based molecular-visualization libraries (e.g., 3Dmol.js, NGL Viewer, etc.). An easy-to-install WordPress plugin enables “plug-and-play” trajectory visualization. PCAViz will appeal to a broad audience of researchers and educators. The source code is available at http://durrantlab.com/pcaviz/, and the WordPress plugin is available via the official WordPress Plugin Directory.

q-Canonical Monte Carlo Sampling for Modeling the Linkage between Charge Regulation and Conformational Equilibria of Peptides

The overall charge content and the patterning of charged residues have a profound impact on the conformational ensembles adopted by intrinsically disordered proteins. These parameters can be altered by charge regulation, which refers to the effects of post-translational modifications, pH-dependent changes to charge, and conformational fluctuations that modify the pKa values of ionizable residues. Although atomistic simulations have played a prominent role in uncovering the major sequence-ensemble relationships of IDPs, most simulations assume fixed charge states for ionizable residues. This may lead to erroneous estimates for conformational equilibria if they are linked to charge regulation. Here, we report the development of a new method we term q-canonical Monte Carlo sampling for modeling the linkage between charge regulation and conformational equilibria. The method, which is designed to be interoperable with the ABSINTH implicit solvation model, operates as follows: For a protein sequence with n ionizable residues, we start with all 2n charge microstates and use a criterion based on model compound pKa values to prune down to a subset of thermodynamically relevant charge microstates. This subset is then grouped into mesostates, where all microstates that belong to a mesostate have the same net charge. Conformational distributions, drawn from a canonical ensemble, are generated for each of the charge microstates that make up a mesostate using a method we designate as proton walk sampling. This method combines Metropolis Monte Carlo sampling in conformational space with an auxiliary Markov process that enables interconversions between charge microstates along a mesostate. Proton walk sampling helps identify the most likely charge microstate per mesostate. We then use thermodynamic integration aided by the multistate Bennett acceptance ratio method to estimate the free energies for converting between mesostates. These free energies are then combined with the per-microstate weights along each mesostate to estimate standard state free energies and pH-dependent free energies for all thermodynamically relevant charge microstates. The results provide quantitative estimates of the probabilities and preferred conformations associated with every thermodynamically accessible charge microstate. We showcase the application of q-canonical sampling using two model systems. The results establish the soundness of the method and the importance of charge regulation in systems characterized by conformational heterogeneity.

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.

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.

Stress-Activated Chaperones: A First Line of Defense

Proteins are constantly challenged by environmental stress conditions that threaten their structure and function. Especially problematic are oxidative, acid, and severe heat stress which induce very rapid and widespread protein unfolding and generate conditions that make canonical chaperones and/or transcriptional responses inadequate to protect the proteome. We review here recent advances in identifying and characterizing stress-activated chaperones which are inactive under non-stress conditions but become potent chaperones under specific protein-unfolding stress conditions. We discuss the post-translational mechanisms by which these chaperones sense stress, and consider the role that intrinsic disorder plays in their regulation and function. We examine their physiological roles under both non-stress and stress conditions, their integration into the cellular proteostasis network, and their potential as novel therapeutic targets.

CHARMM-GUI 10 years for biomolecular modeling and simulation

CHARMM-GUI, http://www.charmm-gui.org, is a web-based graphical user interface that prepares complex biomolecular systems for molecular simulations. CHARMM-GUI creates input files for a number of programs including CHARMM, NAMD, GROMACS, AMBER, GENESIS, LAMMPS, Desmond, OpenMM, and CHARMM/OpenMM. Since its original development in 2006, CHARMM-GUI has been widely adopted for various purposes and now contains a number of different modules designed to set up a broad range of simulations: (1) PDB Reader & Manipulator, Glycan Reader, and Ligand Reader & Modeler for reading and modifying molecules; (2) Quick MD Simulator, Membrane Builder, Nanodisc Builder, HMMM Builder, Monolayer Builder, Micelle Builder, and Hex Phase Builder for building all-atom simulation systems in various environments; (3) PACE CG Builder and Martini Maker for building coarse-grained simulation systems; (4) DEER Facilitator and MDFF/xMDFF Utilizer for experimentally guided simulations; (5) Implicit Solvent Modeler, PBEQ-Solver, and GCMC/BD Ion Simulator for implicit solvent related calculations; (6) Ligand Binder for ligand solvation and binding free energy simulations; and (7) Drude Prepper for preparation of simulations with the CHARMM Drude polarizable force field. Recently, new modules have been integrated into CHARMM-GUI, such as Glycolipid Modeler for generation of various glycolipid structures, and LPS Modeler for generation of lipopolysaccharide structures from various Gram-negative bacteria. These new features together with existing modules are expected to facilitate advanced molecular modeling and simulation thereby leading to an improved understanding of the structure and dynamics of complex biomolecular systems. Here, we briefly review these capabilities and discuss potential future directions in the CHARMM-GUI development project. © 2016 Wiley Periodicals, Inc.

Conditional Chaperone-Client Interactions Revealed by Genetically Encoded Photo-cross-linkers

The cell envelope is an integral and essential component of Gram-negative bacteria. As the front line during host-pathogen interactions, it is directly challenged by host immune responses as well as other harsh extracellular stimuli. The high permeability of the outer-membrane and the lack of ATP energy system render it difficult to maintain important biological activities within the periplasmic space under stress conditions. The HdeA/B chaperone machinery is the only known acid resistant system found in bacterial periplasm, enabling enteric pathogens to survive through the highly acidic human stomach and establish infections in the intestine. These two homologous chaperones belong to a fast growing family of conditionally disordered chaperones that conditionally lose their well-defined three-dimensional structures to exert biological activities. Upon losing ordered structures, these proteins commit promiscuous binding of diverse clients in response to environmental stimulation. For example, HdeA and HdeB are well-folded inactive dimers at neutral pH but become partially unfolded to protect a wide array of acid-denatured proteins upon acid stress. Whether these conditionally disordered chaperones possess client specificities remains unclear. This is in part due to the lack of efficient tools to investigate such versatile and heterogeneous protein-protein interactions under living conditions. Genetically encoded protein photo-cross-linkers have offered a powerful strategy to capture protein-protein interactions, showing great potential in profiling protein interaction networks, mapping binding interfaces, and probing dynamic changes in both physiological and pathological settings. Despite great success, photo-cross-linkers that can simultaneously capture the promiscuous binding partners and directly identify the interaction interfaces remain technically challenging. Furthermore, methods for side-by-side profiling and comparing the condition-dependent client pools from two homologous chaperones are lacking. Herein, we introduce our recent efforts in developing a panel of versatile genetically encoded photo-cross-linkers to study the disorder-mediated chaperone-client interactions in living cells. In particular, we have developed a series of proteomic-based strategies relying on these new photo-cross-linkers to systematically compare the client profiles of HdeA and HdeB, as well as to map their interaction interfaces. These studies revealed the mode-of-action, particularly the client specificity, of these two conditionally disordered chaperones. In the end, some recent elegant work from other groups that applied the genetically encoded photo-cross-linking strategy to illuminate important protein-protein interactions within bacterial cell envelope is also 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.

Detection of the pH-dependent Activity of Escherichia coli Chaperone HdeB In Vitro and In Vivo

Bacteria are frequently exposed to environmental changes, such as alterations in pH, temperature, redox status, light exposure or mechanical force. Many of these conditions cause protein unfolding in the cell and have detrimental impact on the survival of the organism. A group of unrelated, stress-specific molecular chaperones have been shown to play essential roles in the survival of these stress conditions. While fully folded and chaperone-inactive before stress, these proteins rapidly unfold and become chaperone-active under specific stress conditions. Once activated, these conditionally disordered chaperones bind to a large number of different aggregation-prone proteins, prevent their aggregation and either directly or indirectly facilitate protein refolding upon return to non-stress conditions. The primary approach for gaining a more detailed understanding about the mechanism of their activation and client recognition involves the purification and subsequent characterization of these proteins using in vitro chaperone assays. Follow-up in vivo stress assays are absolutely essential to independently confirm the obtained in vitro results. This protocol describes in vitro and in vivo methods to characterize the chaperone activity of E. coli HdeB, an acid-activated chaperone. Light scattering measurements were used as a convenient read-out for HdeB's capacity to prevent acid-induced aggregation of an established model client protein, MDH, in vitro. Analytical ultracentrifugation experiments were applied to reveal complex formation between HdeB and its client protein LDH, to shed light into the fate of client proteins upon their return to non-stress conditions. Enzymatic activity assays of the client proteins were conducted to monitor the effects of HdeB on pH-induced client inactivation and reactivation. Finally, survival studies were used to monitor the influence of HdeB's chaperone function in vivo.

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.

Light-induced Changes in the Dimerization Interface of Bacteriophytochromes

Phytochromes are dimeric photoreceptor proteins that sense red light levels in plants, fungi, and bacteria. The proteins are structurally divided into a light-sensing photosensory module consisting of PAS, GAF, and PHY domains and a signaling output module, which in bacteriophytochromes typically is a histidine kinase (HK) domain. Existing structural data suggest that two dimerization interfaces exist between the GAF and HK domains, but their functional roles remain unclear. Using mutational, biochemical, and computational analyses of the Deinococcus radiodurans phytochrome, we demonstrate that two dimerization interfaces between sister GAF and HK domains stabilize the dimer with approximately equal contributions. The existence of both dimerization interfaces is critical for thermal reversion back to the resting state. We also find that a mutant in which the interactions between the GAF domains were removed monomerizes under red light. This implies that the interactions between the HK domains are significantly altered by photoconversion. The results suggest functional importance of the dimerization interfaces in bacteriophytochromes.

Membrane environment modulates the pKa values of transmembrane helices

In this work, we apply the recently developed constant pH molecular dynamics technique to study protonation equilibria of titratable side chains in the context of simple transmembrane (TM) helices and explore the effect of pH on their configurations in membrane bilayers. We observe that, despite a significant shift toward neutral states, considerable population of different side chains stay in the charged state that give rise to pK_a values around 9.6 for Asp and Glu and 4.5 to 6 for His and Lys side chains, respectively. These charged states are highly stabilized by favorable interactions between head groups, water molecules, and the charged side chains that are facilitated by substantial changes in the configuration of the peptides. The pH dependent configurations and the measured pK_a values are in good agreement with relatively recent solid state NMR measurements. Our results presented here demonstrate that all-atom constant pH molecular dynamics can be applied to membrane proteins and peptides to obtain reliable pK_a values and pH dependent behavior for these systems.

Contribution of the hdeB-like gene (SEN1493) to survival of Salmonella enterica enteritidis Nal(R) at pH 2

Periplasmic proteins are particularly vulnerable to denaturation upon entry into a highly acid environment. In Escherichia coli, a level of protection of these proteins is afforded by acid-inducible chaperonins encoded by hdeAB. In contrast, Salmonella enterica only harbors an hdeB-like gene and it is currently not known what function it plays in this genus. In the present study, the hdeB-like gene was deleted in Salmonella enterica Enteritidis Nal(R) and Salmonella enterica Kentucky 3795. When grown overnight in tryptic soy broth (TSB) medium buffered at pH 5.5 and then exposed to TSB pH 2 for 20 min, Enteritidis wild-type strain experienced a 0.5-log10 reduction in colony-forming units, whereas the deletion strain's surviving cells were reduced by 1.6 log10. No difference in survival was observed in the corresponding Salmonella enterica Kentucky 3795 strains treated the same way. Exposure of the strains to pH 2.5 or 3 resulted in the same log reduction regardless of the presence of the hdeB-like gene. When wild-type and deletion strains of both serovars were grown in medium buffered at pH 7 prior to exposure to the acidic pHs, no difference in survival with respect to serovar or presence/absence of the hdeB-like gene was found. Salmonella enterica Enteritidis Nal(R) carrying its own or the intragenic region upstream of the hdeB-like from Salmonella enterica Kentucky 3795 cloned in front of the gfp gene from pFPV25 showed maximum fluorescence when grown at pH 5.5, whereas the corresponding plasmid-carrying Salmonella enterica Kentucky strains did not exhibit fluorescence regardless of the pH of the growth medium. Therefore, the hdeB-like gene in Salmonella enterica Enteritidis, but not in Salmonella enterica Kentucky 3795, contributed to survival at pH 2 and its expression is responsive to the pH of the medium.

Residue-level resolution of alphavirus envelope protein interactions in pH-dependent fusion

Alphavirus envelope proteins, organized as trimers of E2-E1 heterodimers on the surface of the pathogenic alphavirus, mediate the low pH-triggered fusion of viral and endosomal membranes in human cells. The lack of specific treatment for alphaviral infections motivates our exploration of potential antiviral approaches by inhibiting one or more fusion steps in the common endocytic viral entry pathway. In this work, we performed constant pH molecular dynamics based on an atomic model of the alphavirus envelope with icosahedral symmetry. We have identified pH-sensitive residues that cause the largest shifts in thermodynamic driving forces under neutral and acidic pH conditions for various fusion steps. A series of conserved interdomain His residues is identified to be responsible for the pH-dependent conformational changes in the fusion process, and ligand binding sites in their vicinity are anticipated to be potential drug targets aimed at inhibiting viral infections.

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.

HdeB functions as an acid-protective chaperone in bacteria

Enteric bacteria such as Escherichia coli utilize various acid response systems to counteract the acidic environment of the mammalian stomach. To protect their periplasmic proteome against rapid acid-mediated damage, bacteria contain the acid-activated periplasmic chaperones HdeA and HdeB. Activation of HdeA at pH 2 was shown to correlate with its acid-induced dissociation into partially unfolded monomers. In contrast, HdeB, which has high structural similarities to HdeA, shows negligible chaperone activity at pH 2 and only modest chaperone activity at pH 3. These results raised intriguing questions concerning the physiological role of HdeB in bacteria, its activation mechanism, and the structural requirements for its function as a molecular chaperone. In this study, we conducted structural and biochemical studies that revealed that HdeB indeed works as an effective molecular chaperone. However, in contrast to HdeA, whose chaperone function is optimal at pH 2, the chaperone function of HdeB is optimal at pH 4, at which HdeB is still fully dimeric and largely folded. NMR, analytical ultracentrifugation, and fluorescence studies suggest that the highly dynamic nature of HdeB at pH 4 alleviates the need for monomerization and partial unfolding. Once activated, HdeB binds various unfolding client proteins, prevents their aggregation, and supports their refolding upon subsequent neutralization. Overexpression of HdeA promotes bacterial survival at pH 2 and 3, whereas overexpression of HdeB positively affects bacterial growth at pH 4. These studies demonstrate how two structurally homologous proteins with seemingly identical in vivo functions have evolved to provide bacteria with the means for surviving a range of acidic protein-unfolding conditions.

In silico insights into protein-protein interactions and folding dynamics of the saposin-like domain of Solanum tuberosum aspartic protease

The plant-specific insert is an approximately 100-residue domain found exclusively within the C-terminal lobe of some plant aspartic proteases. Structurally, this domain is a member of the saposin-like protein family, and is involved in plant pathogen defense as well as vacuolar targeting of the parent protease molecule. Similar to other members of the saposin-like protein family, most notably saposins A and C, the recently resolved crystal structure of potato (Solanum tuberosum) plant-specific insert has been shown to exist in a substrate-bound open conformation in which the plant-specific insert oligomerizes to form homodimers. In addition to the open structure, a closed conformation also exists having the classic saposin fold of the saposin-like protein family as observed in the crystal structure of barley (Hordeum vulgare L.) plant-specific insert. In the present study, the mechanisms of tertiary and quaternary conformation changes of potato plant-specific insert were investigated in silico as a function of pH. Umbrella sampling and determination of the free energy change of dissociation of the plant-specific insert homodimer revealed that increasing the pH of the system to near physiological levels reduced the free energy barrier to dissociation. Furthermore, principal component analysis was used to characterize conformational changes at both acidic and neutral pH. The results indicated that the plant-specific insert may adopt a tertiary structure similar to the characteristic saposin fold and suggest a potential new structural motif among saposin-like proteins. To our knowledge, this acidified PSI structure presents the first example of an alternative saposin-fold motif for any member of the large and diverse SAPLIP family.

Uncovering pH-dependent transient states of proteins with buried ionizable residues

The role of pH in regulating biological activity is ubiquitous, and understanding pH-mediated activity has traditionally relied on analyzing static biomolecular structures of highly populated ground states solved near physiological pH. However, recent advances have shown the increasing importance of transiently populated states, the characterization of which is extremely challenging but made plausible with the development of techniques such as relaxation dispersion NMR spectroscopy. To unlock the pH dependence of these transient states with atomistic-level details, we applied the recently developed explicit solvent constant pH molecular dynamics (CPHMD(MSλD)) framework to a series of staphylococcal nuclease (SNase) mutants with buried ionizable residues and probed their dynamics in different pH environments. Among our key findings is the existence of open states in all SNase mutants containing "buried" residues with highly shifted pKa's, where local solvation around the protonation site was observed. The calculated pKa demonstrated good agreement with experimental pKa's, with a low average unsigned error of 1.3 pKa units and correlation coefficient R(2) = 0.78. Sampling both open and closed states in their respective pH range, where they are expected to be dominant, was necessary to reproduce experimental pKa's, and in the most extreme examples of pKa shifts measured, it can be interpreted that the open-state structures are transient at physiological pH, contributing a small population of 1-2%. This suggests that buried ionizable residues can trigger conformational fluctuations that may be observed as transient-state structures at physiological pH. Furthermore, the coupled relationship of both open and closed states and their role in recapitulating macroscopic experimental observables suggest that structural analysis of buried residues may benefit from looking at structural pairs, as opposed to the conventional approach of looking at a single static ground-state conformation.

Recent development and application of constant pH molecular dynamics

Solution pH is a critical environmental factor for chemical and biological processes. Over the last decade, significant efforts have been made in the development of constant pH molecular dynamics (pHMD) techniques for gaining detailed insights into pH-coupled dynamical phenomena. In this article we review the advancement of this field in the past five years, placing a special emphasis on the development of the all-atom continuous pHMD technique. We discuss various applications, including the prediction of pK_a shifts for proteins, nucleic acids and surfactant assemblies, elucidation of pH-dependent population shifts, protein-protein and protein-RNA binding, as well as the mechanisms of pH-dependent self-assembly and phase transitions of surfactants and peptides. We also discuss future directions for the further improvement of the pHMD techniques.

Chaperone-like activity of monomeric human 14-3-3ζ on different protein substrates

Members of the 14-3-3 protein family interact with hundreds of different, predominantly phosphorylated, proteins. 14-3-3 dimers are prevalent but exist at the equilibrium with the monomers. Our previous studies using the engineered monomeric 14-3-3ζ (14-3-3ζm) showed that 14-3-3ζ monomer retained binding activity towards selected phosphorylated partners and, in addition, it prevented heat-induced aggregation of myosin subfragment 1. Since the chaperone-like activity of 14-3-3 monomers has been insufficiently studied, here we have analyzed the effect of 14-3-3ζm on the aggregation of different model proteins. We found that 14-3-3ζm demonstrated considerable chaperone-like activity by inhibiting the DTT-induced aggregation of insulin and thermally-induced aggregation of alcohol dehydrogenase and phosphorylase kinase. Importantly, the anti-aggregating activity of 14-3-3ζm was concentration-dependent and overall, was more pronounced than that of its dimeric counterpart. In some cases, the chaperone-like effect of 14-3-3ζm was comparable, or even higher, than that of the small heat shock proteins, HspB6 and HspB5. We suggest that 14-3-3s not only can bind and regulate the activity of multiple phosphoproteins, but also possess moonlighting chaperone-like activity, which is especially pronounced in the case of monomeric forms of 14-3-3 which can be present under certain stress conditions.

Constant pH molecular dynamics of proteins in explicit solvent with proton tautomerism

pH is a ubiquitous regulator of biological activity, including protein-folding, protein-protein interactions, and enzymatic activity. Existing constant pH molecular dynamics (CPHMD) models that were developed to address questions related to the pH-dependent properties of proteins are largely based on implicit solvent models. However, implicit solvent models are known to underestimate the desolvation energy of buried charged residues, increasing the error associated with predictions that involve internal ionizable residue that are important in processes like hydrogen transport and electron transfer. Furthermore, discrete water and ions cannot be modeled in implicit solvent, which are important in systems like membrane proteins and ion channels. We report on an explicit solvent constant pH molecular dynamics framework based on multi-site λ-dynamics (CPHMD(MSλD)). In the CPHMD(MSλD) framework, we performed seamless alchemical transitions between protonation and tautomeric states using multi-site λ-dynamics, and designed novel biasing potentials to ensure that the physical end-states are predominantly sampled. We show that explicit solvent CPHMD(MSλD) simulations model realistic pH-dependent properties of proteins such as the Hen-Egg White Lysozyme (HEWL), binding domain of 2-oxoglutarate dehydrogenase (BBL) and N-terminal domain of ribosomal protein L9 (NTL9), and the pKa predictions are in excellent agreement with experimental values, with a RMSE ranging from 0.72 to 0.84 pKa units. With the recent development of the explicit solvent CPHMD(MSλD) framework for nucleic acids, accurate modeling of pH-dependent properties of both major class of biomolecules-proteins and nucleic acids is now possible.

NMR-monitored titration of acid-stress bacterial chaperone HdeA reveals that Asp and Glu charge neutralization produces a loosened dimer structure in preparation for protein unfolding and chaperone activation

HdeA is a periplasmic chaperone found in several gram-negative pathogenic bacteria that are linked to millions of cases of dysentery per year worldwide. After the protein becomes activated at low pH, it can bind to other periplasmic proteins, protecting them from aggregation when the bacteria travel through the stomach on their way to colonize the intestines. It has been argued that one of the major driving forces for HdeA activation is the protonation of aspartate and glutamate side chains. The goal for this study, therefore, was to investigate, at the atomic level, the structural impact of this charge neutralization on HdeA during the transition from near-neutral conditions to pH 3.0, in preparation for unfolding and activation of its chaperone capabilities. NMR spectroscopy was used to measure pKa values of Asp and Glu residues and monitor chemical shift changes. Measurements of R2/R1 ratios from relaxation experiments confirm that the protein maintains its dimer structure between pH 6.0 and 3.0. However, calculated correlation times and changes in amide protection from hydrogen/deuterium exchange experiments provide evidence for a loosening of the tertiary and quaternary structures of HdeA; in particular, the data indicate that the dimer structure becomes progressively weakened as the pH decreases. Taken together, these results provide insight into the process by which HdeA is primed to unfold and carry out its chaperone duties below pH 3.0, and it also demonstrates that neutralization of aspartate and glutamate residues is not likely to be the sole trigger for HdeA dissociation and unfolding.

Deconstructing activation events in rhodopsin

Activation of class-A G-protein-coupled receptors (GPCRs) involves large-scale reorganization of the H3/H6 interhelical network. In rhodopsin (Rh), this process is coupled to a change in the protonation state of a key residue, E134, whose exact role in activation is not well understood. Capturing this millisecond pH-dependent process is a well-appreciated challenge. We have developed a scheme combining the harmonic Fourier beads (HFB) method and constant-pH molecular dynamics with pH-based replica exchange (pH-REX) to gain insight into the structural changes that occur along the activation pathway as a function of the protonation state of E134. Our results indicate that E134 is protonated as a consequence of tilting of H6 by ca. 4.0° with respect to its initial position and simultaneous rotation by ca. 23° along its principal axis. The movement of H6 is associated with breakage of the E247-R135 and R135-E134 salt bridges and concomitant release of the E134 side chain, which results in an increase in its pKa value above physiological pH. An increase in the hydrophobicity of the environment surrounding E134 leads to further tilting and rotation of H6 and upshift of the E134 pKa. Such atomic-level information, which is not accessible through experiments, refines the earlier proposed sequential model of Rh activation (see: Zaitseva, E.; et al. Sequential Rearrangement of Interhelical Networks Upon Rhodopsin Activation in Membranes: The Meta IIa Conformational Substate . J. Am. Chem. Soc. 2010, 132, 4815) and argues that the E134 protonation switch is both a cause and a consequence of the H6 motion.

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.

pH-sensitive residues in the p19 RNA silencing suppressor protein from carnation Italian ringspot virus affect siRNA binding stability

Tombusviruses, such as Carnation Italian ringspot virus (CIRV), encode a protein homodimer called p19 that is capable of suppressing RNA silencing in their infected hosts by binding to and sequestering short-interfering RNA (siRNA) away from the RNA silencing pathway. P19 binding stability has been shown to be sensitive to changes in pH but the specific amino acid residues involved have remained unclear. Using constant pH molecular dynamics simulations, we have identified key pH-dependent residues that affect CIRV p19-siRNA binding stability at various pH ranges based on calculated changes in the free energy contribution from each titratable residue. At high pH, the deprotonation of Lys60, Lys67, Lys71, and Cys134 has the largest effect on the binding stability. Similarly, deprotonation of several acidic residues (Asp9, Glu12, Asp20, Glu35, and/or Glu41) at low pH results in a decrease in binding stability. At neutral pH, residues Glu17 and His132 provide a small increase in the binding stability and we find that the optimal pH range for siRNA binding is between 7.0 and 10.0. Overall, our findings further inform recent experiments and are in excellent agreement with data on the pH-dependent binding profile.

Towards Accurate Prediction of Protonation Equilibrium of Nucleic Acids

The role of protonated nucleotides in modulating the pH-dependent properties of nucleic acids is one of the emerging frontiers in the field of nucleic acid biology. The recent development of a constant pH molecular dynamics simulation (CPHMD^MSλD) framework for simulating nucleic acids has provided a tool for realistic simulations of pH-dependent dynamics. We enhanced the CPHMD^MSλD framework with pH-based replica exchange (pH-REX), which significantly improves the sampling of both titration and spatial coordinates. The results from our pK_a calculations for the GAAA tetraloop, which was predicted with lower accuracy previously due to sampling challenges, demonstrates that pH-REX reduces the average unsigned error (AUE) to 0.7 pK_a units, and the error of the most poorly predicted residue A17 was drastically reduced from 2.9 to 1.2 pK_a unit. Lastly, we show that pH-REX CPHMD^MSλD simulations can be used to identify the dominant conformation of nucleic acid structures in alternate pH environments. This work suggests that pH-REX CPHMD^MSλD simulations provide a practical tool for predicting nucleic acid protonation equilibrium from first-principles, and offering structural and mechanistic insight into the study of pH-dependent properties of nucleic acids.

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.

pH-dependent dynamics of complex RNA macromolecules

The role of pH-dependent protonation equilibrium in modulating RNA dynamics and function is one of the key unanswered questions in RNA biology. Molecular dynamics (MD) simulations can provide insight into the mechanistic roles of protonated nucleotides, but it is only capable of modeling fixed protonation states and requires prior knowledge of the key residue's protonation state. Recently, we developed a framework for constant pH molecular dynamics simulations (CPHMDMSλD) of nucleic acids, where the nucleotides' protonation states are modeled as dynamic variables that are coupled to the structural dynamics of the RNA. In the present study, we demonstrate the application of CPHMDMSλD to the lead-dependent ribozyme; establishing the validity of this approach for modeling complex RNA structures. We show that CPHMDMSλD accurately predicts the direction of the pKa shifts and reproduces experimentally-measured microscopic pKa values with an average unsigned error of 1.3 pKa units. The effects of coupled titration states in RNA structures are modeled, and the importance of conformation sampling is highlighted. The general accuracy of CPHMDMSλD simulations in reproducing pH-dependent observables reported in this work demonstrates that constant pH simulations provides a powerful tool to investigate pH-dependent processes in nucleic acids.