Search PubMed⌕ Search

Biomedical subjects

B Honig

Publications and source records attributed to B Honig.

At least 55 records · Page 3Linked to original sources

Free energy determinants of secondary structure formation: II. Antiparallel beta-sheets.

The factors that determine the stability of antiparallel beta-sheets are considered via a theoretical analysis of conformational free energies. A series of idealized model polyalanine beta-sheets are built with constraints such that the angular geometry of hydrogen bonding varies in the range observed in proteins while hydrogen bonding distance remains fixed. The conformations of the sheets generated in this way have a broad distribution of twist angles ranging from highly twisted left-handed to highly twisted right-handed orientations. The association free energies of the sheets are calculated with a gas phase CHARMM potential and FDPB/gamma solvation models. Left-handed structures are found to be less stable than right handed structures due to intrachain steric hindrance in isolated left-handed strands. This explains why antiparallel beta-sheets in proteins are invariably twisted in the right-handed direction. The free energy surface for right-handed sheets shows particular preference for conformations ranging from flat to those that exhibit a pronounced right-handed twist. This suggests that antiparallel beta-sheets can adopt a variety of right-handed conformations, a result that is consistent with observations on known proteins. In parallel with our study of alpha-helices we find that van der Waals and hydrophobic interactions are the primary factor stabilizing polyalanine beta-sheets, while electrostatic interactions including hydrogen bonding are found to be destabilizing. However, in contrast to helices, the net change in conformational free energy involving only backbone-backbone interactions (including beta-carbons) is not sufficient to overcome the loss in configurational entropy that accompanies sheet formation. Rather we suggest that cross-strand non-polar side-chain-side-chain interactions are essential for sheet formation, explaining why large non-polar amino acids have the greatest sheet forming propensities. Thus, sheet propensities involve pairwise interactions and are expected to be context dependent, as has been observed in recent experiments.

Amino Acid Sequence↗

Classical electrostatics in biology and chemistry.

A major revival in the use of classical electrostatics as an approach to the study of charged and polar molecules in aqueous solution has been made possible through the development of fast numerical and computational methods to solve the Poisson-Boltzmann equation for solute molecules that have complex shapes and charge distributions. Graphical visualization of the calculated electrostatic potentials generated by proteins and nucleic acids has revealed insights into the role of electrostatic interactions in a wide range of biological phenomena. Classical electrostatics has also proved to be successful quantitative tool yielding accurate descriptions of electrical potentials, diffusion limited processes, pH-dependent properties of proteins, ionic strength-dependent phenomena, and the solvation free energies of organic molecules.

Biochemistry↗

On the magnitude of the electrostatic contribution to ligand-DNA interactions.

A model based on the nonlinear Poisson-Boltzmann equation is used to study the electrostatic contribution to the binding free energy of a simple intercalating ligand, 3,8-diamino-6-phenylphenanthridine, to DNA. We find that the nonlinear Poisson-Boltzmann model accurately describes both the absolute magnitude of the pKa shift of 3,8-diamino-6-phenylphenanthridine observed upon intercalation and its variation with bulk salt concentration. Since the pKa shift is directly related to the total electrostatic binding free energy of the charged and neutral forms of the ligand, the accuracy of the calculations implies that the electrostatic contributions to binding are accurately predicted as well. Based on our results, we have developed a general physical description of the electrostatic contribution to ligand-DNA binding in which the electrostatic binding free energy is described as a balance between the coulombic attraction of a ligand to DNA and the disruption of solvent upon binding. Long-range coulombic forces associated with highly charged nucleic acids provide a strong driving force for the interaction of cationic ligands with DNA. These favorable electrostatic interactions are, however, largely compensated for by unfavorable changes in the solvation of both the ligand and the DNA upon binding. The formation of a ligand-DNA complex removes both charged and polar groups at the binding interface from pure solvent while it displaces salt from around the nucleic acid. As a result, the total electrostatic binding free energy is quite small. Consequently, nonpolar interactions, such as tight packing and hydrophobic forces, must play a significant role in ligand-DNA stability.

DNA↗

Salt effects on polyelectrolyte-ligand binding: comparison of Poisson-Boltzmann, and limiting law/counterion binding models.

The theory for salt dependence of the free energy, entropy, and enthalpy of a polyelectrolyte in the PB (PB) model is extended to treat the nonspecific salt dependence of polyelectrolyte-ligand binding reactions. The salt dependence of the binding constant (K) is given by the difference in osmotic pressure terms between the reactants and products. For simple 1-1 salts it is shown that this treatment is equivalent to the general preferential interaction model for the salt dependence of binding [C. Anderson and M. Record (1993) Journal of Physical Chemistry, Vol. 97, pp. 7116-7126]. The salt dependence, entropy, and enthalpy are compared for the PB model and one specific form of the preferential interaction coefficient model that uses counterion condensation/limiting law (LL) behavior. The PB and LL models are applied to three ligand-polyelectrolyte systems with the same net ligand charge: a model sphere-cylinder binding reaction, a drug-DNA binding reaction, and a protein-DNA binding reaction. For the small ligands both the PB and limiting law models give (In K vs. In[salt]) slopes close in magnitude to the net ligand charge. However, the enthalpy/entropy breakdown of the salt dependence is quite different. In the PB model there are considerable contributions from electrostatic enthalpy and dielectric (water reorientation) entropy, compared to the predominant ion cratic (release) entropy in the limiting law model. The relative contributions of these three terms in the PB model depends on the ligand: For the protein, ion release entropy is the smallest contribution to the salt dependence of binding. The effect of three approximations made in the LL model is examined: These approximations are (1) the ligand behaves ideally, (2) the preferential interaction coefficient of the polyelectrolyte is unchanged upon ligand binding, and (3) the polyelectrolyte preferential interaction coefficient is given by the limiting law/counterion-condensation value. Analysis of the PB model shows that assumptions 2 and 3 break down at finite salt concentrations. For the small ligands the effects on the slope cancel, however, giving net slopes that are similar in the PB and LL models, but with a different entropy/enthalpy breakdown. For the protein ligand the errors from assumptions 2 and 3 in the LL models do not cancel.

Computer Graphics↗

Salt effects on nucleic acids.

Salt-dependent electrostatic effects are a major factor in determining the stability, structure, reactivity, and binding behavior of nucleic acids. Increasingly detailed theoretical methods, especially those based on Monte Carlo and Poisson-Boltzmann methodologies, combined with powerful computational algorithms are being used to examine how the shape, charge distribution and dielectric properties of the molecules affect the ion distribution in the surrounding aqueous solution, and how they play a role in ligand binding, structural transitions and other biologically important reactions. These studies indicate that inclusion of detailed structural information about the nucleic acid and its ligands is crucial for improving models of nucleic acid electrostatics, and that better treatment of the ion atmosphere and dielectric effects is also of major importance.

Algorithms↗

A free energy analysis of nucleic acid base stacking in aqueous solution.

This paper reports a theoretical study of the free energy contributions to nucleic acid base stacking in aqueous solution. Electrostatic interactions are treated by using the finite difference Poisson-Boltzmann method and nonpolar effects are treated with explicit calculation of van der Waals interactions and/or free energy-surface area relationships. Although for some pairs of bases there is a favorable Coulombic interaction in the stacked conformation, generally the net effect of electrostatic interactions is to oppose stacking. This result is caused by the loss of favorable base-solvent electrostatic interactions, that accompany the partial removal of polar atoms from water in the stacked conformation. Nonpolar interactions, involving the hydrophobic effect and enhancement of van der Waals interactions caused by close-packing, drive stacking. The calculations qualitatively reproduce the experimental dependence of stacking free energy on purine-pyrimidine composition.

Adenine↗

Adding backbone to protein folding: why proteins are polypeptides.

It is argued that the chemical nature of the polypeptide backbone is the central determinant of the three-dimensional structures of proteins. The requirement that buried polar groups form intramolecular hydrogen bonds limits the fold of the backbone to the well known units of secondary structure while the amino acid sequence chooses among the set of conformations available to the backbone. 'Sidechain-only' models, based for example on hydrophobicity patterns, fail to account for the properties of the backbone and thus will have difficulty capturing essential features of a folding pathway. This is evident from the incorrect predictions they make for the conformations of the limiting cases of all-hydrophobic or all-polar sequences.

Journal Article↗

An algorithm to generate low-resolution protein tertiary structures from knowledge of secondary structure.

An algorithm is described to assemble the three-dimensional fold of a protein starting from its secondary structure. A reduced representation of the polypeptide chain is used together with a crude potential based on pair hydrophobicities. The method is shown to be successful in locating the native topology for two 4-alpha-helix bundles, myohemerythrin and cytochrome b-562.

Algorithms↗

Salt effects on ligand-DNA binding. Minor groove binding antibiotics.

Salt dependent electrostatic effects play a central role in intermolecular interactions involving nucleic acids. In this paper, the finite-difference solution to the nonlinear Poisson-Boltzmann (NLPB) equation is used to evaluate the salt dependent contribution to the electrostatic binding free energy of the minor groove binding antibiotics DAPI, Hoechst 33258 and netropsin to DNA using detailed molecular structures of the complexes. For each of these systems, a treatment based on the NLPB equation accurately describes the variation of the experimentally observed binding constant with bulk salt concentration. A solvation formalism is developed in which salt effects are described in terms of three free energy contributions: the electrostatic ion-molecule interaction free energy, delta delta G degrees im; the electrostatic ion-ion interaction free energy, delta delta G degrees ii; and the entropic ion organization free energy, delta delta G degrees org. The electrostatic terms, delta delta G degrees im and delta delta G degrees ii, have both enthalpic and entropic components, while the term delta delta G degrees org is purely a cratic entropy. Each of these terms depends significantly on salt dependent changes in the counterion and coion concentrations around the DNA. In each of the systems studied, univalent ions substantially destabilize charged ligand-DNA complexes at physiological salt concentrations. This effect involves a salt dependent redistribution of counterions near the DNA. The free energy associated with the redistribution of counterions upon binding is dominated by the unfavorable change in the electrostatic ion-molecule interactions, delta delta G degrees im, rather than the change in the cratic entropy of ion organization, delta delta G degrees org. In addition, the observed slope of the salt dependence of the free energy is determined by electrostatic ion-molecule and ion-ion interactions as well as the cratic entropy of ion release. These findings are in contrast to models in which the cratic entropy of counterion release drives binding.

Anti-Bacterial Agents↗

Salt effects on protein-DNA interactions. The lambda cI repressor and EcoRI endonuclease.

In this paper, finite-difference solutions to the nonlinear Poisson-Boltzmann (NLPB) equation are used to calculate the salt dependent contribution to the electrostatic DNA binding free energy for both the lambda cI repressor and the EcoRI endonuclease. For the protein-DNA systems studied, the NLPB method describes nonspecific univalent salt dependent effects on the binding free energy which are in excellent agreement with experimental results. In these systems, the contribution of the ion atmosphere to the binding free energy substantially destabilizes the protein-DNA complexes. The magnitude of this effect involves a macromolecular structure dependent redistribution of both cations and anions around the protein and the DNA which is dominated by long range electrostatic interactions. We find that the free energy associated with global ion redistribution upon binding is more important than changes associated with local protein-DNA interactions (ion-pairs) in determining salt effects. The NLPB model reveals how long range salt effects can play a significant role in the relative stability of protein-DNA complexes with different structures.

Bacteriophage lambda↗

Structural origins of pH and ionic strength effects on protein stability. Acid denaturation of sperm whale apomyoglobin.

A recently developed approach to calculate the pH dependence of protein stability from three-dimensional structure information is applied to the analysis of acid denaturation of sperm whale apomyoglobin. The finite difference Poisson-Boltzmann method is used to calculate pKa values and these are used to obtain titration curves for the folded protein as well as for compact intermediates. The total electrostatic free energy change involved in apomyoglobin unfolding is then evaluated. Calculations are carried out of the unfolding free energy of the native (N) and the compact intermediate (I) of apomyoglobin relative to the unfolded state (U) over a range of pH at various ionic strengths. The contributions from key ionizable groups to the unfolding process are discussed. For the acid-induced partial unfolding of apomyoglobin near pH 5, the transition from N to I is found to be driven by three histidines that are exposed when the B, C, D and E helices unfold. Similarly, the unfolding of the compact intermediate I consisting of the A, G and H helices is driven primarily by a few carboxylic acids with low pKa values in the compact state. This picture is in contrast to the view which attributes acid denaturation to electrostatic repulsion resulting from the build up of positive charge. In fact, charge-charge interactions in myoglobin are found to be attractive at all pH values where the protein unfolds. pH-dependent changes in these interactions contribute to acid denaturation but other electrostatic effects, such as hydrogen bonding and solvation, are important as well. The effect of increasing ionic strength on unfolding is attributed to the decrease of attractive charge-charge interactions which destabilize the N state relative to I, but stabilize the I state relative to U by reducing the pKa shifts of a few critical carboxylic acids. The I state is found to be more stable than U at neutral pH thus accounting for its presence as an intermediate on the protein folding pathway. Our results have implications for the origins of compact intermediates or "molten globule" states.

Animals↗

Evaluation of the conformational free energies of loops in proteins.

In this paper we discuss the problem of including solvation free energies in evaluating the relative stabilities of loops in proteins. A conformational search based on a gas-phase potential function is used to generate a large number of trial conformations. As has been found previously, the energy minimization step in this process tends to pack charged and polar side chains against the protein surface, resulting in conformations which are unstable in the aqueous phase. Various solvation models can easily identify such structures. In order to provide a more severe test of solvation models, gas-phase conformations were generated in which side chains were kept extended so as to maximize their interaction with the solvent. The free energies of these conformations were compared to that calculated for the crystal structure in three loops of the protein E. coli RNase H, with lengths of 7, 8, and 9 residues. Free energies were evaluated with a finite difference Poisson-Boltzmann (FDPB) calculation for electrostatics and a surface area-based term for nonpolar contributions. These were added to a gas-phase potential function. A free energy function based on atomic solvation parameters was also tested. Both functions were quite successful in selecting, based on a free energy criterion, conformations quite close to the crystal structure for two of the three loops. For one loop, which is involved in crystal contacts, conformations that are quite different from the crystal structure were also selected. A method to avoid precision problems associated with using the FDPB method to evaluate conformational free energies in proteins is described.

Amino Acid Sequence↗

Correlating solvation free energies and surface tensions of hydrocarbon solutes.

A simple equation relating ratios of transfer free energies and solvation free energies to surface tension is derived. When applied to hydrocarbons in water, experimental values for macroscopic surface tension yield remarkably accurate predictions of a ratio involving microscopic quantities, if one uses transfer free energies which have been adjusted for the effects of solute/solvent volume differences. The results support the validity of applying macroscopic concepts such as interfacial free energy at the molecular level. They further suggest that molecular volume as well as surface area contributes to the solubility of hydrocarbons in water.

Alkanes↗

The electrostatic basis for the interfacial binding of secretory phospholipases A2.

Biochemical and structural data suggest that electrostatic forces play a critical role in the binding of secretory phospholipases A2 to substrate aggregates (micelles, vesicles, monolayers, and membranes). This initial binding (adsorption) of the enzyme to the interface is kinetically distinct from the subsequent binding of substrate to the buried active site. Thus, in the absence of specific active-site interactions, electrostatic forces operating at the molecular surface may orient and hold the enzyme at the interface. We have calculated the electrostatic potentials for 10 species of secretory phospholipases A2 whose atomic coordinates have been determined by x-ray crystallography. Most of these enzymes show a marked electrostatic sidedness that is accentuated to a variable degree by the presence of the essential cofactor calcium ion. This asymmetry suggests a discrete interfacial binding region on the protein's surface, the location of which is in general agreement with proposals derived from the results of chemical modification, mutational, and crystallographic experiments.

Amino Acid Sequence↗

Calculation of electrostatic effects at the amino terminus of an alpha helix.

It is generally believed that the electrostatic field arising from the dipolar charge distribution in alpha helices is important for protein structure and function. We report a calculation of the electrostatic potential and field at the amino terminus of an alpha helix in water, obtained from a finite difference solution to the Poisson-Boltzmann equation. This method takes into account the detailed helix shape and charge distribution, as well as solvent, and generalized ionic strength effects. The calculated potential and field are found to be in good agreement with the experimentally observed helix-induced Stark effect and pKa shifts of a probe at the N-terminus of a stable, monomeric alpha-helical peptide (Lockhart and Kim, 1992, 1993). Ionic screening effects are reproduced at low salt concentrations. Deviations at higher salt concentrations may result from specific ion effects (specific ion-solute and/or ion-solvent interactions). The FDPB method was used to analyze the contributions from each residue, charged side chains, and solvent to the helix potential and field. Backbone contributions come primarily from the first one to two helical turns. Charged side chains contribute to helix-induced pKa shifts for certain probe-peptide combinations, even at relatively large distances from the probe (> 14 A).

Amino Acid Sequence↗

Environmental effects on the protonation states of active site residues in bacteriorhodopsin.

Finite difference solutions of the Poisson-Boltzmann equation are used to calculate the pKa values of the functionally important ionizable groups in bacteriorhodopsin. There are strong charge-charge interactions between the residues in the binding site leading to the possibility of complex titration behavior. Structured water molecules, if they exist in the binding site, can have significant effects on the calculated pKa by strongly stabilizing ionized species. The ionization states of the Schiff base and Asp-85 are found to be strongly coupled. Small environmental changes, which might occur as a consequence of trans-cis isomerization, are capable of causing large shifts in the relative pKa values of these two groups. This provides an explanation for the protonation of Asp-85 and the deprotonation of the Schiff base in the M state of bacteriorhodopsin. The different behavior of Asp-85 and Asp-212 is discussed in this regard.

Aspartic Acid↗

Photochemotherapy beyond psoriasis.

Photochemotherapy involves the therapeutic use of nonionizing radiation in combination with a photosensitizing chemical to trigger a photochemical reaction that mediates a beneficial effect. The successful introduction and widespread use of psoralen photochemotherapy (PUVA) in the management of psoriasis was the chief stimulus for recent interest in the therapeutic use of nonionizing radiation in various other dermatoses. This article discusses the expanding spectrum of diseases responding to PUVA therapy. More than 30 conditions such as atopic dermatitis, mycosis fungoides, vitiligo, the photodermatoses, chronic graft-versus-host disease, and granuloma annulare have been successfully treated with oral psoralen photochemotherapy. Various mechanisms of response to treatment are discussed including photoimmunologic effects, selective cytotoxicity, alterations of cell function, and stimulation of melanocytes. Finally, the limitations to the use of PUVA therapy are identified and its future use in other cutaneous and systemic diseases are discussed.

Humans↗