Unraveling the molecular mechanism of interactions of the Rho GTPases Cdc42 and Rac1 with the scaffolding protein IQGAP2

Disruption of Ca2+/calmodulin:KSR1 interaction lowers ERK activation

KSR1, a key scaffold protein for the MAPK pathway, facilitates ERK activation upon growth factor stimulation. We recently demonstrated that KSR1 binds the Ca2+-binding protein calmodulin (CaM), thereby providing an intersection between KSR1-mediated and Ca2+ signaling. In this study, we set out to generate a KSR1 point mutant with reduced Ca2+/CaM binding in order to unravel the functional implications of their interaction. To do so, we solved the structural determinants of complex formation. Using purified fragments of KSR1, we showed that Ca2+/CaM binds to the CA3 domain of KSR1. We then used in silico molecular modeling to predict contact residues for binding. This approach identified two possible modes of interaction: (1) binding of extended Ca2+/CaM to a globular conformation of KSR1-CA3 via electrostatic interactions or (2) binding of collapsed Ca2+/CaM to α-helical KSR1-CA3 via hydrophobic interactions. Experimentally, site-directed mutagenesis of the predicted contact residues for the two binding models favored that where collapsed Ca2+/CaM binds to the α-helical conformation of KSR1-CA3. Importantly, replacing KSR1-Phe355 with Asp reduces Ca2+/CaM binding by 76%. The KSR1-F355D mutation also significantly impairs the ability of EGF to activate ERK, which reveals that Ca2+/CaM binding promotes KSR1-mediated MAPK signaling. This work, by uncovering structural insight into the binding of KSR1 to Ca2+/CaM, identifies a KSR1 single-point mutant as a bioreagent to selectively study the crosstalk between Ca2+ and KSR1-mediated signaling.

Altered binding affinity of SIX1-Q177R correlates with enhanced WNT5A and WNT pathway effector expression in Wilms tumor

Wilms tumors present as an amalgam of varying proportions of tissues located within the developing kidney, one being the nephrogenic blastema comprising multipotent nephron progenitor cells (NPCs). The recurring missense mutation Q177R in NPC transcription factors SIX1 and SIX2 is most correlated with tumors of blastemal histology and is significantly associated with relapse. Yet, the transcriptional regulatory consequences of SIX1/2-Q177R that might promote tumor progression and recurrence have not been investigated extensively. Utilizing multiple Wilms tumor transcriptomic datasets, we identified upregulation of the gene encoding non-canonical WNT ligand WNT5A in addition to other WNT pathway effectors in SIX1/2-Q177R mutant tumors. SIX1 ChIP-seq datasets from Wilms tumors revealed shared binding sites for SIX1/SIX1-Q177R within a promoter of WNT5A and at putative distal cis-regulatory elements (CREs). We demonstrate colocalization of SIX1 and WNT5A in Wilms tumor tissue and utilize in vitro assays that support SIX1 and SIX1-Q177R activation of expression from the WNT5A CREs, as well as enhanced binding affinity within the WNT5A promoter that may promote the differential expression of WNT5A and other WNT pathway effectors associated with SIX1-Q177R tumors.

Computational approaches for the design of modulators targeting protein-protein interactions

BACKGROUND

Protein-protein interactions (PPIs) are intriguing targets for designing novel small-molecule inhibitors. The role of PPIs in various infectious and neurodegenerative disorders makes them potential therapeutic targets . Despite being portrayed as undruggable targets, due to their flat surfaces, disorderedness, and lack of grooves. Recent progresses in computational biology have led researchers to reconsider PPIs in drug discovery.

AREAS COVERED

In this review, we introduce in-silico methods used to identify PPI interfaces and present an in-depth overview of various computational methodologies that are successfully applied to annotate the PPIs. We also discuss several successful case studies that use computational tools to understand PPIs modulation and their key roles in various physiological processes.

EXPERT OPINION

Computational methods face challenges due to the inherent flexibility of proteins, which makes them expensive, and result in the use of rigid models. This problem becomes more significant in PPIs due to their flexible and flat interfaces. Computational methods like molecular dynamics (MD) simulation and machine learning can integrate the chemical structure data into biochemical and can be used for target identification and modulation. These computational methodologies have been crucial in understanding the structure of PPIs, designing PPI modulators, discovering new drug targets, and predicting treatment outcomes.

The Antithetic Roles of IQGAP2 and IQGAP3 in Cancers

The scaffold protein family of IQ motif-containing GTPase-activating proteins (IQGAP1, 2, and 3) share a high degree of homology and comprise six functional domains. IQGAPs bind and regulate the cytoskeleton, interact with MAP kinases and calmodulin, and have GTPase-related activity, as well as a RasGAP domain. Thus, IQGAPs regulate multiple cellular processes and pathways, affecting cell division, growth, cell-cell interactions, migration, and invasion. In the past decade, significant evidence on the function of IQGAPs in signal transduction during carcinogenesis has emerged. Compared with IQGAP1, IQGAP2 and IQGAP3 were less analyzed. In this review, we summarize the different signaling pathways affected by IQGAP2 and IQGAP3, and the antithetic roles of IQGAP2 and IQGAP3 in different types of cancer. IQGAP2 expression is reduced and plays a tumor suppressor role in most solid cancer types, while IQGAP3 is overexpressed and acts as an oncogene. In lymphoma, for example, IQGAPs have partially opposite functions. There is considerable evidence that IQGAPs regulate a multitude of pathways to modulate cancer processes and chemoresistance, but some questions, such as how they trigger this signaling, through which domains, and why they play opposite roles on the same pathways, are still unanswered.

Transcriptome analysis of PBMCs isolated from piglets treated with a miR-124 sponge construct identified miR124/IQGAP2/Rho GTPase as a target pathway support Salmonella Typhimurium infection

miR-124 is a significantly up-regulated miRNA in peripheral blood collected from piglets infected with Salmonella Typhimurium, suggesting that it may play an important role in Salmonella pathogenesis. This study focused on the transcriptomic analysis of peripheral blood mononuclear cells (PBMCs) isolated from miR-124 sponge and Salmonella Typhimurium-treated piglets, and trying to investigate the function of miR-124 in Salmonella infection. The transcriptome profiling analysis revealed that 2778 genes in miR-124 sponge + Salmonella Typhimurium treatment versus control, 2271 genes in Salmonella Typhimurium treatment versus control, and 1301 genes in miR-124 sponge + Salmonella Typhimurium versus Salmonella Typhimurium treatment, were differentially expressed, respectively (FDR < 0.05 and fold change > 2.0). Pathway analysis indicated that the MAPK signaling pathway, Ribosome pathway, and T-cell receptor signaling pathway were the most significantly enriched pathway in differentially expressed genes between miR-124 sponge + Salmonella Typhimurium and Salmonella Typhimurium along treatment (FDR < 0.05). Reporter assays and electrophoretic mobility shift assays showed that miR-124 is a crucial regulatory factor that targets IQ motif containing GTPase-activating protein 2 (IQGAP2). Cell culture experiment indicated that miR-124 attenuated the Salmonella Typhimurium-mediated activation of CDC42 and RAC1 (P < 0.05). Cultured PBMCs treated with miR-124 and IQGAP2-siRNA had higher intracellular Salmonella count than control samples, particularly 12 h post-infection (P < 0.05). Immunofluorescence analysis revealed that miR-124 treatment reduced the percentage of LAMP-1-positive phagosomes. The miR-124 could be an important regulator for IQGAP2/Rho GTPase pathway in Salmonella Typhimurium-infected PBMCs, and this pathway could be a target for Salmonella that support its infection in PBMCs in piglets.

Higher-order interactions of Bcr-Abl can broaden chronic myeloid leukemia (CML) drug repertoire

Bcr-Abl, a nonreceptor tyrosine kinase, is associated with leukemias, especially chronic myeloid leukemia (CML). Deletion of Abl's N-terminal region, to which myristoyl is linked, renders the Bcr-Abl fusion oncoprotein constitutively active. The substitution of Abl's N-terminal region by Bcr enables Bcr-Abl oligomerization. Oligomerization is critical: it promotes clustering on the membrane, which is essential for potent MAPK signaling and cell proliferation. Here we decipher the Bcr-Abl specific, step-by-step oligomerization process, identify a specific packing surface, determine exactly how the process is structured and identify its key elements. Bcr's coiled coil (CC) domain at the N-terminal controls Bcr-Abl oligomerization. Crystallography validated oligomerization via Bcr-Abl dimerization between two Bcr CC domains, with tetramerization via tight packing between two binary assemblies. However, the structural principles guiding Bcr CC domain oligomerization are unknown, hindering mechanistic understanding and drugs exploiting it. Using molecular dynamics (MD) simulations, we determine that the binary complex of the Bcr CC domain serves as a basic unit in the quaternary complex providing a specific surface for dimer-dimer packing and higher-order oligomerization. We discover that the small α1-helix is the key. In the binary assembly, the helix forms interchain aromatic dimeric packing, and in the quaternary assembly, it contributes to the specific dimer-dimer packing. Our mechanism is supported by the experimental literature. It offers the key elements controlling this process which can expand the drug discovery strategy, including by Bcr CC-derived peptides, and candidate residues for small covalent drugs, toward quenching oligomerization, supplementing competitive and allosteric tyrosine kinase inhibitors.

Discovery of small molecule inhibitors that effectively disrupt IQGAP1-Cdc42 interaction in breast cancer cells

The small GTPase Cdc42 is an integral component of the cytoskeleton, and its dysregulation leads to pathophysiological conditions, such as cancer. Binding of Cdc42 to the scaffold protein IQGAP1 stabilizes Cdc42 in its active form. The interaction between Cdc42 and IQGAP1 enhances migration and invasion of cancer cells. Disrupting this association could impair neoplastic progression and metastasis; however, no effective means to achieve this has been described. Here, we screened 78,500 compounds using a homogeneous time resolved fluorescence-based assay to identify small molecules that disrupt the binding of Cdc42 to IQGAP1. From the combined results of the validation assay and counter-screens, we selected 44 potent compounds for cell-based experiments. Immunoprecipitation and cell viability analysis rendered four lead compounds, namely NCGC00131308, NCGC00098561, MLS000332963 and NCGC00138812, three of which inhibited proliferation and migration of breast carcinoma cells. Microscale thermophoresis revealed that two compounds bind directly to Cdc42. One compound reduced the amount of active Cdc42 in cells and effectively impaired filopodia formation. Docking analysis provided plausible models of the compounds binding to the hydrophobic pocket adjacent to the GTP binding site of Cdc42. In conclusion, we identified small molecules that inhibit binding between Cdc42 and IQGAP1, which could potentially yield chemotherapeutic agents.

Tension of plus-end tracking protein Clip170 confers directionality and aggressiveness during breast cancer migration

The microtubule (MT) plus-end binding protein Clip170 is associated closely with breast cancer invasion and migration. In this study, Clip170 tension observed by a newly designed cpstFRET tension probe was suggested to be positive related to breast cancer aggressiveness, which could be regulated by α-tubulin detyrosination-induced MT disassembly. Clip170 phosphorylation induced by Ribosomal protein S6 kinase (RSK) could also increase its tension and promote the conversion of a discrete comet-like Clip-170 distribution into a spotty pattern during cancer metastasis. Heightened Clip170 tension was correlated with the formation of cortactin-associated filopodia and lamellipodia, and then promoted invasion and metastasis both in vitro and in vivo. Meanwhile, Clip170 tension enhanced at the leading edge in directional migration, accompanying with IQGAP1 subcellular distribution variation. Our work indicates that the malignancy and directionality during breast cancer migration depend on the magnitude and polarization of Clip170 tension, and we suggest Clip170 tension as a new potential drug target for breast cancer therapy.

CDC42-IQGAP Interactions Scrutinized: New Insights into the Binding Properties of the GAP-Related Domain

The IQ motif-containing GTPase-activating protein (IQGAP) family composes of three highly-related and evolutionarily conserved paralogs (IQGAP1, IQGAP2 and IQGAP3), which fine tune as scaffolding proteins numerous fundamental cellular processes. IQGAP1 is described as an effector of CDC42, although its effector function yet re-mains unclear. Biophysical, biochemical and molecular dynamic simulation studies have proposed that IQGAP RASGAP-related domains (GRDs) bind to the switch regions and the insert helix of CDC42 in a GTP-dependent manner. Our kinetic and equilibrium studies have shown that IQGAP1 GRD binds, in contrast to its C-terminal 794 amino acids (called C794), CDC42 in a nucleotide-independent manner indicating a binding outside the switch regions. To resolve this discrepancy and move beyond the one-sided view of GRD, we carried out affinity measurements and a systematic mutational analysis of the interfacing residues between GRD and CDC42 based on the crystal structure of the IQGAP2 GRD-CDC42Q61L GTP complex. We determined a 100-fold lower affinity of the GRD1 of IQGAP1 and of GRD2 of IQGAP2 for CDC42 mGppNHp in comparison to C794/C795 proteins. Moreover, partial and major mutation of CDC42 switch regions substantially affected C794/C795 binding but only a little GRD1 and remarkably not at all the GRD2 binding. However, we clearly showed that GRD2 contributes to the overall affinity of C795 by using a 11 amino acid mutated GRD variant. Furthermore, the GRD1 binding to the CDC42 was abolished using specific point mutations within the insert helix of CDC42 clearly supporting the notion that CDC42 binding site(s) of IQGAP GRD lies outside the switch regions among others in the insert helix. Collectively, this study provides further evidence for a mechanistic framework model that is based on a multi-step binding process, in which IQGAP GRD might act as a ‘scaffolding domain’ by binding CDC42 irrespective of its nucleotide-bound forms, followed by other IQGAP domains downstream of GRD that act as an effector domain and is in charge for a GTP-dependent interaction with CDC42.

The structural basis of BCR-ABL recruitment of GRB2 in chronic myelogenous leukemia

BCR-ABL drives chronic myeloid leukemia (CML). BCR binding to GRB2 transduces signaling via the Ras/MAPK pathway. Despite considerable data confirming the binding, molecular-level understanding of exactly how the two proteins interact, and, especially, what are the determinants of the specificity of the SH2GRB2 domain-phosphorylated BCR (pBCR) recognition are still open questions. Yet, this is vastly important for understanding binding selectivity, and for predicting the phosphorylated receptors, or peptides, that are likely to bind. Here, we uncover these determinants and ascertain to what extent they relate to the affinity of the interaction. Toward this end, we modeled the complexes of the pBCR and SH2GRB2 and other pY/Y-peptide-SH2 complexes and compared their specificity and affinity. We observed that pBCR's 176FpYVNV180 motif is favorable and specific to SH2GRB2, similar to pEGFR, but not other complexes. SH2GRB2 contains two binding pockets: pY-binding recognition pocket triggers binding, and the specificity pocket whose interaction is governed by N179 in pBCR and W121 in SH2GRB2. Our proposed motif with optimal affinity to SH2GRB2 is E/D-pY-E/V-N-I/L. Collectively, we provide the structural basis of BCR-ABL recruitment of GRB2, outline its specificity hallmarks, and delineate a blueprint for prediction of BCR-binding scaffolds and for therapeutic peptide design.

Role of IQ Motif-Containing GTPase-Activating Proteins in Hepatocellular Carcinoma

IQ motif-containing GTPase-activating proteins (IQGAPs) are a class of scaffolding proteins, including IQGAP1, IQGAP2, and IQGAP3, which govern multiple cellular activities by facilitating cytoskeletal remodeling and cellular signal transduction. The role of IQGAPs in cancer initiation and progression has received increasing attention in recent years, especially in hepatocellular carcinoma (HCC), where the aberrant expression of IQGAPs is closely related to patient prognosis. IQGAP1 and 3 are upregulated and are considered oncogenes in HCC, while IQGAP2 is downregulated and functions as a tumor suppressor. This review details the three IQGAP isoforms and their respective structures. The expression and role of each protein in different liver diseases and mainly in HCC, as well as the underlying mechanisms, are also presented. This review also provides a reference for further studies on IQGAPs in HCC.

The mechanism of activation of MEK1 by B-Raf and KSR1

MEK1 interactions with B-Raf and KSR1 are key steps in Ras/Raf/MEK/ERK signaling. Despite this, vital mechanistic details of how these execute signal transduction are still enigmatic. Among these is why, despite B-Raf and KSR1 kinase domains similarity, the B-Raf/MEK1 and KSR1/MEK1 complexes have distinct contributions to MEK1 activation, and broadly, what is KSR1's role. Our molecular dynamics simulations clarify these still unresolved ambiguities. Our results reveal that the proline-rich (P-rich) loop of MEK1 plays a decisive role in MEK1 activation loop (A-loop) phosphorylation. In the inactive B-Raf/MEK1 heterodimer, the collapsed A-loop of B-Raf interacts with the P-rich loop and A-loop of MEK1, minimizing MEK1 A-loop fluctuation and preventing it from phosphorylation. In the active B-Raf/MEK1 heterodimer, the P-rich loop moves in concert with the A-loop of B-Raf as it extends. This reduces the number of residues interacting with MEK1 A-loop, allowing increased A-loop fluctuation, and bringing Ser222 closer to ATP for phosphorylation. B-Raf αG-helix Arg662 promotes MEK1 activation by orienting Ser218 towards ATP. In KSR1/MEK1, the KSR1 αG-helix has Ala826 in place of B-Raf Arg662. This difference results in much fewer interactions between KSR1 αG-helix and MEK1 A-loop, thus a more flexible A-loop. We postulate that if KSR1 were to adopt an active configuration with an extended A-loop as seen in other protein kinases, then the MEK1 P-rich loop would extend in a similar manner, as seen in the active B-Raf/MEK1 heterodimer. This would result in highly flexible MEK1 A-loop, and KSR1 functioning as an active, B-Raf-like, kinase.

Structural basis of activation of the tumor suppressor protein neurofibromin

Mutations in the NF1 gene cause the familial genetic disease neurofibromatosis type I, as well as predisposition to cancer. The NF1 gene product, neurofibromin, is a GTPase-activating protein and acts as a tumor suppressor by negatively regulating the small GTPase, Ras. However, structural insights into neurofibromin activation remain incompletely defined. Here, we provide cryoelectron microscopy (cryo-EM) structures that reveal an extended neurofibromin homodimer in two functional states: an auto-inhibited state with occluded Ras-binding site and an asymmetric open state with an exposed Ras-binding site. Mechanistically, the transition to the active conformation is stimulated by nucleotide binding, which releases a lock that tethers the catalytic domain to an extended helical repeat scaffold in the occluded state. Structure-guided mutational analysis supports functional relevance of allosteric control. Disease-causing mutations are mapped and primarily impact neurofibromin stability. Our findings suggest a role for nucleotides in neurofibromin regulation and may lead to therapeutic modulation of Ras signaling.

Prediction of Protein–Protein Interactions Between Alsin DH/PH and Rac1 and Resulting Protein Dynamics

Alsin is a protein of 1,657 amino acids known for its crucial role in vesicular trafficking in neurons thanks to its ability to interact with two guanosine triphosphatases, Rac1 and Rab5. Evidence suggests that Rac1 can bind Alsin central region, composed by a Dbl Homology (DH) domain followed by a Pleckstrin Homology (PH) domain, leading to Alsin relocalization. However, Alsin three-dimensional structure and its relationship with known biological functions of this protein are still unknown. In this work, a homology model of the Alsin DH/PH domain was developed and studied through molecular dynamics both in the presence and in the absence of its binding partner, Rac1. Due to different conformations of DH domain, the presence of Rac1 seems to stabilize an open state of the protein, while the absence of its binding partner results in closed conformations. Furthermore, Rac1 interaction was able to reduce the fluctuations in the second conserved region of DH motif, which may be involved in the formation of a homodimer. Moreover, the dynamics of DH/PH was described through a Markov State Model to study the pathways linking the open and closed states. In conclusion, this work provided an all-atom model for the DH/PH domain of Alsin protein; moreover, molecular dynamics investigations suggested underlying molecular mechanisms in the signal transduction between Rac1 and Alsin, providing the basis for a deeper understanding of the whole structure–function relationship for Alsin protein.

Selectivity Determinants of RHO GTPase Binding to IQGAPs

IQ motif-containing GTPase-activating proteins (IQGAPs) modulate a wide range of cellular processes by acting as scaffolds and driving protein components into distinct signaling networks. Their functional states have been proposed to be controlled by members of the RHO family of GTPases, among other regulators. In this study, we show that IQGAP1 and IQGAP2 can associate with CDC42 and RAC1-like proteins but not with RIF, RHOD, or RHO-like proteins, including RHOA. This seems to be based on the distribution of charged surface residues, which varies significantly among RHO GTPases despite their high sequence homology. Although effector proteins bind first to the highly flexible switch regions of RHO GTPases, additional contacts outside are required for effector activation. Sequence alignment and structural, mutational, and competitive biochemical analyses revealed that RHO GTPases possess paralog-specific residues outside the two highly conserved switch regions that essentially determine the selectivity of RHO GTPase binding to IQGAPs. Amino acid substitution of these specific residues in RHOA to the corresponding residues in RAC1 resulted in RHOA association with IQGAP1. Thus, electrostatics most likely plays a decisive role in these interactions.

The small GTPase KlRho5 responds to oxidative stress and affects cytokinesis

Rho5 is the yeast homolog of the human small GTPase Rac1. We characterized the genes encoding Rho5 and the subunits of its dimeric activating guanine-nucleotide-exchange factor (GEF), Dck1 and Lmo1, in the yeast Kluyveromyces lactis. Rapid translocation of the three GFP-tagged components to mitochondria upon oxidative stress and carbon starvation indicate a similar function of KlRho5 in energy metabolism and mitochondrial dynamics as described for its Saccharomyces cerevisiae homolog. Accordingly, Klrho5 deletion mutants are hyper-resistant towards hydrogen peroxide. Moreover, synthetic lethalities of rho5 deletions with key components in nutrient sensing, such as sch9 and gpr1, are not conserved in K. lactis. Instead, Klrho5 deletion mutants display morphological defects with strengthened lateral cell walls and protruding bud scars. The latter result from aberrant cytokinesis, as observed by following the budding process in vivo and by transmission electron microscopy of the bud neck region. This phenotype can be suppressed by KlCDC42G12V, which encodes a hyper-active variant. Data from live-cell fluorescence microscopy support the notion that KlRho5 interferes with the actin moiety of the contractile actomyosin ring, with consequences different from those previously reported for mutants lacking myosin.

The mechanism of activation of monomeric B-Raf V600E

Oncogenic mutations in the serine/threonine kinase B-Raf, particularly the V600E mutation, are frequent in cancer, making it a major drug target. Although much is known about B-Raf's active and inactive states, questions remain about the mechanism by which the protein changes between these two states. Here, we utilize molecular dynamics to investigate both wild-type and V600E B-Raf to gain mechanistic insights into the impact of the Val to Glu mutation. The results show that the wild-type and mutant follow similar activation pathways involving an extension of the activation loop and an inward motion of the αC-helix. The V600E mutation, however, destabilizes the inactive state by disrupting hydrophobic interactions present in the wild-type structure while the active state is stabilized through the formation of a salt bridge between Glu600 and Lys507. Additionally, when the activation loop is extended, the αC-helix is able to move between an inward and outward orientation as long as the DFG motif adopts a specific orientation. In that orientation Phe595 rotates away from the αC-helix, allowing the formation of a salt bridge between Lys483 and Glu501. These mechanistic insights have implications for the development of new Raf inhibitors.

The structural basis of Akt PH domain interaction with calmodulin

Akt plays a key role in the Ras/PI3K/Akt/mTOR signaling pathway. In breast cancer, Akt translocation to the plasma membrane is enabled by the interaction of its pleckstrin homology domain (PHD) with calmodulin (CaM). At the membrane, the conformational change promoted by PIP3 releases CaM and facilitates Thr308 and Ser473 phosphorylation and activation. Here, using modeling and molecular dynamics simulations, we aim to figure out how CaM interacts with Akt's PHD at the atomic level. Our simulations show that CaM-PHD interaction is thermodynamically stable and involves a β-strand rather than an α-helix, in agreement with NMR data, and that electrostatic and hydrophobic interactions are critical. The PHD interacts with CaM lobes; however, multiple modes are possible. IP4, the polar head of PIP3, weakens the CaM-PHD interaction, implicating the release mechanism at the plasma membrane. Recently, we unraveled the mechanism of PI3Kα activation at the atomistic level and the structural basis for Ras role in the activation. Here, our atomistic structural data clarify the mechanism of how CaM interacts, delivers, and releases Akt-the next node in the Ras/PI3K pathway-at the plasma membrane.

Mechanistic Differences of Activation of Rac1P29S and Rac1A159V

Rac1 is a small GTPase that plays key roles in actin reorganization, cell motility, and cell survival/growth as well as in various cancer types and neurodegenerative diseases. Similar to other Ras superfamily GTPases, Rac1 switches between active GTP-bound and inactive GDP-bound states. Switch I and II regions open and close during GDP/GTP exchange. P29S and A159V (paralogous to K-RasA146) mutations are the two most common somatic mutations of Rac1. Rac1P29S is a known hotspot for melanoma, whereas Rac1A159V most commonly occurs in head and neck cancer. To investigate how these substitutions induce the Rac1 dynamics, we used atomistic molecular dynamics simulations on the wild-type Rac1 and two mutant systems (P29S and A159V) in the GTP bound state, and on the wild-type Rac1 and P29S mutated system in the GDP bound state. Here, we show that P29S and A159V mutations activate Rac1 with different mechanisms. In Rac1P29S-GTP, the substitution increases the flexibility of Switch I based on RMSF and dihedral angle calculations and leads to an open conformation. We propose that the open Switch I conformation is one of the underlying reasons for rapid GDP/GTP exchange of Rac1P29S. On the other hand, in Rac1A159V-GTP, some of the contacts of the guanosine ring of GTP with Rac1 are temporarily lost, enabling the guanosine ring to move toward Switch I and subsequently close the switch. Rac1A159V-GTP adopts a Ras state 2 like conformation, where both switch regions are in closed conformation and Thr35 forms a hydrogen bond with the nucleotide.

Active and Inactive Cdc42 Differ in Their Insert Region Conformational Dynamics

Cell division control protein 42 homolog (Cdc42) protein, a Ras superfamily GTPase, regulates cellular activities, including cancer progression. Using all-atom molecular dynamics (MD) simulations and essential dynamic analysis, we investigated the structure and dynamics of the catalytic domains of GDP-bound (inactive) and GTP-bound (active) Cdc42 in solution. We discovered substantial differences in the dynamics of the inactive and active forms, particularly in the "insert region" (residues 122-135), which plays a role in Cdc42 activation and binding to effectors. The insert region has larger conformational flexibility in the GDP-bound Cdc42 than in the GTP-bound Cdc42. The G2 loop and switch I at the effector lobe of the catalytic domain exhibit large conformational changes in both the GDP- and the GTP-bound systems, but in the GTP-bound Cdc42, the switch I interactions with GTP are retained. Oncogenic mutations were identified in the Ras superfamily. In Cdc42, the G12V and Q61L mutations decrease the GTPase activity. We simulated these mutations in both GDP- and GTP-bound Cdc42. Although the overall structural organization is quite similar between the wild type and the mutants, there are small differences in the conformational dynamics, especially in the two switch regions. Taken together, the G12V and Q61L mutations may play a role similar to their K-Ras counterparts in nucleotide binding and activation. The conformational differences, which are mainly in the insert region and, to a lesser extent, in the switch regions flanking the nucleotide binding site, can shed light on binding and activation. We propose that the differences are due to a network of hydrogen bonds that gets disrupted when Cdc42 is bound to GDP, a disruption that does not exist in other Rho GTPases. The differences in the dynamics between the two Cdc42 states suggest that the inactive conformation has reduced ability to bind to effectors.

The Arg-293 of Cryptochrome1 is responsible for the allosteric regulation of CLOCK-CRY1 binding in circadian rhythm

Mammalian circadian clocks are driven by transcription/translation feedback loops composed of positive transcriptional activators (BMAL1 and CLOCK) and negative repressors (CRYPTOCHROMEs (CRYs) and PERIODs (PERs)). CRYs, in complex with PERs, bind to the BMAL1/CLOCK complex and repress E-box-driven transcription of clock-associated genes. There are two individual CRYs, with CRY1 exhibiting higher affinity to the BMAL1/CLOCK complex than CRY2. It is known that this differential binding is regulated by a dynamic serine-rich loop adjacent to the secondary pocket of both CRYs, but the underlying features controlling loop dynamics are not known. Here we report that allosteric regulation of the serine-rich loop is mediated by Arg-293 of CRY1, identified as a rare CRY1 SNP in the Ensembl and 1000 Genomes databases. The p.Arg293His CRY1 variant caused a shortened circadian period in a Cry1-/-Cry2-/- double knockout mouse embryonic fibroblast cell line. Moreover, the variant displayed reduced repressor activity on BMAL1/CLOCK driven transcription, which is explained by reduced affinity to BMAL1/CLOCK in the absence of PER2 compared with CRY1. Molecular dynamics simulations revealed that the p.Arg293His CRY1 variant altered a communication pathway between Arg-293 and the serine loop by reducing its dynamicity. Collectively, this study provides direct evidence that allosterism in CRY1 is critical for the regulation of circadian rhythm.

The quaternary assembly of KRas4B with Raf-1 at the membrane

Proximally located in the membrane, oncogenic Ras dimers (or nanoclusters) can recruit and promote Raf dimerization and MAPK (Raf/MEK/ERK) signaling. Among Ras isoforms, KRas4B is the most frequently mutated. Recent data on the binary KRas4B–Raf-1 complex suggested that Raf-1 CRD not only executes membrane anchorage, but also supports the high-affinity interaction of Raf-1 RBD with KRas4B catalytic domain. For a detailed mechanistic picture of Raf activation at the membrane, we employ explicit MD simulations of the quaternary KRas4B–Raf-1 complex. The complex contains two active GTP-bound KRas4B proteins forming a dimer through the allosteric lobe interface and two tandem RBD-CRD segments of Raf-1 interacting with the effector lobes at both ends of the KRas4B dimer. We show that Raf-1 RBD-CRD supports stable KRas4B dimer at preferred interface and orientation at the membrane, thereby cooperatively enhancing the affinity of the KRas4B–Raf-1 interaction. We propose that a Ras dimer at the membrane can increase the population of proximal Raf kinase domains, promoting kinase domain dimerization in the cytoplasm. Collectively, the dynamic Ras–Raf assembly promotes Raf activation not by allostery; instead, Ras activates Raf by shifting its ensemble toward kinase domain-accessible states through enhanced affinity at the membrane.

Ubiquitination of the scaffold protein IQGAP1 diminishes its interaction with and activation of the Rho GTPase CDC42

IQ motif-containing GTPase-activating protein 1 (IQGAP1) is a scaffold protein that interacts with numerous binding partners and thereby regulates fundamental biological processes. The functions of IQGAP1 are modulated by several mechanisms, including protein binding, self-association, subcellular localization, and phosphorylation. Proteome-wide screens have indicated that IQGAP1 is ubiquitinated, but the possible effects of this post-translational modification on its function are unknown. Here we characterized and evaluated the function of IQGAP1 ubiquitination. Using MS-based analysis in HEK293 cells, we identified six lysine residues (Lys-556, -1155, -1230, -1465, -1475, and -1528) as ubiquitination sites in IQGAP1. To elucidate the biological consequences of IQGAP1 ubiquitination, we converted each of these lysines to arginine and found that replacing two of these residues, Lys-1155 and Lys-1230, in the GAP-related domain of IQGAP1 (termed IQGAP1 GRD-2K) reduces its ubiquitination. Moreover, IQGAP1 GRD-2K bound a significantly greater proportion of the two Rho GTPases cell division cycle 42 (CDC42) and Rac family small GTPase 1 (RAC1) than did WT IQGAP1. Consistent with this observation, reconstitution of IQGAP1-null cells with IQGAP1 GRD-2K significantly increased the amount of active CDC42 and enhanced cell migration significantly more than WT IQGAP1. Our results reveal that ubiquitination of the CDC42 regulator IQGAP1 alters its ability to bind to and activate this GTPase, leading to physiological effects. Collectively, these findings expand our view of the role of ubiquitination in cell signaling and provide additional insight into CDC42 regulation.

Ca2+-Dependent Switch of Calmodulin Interaction Mode with Tandem IQ Motifs in the Scaffolding Protein IQGAP1

IQ domain GTPase-activating scaffolding protein 1 (IQGAP1) mediates cytoskeleton, cell migration, proliferation, and apoptosis events. Calmodulin (CaM) modulates IQGAP1 functions by binding to its four tandem IQ motifs. Exactly how CaM binds the IQ motifs and which functions of IQGAP1 CaM regulates and how are fundamental mechanistic questions. We combine experimental pull-down assays, mutational data, and molecular dynamics simulations to understand the IQ-CaM complexes with and without Ca2+ at the atomic level. Apo-CaM favors the IQ3 and IQ4 motifs but not the IQ1 and IQ2 motifs that lack two hydrophobic residues for interactions with apo-CaM's hydrophobic pocket. Ca2+-CaM binds all four IQ motifs, with both N- and C-lobes tightly wrapped around each motif. Ca2+ promotes IQ-CaM interactions and increases the amount of IQGAP1-loaded CaM for IQGAP1-mediated signaling. Collectively, we describe IQ-CaM binding in atomistic detail and feature the emergence of Ca2+ as a key modulator of the CaM-IQGAP1 interactions.

The Structural Basis of the Farnesylated and Methylated KRas4B Interaction with Calmodulin

Ca²⁺-calmodulin (CaM) extracts KRas4B from the plasma membrane, suggesting that KRas4B/CaM interaction plays a role in regulating Ras signaling. To gain mechanistic insight, we provide a computational model, supported by experimental structural data, of farnesylated/methylated KRas4B_1-185 interacting with CaM in solution and at anionic membranes including signaling lipids. Due to multiple interaction modes, we observe diverse conformational ensembles of the KRas4B-CaM complex. A highly populated conformation reveals the catalytic domain interacting with the N-lobe and the hypervariable region (HVR) wrapping around the linker with the farnesyl docking to the extended CaM's C-lobe pocket. Alternatively, KRas4B can interact with collapsed CaM with the farnesyl penetrating CaM's center. At anionic membranes, CaM interacts with the catalytic domain with large fluctuations, drawing the HVR. Signaling lipids establishing strong salt bridges with CaM prevent membrane departure. Membrane-interacting KRas4B-CaM complex can productively recruit phosphatidylinositol 3-kinase α (PI3Kα) to the plasma membrane, serving as a coagent in activating PI3Kα/Akt signaling.

Emerging roles of allosteric modulators in the regulation of protein-protein interactions (PPIs): A new paradigm for PPI drug discovery

Protein-protein interactions (PPIs) are closely implicated in various types of cellular activities and are thus pivotal to health and disease states. Given their fundamental roles in a wide range of biological processes, the modulation of PPIs has enormous potential in drug discovery. However, owing to the general properties of large, flat, and featureless interfaces of PPIs, previous attempts have demonstrated that the generation of therapeutic agents targeting PPI interfaces is challenging, rendering them almost "undruggable" for decades. To date, rapid progress in chemical and structural biology techniques has promoted the exploitation of allostery as a novel approach in drug discovery. By attaching to allosteric sites that are topologically and spatially distinct from PPI interfaces, allosteric modulators can achieve improved physiochemical properties. Thus, allosteric modulators may represent an alternative strategy to target intractable PPIs and have attracted intense pharmaceutical interest. In this review, we first briefly introduce the characteristics of PPIs and then present different approaches for investigating PPIs, as well as the latest methods for modulating PPIs. Importantly, we comprehensively review the recent progress in the development of allosteric modulators to inhibit or stabilize PPIs. Finally, we conclude with future perspectives on the discovery of allosteric PPI modulators, especially the application of computational methods to aid in allosteric PPI drug discovery.

Computational Structural Biology: Successes, Future Directions, and Challenges

Computational biology has made powerful advances. Among these, trends in human health have been uncovered through heterogeneous 'big data' integration, and disease-associated genes were identified and classified. Along a different front, the dynamic organization of chromatin is being elucidated to gain insight into the fundamental question of genome regulation. Powerful conformational sampling methods have also been developed to yield a detailed molecular view of cellular processes. when combining these methods with the advancements in the modeling of supramolecular assemblies, including those at the membrane, we are finally able to get a glimpse into how cells' actions are regulated. Perhaps most intriguingly, a major thrust is on to decipher the mystery of how the brain is coded. Here, we aim to provide a broad, yet concise, sketch of modern aspects of computational biology, with a special focus on computational structural biology. We attempt to forecast the areas that computational structural biology will embrace in the future and the challenges that it may face. We skirt details, highlight successes, note failures, and map directions.

Methods for Discovering and Targeting Druggable Protein-Protein Interfaces and Their Application to Repurposing

Drug repurposing is a creative and resourceful approach to increase the number of therapies by exploiting available and approved drugs. However, identifying new protein targets for previously approved drugs is challenging. Although new strategies have been developed for drug repurposing, there is broad agreement that there is room for further improvements. In this chapter, we review protein-protein interaction (PPI) interface-targeting strategies for drug repurposing applications. We discuss certain features, such as hot spot residue and hot region prediction and their importance in drug repurposing, and illustrate common methods used in PPI networks to identify drug off-targets. We also collect available online resources for hot spot prediction, binding pocket identification, and interface clustering which are effective resources in polypharmacology. Finally, we provide case studies showing the significance of protein interfaces and hot spots in drug repurposing.

Autoinhibition in Ras effectors Raf, PI3Kα, and RASSF5: a comprehensive review underscoring the challenges in pharmacological intervention

Autoinhibition is an effective mechanism that guards proteins against spurious activation. Despite its ubiquity, the distinct organizations of the autoinhibited states and their release mechanisms differ. Signaling is most responsive to the cell environment only if a small shift in the equilibrium is required to switch the system from an inactive (occluded) to an active (exposed) state. Ras signaling follows this paradigm. This underscores the challenge in pharmacological intervention to exploit and enhance autoinhibited states. Here, we review autoinhibition and release mechanisms at the membrane focusing on three representative Ras effectors, Raf protein kinase, PI3Kα lipid kinase, and NORE1A (RASSF5) tumor suppressor, and point to the ramifications to drug discovery. We further touch on Ras upstream and downstream signaling, Ras activation, and the Ras superfamily in this light, altogether providing a broad outlook of the principles and complexities of autoinhibition.

Evolution of In Silico Strategies for Protein-Protein Interaction Drug Discovery

The advent of advanced molecular modeling software, big data analytics, and high-speed processing units has led to the exponential evolution of modern drug discovery and better insights into complex biological processes and disease networks. This has progressively steered current research interests to understanding protein-protein interaction (PPI) systems that are related to a number of relevant diseases, such as cancer, neurological illnesses, metabolic disorders, etc. However, targeting PPIs are challenging due to their "undruggable" binding interfaces. In this review, we focus on the current obstacles that impede PPI drug discovery, and how recent discoveries and advances in in silico approaches can alleviate these barriers to expedite the search for potential leads, as shown in several exemplary studies. We will also discuss about currently available information on PPI compounds and systems, along with their usefulness in molecular modeling. Finally, we conclude by presenting the limits of in silico application in drug discovery and offer a perspective in the field of computer-aided PPI drug discovery.

Arl2-Mediated Allosteric Release of Farnesylated KRas4B from Shuttling Factor PDEδ

Proper localization of Ras proteins at the plasma membrane (PM) is crucial for their functions. To get to the PM, KRas4B and some other Ras family proteins bind to the PDEδ shuttling protein through their farnesylated hypervariable regions (HVRs). The docking of their farnesyl (and to a lesser extent geranylgeranyl) in the hydrophobic pocket of PDEδ's stabilizes the interaction. At the PM, guanosine 5'-triphosphate (GTP)-bound Arf-like protein 2 (Arl2) assists in the release of Ras from the PDEδ. However, exactly how is still unclear. Using all-atom molecular dynamics simulations, we unraveled the detailed mechanism of Arl2-mediated release of KRas4B, the most abundant oncogenic Ras isoform, from PDEδ. We simulated ternary Arl2-PDEδ-KRas4B HVR complexes and observed that Arl2 binding weakens the PDEδ-farnesylated HVR interaction. Our detailed analysis showed that allosteric changes (involving β6 of PDEδ and additional PDEδ residues) compress the hydrophobic PDEδ pocket and push the HVR out. Mutating PDEδ residues that mediate allosteric changes in PDEδ terminates the release process. Mutant Ras proteins are enriched in human cancers, with currently no drugs in the clinics. This mechanistic account may inspire efforts to develop drugs suppressing oncogenic KRas4B release.

Membrane Thickness as a Key Factor Contributing to the Activation of Osmosensors and Essential Ras Signaling Pathways

The cell membrane provides a functional link between the external environment and the replicating DNA genome by using ligand-gated receptors and chemical signals to activate signaling transduction pathways. However, increasing evidence has also indicated that the phospholipid bilayer itself by altering various physical parameters serves as a sensor that regulate membrane proteins in a specific manner. Changes in thickness and/or curvature of the membrane have been shown to be induced by mechanical forces and transmitted through the transmembrane helices of several types of mechanosensitive (MS) ion channels underlying functions such as osmoregulation in bacteria and sensory processing in mammalian cells. This review focus on recent protein functional and structural data indicating that the activation of bacterial and yeast osmosensors is consistent with thickness-induced tilting changes of the transmembrane domains of these proteins. Membrane thinning in combination with curvature changes may also lead to the lateral transfer of the small lipid-anchored GTPases Ras1 and H-Ras out of lipid rafts for clustering and signaling. The modulation of signaling pathways by amphiphilic peptides and the membrane-active antibiotics colistin and Amphotericin B is also discussed.