Blind prediction of interfacial water positions in CAPRI.

We report the first assessment of blind predictions of water positions at protein-protein interfaces, performed as part of the critical assessment of predicted interactions (CAPRI) community-wide experiment. Groups submitting docking predictions for the complex of the DNase domain of colicin E2 and Im2 immunity protein (CAPRI Target 47), were invited to predict the positions of interfacial water molecules using the method of their choice. The predictions-20 groups submitted a total of 195 models-were assessed by measuring the recall fraction of water-mediated protein contacts. Of the 176 high- or medium-quality docking models-a very good docking performance per se-only 44% had a recall fraction above 0.3, and a mere 6% above 0.5. The actual water positions were in general predicted to an accuracy level no better than 1.5 Å, and even in good models about half of the contacts represented false positives. This notwithstanding, three hotspot interface water positions were quite well predicted, and so was one of the water positions that is believed to stabilize the loop that confers specificity in these complexes. Overall the best interface water predictions was achieved by groups that also produced high-quality docking models, indicating that accurate modelling of the protein portion is a determinant factor. The use of established molecular mechanics force fields, coupled to sampling and optimization procedures also seemed to confer an advantage. Insights gained from this analysis should help improve the prediction of protein-water interactions and their role in stabilizing protein complexes.

Diameter-selective dispersion of carbon nanotubes by ß-lactoglobulin whey protein.

The ß-lactoglobulin (ß-LG) protein was discovered to be an efficient and selective dispersant for carbon nanotubes (CTNs) with certain diameters. A dispersion process of CTNs by the ß-LG was studied, focusing on the relationships between the surface curvature of the CNT and the ß-LG's efficiency in dispersing them, using cryogenic-transmission electron microscopy (cryo-TEM) and optical spectroscopy. Plausible binding sites of the ß-LG, responsible for the interaction of the protein with CNTs of various diameters (surface curvatures) were also investigated and were found to be in good agreement with corresponding docking calculations.

Structure-based in silico identification of ubiquitin-binding domains provides insights into the ALIX-V:ubiquitin complex and retrovirus budding.

The ubiquitylation signal promotes trafficking of endogenous and retroviral transmembrane proteins. The signal is decoded by a large set of ubiquitin (Ub) receptors that tether Ub-binding domains (UBDs) to the trafficking machinery. We developed a structure-based procedure to scan the protein data bank for hidden UBDs. The screen retrieved many of the known UBDs. Intriguingly, new potential UBDs were identified, including the ALIX-V domain. Pull-down, cross-linking and E3-independent ubiquitylation assays biochemically corroborated the in silico findings. Guided by the output model, we designed mutations at the postulated ALIX-V:Ub interface. Biophysical affinity measurements using microscale-thermophoresis of wild-type and mutant proteins revealed some of the interacting residues of the complex. ALIX-V binds mono-Ub with a K(d) of 119 µM. We show that ALIX-V oligomerizes with a Hill coefficient of 5.4 and IC(50) of 27.6 µM and that mono-Ub induces ALIX-V oligomerization. Moreover, we show that ALIX-V preferentially binds K63 di-Ub compared with mono-Ub and K48 di-Ub. Finally, an in vivo functionality assay demonstrates the significance of ALIX-V:Ub interaction in equine infectious anaemia virus budding. These results not only validate the new procedure, but also demonstrate that ALIX-V directly interacts with Ub in vivo and that this interaction can influence retroviral budding.

Computer-based identification of a novel LIMK1/2 inhibitor that synergizes with salirasib to destabilize the actin cytoskeleton.

Neurofibromin regulates cell motility via three distinct GTPase pathways acting through two different domains, the Ras GTPase-activating protein-related domain (GRD) and the pre-GRD domain. First, the GRD domain inhibits Ras-dependent changes in cell motility through the mitogen activated protein cascade. Second, it also regulates Rho-dependent (Ras-independent) changes by activating LIM kinase 2 (LIMK2), an enzyme that phosphorylates and inactivates cofilin (an actin-depolymerizing factor). Third, the pre-GRD domain acts through the Rac1 GTPase, that activate the P21 activated kinase 1 (PAK1)-LIMK1-cofilin pathway. We employed molecular modeling to identify a novel inhibitor of LIMK1/2. The active sites of an ephrin-A receptor (EphA3) and LIMK2 showed marked similarity (60%). On testing a known inhibitor of EphA3, we found that it fits to the LIMK1/2-ATP binding site and to the latter's substrate-binding pockets. We identified a similar compound, T56-LIMKi, and found that it inhibits LIMK1/2 kinase activities. It blocked the phosphorylation of cofilin which led to actin severance and inhibition of tumor cell migration, tumor cell growth, and anchorage-independent colony formation in soft agar. Because modulation of LIMK by neurofibromin is not affected by the Ras inhibitor Salirasib, we examined the combined effect of Salirasib and T56-LIMKi each of which can affect cell motility by a distinct pathway. We found that their combined action on cell proliferation and stress-fiber formation in neurofibromin-deficient cells was synergistic. We suggest that this drug combination may be developed for treatment of neurofibromatosis and cancer.

SymmRef: a flexible refinement method for symmetric multimers.

Symmetric protein complexes are abundant in the living cell. Predicting their atomic structure can shed light on the mechanism of many important biological processes. Symmetric docking methods aim to predict the structure of these complexes given the unbound structure of a single monomer, or its model. Symmetry constraints reduce the search-space of these methods and make the prediction easier compared to asymmetric protein-protein docking. However, the challenge of modeling the conformational changes that the monomer might undergo is a major obstacle. In this article, we present SymmRef, a novel method for refinement and reranking of symmetric docking solutions. The method models backbone and side-chain movements and optimizes the rigid-body orientations of the monomers. The backbone movements are modeled by normal modes minimization and the conformations of the side-chains are modeled by selecting optimal rotamers. Since solved structures of symmetric multimers show asymmetric side-chain conformations, we do not use symmetry constraints in the side-chain optimization procedure. The refined models are re-ranked according to an energy score. We tested the method on a benchmark of unbound docking challenges. The results show that the method significantly improves the accuracy and the ranking of symmetric rigid docking solutions. SymmRef is available for download at http:// bioinfo3d.cs.tau.ac.il/SymmRef/download.html.