Nir Ben-Tal

ConSurf 2005: the projection of evolutionary conservation scores of residues on protein structures

Key amino acid positions that are important for maintaining the 3D structure of a protein and/or its function(s), e.g. catalytic activity, binding to ligand, DNA or other proteins, are often under strong evolutionary constraints. Thus, the biological importance of a residue often correlates with its level of evolutionary conservation within the protein family. ConSurf () is a web-based tool that automatically calculates evolutionary conservation scores and maps them on protein structures via a user-friendly interface. Structurally and functionally important regions in the protein typically appear as patches of evolutionarily conserved residues that are spatially close to each other. We present here version 3.0 of ConSurf. This new version includes an empirical Bayesian method for scoring conservation, which is more accurate than the maximum-likelihood method that was used in the earlier release. Various additional steps in the calculation can now be controlled by a number of advanced options, thus further improving the accuracy of the calculation. Moreover, ConSurf version 3.0 also includes a measure of confidence for the inferred amino acid conservation scores.

Residue frequencies and pairing preferences at protein–protein interfaces

We used a nonredundant set of 621 protein–protein interfaces of known high‐resolution structure to derive residue composition and residue–residue contact preferences. The residue composition at the interfaces, in entire proteins and in whole genomes correlates well, indicating the statistical strength of the data set. Differences between amino acid distributions were observed for interfaces with buried surface area of less than 1,000 Å2 versus interfaces with area of more than 5,000 Å2. Hydrophobic residues were abundant in large interfaces while polar residues were more abundant in small interfaces. The largest residue–residue preferences at the interface were recorded for interactions between pairs of large hydrophobic residues, such as Trp and Leu, and the smallest preferences for pairs of small residues, such as Gly and Ala. On average, contacts between pairs of hydrophobic and polar residues were unfavorable, and the charged residues tended to pair subject to charge complementarity, in agreement with previous reports. A bootstrap procedure, lacking from previous studies, was used for error estimation. It showed that the statistical errors in the set of pairing preferences are generally small; the average standard error is ≈0.2, i.e., about 8% of the average value of the pairwise index (2.9). However, for a few pairs (e.g., Ser–Ser and Glu–Asp) the standard error is larger in magnitude than the pairing index, which makes it impossible to tell whether contact formation is favorable or unfavorable. The results are interpreted using physicochemical factors and their implications for the energetics of complex formation and for protein docking are discussed. Proteins 2001;43:89–102.

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.

ConSurf 2016: an improved methodology to estimate and visualize evolutionary conservation in macromolecules

The degree of evolutionary conservation of an amino acid in a protein or a nucleic acid in DNA/RNA reflects a balance between its natural tendency to mutate and the overall need to retain the structural integrity and function of the macromolecule. The ConSurf web server (http://consurf.tau.ac.il), established over 15 years ago, analyses the evolutionary pattern of the amino/nucleic acids of the macromolecule to reveal regions that are important for structure and/or function. Starting from a query sequence or structure, the server automatically collects homologues, infers their multiple sequence alignment and reconstructs a phylogenetic tree that reflects their evolutionary relations. These data are then used, within a probabilistic framework, to estimate the evolutionary rates of each sequence position. Here we introduce several new features into ConSurf, including automatic selection of the best evolutionary model used to infer the rates, the ability to homology-model query proteins, prediction of the secondary structure of query RNA molecules from sequence, the ability to view the biological assembly of a query (in addition to the single chain), mapping of the conservation grades onto 2D RNA models and an advanced view of the phylogenetic tree that enables interactively rerunning ConSurf with the taxa of a sub-tree.

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.

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.

Mutations in the gene encoding pejvakin, a newly identified protein of the afferent auditory pathway, cause DFNB59 auditory neuropathy

Auditory neuropathy is a particular type of hearing impairment in which neural transmission of the auditory signal is impaired, while cochlear outer hair cells remain functional. Here we report on DFNB59, a newly identified gene on chromosome 2q31.1–q31.3 mutated in four families segregating autosomal recessive auditory neuropathy. DFNB59 encodes pejvakin, a 352-residue protein. Pejvakin is a paralog of DFNA5, a protein of unknown function also involved in deafness. By immunohistofluorescence, pejvakin is detected in the cell bodies of neurons of the afferent auditory pathway. Furthermore, Dfnb59 knock-in mice, homozygous for the R183W variant identified in one DFNB59 family, show abnormal auditory brainstem responses indicative of neuronal dysfunction along the auditory pathway. Unlike previously described sensorineural deafness genes, all of which underlie cochlear cell pathologies, DFNB59 is the first human gene implicated in nonsyndromic deafness due to a neuronal defect.

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.

The ConSurf-DB: pre-calculated evolutionary conservation profiles of protein structures.

ConSurf-DB is a repository for evolutionary conservation analysis of the proteins of known structures in the Protein Data Bank (PDB). Sequence homologues of each of the PDB entries were collected and aligned using standard methods. The evolutionary conservation of each amino acid position in the alignment was calculated using the Rate4Site algorithm, implemented in the ConSurf web server. The algorithm takes into account the phylogenetic relations between the aligned proteins and the stochastic nature of the evolutionary process explicitly. Rate4Site assigns a conservation level for each position in the multiple sequence alignment using an empirical Bayesian inference. Visual inspection of the conservation patterns on the 3D structure often enables the identification of key residues that comprise the functionally important regions of the protein. The repository is updated with the latest PDB entries on a monthly basis and will be rebuilt annually. ConSurf-DB is available online at http://consurfdb.tau.ac.il/

The ConSurf‐HSSP database: The mapping of evolutionary conservation among homologs onto PDB structures

The HSSP (Homology‐Derived Secondary Structure of Proteins) database provides multiple sequence alignments (MSAs) for proteins of known three‐dimensional (3D) structure in the Protein Data Bank (PDB). The database also contains an estimate of the degree of evolutionary conservation at each amino acid position. This estimate, which is based on the relative entropy, correlates with the functional importance of the position; evolutionarily conserved positions (i.e., positions with limited variability and low entropy) are occasionally important to maintain the 3D structure and biological function(s) of the protein. We recently developed the Rate4Site algorithm for scoring amino acid conservation based on their calculated evolutionary rate. This algorithm takes into account the phylogenetic relationships between the homologs and the stochastic nature of the evolutionary process. Here we present the ConSurf‐HSSP database of Rate4Site estimates of the evolutionary rates of the amino acid positions, calculated using HSSPs MSAs. The database provides precalculated evolutionary rates for nearly all of the PDB. These rates are projected, using a color code, onto the protein structure, and can be viewed online using the ConSurf server interface. To exemplify the database, we analyzed in detail the conservation pattern obtained for pyruvate kinase and compared the results with those observed using the relative entropy scores of the HSSP database. It is reassuring to know that the main functional region of the enzyme is detectable using both conservation scores. Interestingly, the ConSurf‐HSSP calculations mapped additional functionally important regions, which are moderately conserved and were overlooked by the original HSSP estimate. The ConSurf‐HSSP database is available online (http://consurf‐hssp.tau.ac.il). Proteins 2005.