Search PubMed⌕ Search

Biomedical subjects

Sarel J Fleishman

Publications and source records attributed to Sarel J Fleishman.

16 recordsLinked to original sources
PubMed2006

Quasi-symmetry in the cryo-EM structure of EmrE provides the key to modeling its transmembrane domain.

Small multidrug resistance (SMR) transporters contribute to bacterial resistance by coupling the efflux of a wide range of toxic aromatic cations, some of which are commonly used as antibiotics and antiseptics, to proton influx. EmrE is a prototypical small multidrug resistance transporter comprising four transmembrane segments (M1-M4) that forms dimers. It was suggested recently that EmrE molecules in the dimer have different topologies, i.e. monomers have opposite orientations with respect to the membrane plane. A 3-D structure of EmrE acquired by electron cryo-microscopy (cryo-EM) at 7.5 Angstroms resolution in the membrane plane showed that parts of the structure are related by quasi-symmetry. We used this symmetry relationship, combined with sequence conservation data, to assign the transmembrane segments in EmrE to the densities seen in the cryo-EM structure. A C alpha model of the transmembrane region was constructed by considering the evolutionary conservation pattern of each helix. The model is validated by much of the biochemical data on EmrE with most of the positions that were identified as affecting substrate translocation being located around the substrate-binding cavity. A suggested mechanism for proton-coupled substrate translocation in small multidrug resistance antiporters provides a mechanistic rationale to the experimentally observed inverted topology.

Amino Acid Sequence↗
PubMed2006

The structural context of disease-causing mutations in gap junctions.

Gap junctions form intercellular channels that mediate metabolic and electrical signaling between neighboring cells in a tissue. Lack of an atomic resolution structure of the gap junction has made it difficult to identify interactions that stabilize its transmembrane domain. Using a recently computed model of this domain, which specifies the locations of each amino acid, we postulated the existence of several interactions and tested them experimentally. We introduced mutations within the transmembrane domain of the gap junction-forming protein connexin that were previously implicated in genetic diseases and that apparently destabilized the gap junction, as evidenced here by the absence of the protein from the sites of cell-cell apposition. The model structure helped identify positions on adjacent helices where second-site mutations restored membrane localization, revealing possible interactions between residue pairs. We thus identified two putative salt bridges and one pair involved in packing interactions in which one disease-causing mutation suppressed the effects of another. These results seem to reveal some of the physical forces that underlie the structural stability of the gap junction transmembrane domain and suggest that abrogation of such interactions bring about some of the effects of disease-causing mutations.

Cell Communication↗
PubMed2006

Progress in structure prediction of alpha-helical membrane proteins.

Transmembrane (TM) proteins comprise 20-30% of the genome but, because of experimental difficulties, they represent less than 1% of the Protein Data Bank. The dearth of membrane protein structures makes computational prediction a potentially important means of obtaining novel structures. Recent advances in computational methods have been combined with experimental data to constrain the modeling of three-dimensional structures. Furthermore, threading and ab initio modeling approaches that were effective for soluble proteins have been applied to TM domains. Surprisingly, experimental structures, proteomic analyses and bioinformatics have revealed unexpected architectures that counter long-held views on TM protein structure and stability. Future computational and experimental studies aimed at understanding the thermodynamic and evolutionary bases of these architectural details will greatly enhance predictive capabilities.

Animals↗
PubMed2006

Intrinsically disordered C-terminal segments of voltage-activated potassium channels: a possible fishing rod-like mechanism for channel binding to scaffold proteins.

Membrane-embedded voltage-activated potassium channels (Kv) bind intracellular scaffold proteins, such as the Post Synaptic Density 95 (PSD-95) protein, using a conserved PDZ-binding motif located at the channels' C-terminal tip. This interaction underlies Kv-channel clustering, and is important for the proper assembly and functioning of the synapse. Here we demonstrate that the C-terminal segments of Kv channels adjacent to the PDZ-binding motif are intrinsically disordered. Phylogenetic analysis of the Kv channel family reveals a cluster of channel sequences belonging to three out of the four main channel families, for which an association is demonstrated between the presence of the consensus terminal PDZ-binding motif and the intrinsically disordered nature of the immediately adjacent C-terminal segment. Our observations, combined with a structural analogy to the N-terminal intra-molecular ball-and-chain mechanism for Kv channel inactivation, suggest that the C-terminal disordered segments of these channel families encode an inter-molecular fishing rod-like mechanism for K(+) channel binding to scaffold proteins.

Amino Acid Motifs↗
PubMed2006

Has the code for protein translocation been broken?

Polypeptides chains are segregated by the translocon channel into secreted or membrane-inserted proteins. Recent reports claim that an in vivo system has been used to break the "amino acid code" used by translocons to make the determination of protein type (i.e. secreted or membrane-inserted). However, the experimental setup used in these studies could have confused the derivation of this code, in particular for polar amino acids. These residues are likely to undergo stabilizing interactions with other protein components in the experiment, shielding them from direct contact with the inhospitable membrane. Hence, it is our view that the "code" for protein translocation has not yet been deciphered and that further experiments are required for teasing apart the various energetic factors contributing to protein translocation.

Animals↗
PubMed2006

Transmembrane protein structures without X-rays.

Transmembrane (TM) proteins constitute 15-30% of the genome, but <1% of the structures in the Protein Data Bank. This discrepancy is disturbing, and emphasizes that structure determination of TM proteins remains challenging. The challenge is greatest for proteins from eukaryotes, the structures of which remain intractable despite tremendous advances that have been made towards structure determination of bacterial TM proteins. Notably, >50% of the membrane protein families in eukaryotes lack bacterial homologs. Therefore, it is conceivable that many more years will elapse before high-resolution structures of eukaryotic TM proteins emerge. Until then, integrated approaches that combine biochemical and computational analyses with low-resolution structures are likely to have increasingly important roles in providing frameworks for the mechanistic understanding of membrane-protein structure and function.

Bacterial Proteins↗
PubMed2004

A Calpha model for the transmembrane alpha helices of gap junction intercellular channels.

Gap junction channels connect the cytoplasms of apposed cells via an intercellular conduit formed by the end-to-end docking of two hexameric hemichannels called connexons. We used electron cryomicroscopy to derive a three-dimensional density map at 5.7 angstroms in-plane and 19.8 angstroms vertical resolution, allowing us to identify the positions and tilt angles for the 24 alpha helices within each hemichannel. The four hydrophobic segments in connexin sequences were assigned to the alpha helices in the map based on biochemical and phylogenetic data. Analyses of evolutionary conservation and compensatory mutations in connexin evolution identified the packing interfaces between the helices. The final model, which specifies the coordinates of Calpha atoms in the transmembrane domain, provides a structural basis for understanding the different physiological effects of almost 30 mutations and polymorphisms in terms of structural deformations at the interfaces between helices, revealing an intimate connection between molecular structure and disease.

Amino Acid Sequence↗
PubMed2004

An automatic method for predicting transmembrane protein structures using cryo-EM and evolutionary data.

The transmembrane (TM) domains of many integral membrane proteins are composed of alpha-helix bundles. Structure determination at high resolution (<4 A) of TM domains is still exceedingly difficult experimentally. Hence, some TM-protein structures have only been solved at intermediate (5-10 A) or low (>10 A) resolutions using, for example, cryo-electron microscopy (cryo-EM). These structures reveal the packing arrangement of the TM domain, but cannot be used to determine the positions of individual amino acids. The observation that typically, the lipid-exposed faces of TM proteins are evolutionarily more variable and less charged than their core provides a simple rule for orienting their constituent helices. Based on this rule, we developed score functions and automated methods for orienting TM helices, for which locations and tilt angles have been determined using, e.g., cryo-EM data. The method was parameterized with the aim of retrieving the native structure of bacteriorhodopsin among near- and far-from-native templates. It was then tested on proteins that differ from bacteriorhodopsin in their sequences, architectures, and functions, such as the acetylcholine receptor and rhodopsin. The predicted structures were within 1.5-3.5 A from the native state in all cases. We conclude that the computational method can be used in conjunction with cryo-EM data to obtain approximate model structures of TM domains of proteins for which a sufficiently heterogeneous set of homologs is available. We also show that in those proteins in which relatively short loops connect neighboring helices, the scoring functions can discriminate between near- and far-from-native conformations even without the constraints imposed on helix locations and tilt angles that are derived from cryo-EM.

Algorithms↗
PubMed2004

Assigning transmembrane segments to helices in intermediate-resolution structures.

MOTIVATION: Transmembrane (TM) proteins that form alpha-helix bundles constitute approximately 50% of contemporary drug targets. Yet, it is difficult to determine their high-resolution (< 4 A) structures. Some TM proteins yield more easily to structure determination using cryo electron microscopy (cryo-EM), though this technique most often results in lower resolution structures, precluding an unambiguous assignment of TM amino acid sequences to the helices seen in the structure. We present computational tools for assigning the TM segments in the protein's sequence to the helices seen in cryo-EM structures. RESULTS: The method examines all feasible TM helix assignments and ranks each one based on a score function that was derived from loops in the structures of soluble alpha-helix bundles. A set of the most likely assignments is then suggested. We tested the method on eight TM chains of known structures, such as bacteriorhodopsin and the lactose permease. Our results indicate that many assignments can be rejected at the outset, since they involve the connection of pairs of remotely placed TM helices. The correct assignment received a high score, and was ranked highly among the remaining assignments. For example, in the lactose permease, which contains 12 TM helices, most of which are connected by short loops, only 12 out of 479 million assignments were found to be feasible, and the native one was ranked first. AVAILABILITY: The program and the non-redundant set of protein structures used here are available at http://www.cs.tau.ac.il/~angela

Amino Acid Sequence↗
PubMed2004

An evolutionarily conserved network of amino acids mediates gating in voltage-dependent potassium channels.

A novel sequence-analysis technique for detecting correlated amino acid positions in intermediate-size protein families (50-100 sequences) was developed, and applied to study voltage-dependent gating of potassium channels. Most contemporary methods for detecting amino acid correlations within proteins use very large sets of data, typically comprising hundreds or thousands of evolutionarily related sequences, to overcome the relatively low signal-to-noise ratio in the analysis of co-variations between pairs of amino acid positions. Such methods are impractical for voltage-gated potassium (Kv) channels and for many other protein families that have not yet been sequenced to that extent. Here, we used a phylogenetic reconstruction of paralogous Kv channels to follow the evolutionary history of every pair of amino acid positions within this family, thus increasing detection accuracy of correlated amino acids relative to contemporary methods. In addition, we used a bootstrapping procedure to eliminate correlations that were statistically insignificant. These and other measures allowed us to increase the method's sensitivity, and opened the way to reliable identification of correlated positions even in intermediate-size protein families. Principal-component analysis applied to the set of correlated amino acid positions in Kv channels detected a network of inter-correlated residues, a large fraction of which were identified as gating-sensitive upon mutation. Mapping the network of correlated residues onto the 3D structure of the Kv channel from Aeropyrum pernix disclosed correlations between residues in the voltage-sensor paddle and the pore region, including regions that are involved in the gating transition. We discuss these findings with respect to the evolutionary constraints acting on the channel's various domains. The software is available on our website

Amino Acid Sequence↗
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

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

pANT: a method for the pairwise assessment of nonfunctionalization times of processed pseudogenes.

We present a method for pairwise Assessment of Nonfunctionalization Times (pANT) in processed pseudogenes. Contrary to existing methods for estimating nonfunctionalization times, pANT utilizes previously calculated probabilities of nucleotide substitution as explicit rate measurements, rather than assume that the substitution rates are the same for all nucleotides. Thus, the method allows a more accurate computation of the time that has elapsed since the nonfunctionalization of a pseudogene. Whereas existing methods require the sequence of an orthologous functional gene, which is not always at hand, pANT only uses the pairwise alignment of the gene/pseudogene pair, thus expanding the range of problems that can be tackled. To estimate evolutionary times in nonfunctional sequences, pANT measures the differences in the pairwise alignment of a gene and its paralogous processed pseudogene, using only the first and second codon positions. It assumes that, because of functional constraints, these positions in the sequence of the functional homolog have not changed since the time of nonfunctionalization of the pseudogene. Hence, the sequence of the gene may be used as the ancestor of the pseudogene. We show that the method's reliance on a detailed substitution matrix, which is derived separately for each species, makes it more accurate than existing methods. We applied pANT to the case of the unitary alpha-1,3-galactosyltransferase human pseudogene and found that our estimate of the nonfunctionalization time was in agreement with that obtained by taxonomic and paleontological considerations pertaining to the divergence between platyrrhines (New World monkeys) and cattarhines (Old World monkeys).

Animals↗
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

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↗