Biomedical subjects

Dror Tobi

Publications and source records attributed to Dror Tobi.

7 recordsLinked to original sources

PubMed2006

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 A 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.

Antigens↗

PubMed2005

Structural changes involved in protein binding correlate with intrinsic motions of proteins in the unbound state.

Protein-protein binding usually involves structural changes that may extend beyond the rearrangements on a local scale, and cannot be explained by a classical lock-and-key mechanism. Several models have been advanced to explain the flexible binding of proteins such as the induced fit mechanism where the ligand is postulated to induce a conformational change at the interaction site upon binding, or the preexisting equilibrium hypothesis that assumes that protein samples an ensemble of conformations at equilibrium conditions and that the ligand binds selectively to an active conformation. We explored the equilibrium motions of proteins that exhibit relatively large (nonlocal) conformational changes upon protein binding using the Gaussian network model and the anisotropic network model of protein dynamics. For four complexes, LIR-1/HLA-A2, Actin/DNase I, CDK2/cyclin, and CDK6/p16(INK4a), the motions calculated for the monomer exhibiting the largest conformational change, in its unbound (free) form, correlate with the experimentally observed structural changes upon binding. This study emphasizes the preexisting equilibrium/conformational selection as a mechanism for protein-protein interaction and lends support the concept that proteins, in their native conformation, are predisposed to undergo conformational fluctuations that are relevant to, or even required for, their biological functions.

Actins↗

PubMed2004

Discovery and characterization of novel small molecule inhibitors of human Cdc25B dual specificity phosphatase.

Cdc25A and Cdc25B dual-specificity phosphatases are key regulators of cell cycle transition and proliferation. They have oncogenic properties and are overexpressed in many human tumors. Because selective Cdc25 phosphatase inhibitors would be valuable biological tools and possible therapeutic agents, we have assayed a small molecule library for in vitro inhibition of Cdc25. We now report the identification of two new structurally distinct classes of Cdc25 inhibitors with cellular activity. The cyclopentaquinoline 3a,4,5,9b-tetrahydro-3H-cyclopenta[c]quinoline-4,8-dicarboxylic acid (5661118) and the naphthofurandione 3-benzoyl-naphtho[1,2-b]furan-4,5-dione (5169131) had in vitro IC50 values of 2.5 to 11 microM against recombinant Cdc25 and were less potent inhibitors of other phosphatases. Unlike 5661118, 5169131 caused reversible inhibition of Cdc25B and displayed competitive inhibitor kinetics. No growth inhibitory activity was seen with 5661118, whereas 10 to 30 microM 5169131 caused G1/S and G2/M arrest. We also found that 5169131 inhibited human PC-3 prostate and MDA-MB-435 breast cancer cell proliferation. Concentration-dependent Tyr15 hyperphosphorylation was seen on cyclin-dependent kinase with a 1-h 5169131 treatment, consistent with Cdc25 inhibition. Cells resistant to DNA toposiomerase II inhibitors were as sensitive to 5169131 as parental cells, indicating that this quinone compound does not inhibit topoisomerase II in vivo. Molecular modeling was used to predict a potential interaction site between the inhibitor and Cdc25B and to provide insights as to the molecular origins of the experimental observations. Based on its kinetic profile and cellular activity, we suggest that 5169131 could be an excellent tool for further studies on the cellular roles of Cdc25.

Cell Cycle↗

PubMed2004

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 (20(4)) 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.

Algorithms↗

PubMed2003

Sequence-based design of kinase inhibitors applicable for therapeutics and target identification.

A platform for specifically modulating kinase-dependent signaling using peptides derived from the catalytic domain of the kinase is presented. This technology, termed KinAce, utilizes the canonical structure of protein kinases. The targeted regions (subdomain V and subdomains IX and X) are analyzed and their sequence, three-dimensional structure, and involvement in protein-protein interaction are highlighted. Short myristoylated peptides were derived from the target regions of the tyrosine kinases c-Kit and Lyn and the serine/threonine kinases 3-phosphoinositide-dependent kinase-1 (PDK1) and Akt/protein kinase B (PKB). For each kinase an active designer peptide is shown to selectively inhibit the signaling of the kinase from which it is derived, and to inhibit cancer cell proliferation in the micromolar range. This technology emerges as an applicable tool for deriving sequence-based selective inhibitors for a broad range of protein kinases as hits that may be further developed into drugs. Moreover, it enables identification of novel kinase targets for selected therapeutic indications as demonstrated in the KinScreen application.

3-Phosphoinositide-Dependent Protein Kinases↗

PubMed2003

Allosteric changes in protein structure computed by a simple mechanical model: hemoglobin T<-->R2 transition.

Information on protein dynamics has been usually inferred from spectroscopic studies of parts of the proteins, or indirectly from the comparison of the conformations assumed in the presence of different substrates or ligands. While molecular simulations also provide information on protein dynamics, they usually suffer from incomplete sampling of conformational space, and become prohibitively expensive when exploring the collective dynamics of large macromolecular structures. Here, we explore the dynamics of a well-studied allosteric protein, hemoglobin (Hb), to show that a simple mechanical model based on Gaussian fluctuations of residues can efficiently predict the transition between the tense (T, unliganded) and relaxed (R or R2, O(2) or CO-bound) forms of Hb. The passage from T into R2 is shown to be favored by the global mode of motion, which, in turn is driven by entropic effects. The major difference between the dynamics of the T and R2 forms is the loss of the hinge-bending role of alpha(1)-beta(2) (or alpha(2)-beta(1)) interfacial residues at alpha Phe36-His45 and beta Thr87-Asn102 in the R2 form, which implies a decreased cooperativity in the higher affinity (R2) form of Hb, consistent with many experimental studies. The involvement of the proximal histidine beta His92 in this hinge region suggests that the allosteric propagation of the local structural changes (induced upon O(2) binding) into global ones occur via hinge regions. This is the first demonstration that there is an intrinsic tendency of Hb to undergo T-->R2 transition, induced by purely elastic forces of entropic origin that are uniquely defined for the particular contact topology of the T form.

Animals↗

PubMed2003

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.

Computer Simulation↗