Discovery of Rab1 binding sites using an ensemble of clustering methods

Focus on the Small GTPase Rab1: A Key Player in the Pathogenesis of Parkinson's Disease

Parkinson’s disease (PD) is the second most frequent neurodegenerative disease. It is characterized by the loss of dopaminergic neurons in the substantia nigra and the formation of large aggregates in the survival neurons called Lewy bodies, which mainly contain α-synuclein (α-syn). The cause of cell death is not known but could be due to mitochondrial dysfunction, protein homeostasis failure, and alterations in the secretory/endolysosomal/autophagic pathways. Survival nigral neurons overexpress the small GTPase Rab1. This protein is considered a housekeeping Rab that is necessary to support the secretory pathway, the maintenance of the Golgi complex structure, and the regulation of macroautophagy from yeast to humans. It is also involved in signaling, carcinogenesis, and infection for some pathogens. It has been shown that it is directly linked to the pathogenesis of PD and other neurodegenerative diseases. It has a protective effect against α–σψν toxicity and has recently been shown to be a substrate of LRRK2, which is the most common cause of familial PD and the risk of sporadic disease. In this review, we analyze the key aspects of Rab1 function in dopamine neurons and its implications in PD neurodegeneration/restauration. The results of the current and former research support the notion that this GTPase is a good candidate for therapeutic strategies.

Structural dynamics and allostery of Rab proteins: strategies for drug discovery and design

Rab proteins represent the largest family of the Rab superfamily guanosine triphosphatase (GTPase). Aberrant human Rab proteins are associated with multiple diseases, including cancers and neurological disorders. Rab subfamily members display subtle conformational variations that render specificity in their physiological functions and can be targeted for subfamily-specific drug design. However, drug discovery efforts have not focused much on targeting Rab allosteric non-nucleotide binding sites which are subjected to less evolutionary pressures to be conserved, hence are likely to offer subfamily specificity and may be less prone to undesirable off-target interactions and side effects. To discover druggable allosteric binding sites, Rab structural dynamics need to be first incorporated using multiple experimentally and computationally obtained structures. The high-dimensional structural data may necessitate feature extraction methods to identify manageable representative structures for subsequent analyses. We have detailed state-of-the-art computational methods to (i) identify binding sites using data on sequence, shape, energy, etc., (ii) determine the allosteric nature of these binding sites based on structural ensembles, residue networks and correlated motions and (iii) identify small molecule binders through structure- and ligand-based virtual screening. To benefit future studies for targeting Rab allosteric sites, we herein detail a refined workflow comprising multiple available computational methods, which have been successfully used alone or in combinations. This workflow is also applicable for drug discovery efforts targeting other medically important proteins. Depending on the structural dynamics of proteins of interest, researchers can select suitable strategies for allosteric drug discovery and design, from the resources of computational methods and tools enlisted in the workflow.

Identification of native protein structures captured by principal interactions

Background

Evaluation of protein structure is based on trustworthy potential function. The total potential of a protein structure is approximated as the summation of all pair-wise interaction potentials. Knowledge-based potentials (KBP) are one type of potential functions derived by known experimentally determined protein structures. Although several KBP functions with different methods have been introduced, the key interactions that capture the total potential have not studied yet.

Results

In this study, we seek the interaction types that preserve as much of the total potential as possible. We employ a procedure based on the principal component analysis (PCA) to extract the significant and key interactions in native protein structures. We call these interactions as principal interactions and show that the results of the model that considers only these interactions are very close to the full interaction model that considers all interactions in protein fold recognition. In fact, the principal interactions maintain the discriminative power of the full interaction model. This method was evaluated on 3 KBPs with different contact definitions and thresholds of distance and revealed that their corresponding principal interactions are very similar and have a lot in common. Additionally, the principal interactions consisted of 20 % of the full interactions on average, and they are between residues, which are considered important in protein folding.

Conclusions

This work shows that all interaction types are not equally important in discrimination of native structure. The results of the reduced model based on principal interactions that were very close to the full interaction model suggest that a new strategy is needed to capture the role of remaining interactions (non-principal interactions) to improve the power of knowledge-based potential functions.

Allosteric binding sites in Rab11 for potential drug candidates

Rab11 is an important protein subfamily in the RabGTPase family. These proteins physiologically function as key regulators of intracellular membrane trafficking processes. Pathologically, Rab11 proteins are implicated in many diseases including cancers, neurodegenerative diseases and type 2 diabetes. Although they are medically important, no previous study has found Rab11 allosteric binding sites where potential drug candidates can bind to. In this study, by employing multiple clustering approaches integrating principal component analysis, independent component analysis and locally linear embedding, we performed structural analyses of Rab11 and identified eight representative structures. Using these representatives to perform binding site mapping and virtual screening, we identified two novel binding sites in Rab11 and small molecules that can preferentially bind to different conformations of these sites with high affinities. After identifying the binding sites and the residue interaction networks in the representatives, we computationally showed that these binding sites may allosterically regulate Rab11, as these sites communicate with switch 2 region that binds to GTP/GDP. These two allosteric binding sites in Rab11 are also similar to two allosteric pockets in Ras that we discovered previously.

Structural analysis of protein tyrosine phosphatase 1B reveals potentially druggable allosteric binding sites

Catalytic proteins such as human protein tyrosine phosphatase 1B (PTP1B), with conserved and highly polar active sites, warrant the discovery of druggable nonactive sites, such as allosteric sites, and potentially, therapeutic small molecules that can bind to these sites. Catalyzing the dephosphorylation of numerous substrates, PTP1B is physiologically important in intracellular signal transduction pathways in diverse cell types and tissues. Aberrant PTP1B is associated with obesity, diabetes, cancers, and neurodegenerative disorders. Utilizing clustering methods (based on root mean square deviation, principal component analysis, nonnegative matrix factorization, and independent component analysis), we have examined multiple PTP1B structures. Using the resulting representative structures in different conformational states, we determined consensus clustroids and used them to identify both known and novel binding sites, some of which are potentially allosteric. We report several lead compounds that could potentially bind to the novel PTP1B binding sites and can be further optimized. Considering the possibility for drug repurposing, we discovered homologous binding sites in other proteins, with ligands that could potentially bind to the novel PTP1B binding sites.

Molecular Dynamics Simulation and Prediction of Druggable Binding Sites

Binding site identification and druggability evaluation are two essential steps in structure-based drug design. A druggable binding site tends to have high binding affinity to drug-like molecules. Predicting such sites can have a significant impact on a drug design campaign. This chapter focuses on summarizing the different methods that are used to predict druggable binding sites. The chapter also discusses the importance of including protein flexibility in the search process and the use of molecular dynamics simulations to address this aspect. Case studies from the literature are also summarized and discussed. We hope that this chapter would provide an overview on the different methods employed in binding site identification evaluation.