PubMed Health⌕ Search

Biomedical subjects

N Ben-Tal

Publications and source records attributed to N Ben-Tal.

16 recordsLinked to original sources
PubMed2000

Theoretical calculations of the permeability of monensin-cation complexes in model bio-membranes.

Monensin is one of the best-characterized ionophores; it functions in the electroneutral exchange of cations between the extracellular and cytoplasmic sides of cell membranes. The X-ray crystal structures of monensin in free acid form and in complex with Na(+), K(+) and Ag(+) are known and we have recently measured the diffusion rates of monensin in free acid form (Mo-H) and in complex with Na(+) (Mo-Na) and with K(+) (Mo-K) using laser pulse techniques. The results have shown that Mo-H diffuses across the membrane one order of magnitude faster than Mo-Na and two orders of magnitude faster than Mo-K. Here, we report calculations of the translocation free energy of these complexes across the membrane along the most favorable path, i.e. the lowest free energy path. The calculations show that the most favorable orientation of monensin is with its hydrophobic furanyl and pyranyl moieties in the hydrocarbon region of the membrane and the carboxyl group and the cation at the water-membrane interface. Further, the calculations show that Mo-H is likely to be inserted deeper than Mo-Na into the bilayer, and that the free energy barrier for transfer of Mo-H across the membrane is approximately 1 kcal/mol lower than for Mo-Na, in good agreement with our measurements. Our results show that the Mo-K complex is unlikely to diffuse across lipid bilayers in its X-ray crystal structure, in contrast to the Mo-H and Mo-Na complexes. Apparently, when diffusing across the membrane, the Mo-K complex assumes a different conformation and/or thinning defects in the bilayer lower significantly the free energy barrier for the process. The suitability of the model for treating the membrane association of small molecules is discussed in view of the successes and failures observed for the monensin system.

Cations, Monovalent↗
PubMed2000

Association entropy in adsorption processes.

The association of two species to form a bound complex, e.g., the binding of a ligand to a protein or the adsorption of a peptide on a lipid membrane, involves an entropy loss, reflecting the conversion of free translational and rotational degrees of freedom into bound motions. Previous theoretical estimates of the standard entropy change in bimolecular binding processes, DeltaS(o), have been derived from the root-mean-square fluctuations in protein crystals, suggesting DeltaS(o) approximately -50 e.u., i.e., TDeltaS degrees approximately -25 kT = -15 kcal/mol. In this work we focus on adsorption, rather than binding processes. We first present a simple statistical-thermodynamic scheme for calculating the adsorption entropy, including its resolution into translational and rotational contributions, using the known distance-orientation dependent binding (adsorption) potential. We then utilize this scheme to calculate the free energy of interaction and entropy of pentalysine adsorption onto a lipid membrane, obtaining TDeltaS(o) approximately -1.7 kT approximately -1.3 kcal/mol. Most of this entropy change is due to the conversion of one free translation into a bound motion, the rest arising from the confinement of two rotational degrees of freedom. The smaller entropy loss in adsorption compared to binding processes arises partly because a smaller number of degrees of freedom become restricted, but mainly due to the fact that the binding potential is much "softer."

Adsorption↗
PubMed2000

Calculations suggest a pathway for the transverse diffusion of a hydrophobic peptide across a lipid bilayer.

Alamethicin is a hydrophobic antibiotic peptide 20 amino acids in length. It is predominantly helical and partitions into lipid bilayers mostly in transmembrane orientations. The rate of the peptide transverse diffusion (flip-flop) in palmitoyl-oleyl-phosphatidylcholine vesicles has been measured recently and the results suggest that it involves an energy barrier, presumably due to the free energy of transfer of the peptide termini across the bilayer. We used continuum-solvent model calculations, the known x-ray crystal structure of alamethicin and a simplified representation of the lipid bilayer as a slab of low dielectric constant to calculate the flip-flop rate. We assumed that the lipids adjust rapidly to each configuration of alamethicin in the bilayer because their motions are significantly faster than the average peptide flip-flop time. Thus, we considered the process as a sequence of discrete peptide-membrane configurations, representing critical steps in the diffusion, and estimated the transmembrane flip-flop rate from the calculated free energy of the system in each configuration. Our calculations indicate that the simplest possible pathway, i.e., the rotation of the helix around the bilayer midplane, involving the simultaneous burial of the two termini in the membrane, is energetically unfavorable. The most plausible alternative is a two-step process, comprised of a rotation of alamethicin around its C-terminus residue from the initial transmembrane orientation to a surface orientation, followed by a rotation around the N-terminus residue from the surface to the final reversed transmembrane orientation. This process involves the burial of one terminus at a time and is much more likely than the rotation of the helix around the bilayer midplane. Our calculations give flip-flop rates of approximately 10(-7)/s for this pathway, in accord with the measured value of 1.7 x 10(-6)/s.

Alamethicin↗
PubMed2000

Continuum solvent model calculations of alamethicin-membrane interactions: thermodynamic aspects.

Alamethicin is a 20-amino acid antibiotic peptide that forms voltage-gated ion channels in lipid bilayers. Here we report calculations of its association free energy with membranes. The calculations take into account the various free-energy terms that contribute to the transfer of the peptide from the aqueous phase into bilayers of different widths. The electrostatic and nonpolar contributions to the solvation free energy are calculated using continuum solvent models. The contributions from the lipid perturbation and membrane deformation effects and the entropy loss associated with peptide immobilization in the bilayer are estimated from a statistical thermodynamic model. The calculations were carried out using two classes of experimentally observed conformations, both of which are helical: the NMR and the x-ray crystal structures. Our calculations show that alamethicin is unlikely to partition into bilayers in any of the NMR conformations because they have uncompensated backbone hydrogen bonds and their association with the membrane involves a large electrostatic solvation free energy penalty. In contrast, the x-ray conformations provide enough backbone hydrogen bonds for the peptide to associate with bilayers. We tested numerous transmembrane and surface orientations of the peptide in bilayers, and our calculations indicate that the most favorable orientation is transmembrane, where the peptide protrudes approximately 4 A into the water-membrane interface, in very good agreement with electron paramagnetic resonance and oriented circular dichroism measurements. The calculations were carried out using two alamethicin isoforms: one with glutamine and the other with glutamate in the 18th position. The calculations indicate that the two isoforms have similar membrane orientations and that their insertion into the membrane is likely to involve a 2-A deformation of the bilayer, again, in good agreement with experimental data. The implications of the results for the biological function of alamethicin and its capacity to oligomerize and form ion channels are discussed.

Alamethicin↗
PubMed1999

kPROT: a knowledge-based scale for the propensity of residue orientation in transmembrane segments. Application to membrane protein structure prediction.

Modeling of integral membrane proteins and the prediction of their functional sites requires the identification of transmembrane (TM) segments and the determination of their angular orientations. Hydrophobicity scales predict accurately the location of TM helices, but are less accurate in computing angular disposition. Estimating lipid-exposure propensities of the residues from statistics of solved membrane protein structures has the disadvantage of relying on relatively few proteins. As an alternative, we propose here a scale of knowledge-based Propensities for Residue Orientation in Transmembrane segments (kPROT), derived from the analysis of more than 5000 non-redundant protein sequences. We assume that residues that tend to be exposed to the membrane are more frequent in TM segments of single-span proteins, while residues that prefer to be buried in the transmembrane bundle interior are present mainly in multi-span TMs. The kPROT value for each residue is thus defined as the logarithm of the ratio of its proportions in single and multiple TM spans. The scale is refined further by defining it for three discrete sections of the TM segment; namely, extracellular, central, and intracellular. The capacity of the kPROT scale to predict angular helical orientation was compared to that of alternative methods in a benchmark test, using a diversity of multi-span alpha-helical transmembrane proteins with a solved 3D structure. kPROT yielded an average angular error of 41 degrees, significantly lower than that of alternative scales (62 degrees -68 degrees ). The new scale thus provides a useful general tool for modeling and prediction of functional residues in membrane proteins. A WWW server (http://bioinfo.weizmann.ac.il/kPROT) is available for automatic helix orientation prediction with kPROT.

Amino Acid Sequence↗
PubMed1999

Electrostatic properties of membranes containing acidic lipids and adsorbed basic peptides: theory and experiment.

The interaction of heptalysine with vesicles formed from mixtures of the acidic lipid phosphatidylserine (PS) and the zwitterionic lipid phosphatidylcholine (PC) was examined experimentally and theoretically. Three types of experiments showed that smeared charge theories (e.g., Gouy-Chapman-Stern) underestimate the membrane association when the peptide concentration is high. First, the zeta potential of PC/PS vesicles in 100 mM KCl solution increased more rapidly with heptalysine concentration (14.5 mV per decade) than predicted by a smeared charge theory (6.0 mV per decade). Second, changing the net surface charge density of vesicles by the same amount in two distinct ways produced dramatically different effects: the molar partition coefficient decreased 1000-fold when the mole percentage of PS was decreased from 17% to 4%, but decreased only 10-fold when the peptide concentration was increased to 1 microM. Third, high concentrations of basic peptides reversed the charge on PS and PC/PS vesicles. Calculations based on finite difference solutions to the Poisson-Boltzmann equation applied to atomic models of heptalysine and PC/PS membranes provide a molecular explanation for the observations: a peptide adsorbing to the membrane in the presence of other surface-adsorbed peptides senses a local potential more negative than the average potential. The biological implications of these "discreteness-of-charge" effects are discussed.

2-Naphthylamine↗
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↗
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

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

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

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↗