Dror Tobi

Distance-Dependent, Pair Potential for Protein Folding: Results From Linear Optimization

The results of an optimization of a folding potential are reported. The complete energy function is modeled as a sum of pairwise interactions with a flexible functional form. The relevant distance between two amino acids (2 − 9 Å) is divided into 13 intervals, and the energy of each interval is optimized independently. We show, in accord with a previous publication (Tobi et al., Proteins 2000;40:71–85) that it is impossible to find a pair potential with the above flexible form that recognizes all native folds. Nevertheless, a potential that rates correctly a subset of the decoy structures was constructed and optimized. The resulting potential is compared with a distance‐dependent statistical potential of Bahar and Jernigan. It is further tested against decoy structures that were created in the Levitts group. On average, the new potential places native shapes lower in energy and provides higher Z scores than other potentials. Proteins 2000;41:40–46.

On the design and analysis of protein folding potentials

Pairwise interaction models to recognize native folds are designed and analyzed. Different sets of parameters are considered but the focus was on 20 × 20 contact matrices. Simultaneous solution of inequalities and minimization of the variance of the energy find matrices that recognize exactly the native folds of 572 sequences and structures from the protein data bank (PDB). The set includes many homologous pairs, which present a difficult recognition problem. Significant recognition ability is recovered with a small number of parameters (e.g., the H/P model). However, full recognition requires a complete set of amino acids. In addition to structures from the PDB, a folding program (MONSSTER) was used to generate decoy structures for 75 proteins. It is impossible to recognize all the native structures of the extended set by contact potentials. We therefore searched for a new functional form. An energy function U, which is based on a sum of general pairwise interactions limited to a resolution of 1 angstrom, is considered. This set was infeasible too. We therefore conjecture that it is not possible to find a folding potential, resolved to 1 angstrom, which is a sum of pair interactions. Proteins 2000;40:71–85.

The dominant interaction between peptide and urea is electrostatic in nature: A molecular dynamics simulation study

The conformational equilibrium of a blocked valine peptide in water and aqueous urea solution is studied using molecular dynamics simulations. Pair correlation functions indicate enhanced concentration of urea near the peptide. Stronger hydrogen bonding of urea–peptide compared to water–peptide is observed with preference for helical conformation. The potential of mean force, computed using umbrella sampling, shows only small differences between urea and water solvation that are difficult to quantify. The changes in solvent structure around the peptide are explained by favorable electrostatic interactions (hydrogen bonds) of urea with the peptide backbone. There is no evidence for significant changes in hydrophobic interactions in the two conformations of the peptide in urea solution. Our simulations suggest that urea denatures proteins by preferentially forming hydrogen bonds to the peptide backbone, reducing the barrier for exposing protein residues to the solvent, and reaching the unfolded state.

Synaptotagmin restores kinetic properties of a syntaxin‐associated N‐type voltage sensitive calcium channel

© 1997 Federation of European Biochemical Societies.

The α2/δ subunit of voltage sensitive Ca2+ channels is a single transmembrane extracellular protein which is involved in regulated secretion

The membrane topology of α2/δ subunit was investigated utilizing electrophysiological functional assay and specific anti‐α2 antibodies. (a) cRNA encoding a deleted α2/δ subunit was coinjected with α1C subunit of the L‐type calcium channel into Xenopus oocytes. The truncated form, lacking the third putative TM domain (α2/δΔTMIII), failed to amplify the expressed inward currents, normally induced by α 1c coinjected with intact α2/δ subunit. Western blot analysis of α2/δΔTMIII shows the appearance of a degraded α2 protein and no expression of the full‐size two‐TM truncated‐protein. The improper processing of α2/δΔTMIII suggests that the α2/δ is a single TM domain protein and the TM region is positioned at the δ subunit. (b) External application of anti‐α2 antibodies, prepared for an epitope within the alternatively spliced and ‘intracellular’ region, inhibits depolarization induced secretion in PC12, further supporting an external location of the α2 subunit and establishing δ subunit as the only membrane anchor for the extracellular α2 subunit.

Optimal design of protein docking potentials : Efficiency and limitations

Protein–protein docking is a challenging computational problem in functional genomics, particularly when one or both proteins undergo conformational change(s) upon binding. The major challenge is to define scoring function soft enough to tolerate these changes and specific enough to distinguish between near‐native and “misdocked” conformations. Using a linear programming technique, we derived protein docking potentials (PDPs) that comply with this requirement. We considered a set of 63 nonredundant complexes to this aim, and generated 400,000 putative docked complexes (decoys) based on shape complementarity criterion for each complex. The PDPs were required to yield for the native (correctly docked) structure a potential energy lower than those of all the nonnative (misdocked) structures. The energy constraints applied to all complexes led to ca. 25 million inequalities, the simultaneous solution of which yielded an optimal set of PDPs that discriminated the correctly docked (up to 4.0 Å root‐mean‐square deviation from known complex structure) structure among the 85 top‐ranking (0.02%) decoys in 59/63 examined bound–bound cases. The high performance of the potentials was further verified in jackknife tests and by ranking putative docked conformation submitted to CAPRI. In addition to their utility in identifying correctly folded complexes, the PDPs reveal biologically meaningful features that distinguish docking potentials from folding potentials. Proteins 2006.

Docking of antizyme to ornithine decarboxylase and antizyme inhibitor using experimental mutant and double-mutant cycle data.

Antizyme (Az) is a highly conserved key regulatory protein bearing a major role in regulating polyamine levels in the cell. It has the ability to bind and inhibit ornithine decarboxylase (ODC), targeting it for degradation. Az inhibitor (AzI) impairs the activity of Az. In this study, we mapped the binding sites of ODC and AzI on Az using Ala scan mutagenesis and generated models of the two complexes by constrained computational docking. In order to scan a large number of mutants in a short time, we developed a workflow combining high-throughput mutagenesis, small-scale parallel partial purification of His-tagged proteins and their immobilization on a tris-nitrilotriacetic-acid-coated surface plasmon resonance chip. This combination of techniques resulted in a significant reduction in time for production and measurement of large numbers of mutant proteins. The data-driven docking results suggest that both proteins occupy the same binding site on Az, with Az binding within a large groove in AzI and ODC. However, single-mutant data provide information concerning the location of the binding sites only, not on their relative orientations. Therefore, we generated a large number of double-mutant cycles between residues on Az and ODC and used the resulting interaction energies to restrict docking. The model of the complex is well defined and accounts for the mutant data generated here, and previously determined biochemical data for this system. Insights on the structure and function of the complexes, as well as general aspects of the method, are discussed.

A Sequence Alignment-Independent Method for Protein Classification

Annotation of the rapidly accumulating body of sequence data relies heavily on the detection of remote homologues and functional motifs in protein families. The most popular methods rely on sequence alignment. These include programs that use a scoring matrix to compare the probability of a potential alignment with random chance and programs that use curated multiple alignments to train profile hidden Markov models (HMMs). Related approaches depend on bootstrapping multiple alignments from a single sequence. However, alignment-based programs have limitations. They make the assumption that contiguity is conserved between homologous segments, which may not be true in genetic recombination or horizontal transfer. Alignments also become ambiguous when sequence similarity drops below 40%. This has kindled interest in classification methods that do not rely on alignment. An approach to classification without alignment based on the distribution of contiguous sequences of four amino acids (4-grams) was developed. Interest in 4-grams stemmed from the observation that almost all theoretically possible 4-grams (204) occur in natural sequences and the majority of 4-grams are uniformly distributed. This implies that the probability of finding identical 4-grams by random chance in unrelated sequences is low. A Bayesian probabilistic model was developed to test this hypothesis. For each protein family in Pfam-A and PIR-PSD, a feature vector called a probe was constructed from the set of 4-grams that best characterised the family. In rigorous jackknife tests, unknown sequences from Pfam-A and PIR-PSD were compared with the probes for each family. A classification result was deemed a true positive if the probe match with the highest probability was in first place in a rank-ordered list. This was achieved in 70% of cases. Analysis of false positives suggested that the precision might approach 85% if selected families were clustered into subsets. Case studies indicated that the 4-grams in common between an unknown and the best matching probe correlated with functional motifs from PRINTS. The results showed that remote homologues and functional motifs could be identified from an analysis of 4-gram patterns.

Dynamics alignment: Comparison of protein dynamics in the scop database

A novel methodology for comparison of protein dynamics is presented. Protein dynamics is calculated using the Gaussian network model and the modes of motion are globally aligned using the dynamic programming algorithm of Needleman and Wunsch, commonly used for sequence alignment. The alignment is fast and can be used to analyze large sets of proteins. The methodology is applied to the four major classes of the SCOP database: “all alpha proteins,” “all beta proteins,” “alpha and beta proteins,” and “alpha/beta proteins”. We show that different domains may have similar global dynamics. In addition, we report that the dynamics of “all alpha proteins” domains are less specific to structural variations within a given fold or superfamily compared with the other classes. We report that domain pairs with the most similar and the least similar global dynamics tend to be of similar length. The significance of the methodology is that it suggests a new and efficient way of mapping between the global structural features of protein families/subfamilies and their encoded dynamics. Proteins 2012;

Toward the development of a novel non-RGD cyclic peptide drug conjugate for treatment of human metastatic melanoma

The newly discovered short (9 amino acid) non-RGD S-S bridged cyclic peptide ALOS-4 (H-cycl(Cys-Ser-Ser-Ala-Gly-Ser-Leu-Phe-Cys)-OH), which binds to integrin αvβ3 is investigated as peptide carrier for targeted drug delivery against human metastatic melanoma. ALOS4 binds specifically the αvβ3 overexpressing human metastatic melanoma WM-266-4 cell line both in vitro and in ex vivo assays. Coupling ALOS4 to the topoisomerase I inhibitor Camptothecin (ALOS4-CPT) increases the cytotoxicity of CPT against human metastatic melanoma cells while reduces dramatically the cytotoxicity against non-cancerous cells as measured by the levels of γH2A.X, active caspase 3 and cell viability. Moreover, conjugating ALOS4 to CPT even increases the chemo-stability of CPT under physiological pH. Bioinformatic analysis using Rosetta platform revealed potential docking sites of ALOS4 on the αvβ3 integrin which are distinct from the RGD binding sites. We propose to use this specific non-RGD cyclic peptide as the therapeutic carrier for conjugation of drugs in order to improve efficacy and reduce toxicity of currently available treatments of human malignant melanoma.