PubMed Health⌕ Search

Biomedical subjects

Beatriz Ibarra-Molero

Publications and source records attributed to Beatriz Ibarra-Molero.

12 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↗
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↗
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

Salt-bridges can stabilize but do not accelerate the folding of the homodimeric coiled-coil peptide GCN4-p1.

Double mutant cycle analysis was employed to ascertain the role of intra- and interchain salt-bridges in the folding and stability of the dimeric coiled-coil peptide, GCN4-p1, the 33-residue leucine zipper domain of the transcriptional activator GCN4. Equilibrium circular dichroism studies of the urea-induced unfolding reaction at neutral pH revealed that both types of ionic interactions, localized primarily in the N-terminal portion of the molecule, enhance the stability of the native coiled-coil. By contrast, comparable stopped-flow circular dichroism studies indicate that the salt-bridge interactions, with one possible exception, are not well formed in the transition state for folding. Although the E22Q/R25A double mutant failed to fold, fragmentation studies suggest that the E22/R25 intramolecular salt-bridge may play a critical role in stabilizing C-terminal nascent helices that drive the association reaction. The remaining salt-bridges appear to stabilize the parallel-stranded coiled-coil architecture of GCN4-p1 only after the peptide traverses the rate-limiting, dimeric transition state.

Amino Acid Substitution↗
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↗
PubMed2003

Dissecting structural and electrostatic interactions of charged groups in alpha-sarcin. An NMR study of some mutants involving the catalytic residues.

The cytotoxic ribonuclease alpha-sarcin is the best characterized member of the ribotoxin family. Ribotoxins share a common structural core, catalytic residues, and active site topology with members of the broader family of nontoxic microbial extracellular RNases. They are, however, much more specific in their biological action. To shed light on the highly specific alpha-sarcin activity, we have evaluated the structural and electrostatic interactions of its charged groups, by combining the structural and pK(a) characterization by NMR of several variants with theoretical calculations based on the Tanford-Kirkwood and Poisson-Boltzmann models. The NMR data reveal that the global conformation of wild-type alpha-sarcin is preserved in the H50Q, E96Q, H137Q, and H50/137Q variants, and that His137 is involved in an H-bond that is crucial in maintaining the active site structure and in reinforcing the stability of the enzyme. The loss of this H-bond in the H137Q and H50/137Q variants modifies the local structure of the active site. The pK(a) values of active site groups H50, E96, and H137 in the four variants have been determined by two-dimensional NMR. The catalytic dyad of E96 and H137 is not sensitive to charge replacements, since their pK(a) values vary less than +/-0.3 pH unit with respect to those of the wild type. On the contrary, the pK(a) of His50 undergoes drastic changes when compared to its value in the intact protein. These amount to an increase of 0.5 pH unit or a decrease of 1.1 pH units depending on whether a positive or negative charge is substituted at the active site. The main determinants of the pK(a) values of most of the charged groups in alpha-sarcin have been established by considering the NMR results in conjunction with those derived from theoretical pK(a) calculations. With regard to the active site residues, the H50 pK(a) is chiefly influenced by electrostatic interactions with E96 and H137, whereas the effect of the low dielectric constant and the interaction with R121 appear to be the main determinants of the altered pK(a) value of E96 and H137. Charge-charge interactions and an increased level of burial perturb the pK(a) values of the active site residues of alpha-sarcin, which can account for its reduced ribonucleolytic activity and its high specificity.

Amino Acid Substitution↗
PubMed2002

L-phenylalanine binding and domain organization in human phenylalanine hydroxylase: a differential scanning calorimetry study.

Human phenylalanine hydroxylase (hPAH) is a tetrameric enzyme that catalyzes the hydroxylation of L-phenylalanine (L-Phe) to L-tyrosine; a dysfunction of this enzyme causes phenylketonuria. Each subunit in hPAH contains an N-terminal regulatory domain (Ser2-Ser110), a catalytic domain (Asp112-Arg410), and an oligomerization domain (Ser411-Lys452) including dimerization and tetramerization motifs. Two partially overlapping transitions are seen in differential scanning calorimetry (DSC) thermograms for wild-type hPAH in 0.1 M Na-Hepes buffer, 0.1 M NaCl, pH 7.0. Although these transitions are irreversible, studies on their scan-rate dependence support that the equilibrium thermodynamics analysis is permissible in this case. Comparison with the DSC thermograms for truncated forms of the enzyme, studies on the protein and L-Phe concentration effects on the transitions, and structure-energetic calculations based on a modeled structure support that the thermal denaturation of hPAH occurs in three stages: (i) unfolding of the four regulatory domains, which is responsible for the low-temperature calorimetric transition; (ii) unfolding of two (out of the four) catalytic domains, which is responsible for the high-temperature transition; and (iii) irreversible protein denaturation, which is likely responsible for the observed exothermic distortion in the high-temperature side of the high-temperature transition. Stages 1 and 2 do not appear to be two-state processes. We present an approach to the analysis of ligand effects on DSC transition temperatures, which is based on the general binding polynomial formalism and is not restricted to two-state transitions. Application of this approach to the L-Phe effect on the DSC thermograms for hPAH suggests that (i) there are no binding sites for L-Phe in the regulatory domains; therefore, contrary to the common belief, the activation of PAH by L-Phe seems to be the result of its homotropic cooperative binding to the active sites. (ii) The regulatory domain appears to be involved in cooperativity through its interactions with the catalytic and oligomerization domains; thus, upon regulatory domain unfolding, the cooperativity in the binding of L-Phe to the catalytic domains seems to be lost and the value of the L-Phe concentration corresponding to half-saturation is increased. Overall, our results contribute to the understanding of the conformational stability and the substrate-induced cooperative activation of this important enzyme.

Calorimetry, Differential Scanning↗
PubMed2002

Electrostatic interactions in ubiquitin: stabilization of carboxylates by lysine amino groups.

To explore electrostatic interactions in ubiquitin, pK(a) values have been determined by NMR for all 12 carboxyl groups in wild-type ubiquitin and in variants where single lysines have been replaced by neutral residues. Aspartate pK(a) values in ubiquitin range from 3.1 to 3.8 and are generally less than model compound values. Most aspartate pK(a) values are within 0.2 pH unit of those predicted with a simple Tanford-Kirkwood model. Glutamate pK(a) values range from 3.8 to 4.5, close to model compound values and differing by 0.1-0.8 pH unit from calculated values. To determine the role of positive charges in modulating carboxyl pK(a) values, we mutated lysines at positions 11, 29, and 33 to glutamine and threonine. NMR studies with these six single-site mutants reveal significant interactions of Lys 11 and Lys 29 with Glu 34 and Asp 21, respectively: pK(a) values for Glu 34 and Asp 21 increase by approximately 0.5-0.8 pH unit, similar to predicted values, when the lysines are replaced by neutral residues. In contrast, the predicted interaction between Lys 33 and Glu 34 is not observed experimentally. In some instances, substitution of lysine by glutamine and threonine did not lead to the same changes in carboxyl pK(a) values. These may reflect new short-range interactions between the mutated residues and the carboxyl groups. Carboxyl pK(a) shifts > 0.5 pH unit result from mutations at groups that are <5 A from the carboxyl group. No interactions are observed at >10 A.

Amides↗