Search PubMed⌕ Search

Biomedical subjects

A Warshel

Publications and source records attributed to A Warshel.

At least 55 records · Page 3Linked to original sources
PubMed1991

Dipoles localized at helix termini of proteins stabilize charges.

The charge-stabilization effect associated with alpha-helices in proteins has been reexamined by microscopic calculations without any a priori assumptions about the dielectric constant of the protein. The calculations reproduce the observed charge stabilization effect of a helix in two well-defined test cases: a histidine residue situated at the C-terminal end of a helix in barnase and the sulfate ligand located near the N-terminal end of a helix of the sulfate-binding protein. They also show that the effective dielectric constant for helix-charge interactions is much larger than previously assumed and that the stabilizing effect of the helix is not associated with the helix macrodipole but rather with a few localized dipoles confined mostly to the first turn of the helix. It is predicted that mutations at one end of the helix should have very small effects on the stabilization of charges at the opposite terminus. It is pointed out that the relatively short-ranged effect of the helix is essentially similar to other cases in which localized dipoles play key roles in electrostatic stabilization.

Bacterial Proteins↗
PubMed1990

Electrostatic control of charge separation in bacterial photosynthesis.

Electrostatic interaction energies of the electron carriers with their surroundings in a photosynthetic bacterial reaction center are calculated. The calculations are based on the detailed crystal structure of reaction centers from Rhodopseu-domonas viridis, and use an iterative, self-consistent procedure to evaluate the effects of induced dipoles in the protein and the surrounding membrane. To obtain the free energies of radical-pair states, the calculated electrostatic interaction energies are combined with the experimentally measured midpoint redox potentials of the electron carriers and of bacteriochlorophyll (BChl) and bacteriopheophytin (BPh) in vitro. The P+HL- radical-pair, in which an electron has moved from the primary electron donor (P) to a BPh on the 'L' side of the reaction center (HL), is found to lie approx. 2.0 kcal/mol below the lowest excited singlet state (P*), when the radical-pair is formed in the static crystallographic structure. The reorganization energy for the subsequent relaxation of P+HL- is calculated to be 5.0 kcal/mol, so that the relaxed radical-pair lies about 7 kcal/mol below P*. The unrelaxed P+BL- radical-pair, in which the electron acceptor is the accessory BChl located between P and HL, appears to be essentially isoenergetic with P*.P+BM-, in which an electron moves to the BChl on the 'M' side, is calculated to lie about 5.5 kcal/mol above P*. These results have an estimated error range of +/- 2.5 kcal/mol. They are shown to be relatively insensitive to various details of the model, including the charge distribution in P+, the atomic charges used for the amino acid residues, the boundaries of the structural region that is considered microscopically and the treatments of the histidyl ligands of P and of potentially ionizable amino acids. The calculated free energies are consistent with rapid electron transfer from P* to HL by way of BL, and with a much slower electron transfer to the pigments on the M side. Tyrosine M208 appears to play a particularly important role in lowering the energy of P+BL-. Electrostatic interactions with the protein favor localization of the positive charge of P+ on PM, one of the two BChl molecules that make up the electron donor.

Bacterial Proteins↗
PubMed1989

Dispersed polaron simulations of electron transfer in photosynthetic reaction centers.

A microscopic method for simulating quantum mechanical, nuclear tunneling effects in biological electron transfer reactions is presented and applied to several electron transfer steps in photosynthetic bacterial reaction centers. In this "dispersed polaron" method the fluctuations of the protein and the electron carriers are projected as effective normal modes onto an appropriate reaction coordinate and used to evaluate the quantum mechanical rate constant. The simulations, based on the crystallographic structure of the reaction center from Rhodopseudomonas viridis, focus on electron transfer from a bacteriopheophytin to a quinone and the subsequent back-reaction. The rates of both of these reactions are almost independent of temperature or even increase with decreasing temperature. The simulations reproduce this unusual temperature dependence in a qualitative way, without the use of adjustable parameters for the protein's Franck-Condon factors. The observed dependence of the back-reaction on the free energy of the reaction also is reproduced, including the special behavior in the "inverted region."

Bacterial Proteins↗
PubMed1989

Calculations of free energy profiles for the staphylococcal nuclease catalyzed reaction.

Calculations of the free energy profile for the first two (rate-limiting) steps of the staphylococcal nuclease catalyzed reaction are reported. The calculations are based on the empirical valence bond method in combination with free energy perturbation molecular dynamics simulations. The calculated activation free energy is in good agreement with experimental kinetic data, and the catalytic effect of the enzyme is reproduced without any arbitrary adjustment of parameters. The enormous reduction of the activation barrier (relative to the reference reaction in water) appears to be largely associated with the strong electrostatic effect of the Ca2+ ion and the two arginine residues in the active site. This favorable electrostatic environment reduces the cost of the general-base catalysis step by almost 15 kcal/mol (by stabilizing the OH- nucleophile) and then stabilizes the developing negative charge on the 5'-phosphate group in the second step of the reaction by about 19 kcal/mol. The basic features of the originally postulated enzyme mechanism (Cotton et al., 1979) are found to be compatible with the observed activation free energy. However, the proposed modification of the mechanism (Sepersu et al., 1987), in which Arg 87 interacts only with the pentacoordinated transition state, is supported by the simulations. Further calculations on the D21E mutant also give results in good agreement with kinetic data.

Binding Sites↗
PubMed1989

How do serine proteases really work?

Recent advances in genetic engineering have led to a growing acceptance of the fact that enzymes work like other catalysts by reducing the activation barriers of the corresponding reactions. However, the key question about the action of enzymes is not related to the fact that they stabilize transition states but to the question to how they accomplish this task. This work considers the catalytic reaction of serine proteases and demonstrates how one can use a combination of calculations and experimental information to elucidate the key contributions to the catalytic free energy. Recent reports about genetic modifications of the buried aspartic group in serine proteases, which established the large effect of this group (but could not determine its origin), are analyzed. Two independent methods indicate that the buried aspartic group in serine proteases stabilizes the transition state by electrostatic interactions rather than by alternative mechanisms. Simple free energy considerations are used to eliminate the double proton-transfer mechanism (which is depicted in many textbooks as the key catalytic factor in serine proteases). The electrostatic stabilization of the oxyanion side of the transition state is also considered. It is argued that serine proteases and other enzymes work by providing electrostatic complementarity to the changes in charge distribution occurring during the reactions they catalyze.

Kinetics↗
PubMed1989

Role of arginine-38 in regulation of the cytochrome c oxidation-reduction equilibrium.

Arg-38 is an internal residue of mitochondrial cytochrome c that is close to heme propionate-7. Previous work comparing the behavior of cytochromes c from several species [Moore, G. R., Harris, D. E., Leitch, F. A., & Pettigrew, G. W. (1984) Biochim. Biophys. Acta 764, 331-342] has suggested that Arg-38 lowers the pKa of this propionate group and thereby accounts for the relative pH independence of the cytochrome c reduction potential from pH 5 to pH 8. The influence of Arg-38 on the oxidation-reduction equilibrium of yeast iso-1-cytochrome c has now been investigated by electrochemical, NMR, and theoretical analysis of six specifically mutated forms of this protein in which Arg has been replaced by Lys, His, Gln, Asn, Leu, or Ala. As the electron-withdrawing character of the residue at position 38 decreases, the reduction potential of the protein also decreases, with the largest decrease (ca. 50 mV) observed for the Ala variant. However, the variation in the reduction potentials of the mutants as a function of pH was similar to that observed for the wild-type protein. The effects of some of these mutations on the pKa values of His-33 and His-39 have been determined by NMR spectroscopy and found to be minimal. Calculations of the electrostatic free energy for the Leu-38 variant predict a decrease in the reduction potential of this mutant that is remarkably close to that observed experimentally. This work establishes that while Arg-38 contributes to the relatively high reduction potential of cytochrome c, this residue does not appear to be the sole functionality responsible for lowering the heme propionate-7 pKa.

Amino Acids↗
PubMed1989

Energetics of ion permeation through membrane channels. Solvation of Na+ by gramicidin A.

Calculations of the solvation energetics for a Na+ ion inside the Gramicidin A channel and in water are presented. The protein dipoles Langevin dipoles (PDLD) method is used to obtain an electrostatic free energy profile for ion permeation through the channel. To gauge the quality of the PDLD results the solvation free energy of a Na+ ion in water and in the center of the channel is also calculated using free energy perturbation (FEP) simulations. The effect of the polarisability of the surrounding lipid membrane is taken into account by representing the membrane by a large grid of polarisable point dipoles. The two methods give similar solvation energies in the interior of the channel and these are less than 5 kcal/mol above the solvation free energy for Na+ in water, in good agreement with experimental data on the activation barriers for ion permeation. It appears that the problems associated with previous calculations of energy profiles in membrane channels can be overcome by a consistent treatment of all the relevant electrostatic contributions. In particular, we find that the induced dipoles of the membrane and the protein contributes with approximately 10 kcal/mol to the solvation energy inside the channel and can therefore not be discarded in a realistic description of ion solvation in the Gramicidin channel.

Gramicidin↗
PubMed1989

Enzymes work by solvation substitution rather than by desolvation.

Considerable attention has recently been drawn to the hypothesis that enzymes catalyze their reactions by displacing solvent and creating an environment similar to the gas phase for the reacting substrates. This "desolvation hypothesis" is reexamined in this paper by defining a common reference energy for reactions in various environments. It is argued that consistent attempts to describe the actual energetics of enzymatic reactions, taking either gas phase or solution as a reference, would contradict the above hypothesis. That is, the enzyme does remove water molecules from its substrate, but substitutes these molecules for another polar environment (namely, its active site). By taking amide hydrolysis as an example, we use experimentally estimated solvation energies and analyze the reaction profile in the gas phase, in solution, and in enzyme active sites. We show that the gas-phase reaction is characterized by an enormous activation barrier (associated with forming the charged nucleophile from neutral fragments), although the nucleophilic attack is essentially barrierless. On the other hand, the enzyme and solution reactions are found to have similar reaction profiles, with a lower activation barrier for the enzymatic reaction. Presumably, the fact that previous analyses of this problem did not involve the construction of the relevant thermodynamic cycles (and quantitative estimates of the corresponding solvation energies) led to the desolvation hypothesis. Our conclusion is that enzyme active sites provide specific polar environments that do not resemble the gas phase but that are designed for electrostatic stabilization of ionic transition states and that "solvate" these states more than water does.

Binding Sites↗
PubMed1988

Why ion pair reversal by protein engineering is unlikely to succeed.

Genetic engineering is a powerful tool for exploring correlations between structure and function in proteins, but as yet we are unable to use it for effective protein design. One of the most interesting examples, which would seem to be obvious, is reversing the polarity of an ion pair. Changing a positively charged protein group, that provides a strong binding for negative substrates, to a negative group is expected to provide an effective binding site for a positively charged substrate. But several recent experiments on aspartate aminotransferase, trypsin and aspartate transcarbamoylase (Schachman, H. K. personal communication) have indicated that polarity reversal is not so successful. Here we argue that the same factors that make the enzyme an effective system for the (-+) pair will make it a much less effective system for the (+-) pair. We also point out that the unusually low effective dielectric constant (epsilon approximately equal to 13) for the (-+) interaction is due to its microenvironment and this will destabilize a (+-) arrangement having an entirely different dielectric constant (epsilon approximately equal to 80). The calculations presented here evaluate the energetics of ion pairs in protein active sites on a semiquantitative level. This is particularly important when dealing with strong, functionally important interactions that are difficult to evaluate with macroscopic models.

Aspartate Aminotransferases↗
PubMed1988

Evaluation of catalytic free energies in genetically modified proteins.

A combination of the empirical valence bond method and a free energy perturbation approach is used to simulate the activity of genetically modified enzymes. The simulations reproduce in a semiquantitative way the observed effects of mutations on the activity and binding free energies of trypsin and subtilisin. This suggests that we are approaching a stage of quantitative structure-function correlation of enzymes. The analysis of the calculations points towards the electrostatic energy of the reacting system as the key factor in enzyme catalysis. The changes in the charges of the reacting system and the corresponding changes in "solvation" free energy (generalized here as the interaction between the charges and the given microenvironment) are emphasized. It is argued that a reliable evaluation of these changes might be sufficient for correlating structure and catalysis. The use of free energy perturbation methods and thermodynamic cycles for evaluation of solvation energies and reactivity is discussed, pointing out our early contributions. The apparent elaborated nature of our treatment is clarified, explaining that such a treatment is essential for consistent calculations of chemical reactions in polar environments. The problems associated with seemingly more rigorous quantum mechanical methods are discussed, emphasizing the inconsistency associated with using gas phase charge distributions. The importance of dynamic aspects is examined by evaluating the autocorrelation of the protein "reaction field" on the reacting substrate. It is found that, at least in the present case, dynamic effects are not important. The nature of the catalytic free energy is considered, arguing that the protein provides preoriented dipoles (polarized to stabilize the transition state charge distribution) and small reorganization energy, thus reducing the activation free energy. The corresponding catalytic free energy is related to the folding free energy, which is being invested in aligning the active site dipoles.

Amino Acid Sequence↗
PubMed1987

Semiquantitative calculations of catalytic free energies in genetically modified enzymes.

The catalytic free energy and binding free energies of the native and the Asn-155----Thr, Asn-155----Leu, and Asn-155----Ala mutants of subtilisin are calculated by the empirical valence bond method and a free energy perturbation method. Two simple procedures are used; one "mutates" the substrate, and the other "mutates" the enzyme. The calculated changes in free energies (delta delta G not equal to cat and delta delta Gbind) between the mutant and native enzymes are within 1 kcal/mol of the corresponding observed values. This indicates that we are approaching a quantitative structure-function correlation. The calculated changes in catalytic free energies are almost entirely due to the electrostatic interaction between the enzyme-water system and the charges of the reacting system. This supports the idea that the electrostatic free energy associated with the changes of charges of the reacting system is the key factor in enzyme catalysis.

Amino Acid Sequence↗
PubMed1986

Free energy of charges in solvated proteins: microscopic calculations using a reversible charging process.

Evaluation of the free energy of ionization of acidic groups in proteins may be used as a powerful and general test case for determining the reliability of calculations of electrostatic energies in macromolecules. This work attacks this test case by using an adiabatic charging process that evaluates the changes in free energies associated with ionizing the acidic groups Asp-3 and Glu-7 in bovine pancreatic trypsin inhibitor and aspartic acid in solution. The results of these free energy calculations are very encouraging; the error range is about 1 kcal/mol for these free energy changes of about-70 kcal/mol. This indicates that we are finally approaching the stage of obtaining quantitative results in modeling the energetics of solvated proteins.

Aprotinin↗
PubMed1986

Control of the redox potential of cytochrome c and microscopic dielectric effects in proteins.

X-ray structural information provides the opportunity to explore quantitatively the relation between the microenvironments of heme proteins and their redox potentials. This can be done by considering the protein as a "solvent" for its redox center and calculating the difference between the electrostatic energy of the reduced and oxidized heme. Such calculations are presented here, applying the protein dipoles-Langevin dipoles (PDLD) model to cytochrome c. The calculations focus on an evaluation of the difference between the redox potentials of cytochrome c and the octapeptide-methionine complex formed by hydrolysis of cytochrome c. The corresponding difference (approximately 7 kcal/mol) is accounted for by the PDLD calculations. It is found that the protein provides basically a low dielectric environment for the heme, which destabilizes the oxidized heme (relative to its energy in water). The effect of the charged propionic acids on the heme is examined in a preliminary way. It is found that the negative charges of these groups are in a hydrophilic rather than a hydrophobic environment and that the protein-water system provides an effective high dielectric constant for their interaction with the heme. The dual nature of the dielectric effect of the cytochrome (a low dielectric constant for the self-energy of the heme and a high dielectric constant for charge-charge interactions) is discussed. The findings of this work are consistent with the difference between the folding energies of the reduced and oxidized cytochrome c.

Cytochrome c Group↗
PubMed1986

Toward computer-aided site-directed mutagenesis of enzymes.

A preliminary attempt to simulate the observed effect of a site-directed mutagenesis of rat trypsin gives encouraging results. The calculations reproduce in a semiquantitative way the observed change in the activation barrier of the rate-limiting step of amide hydrolysis. This result, which did not require any adjustable parameters, indicates that our method may provide a reliable basis for computer-aided enzyme design. In addition to the potentially practical value of the calculations, they provide important mechanistic information--that is, the change in the catalytic effect in trypsin appears to be almost exclusively due to the change in the electrostatic stabilization of the ionic configurations. This supports the view that electrostatic effects are the major factor in enzyme catalysis.

Animals↗
PubMed1985

Calculations of electrostatic energies in proteins. The energetics of ionized groups in bovine pancreatic trypsin inhibitor.

The importance of including different energy contributions in calculations of electrostatic energies in proteins is examined by calculating the intrinsic pKa values of the acidic groups of bovine pancreatic trypsin inhibitor. It appears that such calculations provide a powerful and revealing test; the relevant solvation energies of the ionized acids are of the order of -70 kcal/mol (1 cal = 4.184 J), and microscopic calculations that do not attempt to simulate the complete protein dielectric effect (including the surrounding solvent) can underestimate the solvation energy by as much as 50 kcal/mol. Reproducing correctly, by the same set of parameters, the solvation energies of ionized acids in different sites of a protein cannot be accomplished by including only part of the key energy contributions. The problems associated with macroscopic calculations are also considered and illustrated by the specific case of bovine pancreatic trypsin inhibitor. A promising approach is shown to be provided by a refinement of the previously developed Protein Dipoles Langevin Dipoles model. This model seems to represent consistently the microscopic dielectric of the protein and the surrounding water molecules. The model overcomes the problems associated with the macroscopic models (by treating explicitly the solvent molecules) and avoids the convergence problems associated with all-atom solvent models (by treating the average solvent polarization rather than averaging the actual polarization energy). This paper describes in detail the actual implementation of the model and examines its performance in evaluating intrinsic pKa values. Preliminary microscopic considerations of charge-charge interactions are presented.

Animals↗