Exploring the mechanism of the PTP1B inhibitors by molecular dynamics and experimental study

Protein tyrosine phosphatase 1B (PTP1B) has proven to be an attractive target for the treatment of cancer, diabetes and other diseases. Although many PTP1B inhibitors with various scaffolds have been developed, there is still a lack of PTP1B inhibitor with high specificity and acceptable pharmacological properties. Therefore, it is urgent to develop more methods to explore complex action mode of PTP1B and ligands for designing ideal PTP1B modulators. In this work, we developed a potential molecular dynamics (MD) analytic mode to analyze the mechanism of active compounds 6a and 6e against PTP1B from different perspectives, including the stable ability, interactions and binding site of ligand and protein, the binding energy, relative movement between residues and changes in protein internal interactions. The simulated results demonstrated that compound 6a bound more stably to the active pocket of PTP1B than 6e due to its smaller molecular volume (326 Å3), matched electronegativity, and enhanced the positive correlation motion of residues, especially for WPD loop and P loop. Lastly, compound 6a as a competitive inhibitor for PTP1B was verified by enzyme kinetic assay. This work successfully studied the mechanism of compound 6a against PTP1B from various aspects, enriched the analysis of interaction mode between PTP1B and inhibitors. In summary, we hope that this work could provide more theoretical information for designing and developing more novel and ideal PTP1B inhibitors in the future.

Challenges and Opportunities in the Crusade of BRAF Inhibitors: From 2002 to 2022

Serine/threonine-protein kinase B-Raf (BRAF; RAF = rapidly accelerated fibrosarcoma) plays an important role in the mitogen-activated protein kinase (MAPK) signaling cascade. Somatic mutations in the BRAF gene were first discovered in 2002 by Davies et al., which was a major breakthrough in cancer research. Subsequently, three different classes of BRAF mutants have been discovered. This class includes class I monomeric mutants (BRAFV600), class II BRAF homodimer mutants (non-V600), and class III BRAF heterodimers (non-V600). Cancers caused by these include melanoma, thyroid cancer, ovarian cancer, colorectal cancer, nonsmall cell lung cancer, and others. In this study, we have highlighted the major binding pockets in BRAF protein, their active and inactive conformations with inhibitors, and BRAF dimerization and its importance in paradoxical activation and BRAF mutation. We have discussed the first-, second-, and third-generation drugs approved by the Food and Drug Administration and drugs under clinical trials with all four different binding approaches with DFG-IN/OUT and αC-IN/OUT for BRAF protein. We have investigated particular aspects and difficulties with all three generations of inhibitors. Finally, this study has also covered recent developments in synthetic BRAF inhibitors (from their discovery in 2002 to 2022), their unique properties, and importance in inhibiting BRAF mutants.

Experimental and theoretical study of the physicochemical properties of the novel imidazole-based eutectic solvent

Novel solvents and their applications are experiencing an increasing interest by the scientific community. Imidazole has been utilized as a major component in many successful ionic liquids. However, very limited studies were reported for using it as a hydrogen bond acceptor in the synthesis of eutectic solvents. In this work, a novel eutectic solvent composed of Imidazole and Monoethanolamine (MEA) is synthesized at different molar ratios. The basic physicochemical properties such as melting point, density, viscosity, and refractive index were measured at different temperatures and modeled as a function of molar composition and temperature. FTIR and ¹H NMR analyses were conducted and, the nature and strength of the molecular interaction between the two solvent molecules were investigated by conducting combined molecular dynamics (MD) simulations and density functional theory (DFT) calculations. The study revealed the electrostatic H-bonding nature of interaction with strength related to their bond distances. The binding energy between the two DES ingredients is proportional to the amount of MEA in the DES due to increasing the H-bonding interactions between Imidazole and MEA molecules. These findings suggest that DES might be used in a variety of chemical and industrial applications.

Gaussian field-based 3D-QSAR and molecular simulation studies to design potent pyrimidine-sulfonamide hybrids as selective BRAFV600E inhibitors

The "RAS-RAF-MEK-ERK" pathway is an important signaling pathway in melanoma. BRAFV600E (70-90%) is the most common mutation in this pathway. BRAF inhibitors have four types of conformers: type I (αC-IN/DFG-IN), type II (αC-IN/DFG-OUT), type I1/2 (αC-OUT/DFG-IN), and type I/II (αC-OUT/DFG-OUT). First- and second-generation BRAF inhibitors show resistance to BRAFV600E and are ineffective against malignancies induced by dimer BRAF mutants causing 'paradoxical' activation. In the present study, we performed molecular modeling of pyrimidine-sulfonamide hybrids inhibitors using 3D-QSAR, molecular docking, and molecular dynamics simulations. Previous reports reveal the importance of pyrimidine and sulfonamide moieties in the development of BRAFV600E inhibitors. Analysis of 3D-QSAR models provided novel pyrimidine sulfonamide hybrid BRAFV600E inhibitors. The designed compounds share similarities with several structural moieties present in first- and second-generation BRAF inhibitors. A total library of 88 designed compounds was generated and molecular docking studies were performed with them. Four molecules (T109, T183, T160, and T126) were identified as hits and selected for detailed studies. Molecular dynamics simulations were performed at 900 ns and binding was calculated. Based on molecular docking and simulation studies, it was found that the designed compounds have better interactions with the core active site [the nucleotide (ADP or ATP) binding site, DFG motif, and the phospho-acceptor site (activation segment) of BRAFV600E protein than previous inhibitors. Similar to the FDA-approved BRAFV600E inhibitors the developed compounds have [αC-OUT/DFG-IN] conformation. Compounds T126, T160 and T183 interacted with DIF (Leu505), making them potentially useful against BRAFV600E resistance and malignancies induced by dimer BRAF mutants. The synthesis and biological evaluation of the designed molecules is in progress, which may lead to some potent BRAFV600E selective inhibitors.

Ras isoform-specific expression, chromatin accessibility, and signaling

The anchorage of Ras isoforms in the membrane and their nanocluster formations have been studied extensively, including their detailed interactions, sizes, preferred membrane environments, chemistry, and geometry. However, the staggering challenge of their epigenetics and chromatin accessibility in distinct cell states and types, which we propose is a major factor determining their specific expression, still awaits unraveling. Ras isoforms are distinguished by their C-terminal hypervariable region (HVR) which acts in intracellular transport, regulation, and membrane anchorage. Here, we review some isoform-specific activities at the plasma membrane from a structural dynamic standpoint. Inspired by physics and chemistry, we recognize that understanding functional specificity requires insight into how biomolecules can organize themselves in different cellular environments. Within this framework, we suggest that isoform-specific expression may largely be controlled by the chromatin density and physical compaction, which allow (or curb) access to "chromatinized DNA." Genes are preferentially expressed in tissues: proteins expressed in pancreatic cells may not be equally expressed in lung cells. It is the rule-not an exception, and it can be at least partly understood in terms of chromatin organization and accessibility state. Genes are expressed when they can be sufficiently exposed to the transcription machinery, and they are less so when they are persistently buried in dense chromatin. Notably, chromatin accessibility can similarly determine expression of drug resistance genes.

Probing the Relation Between Community Evolution in Dynamic Residue Interaction Networks and Xylanase Thermostability

Residue-residue interactions are the basis of protein thermostability. The molecular conformations of Streptomyces lividans xylanase (xyna_strli) and Thermoascus aurantiacus xylanase (xyna_theau) at 300K, 325K and 350K were obtained by Molecular Dynamics (MD) simulations. Dynamic weighted residue interaction networks were constructed and the rigid-communities were detected using the ESPRA algorithm and the Evolving Graph+Fast-Newman algorithm. The residues in the rigid-communities are primarily located in loop2, short helixes α2^', α3^', α4^' and helixes α3 and α4. Thus, the rigid-community is close to the N-terminus of xylanase, which is usually stabilized to increase thermostability using site-directed mutagenesis. The evolution of the rigid-community with increasing temperature shows a stable synergistic interaction between loop2, α2^', α3^' and α4^' in xyna_theau. In particular, the short helixes α2^' and α3^' form a "thermo helix" to promote thermostability. In addition, tight global interactions between loop2, α2^', α3^', α3, α4^' and α4 of xyna_theau are identified, consisting mainly of hydrogen bonds, van der Waals forces and π-π stacking. These residue interactions are more resistant to high temperatures than those in xyna_strli. Robust residue interactions within these secondary structures are key factors influencing xyna_strli and xyna_theau thermostability. Analyzing the rigid-community can elucidate the cooperation of secondary structures, which cannot be discovered from sequence and 3D structure alone.

Exploring the dynamic mechanism of allosteric drug SHP099 inhibiting SHP2E69K

The E69K mutation is one of the most frequent protein tyrosine phosphatase-2 (SHP2) mutations in leukemia, and it can cause the increase in the protein activity. Recent studies have shown that the E69K mutation was fairly sensitive to the allosteric inhibitor of SHP2 (SHP099). However, the molecular mechanism of the allosteric drug SHP099 inhibiting SHP2^E69K remains unclear. Thus, the molecular dynamic simulations and the post-dynamics analyses (RMSF, PCA, DCCM, RIN and the binding free energies) for SHP2^WT, SHP2^WT-SHP099, SHP2^E69K and SHP2^E69K-SHP099 were carried out, respectively. Owing to the strong binding affinity of SHP099 to residues Thr219 and Arg220, the flexibility of linker region (residues Val209-Arg231) was reduced. Moreover, the presence of SHP099 kept the autoinhibition state of the SHP2 protein through enhancing the interactions between the linker region and Q loop in PTP domain, such as Thr219/Val490, Thr219/Asn491, Arg220/Ile488 and Leu254/Asn491. In addition, it was found that the residues (Thr219, Arg220, Leu254 and Asn491) might be the key residues responsible for the conformational changes of protein. Overall, this study may provide an important basis for understanding how the SHP099 effectively inhibited the SHP2^E69K activity at the molecular level.

Drug Research Meets Network Science: Where Are We?

Network theory provides one of the most potent analysis tools for the study of complex systems. In this paper, we illustrate the network-based perspective in drug research and how it is coherent with the new paradigm of drug discovery. We first present data sources from which networks are built, then show some examples of how the networks can be used to investigate drug-related systems. A section is devoted to network-based inference applications, i.e., prediction methods based on interactomes, that can be used to identify putative drug-target interactions without resorting to 3D modeling. Finally, we present some aspects of Boolean networks dynamics, anticipating that it might become a very potent modeling framework to develop in silico screening protocols able to simulate phenotypic screening experiments. We conclude that network applications integrated with machine learning and 3D modeling methods will become an indispensable tool for computational drug discovery in the next years.

The Mystery of Rap1 Suppression of Oncogenic Ras

Decades ago, Rap1, a small GTPase very similar to Ras, was observed to suppress oncogenic Ras phenotype, reverting its transformation. The proposed reason, persisting since, has been competition between Ras and Rap1 for a common target. Yet, none was found. There was also Rap1's puzzling suppression of Raf-1 versus activation of BRAF. Reemerging interest in Rap1 envisages capturing its Ras suppression action by inhibitors. Here, we review the literature and resolve the enigma. In vivo oncogenic Ras exists in isoform-distinct nanoclusters. The presence of Rap1 within the nanoclusters reduces the number of the clustered oncogenic Ras molecules, thus suppressing Raf-1 activation and mitogen-activated protein kinase (MAPK) signaling. Nanoclustering suggests that Rap1 suppression is Ras isoform dependent. Altogether, a potent Rap1-like inhibitor appears unlikely.

Soft regions of protein surface are potent for stable dimer formation

By having knowledge about the characteristics of protein interaction interfaces, we will be able to manipulate protein complexes for therapies. Dimer state is considered as the primary alphabet of the most proteins' quaternary structure. The properties of binding interface between subunits and of noninterface region define the specificity and stability of the intended protein complex. Considering some topological properties and amino acids' affinity for binding in interfaces of protein dimers, we construct the interface-specific recurrence plots. The data obtained from recurrence quantitative analysis, and accessibility-related metrics help us to classify the protein dimers into four distinct classes. Some mechanical properties of binding interfaces are computed for each predefined class of the dimers. The computed mechanical characteristics of binding patch region are compared with those of nonbinding region of proteins. Our observations indicate that the mechanical properties of protein binding sites have a decisive impact on determining the dimer stability. We introduce a new concept in analyzing protein structure by considering mechanical properties of protein structure. We conclude that the interface region between subunits of stable dimers is usually mechanically softer than the interface of unstable protein dimers. AbbreviationsAABaverage affinity for bindingANManisotropic network modelAPCaffinity propagation clusteringASAaccessible surface areaCCDinter residues distanceCSCcomplex stability codeDMdistance matrixΔG_dissPISA-computed dissociation free energyGNMGaussian normal mode analysisNMAnormal mode analysisPBPprotein binding patchPISAproteins, interfaces, structures and assembliesrASArelative accessible area in respect to unfolded state of residuesRMrecurrence matrixrPrelative protrusionRPrecurrence plotRQArecurrence quantitative analysisSEMstandard error of meanCommunicated by Ramaswamy H. Sarma.

A patent review of BRAF inhibitors: 2013-2018

Introduction: As a key element in arguably the most important pathway MAPK signaling, the BRAF kinase gives rise to severe diseases including cancers when pathologically activated. Extensive research on BRAFi (BRAF inhibitor) has been carried out to profile the characters for optimized agents and to elaborate the therapeutic strategies for the related cancer treatment. Areas covered: This review gives an overview of recently approved BRAF agents on function mode, therapeutic efficacy, and deficiency, based on which current challenges and corresponding strategies were presented. New entities as BRAFi for medical purpose in patent literature during the period 2013-2018 were also briefly introduced. Expert opinion: With the disclosure of paradox-breaker BRAFi PLX7904 crystal in complex with BRAF, the rational design for next-generation BRAFi is becoming ever more feasible. Accompanying therapeutic strategies in BRAFi elaboration may also provide flexible choice in the future 'personal medicine'. Further digging in the greatly enriched BRAFi pool will greatly benefit the drug design processes such as FBDD- and SBDD-driven development.

The structural basis for Ras activation of PI3Kα lipid kinase

PI3Kα is a principal Ras effector that phosphorylates PIP2 to PIP3 in the PI3K/Akt/mTOR pathway. How Ras activates PI3K has been unclear: is Ras' role confined to PI3K recruitment to the membrane or does Ras activation also involve allostery? Recently, we determined the mechanism of PI3Kα activation at the atomic level. We showed the vital role and significance of conformational change in PI3Kα activation. Here, by a 'best-match for hydrogen-bonding pair' (BMHP) computational protocol and molecular dynamics (MD) simulations, we model the atomic structure of KRas4B in complex with the Ras binding domain (RBD) of PI3Kα, striving to understand the mechanism of PI3Kα activation by Ras. Point mutations T208D, K210E, and K227E disrupt the KRas4B-RBD interface in the models, in line with the experiments. We identify allosteric signaling pathways connecting Ras to RBD in the p110α subunit. However, the observed weak allosteric signals coupled with the detailed mechanism of PI3Kα activation make us conclude that the dominant mechanistic role of Ras is likely to be recruitment and restriction of the PI3Kα population at the membrane. Thus, RTK recruits the PI3Kα to the membrane and activates it by relieving its autoinhibition exerted by the nSH2 domain, leading to exposure of the kinase domain, which permits PIP2 binding. Ras recruitment can shift the PI3Kα ensemble toward a population where the kinase domain surface and the active site position and orientation favor PIP2 insertion. This work helps elucidate Ras-mediated PI3K activation and explores the structural basis for Ras-PI3Kα drug discovery.

Integration of network models and evolutionary analysis into high-throughput modeling of protein dynamics and allosteric regulation: theory, tools and applications

Proteins are dynamical entities that undergo a plethora of conformational changes, accomplishing their biological functions. Molecular dynamics simulation and normal mode analysis methods have become the gold standard for studying protein dynamics, analyzing molecular mechanism and allosteric regulation of biological systems. The enormous amount of the ensemble-based experimental and computational data on protein structure and dynamics has presented a major challenge for the high-throughput modeling of protein regulation and molecular mechanisms. In parallel, bioinformatics and systems biology approaches including genomic analysis, coevolution and network-based modeling have provided an array of powerful tools that complemented and enriched biophysical insights by enabling high-throughput analysis of biological data and dissection of global molecular signatures underlying mechanisms of protein function and interactions in the cellular environment. These developments have provided a powerful interdisciplinary framework for quantifying the relationships between protein dynamics and allosteric regulation, allowing for high-throughput modeling and engineering of molecular mechanisms. Here, we review fundamental advances in protein dynamics, network theory and coevolutionary analysis that have provided foundation for rapidly growing computational tools for modeling of allosteric regulation. We discuss recent developments in these interdisciplinary areas bridging computational biophysics and network biology, focusing on promising applications in allosteric regulations, including the investigation of allosteric communication pathways, protein–DNA/RNA interactions and disease mutations in genomic medicine. We conclude by formulating and discussing future directions and potential challenges facing quantitative computational investigations of allosteric regulatory mechanisms in protein systems.

Covalent vs. Non-Covalent Inhibition: Tackling Drug Resistance in EGFR - A Thorough Dynamic Perspective

A persistent challenge in the treatment of non-small cell lung cancer (NSCLC) with EGFR is the emergence of drug-resistant caused by somatic mutations. The EGFR L858R/T790 M double mutant (EGFR^DM ) was found to be the most alarming variant. Despite the development of a wide range of inhibitors, none of them could inhibit EGFR^DM effectively. Recently, 11h and 45a, have been found to be potent inhibitors against EGFR^DM through two distinctive mechanisms, non-covalent and covalent binding, respectively. However, the structural and dynamic implications of the two modes of inhibitions remain unexplored. Herein, two molecular dynamics simulation protocols, coupled with free-energy calculations, were applied to gain insight into the atomistic nature of each binding mode. The comparative analysis confirmed that there is a significant difference in the binding free energy between 11h and 45a (ΔΔG_bind =-21.17 kcal/mol). The main binding force that governs the binding of both inhibitors is vdW, with a higher contribution for 45a. Two residues ARG841 and THR854 were found to have curtailed role in the binding of 45a to EGFR^DM by stabilizing its flexible alcohol chain. The 45a binding to EGFR^DM induces structural rearrangement in the active site to allow easier accessibility of 45a to target residue CYS797. The findings of this work can substantially shed light on new strategies for developing novel classes of covalent and non-covalent inhibitors with increased specificity and potency.

Identification and Biological Evaluation of Novel Type II B-RafV600E Inhibitors

The mitogen-activated protein kinase (MAPK) pathway plays a vital role in signal transduction networks. Severe diseases may be triggered if it is disturbed by mutated components, especially the kinase B-Raf^V600E . New inhibitors of the kinase are needed as cases of relapse and resistance with the known drugs have been widely reported in the clinic. In the present work, a new class of B-Raf^V600E inhibitors was identified by fragment linking. In vitro and in vivo assays were used to demonstrate the pharmacological properties of the compounds. 3-{3-[4-Chloro-3-(trifluoromethyl)phenyl]ureido}-N-[1-(4-methoxyphenyl)-1H-pyrazol-4-yl]benzamide was the most potent agent with IC₅₀ values of 0.035±0.004 μm (B-Raf^V600E kinase) and 0.39±0.04 μm (A375 cells). Furthermore, no obvious toxicity was observed. Collectively, the results favored justified the design rationale and hinted that this new chemotype might be worth studying further to develop novel B-Raf inhibitor candidates.

Using molecular simulation to explore the nanoscale dynamics of the plant kinome

Eukaryotic protein kinases (PKs) are a large family of proteins critical for cellular response to external signals, acting as molecular switches. PKs propagate biochemical signals by catalyzing phosphorylation of other proteins, including other PKs, which can undergo conformational changes upon phosphorylation and catalyze further phosphorylations. Although PKs have been studied thoroughly across the domains of life, the structures of these proteins are sparsely understood in numerous groups of organisms, including plants. In addition to efforts towards determining crystal structures of PKs, research on human PKs has incorporated molecular dynamics (MD) simulations to study the conformational dynamics underlying the switching of PK function. This approach of experimental structural biology coupled with computational biophysics has led to improved understanding of how PKs become catalytically active and why mutations cause pathological PK behavior, at spatial and temporal resolutions inaccessible to current experimental methods alone. In this review, we argue for the value of applying MD simulation to plant PKs. We review the basics of MD simulation methodology, the successes achieved through MD simulation in animal PKs, and current work on plant PKs using MD simulation. We conclude with a discussion of the future of MD simulations and plant PKs, arguing for the importance of molecular simulation in the future of plant PK research.

Development of the first model of a phosphorylated, ATP/Mg2+-containing B-Raf monomer by molecular dynamics simulations: a tool for structure-based design

A model of phosphorylated and ATP-containing B-Raf protein kinase is needed as a tool for the structure-based design of new allosteric inhibitors, since no crystal structure of such a system has been resolved. Here, we present the development of such a model as well as a thorough analysis of its structural features. This model was prepared using a systematic molecular dynamics approach considering the presence or absence of both the phosphate group at the Thr599 site and the ATP molecule. Then, different structural features (i.e. DFG motif, Mg²⁺ binding loop, activation loop, phosphorylation site and αC-helix region) were analysed for each trajectory to validate the aimed 2pBRAF_ATP model. Moreover, the structure and activating interactions of this 2pBRAF_ATP model were found to be in agreement with previously reported information. Finally, the model was further validated by means of a molecular docking study with our previously developed lead compound I confirming that this ATP-containing, phosphorylated protein model is suitable for further structure-based design studies.

The molecular mechanism behind resistance of the kinase FLT3 to the inhibitor quizartinib

Fms-like tyrosine kinase 3 (FLT3) is a receptor tyrosine kinase that is a drug target for leukemias. Several potent inhibitors of FLT3 exist, and bind to the inactive form of the enzyme. Unfortunately, resistance due to mutations in the kinase domain of FLT3 limits the therapeutic effects of these inhibitors. As in many other cases, it is not straightforward to explain why certain mutations lead to drug resistance. Extensive fully atomistic molecular dynamics (MD) simulations of FLT3 were carried out with an inhibited form (FLT-quizartinib complex), a free (apo) form, and an active conformation. In all cases, both the wild type (wt) proteins and two resistant mutants (D835F and Y842H) were studied. Analysis of the simulations revealed that impairment of protein-drug interactions cannot explain the resistance mutations in question. Rather, it appears that the active state of the mutant forms is perturbed by the mutations. It is therefore likely that perturbation of deactivation of the protein (which is necessary for drug binding) is responsible for the reduced affinity of the drug to the mutants. Importantly, this study suggests that it is possible to explain the source of resistance by mutations in FLT3 by an analysis of unbiased MD simulations.

Computational Analysis of Residue Interaction Networks and Coevolutionary Relationships in the Hsp70 Chaperones: A Community-Hopping Model of Allosteric Regulation and Communication

Allosteric interactions in the Hsp70 proteins are linked with their regulatory mechanisms and cellular functions. Despite significant progress in structural and functional characterization of the Hsp70 proteins fundamental questions concerning modularity of the allosteric interaction networks and hierarchy of signaling pathways in the Hsp70 chaperones remained largely unexplored and poorly understood. In this work, we proposed an integrated computational strategy that combined atomistic and coarse-grained simulations with coevolutionary analysis and network modeling of the residue interactions. A novel aspect of this work is the incorporation of dynamic residue correlations and coevolutionary residue dependencies in the construction of allosteric interaction networks and signaling pathways. We found that functional sites involved in allosteric regulation of Hsp70 may be characterized by structural stability, proximity to global hinge centers and local structural environment that is enriched by highly coevolving flexible residues. These specific characteristics may be necessary for regulation of allosteric structural transitions and could distinguish regulatory sites from nonfunctional conserved residues. The observed confluence of dynamics correlations and coevolutionary residue couplings with global networking features may determine modular organization of allosteric interactions and dictate localization of key mediating sites. Community analysis of the residue interaction networks revealed that concerted rearrangements of local interacting modules at the inter-domain interface may be responsible for global structural changes and a population shift in the DnaK chaperone. The inter-domain communities in the Hsp70 structures harbor the majority of regulatory residues involved in allosteric signaling, suggesting that these sites could be integral to the network organization and coordination of structural changes. Using a network-based formalism of allostery, we introduced a community-hopping model of allosteric communication. Atomistic reconstruction of signaling pathways in the DnaK structures captured a direction-specific mechanism and molecular details of signal transmission that are fully consistent with the mutagenesis experiments. The results of our study reconciled structural and functional experiments from a network-centric perspective by showing that global properties of the residue interaction networks and coevolutionary signatures may be linked with specificity and diversity of allosteric regulation mechanisms.

The diversity of allosteric mechanisms in the Hsp70 proteins could range from modulation of the inter-domain interactions and conformational dynamics to fine-tuning of the Hsp70 interactions with co-chaperones. The goal of this study is to present a systematic computational analysis of the dynamic and evolutionary factors underlying allosteric structural transformations of the Hsp70 proteins. We investigated the relationship between functional dynamics, residue coevolution, and network organization of residue interactions in the Hsp70 proteins. The results of this study revealed that conformational dynamics of the Hsp70 proteins may be linked with coevolutionary propensities and mutual information dependencies of the protein residues. Modularity and connectivity of allosteric interactions in the Hsp70 chaperones are coordinated by stable functional sites that feature unique coevolutionary signatures and high network centrality. The emergence of the inter-domain communities that are coordinated by functional centers and include highly coevolving residues could facilitate structural transitions through cooperative reorganization of the local interacting modules. We determined that the differences in the modularity of the residue interactions and organization of coevolutionary networks in DnaK may be associated with variations in their allosteric mechanisms. The network signatures of the DnaK structures are characteristic of a population-shift allostery that allows for coordinated structural rearrangements of local communities. A dislocation of mediating centers and insufficient coevolutionary coupling between functional regions may render a reduced cooperativity and promote a limited entropy-driven allostery in the Sse1 chaperone that occurs without structural changes. The results of this study showed that a network-centric framework and a community-hopping model of allosteric communication pathways may provide novel insights into molecular and evolutionary principles of allosteric regulation in the Hsp70 proteins.

Exploring Molecular Mechanisms of Paradoxical Activation in the BRAF Kinase Dimers: Atomistic Simulations of Conformational Dynamics and Modeling of Allosteric Communication Networks and Signaling Pathways

The recent studies have revealed that most BRAF inhibitors can paradoxically induce kinase activation by promoting dimerization and enzyme transactivation. Despite rapidly growing number of structural and functional studies about the BRAF dimer complexes, the molecular basis of paradoxical activation phenomenon is poorly understood and remains largely hypothetical. In this work, we have explored the relationships between inhibitor binding, protein dynamics and allosteric signaling in the BRAF dimers using a network-centric approach. Using this theoretical framework, we have combined molecular dynamics simulations with coevolutionary analysis and modeling of the residue interaction networks to determine molecular determinants of paradoxical activation. We have investigated functional effects produced by paradox inducer inhibitors PLX4720, Dabrafenib, Vemurafenib and a paradox breaker inhibitor PLX7904. Functional dynamics and binding free energy analyses of the BRAF dimer complexes have suggested that negative cooperativity effect and dimer-promoting potential of the inhibitors could be important drivers of paradoxical activation. We have introduced a protein structure network model in which coevolutionary residue dependencies and dynamic maps of residue correlations are integrated in the construction and analysis of the residue interaction networks. The results have shown that coevolutionary residues in the BRAF structures could assemble into independent structural modules and form a global interaction network that may promote dimerization. We have also found that BRAF inhibitors could modulate centrality and communication propensities of global mediating centers in the residue interaction networks. By simulating allosteric communication pathways in the BRAF structures, we have determined that paradox inducer and breaker inhibitors may activate specific signaling routes that correlate with the extent of paradoxical activation. While paradox inducer inhibitors may facilitate a rapid and efficient communication via an optimal single pathway, the paradox breaker may induce a broader ensemble of suboptimal and less efficient communication routes. The central finding of our study is that paradox breaker PLX7904 could mimic structural, dynamic and network features of the inactive BRAF-WT monomer that may be required for evading paradoxical activation. The results of this study rationalize the existing structure-functional experiments by offering a network-centric rationale of the paradoxical activation phenomenon. We argue that BRAF inhibitors that amplify dynamic features of the inactive BRAF-WT monomer and intervene with the allosteric interaction networks may serve as effective paradox breakers in cellular environment.