Biomedical subjects

Nir Ben-Tal

Publications and source records attributed to Nir Ben-Tal.

36 records · Page 2Linked to original sources

PubMed2004

ConSeq: the identification of functionally and structurally important residues in protein sequences.

MOTIVATION: ConSeq is a web server for the identification of biologically important residues in protein sequences. Functionally important residues that take part, e.g. in ligand binding and protein-protein interactions, are often evolutionarily conserved and are most likely to be solvent-accessible, whereas conserved residues within the protein core most probably have an important structural role in maintaining the protein's fold. Thus, estimated evolutionary rates, as well as relative solvent accessibility predictions, are assigned to each amino acid in the sequence; both are subsequently used to indicate residues that have potential structural or functional importance. AVAILABILITY: The ConSeq web server is available at http://conseq.bioinfo.tau.ac.il/ SUPPLEMENTARY INFORMATION: The ConSeq methodology, a description of its performance in a set of five well-documented proteins, a comparison to other methods, and the outcome of its application to a set of 111 proteins of unknown function, are presented at http://conseq.bioinfo.tau.ac.il/ under 'OVERVIEW', 'VALIDATION', 'COMPARISON' and 'PREDICTIONS', respectively.

Algorithms↗

PubMed2004

A putative mechanism for downregulation of the catalytic activity of the EGF receptor via direct contact between its kinase and C-terminal domains.

Tyrosine kinase receptors of the EGFR family play a significant role in vital cellular processes and in various cancers. EGFR members are unique among kinases, as the regulatory elements of their kinase domains are constitutively ready for catalysis. Nevertheless, the receptors are not constantly active. This apparent paradox has prompted us to seek mechanisms of regulation in EGFR's cytoplasmic domain that do not involve conformational changes of the kinase domain. Our computational analyses, based on the three-dimensional structure of EGFR's kinase domain suggest that direct contact between the kinase and a segment from the C-terminal regulatory domains inhibits enzymatic activity. EGFR activation would then involve temporal dissociation of this stable complex, for example, via ligand-induced contact formation between the extracellular domains, leading to the reorientation of the transmembrane and intracellular domains. The model provides an explanation at the molecular level for the effects of several cancer-causing EGFR mutations.

Amino Acid Sequence↗

PubMed2004

Increased concentration of polyvalent phospholipids in the adsorption domain of a charged protein.

We studied the adsorption of a charged protein onto an oppositely charged membrane, composed of mobile phospholipids of differing valence, using a statistical-thermodynamical approach. A two-block model was employed, one block corresponding to the protein-affected region on the membrane, referred to as the adsorption domain, and the other to the unaffected remainder of the membrane. We calculated the protein-induced lipid rearrangement in the adsorption domain as arising from the interplay between the electrostatic interactions in the system and the mixing entropy of the lipids. Equating the electrochemical potentials of the lipids in the two blocks yields an expression for the relations among the various lipid fractions in the adsorption domain, indicating a sensitive dependence of lipid fraction on valence. This expression is a result of the two-block picture but does not depend on further details of the protein-membrane interaction. We subsequently calculated the lipid fractions themselves using the Poisson-Boltzmann theory. We examined the dependence of lipid enrichment, i.e., the ratio between the lipid fractions inside and outside the adsorption domain, on various parameters such as ionic strength and lipid valence. Maximum enrichment was found for lipid valence in the range between -3 and -4 in physiological conditions. Our results are in qualitative agreement with recent experimental studies on the interactions between peptides having a domain of basic residues and membranes containing a small fraction of the polyvalent phosphatidylinositol 4,5-bisphosphate (PIP2). This study provides theoretical support for the suggestion that proteins adsorbed onto membranes through a cluster of basic residues may sequester PIP2 and other polyvalent lipids.

Lipid Bilayers↗

PubMed2004

The conserved carboxy terminus of the capsid domain of human immunodeficiency virus type 1 gag protein is important for virion assembly and release.

The retroviral Gag precursor plays an important role in the assembly of virion particles. The capsid (CA) protein of the Gag molecule makes a major contribution to this process. In the crystal structure of the free CA protein of the human immunodeficiency virus type 1 (HIV-1), 11 residues of the C terminus were found to be unstructured, and to date no information exists on the structure of these residues in the context of the Gag precursor molecule. We performed phylogenetic analysis and demonstrated a high degree of conservation of these 11 amino acids. Deletion of this cluster or introduction of various point mutations into these residues resulted in significant impairment of particle infectivity. In this cluster, two putative structural regions were identified, residues that form a hinge region (353-VGGP-356) and those that contribute to an alpha-helix (357-GHKARVL-363). Overall, mutations in these regions resulted in inhibition of virion production, but mutations in the hinge region demonstrated the most significant reduction. Although all the Gag mutants appeared to have normal Gag-Gag and Gag-RNA interactions, the hinge mutants were characterized by abnormal formation of cytoplasmic Gag complexes. Gag proteins with mutations in the hinge region demonstrated normal membrane association but aberrant rod-like membrane structures. More detailed analysis of these structures in one of the mutants demonstrated abnormal trapped Gag assemblies. These data suggest that the conserved CA C terminus is important for HIV-1 virion assembly and release and define a putative target for drug design geared to inhibit the HIV-1 assembly process.

Amino Acid Sequence↗

PubMed2004

Free diffusion of steroid hormones across biomembranes: a simplex search with implicit solvent model calculations.

Steroid hormones such as progesterone, testosterone, and estradiol are derived from cholesterol, a major constituent of biomembranes. Although the hormones might be expected to associate with the bilayer in a fashion similar to that of cholesterol, their biological action in regulating transcription of target genes involves transbilayer transfer by free diffusion, which is not observed for cholesterol. We used a novel combination of a continuum-solvent model and the downhill simplex search method for the calculation of the free energy of interaction of these hormones with lipid membranes, and compared these values to that of cholesterol-membrane interaction. The hormones were represented in atomic detail and the membrane as a structureless hydrophobic slab embedded in implicit water. A deep free-energy minimum of approximately -15 kcal/mol was obtained for cholesterol at its most favorable location in the membrane, whereas the most favorable locations for the hormones were associated with shallower minima of -5.0 kcal/mol or higher. The free-energy difference, which is predominantly due to the substitution of cholesterol's hydrophobic tail with polar groups, explains the different manner in which cholesterol and the hormones interact with the membrane. Further calculations were conducted to estimate the rate of transfer of the hormones from the aqueous phase into hexane, and from hexane back into the aqueous phase. The calculated rates agreed reasonably well with measurements in closely related systems. Based on these calculations, we suggest putative pathways for the free diffusion of the hormones across biomembranes. Overall, the calculations imply that the hormones may rapidly cross biomembrane barriers. Implications for gastrointestinal absorption and transfer across the blood-brain barrier and for therapeutic uses are discussed.

Algorithms↗

PubMed2003

Evolutionary analysis reveals collective properties and specificity in the C-type lectin and lectin-like domain superfamily.

Members of the C-type lectin/C-type lectin-like domain (CTL/CTLD) superfamily share a common fold and are involved in a variety of functions, such as generalized defense mechanisms against foreign agents, discrimination between healthy and pathogen-infected cells, and endocytosis and blood coagulation. In this work we used ConSurf, a computer program recently developed in our lab, to perform an evolutionary analysis of this superfamily in order to further identify characteristics of all or part of its members. Given a set of homologous proteins in the form of multiple sequence alignment (MSA) and an inferred phylogenetic tree, ConSurf calculates the conservation score in every alignment position, taking into account the relationships between the sequences and the physicochemical similarity between the amino acids. The scores are then color-coded onto the three-dimensional structure of one of the homologous proteins. We provide here and at http://ashtoret.tau.ac.il/ approximately sharon a detailed analysis of the conservation pattern obtained for the entire superfamily and for two subgroups of proteins: (a) 21 CTLs and (b) 11 heterodimeric CTLD toxins. We show that, in general, proteins of the superfamily have one face that is constructed mostly of conserved residues and another that is not, and we suggest that the former face is involved in binding to other proteins or domains. In the CTLs examined we detected a region of highly conserved residues, corresponding to the known calcium- and carbohydrate-binding site of the family, which is not conserved throughout the entire superfamily, and in the CTLD toxins we found a patch of highly conserved residues, corresponding to the known dimerization region of these proteins. Our analysis also detected patches of conserved residues with yet unknown function(s).

Amino Acid Sequence↗

PubMed2003

Implicit solvent model estimates of the stability of model structures of the alamethicin channel.

Alamethicin is a hydrophobic helical peptide of 20 residues, which oligomerizes to form ion-conducting channels in membranes. The behavior of an intact alamethicin channel in POPC bilayers was recently studied, using 2 ns molecular dynamics (MD) simulations of a model hexameric channel. These simulations produced numerous conformations of the channel. In the present study, we used 11 of these channel conformations and carried out continuum-solvent model calculations, similar to those used for the monomers in our previous studies, to investigate the energetics of the channel inside the lipid bilayer. Our results suggest that, out of the 11 channel conformations produced by the MD simulations, only four are stable inside the lipid bilayer, with water-to-membrane free energies of transfer ranging from approximately -6 to approximately -10 kcal/mol. Analysis of the results suggests two causes for the apparent instability of the remainder of the structures inside the lipid bilayer, both resulting from the desolvation of channel polar groups (i.e. their transfer from the aqueous phase into the bilayer). The first is specific, uncompensated backbone hydrogen bonds, which exist in the region of the channel exposed to the hydrocarbon of the lipid bilayer. The second is exposure of intra-pore water molecules to the surrounding lipid. Thus, the association of these structures with the membrane involves a large electrostatic desolvation free-energy penalty. The apparent conflict between continuum-solvent and MD calculations, and its significance for the interpretation of membrane proteins simulations, are discussed.

Alamethicin↗

PubMed2003

Interactions of hydrophobic peptides with lipid bilayers: Monte Carlo simulations with M2delta.

We introduce here a novel Monte Carlo simulation method for studying the interactions of hydrophobic peptides with lipid membranes. Each of the peptide's amino acids is represented as two interaction sites: one corresponding to the backbone alpha-carbon and the other to the side chain, with the membrane represented as a hydrophobic profile. Peptide conformations and locations in the membrane and changes in the membrane width are sampled using the Metropolis criterion, taking into account the underlying energetics. Using this method we investigate the interactions between the hydrophobic peptide M2delta and a model membrane. The simulations show that starting from an extended conformation in the aqueous phase, the peptide first adsorbs onto the membrane surface, while acquiring an ordered helical structure. This is followed by formation of a helical-hairpin and insertion into the membrane. The observed path is in agreement with contemporary understanding of peptide insertion into biological membranes. Two stable orientations of membrane-associated M2delta were obtained: transmembrane (TM) and surface, and the value of the water-to-membrane transfer free energy of each of them is in agreement with calculations and measurements on similar cases. M2delta is most stable in the TM orientation, where it assumes a helical conformation with a tilt of 14 degrees between the helix principal axis and the membrane normal. The peptide conformation agrees well with the experimental data; average root-mean-square deviations of 2.1 A compared to nuclear magnetic resonance structures obtained in detergent micelles and supported lipid bilayers. The average orientation of the peptide in the membrane in the most stable configurations reported here, and in particular the value of the tilt angle, are in excellent agreement with the ones calculated using the continuum-solvent model and the ones observed in the nuclear magnetic resonance studies. This suggests that the method may be used to predict the three-dimensional structure of TM peptides.

Amino Acid Sequence↗

PubMed2003

Interactions of the M2delta segment of the acetylcholine receptor with lipid bilayers: a continuum-solvent model study.

M2delta, one of the transmembrane segments of the nicotinic acetylcholine receptor, is a 23-amino-acid peptide, frequently used as a model for peptide-membrane interactions. In this and the companion article we describe studies of M2delta-membrane interactions, using two different computational approaches. In the present work, we used continuum-solvent model calculations to investigate key thermodynamic aspects of its interactions with lipid bilayers. M2delta was represented in atomic detail and the bilayer was represented as a hydrophobic slab embedded in a structureless aqueous phase. Our calculations show that the transmembrane orientation is the most favorable orientation of the peptide in the bilayer, in good agreement with both experimental and computational data. Moreover, our calculations produced the free energy of association of M2delta with the lipid bilayer, which, to our knowledge, has not been reported to date. The calculations included 10 structures of M2delta, determined by nuclear magnetic resonance in dodecylphosphocholine micelles. All the structures were found to be stable inside the lipid bilayer, although their water-to-membrane transfer free energies differed by as much as 12 kT. Although most of the structures were roughly linear, a single structure had a kink in its central region. Interestingly, this structure was found to be the most stable inside the lipid bilayer, in agreement with molecular dynamics simulations of the peptide and with the recently determined structure of the intact receptor. Our analysis showed that the kink reduced the polarity of the peptide in its central region by allowing the electrostatic masking of the Gln13 side chain in that area. Our calculations also showed a tendency for the membrane to deform in response to peptide insertion, as has been previously found for the membrane-active peptides alamethicin and gramicidin. The results are compared to Monte Carlo simulations of the peptide-membrane system, as presented in the accompanying article.

Alamethicin↗

PubMed2003

ConSurf: identification of functional regions in proteins by surface-mapping of phylogenetic information.

UNLABELLED: We recently developed algorithmic tools for the identification of functionally important regions in proteins of known three dimensional structure by estimating the degree of conservation of the amino-acid sites among their close sequence homologues. Projecting the conservation grades onto the molecular surface of these proteins reveals patches of highly conserved (or occasionally highly variable) residues that are often of important biological function. We present a new web server, ConSurf, which automates these algorithmic tools. ConSurf may be used for high-throughput characterization of functional regions in proteins. AVAILABILITY: The ConSurf web server is available at:http://consurf.tau.ac.il. SUPPLEMENTARY INFORMATION: A set of examples is available at http://consurf.tau.ac.il under 'GALLERY'.

Amino Acids↗

PubMed2002

A putative molecular-activation switch in the transmembrane domain of erbB2.

Overexpression of the receptor tyrosine kinase (RTK) erbB2 (also designated neu or HER2) was implicated in causing a variety of human cancers, including mammary and ovarian carcinomas. Ligand-induced receptor dimerization is critical for stimulation of the intrinsic protein tyrosine kinase (PTK) of RTKs. It was therefore proposed that PTK activity is stimulated as a result of the reorientation of the cytoplasmic domains within receptor dimers, leading to transautophosphorylation and stimulation of enzymatic activity. Here, we propose a molecular mechanism for rotation-coupled activation of the erbB2 receptor. Using a computational exploration of conformation space of the transmembrane (TM) segments of an erbB2 homodimer, we found two stable conformations of the TM domain. We suggest that these conformations correspond to the active and inactive states of erbB2, and that the receptor molecules may switch from one conformation to the other without crossing exceedingly unfavorable states. This model provides an explanation for the biochemical and oncogenic properties of erbB2, such as the effects of erbB2 overexpression on kinase activity and cell transformation. Furthermore, the opposing effects of the neu* activating oncogenic point mutation and the Val-655-->Ile single-nucleotide polymorphism shown to be linked to reduced risk of breast cancer are explained in terms of shifts in the equilibrium between the active and inactive states of erbB2 in vivo.

Allosteric Regulation↗

PubMed2002

Type I chaperonins: not all are created equal.

Type I chaperonins play an essential role in the folding of newly translated and stress-denatured proteins in eubacteria, mitochondria and chloroplasts. Since their discovery, the bacterial chaperonins have provided an excellent model system for investigating the mechanism by which chaperonins mediate protein folding. Due to the high conservation of the primary sequence among Type I chaperonins, it is generally accepted that organellar chaperonins function similar to the bacterial ones. However, recent studies indicate that the chloroplast and mitochondrial chaperonins possess unique structural and functional properties that distinguish them from their bacterial homologs. This review focuses on the unique properties of organellar chaperonins.

Bacteria↗

PubMed2002

Long-range effects on the binding of the influenza HA to receptors are mediated by changes in the stability of a metastable HA conformation.

X-ray studies show that influenza hemagglutinin (HA) forms an elongated structure connecting the influenza virus at one end to cell-surface receptors at the other. At neutral pH, the 20 N-terminal residues of HA2-referred to as the fusion peptide-are buried in a hydrophobic pocket, about 100 A away from the receptor-binding site, and thus seem unlikely to affect HA binding to the receptor. To test this assumption, we mutated residues in the fusion peptide, heterologically expressed the mutated proteins in COS7 cells, and examined their ability to bind fluorescently labeled red blood cells (RBCs). Surprisingly, a significantly reduced binding was recorded for some of the mutants. Ample experimental data indicate that HA has at least two forms: the native structure at neutral pH (N) that is metastable and the fusogenic form (F), observed at low pH, which is stable. Thus, a simple interpretation of our data is that HA can bind to its receptors at the RBC surface in the N form much more effectively than in the F (or in any other stable) form and that the altered binding properties are due to destabilizing effects of the mutations on the N form. That is, some of the mutations involve reduction in the free energy barrier between the N and F forms. This, in turn, leads to reduction in the population of the N form, which is the only form capable of binding to the cell-surface receptors. To explore this possibility, we estimated the stability free energy difference between HA wild-type (wt) and mutants in the N form using an empirical surface tension coefficient. The calculated stability differences correlated well with the measured binding, supporting the above interpretation. Our results are examined taking into account the available experimental data on the affinity of different soluble and membrane-attached forms of HA to its receptors.

Animals↗

PubMed2002

A novel scoring function for predicting the conformations of tightly packed pairs of transmembrane alpha-helices.

Pairs of helices in transmembrane (TM) proteins are often tightly packed. We present a scoring function and a computational methodology for predicting the tertiary fold of a pair of alpha-helices such that its chances of being tightly packed are maximized. Since the number of TM protein structures solved to date is small, it seems unlikely that a reliable scoring function derived statistically from the known set of TM protein structures will be available in the near future. We therefore constructed a scoring function based on the qualitative insights gained in the past two decades from the solved structures of TM and soluble proteins. In brief, we reward the formation of contacts between small amino acid residues such as Gly, Cys, and Ser, that are known to promote dimerization of helices, and penalize the burial of large amino acid residues such as Arg and Trp. As a case study, we show that our method predicts the native structure of the TM homodimer glycophorin A (GpA) to be, in essence, at the global score optimum. In addition, by correlating our results with empirical point mutations on this homodimer, we demonstrate that our method can be a helpful adjunct to mutation analysis. We present a data set of canonical alpha-helices from the solved structures of TM proteins and provide a set of programs for analyzing it (http://ashtoret.tau.ac.il/~sarel). From this data set we derived 11 helix pairs, and conducted searches around their native states as a further test of our method. Approximately 73% of our predictions showed a reasonable fit (RMS deviation <2A) with the native structures compared to the success rate of 8% expected by chance. The search method we employ is less effective for helix pairs that are connected via short loops (<20 amino acid residues), indicating that short loops may play an important role in determining the conformation of alpha-helices in TM proteins.

Cell Membrane↗

PubMed2002

Stability of an ion channel in lipid bilayers: implicit solvent model calculations with gramicidin.

Gramicidin is a helical peptide, 15 residues in length, which dimerizes to form ion-conducting channels in lipid bilayers. Here we report calculations of its free energy of transfer from the aqueous phase into bilayers of different widths. The electrostatic and nonpolar contributions to the desolvation free energy were calculated using implicit solvent models, in which gramicidin was described in atomic detail and the hydrocarbon region of the membrane was described as a slab of hydrophobic medium embedded in water. The free energy penalties from the lipid perturbation and membrane deformation effects, and the entropy loss associated with gramicidin immobilization in the bilayer, were estimated from a statistical thermodynamic model of the bilayer. The calculations were carried out using two classes of experimentally observed conformations: a head-to-head dimer of two single-stranded (SS) beta-helices and a double-stranded (DS) intertwined double helix. The calculations showed that gramicidin is likely to partition into the bilayer in all of these conformations. However, the SS conformation was found to be significantly more stable than the DS in the bilayer, in agreement with most of the experimental data. We tested numerous transmembrane and surface orientations of gramicidin in bilayers of various widths. Our calculations indicate that the most favorable orientation is transmembrane, which is indeed to be expected from a channel-forming peptide. The calculations demonstrate that gramicidin insertion into the membrane is likely to involve a significant deformation of the bilayer to match the hydrophobic width of the peptide (22 A), again in good agreement with experimental data. Interestingly, deformation of the bilayer was induced by all of the gramicidin conformations.

Gramicidin↗

PubMed2002

On the regulatory role of dipeptidyl peptidase IV (=CD=adenosine deaminase complexing protein) on adenosine deaminase activity.

The molecular mechanism controlling the variable activity of the malignancy marker adenosine deaminase (ADA) is enigmatic. ADA activity was found to be modulated by the membrane-bound adenosine deaminase complexing protein (CP=DPPIV=CD26). The role of lipid-protein interactions in this modulation was sought. While direct solubilization of ADA in vesicles resulted in loss of ADA activity, the binding of ADA to CP reconstituted in vesicles restored the specific activity. The activity of ADA, free or bound to CP in solution, resulted in continuous linear Arrhenius plots. However, ADA bound to reconstituted CP exhibited two breaks associated with approximately 30% increased activity, at 25 and 13 degrees C, yielding three lines with similar apparent activation energies (E(a)). Continuum solvent model calculations of the free energy of transfer of the transmembrane helix of CP from the aqueous phase into membranes of various widths show that the most favorable orientations of the helix above and below the main phase transition may be different. We suggest that the 20% change in the thickness of the bilayer below and above the main phase transition may modify the orientation of CP in the membrane, thereby affecting substrate accessibility of ADA. This could account for ADA's reduced activity associated with increased membrane fluidity in transformed vs. normal fibroblasts.

Adenosine Deaminase↗

PubMed2002

A mutation in GJB3 is associated with recessive erythrokeratodermia variabilis (EKV) and leads to defective trafficking of the connexin 31 protein.

Erythrokeratodermia variabilis (EKV) is a skin disorder characterized by variable (transient) erythemas and fixed keratosis. The disorder maps to chromosome 1p34-35, a location that contains the GJB3 gene encoding the gap junction protein connexin 31. Until now, only heterozygote mutations in the form of dominant inheritance have been described in this gene associated with EKV. We report here a homozygote mutation in the connexin 31 gene, found in a family that shows recessive inheritance of the disorder, thus providing the first molecular support for a recessive variant of EKV. The entire GJB3 coding sequence was scanned for mutations by sequencing. We detected a T-->C transition at position 101 of the coding sequence, which replaces a leucine with a proline at residue 34 of the protein (L34P). Evolutionary analysis shows that this mutation is located at a highly conserved region of connexin in the first putative transmembrane helix (TMH). In transfected keratinocytes, L34P connexin 31 had a cytoplasmic distribution, suggesting that the mutant form of this protein will not form normal gap junctions between adjacent cells. The change of leucine to proline is likely to alter the structure of the first TMH of connexin by inducing a kink, thus influencing connexon structure and function.

Amino Acid Sequence↗

PubMed2002

Rate4Site: an algorithmic tool for the identification of functional regions in proteins by surface mapping of evolutionary determinants within their homologues.

MOTIVATION: A number of proteins of known three-dimensional (3D) structure exist, with yet unknown function. In light of the recent progress in structure determination methodology, this number is likely to increase rapidly. A novel method is presented here: 'Rate4Site', which maps the rate of evolution among homologous proteins onto the molecular surface of one of the homologues whose 3D-structure is known. Functionally important regions often correspond to surface patches of slowly evolving residues. RESULTS: Rate4Site estimates the rate of evolution of amino acid sites using the maximum likelihood (ML) principle. The ML estimate of the rates considers the topology and branch lengths of the phylogenetic tree, as well as the underlying stochastic process. To demonstrate its potency, we study the Src SH2 domain. Like previously established methods, Rate4Site detected the SH2 peptide-binding groove. Interestingly, it also detected inter-domain interactions between the SH2 domain and the rest of the Src protein that other methods failed to detect.

Algorithms↗