PubMed Health⌕ Search

Biomedical subjects

B Honig

Publications and source records attributed to B Honig.

At least 19 recordsLinked to original sources
PubMed1992

Analysis of the heat capacity dependence of protein folding.

This paper presents an analysis of plots of enthalpy versus heat capacity change at 25 degrees C for the unfolding of proteins and for the dissolution of gaseous, liquid and solid solutes, first reported by Murphy, Privalov & Gill. The negative slope in the enthalpy plot for proteins is interpreted as arising from a large penalty associated with burying polar groups in the protein interior. The small enthalpy changes that accompany protein unfolding at 25 degrees C are also discussed. It is argued that the combined effects of hydrogen bond formation and close packing predict a large positive enthalpy of unfolding. Electrostatic calculations indicate that the penalty associated with burying polar groups is large enough to effectively cancel these terms, leading to the small net enthalpy changes that are observed. The free energy changes associated with protein folding are also discussed. The free energy cost of burying polar groups largely compensates for the stabilizing contribution of the hydrophobic effect and would appear to account for the fact that proteins are marginally stable, independent of their size and of their relative hydrophobicities.

Hydrogen Bonding↗
PubMed1992

Perturbing the polar environment of Asp102 in trypsin: consequences of replacing conserved Ser214.

Much of the catalytic power of trypsin is derived from the unusual buried, charged side chain of Asp102. A polar cave provides the stabilization for maintaining the buried charge, and it features the conserved amino acid Ser214 adjacent to Asp102. Ser214 has been replaced with Ala, Glu, and Lys in rat anionic trypsin, and the consequences of these changes have been determined. Three-dimensional structures of the Glu and Lys variant trypsins reveal that the new 214 side chains are buried. The 2.2-A crystal structure (R = 0.150) of trypsin S214K shows that Lys214 occupies the position held by Ser214 and a buried water molecule in the buried polar cave. Lys214-N zeta is solvent inaccessible and is less than 5 A from the catalytic Asp102. The side chain of Glu214 (2.8 A, R = 0.168) in trypsin S214E shows two conformations. In the major one, the Glu carboxylate in S214E forms a hydrogen bond with Asp102. Analytical isoelectrofocusing results show that trypsin S214K has a significantly different isoelectric point than trypsin, corresponding to an additional positive charge. The kinetic parameter kcat demonstrates that, compared to trypsin, S214K has 1% of the catalytic activity on a tripeptide amide substrate and S214E is 44% as active. Electrostatic potential calculations provide corroboration of the charge on Lys214 and are consistent with the kinetic results, suggesting that the presence of Lys214 has disturbed the electrostatic potential of Asp102.

Amino Acid Sequence↗
PubMed1992

The electrostatic contribution to DNA base-stacking interactions.

Base-stacking and phosphate-phosphate interactions in B-DNA are studied using the finite difference Poisson-Boltzmann equation. Interaction energies and dielectric constants are calculated and compared to the predictions of simple dielectric models. No extant simple dielectric model adequately describes phosphate-phosphate interactions. Electrostatic effects contribute negligibly to the sequence and conformational dependence of base-stacking interactions. Electrostatic base-stacking interactions can be adequately modeled using the Hingerty screening function. The repulsive and dispersive Lennard-Jones interactions dominate the dependence of the stacking interactions on roll, tilt, twist, and propellor. The Lennard-Jones stacking energy in ideal B-DNA is found to be essentially independent of sequence.

Base Sequence↗
PubMed1991

Electrostatic control of midpoint potentials in the cytochrome subunit of the Rhodopseudomonas viridis reaction center.

The photosynthetic reaction center of Rhodopseudomonas viridis has four hemes in a linear arrangement with alternating high- and low-potential sites. Their midpoints are -60, 20, 310, and 380 mV [Dracheva, S. M., Drachev, L. A., Konstantinov, A. A., Semenov, A. Y., Skulachev, V. P., Arutjunjan, A. M., Shuvalov, V. A. & Zaberezhnaya, S. M. (1988) Eur. J. Biochem. 171, 253-264]. Electrostatic calculations reproduce the 440-mV midpoint spread and assignments of high- and low-potential hemes. When calculations on model compounds to connect the theoretical midpoints to the standard hydrogen electrode are used, the absolute electrochemical midpoints for the reaction center hemes are also in good agreement with experiment. The free energy of oxidation is found to be dependent on pairwise interactions with charged amino acids, heme propionic acids, previously oxidized hemes, and axial ligands and on the reaction field induced by heme oxidation.

Amino Acid Sequence↗
PubMed1991

Extracting hydrophobic free energies from experimental data: relationship to protein folding and theoretical models.

Solubility and vapor pressure measurements of hydrocarbons in water are generally thought to provide estimates of the strength of the hydrophobic effect in the range 20-30 cal/(mol.A2). Our reassessment of the solubility data on the basis of new developments in solution thermodynamics suggests that the hydrophobic surface free energy for hydrocarbon solutes is 46-47 cal/(mol.A2), although the actual value depends strongly on curvature effects [Nicholls et al. (1991) Proteins (in press); Sharp et al. (1991) Science 252, 106-109]. The arguments to support such a significant increase in the estimate of the hydrophobic effect stem partly from theoretical considerations and partly from the experimental results of De Young and Dill [(1990) J. Phys. Chem. 94, 801-809] on benzene partition between water and alkane solvents. Previous estimates of the hydrophobic effect derive from an analysis of solute partition data, which does not fully account for changes in volume entropy. We show here how the ideal gas equations, combined with experimental molar volumes, can account for such changes. Revised solubility scales for the 20 amino acids, based on cyclohexane to water and octanol to water transfer energies, are derived. The agreement between these scales, particularly the octanol scale, and mutant protein stability measurements from Kellis et al. [(1989) Biochemistry 28, 4914-4922] and Shortle et al. [(1990) Biochemistry 29, 8033-8041] is good. The increased strength of the hydrophobic interaction has implications for the energetics of protein folding, substrate binding, and nucleic acid base stacking and the interpretation of computer simulations.

Alkanes↗
PubMed1991

Reconciling the magnitude of the microscopic and macroscopic hydrophobic effects.

The magnitude of the hydrophobic effect, as measured from the surface area dependence of the solubilities of hydrocarbons in water, is generally thought to be about 25 calories per mole per square angstrom (cal mol-1 A-2). However, the surface tension at a hydrocarbon-water interface, which is a "macroscopic" measure of the hydrophobic effect, is approximately 72 cal mol-1 A-2. In an attempt to reconcile these values, alkane solubility data have been reevaluated to account for solute-solvent size differences, leading to a revised "microscopic" hydrophobic effect of 47 cal mol-1 A-2. This value, when used in a simple geometric model for the curvature dependence of the hydrophobic effect, predicts a macroscopic alkane-water surface tension that is close to the macroscopic value.

Hydrocarbons↗
PubMed1991

Protein folding and association: insights from the interfacial and thermodynamic properties of hydrocarbons.

We demonstrate in this work that the surface tension, water-organic solvent, transfer-free energies and the thermodynamics of melting of linear alkanes provide fundamental insights into the nonpolar driving forces for protein folding and protein binding reactions. We first develop a model for the curvature dependence of the hydrophobic effect and find that the macroscopic concept of interfacial free energy is applicable at the molecular level. Application of a well-known relationship involving surface tension and adhesion energies reveals that dispersion forces play little or no net role in hydrophobic interactions; rather, the standard model of disruption of water structure (entropically driven at 25 degrees C) is correct. The hydrophobic interaction is found, in agreement with the classical picture, to provide a major driving force for protein folding. Analysis of the melting behavior of hydrocarbons reveals that close packing of the protein interior makes only a small free energy contribution to folding because the enthalpic gain resulting from increased dispersion interactions (relative to the liquid) is countered by the freezing of side chain motion. The identical effect should occur in association reactions, which may provide an enormous simplification in the evaluation of binding energies. Protein binding reactions, even between nearly planar or concave/convex interfaces, are found to have effective hydrophobicities considerably smaller than the prediction based on macroscopic surface tension. This is due to the formation of a concave collar region that usually accompanies complex formation. This effect may preclude the formation of complexes between convex surfaces.

Algorithms↗
PubMed1991

The inclusion of electrostatic hydration energies in molecular mechanics calculations.

The problem of including solvent effects in molecular mechanics calculations is discussed. It is argued that the neglect of charge-solvent (solvation) interactions can introduce significant errors. The finite difference Poisson-Boltzmann (FDPB) method for calculating electrostatic interactions is summarized and is used as a basis for introducing a new pairwise energy term which accounts for charge-solvent interactions. This term acts between all pairs of atoms usually considered in molecular mechanics calculations and can be easily incorporated into existing force fields. As an example, a parameterization is developed for the CHARMm force field and the results compared to the predictions of the FDPB method. An approach to the realistic incorporation of solvent screening into force fields is also outlined.

Chemical Phenomena↗
PubMed1990

Electrical potential of transfer RNAs: codon-anticodon recognition.

Calculations of the electrostatic potentials were made around yeast elongator phenylalanine, aspartate tRNAs, and yeast initiator methionine tRNA in aqueous solution at physiological ionic strength. The calculations were carried out with a finite difference algorithm for solving the nonlinear Poisson-Boltzmann equation that incorporates the screening effects of the electrolyte, the exclusion of ions by the molecule, the molecular shape, and the different polarizabilities of the solvent and the tRNA. The initiator tRNA is surrounded by uniformly spaced contours of negative potential. The elongator tRNAs are also surrounded by a similar contour pattern except in the anticodon region where there is a pronounced "hole" in the potential surface. This hole is caused by an invagination of the potential contours, which also results in an increase in the local field strength. The effect of this hole is that the anticodon region in the elongator tRNAs is the least negative, or conversely the most positive, region of the molecule. This hole, which is not found when simple Coulombic potentials are used, is due both to the structure of the elongator tRNA anticodon loops and to the different polarizabilities of the solvent and tRNA. The existence of the potential hole in elongator tRNAs may account in part for their ability to associate with other negatively charged macromolecules, in particular mRNA. Moreover, it suggests that the anticodon loop of elongator tRNAs is the energetically most favorable point of approach of mRNA to tRNA.

Anticodon↗
PubMed1990

Octopus photoreceptor membranes. Surface charge density and pK of the Schiff base of the pigments.

The chromophore of octopus rhodopsin is 11-cis retinal, linked via a protonated Schiff base to the protein backbone. Its stable photoproduct, metarhodopsin, has all-trans retinal as its chromphore. The Schiff base of acid metarhodopsin (lambda max = 510 nm) is protonated, whereas that of alkaline metarhodopsin (lambda max = 376 nm) is unprotonated. Metarhodopsin in photoreceptor membranes was titrated and the apparent pK of the Schiff base was measured at different ionic strengths. From these salt-dependent pKs the surface charge density of the octopus photoreceptor membranes and the intrinsic Schiff base pK of metarhodopsin were obtained. The surface charge density is sigma = -1.6 +/- 0.1 electronic charges per 1,000 A2. Comparison of the measured surface charge density with values from octopus rhodopsin model structures suggests that the measured value is for the extracellular surface and so the Schiff base in metarhodopsin is freely accessible to protons from the extracellular side of the membrane. The intrinsic Schiff base pK of metarhodopsin is 8.44 +/- 0.12, whereas that of rhodopsin is found to be 10.65 +/- 0.10 in 4.0 M KCl. These pK values are significantly higher than the pK value around 7.0 for a retinal Schiff base in a polar solvent; we suggest that a plausible mechanism to increase the pK of the retinal pigments is the preorganization of their chromophore-binding sites. The preorganized site stabilizes the protonated Schiff base with respect to the unprotonated one. The difference in the pK for the octopus rhodopsin compared with metarhodopsin is attributed to the relative freedom of the latter's chromophore-binding site to rearrange itself after deprotonation of the Schiff base.

Amino Acid Sequence↗
PubMed1989

Electrical potentials in trypsin isozymes.

Cow and rat trypsin differ in net charge by 12.5 units yet have the same enzymatic mechanism. The role of electrical potentials in the catalytic mechanism of these trypsin isozymes is investigated by using the finite difference Poisson-Boltzmann method. The calculations reveal that the active sites are effectively shielded from surface charge, thus making it possible for the two enzymes to have essentially identical potentials in their catalytically important regions. The potentials in both active sites are dominated by local interactions arising both from partial charges and from the negative charge on Asp-102. The latter is found to stabilize the transition state by about 4 kcal/mol, a value that is consistent with the extent of reduced catalytic activity in the variant Asn-102 trypsin, in which the negative charge is absent. The calculations predict that Asp-102 is ionized and that His-57 is neutral in the resting state of the enzyme. In contrast to their negligible effect on catalytic activity, the cumulative effect of surface charges is found to raise the pK of the N-terminal alpha-amino group of Ile-16 in the rat enzyme by about 1.5 units relative to that of cow trypsin. This charged amino acid forms an ionic bond with Asp-194, which stabilizes the active conformation of the enzyme. An increase in pK of Ile-16 thus provides a possible explanation for the retention of activity of rat trypsin at high pH. The results of this study could not have been obtained from an electrostatic model based on Coulombic potentials.

Amino Acids↗
PubMed1989

The electrostatic potential of B-DNA.

Electrostatic potentials around DNA are obtained by solving the nonlinear Poisson-Boltzmann (PB) equation. The detailed charge distribution of the DNA and the different polarizabilities of the macromolecule and solvent are included explicitly in the calculations. The PB equation is solved using extensions of a finite difference approach applied previously to proteins. Electrical potentials and ion concentrations are compared to those obtained with simpler models. It is found that the shape of the dielectric boundary between the macromolecule and solvent has significant effects on the calculated potentials near the surface, particularly in the grooves. Sequence-specific patterns are found, the most surprising result being the existence of positive regions of potential near the bases in both the major and minor grooves. The effect of solvent and ionic atmosphere screening of phosphate-phosphate repulsions is studied, and an effective dielectric function, appropriate for molecular mechanics simulations, is derived.

DNA↗
PubMed1989

Wavelength regulation in iodopsin, a cone pigment.

The opsin shift, the difference in wavenumber between the absorption peak of a visual pigment and the protonated Schiff base of the chromophore, represents the influence of the opsin binding site on the chromophore. The opsin shift for the chicken cone pigment iodopsin is much larger than that for rhodopsin. To understand the origin of this opsin shift and the mechanism of wavelength regulation in iodopsin, a series of synthetic 9-cis and 11-cis dehydro- and dihydro-retinals was used to regenerate iodopsin-based pigments. The opsin shifts of these pigments are quite similar to those found in bacteriorhodopsin-based artificial pigments. On the basis of these studies, a tentative model of wavelength regulation in iodopsin is proposed.

Animals↗
PubMed1989

Destabilization of an alpha-helix-bundle protein by helix dipoles.

The finite difference Poisson-Boltzmann method is used to calculate the electrostatic work of assembling the four alpha-helices of Themiste dyscritum hemerythrin to form the protein's observed antiparallel helical bundle. The calculations account for the interaction of each helix dipole with the high-dielectric solvent as well as for pairwise interactions of the dipoles with each other. We find that the electrostatic work of assembly is dominated by unfavorable changes in dipole-solvent interactions rather than by favorable interactions between antiparallel helices. Furthermore, the electrostatic energy difference between the observed arrangement of helices in hemerythrin and at least one other possible helical arrangement is less than 1 kT. These results suggest that the helix dipole actually destabilizes the helical bundle and that it plays little or no role in producing the observed bundle geometry.

Calorimetry↗
PubMed1988

Calculation of the total electrostatic energy of a macromolecular system: solvation energies, binding energies, and conformational analysis.

In this report we describe an accurate numerical method for calculating the total electrostatic energy of molecules of arbitrary shape and charge distribution, accounting for both Coulombic and solvent polarization terms. In addition to the solvation energies of individual molecules, the method can be used to calculate the electrostatic energy associated with conformational changes in proteins as well as changes in solvation energy that accompany the binding of charged substrates. The validity of the method is examined by calculating the hydration energies of acetate, methyl ammonium, ammonium, and methanol. The method is then used to study the relationship between the depth of a charge within a protein and its interaction with the solvent. Calculations of the relative electrostatic energies of crystal and misfolded conformations of Themiste dyscritum hemerythrin and the VL domain of an antibody are also presented. The results indicate that electrostatic charge-solvent interactions strongly favor the crystal structures. More generally, it is found that charge-solvent interactions, which are frequently neglected in protein structure analysis, can make large contributions to the total energy of a macromolecular system.

Binding, Competitive↗
PubMed1988

Chromophore/protein and chromophore/anion interactions in halorhodopsin.

Halorhodopsin (HR), the light-driven chloride transport pigment of Halobacterium halobium, was bleached and reconstituted with retinal analogues with the pi electron system interrupted at different locations (dihydroretinals). The absorption maxima of the artificial pigments formed with the dihydroretinals are found to be very similar to those of the corresponding pigments formed by reconstitution of bacteriorhodopsin (BR) and sensory rhodopsin (SR). This strongly suggests that the distribution of charges around the retinal is similar in all three bacterial rhodopsins. Comparison of the primary, and proposed secondary, structures for HR and BR reveal conserved asparagine (asp) and arginine (arg) residues, which are likely candidates for the ionizable amino acids that interact with the retinal. In a second set of experiments absorption shifts due to the binding of anions to Sites I and II in HR, reconstituted with different retinal analogues, were used to estimate the locations of these binding sites relative to the retinal. Site I is localized near the Schiff base, and Site II near the ionone ring. On the basis of these results a structural model for HR is proposed, which accounts for the spectroscopic properties of HR in terms of the three buried arg residues and two of the buried asp residues in the protein.

Bacteriorhodopsins↗