Search PubMed⌕ Search

Biomedical subjects

V J Hilser

Publications and source records attributed to V J Hilser.

16 recordsLinked to original sources

The origin of pH-dependent changes in m-values for the denaturant-induced unfolding of proteins.

Denaturant-induced unfolding is one of the most prevalent means of evaluating the structural stability of proteins and of determining the energetic consequences of mutations or changes in solution conditions. In spite of the widespread use of this approach, controversies and inconsistencies still persist with regard to the interpretation of the results of such studies. For example, most proteins show either a significant increase or a decrease (as much as 100 %) in the denaturant-dependence of the free energy of unfolding (i.e. the m-value) under increasingly acidic conditions. The pH dependence of the m-value is given different interpretations depending on whether the m-values increase or decrease with decreasing pH. In cases where m-values decrease, the decrease is attributed to the presence of an intermediate that becomes transiently stabilized during the unfolding transition at low pH. Cases where m-values increase as pH is lowered are usually interpreted in terms of an increase in the amount of surface area exposed by the denatured state at low pH. We have developed a general thermodynamic model that accounts for both types of behavior in terms of an intermediate that is populated throughout the unfolding transition. The model provides a unified framework for explaining both types of observed behavior, and the validity of the model was tested through the analysis of the pH dependence of m-values of staphylococcal nuclease. According to the model, the observed increase in m-values with decreasing pH is consistent with the existence of an intermediate that is populated during urea and guanidine unfolding. The intermediate becomes less populated during the unfolding transition at lower pH values giving rise to the apparent increase in m-values. These results argue that the prevailing interpretation need not apply to all proteins.

Acids↗

Thermodynamic propensities of amino acids in the native state ensemble: implications for fold recognition.

An amino acid sequence, in the context of the solvent environment, contains all of the thermodynamic information necessary to encode a three-dimensional protein structure. To investigate the relationship between an amino acid sequence and its corresponding protein fold, a database of thermodynamic stability information was assembled that spanned 2951 residues from 44 nonhomologous proteins. This information was obtained using the COREX algorithm, which computes an ensemble-based description of the native state of a protein. It was observed that amino acid types partitioned unequally into high, medium, and low thermodynamic stability environments. Furthermore, these distributions were reproducible and were significantly different than those expected from random partitioning. To assess the structural importance of the distributions, simple fold-recognition experiments were performed based on a 3D-1D scoring matrix containing only COREX residue stability information. This procedure was able to recover amino acid sequences corresponding to correct target structures more effectively than scoring matrices derived from randomized data. High-scoring sequences were often aligned correctly with their corresponding target profiles, suggesting that calculated thermodynamic stability profiles have the potential to encode sequence information. As a control, identical fold-recognition experiments were performed on the same database of proteins using DSSP secondary structure information in the scoring matrix, instead of COREX residue stability information. The comparable performance of both approaches suggested that COREX residue stability information and secondary structure information could be of equivalent utility in more sophisticated fold-recognition techniques. The results of this work are a consequence of the idea that amino acid sequences fold not into single, rigidly stable structures but rather into thermodynamic ensembles best represented by a time-averaged structure.

Algorithms↗

Binding sites in Escherichia coli dihydrofolate reductase communicate by modulating the conformational ensemble.

To explore how distal mutations affect binding sites and how binding sites in proteins communicate, an ensemble-based model of the native state was used to define the energetic connectivities between the different structural elements of Escherichia coli dihydrofolate reductase. Analysis of this model protein has allowed us to identify two important aspects of intramolecular communication. First, within a protein, pair-wise couplings exist that define the magnitude and extent to which mutational effects propagate from the point of origin. These pair-wise couplings can be identified from a quantity we define as the residue-specific connectivity. Second, in addition to the pair-wise energetic coupling between residues, there exists functional connectivity, which identifies energetic coupling between entire functional elements (i.e., binding sites) and the rest of the protein. Analysis of the energetic couplings provides access to the thermodynamic domain structure in dihydrofolate reductase as well as the susceptibility of the different regions of the protein to both small-scale (e.g., point mutations) and large-scale perturbations (e. g., binding ligand). The results point toward a view of allosterism and signal transduction wherein perturbations do not necessarily propagate through structure via a series of conformational distortions that extend from one active site to another. Instead, the observed behavior is a manifestation of the distribution of states in the ensemble and how the distribution is affected by the perturbation.

Binding Sites↗

Ensemble modulation as an origin of denaturant-independent hydrogen exchange in proteins.

Native state hydrogen exchange (HX) has become a powerful tool for the analysis of conformational states that exist under native conditions. However, the interpretation of HX data in terms of conformational fluctuations is still controversial. In particular, it has been shown that many residues display exchange behavior that is independent of denaturant concentration. It has been postulated that this lack of denaturant dependence results from local fluctuations that do not expose appreciable amounts of buried surface area. Here, we use a general thermodynamic description of HX to explore the different possibilities for this behavior. We find that the denaturant dependence seen in HX experiments under native conditions is not a de facto indication of the amount of surface area exposure required for exchange. Instead, this behavior results from the relatively homogenous character of the conformational ensemble that exists under native conditions and the non-specific nature of denaturant effects. Furthermore, a comparison of the HX behavior from a stabilized mutant of Staphylococcal nuclease (SNase) with that predicted for the wild-type SNase from the COREX algorithm suggests that denaturant-independent exchange of many residues is consistent with significant (approximately 10 %) surface area exposure for this protein.

Algorithms↗

The structural distribution of cooperative interactions in proteins: analysis of the native state ensemble.

Cooperative interactions link the behavior of different amino acid residues within a protein molecule. As a result, the effects of chemical or physical perturbations to any given residue are propagated to other residues by an intricate network of interactions. Very often, amino acids "sense" the effects of perturbations occurring at very distant locations in the protein molecule. In these studies, we have investigated by computer simulation the structural distribution of those interactions. We show here that cooperative interactions are not intrinsically bi-directional and that different residues play different roles within the intricate network of interactions existing in a protein. The effect of a perturbation to residue j on residue k is not necessarily equal to the effect of the same perturbation to residue k on residue j. In this paper, we introduce a computer algorithm aimed at mapping the network of cooperative interactions within a protein. This algorithm exhaustively performs single site thermodynamic mutations to each residue in the protein and examines the effects of those mutations on the distribution of conformational states. The algorithm has been applied to three different proteins (lambda repressor fragment 6-85, chymotrypsin inhibitor 2, and barnase). This algorithm accounts well for the observed behavior of these proteins.

Algorithms↗

Structure-based statistical thermodynamic analysis of T4 lysozyme mutants: structural mapping of cooperative interactions.

The recent development of a structural parameterization of the energetics of protein folding has permitted the incorporation of the functions that describe the enthalpy, entropy and heat capacity changes, i.e. the individual components of the Gibbs energy, into a statistical thermodynamic formalism that describes the distribution of conformational states under equilibrium conditions. The goal of this approach is to construct with the computer a large ensemble of conformational states, and then to derive the most probable population distribution, i.e. the distribution of states that best accounts for a wide array of experimental observables. This analysis has been applied to four different mutants of T4 lysozyme (S44A, S44G, V131A, V131G). It is shown that the structural parameterization predicts well the stability of the protein and the effects of the mutations. The entire set of folding constants per residue has been calculated for the four mutants. In all cases, the effect of the mutations propagates beyond the mutation site itself through sequence and three-dimensional space. This phenomenon occurs despite the fact that the mutations are at solvent-exposed locations and do not directly affect other interactions in the protein. These results suggest that single amino acid mutations at solvent-exposed locations, or other locations that cause a minimal perturbation, can be used to identify the extent of cooperative interactions. The magnitude and extent of these effects and the accuracy of the algorithm can be tested by means of NMR-detected hydrogen exchange.

Bacteriophage T4↗

Predicting the equilibrium protein folding pathway: structure-based analysis of staphylococcal nuclease.

The equilibrium folding pathway of staphylococcal nuclease (SNase) has been approximated using a statistical thermodynamic formalism that utilizes the high-resolution structure of the native state as a template to generate a large ensemble of partially folded states. Close to 400,000 different states ranging from the native to the completely unfolded states were included in the analysis. The probability of each state was estimated using an empirical structural parametrization of the folding energetics. It is shown that this formalism predicts accurately the stability of the protein, the cooperativity of the folding/unfolding transition observed by differential scanning calorimetry (DSC) or urea denaturation and the thermodynamic parameters for unfolding. More importantly, this formalism provides a quantitative account of the experimental hydrogen exchange protection factors measured under native conditions for SNase. These results suggest that the computer-generated distribution of states approximates well the ensemble of conformations existing in solution. Furthermore, this formalism represents the first model capable of quantitatively predicting within a unified framework the probability distribution of states seen under native conditions and its change upon unfolding.

Enzyme Stability↗

Structure-based calculation of the equilibrium folding pathway of proteins. Correlation with hydrogen exchange protection factors.

A new statistical thermodynamic formalism has been developed in order to describe the equilibrium folding pathway of proteins. The resulting formalism allows calculation of the probabilities that individual amino acid residues will be in a native or native-like conformation for any given degree of folding of the protein molecule. The residue probabilities are defined by the probability distribution of conformational states and can be used to calculate experimental quantities like native-state, hydrogen exchange protection factors. A combinatorial algorithm aimed at generating a large ensemble of conformational states (10(4) to 10(6)) using the native structure as a template has been developed. The Gibbs energy and corresponding probability of each conformational state is estimated by using a previously developed structural parametrization of the energetics. The approach has been applied to five different proteins: hen egg-white lysozyme, equine lysozyme, bovine pancreatic trypsin inhibitor, staphylococcal nuclease and turkey ovomucoid third domain. The validity of the approach has been tested by comparing predicted and experimental hydrogen exchange protection factors. It is shown that for the above proteins 76%, 73%, 74%, 78% and 81% of all observed protection factors are predicted correctly. Furthermore, on average, the magnitude of the predicted protection factors, expressed as apparent free energies per residue deviate less than 1 kcal/mol from those obtained experimentally. These results represent the first attempt at predicting both the location and magnitude of hydrogen exchange protection factors from the high-resolution structure of a protein. The good agreement between experimental and predicted values has permitted a close examination of the nature of the equilibrium folding intermediates existing under conditions of maximal stability of the native state.

Animals↗

The magnitude of the backbone conformational entropy change in protein folding.

The magnitude of the conformational entropy change experienced by the peptide backbone upon protein folding was investigated experimentally and by computational analysis. Experimentally, two different pairs of mutants of a 33 amino acid peptide corresponding to the leucine zipper region of GCN4 were used for high-sensitivity microcalorimetric analysis. Each pair of mutants differed only by having alanine or glycine at a specific solvent-exposed position under conditions in which the differences in stability could be attributed to differences in the conformational entropy of the unfolded state. The mutants studied were characterized by different stabilities but had identical heat capacity changes of unfolding (delta Cp), identical solvent-related entropies of unfolding (delta Ssolv), and identical enthalpies of unfolding (delta H) at equivalent temperatures. Accordingly, the differences in stability between the different mutants could be attributed to differences in conformational entropy. The computational studies were aimed at generating the energy profile of backbone conformations as a function of the main chain dihedral angles phi and phi. The energy profiles permit a direct calculation of the probability distribution of different conformers and therefore of the conformational entropy of the backbone. The experimental results presented in this paper indicate that the presence of the methyl group in alanine reduces the conformational entropy of the peptide backbone by 2.46 +/- 0.2 cal/K. mol with respect to that of glycine, consistent with a 3.4-fold reduction in the number of allowed conformations in the alanine-containing peptides. Similar results were obtained from the energy profiles. The computational analysis also indicates that the addition of further carbon atoms to the side chain had only a small effect as long as the side chains were unbranched at position beta. A further reduction with respect to Ala of only 0.61 and 0.81 cal/K. mol in the backbone entropy was obtained for leucine and lysine, respectively. beta-branching (Val) produces the largest decrease in conformational entropy (1.92 cal/K.mol less than Ala). Finally, the backbone entropy change associated with the unfolding of an alpha-helix is 6.51 cal/K.mol for glycine. These and previous results have allowed a complete estimation of the conformational entropy changes associated with protein folding.

Alanine↗

The enthalpy change in protein folding and binding: refinement of parameters for structure-based calculations.

Two effects are mainly responsible for the observed enthalpy change in protein unfolding: the disruption of internal interactions within the protein molecule (van der Waals, hydrogen bonds, etc.) and the hydration of the groups that are buried in the native state and become exposed to the solvent on unfolding. In the traditional thermodynamic analysis, the effects of hydration have usually been evaluated using the thermodynamic data for the transfer of small model compounds from the gas phase to water. The contribution of internal interactions, on the other hand, are usually estimated by subtracting the hydration effects from the experimental enthalpy of unfolding. The main drawback of this approach is that the enthalpic contributions of hydration, and those due to the disruption of internal interactions, are more than one order of magnitude larger than the experimental enthalpy value. The enthalpy contributions of hydration and disruption of internal interactions have opposite signs and cancel each other almost completely resulting in a final value that is over 10 times smaller than the individual terms. For this reason, the classical approach cannot be used to accurately predict unfolding enthalpies from structure: any error in the estimation of the hydration enthalpy will be amplified by a factor of 10 or more in the estimation of the unfolding enthalpy. Recently, it has been shown that simple parametric equations that relate the enthalpy change with certain structural parameters, especially changes in solvent accessible surface areas have considerable predictive power. In this paper, we provide a physical foundation to that parametrization and in the process we present a system of equations that explicitly includes the enthalpic effects of the packing density between the different atoms within the protein molecule. Using this approach, the error in the prediction of folding/unfolding enthalpies at 60 degrees C, the median temperature for thermal unfolding, is better than +/- 3% (standard deviation = 4 kcal/mol).

Protein Conformation↗

The heat capacity of proteins.

The heat capacity plays a major role in the determination of the energetics of protein folding and molecular recognition. As such, a better understanding of this thermodynamic parameter and its structural origin will provide new insights for the development of better molecular design strategies. In this paper we have analyzed the absolute heat capacity of proteins in different conformations. The results of these studies indicate that three major terms account for the absolute heat capacity of a protein: (1) one term that depends only on the primary or covalent structure of a protein and contains contributions from vibrational frequencies arising from the stretching and bending modes of each valence bond and internal rotations; (2) a term that contains the contributions of noncovalent interactions arising from secondary and tertiary structure; and (3) a term that contains the contributions of hydration. For a typical globular protein in solution the bulk of the heat capacity at 25 degrees C is given by the covalent structure term (close to 85% of the total). The hydration term contributes about 15 and 40% to the total heat capacity of the native and unfolded states, respectively. The contribution of non-covalent structure to the total heat capacity of the native state is positive but very small and does not amount to more than 3% at 25 degrees C. The change in heat capacity upon unfolding is primarily given by the increase in the hydration term (about 95%) and to a much lesser extent by the loss of noncovalent interactions (up to approximately 5%).(ABSTRACT TRUNCATED AT 250 WORDS)

Amino Acids↗

Comparison of the utility of capillary zone electrophoresis and high-performance liquid chromatography in peptide mapping and separation.

Capillary zone electrophoresis (CZE) and high-performance liquid chromatography (HPLC) have been used in the analysis of the primary structure of recombinant human insulin-like growth factor I (rhIGF-I). CZE both complements and supplements HPLC separations. CZE has been used to resolve peaks which co-elute on HPLC, as well as to help establish the identity of tryptic fragments in peptide mapping experiments.

Amino Acid Sequence↗

Protein and peptide mobility in capillary zone electrophoresis. A comparison of existing models and further analysis.

Capillary zone electrophoresis of peptide fragments from the tryptic digest of human recombinant insulin-like growth factor I (rhIGF-I) has been carried out and the observed mobilities used to compare the relative applicability of existing mobility models. In addition, the physical forces affecting electromigration have been systematically analyzed in order to more accurately describe the physical chemistry involved. Such an approach should further improve the ability to predict electrophoretic mobility in capillary zone electrophoresis.

Amino Acid Sequence↗

Analysis of thermally induced protein folding/unfolding transitions using free solution capillary electrophoresis.

It is shown that free solution capillary electrophoresis (FSCE) can be used to monitor the temperature-dependent folding/unfolding transitions of proteins. Furthermore, analysis of the data obtained by FSCE can be used to estimate the apparent thermodynamic parameters (enthalpy change (delta HvH), entropy change (delta S), and transition temperature (Tm)) associated with the folding/unfolding transition. In addition to mobility changes associated with the transition, FSCE analysis is unique in its ability to provide access to the population distribution of mobility states. This is demonstrated by the temperature-dependent change in the electrophoretic peak width and by the appearance of multiple peaks for very slow equilibrium or irreversible processes. Moreover, by comparing the mobility of the denatured state to that of unstructured model peptides, it is possible to characterize the relative degree of structure present in the unfolded state of a protein. This methodology has been applied to the analysis of the thermally induced unfolding of lysozyme at low pH. It is shown that the mobility of thermally denatured lysozyme can be described by the same function that describes unstructured, fully solvated peptides. On the contrary, the mobility of the native lysozyme is significantly higher than the value predicted by the same function. The accuracy of the apparent thermodynamic parameters obtained by this methodology compare within error with values obtained by direct calorimetric measurements using differential scanning calorimetry.

Amino Acid Sequence↗