Biomedical subjects

Jose M Sanchez-Ruiz

Publications and source records attributed to Jose M Sanchez-Ruiz.

At least 19 recordsLinked to original sources

PubMed2006

Natural selection for kinetic stability is a likely origin of correlations between mutational effects on protein energetics and frequencies of amino acid occurrences in sequence alignments.

It appears plausible that natural selection constrains, to some extent at least, the stability in many natural proteins. If, during protein evolution, stability fluctuates within a comparatively narrow range, then mutations are expected to be fixed with frequencies that reflect mutational effects on stability. Indeed, we recently reported a robust correlation between the effect of 27 conservative mutations on the thermodynamic stability (unfolding free energy) of Escherichia coli thioredoxin and the frequencies of residues occurrences in sequence alignments. We show here that this correlation likely implies a lower limit to thermodynamic stability of only a few kJ/mol below the unfolding free energy of the wild-type (WT) protein. We suggest, therefore, that the correlation does not reflect natural selection of thermodynamic stability by itself, but of some other factor which is linked to thermodynamic stability for the mutations under study. We propose that this other factor is the kinetic stability of thioredoxin in vivo, since( i) kinetic stability relates to irreversible denaturation, (ii) the rate of irreversible denaturation in a crowded cellular environment (or in a harsh extracellular environment) is probably determined by the rate of unfolding, and (iii) the half-life for unfolding changes in an exponential manner with activation free energy and, consequently, comparatively small free energy effects can have deleterious consequences for kinetic stability. This proposal is supported by the results of a kinetic study of the WT form and the 27 single-mutant variants of E. coli thioredoxin based on the global analyses of chevron plots and equilibrium unfolding profiles determined from double-jump unfolding assays. This kinetic study suggests, furthermore, one of the factors that may contribute to the high activation free energy for unfolding in thioredoxin (required for kinetic stability), namely the energetic optimization of native-state residue environments in regions, which become disrupted in the transition state for unfolding.

Amino Acid Sequence↗

PubMed2006

Energetic and structural consequences of desolvation/solvation barriers to protein folding/unfolding assessed from experimental unfolding rates.

Theoretical work has suggested the existence of solvation/desolvation barriers in protein folding/unfolding processes. We propose that the energetic and structural consequences of such barriers for the folding transition state can be assessed from experimental unfolding rates using well-established structure-energetics relationships. For a set of proteins of size within the 60-130 number-of-residues range, we find energetic effects associated to solvation/desolvation on the order of 10(2) kJ/mol. This supports that the folding transition states may be characterized by large networks of water-unsatisfied, broken internal contacts. In terms of buried surface, we estimate the typical network size to be on the order of several thousands of A2, or approximately 50% of the total change in accessible surface area upon unfolding. The analyses reported here thus suggest a clear structural picture for the different energetic balance of native and folding transition states.

Computer Simulation↗

PubMed2006

Ligand effects on protein thermodynamic stability.

A simple, partition-function formalism is used to describe the coupling between ligand binding and protein equilibrium unfolding. This general theoretical framework is shown to provide an adequate basis for the analysis of experimental ligand effects on the unfolding of complex protein systems. Nevertheless, the most important consequences of ligand binding for protein thermodynamic stability, as exposed by the partition-function approach, are found to be those demonstrated by Julian Sturtevant about 20 years ago.

Humans↗

PubMed2006

Role of solvation barriers in protein kinetic stability.

The stability of several protein systems of interest has been shown to have a kinetic basis. Besides the obvious biotechnological implications, the general interest of understanding protein kinetic stability is emphasized by the fact that some emerging molecular approaches to the inhibition of amyloidogenesis focus on the increase of the kinetic stability of protein native states. Lipases are among the most important industrial enzymes. Here, we have studied the thermal denaturation of the wild-type form, four single-mutant variants and two highly stable, multiple-mutant variants of lipase from Thermomyces lanuginosa. In all cases, thermal denaturation was irreversible, kinetically controlled and conformed to the two-state irreversible model. This result supports that the novel molecular-dynamics-focused, directed-evolution approach involved in the preparation of the highly stable variants is successful likely because it addresses kinetic stability and, in particular, because heated molecular dynamics simulations possibly identify regions of disrupted native interactions in the transition state for irreversible denaturation. Furthermore, we find very large mutation effects on activation enthalpy and entropy, which were not accompanied by similarly large changes in kinetic urea m-value. From this we are led to conclude that these mutation effects are associated to some structural feature of the transition state for the irreversible denaturation process that is not linked to large changes in solvent accessibility. Recent computational studies have suggested the existence of solvation/desolvation barriers in at least some protein folding/unfolding processes. We thus propose that a solvation barrier (arising from the asynchrony between breaking of internal contacts and water penetration) may contribute to the kinetic stability of lipase from T. lanuginosa (and, possibly, to the kinetic stability of other proteins as well).

Ascomycota↗

PubMed2005

Direct measurement of barrier heights in protein folding.

A serious limitation in the study of protein folding reactions resides in the lack of experimental methods to measure the absolute height of the free energy barrier. This is particularly unfortunate given that if folding barriers are small, as theory predicts, it might be possible to resolve folding mechanisms directly. Here we explore the performance of a recently developed method to extract folding barriers from equilibrium differential scanning calorimetry (DSC) experiments. To test the method, we compare the thermodynamic barrier heights for 15 proteins obtained from available DSC data with the folding rates measured in kinetic experiments. The correlation between these two parameters is very good (r2 = 0.9) and has a slope consistent with the same energy scale. These results confirm that it is possible to measure free energy barriers for natural proteins from equilibrium DSC experiments. Furthermore, the measured barrier heights are small (<8 RT), in general, and marginal or nonexisting for fast-folding proteins.

Calorimetry, Differential Scanning↗

PubMed2005

A simple tool to explore the distance distribution of correlated mutations in proteins.

The analysis of correlated mutations in protein sequence alignments is of considerable interest, since it may provide useful energetic and even structural information (ideally, residue contacts). However, a number of recent experimental studies support the existence of long-distance communication in proteins, a fact that may lead to correlation between distant residues. We introduce in this work a simple statistical procedure to describe the relation structure--alignments on the basis of the residue--residue distance dependence of the number of residue couples over given thresholds of a correlation measure (such as a covariance value). This procedure may lead to clear pictures of the distance distribution of correlated mutations and may provide a simple but efficient tool to explore the different structural features that are reflected in the sequence alignments.

Amino Acids↗

PubMed2005

A stability pattern of protein hydrophobic mutations that reflects evolutionary structural optimization.

We have determined the effect of mutations involving isoleucine and valine (i.e., mutations I-->V and V-->I) on the stability of Escherichia coli thioredoxin. Despite the similarity in chemical structure (V and I differ only in a methyl group), we find that many environments are optimized to a significant extent for either V or I. We find, furthermore, that a plot of effect of hydrophobic mutations on stability versus packing density shows a strikingly simple pattern that clearly reflects evolutionary structural optimization. The existence of such patterns suggests the possibility of rationalizing (and perhaps even predicting) mutation effects on protein stability on the basis of evolutionary models. By "evolutionary model" we specifically refer in this context to a model for mutation effects on stability in which certain physical features of the mutated residue environments are evaluated from an assumption regarding how such environments have been selected during protein evolution (as opposed to a purely "physical model" in which those features would be derived from some kind of energetics analysis of the protein structural characteristics). To illustrate this novel approach and provide general guidelines for its application, we develop here a simple evolutionary model that successfully explains the effect of the I<-->V mutations on thioredoxin stability.

Algorithms↗

PubMed2005

Robustness of downhill folding: guidelines for the analysis of equilibrium folding experiments on small proteins.

Previously, we identified the protein BBL as a downhill folder. This conclusion was based on the statistical mechanical analysis of equilibrium experiments performed in two variants of BBL, one with a fluorescent label at the N-terminus, and another one labeled at both ends. A recent report has claimed that our results are an artifact of label-induced aggregation and that BBL with no fluorescent labels and a longer N-terminal tail folds in a two-state fashion. Here, we show that singly and doubly labeled BBL do not aggregate, unfold reversibly, and have the same thermodynamic properties when studied under appropriate experimental conditions (e.g., our original conditions (1)). With an elementary analysis of the available data on the nonlabeled BBL (2), we also show that this slightly more stable BBL variant is not a two-state folder. We discuss the problems that led to its previous misclassification and how they can be avoided. Finally, we investigate the equilibrium unfolding of the singly labeled BBL with both ends protected by acetylation and amidation. This variant has the same thermodynamic stability of the nonlabeled BBL and displays all the equilibrium signatures of downhill folding. From all these observations, we conclude that fluorescent labels do not perturb the thermodynamic properties of BBL, which consistently folds downhill regardless of its stability and specific protein tails. The work on BBL illustrates the shortcomings of applying conventional procedures intended to distinguish between two-state and three-state folding models to small fast-folding proteins.

Amino Acid Sequence↗

PubMed2004

The effect of charge-introduction mutations on E. coli thioredoxin stability.

Technological applications of proteins are often hampered by their low-stability and, consequently, the development of procedures for protein stabilization is of considerable biotechnological interest. Here, we use simple electrostatics to determine positions in E. coli thioredoxin at which mutations that introduce new charged residues are expected to lead to stability enhancement. We also obtain the corresponding mutants and characterize their stability using differential scanning calorimetry. The results are interpreted in terms of the accessibility in the native structure of the mutated residues and the potential effect of the mutations on the residual structure of the denatured state.

Calorimetry, Differential Scanning↗

PubMed2004

Empirical parametrization of pK values for carboxylic acids in proteins using a genetic algorithm.

Considerable effort has been devoted to the development of theoretical electrostatic methods to predict the pK values of ionizable residues in proteins. However, predictions appear often to be still at the qualitative or semi-quantitative level. We believe that, with the increasing number experimentally available pK values for proteins of known structure, an alternative approach becomes feasible: the empirical parametrization of the experimental protein pK database. Of course, in the long term, this empirical approach is no substitute for rigorous electrostatic analysis but, in the short term, it may prove to have useful predictive power and it may help to pinpoint the main structural determinants of pK values in proteins. Here we demonstrate the feasibility of the parametrization approach by fitting (using a genetic algorithm as fitting tool) the database for carboxylic acid pK values in proteins on the basis of an empirical equation that takes into account the two following kinds of effects: (1) long-range charge-charge interactions; (2) interactions of the given carboxylic acid group with its environment in the protein, which are described in terms of contributions from the different kind of atoms present in the protein (atomic contributions).

Algorithms↗

PubMed2004

Exploring protein-folding ensembles: a variable-barrier model for the analysis of equilibrium unfolding experiments.

Recent theoretical and experimental results point to the existence of small barriers to protein folding. These barriers can even be absent altogether, resulting in a continuous folding transition (i.e., downhill folding). With small barriers, the detailed properties of folding ensembles may become accessible to equilibrium experiments. However, further progress is hampered because folding experiments are interpreted with chemical models (e.g., the two-state model), which assume the existence of well defined macrostates separated by arbitrarily high barriers. Here we introduce a phenomenological model based on the classical Landau theory for critical transitions. In this physical model the height of the thermodynamic free energy barrier and the general properties of the folding ensemble are directly obtained from the experimental data. From the analysis of differential scanning calorimetry data alone, our model identifies the presence of a significant (>35 kJ/mol) barrier for the two-state protein thioredoxin and the absence of a barrier for BBL, a previously characterized downhill folding protein. These results illustrate the potential of our approach for extracting the general features of protein ensembles from equilibrium folding experiments.

Calorimetry, Differential Scanning↗

PubMed2004

Do proteins always benefit from a stability increase? Relevant and residual stabilisation in a three-state protein by charge optimisation.

The vast majority of our knowledge on protein stability arises from the study of simple two-state models. However, proteins displaying equilibrium intermediates under certain conditions abound and it is unclear whether the energetics of native/intermediate equilibria is well represented in current knowledge. We consider here that the overall conformational stability of three-state proteins is made of a "relevant" term and a "residual" one, corresponding to the free energy differences of the native to intermediate (N-to-I) and intermediate to denatured (I-to-D) equilibria, respectively. The N-to-I free energy difference is considered to be the relevant stability because protein-unfolding intermediates are likely devoid of biological activity. We use surface charge optimisation to first increase the overall (N-to-D) stability of a model three-state protein (apoflavodoxin) and then investigate whether the stabilisation obtained is realised into relevant or into residual stability. Most of the mutations designed from electrostatic calculations or from simple sequence conservation analysis produce large increases in the overall stability of the protein. However, in most cases, this simply leads to similarly large increases of the residual stability. Two mutations, nevertheless, show a different trend and increase the relevant stability of the protein substantially. When all the mutations are mapped onto the structure of the apoflavodoxin thermal-unfolding intermediate (obtained independently by equilibrium phi-analysis and NMR) they cluster perfectly so that the mutations increasing the relevant stability appear in the small unstructured region of the intermediate and the others in the native-like region. This illustrates the need for specific investigation of N-to-I equilibria and the structure of protein intermediates, and indicates that it is possible to rationally stabilise a protein against partial unfolding once the structure of the intermediate conformation is known, even if at low resolution.

Anabaena↗

PubMed2004

Relation between protein stability, evolution and structure, as probed by carboxylic acid mutations.

Native proteins are marginally stable. Low thermodynamic stability may actually be advantageous, although the accumulation of neutral, destabilizing mutations may have also contributed to it. In any case, once marginal stability has been reached, it appears plausible that mutations at non-constrained positions become fixed in the course of evolution (due to random drift) with frequencies that roughly reflect the mutation effects on stability ("pseudo-equilibrium hypothesis"). We have found that all glutamate-->aspartate mutations in wild-type Escherichia coli thioredoxin are destabilizing, as well as most of the aspartate-->glutamate mutations. Furthermore, the effect of these mutations on thioredoxin thermodynamic stability shows a robust correlation with the frequencies of occurrence of the involved residues in several-hundred sequence alignments derived from a BLAST search. These results provide direct and quantitative experimental evidence for the pseudo-equilibrium hypothesis and should have general consequences for the interpretation of mutation effects on protein stability, as they suggest that residue environments in proteins may be optimized for stabilizing interactions to a remarkable degree of specificity. We also provide evidence that such stabilizing interactions may be detected in sequence alignments, and briefly discuss the implications of this possibility for the derivation of structural information (on native and denatured states) from comparative sequence analyses.

Aspartic Acid↗

PubMed2004

The efficiency of different salts to screen charge interactions in proteins: a Hofmeister effect?

Understanding the screening by salts of charge-charge interactions in proteins is important for at least two reasons: a), screening by intracellular salt concentration may modulate the stability and interactions of proteins in vivo; and b), the in vitro experimental estimation of the contributions from charge-charge interactions to molecular processes involving proteins is generally carried out on the basis of the salt effect on process energetics, under the assumption that these interactions are screened out by moderate salt concentrations. Here, we explore experimentally the extent to which the screening efficiency depends on the nature of the salt. To this end, we have carried out an energetic characterization of the effect of NaCl (a nondenaturing salt), guanidinium chloride (a denaturing salt), and guanidinium thiocyanate (a stronger denaturant) on the stability of the wild-type form and a T14K variant of Escherichia coli thioredoxin. Our results suggest that the efficiency of different salts to screen charge-charge interactions correlates with their denaturing strength and with the position of the constituent ions in the Hofmeister rankings. This result appears consistent with the plausible relation of the Hofmeister rankings with the extent of solute accumulation/exclusion from protein surfaces.

Algorithms↗

PubMed2004

Structural and stability effects of phosphorylation: Localized structural changes in phenylalanine hydroxylase.

Phosphorylation of phenylalanine hydroxylase (PAH) at Ser16 by cAMP-dependent protein kinase increases the basal activity of the enzyme and its resistance to tryptic proteolysis. The modeled structures of the full-length phosphorylated and unphosphorylated enzyme were subjected to molecular dynamics simulations, and we analyzed the energy of charge-charge interactions for individual ionizable residues in the final structures. These calculations showed that the conformational changes induced by incorporation of phosphate were localized and limited mostly to the region around the phosphoserine (Arg13-Asp17) and a region around the active site in the catalytic domain that includes residues involved in the binding of the iron and the substrate L-Phe (Arg270 and His285). The absence of a generalized conformational change was confirmed by differential scanning calorimetry, thermal-dependent circular dichroism, fluorescence spectroscopy, and limited chymotryptic proteolysis of the phosphorylated and unphosphorylated PAH. Our results explain the effect of phosphorylation of PAH on both the resistance to proteolysis specifically by trypsin-like enzymes and on the increase in catalytic efficiency.

Calorimetry, Differential Scanning↗

PubMed2003

Role of residual structure in the unfolded state of a thermophilic protein.

Ribonucleases H from the thermophilic bacterium Thermus thermophilus and the mesophile Escherichia coli demonstrate a dramatic and surprising difference in their change in heat capacity upon unfolding (DeltaCp degrees ). The lower DeltaCp degrees of the thermophilic protein directly contributes to its higher thermal denaturation temperature (Tm). We propose that this DeltaCp degrees difference originates from residual structure in the unfolded state of the thermophilic protein; we verify this hypothesis by using a mutagenic approach. Residual structure in the unfolded state may provide a mechanism for balancing a high Tm with the optimal thermodynamic stability for a protein's function. Structure in the unfolded state is shown to differentially affect the thermodynamic profiles of thermophilic and mesophilic proteins.

Escherichia coli↗

PubMed2003

Energetic evidence for formation of a pH-dependent hydrophobic cluster in the denatured state of Thermus thermophilus ribonuclease H.

NMR studies on the denatured states of proteins indicate that residual structure often resides predominantly in hydrophobic clusters. Such hydrophobic cluster formation implies burial of apolar surface and, consequently, is expected to cause a decrease in heat capacity. We report here that, in the case of ribonuclease H from the thermophile Thermus thermophilus, a sharp decrease in denatured-state heat capacity occurs at about pH 3.8; this result points to the formation of hydrophobic clusters triggered by the protonation of several (about four) carboxylic acid groups, and indicates that the burial of apolar surface is favored by the less hydrophilic character of the uncharged forms of Asp and Glu side-chains. The process is not accompanied by large changes in optically active structure, but appears to be highly cooperative, as indicated by the sharpness of the pH-induced transition in the heat capacity. This acid-induced hydrophobic burial in denatured T.thermophilus ribonuclease H is clearly reflected in the pH dependence of the denaturation temperature (i.e. an abrupt change of slope at about pH 3.8 is seen in the plot of denaturation temperature versus pH), supporting a role for such denatured-state hydrophobic clusters in protein stability. The finding of cooperative protonation of several groups coupled to surface burial in denatured T.thermophilus ribonuclease H emphasizes the potential complexity of denatured-state electrostatics and advises caution when attempting to predict denatured-state properties on the basis of simple electrostatic models. Finally, our results suggest a higher propensity for hydrophobic cluster formation in the denatured state of T.thermophilus ribonuclease H as compared with that of its mesophilic counterpart from Escherichia coli.

Anilino Naphthalenesulfonates↗

PubMed2002

Experimental identification of downhill protein folding.

Theory predicts the existence of barrierless protein folding. Without barriers, folding should be noncooperative and the degree of native structure should be coupled to overall protein stability. We investigated the thermal unfolding of the peripheral subunit binding domain from Escherichia coli's 2-oxoglutarate dehydrogenase multienzyme complex (termed BBL) with a combination of spectroscopic techniques and calorimetry. Each technique probed a different feature of protein structure. BBL has a defined three-dimensional structure at low temperatures. However, each technique showed a distinct unfolding transition. Global analysis with a statistical mechanical model identified BBL as a downhill-folding protein. Because of BBL's biological function, we propose that downhill folders may be molecular rheostats, in which effects could be modulated by altering the distribution of an ensemble of structures.

Acyltransferases↗