Ensemble-based modeling and rigidity decomposition of allosteric interaction networks and communication pathways in cyclin-dependent kinases: Differentiating kinase clients of the Hsp90-Cdc37 chaperone

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

Conformational Dynamics and Mechanisms of Client Protein Integration into the Hsp90 Chaperone Controlled by Allosteric Interactions of Regulatory Switches: Perturbation-Based Network Approach for Mutational Profiling of the Hsp90 Binding and Allostery

Understanding the allosteric mechanisms of the Hsp90 chaperone interactions with cochaperones and client protein clientele is fundamental to dissect activation and regulation of many proteins. In this work, atomistic simulations are combined with perturbation-based approaches and dynamic network modeling for a comparative mutational profiling of the Hsp90 binding and allosteric interaction networks in the three Hsp90 maturation complexes with FKBP51 and P23 cochaperones and the glucocorticoid receptor (GR) client. The conformational dynamics signatures of the Hsp90 complexes and dynamics fluctuation analysis revealed how the intrinsic plasticity of the Hsp90 dimer can be modulated by cochaperones and client proteins to stabilize the closed dimer state required at the maturation stage of the ATPase cycle. In silico deep mutational scanning of the protein residues characterized the hot spots of protein stability and binding affinity in the Hsp90 complexes, showing that binding hot spots may often coincide with the regulatory centers that modulate dynamic allostery in the Hsp90 dimer. We introduce a perturbation-based network approach for mutational scanning of allosteric residue potentials and characterize allosteric switch clusters that control mechanism of cochaperone-dependent client recognition and remodeling by the Hsp90 chaperone. The results revealed a conserved network of allosteric switches in the Hsp90 complexes that allow cochaperones and GR protein to become integrated into the Hsp90 system by anchoring to the conformational switch points in the functional Hsp90 regions. This study suggests that the Hsp90 binding and allostery may operate under a regulatory mechanism in which activation or repression of the Hsp90 activity can be pre-encoded in the allosterically regulated Hsp90 dimer motions. By binding directly to the conformational switch centers on the Hsp90, cochaperones and interacting proteins can efficiently modulate the allosteric interactions and long-range communications required for client remodeling and activation.

Identification of Allosteric Effects in Proteins by Elastic Network Models

Allostery is a fundamental regulatory mechanism in the majority of biological processes of molecular machines. Allostery is well-known as a dynamic-driven process, and thus, the molecular mechanism of allosteric signal transmission needs to be established. Elastic network models (ENMs) provide efficient methods for investigating the intrinsic dynamics and allosteric communication pathways in proteins. In this chapter, two ENM methods including Gaussian network model (GNM) coupled with Markovian stochastic model, as well as the anisotropic network model (ANM), were introduced to identify allosteric effects in hemoglobins. Techniques on model parameters, scripting and calculation, analysis, and visualization are shown step by step.

Allosteric Regulation at the Crossroads of New Technologies: Multiscale Modeling, Networks, and Machine Learning

Allosteric regulation is a common mechanism employed by complex biomolecular systems for regulation of activity and adaptability in the cellular environment, serving as an effective molecular tool for cellular communication. As an intrinsic but elusive property, allostery is a ubiquitous phenomenon where binding or disturbing of a distal site in a protein can functionally control its activity and is considered as the "second secret of life." The fundamental biological importance and complexity of these processes require a multi-faceted platform of synergistically integrated approaches for prediction and characterization of allosteric functional states, atomistic reconstruction of allosteric regulatory mechanisms and discovery of allosteric modulators. The unifying theme and overarching goal of allosteric regulation studies in recent years have been integration between emerging experiment and computational approaches and technologies to advance quantitative characterization of allosteric mechanisms in proteins. Despite significant advances, the quantitative characterization and reliable prediction of functional allosteric states, interactions, and mechanisms continue to present highly challenging problems in the field. In this review, we discuss simulation-based multiscale approaches, experiment-informed Markovian models, and network modeling of allostery and information-theoretical approaches that can describe the thermodynamics and hierarchy allosteric states and the molecular basis of allosteric mechanisms. The wealth of structural and functional information along with diversity and complexity of allosteric mechanisms in therapeutically important protein families have provided a well-suited platform for development of data-driven research strategies. Data-centric integration of chemistry, biology and computer science using artificial intelligence technologies has gained a significant momentum and at the forefront of many cross-disciplinary efforts. We discuss new developments in the machine learning field and the emergence of deep learning and deep reinforcement learning applications in modeling of molecular mechanisms and allosteric proteins. The experiment-guided integrated approaches empowered by recent advances in multiscale modeling, network science, and machine learning can lead to more reliable prediction of allosteric regulatory mechanisms and discovery of allosteric modulators for therapeutically important protein targets.

Perturb-Scan-Pull: A Novel Method Facilitating Conformational Transitions in Proteins

Conformational transitions in proteins facilitate precise physiological functions. Therefore, it is crucial to understand the mechanisms underlying these processes to modulate protein function. Yet, studying structural and dynamical properties of proteins is notoriously challenging due to the complexity of the underlying potential energy surfaces (PES). We have previously developed the perturbation-response scanning (PRS) method to identify key residues that participate in the communication network responsible for specific conformational transitions. PRS is based on a residue-by-residue scan of the protein to determine the subset of residues/forces which provide the closest conformational change leading to a target conformational state, inasmuch as linear response theory applies to these motions. Here, we develop a novel method to further evaluate if conformational transitions may be triggered on the PES. We aim to study functionally relevant conformational transitions in proteins by using results obtained from PRS and feeding them as inputs to steered molecular dynamics simulations. The success and the transferability of the method are evaluated on three protein systems having different complexities of motion on the PES: calmodulin, adenylate kinase, and bacterial ferric binding protein. We find that the method captures the target conformation, while providing key residues and the optimum paths with relatively low free energy profiles.

Data-driven computational analysis of allosteric proteins by exploring protein dynamics, residue coevolution and residue interaction networks

BACKGROUND

Computational studies of allosteric interactions have witnessed a recent renaissance fueled by the growing interest in modeling of the complex molecular assemblies and biological networks. Allosteric interactions in protein structures allow for molecular communication in signal transduction networks.

METHODS

In this work, we performed a large scale comprehensive and multi-faceted analysis of >300 diverse allosteric proteins and complexes with allosteric modulators. By modeling and exploring coarse-grained dynamics, residue coevolution, and residue interaction networks for allosteric proteins, we have determined unifying molecular signatures shared by allosteric systems.

RESULTS

The results of this study have suggested that allosteric inhibitors and allosteric activators may differentially affect global dynamics and network organization of protein systems, leading to diverse allosteric mechanisms. By using structural and functional data on protein kinases, we present a detailed case study that that included atomic-level analysis of coevolutionary networks in kinases bound with allosteric inhibitors and activators.

CONCLUSIONS

We have found that coevolutionary networks can form direct communication pathways connecting functional regions and can recapitulate key regulatory sites and interactions responsible for allosteric signaling in the studied protein systems. The results of this computational investigation are compared with the experimental studies and reveal molecular signatures of known regulatory hotspots in protein kinases.

GENERAL SIGNIFICANCE

This study has shown that allosteric inhibitors and allosteric activators can have a different effect on residue interaction networks and can exploit distinct regulatory mechanisms, which could open up opportunities for probing allostery and new drug combinations with broad range of activities.

Integration of Random Forest Classifiers and Deep Convolutional Neural Networks for Classification and Biomolecular Modeling of Cancer Driver Mutations

Development of machine learning solutions for prediction of functional and clinical significance of cancer driver genes and mutations are paramount in modern biomedical research and have gained a significant momentum in a recent decade. In this work, we integrate different machine learning approaches, including tree based methods, random forest and gradient boosted tree (GBT) classifiers along with deep convolutional neural networks (CNN) for prediction of cancer driver mutations in the genomic datasets. The feasibility of CNN in using raw nucleotide sequences for classification of cancer driver mutations was initially explored by employing label encoding, one hot encoding, and embedding to preprocess the DNA information. These classifiers were benchmarked against their tree-based alternatives in order to evaluate the performance on a relative scale. We then integrated DNA-based scores generated by CNN with various categories of conservational, evolutionary and functional features into a generalized random forest classifier. The results of this study have demonstrated that CNN can learn high level features from genomic information that are complementary to the ensemble-based predictors often employed for classification of cancer mutations. By combining deep learning-generated score with only two main ensemble-based functional features, we can achieve a superior performance of various machine learning classifiers. Our findings have also suggested that synergy of nucleotide-based deep learning scores and integrated metrics derived from protein sequence conservation scores can allow for robust classification of cancer driver mutations with a limited number of highly informative features. Machine learning predictions are leveraged in molecular simulations, protein stability, and network-based analysis of cancer mutations in the protein kinase genes to obtain insights about molecular signatures of driver mutations and enhance the interpretability of cancer-specific classification models.

Atomistic Modeling of the ABL Kinase Regulation by Allosteric Modulators Using Structural Perturbation Analysis and Community-Based Network Reconstruction of Allosteric Communications

In this work, we have examined the molecular mechanisms of allosteric regulation of the ABL tyrosine kinase at the atomic level. Atomistic modeling of the ABL complexes with a panel of allosteric modulators has been performed using a combination of molecular dynamics simulations, structural residue perturbation scanning, and a novel community analysis of the residue interaction networks. Our results have indicated that allosteric inhibitors and activators may exert a differential control on allosteric signaling between the kinase binding sites and functional regions. While the inhibitor binding can strengthen the closed ABL state and induce allosteric communications directed from the allosteric pocket to the ATP binding site, the DPH activator may induce a more dynamic open form and activate allosteric couplings between the ATP and substrate binding sites. By leveraging a network-centric theoretical framework, we have introduced a novel community analysis method and global topological parameters that have unveiled the hierarchical modularity and the intercommunity bridging sites in the residue interaction network. We have found that allosteric functional hotspots responsible for the kinase regulation may serve the intermodular bridges in the global interaction network. The central conclusion from this analysis is that the regulatory switch centers play a fundamental role in the modular network organization of ABL as the unique intercommunity bridges that connect the SH2 and SH3 domains with the catalytic core into a functional kinase assembly. The hierarchy of network organization in the ABL regulatory complexes may allow for the synergistic action of dense intercommunity links required for the robust signal transfer in the catalytic core and sparse network bridges acting as the regulatory control points that orchestrate allosteric transitions between the inhibited and active kinase forms.

The Role of Heat Shock Proteins in Regulating Receptor Signal Transduction

Heat shock proteins (Hsp) are a class of stress-inducible proteins that mainly act as molecular protein chaperones. This chaperone activity is diverse, including assisting in nascent protein folding and regulating client protein location and translocation within the cell. The main proteins within the Hsp family, particularly Hsp70 and Hsp90, also have a highly diverse and numerous set of protein clients, which when combined with the high expression levels of Hsp proteins (2%-6% of total protein content) establishes these molecules as "central regulators" of cell protein physiology. Among the client proteins, Hsps regulate numerous signal-transduction and receptor-regulatory kinases, and indeed directly regulate some receptors themselves. This also makes the Hsps, particularly Hsp90, central regulators of signal-transduction machinery, with important impacts on endogenous and drug ligand responses. Among these roles, Hsp90 in particular acts to maintain mature signaling kinases in a metastable conformation permissive for signaling activation. In this review, we will focus on the roles of the Hsps, with a special focus on Hsp90, in regulating receptor signaling and subsequent physiologic responses. We will also explore potential means to manipulate Hsp function to improve receptor-targeted therapies. Overall, Hsps are important regulators of receptor signaling that are receiving increasing interest and exploration, particularly as Hsp90 inhibitors progress toward clinical approval for the treatment of cancer. Understanding the complex interplay of Hsp regulation of receptor signaling may provide important avenues to improve patient treatment.

Dissecting Structure-Encoded Determinants of Allosteric Cross-Talk between Post-Translational Modification Sites in the Hsp90 Chaperones

Post-translational modifications (PTMs) represent an important regulatory instrument that modulates structure, dynamics and function of proteins. The large number of PTM sites in the Hsp90 proteins that are scattered throughout different domains indicated that synchronization of multiple PTMs through a combinatorial code can be invoked as an important mechanism to orchestrate diverse chaperone functions and recognize multiple client proteins. In this study, we have combined structural and coevolutionary analysis with molecular simulations and perturbation response scanning analysis of the Hsp90 structures to characterize functional role of PTM sites in allosteric regulation. The results reveal a small group of conserved PTMs that act as global mediators of collective dynamics and allosteric communications in the Hsp90 structures, while the majority of flexible PTM sites serve as sensors and carriers of the allosteric structural changes. This study provides a comprehensive structural, dynamic and network analysis of PTM sites across Hsp90 proteins, identifying specific role of regulatory PTM hotspots in the allosteric mechanism of the Hsp90 cycle. We argue that plasticity of a combinatorial PTM code in the Hsp90 may be enacted through allosteric coupling between effector and sensor PTM residues, which would allow for timely response to structural requirements of multiple modified enzymes.