Search PubMed⌕ Search

Biomedical subjects

B Honig

Publications and source records attributed to B Honig.

At least 37 records · Page 2Linked to original sources
PubMed1999

Structural determinants of trypsin affinity and specificity for cationic inhibitors.

The binding free energies of four inhibitors to bovine beta-trypsin are calculated. The inhibitors use either ornithine, lysine, or arginine to bind to the S1 specificity site. The electrostatic contribution to binding free energy is calculated by solving the finite difference Poisson-Boltzmann equation, the contribution of nonpolar interactions is calculated using a free energy-surface area relationship and the loss of conformational entropy is estimated both for trypsin and ligand side chains. Binding free energy values are of a reasonable magnitude and the relative affinity of the four inhibitors for trypsin is correctly predicted. Electrostatic interactions are found to oppose binding in all cases. However, in the case of ornithine- and lysine-based inhibitors, the salt bridge formed between their charged group and the partially buried carboxylate of Asp189 is found to stabilize the complex. Our analysis reveals how the molecular architecture of the trypsin binding site results in highly specific recognition of substrates and inhibitors. Specifically, partially burying Asp189 in the inhibitor-free enzyme decreases the penalty for desolvation of this group upon complexation. Water molecules trapped in the binding interface further stabilize the buried ion pair, resulting in a favorable electrostatic contribution of the ion pair formed with ornithine and lysine side chains. Moreover, all side chains that form the trypsin specificity site are partially buried, and hence, relatively immobile in the inhibitor-free state, thus reducing the entropic cost of complexation. The implications of the results for the general problem of recognition and binding are considered. A novel finding in this regard is that like charged molecules can have electrostatic contributions to binding that are more favorable than oppositely charged molecules due to enhanced interactions with the solvent in the highly charged complex that is formed.

Animals↗
PubMed1999

GRASS: a server for the graphical representation and analysis of structures.

GRASS (Graphical Representation and Analysis of Structures Server), a new web-based server, is described. GRASS exploits many of the features of the GRASP program and is designed to provide interactive molecular graphics and quantitative analysis tools with a simple interface over the World-Wide Web. Using GRASS, it is now possible to view many surface features of biological macromolecules on either standard workstations used in macromolecular analysis or personal computers. The result is a World-Wide Web-based, platform-independent, easily used tool for macromolecular visualization and structure analysis.

Computer Graphics↗
PubMed1998

Electrostatics and the membrane association of Src: theory and experiment.

The binding of Src to phospholipid membranes requires both hydrophobic insertion of its myristate into the hydrocarbon interior of the membrane and nonspecific electrostatic interaction of its N-terminal cluster of basic residues with acidic phospholipids. We provide a theoretical description of the electrostatic partitioning of Src onto phospholipid membranes. Specifically, we use molecular models to represent a nonmyristoylated peptide corresponding to residues 2-19 of Src [nonmyr-Src(2-19); GSSKSKPKDPSQRRRSLE-NH2] and a phospholipid bilayer, calculate the electrostatic interaction by solving the nonlinear Poisson-Boltzmann equation, and predict the molar partition coefficient using statistical thermodynamics. The theoretical predictions agree with experimental data obtained by measuring the partitioning of nonmyr-Src(2-19) onto phospholipid vesicles: membrane binding increases as the mole percent of acidic lipid in the vesicles is increased, the ionic strength of the solution is decreased, or the net positive charge of the peptide is increased. The theoretical model also correctly predicts the measured partitioning of the myristoylated peptide, myr-Src(2-19); for example, adding 33% acidic lipid to electrically neutral vesicles increases the partitioning of myr-Src(2-19) 100-fold. Phosphorylating either serine 12 (by protein kinase C) or serine 17 (by cAMP-dependent protein kinase) decreases the partitioning of myr-Src(2-19) onto vesicles containing acidic lipid 10-fold. We investigated the effect of phosphorylation on the localization of Src to biological membranes by expressing fusion constructs of Src's N terminus with a soluble carrier protein in COS-1 cells; phosphorylation produces a small shift in the distribution of the Src chimeras from the plasma membrane to the cytosol.

Amino Acid Sequence↗
PubMed1998

Electrostatic contributions to the binding free energy of the lambdacI repressor to DNA.

A model based on the nonlinear Poisson-Boltzmann (NLPB) equation is used to study the electrostatic contribution to the binding free energy of the lambdacI repressor to its operator DNA. In particular, we use the Poisson-Boltzmann model to calculate the pKa shift of individual ionizable amino acids upon binding. We find that three residues on each monomer, Glu34, Glu83, and the amino terminus, have significant changes in their pKa and titrate between pH 4 and 9. This information is then used to calculate the pH dependence of the binding free energy. We find that the calculated pH dependence of binding accurately reproduces the available experimental data over a range of physiological pH values. The NLPB equation is then used to develop an overall picture of the electrostatics of the lambdacI repressor-operator interaction. We find that long-range Coulombic forces associated with the highly charged nucleic acid provide a strong driving force for the interaction of the protein with the DNA. These favorable electrostatic interactions are opposed, however, by unfavorable changes in the solvation of both the protein and the DNA upon binding. Specifically, the formation of a protein-DNA complex removes both charged and polar groups at the binding interface from solvent while it displaces salt from around the nucleic acid. As a result, the electrostatic contribution to the lambdacI repressor-operator interaction opposes binding by approximately 73 kcal/mol at physiological salt concentrations and neutral pH. A variety of entropic terms also oppose binding. The major force driving the binding process appears to be release of interfacial water from the protein and DNA surfaces upon complexation and, possibly, enhanced packing interactions between the protein and DNA in the interface. When the various nonelectrostatic terms are described with simple models that have been applied previously to other binding processes, a general picture of protein/DNA association emerges in which binding is driven by the nonpolar interactions, whereas specificity results from electrostatic interactions that weaken binding but are necessary components of any protein/DNA complex.

DNA, Viral↗
PubMed1997

Electrostatic interaction of myristoylated proteins with membranes: simple physics, complicated biology.

Cell membrane association by several important peripheral proteins, such as Src, MARCKS, HIV-1 Gag, and K-Ras, requires nonspecific electrostatic interactions between a cluster of basic residues on the protein and acidic phospholipids in the plasma membrane. A simple theoretical model based on the nonlinear Poisson-Boltzmann equation describes well the experimentally measured electrostatic association between such proteins and the cell membrane.

Amino Acid Sequence↗
PubMed1997

On the calculation of binding free energies using continuum methods: application to MHC class I protein-peptide interactions.

This paper describes a methodology to calculate the binding free energy (delta G) of a protein-ligand complex using a continuum model of the solvent. A formal thermodynamic cycle is used to decompose the binding free energy into electrostatic and non-electrostatic contributions. In this cycle, the reactants are discharged in water, associated as purely nonpolar entities, and the final complex is then recharged. The total electrostatic free energies of the protein, the ligand, and the complex in water are calculated with the finite difference Poisson-Boltzmann (FDPB) method. The nonpolar (hydrophobic) binding free energy is calculated using a free energy-surface area relationship, with a single alkane/water surface tension coefficient (gamma aw). The loss in backbone and side-chain configurational entropy upon binding is estimated and added to the electrostatic and the nonpolar components of delta G. The methodology is applied to the binding of the murine MHC class I protein H-2Kb with three distinct peptides, and to the human MHC class I protein HLA-A2 in complex with five different peptides. Despite significant differences in the amino acid sequences of the different peptides, the experimental binding free energy differences (delta delta Gexp) are quite small (< 0.3 and < 2.7 kcal/mol for the H-2Kb and HLA-A2 complexes, respectively). For each protein, the calculations are successful in reproducing a fairly small range of values for delta delta Gcalc (< 4.4 and < 5.2 kcal/mol, respectively) although the relative peptide binding affinities of H-2Kb and HLA-A2 are not reproduced. For all protein-peptide complexes that were treated, it was found that electrostatic interactions oppose binding whereas nonpolar interactions drive complex formation. The two types of interactions appear to be correlated in that larger nonpolar contributions to binding are generally opposed by increased electrostatic contributions favoring dissociation. The factors that drive the binding of peptides to MHC proteins are discussed in light of our results.

Animals↗
PubMed1997

Electrostatic binding of proteins to membranes. Theoretical predictions and experimental results with charybdotoxin and phospholipid vesicles.

We previously applied the Poisson-Boltzmann equation to atomic models of phospholipid bilayers and basic peptides to calculate their electrostatic interactions from first principles (Ben-Tal, N., B. Honig, R. M. Peitzsch, G. Denisov, and S. McLaughlan. 1996. Binding of small basic peptides to membranes containing acidic lipids. Theoretical models and experimental results. Biophys. J. 71:561-575). Specifically, we calculated the molar partition coefficient, K (the reciprocal of the lipid concentration at which 1/2 the peptide is bound), of simple basic peptides (e.g., pentalysine) with phospholipid vesicles. The theoretical predictions agreed well with experimental measurements of the binding, but the agreement could have been fortuitous because the structure(s) of these flexible peptides is not known. Here we use the same theoretical approach to calculate the membrane binding of two small proteins of known structure: charybdotoxin (CTx) and iberiotoxin (IbTx); we also measure the binding of these proteins to phospholipid vesicles. The theoretical model describes accurately the dependence of K on the ionic strength and mol % acidic lipid in the membrane for both CTx (net charge +4) and IbTx (net charge +2). For example, the theory correctly predicts that the value of K for the binding of CTx to a membrane containing 33% acidic lipid should decrease by a factor of 10(5) when the salt concentration increases from 10 to 200 mM. We discuss the limitations of the theoretical approach and also consider a simple extension of the theory that incorporates nonpolar interactions.

Biophysical Phenomena↗
PubMed1996

Free energy determinants of secondary structure formation: III. beta-turns and their role in protein folding.

The stability of beta-turns is calculated as a function of sequence and turn type with a Monte Carlo sampling technique. The conformational energy of four internal hydrogen-bonded turn types, I, I', II and II', is obtained by evaluating their gas phase energy with the CHARMM force field and accounting for solvation effects with the Finite Difference Poisson-Boltzmann (FDPB) method. All four turn types are found to be less stable than the coil state, independent of the sequence in the turn. The free-energy penalties associated with turn formation vary between 1.6 kcal/mol and 7.7 kcal/mol, depending on the sequence and turn type. Differences in turn stability arise mainly from intraresidue interactions within the two central residues of the turn. For each combination of the two central residues, except for -Gly-Gly-, the most stable beta-turn type is always found to occur most commonly in native proteins. The fact that a model based on local interactions accounts for the observed preference of specific sequences suggests that long-range tertiary interactions tend to play a secondary role in determining turn conformation. In contrast, for beta-hairpins, long-range interactions appear to dominate. Specifically, due to the right-handed twist of beta-strands, type I' turns for -Gly-Gly- are found to occur with high frequency, even when local energetics would dictate otherwise. The fact that any combination of two residues is found able to adopt a relatively low-energy turn structure explains why the amino acid sequence in turns is highly variable. The calculated free-energy cost of turn formation, when combined with related numbers obtained for alpha-helices and beta-sheets, suggests a model for the initiation of protein folding based on metastable fragments of secondary structure.

Data Interpretation, Statistical↗
PubMed1996

The electrostatic contribution to the B to Z transition of DNA.

In this paper, the finite difference nonlinear Poisson-Boltzmann (NLPB) equation is used to calculate the electrostatic contribution to the B to Z transition of DNA using detailed molecular structures of each DNA form. The electrostatic transition free energy is described as a balance between the change in intramolecular Coulombic interactions and charge-dependent interactions between the DNA and the solvent. As in many prior studies, we find that the larger electrostatic repulsions among the more closely spaced Z-DNA phosphates destabilize this form compared to B-DNA in the absence of solvent. However, as a result of the more compact three-dimensional geometry of Z-DNA, both water and salt are found to strongly stabilize this conformation to the extent that the total electrostatic free energy favors the B to Z transition in aqueous solution. Water acts not only by screening the inter-phosphate repulsions but also by solvating both charged and polar groups on Z-DNA more favorably than B-DNA. In addition, Z-DNA is stabilized by a substantially higher concentration of nearby counterions than B-DNA. The relative stabilization of Z-DNA by salt increases with increasing bulk salt concentration, leading to the high-salt B to Z transition. We find that the salt dependence of the B to Z transition free energy calculated with the NLPB equation agrees reasonably well with experimental results. Since electrostatic interactions are found to favor the Z-form, nonelectrostatic forces must be responsible for the relative stability of B-DNA in solution. An analysis of these forces suggests that the conformational entropy may play an important role.

Computer Simulation↗
PubMed1996

Statistical thermodynamic analysis of peptide and protein insertion into lipid membranes.

A statistical thermodynamic approach is used to analyze the various contributions to the free energy change associated with the insertion of proteins and protein fragments into lipid bilayers. The partition coefficient that determines the equilibrium distribution of proteins between the membrane and the solution is expressed as the ratio between the partition functions of the protein in the two phases. It is shown that when all of the relevant degrees of freedom (i.e., those that change their character upon insertion into the membrane) can be treated classically, the partition coefficient is fully determined by the ratio of the configurational integrals and thus does not involve any mass-dependent factors, a conclusion that is also valid for related processes such as protein adsorption on a membrane surface or substrate binding to proteins. The partition coefficient, and hence the transfer free energy, depend only on the potential energy of the protein in the membrane. Expressing this potential as a sum of a "static" term, corresponding to the equilibrium (minimal free energy) configuration of the protein in the membrane, and a "dynamical" term representing fluctuations around the equilibrium configuration, we show that the static term contains the "solvation" and "lipid perturbation" contributions to the transfer free energy. The dynamical term is responsible for the "immobilization" free energy, reflecting the loss of translational and rotational entropy of the protein upon incorporation into the membrane. Based on a recent molecular theory of lipid-protein interactions, the lipid perturbation and immobilization contributions are then expressed in terms of the elastic deformation free energy resulting from the perturbation of the lipid environment by the foreign (protein) inclusion. The model is formulated for cylindrically shaped proteins, and numerical estimates are given for the insertion of an alpha-helical peptide into a lipid bilayer. The immobilization free energy is shown to be considerably smaller than in previous estimates of this quantity, and the origin of the difference is discussed in detail.

Biophysical Phenomena↗
PubMed1996

Binding of small basic peptides to membranes containing acidic lipids: theoretical models and experimental results.

We measured directly the binding of Lys3, Lys5, and Lys7 to vesicles containing acidic phospholipids. When the vesicles contain 33% acidic lipids and the aqueous solution contains 100 mM monovalent salt, the standard Gibbs free energy for the binding of these peptides is 3, 5, and 7 kcal/mol, respectively. The binding energies decrease as the mol% of acidic lipids in the membrane decreases and/or as the salt concentration increases. Several lines of evidence suggest that these hydrophilic peptides do not penetrate the polar headgroup region of the membrane and that the binding is mainly due to electrostatic interactions. To calculate the binding energies from classical electrostatics, we applied the nonlinear Poisson-Boltzmann equation to atomic models of the phospholipid bilayers and the basic peptides in aqueous solution. The electrostatic free energy of interaction, which arises from both a long-range coulombic attraction between the positively charged peptide and the negatively charged lipid bilayer, and a short-range Born or image charge repulsion, is a minimum when approximately 2.5 A (i.e., one layer of water) exists between the van der Waals surfaces of the peptide and the lipid bilayer. The calculated molar association constants, K, agree well with the measured values: K is typically about 10-fold smaller than the experimental value (i.e., a difference of about 1.5 kcal/mol in the free energy of binding). The predicted dependence of K (or the binding free energies) on the ionic strength of the solution, the mol% of acidic lipids in the membrane, and the number of basic residues in the peptide agree very well with the experimental measurements. These calculations are relevant to the membrane binding of a number of important proteins that contain clusters of basic residues.

Amino Acid Sequence↗
PubMed1996

Electrostatic coupling between retinal isomerization and the ionization state of Glu-204: a general mechanism for proton release in bacteriorhodopsin.

The pKa values of ionizable groups that lie between the active site region of bacteriorhodopsin (bR) and the extracellular surface of the protein are reported. Glu-204 is found to have an elevated pKa in the resting state of bR, suggesting that it corresponds to the proton-releasing group in bR. Its elevated pKa is predicted to be due in part to strong repulsive interactions with Glu-9. Following trans-cis isomerization of the retinal chromophore and the transfer of a proton to Asp-85, polar groups on the protein are able to interact more strongly with the ionized state of Glu-204, leading to a substantial reduction of its pKa. This suggests a general mechanism for proton release in which isomerization and subsequent charge separation initially produce a new electrostatic balance in the active site of bR. Here it is proposed that those events in turn drives a conformational change in the protein in which the ionized state of Glu-204 can be stabilized through interactions with groups that were previously inaccessible. Whether these groups should be identified with polar moieties in the protein, bound waters, or Arg-82 is an important mechanistic question whose elucidation will require further study.

Bacteriorhodopsins↗
PubMed1996

Helix-helix interactions in lipid bilayers.

Using a continuum model, we calculated the electrostatic interaction free energy between two alpha-helices in three environments: the aqueous phase, a low dielectric alkane phase, and a simple representation of a lipid bilayer. As was found in previous work, helix-helix interactions in the aqueous phase are quite weak, because of solvent screening, and slightly repulsive, because of desolvation effects that accompany helix assembly. In contrast, the interactions can be quite strong in a hypothetical alkane phase because desolvation effects are essentially nonexistent and because helix-helix interactions are not well screened. In this type of environment, the antiparallel helix orientation is strongly favored over the parallel orientation. In previous work we found that the free energy penalty associated with burying helix termini in a bilayer is quite high, which is why the termini tend to protrude into the solvent. Under these conditions the electrostatic interaction is strongly screened by solvent; indeed, it is sufficient for the termini to protrude a few angstroms from the two surfaces of the bilayer for their interaction to diminish almost completely. The effect is consistent with the classical model of the helix dipole in which the dipole moment is represented by point charges located at either terminus. Our results suggest, in agreement with previous models, that there is no significant nonspecific driving force for helix aggregation and, hence, that membrane protein folding must be driven by specific interactions such as close packing and salt-bridge and hydrogen bond formation.

Computer Simulation↗
PubMed1996

Free-energy determinants of alpha-helix insertion into lipid bilayers.

A detailed treatment is provided of the various free-energy terms that contribute to the transfer of a polyalanine alpha-helix from the aqueous phase into lipid bilayers. In agreement with previous work, the hydrophobic effect is found to provide the major driving force for helix insertion. However, an opposing effect of comparable magnitude is also identified and is attributed to the large free-energy penalty associated with the desolvation of peptide hydrogen bonds on transfer to the low dielectric environment of the bilayer. Lipid perturbation effects as well as the entropy loss associated with helix immobilization in the bilayer are also evaluated. Two configurations of a membrane-bound 25mer polyalanine helix were found to be lower in free energy than the isolated helix in the aqueous phase. The first corresponds to the case of vertical insertion, in which a helix terminus protrudes from each side of the bilayer. The second minimum is for the case of horizontal insertion, for which the helix is adsorbed upon the surface of the bilayer. The calculated free-energy minima are found to be in good agreement with recent measurements of related systems. Large free-energy barriers resulting from desolvation of unsatisfied hydrogen-bonding groups at the helix termini are obtained for both insertion processes. The barriers for insertion are significantly reduced if the helix termini are assumed to be "capped" through the formation of hydrogen bonds with polar sidechains. For uncapped helices, our results support recently proposed models in which helices are inserted by first adsorbing on the membrane surface and then having one terminus "swing around" so as to penetrate the bilayer.

Biophysical Phenomena↗
PubMed1996

Calculated coupling of electron and proton transfer in the photosynthetic reaction center of Rhodopseudomonas viridis.

Based on new Rhodopseudomonas (Rp.) viridis reaction center (RC) coordinates with a reliable structure of the secondary acceptor quinone (QB) site, a continuum dielectric model and finite difference technique have been used to identify clusters of electrostatically interacting ionizable residues. Twenty-three residues within a distance of 25 A from QB (QB cluster) have been shown to be strongly electrostatically coupled to QB, either directly or indirectly. An analogous cluster of 24 residues is found to interact with QA (QA cluster). Both clusters extend to the cytoplasmic surface in at least two directions. However, the QB cluster differs from the QA cluster in that it has a surplus of acidic residues, more strong electrostatic interactions, is less solvated, and experiences a strong positive electrostatic field arising from the polypeptide backbone. Consequently, upon reduction of QA or QB, it is the QB cluster, and not the QA cluster, which is responsible for substoichiometric proton uptake at neutral pH. The bulk of the changes in the QB cluster are calculated to be due to the protonation of a tightly coupled cluster of the three Glu residues (L212, H177, and M234) within the QB cluster. If the lifetime of the doubly reduced state QB2- is long enough, Asp M43 and Ser L223 are predicted to also become protonated. The calculated complex titration behavior of the strongly interacting residues of the QB cluster and the resulting electrostatic response to electron transfer may be a common feature in proton-transferring membrane protein complexes.

Amino Acid Sequence↗
PubMed1996

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.

Hydrogen Bonding↗
PubMed1995

Free energy determinants of secondary structure formation: I. alpha-Helices.

The Zimm-Bragg parameters s and sigma are calculated for the helix-coil transition of poly-L-alanine. The theoretical approach involves evaluating gas phase conformational energies for both coil and helical states using the CHARMM potential function and accounting for solvation effects with various continuum solvation models. Conformational free energies are then incorporated into a formalism developed by Go et al. for the calculation of s and sigma. Calculated values for both s and sigma as well as the enthalpy change associated with helix formation are in good agreement with experimental data when the Finite Difference Poisson-Boltzmann (FDPB) method is used to treat solvent effects. The driving force for the helix-coil transition is analyzed in terms of individual free energy components. Hydrogen bond formation is found to contribute little to helix stability because the internal hydrogen bonding energy is largely canceled by the large free energy cost associated with removing polar groups from water. The entropic cost associated with fixing backbone dihedral angles in the helical conformation is found to be approximately 7 e.u./residue (about 2 kcal/mol at room temperature). The major driving force favoring helix formation can be associated with interactions including enhanced van der Waals interactions in the close-packed helix conformation and the hydrophobic effect. These contribute about 2 kcal/mol favoring the helical state. The differences in helical propensities between alanine and glycine are attributed primarily to hydrophobic and packing interactions involving the C beta with a smaller contribution arising from increased conformational freedom for glycine in the coil state. The description of helix formation presented here is consistent with previous conclusions regarding tertiary structure formation which suggest that hydrophobic and close-packed interactions provide stability while hydrogen bond formation constitutes a structural constraint imposed by the high free energy cost associated with burying unsatisfied hydrogen bonding groups. alpha-Helix formation may thus be viewed as a form of hydrophobic collapse constrained by the requirement that polar groups be either exposed to solvent or form hydrogen bonds. More generally it appears from this study that for a folding model to be a realistic, it must properly account for the chemical nature of the polypeptide chain, particularly the solvation energetics of amide groups.

Amino Acid Sequence↗