Mutational analysis of the energetics of the GrpE.DnaK binding interface: equilibrium association constants by sedimentation velocity analytical ultracentrifugation

New insights into the structure and function of the complex between the Escherichia coli Hsp70, DnaK, and its nucleotide-exchange factor, GrpE

The 70 kDa heat shock proteins (Hsp70s) play a pivotal role in many cellular functions using allosteric communication between their nucleotide-binding domain (NBD) and substrate-binding domain, mediated by an interdomain linker, to modulate their affinity for protein clients. Critical to modulation of the Hsp70 allosteric cycle, nucleotide-exchange factors (NEFs) act by a conserved mechanism involving binding to the ADP-bound NBD and opening of the nucleotide-binding cleft to accelerate the release of ADP and binding of ATP. The crystal structure of the complex between the NBD of the Escherichia coli Hsp70, DnaK, and its NEF, GrpE, was reported previously, but the GrpE in the complex carried a point mutation (G122D). Both the functional impact of this mutation and its location on the NEF led us to revisit the DnaK NBD/GrpE complex structurally using AlphaFold modeling and validation by solution methods that report on protein conformation and mutagenesis. This work resulted in a new model for the DnaK NBD in complex with GrpE in which subdomain IIB of the NBD rotates more than in the crystal structure, resulting in an open conformation of the nucleotide-binding cleft, which now resembles more closely what is seen in other Hsp/NEF complexes. Moreover, the new model is consistent with the increased ADP off-rate accompanying GrpE binding. Excitingly, our findings point to an interdomain allosteric signal in DnaK triggered by GrpE binding.

Functional analysis of the Chloroplast GrpE (CGE) proteins from Arabidopsis thaliana

The function of proteins depends on specific partners that regulate protein folding, degradation and protein-protein interactions, such partners are the chaperones and cochaperones. In chloroplasts, proteins belonging to several families of chaperones have been identified: chaperonins (Cpn60s), Hsp90s (Hsp90-5/Hsp90C), Hsp100s (Hsp93/ClpC) and Hsp70s (cpHsc70s). Several lines of evidence have demonstrated that cpHsc70 chaperones are involved in molecular processes like protein import, protein folding and oligomer formation that impact important physiological aspects in plants such as thermotolerance and thylakoid biogenesis. Despite the vast amount of data existing around the function of cpHcp70s chaperones, very little attention has been paid to the roles of DnaJ and GrpE cochaperones in the chloroplast. In this study, we performed a phylogenetic analysis of the chloroplastic GrpE (CGE) proteins from 71 species. Based on their phylogenetic relationships and on a motif enrichment analysis, we propose a classification system for land plants' CGEs, which include two independent groups with specific primary structure traits. Furthermore, using in vivo assays we determined that the two CGEs from A. thaliana (AtCGEs) complement the mutant phenotype displayed by a knockout E. coli strain defective in the bacterial grpE gene. Moreover, we determined in planta that the two AtCGEs are bona fide chloroplastic proteins, which form the essential homodimers needed to establish direct physical interactions with the cpHsc70-1 chaperone. Finally, we found evidence suggesting that AtCGE1 is involved in specific physiological phenomena in A. thaliana, such as the chloroplastic response to heat stress, and the correct oligomerization of the photosynthesis-related LHCII complex.

Two Arabidopsis Chloroplast GrpE Homologues Exhibit Distinct Biological Activities and Can Form Homo- and Hetero-Oligomers

Flowering plants have evolved two distinct clades of chloroplast GrpE homologues (CGEs), which are the nucleotide exchange factor for Hsp70. In Arabidopsis, they are named AtCGE1 (At5g17710) and AtCGE2 (At1g36390). Characterization of their corresponding T-DNA insertion mutants revealed that there is no visible change in phenotype except a defect in protein import in an AtCGE2-knockout mutant under normal growth conditions. However, the embryo development of an AtCGE1-knockout mutant was arrested early at the globular stage. An AtCGE1-knockdown mutant, harboring a T-DNA insertion in the 5'-UTR region, exhibited growth retardation and protein import defect, and its mutant phenotypes became more severe when AtCGE2 was further knocked out. Sub-organellar distribution implied that AtCGE2 might be important for membrane biology due to its preferential association with chloroplast membranes. Biochemical studies and complementation tests showed that only AtCGE1, but not AtCGE2, can effectively rescue the heat-sensitive phenotype of Escherichia coli grpE mutant and robustly stimulate the refolding of denatured luciferase by DnaK. Interestingly, AtCGE1 and AtCGE2 are tending to form heterocomplexes, which exhibit comparable co-chaperone activity to AtCGE1 homocomplexes. Our data indicate that AtCGE1 is the principle functional homologue of GrpE. The possibility that AtCGE2 has a subsidiary or regulatory function through homo- and/or hetero-oligomerization is discussed.

The use of analytical sedimentation velocity to extract thermodynamic linkage

For 25 years, the Gibbs Conference on Biothermodynamics has focused on the use of thermodynamics to extract information about the mechanism and regulation of biological processes. This includes the determination of equilibrium constants for macromolecular interactions by high precision physical measurements. These approaches further reveal thermodynamic linkages to ligand binding events. Analytical ultracentrifugation has been a fundamental technique in the determination of macromolecular reaction stoichiometry and energetics for 85 years. This approach is highly amenable to the extraction of thermodynamic couplings to small molecule binding in the overall reaction pathway. In the 1980s this approach was extended to the use of sedimentation velocity techniques, primarily by the analysis of tubulin-drug interactions by Na and Timasheff. This transport method necessarily incorporates the complexity of both hydrodynamic and thermodynamic nonideality. The advent of modern computational methods in the last 20 years has subsequently made the analysis of sedimentation velocity data for interacting systems more robust and rigorous. Here we review three examples where sedimentation velocity has been useful at extracting thermodynamic information about reaction stoichiometry and energetics. Approaches to extract linkage to small molecule binding and the influence of hydrodynamic nonideality are emphasized. These methods are shown to also apply to the collection of fluorescence data with the new Aviv FDS.

Rigorous assessment and integration of the sequence and structure based features to predict hot spots

Background

Systematic mutagenesis studies have shown that only a few interface residues termed hot spots contribute significantly to the binding free energy of protein-protein interactions. Therefore, hot spots prediction becomes increasingly important for well understanding the essence of proteins interactions and helping narrow down the search space for drug design. Currently many computational methods have been developed by proposing different features. However comparative assessment of these features and furthermore effective and accurate methods are still in pressing need.

Results

In this study, we first comprehensively collect the features to discriminate hot spots and non-hot spots and analyze their distributions. We find that hot spots have lower relASA and larger relative change in ASA, suggesting hot spots tend to be protected from bulk solvent. In addition, hot spots have more contacts including hydrogen bonds, salt bridges, and atomic contacts, which favor complexes formation. Interestingly, we find that conservation score and sequence entropy are not significantly different between hot spots and non-hot spots in Ab+ dataset (all complexes). While in Ab- dataset (antigen-antibody complexes are excluded), there are significant differences in two features between hot pots and non-hot spots. Secondly, we explore the predictive ability for each feature and the combinations of features by support vector machines (SVMs). The results indicate that sequence-based feature outperforms other combinations of features with reasonable accuracy, with a precision of 0.69, a recall of 0.68, an F1 score of 0.68, and an AUC of 0.68 on independent test set. Compared with other machine learning methods and two energy-based approaches, our approach achieves the best performance. Moreover, we demonstrate the applicability of our method to predict hot spots of two protein complexes.

Conclusion

Experimental results show that support vector machine classifiers are quite effective in predicting hot spots based on sequence features. Hot spots cannot be fully predicted through simple analysis based on physicochemical characteristics, but there is reason to believe that integration of features and machine learning methods can remarkably improve the predictive performance for hot spots.

Prediction of protein interaction hot spots using rough set-based multiple criteria linear programming

Protein-protein interactions are fundamentally important in many biological processes and it is in pressing need to understand the principles of protein-protein interactions. Mutagenesis studies have found that only a small fraction of surface residues, known as hot spots, are responsible for the physical binding in protein complexes. However, revealing hot spots by mutagenesis experiments are usually time consuming and expensive. In order to complement the experimental efforts, we propose a new computational approach in this paper to predict hot spots. Our method, Rough Set-based Multiple Criteria Linear Programming (RS-MCLP), integrates rough sets theory and multiple criteria linear programming to choose dominant features and computationally predict hot spots. Our approach is benchmarked by a dataset of 904 alanine-mutated residues and the results show that our RS-MCLP method performs better than other methods, e.g., MCLP, Decision Tree, Bayes Net, and the existing HotSprint database. In addition, we reveal several biological insights based on our analysis. We find that four features (the change of accessible surface area, percentage of the change of accessible surface area, size of a residue, and atomic contacts) are critical in predicting hot spots. Furthermore, we find that three residues (Tyr, Trp, and Phe) are abundant in hot spots through analyzing the distribution of amino acids.

Effect of kinetics on sedimentation velocity profiles and the role of intermediates

We have previously presented a tutorial on direct boundary fitting of sedimentation velocity data for kinetically mediated monomer-dimer systems [Correia and Stafford, 2009]. We emphasized the ability of Sedanal to fit for the k(off) values and measure their uncertainty at the 95% confidence interval. We concluded for a monomer-dimer system the range of well-determined k(off) values is limited to 0.005-10(-5) s(-1) corresponding to relaxation times of approximately 70 to approximately 33,000 s. More complicated reaction schemes introduce the potential complexity of low concentrations of an intermediate that may also influence the kinetic behavior during sedimentation. This can be seen in a cooperative ABCD system (A+B --> C; B+C --> D) where C, the 1:1 complex, is sparsely populated (K(1)=10(4) M(-1), K(2)=10(8) M(-1)). Under these conditions a k(1,off)<0.01 s(-1) produces slow kinetic features. The low concentration of species C contributes to this effect while still allowing the accurate estimation of k(1,off) (although k(2,off) can readily compensate and contribute to the kinetics). More complex reactions involving concerted assembly or cooperative ring formation with low concentrations of intermediate species also display kinetic effects due to a slow flux of material through the sparsely populated intermediate states. This produces a kinetically limited reaction boundary that produces partial resolution of individual species during sedimentation. Cooperativity of ring formation drives the reaction and thus separation of these two effects, kinetics and energetics, can be challenging. This situation is experimentally exhibited by systems that form large oligomers or rings and may especially contribute to formation of micelles and various protein aggregation diseases including formation of beta-amyloid and tau aggregates. Simulations, quantitative parameter estimation by direct boundary fitting and diagnostic features for these systems are presented with an emphasis on the features available in Sedanal to simulate and analyze kinetically mediated systems.

Sedimentation velocity ultracentrifugation analysis for hydrodynamic characterization of G-quadruplex structures

Analytical ultracentrifugation (AUC) is a powerful technique for the characterization of hydrodynamic and thermodynamic properties. The intent of this article is to demonstrate the utility of sedimentation velocity (SV) studies to obtain hydrodynamic information for G-quadruplex (GQ) systems and to provide insights into one part of this process, namely, data analysis of existing SV data. An array of data analysis software is available, mostly written and continually developed by established researchers in the AUC field, with particularly rapid advances in the analysis of SV data. Each program has its own learning curve, and this article is intended as a resource in the data analysis process for beginning researchers in the field. We discuss the application of three of the most commonly used data analysis programs, DCDT+, Sedfit, and SedAnal, to the interpretation of SV data obtained in our laboratory on two GQ systems.

Extracting equilibrium constants from kinetically limited reacting systems

It has been known for some time that slow kinetics will distort the shape of a reversible reaction boundary. Here we present a tutorial on direct boundary fitting of sedimentation velocity data for a monomer-dimer system that exhibits kinetic effects. Previous analysis of a monomer-dimer system suggested that rapid reaction behavior will persist until the relaxation time of the system exceeds 100 s (reviewed in Kegeles and Cann, 1978). Utilizing a kinetic integrator feature in Sedanal (Stafford and Sherwood, 2004), we can now fit for the k(off) values and measure the uncertainty at the 95% confidence interval. For the monomer-dimer system the range of well determined k(off) values is limited to 0.005 to 10(-5) s(-1) corresponding to relaxation times (at a loading concentration of the Kd) of approximately 70 to approximately 33,000 s. For shorter relaxation times the system is fast and only the equilibrium constant K but not k(off) can be uniquely determined. For longer relaxation times the system is irreversibly slow, and assuming the system was at initial equilibrium before the start of the run, only the equilibrium constant K but not k(off) can be uniquely determined.

Analytical ultracentrifugation: sedimentation velocity and sedimentation equilibrium

Analytical ultracentrifugation (AUC) is a versatile and powerful method for the quantitative analysis of macromolecules in solution. AUC has broad applications for the study of biomacromolecules in a wide range of solvents and over a wide range of solute concentrations. Three optical systems are available for the analytical ultracentrifuge (absorbance, interference, and fluorescence) that permit precise and selective observation of sedimentation in real time. In particular, the fluorescence system provides a new way to extend the scope of AUC to probe the behavior of biological molecules in complex mixtures and at high solute concentrations. In sedimentation velocity (SV), the movement of solutes in high centrifugal fields is interpreted using hydrodynamic theory to define the size, shape, and interactions of macromolecules. Sedimentation equilibrium (SE) is a thermodynamic method where equilibrium concentration gradients at lower centrifugal fields are analyzed to define molecule mass, assembly stoichiometry, association constants, and solution nonideality. Using specialized sample cells and modern analysis software, researchers can use SV to determine the homogeneity of a sample and define whether it undergoes concentration-dependent association reactions. Subsequently, more thorough model-dependent analysis of velocity and equilibrium experiments can provide a detailed picture of the nature of the species present in solution and their interactions.

GrpE N-terminal domain contributes to the interaction with Dnak and modulates the dynamics of the chaperone substrate binding domain

GrpE acts as a nucleotide exchange factor for DnaK, the main Hsp70 protein in bacteria, accelerating ADP/ATP exchange by several orders of magnitude. GrpE is a homodimer, each subunit containing three structural domains: a N-terminal unordered segment, two long coils and a C-terminal globular domain formed by a four-helix bundle, and a beta-subdomain. GrpE association to DnaK nucleotide-binding domain involves side-chain and backbone interactions located within the "headpiece" of the cochaperone, which consists of the C-terminal half of the coils, the four-helix bundle and the beta-subdomain. However, the role of the GrpE N-terminal region in the interaction with DnaK and the activity of the cochaperone remain controversial. In this study we explore the contribution of this domain to the binding reaction, using the wild-type proteins, two deletion mutants of GrpE (GrpE(34-197) and GrpE(69-197)) and the isolated DnaK nucleotide-binding domain. Analysis of the thermodynamic binding parameters obtained by isothermal titration calorimetry shows that both GrpE N-terminal segments, 1-33 and 34-68, contribute to the binding reaction. Partial proteolysis and substrate dissociation kinetics also suggest that the N-terminal half of GrpE coils (residues 34-68) interacts with DnaK interdomain linker, regulates the nucleotide exchange activity of the cochaperone and is required to stabilize DnaK-substrate complexes in the ADP-bound conformation.

Low resolution structure and stability studies of human GrpE#2, a mitochondrial nucleotide exchange factor

GrpE acts as a nucleotide exchange factor for the Hsp70 chaperone system. Only one GrpE isoform is present in Escherichia coli, but for reasons not yet well understood, two GrpE isoforms have been found in mammalian mitochondria.Therefore, studies aimed at evaluating the physico-chemical characteristics of these proteins are important for the comprehension of the function of the Hsp70 chaperone system in different organisms. Here we report biophysical studies on human mitochondrial GrpE isoform 2. Small angle X-ray scattering measurements of human GrpE isoform 2 showed that this protein has a quaternary structure which is similar to those of human GrpE isoform 1 and E. coli GrpE: a dimer with a cruciform elongated shape. However, mitochondrial isoforms differed from each other regarding chemical and thermal denaturation profiles. This fact, combined with results of distinct expression patterns previously reported, point out that these proteins may have different response to external stimuli. Our results also indicate that human GrpE isoform 2 is more similar to the GrpE from E. coli than to human GrpE isoform 1. These results are relevant because differences in the conformation of Hsp70 co-chaperones are considered to be one of the reasons for functional diversity of this system.

Mutational Analysis of Target Enzyme Recognition of the β-Trefoil Fold Barley α-Amylase/Subtilisin Inhibitor

The barley alpha-amylase/subtilisin inhibitor (BASI) inhibits alpha-amylase 2 (AMY2) with subnanomolar affinity. The contribution of selected side chains of BASI to this high affinity is discerned in this study, and binding to other targets is investigated. Seven BASI residues along the AMY2-BASI interface and four residues in the putative protease-binding loop on the opposite side of the inhibitor were mutated. A total of 15 variants were compared with the wild type by monitoring the alpha-amylase and protease inhibitory activities using Blue Starch and azoalbumin, respectively, and the kinetics of binding to target enzymes by surface plasmon resonance. Generally, the mutations had little effect on k(on), whereas the k(off) values were increased up to 67-fold. The effects on the inhibitory activity, however, were far more pronounced, and the K(i) values of some mutants on the AMY2-binding side increased 2-3 orders of magnitude, whereas mutations on the other side of the inhibitor had virtually no effect. The mutants K140L, D150N, and E168T lost inhibitory activity, revealing the pivotal role of charge interactions for BASI activity on AMY2. A fully hydrated Ca(2+) at the AMY2-BASI interface mediates contacts to the catalytic residues of AMY2. Mutations involving residues contacting the solvent ligands of this Ca(2+) had weaker affinity for AMY2 and reduced sensitivity to the Ca(2+) modulation of the affinity. These results suggest that the Ca(2+) and its solvation sphere are integral components of the AMY2-BASI complex, thus illuminating a novel mode of inhibition and a novel role for calcium in relation to glycoside hydrolases.