GalaxyHomomer: a web server for protein homo-oligomer structure prediction from a monomer sequence or structure.

Homo-oligomerization of proteins is abundant in nature, and is often intimately related with the physiological functions of proteins, such as in metabolism, signal transduction or immunity. Information on the homo-oligomer structure is therefore important to obtain a molecular-level understanding of protein functions and their regulation. Currently available web servers predict protein homo-oligomer structures either by template-based modeling using homo-oligomer templates selected from the protein structure database or by ab initio docking of monomer structures resolved by experiment or predicted by computation. The GalaxyHomomer server, freely accessible at http://galaxy.seoklab.org/homomer, carries out template-based modeling, ab initio docking or both depending on the availability of proper oligomer templates. It also incorporates recently developed model refinement methods that can consistently improve model quality. Moreover, the server provides additional options that can be chosen by the user depending on the availability of information on the monomer structure, oligomeric state and locations of unreliable/flexible loops or termini. The performance of the server was better than or comparable to that of other available methods when tested on benchmark sets and in a recent CASP performed in a blind fashion.

A Cooperative Folding Unit as the Structural Link for Energetic Coupling within a Protein.

Previously, we demonstrated that binding of a ligand to Escherichia coli cofactor-dependent phosphoglycerate mutase (dPGM), a homodimeric protein, is energetically coupled with dimerization. The equilibrium unfolding of dPGM occurs with a stable, monomeric intermediate. Binding of several nonsubstrate metabolites stabilizes the dimeric native form over the monomeric intermediate, reducing the population of the intermediate. Both the active site and the dimer interface appear to be unfolded in the intermediate. We hypothesized that a loop containing residues 118-152 was responsible for the energetic coupling between the dimer interface and the distal active site and was unfolded in the intermediate. Here, we investigated the structure of the dPGM intermediate by probing side-chain interactions and solvent accessibility of the peptide backbone. By comparing the effect of a mutation on the global stability and the stability of the intermediate, we determine an equilibrium φ value (φeq value), which provides information about whether side-chain interactions are retained or lost in the intermediate. Hydrogen/deuterium exchange coupled with mass spectrometry (HDX-MS) was used to investigate differences in the solvent accessibility of the peptide backbone in the intermediate and native forms of dPGM. The results of φeq value analysis and HDX-MS reveal the least stable folding unit of dPGM, which is unfolded in the intermediate and links the active site to the dimer interface. The structure of the intermediate reveals how the cooperative network of residues in dPGM gives rise to the observed energetic coupling between dimerization and ligand binding.

Discovery of Nicotinamide Adenine Dinucleotide Binding Proteins in the Escherichia coli Proteome Using a Combined Energetic- and Structural-Bioinformatics-Based Approach.

Protein-ligand interaction plays a critical role in regulating the biochemical functions of proteins. Discovering protein targets for ligands is vital to new drug development. Here, we present a strategy that combines experimental and computational approaches to identify ligand-binding proteins in a proteomic scale. For the experimental part, we coupled pulse proteolysis with filter-assisted sample preparation (FASP) and quantitative mass spectrometry. Under denaturing conditions, ligand binding affected protein stability, which resulted in altered protein abundance after pulse proteolysis. For the computational part, we used the software Patch-Surfer2.0. We applied the integrated approach to identify nicotinamide adenine dinucleotide (NAD)-binding proteins in the Escherichia coli proteome, which has over 4200 proteins. Pulse proteolysis and Patch-Surfer2.0 identified 78 and 36 potential NAD-binding proteins, respectively, including 12 proteins that were consistently detected by the two approaches. Interestingly, the 12 proteins included 8 that are not previously known as NAD binders. Further validation of these eight proteins showed that their binding affinities to NAD computed by AutoDock Vina are higher than their cognate ligands and also that their protein ratios in the pulse proteolysis are consistent with known NAD-binding proteins. These results strongly suggest that these eight proteins are indeed newly identified NAD binders.

Energetic Coupling between Ligand Binding and Dimerization in Escherichia coli Phosphoglycerate Mutase.

Energetic coupling of two molecular events in a protein molecule is ubiquitous in biochemical reactions mediated by proteins, such as catalysis and signal transduction. Here, we investigate energetic coupling between ligand binding and folding of a dimer using a model system that shows three-state equilibrium unfolding of an exceptional quality. The homodimeric Escherichia coli cofactor-dependent phosphoglycerate mutase (dPGM) was found to be stabilized by ATP in a proteome-wide screen, although dPGM does not require or utilize ATP for enzymatic function. We investigated the effect of ATP on the thermodynamic stability of dPGM using equilibrium unfolding. We found that, in the absence of ATP, dPGM populates a partially unfolded, monomeric intermediate during equilibrium unfolding. However, addition of 1.0 mM ATP drastically reduces the population of the intermediate by selectively stabilizing the native dimer. Using a computational ligand docking method, we predicted ATP binds to the active site of the enzyme using the triphosphate group. By performing equilibrium unfolding and isothermal titration calorimetry with active-site variants of dPGM, we confirmed that active-site residues are involved in ATP binding. Our findings show that ATP promotes dimerization of the protein by binding to the active site, which is distal from the dimer interface. This cooperativity suggests an energetic coupling between the active site and the dimer interface. We also propose a structural link to explain how ligand binding to the active site is energetically coupled with dimerization.

Bacteriorhodopsin folds through a poorly organized transition state.

The folding mechanisms of helical membrane proteins remain largely uncharted. Here we characterize the kinetics of bacteriorhodopsin folding and employ φ-value analysis to explore the folding transition state. First, we developed and confirmed a kinetic model that allowed us to assess the rate of folding from SDS-denatured bacteriorhodopsin (bRU) and provides accurate thermodynamic information even under influence of retinal hydrolysis. Next, we obtained reliable φ-values for 16 mutants of bacteriorhodopsin with good coverage across the protein. Every φ-value was less than 0.4, indicating the transition state is not uniquely structured. We suggest that the transition state is a loosely organized ensemble of conformations.

Thermodynamic stability of bacteriorhodopsin mutants measured relative to the bacterioopsin unfolded state.

The stability of bacteriorhodopsin (bR) has often been assessed using SDS unfolding assays that monitor the transition of folded bR (bR(f)) to unfolded (bR(u)). While many criteria suggest that the unfolding curves reflect thermodynamic stability, slow retinal (RET) hydrolysis during refolding makes it impossible to perform the most rigorous test for equilibrium, i.e., superimposable unfolding and refolding curves. Here we made a new equilibrium test by asking whether the refolding rate in the transition zone is faster than RET hydrolysis. We find that under conditions we have used previously, refolding is in fact slower than hydrolysis, strongly suggesting that equilibrium is not achieved. Instead, the apparent free energy values reported previously are dominated by unfolding rates. To assess how different the true equilibrium values are, we employed an alternative method by measuring the transition of bR(f) to unfolded bacterioopsin (bO(u)), the RET-free form of unfolded protein. The bR(f)-to-bO(u) transition is fully reversible, particular when we add excess RET. We compared the difference in unfolding free energies for 13 bR mutants measured by both assays. For 12 of the 13 mutants with a wide range of stabilities, the results are essentially the same within experimental error. The congruence of the results is fortuitous and suggests the energetic effects of most mutations may be focused on the folded state. The bR(f)-to-bO(u) reaction is inconvenient because many days are required to reach equilibrium, but it is the preferable measure of thermodynamic stability. This article is part of a Special Issue entitled: Protein Folding in Membranes.

Chaperone action of a cofactor in protein folding.

Previously, we have reported that ATP accelerates the folding and unfolding of Escherichia coli glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which is a glycolytic enzyme utilizing NAD+ as a cofactor. Because ATP and NAD+ share the ADP moiety, we hypothesized that NAD+ also accelerates the folding of GAPDH and that the common structural motif between ATP and NAD+ is responsible for the chaperone activity. After confirming that NAD+ indeed accelerates the folding of GAPDH, we examined the chaperone activity of the structural fragments of NAD+ (ADP, AMP, adenosine, and nicotinamide monophosphate). Our finding showed that ADP and AMP significantly speed up the folding of GAPDH, while adenosine and nicotinamide monophosphate do not. ADP and AMP also dramatically speed up the unfolding of GAPDH by selectively stabilizing a transition state in which GAPDH has a partially unfolded conformation. Similar to the previously reported effect of ATP on the equilibrium unfolding of GAPDH, a partially unfolded intermediate also accumulates in the presence of ADP and AMP. Based on the effect of the structural fragments of NAD+ on the folding of GAPDH, we identified that AMP is the structural determinant of the chaperone activity of ATP and NAD+ . Also, we propose a plausible model to explain how NAD+ accelerates the folding of GAPDH through a stepwise development of molecular interactions with the protein.

Revisiting absorbance at 230nm as a protein unfolding probe.

Thermodynamic stability and unfolding kinetics of proteins are typically determined by monitoring protein unfolding with spectroscopic probes, such as circular dichroism (CD) and fluorescence. UV absorbance at 230nm (A(230)) is also known to be sensitive to protein conformation. However, its feasibility for quantitative analysis of protein energetics has not been assessed. Here we evaluate A(230) as a structural probe to determine thermodynamic stability and unfolding kinetics of proteins. By using Escherichia coli maltose binding protein (MBP) and E. coli ribonuclease H (RNase H) as our model proteins, we monitored their unfolding in urea and guanidinium chloride with A(230). Significant changes in A(230) were observed with both proteins on unfolding in the chemical denaturants. The global stabilities were successfully determined by measuring the change in A(230) in varying concentrations of denaturants. Also, unfolding kinetics was investigated by monitoring the change in A(230) under denaturing conditions. The results were quite consistent with those determined by CD. Unlike CD, A(230) allowed us to monitor protein unfolding in a 96-well microtiter plate with a UV plate reader. Our finding suggests that A(230) is a valid and convenient structural probe to determine thermodynamic stability and unfolding kinetics of proteins with many potential applications.

Energetics-based protein profiling on a proteomic scale: identification of proteins resistant to proteolysis.

Native states of proteins are flexible, populating more than just the unique native conformation. The energetics and dynamics resulting from this conformational ensemble are inherently linked to protein function and regulation. Proteolytic susceptibility is one feature determined by this conformational energy landscape. As an attempt to investigate energetics of proteins on a proteomic scale, we challenged the Escherichia coli proteome with extensive proteolysis and determined which proteins, if any, have optimized their energy landscape for resistance to proteolysis. To our surprise, multiple soluble proteins survived the challenge. Maltose binding protein, a survivor from thermolysin digestion, was characterized by in vitro biophysical studies to identify the physical origin of proteolytic resistance. This experimental characterization shows that kinetic stability is responsible for the unusual resistance in maltose binding protein. The biochemical functions of the identified survivors suggest that many of these proteins may have evolved extreme proteolytic resistance because of their critical roles under stressed conditions. Our results suggest that under functional selection proteins can evolve extreme proteolysis resistance by modulating their conformational energy landscapes without the need to invent new folds, and that proteins can be profiled on a proteomic scale according to their energetic properties by using proteolysis as a structural probe.

Salt bridge as a gatekeeper against partial unfolding.

Salt bridges are frequently observed in protein structures. Because the energetic contribution of salt bridges is strongly dependent on the environmental context, salt bridges are believed to contribute to the structural specificity rather than the stability. To test the role of salt bridges in enhancing structural specificity, we investigated the contribution of a salt bridge to the energetics of native-state partial unfolding in a cysteine-free version of Escherichia coli ribonuclease H (RNase H*). Thermolysin cleaves a protruding loop of RNase H(*) through transient partial unfolding under native conditions. Lys86 and Asp108 in RNase H(*) form a partially buried salt bridge that tethers the protruding loop. Investigation of the global stability of K86Q/D108N RNase H(*) showed that the salt bridge does not significantly contribute to the global stability. However, K86Q/D108N RNase H(*) is greatly more susceptible to proteolysis by thermolysin than wild-type RNase H(*) is. The free energy for partial unfolding determined by native-state proteolysis indicates that the salt bridge significantly increases the energy for partial unfolding by destabilizing the partially unfolded form. Double mutant cycles with single and double mutations of the salt bridge suggest that the partially unfolded form is destabilized due to a significant decrease in the interaction energy between Lys86 and Asp108 upon partial unfolding. This study demonstrates that, even in the case that a salt bridge does not contribute to the global stability, the salt bridge may function as a gatekeeper against partial unfolding that disturbs the optimal geometry of the salt bridge.

Effect of circular permutations on transient partial unfolding in proteins.

Under native conditions, proteins can undergo transient partial unfolding, which may cause proteins to misfold or aggregate. A change in sequence connectivity by circular permutation may affect the energetics of transient partial unfolding in proteins without altering the three-dimensional structures. Using Escherichia coli dihydrofolate reductase (DHFR) as a model system, we investigated how circular permutation affects transient partial unfolding in proteins. We constructed three circular permutants, CP18, CP37, and CP87, with the new N-termini at residue 18, 37, and 87, respectively, and probed transient partial unfolding by native-state proteolysis. The new termini in CP18, CP37, and CP87 are within, near, and distal to the Met20 loop, which is known to be dynamic and also part of the region that undergoes transient unfolding in wild-type DHFR. The stabilities of both native and partially unfolded forms of CP18 are similar to those of wild-type DHFR, suggesting that the influence of introducing new termini in a dynamic region to the protein is minimal. CP37 has a significantly more accessible partially unfolded form than wild-type DHFR, demonstrating that introducing new termini near a dynamic region may promote transient partial unfolding. CP87 has significantly destabilized native and partially unfolded forms, confirming that modification of the folded region in a partially unfolded form destabilizes the partially unfolded form similar to the native form. Our findings provide valuable guidelines to control transient partial unfolding in designing circular permutants in proteins.

Probing the high energy states in proteins by proteolysis.

Unless the native conformation has an unstructured region, proteases cannot effectively digest a protein under native conditions. Digestion must occur from a higher energy form, when at least some part of the protein is exposed to solvent and becomes accessible by proteases. Monitoring the kinetics and denaturant dependence of proteolysis under native conditions yields insight into the mechanism of proteolysis as well as these high-energy conformations. We propose here a generalized approach to exploit proteolysis as a tool to probe high-energy states in proteins. This "native state proteolysis" experiment was carried out on Escherichia coli ribonuclease HI. Mass spectrometry and N-terminal sequencing showed that thermolysin cleaves the peptide bond between Thr92 and Ala93 in an extended loop region of the protein. By comparing the proteolysis rate of the folded protein and a peptidic substrate mimicking the sequence at the cleavage site, the energy required to reach the susceptible state (Delta G(proteolysis)) was determined. From the denaturant dependence of Delta G(proteolysis), we determined that thermolysin digests this protein through a local fluctuation, i.e. localized unfolding with minimal change in solvent assessable surface area. Proteolytic susceptibilities of proteins are discussed based on the finding of this local fluctuation mechanism for proteolysis under native conditions.

Ligand binding to a high-energy partially unfolded protein.

The conformational energy landscape of a protein determines populations of all possible conformations of the protein and also determines the kinetics of the conversion between the conformations. Interaction with ligands influences the conformational energy landscapes of proteins and shifts populations of proteins in different conformational states. To investigate the effect of ligand binding on partial unfolding of a protein, we use Escherichia coli dihydrofolate reductase (DHFR) and its functional ligand NADP(+) as a model system. We previously identified a partially unfolded form of DHFR that is populated under native conditions. In this report, we determined the free energy for partial unfolding of DHFR at varying concentrations of NADP(+) and found that NADP(+) binds to the partially unfolded form as well as the native form. DHFR unfolds partially without releasing the ligand, though the binding affinity for NADP(+) is diminished upon partial unfolding. Based on known crystallographic structures of NADP(+) -bound DHFR and the model of the partially unfolded protein we previously determined, we propose that the adenosine-binding domain of DHFR remains folded in the partially unfolded form and interacts with the adenosine moiety of NADP(+) . Our result demonstrates that ligand binding may affect the conformational free energy of not only native forms but also high-energy non-native forms.

Analysis of the stability of multimeric proteins by effective DeltaG and effective m-values.

Analyzing the stability of a multimeric protein is challenging because of the intrinsic difficulty in handling the mathematical model for the folded multimer-unfolded monomer equilibrium. To circumvent this problem, we introduce the concept of effective stability, DeltaGeff (= -RTlnKeff), where Keff is the equilibrium constant expressed in monomer units. Analysis of the denaturant effect on DeltaGeff gives new insight into the stability of multimeric proteins. When a multimeric protein is mostly folded, the dependence of effective stability on denaturant concentration (effective m-value) is simply the m-value of its monomeric unit. However, when the protein is mostly unfolded, its stability depends on denaturant concentration with the m-value of its multimeric form. We also find that the effective m-value at the Cm is a good approximation of the apparent m-value determined by fitting the equilibrium unfolding data from multimeric proteins with a two-state monomer model. Moreover, when the m-value of a monomeric unit is estimated from its size, the effective stability of a multimeric protein can be determined simply from Cm and this estimated m-value. These simple and intuitive approaches will allow a facile analysis of the stability of multimeric proteins. These analyses are also applicable for high-throughput analysis of protein stability on a proteomic scale.

Product inhibition in native-state proteolysis.

The proteolysis kinetics of intact proteins by nonspecific proteases provides valuable information on transient partial unfolding of proteins under native conditions. Native-state proteolysis is an approach to utilize the proteolysis kinetics to assess the energetics of partial unfolding in a quantitative manner. In native-state proteolysis, folded proteins are incubated with nonspecific proteases, and the rate of proteolysis is determined from the disappearance of the intact protein. We report here that proteolysis of intact proteins by nonspecific proteases, thermolysin and subtilisin deviates from first-order kinetics. First-order kinetics has been assumed for the analysis of native-state proteolysis. By analyzing the kinetics of proteolysis with varying concentrations of substrate proteins and also with cleavage products, we found that the deviation from first-order kinetics results from product inhibition. A kinetic model including competitive product inhibition agrees well with the proteolysis time course and allows us to determine the uninhibited rate constant for proteolysis as well as the apparent inhibition constant. Our finding suggests that the likelihood of product inhibition must be considered for quantitative assessment of proteolysis kinetics.

Structure of a partially unfolded form of Escherichia coli dihydrofolate reductase provides insight into its folding pathway.

Proteins frequently fold via folding intermediates that correspond to local minima on the conformational energy landscape. Probing the structure of the partially unfolded forms in equilibrium under native conditions can provide insight into the properties of folding intermediates. To elucidate the structures of folding intermediates of Escherichia coli dihydrofolate reductase (DHFR), we investigated transient partial unfolding of DHFR under native conditions. We probed the structure of a high-energy conformation susceptible to proteolysis (cleavable form) using native-state proteolysis. The free energy for unfolding to the cleavable form is clearly less than that for global unfolding. The dependence of the free energy on urea concentration (m-value) also confirmed that the cleavable form is a partially unfolded form. By assessing the effect of mutations on the stability of the partially unfolded form, we found that native contacts in a hydrophobic cluster formed by the F-G and Met-20 loops on one face of the central ß-sheet are mostly lost in the partially unfolded form. Also, the folded region of the partially unfolded form is likely to have some degree of structural heterogeneity. The structure of the partially unfolded form is fully consistent with spectroscopic properties of the near-native kinetic intermediate observed in previous folding studies of DHFR. The findings suggest that the last step of the folding of DHFR involves organization in the structure of two large loops, the F-G and Met-20 loops, which is coupled with compaction of the rest of the protein.

Selective stabilization of a partially unfolded protein by a metabolite.

When proteins fold in vivo, the intermediates that exist transiently on their folding pathways are exposed to the potential interactions with a plethora of metabolites within the cell. However, these potential interactions are commonly ignored. Here, we report a case in which a ubiquitous metabolite interacts selectively with a nonnative conformation of a protein and facilitates protein folding and unfolding process. From our previous proteomics study, we have discovered that Escherichia coli glyceraldehyde-3-phosphate dehydrogenase (GAPDH), which is not known to bind ATP under native conditions, is apparently destabilized in the presence of a physiological concentration of ATP. To decipher the origin of this surprising effect, we investigated the thermodynamics and kinetics of folding and unfolding of GAPDH in the presence of ATP. Equilibrium unfolding of the protein in urea showed that a partially unfolded equilibrium intermediate accumulates in the presence of ATP. This intermediate has a quaternary structure distinct from the native protein. Also, ATP significantly accelerates the unfolding of GAPDH by selectively stabilizing a transition state that is distinct from the native state of the protein. Moreover, ATP also significantly accelerates the folding of GAPDH. These results demonstrate that ATP interacts specifically with a partially unfolded form of GAPDH and affects the kinetics of folding and unfolding of this protein. This unusual effect of ATP on the folding of GAPDH implies that endogenous metabolites may facilitate protein folding in vivo by interacting with partially unfolded intermediates.

Revisiting the folding kinetics of bacteriorhodopsin.

The elucidation of the physical principles that govern the folding and stability of membrane proteins is one of the greatest challenges in protein science. Several insights into the folding of α-helical membrane proteins have come from the investigation of the conformational equilibrium of H. halobium bacteriorhodopsin (bR) in mixed micelles using SDS as a denaturant. In an effort to confirm that folded bR and SDS-denatured bR reach the same conformational equilibrium, we found that bR folding is significantly slower than has been previously known. Interrogation of the effect of the experimental variables on folding kinetics reveals that the rate of folding is dependent not only on the mole fraction of SDS but also on the molar concentrations of mixed micelle components, a variable that was not controlled in the previous study of bR folding kinetics. Moreover, when the molar concentrations of mixed micelle components are fixed at the concentrations commonly employed for bR equilibrium studies, conformational relaxation in the transition zone is slower than hydrolysis of the retinal Schiff base. As a result, the conformational equilibrium between folded bR and SDS-denatured bR cannot be achieved under the conventional condition. Our finding suggests that the molar concentrations of mixed micelle components are important experimental variables in the investigation of the kinetics and thermodynamics of bR folding and should be accounted for to ensure the accurate assessment of the conformational equilibrium of bR without the interference of retinal hydrolysis.

Simplified proteomics approach to discover protein-ligand interactions.

Identifying targets of biologically active small molecules is an essential but still challenging task in drug research and chemical genetics. Energetics-based target identification is an approach that utilizes the change in the conformational stabilities of proteins upon ligand binding in order to identify target proteins. Different from traditional affinity-based capture approaches, energetics-based methods do not require any labeling or immobilization of the test molecule. Here, we report a surprisingly simple version of energetics-based target identification, which only requires ion exchange chromatography, SDS PAGE, and minimal use of mass spectrometry. The complexity of a proteome is reduced through fractionation by ion exchange chromatography. Urea-induced unfolding of proteins in each fraction is then monitored by the significant increase in proteolytic susceptibility upon unfolding in the presence and the absence of a ligand. Proteins showing a different degree of unfolding with the ligand are identified by SDS PAGE followed by mass spectrometry. Using this approach, we identified ATP-binding proteins in the Escherichia coli proteome. In addition to known ATP-binding proteins, we also identified a number of proteins that were not previously known to interact with ATP. To validate one such finding, we cloned and purified phosphoglyceromutase, which was not previously known to bind ATP, and confirmed that ATP indeed stabilizes this protein. The combination of fractionation and pulse proteolysis offers an opportunity to investigate protein-drug or protein-metabolite interactions on a proteomic scale with minimal instrumentation and without modification of a molecule of interest.

Energetics-based discovery of protein-ligand interactions on a proteomic scale.

Biochemical functions of proteins in cells frequently involve interactions with various ligands. Proteomic methods for the identification of proteins that interact with specific ligands such as metabolites, signaling molecules, and drugs are valuable in investigating the regulatory mechanisms of cellular metabolism, annotating proteins with unknown functions, and elucidating pharmacological mechanisms. Here we report an energetics-based target identification method in which target proteins in a cell lysate are identified by exploiting the effect of ligand binding on their stabilities. Urea-induced unfolding of proteins in cell lysates is probed by a short pulse of proteolysis, and the effect of a ligand on the amount of folded protein remaining is monitored on a proteomic scale. As proof of principle, we identified proteins that interact with ATP in the Escherichia coli proteome. Literature and database mining confirmed that a majority of the identified proteins are indeed ATP-binding proteins. Four identified proteins that were previously not known to interact with ATP were cloned and expressed to validate the result. Except for one protein, the effects of ATP on urea-induced unfolding were confirmed. Analyses of the protein sequences and structure models were also employed to predict potential ATP binding sites in the identified proteins. Our results demonstrate that this energetics-based target identification approach is a facile method to identify proteins that interact with specific ligands on a proteomic scale.