Interfaces in Molecular Docking

An automated decision-tree approach to predicting protein interaction hot spots

Protein-protein interactions can be altered by mutating one or more "hot spots," the subset of residues that account for most of the interface's binding free energy. The identification of hot spots requires a significant experimental effort, highlighting the practical value of hot spot predictions. We present two knowledge-based models that improve the ability to predict hot spots: K-FADE uses shape specificity features calculated by the Fast Atomic Density Evaluation (FADE) program, and K-CON uses biochemical contact features. The combined K-FADE/CON (KFC) model displays better overall predictive accuracy than computational alanine scanning (Robetta-Ala). In addition, because these methods predict different subsets of known hot spots, a large and significant increase in accuracy is achieved by combining KFC and Robetta-Ala. The KFC analysis is applied to the calmodulin (CaM)/smooth muscle myosin light chain kinase (smMLCK) interface, and to the bone morphogenetic protein-2 (BMP-2)/BMP receptor-type I (BMPR-IA) interface. The results indicate a strong correlation between KFC hot spot predictions and mutations that significantly reduce the binding affinity of the interface.

Virtual Screening, Molecular Interaction Field, Molecular Dynamics, Docking, Density Functional, and ADMET Properties of Novel AChE Inhibitors in Alzheimer's Disease

Alzheimer's disease (AD) affects approximately 10% of the world's population with 65 years of age, being the most common form of dementia in adults and is characterized by senile plaquets and cholinergic deficits. Many drugs currently used for the treatment of the AD are based on the improvement of cholinergic neurotransmission achieved by Acetylcholinesterase (AChE) inhibition, the enzyme responsible for acetylcholine hydrolysis. We have focused in this work on the usage of computer-aided molecular design by virtual screening, molecular dynamics with implicit and explicit water solvation, density functional, molecular interaction field studies, docking procedures, ADMET predictions in order to propose novel potential AChE inhibitor for the treatment of Alzheimer's disease.

TRP1 interacting PDZ-domain protein GIPC forms oligomers and is localized to intracellular vesicles in human melanocytes

PDZ proteins coordinate assembly of protein complexes that participate in diverse biological processes. GIPC is a multifunctional PDZ protein that interacts with several soluble and membrane proteins. Unlike most PDZ proteins, GIPC contains single PDZ domain and the mechanisms by which GIPC mediates its actions remain unclear. We investigated the possibility that in lieu of multiple PDZ domains, GIPC forms multimers. Here, we demonstrate that GIPC can bind to itself and that the PDZ domain is involved in GIPC-GIPC interaction. Gel filtration, sucrose gradient centrifugation and chemical cross-linking showed that whereas bulk of cytosolic GIPC was present as monomer, oligomers with an estimated molecular mass corresponding to GIPC homotrimer were readily detectable in the membrane fraction. Modeling of GIPC PDZ domain showed feasibility of trimerization. Immunogold electron microscopy showed that GIPC is present in clusters near vesicles. Our data suggest that oligomers of GIPC mediate its functions in melanocytes.

Disruption of shape-complementarity markers to create cytotoxic variants of ribonuclease A

Onconase (ONC), an amphibian member of the bovine pancreatic ribonuclease A (RNase A) superfamily, is in phase III clinical trials as a treatment for malignant mesothelioma. RNase A is a far more efficient catalyst of RNA cleavage than ONC but is not cytotoxic. The innate ability of ONC to evade the cytosolic ribonuclease inhibitor protein (RI) is likely to be a primary reason for its cytotoxicity. In contrast, the non-covalent interaction between RNase A and RI is one of the strongest known, with the RI.RNase A complex having a K(d) value in the femtomolar range. Here, we report on the use of the fast atomic density evaluation (FADE) algorithm to identify regions in the molecular interface of the RI.RNase A complex that exhibit a high degree of geometric complementarity. Guided by these "knobs" and "holes", we designed variants of RNase A that evade RI. The D38R/R39D/N67R/G88R substitution increased the K(d) value of the pRI.RNase A complex by 20 x 10(6)-fold (to 1.4 microM) with little change to catalytic activity or conformational stability. This and two related variants of RNase A were more toxic to human cancer cells than was ONC. Notably, these cytotoxic variants exerted their toxic activity on cancer cells selectively, and more selectively than did ONC. Substitutions that further diminish affinity for RI (which has a cytosolic concentration of 4 microM) are unlikely to produce a substantial increase in cytotoxic activity. These results demonstrate the utility of the FADE algorithm in the examination of protein-protein interfaces and represent a landmark towards the goal of developing chemotherapeutics based on mammalian ribonucleases.