Bayesian network models identify cooperative GPCR:G protein interactions that contribute to G protein coupling

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network- like behavior and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces, it is challenging to determine which amino acid pair interactions are cooperative. In this work we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein coupled receptor (GPCR) complexes with Gα subunits. Our results reveal six GPCR:Gα contacts that are common to the different Gα subtypes and show strong cooperativity in the formation of interface. Both the C-terminus helix5 and the core of the G protein are codependent entities and play an important role in GPCR coupling. We show that a promiscuous GPCR coupling to different Gα subtypes, makes all the GPCR:Gα contacts that are specific to each Gα subtype (Gαs, Gαi and Gαq). This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.

Stabilization of interdomain interactions in G protein α subunits as a determinant of Gαi subtype signaling specificity

Highly homologous members of the Gαi family, Gαi1-3, have distinct tissue distributions and physiological functions, yet their biochemical and functional properties are very similar. We recently identified PDZ-RhoGEF (PRG) as a novel Gαi1 effector that is poorly activated by Gαi2. In a proteomic proximity labeling screen we observed a strong preference for Gαi1 relative to Gαi2 with respect to engagement of a broad range of potential targets. We investigated the mechanistic basis for this selectivity using PRG as a representative target. Substitution of either the helical domain (HD) from Gαi1 into Gαi2 or substitution of a single amino acid, A230 in Gαi2 with the corresponding D in Gαi1, largely rescues PRG activation and interactions with other potential Gαi targets. Molecular dynamics simulations combined with Bayesian network models revealed that in the GTP bound state, separation at the HD-Ras-like domain (RLD) interface is more pronounced in Gαi2 than Gαi1. Mutation of A230 to D in Gαi2 stabilizes HD-RLD interactions via ionic interactions with R145 in the HD which in turn modify the conformation of Switch III. These data support a model where D229 in Gαi1 interacts with R144 and stabilizes a network of interactions between HD and RLD to promote protein target recognition. The corresponding A230 in Gαi2 is unable to stabilize this network leading to an overall lower efficacy with respect to target interactions. This study reveals distinct mechanistic properties that could underly differential biological and physiological consequences of activation of Gαi1 or Gαi2 by G protein-coupled receptors.

Oligo-PROTAC strategy for cell-selective and targeted degradation of activated STAT3

Decoy oligodeoxynucleotides (ODNs) allow targeting undruggable transcription factors, such as STAT3, but their limited potency and lack of delivery methods hampered translation. To overcome these challenges, we conjugated a STAT3-specific decoy to thalidomide, a ligand to cereblon in E3 ubiquitin ligase complex, to generate a proteolysis-targeting chimera (STAT3DPROTAC). STAT3DPROTAC downregulated STAT3 in target cells, but not STAT1 or STAT5. Computational modeling of the STAT3DPROTAC ternary complex predicted two surface lysines, K601 and K626, in STAT3 as potential ubiquitination sites. Accordingly, K601/K626 point mutations in STAT3, as well as proteasome inhibition or cereblon deletion, alleviated STAT3DPROTAC effect. Next, we conjugated STAT3DPROTAC to a CpG oligonucleotide targeting Toll-like receptor 9 (TLR9) to generate myeloid/B cell-selective C-STAT3DPROTAC. Naked C-STAT3DPROTAC was spontaneously internalized by TLR9+ myeloid cells, B cells, and human and mouse lymphoma cells but not by T cells. C-STAT3DPROTAC effectively decreased STAT3 protein levels and also STAT3-regulated target genes critical for lymphoma cell proliferation and/or survival (BCL2L1, CCND2, and MYC). Finally, local C-STAT3DPROTAC administration to human Ly3 lymphoma-bearing mice triggered tumor regression, while control C-STAT3D and C-SCR treatments had limited effects. Our results underscore the feasibility of using a PROTAC strategy for cell-selective, decoy oligonucleotide-based STAT3 targeting of and potentially other tumorigenic transcription factors for cancer therapy.

A novel class of inhibitors that disrupts the stability of integrin heterodimers identified by CRISPR-tiling-instructed genetic screens

The plasma membrane is enriched for receptors and signaling proteins that are accessible from the extracellular space for pharmacological intervention. Here we conducted a series of CRISPR screens using human cell surface proteome and integrin family libraries in multiple cancer models. Our results identified ITGAV (integrin αV) and its heterodimer partner ITGB5 (integrin β5) as the essential integrin α/β pair for cancer cell expansion. High-density CRISPR gene tiling further pinpointed the integral pocket within the β-propeller domain of ITGAV for integrin αVβ5 dimerization. Combined with in silico compound docking, we developed a CRISPR-Tiling-Instructed Computer-Aided (CRISPR-TICA) pipeline for drug discovery and identified Cpd_AV2 as a lead inhibitor targeting the β-propeller central pocket of ITGAV. Cpd_AV2 treatment led to rapid uncoupling of integrin αVβ5 and cellular apoptosis, providing a unique class of therapeutic action that eliminates the integrin signaling via heterodimer dissociation. We also foresee the CRISPR-TICA approach to be an accessible method for future drug discovery studies.

Therapeutic targeting Tudor domains in leukemia via CRISPR-Scan Assisted Drug Discovery

Epigenetic dysregulation has been reported in multiple cancers including leukemias. Nonetheless, the roles of the epigenetic reader Tudor domains in leukemia progression and therapy remain unexplored. Here, we conducted a Tudor domain-focused CRISPR screen and identified SGF29, a component of SAGA/ATAC acetyltransferase complexes, as a crucial factor for H3K9 acetylation, ribosomal gene expression, and leukemogenesis. To facilitate drug development, we integrated the CRISPR tiling scan with compound docking and molecular dynamics simulation, presenting a generally applicable strategy called CRISPR-Scan Assisted Drug Discovery (CRISPR-SADD). Using this approach, we identified a lead inhibitor that selectively targets SGF29's Tudor domain and demonstrates efficacy against leukemia. Furthermore, we propose that the structural genetics approach used in our study can be widely applied to diverse fields for de novo drug discovery.

Molecular annotation of G protein variants in a neurological disorder

Heterotrimeric G proteins transduce extracellular chemical messages to generate appropriate intracellular responses. Point mutations in GNAO1, encoding the G protein αo subunit, have been implicated in a pathogenic condition characterized by seizures, movement disorders, intellectual disability, and developmental delay (GNAO1 disorder). However, the effects of these mutations on G protein structure and function are unclear. Here, we report the effects of 55 mutations on Gαo conformation, thermostability, nucleotide binding, and hydrolysis, as well as interaction with Gβγ subunits, receptors, and effectors. Our effort reveals four functionally distinct groups of mutants, including one group that sequesters receptors and another that sequesters Gβγ, both acting in a genetically dominant manner. These findings provide a more comprehensive understanding of disease-relevant mutations and reveal that GNAO1 disorder is likely composed of multiple mechanistically distinct disorders that will likely require multiple therapeutic strategies.

Cellular Lipids Regulate the Conformational Ensembles of the Disordered Intracellular Loop 3 in β2 Adrenergic Receptor

The structurally disordered intracellular loops (ICLs) of G protein-coupled receptors (GPCRs) play a critical role in G protein coupling. In our previous work, we used a combination of FRET-based and computational methodologies to show that the third intracellular loop (ICL3) modulates the activity and G protein coupling selectivity in GPCRs. In the current study, we have uncovered the role of several lipid components in modulating the conformational ensemble of ICL3 of the β2-adrenergic receptor (β2AR). Our findings indicate that phosphatidylinositol 4,5-bisphosphate (PIP2) in the inner leaflet of the membrane bilayer acts as a stabilizing anchor for ICL3, opening the intracellular cavity to facilitate G protein coupling. This interaction between PIP2 and ICL3 causes tilting of β2AR within the cellular membrane. Notably, this tilting of the receptor is supported by ganglioside GM3 stabilizing the extracellular loops on the outer leaflet of the bilayer, thereby exerting an allosteric effect on the orthosteric ligand binding pocket. Our results underscore the significance of lipids in modulating GPCR activity, proposing an allosteric mechanism that occurs through the receptor's orientation within the membrane.

Engineered odorant receptors illuminate structural principles of odor discrimination

A central challenge in olfaction is understanding how the olfactory system detects and distinguishes odorants with diverse physicochemical properties and molecular configurations. Vertebrate animals perceive odors via G protein-coupled odorant receptors (ORs). In humans, ~400 ORs enable the sense of smell. The OR family is composed of two major classes: Class I ORs are tuned to carboxylic acids while Class II ORs, representing the vast majority of the human repertoire, respond to a wide variety of odorants. How ORs recognize chemically diverse odorants remains poorly understood. A fundamental bottleneck is the inability to visualize odorant binding to ORs. Here, we uncover fundamental molecular properties of odorant-OR interactions by employing engineered ORs crafted using a consensus protein design strategy. Because such consensus ORs (consORs) are derived from the 17 major subfamilies of human ORs, they provide a template for modeling individual native ORs with high sequence and structural homology. The biochemical tractability of consORs enabled four cryoEM structures of distinct consORs with unique ligand recognition properties. The structure of a Class I consOR, consOR51, showed high structural similarity to the native human receptor OR51E2 and yielded a homology model of a related member of the human OR51 family with high predictive power. Structures of three Class II consORs revealed distinct modes of odorant-binding and activation mechanisms between Class I and Class II ORs. Thus, the structures of consORs lay the groundwork for understanding molecular recognition of odorants by the OR superfamily.

Bayesian network models identify co-operative GPCR:G protein interactions that contribute to G protein coupling

Cooperative interactions in protein-protein interfaces demonstrate the interdependency or the linked network-like behavior of interface interactions and their effect on the coupling of proteins. Cooperative interactions also could cause ripple or allosteric effects at a distance in protein-protein interfaces. Although they are critically important in protein-protein interfaces it is challenging to determine which amino acid pair interactions are cooperative. In this work we have used Bayesian network modeling, an interpretable machine learning method, combined with molecular dynamics trajectories to identify the residue pairs that show high cooperativity and their allosteric effect in the interface of G protein-coupled receptor (GPCR) complexes with G proteins. Our results reveal a strong co-dependency in the formation of interface GPCR:G protein contacts. This observation indicates that cooperativity of GPCR:G protein interactions is necessary for the coupling and selectivity of G proteins and is thus critical for receptor function. We have identified subnetworks containing polar and hydrophobic interactions that are common among multiple GPCRs coupling to different G protein subtypes (Gs, Gi and Gq). These common subnetworks along with G protein-specific subnetworks together confer selectivity to the G protein coupling. This work underscores the potential of data-driven Bayesian network modeling in elucidating the intricate dependencies and selectivity determinants in GPCR:G protein complexes, offering valuable insights into the dynamic nature of these essential cellular signaling components.

Unraveling allostery within the angiotensin II type 1 receptor for Gαq and β-arrestin coupling

G protein-coupled receptors engage both G proteins and β-arrestins, and their coupling can be biased by ligands and mutations. Here, to resolve structural elements and mechanisms underlying effector coupling to the angiotensin II (AngII) type 1 receptor (AT1R), we combined alanine scanning mutagenesis of the entire sequence of the receptor with pharmacological profiling of Gαq and β-arrestin engagement to mutant receptors and molecular dynamics simulations. We showed that Gαq coupling to AT1R involved a large number of residues spread across the receptor, whereas fewer structural regions of the receptor contributed to β-arrestin coupling regulation. Residue stretches in transmembrane domain 4 conferred β-arrestin bias and represented an important structural element in AT1R for functional selectivity. Furthermore, we identified allosteric small-molecule binding sites that were enclosed by communities of residues that produced biased signaling when mutated. Last, we showed that allosteric communication within AT1R emanating from the Gαq coupling site spread beyond the orthosteric AngII-binding site and across different regions of the receptor, including currently unresolved structural regions. Our findings reveal structural elements and mechanisms within AT1R that bias Gαq and β-arrestin coupling and that could be harnessed to design biased receptors for research purposes and to develop allosteric modulators.

Oligo-PROTAC strategy for cell-selective and targeted degradation of activated STAT3

Decoy-oligodeoxynucleotides (D-ODNs) can target undruggable transcription factors, such as STAT3. However, challenges in D-ODN delivery and potency hampered their translation. To overcome these limitations, we conjugated STAT3-specific D-ODN to thalidomide (Tha), a known ligand to cereblon (CRBN, a component of E3 ubiquitin ligase) to generate a proteolysis-targeting chimera (STAT3D PROTAC ). STAT3D PROTAC downregulated STAT3, but not STAT1 or STAT5, in target cells. Computational modeling of the STAT3D PROTAC ternary complex predicted two surface lysines on STAT3, K601 and K626 as potential ubiquitination sites for the PROTAC bound E3 ligase. Accordingly, K601/K626 point mutations in STAT3, as well as proteasome inhibitors, and CRBN deletion alleviated STAT3D PROTAC effect. Next, we conjugated STAT3D PROTAC to a CpG ligand targeting Toll-like receptor 9 (TLR9) to generate myeloid/B-cell-selective C-STAT3D PROTAC conjugate. Naked C-STAT3D PROTAC was spontaneously internalized by TLR9 + myeloid cells, B cells as well as human Ly18 and mouse A20 lymphoma cells, but not by T cells. C-STAT3D PROTAC decreased STAT3 levels to 50% at 250 nM and over 85% at 2 µM dosing in myeloid cells. We also observed significantly improved downregulation of STAT3 target genes involved in lymphoma cell proliferation and/or survival ( BCL2L1, CCND2, MYC ). Finally, we assessed the antitumor efficacy of C-STAT3D PROTAC compared to C-STAT3D or scrambled control (C-SCR) against human lymphoma xenotransplants. Local C-STAT3D PROTAC administration triggered lymphoma regression while control treatments had limited effects. Our results underscore feasibility of using PROTAC strategy for cell-selective, decoy oligonucleotide-based targeting of STAT3 and potentially other tumorigenic transcription factors for cancer therapy.

Molecular recognition of an aversive odorant by the murine trace amine-associated receptor TAAR7f

There are two main families of G protein-coupled receptors that detect odours in humans, the odorant receptors (ORs) and the trace amine-associated receptors (TAARs). Their amino acid sequences are distinct, with the TAARs being most similar to the aminergic receptors such as those activated by adrenaline, serotonin and histamine. To elucidate the structural determinants of ligand recognition by TAARs, we have determined the cryo-EM structure of a murine receptor, mTAAR7f, coupled to the heterotrimeric G protein Gs and bound to the odorant N,N-dimethylcyclohexylamine (DMCH) to an overall resolution of 2.9 Å. DMCH is bound in a hydrophobic orthosteric binding site primarily through van der Waals interactions and a strong charge-charge interaction between the tertiary amine of the ligand and an aspartic acid residue. This site is distinct and non-overlapping with the binding site for the odorant propionate in the odorant receptor OR51E2. The structure, in combination with mutagenesis data and molecular dynamics simulations suggests that the activation of the receptor follows a similar pathway to that of the β-adrenoceptors, with the significant difference that DMCH interacts directly with one of the main activation microswitch residues.

Stabilization of Interdomain Interactions in G protein αi Subunits Determines Gαi Subtype Signaling Specificity

Highly homologous members of the Gαi family, Gαi1-3, have distinct tissue distributions and physiological functions, yet the functional properties of these proteins with respect to GDP/GTP binding and regulation of adenylate cyclase are very similar. We recently identified PDZ-RhoGEF (PRG) as a novel Gαi1 effector, however, it is poorly activated by Gαi2. Here, in a proteomic proximity labeling screen we observed a strong preference for Gαi1 relative to Gαi2 with respect to engagement of a broad range of potential targets. We investigated the mechanistic basis for this selectivity using PRG as a representative target. Substitution of either the helical domain (HD) from Gαi1 into Gαi2 or substitution of a single amino acid, A230 in Gαi2 to the corresponding D in Gαi1, largely rescues PRG activation and interactions with other Gαi targets. Molecular dynamics simulations combined with Bayesian network models revealed that in the GTP bound state, dynamic separation at the HD-Ras-like domain (RLD) interface is prevalent in Gαi2 relative to Gαi1 and that mutation of A230s4h3.3 to D in Gαi2 stabilizes HD-RLD interactions through formation of an ionic interaction with R145HD.11 in the HD. These interactions in turn modify the conformation of Switch III. These data support a model where D229s4h3.3 in Gαi1 interacts with R144HD.11 stabilizes a network of interactions between HD and RLD to promote protein target recognition. The corresponding A230 in Gαi2 is unable to form the "ionic lock" to stabilize this network leading to an overall lower efficacy with respect to target interactions. This study reveals distinct mechanistic properties that could underly differential biological and physiological consequences of activation of Gαi1 or Gαi2 by GPCRs.

Dynamic phosphatase-recruitment controls B-cell selection and oncogenic signaling

Initiation of B-cell receptor (BCR) ¹ signaling, and subsequent antigen-encounter in germinal centers ^2,3 represent milestones of B-lymphocyte development that are both marked by sharp increases of CD25 surface-expression. Oncogenic signaling in B-cell leukemia (B-ALL) ⁴ and lymphoma ⁵ also induced CD25-surface expression. While CD25 is known as an IL2-receptor chain on T- and NK-cells ^6-9 , the significance of its expression on B-cells was unclear. Our experiments based on genetic mouse models and engineered patient-derived xenografts revealed that, rather than functioning as an IL2-receptor chain, CD25 expressed on B-cells assembled an inhibitory complex including PKCδ and SHIP1 and SHP1 phosphatases for feedback control of BCR-signaling or its oncogenic mimics. Recapitulating phenotypes of genetic ablation of PKCδ ¹⁰ - ¹² , SHIP1 ^13,14 and SHP1 ^{14, 15,16} , conditional CD25-deletion decimated early B-cell subsets but expanded mature B-cell populations and induced autoimmunity. In B-cell malignancies arising from early (B-ALL) and late (lymphoma) stages of B-cell development, CD25-loss induced cell death in the former and accelerated proliferation in the latter. Clinical outcome annotations mirrored opposite effects of CD25-deletion: high CD25 expression levels predicted poor clinical outcomes for patients with B-ALL, in contrast to favorable outcomes for lymphoma-patients. Biochemical and interactome studies revealed a critical role of CD25 in BCR-feedback regulation: BCR-signaling induced PKCδ-mediated phosphorylation of CD25 on its cytoplasmic tail (S ²⁶⁸ ). Genetic rescue experiments identified CD25-S ²⁶⁸ tail-phosphorylation as central structural requirement to recruit SHIP1 and SHP1 phosphatases to curb BCR-signaling. A single point mutation CD25 ^S268A abolished recruitment and activation of SHIP1 and SHP1 to limit duration and strength of BCR-signaling. Loss of phosphatase-function, autonomous BCR-signaling and Ca ²⁺ -oscillations induced anergy and negative selection during early B-cell development, as opposed to excessive proliferation and autoantibody production in mature B-cells. These findings highlight the previously unrecognized role of CD25 in assembling inhibitory phosphatases to control oncogenic signaling in B-cell malignancies and negative selection to prevent autoimmune disease.

A Closer Look at EGFR Inhibitor Resistance in Non-Small Cell Lung Cancer through the Lens of Precision Medicine

The development of EGFR small-molecule inhibitors has provided significant benefit for the affected patient population. Unfortunately, current inhibitors are no curative therapy, and their development has been driven by on-target mutations that interfere with binding and thus inhibitory activity. Genomic studies have revealed that, in addition to these on-target mutations, there are also multiple off-target mechanisms of EGFR inhibitor resistance and novel therapeutics that can overcome these challenges are sought. Resistance to competitive 1st-generation and covalent 2nd- and 3rd-generation EGFR inhibitors is overall more complex than initially thought, and novel 4th-generation allosteric inhibitors are expected to suffer from a similar fate. Additional nongenetic mechanisms of resistance are significant and can include up to 50% of the escape pathways. These potential targets have gained recent interest and are usually not part of cancer panels that look for alterations in resistant patient specimen. We discuss the duality between genetic and nongenetic EGFR inhibitor drug resistance and summarize current team medicine approaches, wherein clinical developments, hand in hand with drug development research, drive potential opportunities for combination therapy.

Autoregulation of GPCR signalling through the third intracellular loop

The third intracellular loop (ICL3) of the G protein-coupled receptor (GPCR) fold is important for the signal transduction process downstream of receptor activation1-3. Despite this, the lack of a defined structure of ICL3, combined with its high sequence divergence among GPCRs, complicates characterization of its involvement in receptor signalling4. Previous studies focusing on the β2 adrenergic receptor (β2AR) suggest that ICL3 is involved in the structural process of receptor activation and signalling5-7. Here we derive mechanistic insights into the role of ICL3 in β2AR signalling, observing that ICL3 autoregulates receptor activity through a dynamic conformational equilibrium between states that block or expose the receptor's G protein-binding site. We demonstrate the importance of this equilibrium for receptor pharmacology, showing that G protein-mimetic effectors bias the exposed states of ICL3 to allosterically activate the receptor. Our findings additionally reveal that ICL3 tunes signalling specificity by inhibiting receptor coupling to G protein subtypes that weakly couple to the receptor. Despite the sequence diversity of ICL3, we demonstrate that this negative G protein-selection mechanism through ICL3 extends to GPCRs across the superfamily, expanding the range of known mechanisms by which receptors mediate G protein subtype selective signalling. Furthermore, our collective findings suggest ICL3 as an allosteric site for receptor- and signalling pathway-specific ligands.

Structural basis of odorant recognition by a human odorant receptor

Our sense of smell enables us to navigate a vast space of chemically diverse odour molecules. This task is accomplished by the combinatorial activation of approximately 400 odorant G protein-coupled receptors encoded in the human genome1-3. How odorants are recognized by odorant receptors remains unclear. Here we provide mechanistic insight into how an odorant binds to a human odorant receptor. Using cryo-electron microscopy, we determined the structure of the active human odorant receptor OR51E2 bound to the fatty acid propionate. Propionate is bound within an occluded pocket in OR51E2 and makes specific contacts critical to receptor activation. Mutation of the odorant-binding pocket in OR51E2 alters the recognition spectrum for fatty acids of varying chain length, suggesting that odorant selectivity is controlled by tight packing interactions between an odorant and an odorant receptor. Molecular dynamics simulations demonstrate that propionate-induced conformational changes in extracellular loop 3 activate OR51E2. Together, our studies provide a high-resolution view of chemical recognition of an odorant by a vertebrate odorant receptor, providing insight into how this large family of G protein-coupled receptors enables our olfactory sense.

Catalytic site mutations confer multiple states of G protein activation

Heterotrimeric guanine nucleotide-binding proteins (G proteins) that function as molecular switches for cellular growth and metabolism are activated by GTP and inactivated by GTP hydrolysis. In uveal melanoma, a conserved glutamine residue critical for GTP hydrolysis in the G protein α subunit is often mutated in Gαq or Gα11 to either leucine or proline. In contrast, other glutamine mutations or mutations in other Gα subtypes are rare. To uncover the mechanism of the genetic selection and the functional role of this glutamine residue, we analyzed all possible substitutions of this residue in multiple Gα isoforms. Through cell-based measurements of activity, we showed that some mutants were further activated and inactivated by G protein-coupled receptors. Through biochemical, molecular dynamics, and nuclear magnetic resonance-based structural studies, we showed that the Gα mutants were functionally distinct and conformationally diverse, despite their shared inability to hydrolyze GTP. Thus, the catalytic glutamine residue contributes to functions beyond GTP hydrolysis, and these functions include subtype-specific, allosteric modulation of receptor-mediated subunit dissociation. We conclude that G proteins do not function as simple on-off switches. Rather, signaling emerges from an ensemble of active states, a subset of which are favored in disease and may be uniquely responsive to receptor-directed ligands.

Biased agonists differentially modulate the receptor conformation ensembles in Angiotensin II type 1 receptor

The structural features that contribute to the efficacy of biased agonists targeting G protein-coupled receptors (GPCRs) towards G proteins or β-arrestin (β-arr) signaling pathways is nebulous, although such knowledge is critical in designing biased ligands. The dynamics of the agonist-GPCR complex is one of the critical factors in determining agonist bias. Angiotensin II type I receptor (AT1R) is an ideal model system to study the molecular basis of bias since it has multiple β-arr2 and Gq protein biased agonists as well as experimentally solved three dimensional structures. Using Molecular Dynamics (MD) simulations for the Angiotensin II type I receptor (AT1R) bound to ten different agonists, we infer that the agonist bound receptor samples conformations with different relative weights, from both the inactive and active state ensembles of the receptor. This concept is perhaps extensible to other class A GPCRs. Such a weighted mixed ensemble recapitulates the inter-residue distance distributions measured for different agonists bound AT1R using DEER experiments. The ratio of the calculated relative strength of the allosteric communication to β-arr2 vs Gq coupling sites scale similarly to the experimentally measured bias factors. Analysis of the inter-residue distance distributions of the activation microswitches involved in class A GPCR activation suggests that β-arr2 biased agonists turn on different combination of microswitches with different relative strengths of activation. We put forth a model that activation microswitches behave like rheostats that tune the relative efficacy of the biased agonists toward the two signaling pathways. Finally, based on our data we propose that the agonist specific residue contacts in the binding site elicit a combinatorial response in the microswitches that in turn differentially modulate the receptor conformation ensembles resulting in differences in coupling to Gq and β-arrestin.

Dynamic spatiotemporal determinants modulate GPCR:G protein coupling selectivity and promiscuity

Recent studies have shown that G protein coupled receptors (GPCRs) show selective and promiscuous coupling to different Gα protein subfamilies and yet the mechanisms of the range of coupling preferences remain unclear. Here, we use Molecular Dynamics (MD) simulations on ten GPCR:G protein complexes and show that the location (spatial) and duration (temporal) of intermolecular contacts at the GPCR:Gα protein interface play a critical role in how GPCRs selectively interact with G proteins. We identify that some GPCR:G protein interface contacts are common across Gα subfamilies and others specific to Gα subfamilies. Using large scale data analysis techniques on the MD simulation snapshots we derive a spatio-temporal code for contacts that confer G protein selective coupling and validated these contacts using G protein activation BRET assays. Our results demonstrate that promiscuous GPCRs show persistent sampling of the common contacts more than G protein specific contacts. These findings suggest that GPCRs maintain contact with G proteins through a common central interface, while the selectivity comes from G protein specific contacts at the periphery of the interface.

Utilization of model-agnostic explainable artificial intelligence frameworks in oncology: a narrative review

Background and Objective

Machine learning (ML) models are increasingly being utilized in oncology research for use in the clinic. However, while more complicated models may provide improvements in predictive or prognostic power, a hurdle to their adoption are limits of model interpretability, wherein the inner workings can be perceived as a "black box". Explainable artificial intelligence (XAI) frameworks including Local Interpretable Model-agnostic Explanations (LIME) and SHapley Additive exPlanations (SHAP) are novel, model-agnostic approaches that aim to provide insight into the inner workings of the "black box" by producing quantitative visualizations of how model predictions are calculated. In doing so, XAI can transform complicated ML models into easily understandable charts and interpretable sets of rules, which can give providers with an intuitive understanding of the knowledge generated, thus facilitating the deployment of such models in routine clinical workflows.

Methods

We performed a comprehensive, non-systematic review of the latest literature to define use cases of model-agnostic XAI frameworks in oncologic research. The examined database was PubMed/MEDLINE. The last search was run on May 1, 2022.

Key Content and Findings

In this review, we identified several fields in oncology research where ML models and XAI were utilized to improve interpretability, including prognostication, diagnosis, radiomics, pathology, treatment selection, radiation treatment workflows, and epidemiology. Within these fields, XAI facilitates determination of feature importance in the overall model, visualization of relationships and/or interactions, evaluation of how individual predictions are produced, feature selection, identification of prognostic and/or predictive thresholds, and overall confidence in the models, among other benefits. These examples provide a basis for future work to expand on, which can facilitate adoption in the clinic when the complexity of such modeling would otherwise be prohibitive.

Conclusions

Model-agnostic XAI frameworks offer an intuitive and effective means of describing oncology ML models, with applications including prognostication and determination of optimal treatment regimens. Using such frameworks presents an opportunity to improve understanding of ML models, which is a critical step to their adoption in the clinic.

Activation mechanism of the class D fungal GPCR dimer Ste2

The fungal class D1 G-protein-coupled receptor (GPCR) Ste2 has a different arrangement of transmembrane helices compared with mammalian GPCRs and a distinct mode of coupling to the heterotrimeric G protein Gpa1-Ste2-Ste181. In addition, Ste2 lacks conserved sequence motifs such as DRY, PIF and NPXXY, which are associated with the activation of class A GPCRs2. This suggested that the activation mechanism of Ste2 may also differ. Here we determined structures of Saccharomyces cerevisiae Ste2 in the absence of G protein in two different conformations bound to the native agonist α-factor, bound to an antagonist and without ligand. These structures revealed that Ste2 is indeed activated differently from other GPCRs. In the inactive state, the cytoplasmic end of transmembrane helix H7 is unstructured and packs between helices H1-H6, blocking the G protein coupling site. Agonist binding results in the outward movement of the extracellular ends of H6 and H7 by 6 Å. On the intracellular surface, the G protein coupling site is formed by a 20 Å outward movement of the unstructured region in H7 that unblocks the site, and a 12 Å inward movement of H6. This is a distinct mechanism in GPCRs, in which the movement of H6 and H7 upon agonist binding facilitates G protein coupling.

Sequence coevolution and structure stabilization modulate olfactory receptor expression

Olfactory receptors (ORs) belong to class A G-protein coupled receptors (GPCRs) and are activated by a variety of odorants. To date, there is no three-dimensional structure of an OR available. One of the major bottlenecks in obtaining purified protein for structural studies of ORs is their poor expression in heterologous cells. To design mutants that enhance expression and thereby enable protein purification, we first identified computable physical properties that recapitulate OR and class A GPCR expression and further conducted an iterative computational prediction-experimental test cycle and generated human OR mutants that express as high as biogenic amine receptors for which structures have been solved. In the process of developing the computational method to recapitulate the expression of ORs in membranes, we identified properties, such as amino acid sequence coevolution, and the strength of the interactions between intracellular loop 1 (ICL1) and the helix 8 region of ORs, to enhance their heterologous expression. We identified mutations that are directly located in these regions as well as other mutations not located in these regions but allosterically strengthen the ICL1-helix 8 enhance expression. These mutants also showed functional responses to known odorants. This method to enhance heterologous expression of mammalian ORs will facilitate high-throughput "deorphanization" of ORs, and enable OR purification for biochemical and structural studies to understand odorant-OR interactions.

Allosteric communication regulates ligand-specific GPCR activity

G protein-coupled receptors (GPCRs) are membrane-bound proteins that are ubiquitously expressed in many cell types and take part in mediating multiple signaling pathways. GPCRs are dynamic proteins and exist in an equilibrium between an ensemble of conformational states such as inactive and fully active states. This dynamic nature of GPCRs is one of the factors that confers their basal activity even in the absence of any ligand-mediated activation. Ligands selectively bind and stabilize a subset of the conformations from the ensemble leading to a shift in the equilibrium toward the inactive or the active state depending on the nature of the ligand. This ligand-selective effect is achieved through allosteric communication between the ligand binding site and G protein or β-arrestin coupling site. Similarly, the G protein coupling to the receptor exerts the allosteric effect on the ligand binding region leading to increased binding affinity for agonists and decreased affinity for antagonists or inverse agonists. In this review, we enumerate the current state of our understanding of the mechanism of allosteric communication in GPCRs with a specific focus on the critical role of computational methods in delineating the residues involved in allosteric communication. Analyzing allosteric communication mechanism using molecular dynamics simulations has revealed (a) a structurally conserved mechanism of allosteric communication that regulates the G protein coupling, (b) a rational structure-based approach to designing selective ligands, and (c) an approach to designing allosteric GPCR mutants that are either ligand and G protein or β-arrestin selective.

A universal allosteric mechanism for G protein activation

G proteins play a central role in signal transduction and pharmacology. Signaling is initiated by cell-surface receptors, which promote guanosine triphosphate (GTP) binding and dissociation of Gα from the Gβγ subunits. Structural studies have revealed the molecular basis of subunit association with receptors, RGS proteins, and downstream effectors. In contrast, the mechanism of subunit dissociation is poorly understood. We use cell signaling assays, molecular dynamics (MD) simulations, and biochemistry and structural analyses to identify a conserved network of amino acids that dictates subunit release. In the presence of the terminal phosphate of GTP, a glycine forms a polar network with an arginine and glutamate, putting torsional strain on the subunit binding interface. This "G-R-E motif" secures GTP and, through an allosteric link, discharges the Gβγ dimer. Replacement of network residues prevents subunit dissociation regardless of agonist or GTP binding. These findings reveal the molecular basis of the final committed step of G protein activation.

Kinase inhibitors allosterically disrupt a regulatory interaction to enhance PKCα membrane translocation

The eukaryotic kinase domain has multiple intrinsically disordered regions whose conformation dictates kinase activity. Small molecule kinase inhibitors (SMKIs) rely on disrupting the active conformations of these disordered regions to inactivate the kinase. While SMKIs are selected for their ability to cause this disruption, the allosteric effects of conformational changes in disordered regions is limited by a lack of dynamic information provided by traditional structural techniques. In this study, we integrated multiscale molecular dynamics simulations with FRET sensors to characterize a novel allosteric mechanism that is selectively triggered by SMKI binding to the protein kinase Cα domain. The indole maleimide inhibitors BimI and sotrastaurin were found to displace the Gly-rich loop (G-loop) that normally shields the ATP-binding site. Displacement of the G-loop interferes with a newly identified, structurally conserved binding pocket for the C1a domain on the N lobe of the kinase domain. This binding pocket, in conjunction with the N-terminal regulatory sequence, masks a diacylglycerol (DAG) binding site on the C1a domain. SMKI-mediated displacement of the G-loop released C1a and exposed the DAG binding site, enhancing protein kinase Cα translocation both to synthetic lipid bilayers and to live cell membranes in the presence of DAG. Inhibitor chemotype determined the extent of the observed allosteric effects on the kinase domain and correlated with the extent of membrane recruitment. Our findings demonstrate the allosteric effects of SMKIs beyond the confines of kinase catalytic conformation and provide an integrated computational-experimental paradigm to investigate parallel mechanisms in other kinases.

Structure of the class D GPCR Ste2 dimer coupled to two G proteins

G-protein-coupled receptors (GPCRs) are divided phylogenetically into six classes^1,2, denoted A to F. More than 370 structures of vertebrate GPCRs (belonging to classes A, B, C and F) have been determined, leading to a substantial understanding of their function³. By contrast, there are no structures of class D GPCRs, which are found exclusively in fungi where they regulate survival and reproduction. Here we determine the structure of a class D GPCR, the Saccharomyces cerevisiae pheromone receptor Ste2, in an active state coupled to the heterotrimeric G protein Gpa1-Ste4-Ste18. Ste2 was purified as a homodimer coupled to two G proteins. The dimer interface of Ste2 is formed by the N terminus, the transmembrane helices H1, H2 and H7, and the first extracellular loop ECL1. We establish a class D1 generic residue numbering system (CD1) to enable comparisons with orthologues and with other GPCR classes. The structure of Ste2 bears similarities in overall topology to class A GPCRs, but the transmembrane helix H4 is shifted by more than 20 Å and the G-protein-binding site is a shallow groove rather than a cleft. The structure provides a template for the design of novel drugs to target fungal GPCRs, which could be used to treat numerous intractable fungal diseases⁴.

IFITM3 functions as a PIP3 scaffold to amplify PI3K signalling in B cells

Interferon-induced transmembrane protein 3 (IFITM3) has previously been identified as an endosomal protein that blocks viral infection^1-3. Here we studied clinical cohorts of patients with B cell leukaemia and lymphoma, and identified IFITM3 as a strong predictor of poor outcome. In normal resting B cells, IFITM3 was minimally expressed and mainly localized in endosomes. However, engagement of the B cell receptor (BCR) induced both expression of IFITM3 and phosphorylation of this protein at Tyr20, which resulted in the accumulation of IFITM3 at the cell surface. In B cell leukaemia, oncogenic kinases phosphorylate IFITM3 at Tyr20, which causes constitutive localization of this protein at the plasma membrane. In a mouse model, Ifitm3^-/- naive B cells developed in normal numbers; however, the formation of germinal centres and the production of antigen-specific antibodies were compromised. Oncogenes that induce the development of leukaemia and lymphoma did not transform Ifitm3^-/- B cells. Conversely, the phosphomimetic IFITM3(Y20E) mutant induced oncogenic PI3K signalling and initiated the transformation of premalignant B cells. Mechanistic experiments revealed that IFITM3 functions as a PIP3 scaffold and central amplifier of PI3K signalling. The amplification of PI3K signals depends on IFITM3 using two lysine residues (Lys83 and Lys104) in its conserved intracellular loop as a scaffold for the accumulation of PIP3. In Ifitm3^-/- B cells, lipid rafts were depleted of PIP3, which resulted in the defective expression of over 60 lipid-raft-associated surface receptors, and impaired BCR signalling and cellular adhesion. We conclude that the phosphorylation of IFITM3 that occurs after B cells encounter antigen induces a dynamic switch from antiviral effector functions in endosomes to a PI3K amplification loop at the cell surface. IFITM3-dependent amplification of PI3K signalling, which in part acts downstream of the BCR, is critical for the rapid expansion of B cells with high affinity to antigen. In addition, multiple oncogenes depend on IFITM3 to assemble PIP3-dependent signalling complexes and amplify PI3K signalling for malignant transformation.

Activation Microswitches in Adenosine Receptor A2A Function as Rheostats in the Cell Membrane

Although multiple components of the cell membrane modulate the stability and activation of G protein-coupled receptors (GPCRs), insights into the dynamics of GPCR structures come from biophysical studies conducted in detergents. This is because of the challenges of studying activation in a multicomponent lipid bilayer. To understand the role of cellular membrane lipids and cations in GPCR activation, we performed multiscale molecular dynamics simulations (56 μs) on three different conformational states of adenosine receptor A2AR, in both the cell membrane-like lipid bilayer and in detergent micelles. Molecular dynamics (MD) simulations show that the phosphatidylinositol bisphosphate (PIP2) interacts with the basic residues in the intracellular regions of A2AR, thereby reducing the flexibility of the receptor in the inactive state and limiting the transition to the active-intermediate state. In the G protein-coupled fully active state, PIP2 stabilizes the GPCR:G protein complex. Such stiffening effects are absent in non-ionic detergent micelles, and therefore, more transitions have been observed in detergents. The inter-residue distances that change significantly upon GPCR activation are known as activation microswitches. The activation microswitches show different levels of activation in the cell membrane, in the pure POPC bilayer, and in detergents. Thus, the temporal heat map of different activation microswitches calculated from the MD simulations suggests a rheostat model of GPCR activation microswitches rather than the binary switch model. These simulation results connect the chemistry of cell membrane lipids to receptor activity, which is useful for the design of detergents mimicking the cell membrane.

Structural characterization of a dimeric complex between the short cytoplasmic domain of CEACAM1 and the pseudo tetramer of S100A10-Annexin A2 using NMR and molecular dynamics

AIIt, a heterotetramer of S100A10 (P11) and Annexin A2, plays a key role in calcium dependent, membrane associations with a variety of proteins. We previously showed that AIIt interacts with the short cytoplasmic domain (12 amino acids) of CEACAM1 (CEACAM1-SF). Since the cytoplasmic domains of CEACAM1 help regulate the formation of cis- or trans-dimers at the cell membrane, we investigated the possible role of their association with AIIt in this process. Using NMR and molecular dynamics, we show that AIIt and its pseudoheterodimer interacts with two molecules of short cytoplasmic domain isoform peptides, and that interaction depends on the binding motif 454-Phe-Gly-Lys-Thr-457 where Phe-454 binds in a hydrophobic pocket of AIIt, the null mutation Phe454Ala reduces binding by 2.5 fold, and the pseudophosphorylation mutant Thr457Glu reduces binding by three fold. Since these two residues in CEACAM1-SF were also found to play a role in the binding of calmodulin and G-actin at the membrane, we hypothesize a sequential set of three interactions are responsible for regulation of cis- to trans-dimerization of CEACAM1. The hydrophobic binding pocket in AIIt corresponds to a previously identified binding pocket for a peptide found in SMARCA3 and AHNAK, suggesting a conserved functional motif in AIIt allowing multiple proteins to reversibly interact with integral membrane proteins in a calcium dependent manner.

How Do Branched Detergents Stabilize GPCRs in Micelles?

The structural and functional properties of G protein-coupled receptors (GPCRs) are often studied in a detergent micellar environment, but many GPCRs tend to denature or aggregate in short alkyl chain detergents. In our previous work [Lee, S., et al. (2016) J. Am. Chem. Soc.138, 15425–15433], we showed that GPCRs in alkyl glucosides were highly dynamic, resulting in the penetration of detergent molecules between transmembrane α-helices, which is the initial step in receptor denaturation. Although this was not observed for GPCRs in dodecyl maltoside (DDM, also known as lauryl maltoside), even this detergent is not mild enough to preserve the integrity of many GPCRs during purification. Lauryl maltose neopentylglycol (LMNG) detergents have been found to have significant advantages for purifying GPCRs in a native state as they impart more stability to the receptor than DDM. To gain insights into how they stabilize GPCRs, we used atomistic molecular dynamics simulations of wild type adenosine A_2A receptor (WT-A_2AR), thermostabilized A_2AR (tA_2AR), and wild type β₂-adrenoceptor (β₂AR) in a variety of detergents (LMNG, DMNG, OGNG, and DDM). Analysis of molecular dynamics simulations of tA_2AR in LMNG, DMNG, and OGNG showed that this series of detergents exhibited behavior very similar to that of an analogous series of detergents DDM, DM, and OG in our previous study. However, there was a striking difference upon comparison of the behavior of LMNG to that of DDM. LMNG showed considerably less motion than DDM, which resulted in the enhanced density of the aliphatic chains around the hydrophobic regions of the receptor and considerably more hydrogen bond formation between the head groups. This contributed to enhanced interaction energies between both detergent molecules and between the receptor and detergent, explaining the enhanced stability of GPCRs purified in this detergent. Branched detergents occlude between transmembrane helices and reduce their flexibility. Our results provide a rational foundation to develop detergent variants for stabilizing membrane proteins.

Structural instability and divergence from conserved residues underlie intracellular retention of mammalian odorant receptors

Mammalian odorant receptors are a diverse and rapidly evolving set of G protein-coupled receptors expressed in olfactory cilia membranes. Most odorant receptors show little to no cell surface expression in nonolfactory cells due to endoplasmic reticulum retention, which has slowed down biochemical studies. Here we provide evidence that structural instability and divergence from conserved residues of individual odorant receptors underlie intracellular retention using a combination of large-scale screening of odorant receptors cell surface expression in heterologous cells, point mutations, structural modeling, and machine learning techniques. We demonstrate the importance of conserved residues by synthesizing consensus odorant receptors that show high levels of cell surface expression similar to conventional G protein-coupled receptors. Furthermore, we associate in silico structural instability with poor cell surface expression using molecular dynamics simulations. We propose an enhanced evolutionary capacitance of olfactory sensory neurons that enable the functional expression of odorant receptors with cryptic mutations.

Machine Learning for Prioritization of Thermostabilizing Mutations for G-Protein Coupled Receptors

Although the three-dimensional structures of G-protein coupled receptors (GPCRs), the largest superfamily of drug targets, have enabled structure-based drug design, there are no structures available for 87% of GPCRs. This is due to the stiff challenge in purifying the inherently flexible GPCRs. Identifying thermostabilized mutant GPCRs via systematic alanine scanning mutations has been a successful strategy in stabilizing GPCRs, but it remains a daunting task for each GPCR. We developed a computational method that combines sequence-, structure-, and dynamics-based molecular properties of GPCRs that recapitulate GPCR stability, with four different machine learning methods to predict thermostable mutations ahead of experiments. This method has been trained on thermostability data for 1231 mutants, the largest publicly available data set. A blind prediction for thermostable mutations of the complement factor C5a receptor 1 retrieved 36% of the thermostable mutants in the top 50 prioritized mutants compared to 3% in the first 50 attempts using systematic alanine scanning.

Prediction of Conformation Specific Thermostabilizing Mutations for Class A G Protein-Coupled Receptors

G protein-coupled receptors (GPCRs) are highly flexible and prone to denaturation during protein extraction in detergents and purification. This poses a huge challenge to purify a conformationally homogeneous solution of GPCRs. Thermostabilizing mutations have been used widely to purify and obtain crystal structures of several GPCRs. However, identifying thermostabilizing mutations for GPCRs remains a tedious and expensive task as they are not transferable even among closely related GPCRs. Additionally, the mutations stabilizing one conformational state of a GPCR do not always stabilize other conformational state(s) of the same GPCR. Previously we developed a computational method, LiticonDesign, for rapid prediction of thermostabilizing mutations for a specific GPCR conformation. In this study, we have used LiticonDesign to predict thermostabilizing mutations for the agonist bound active-intermediate state of the human adenosine receptor (A_2AR) using the structure of the inactive state of the same GPCR and vice versa. Our study shows that the thermostable mutation predictions using LiticonDesign, for an active-intermediate state of a GPCR (A_2AR in our case), requires a homology model that is derived from an active/active-intermediate state GPCR structure as a template. Similarly, the homology models derived from inactive state GPCR conformations are better in predicting the thermostable mutations for the inactive state of A_2AR. Overall, LiticonDesign method is not only efficient in predicting thermostabilizing mutations for a given GPCR sequence but also can recover conformation specific mutations for a state of interest, if a suitable starting structure of desired conformation is chosen.

Generation of dual specific bivalent BiTEs (dbBIspecific T-cell engaging antibodies) for cellular immunotherapy

Background

Bispecific T-cell engaging antibodies (BiTES), comprising dual anti-CD3 and anti-tumor antigen scFv fragments, are important therapeutic agents for the treatment of cancer. The dual scFv construct for BiTES requires proper protein folding while their small molecular size leads to rapid kidney clearance.

Methods

An intact (150 kDa) anti-tumor antigen antibody to CEA was joined in high yield (ca. 30%) to intact (150 kDa) anti-murine and anti-human CD3 antibodies using hinge region specific Click chemistry to form dual-specific, bivalent BiTES (dbBiTES, 300 kDa). dbBiTEs were tested in vitro by EM, flow cytometry and cell cytoxicity and in vivo by PET tumor imaging and redirected T-cell therapy.

Results

The interlocked hinge regions are compatible with a structural model that fits the electron micrographs of 300 kDa particles. Compared to intact anti-CEA antibody, dbBiTES exhibit high in vitro cytotoxicity, high in vivo tumor targeting as demonstrated by PET imaging, and redirected dbBiTE coated T-cells (1 microgram/10 million cells) that kill CEA⁺ target cells in vivo in CEA transgenic mice.

Conclusion

dbBiTE redirected T-cell therapy is a promising, efficient approach for targeting and killing cancer cells.

Electronic supplementary material

The online version of this article (10.1186/s12885-019-6056-8) contains supplementary material, which is available to authorized users.

Conformational plasticity of the intracellular cavity of GPCR-G-protein complexes leads to G-protein promiscuity and selectivity

While the dynamics of the intracellular surface in agonist-stimulated GPCRs is well studied, the impact of GPCR dynamics on G-protein selectivity remains unclear. Here, we combine molecular dynamics simulations with live-cell FRET and secondary messenger measurements, for 21 GPCR-G-protein combinations, to advance a dynamic model of the GPCR-G-protein interface. Our data show C terminus peptides of Gαs, Gαi, and Gαq proteins assume a small ensemble of unique orientations when coupled to their cognate GPCRs, similar to the variations observed in 3D structures of GPCR-G-protein complexes. The noncognate G proteins interface with latent intracellular GPCR cavities but dissociate due to weak and unstable interactions. Three predicted mutations in β2-adrenergic receptor stabilize binding of noncognate Gαq protein in its latent cavity, allowing promiscuous signaling through both Gαs and Gαq in a dose-dependent manner. This demonstrates that latent GPCR cavities can be evolved, by design or nature, to tune G-protein selectivity, giving insights to pluridimensional GPCR signaling.

Efficient and Selective Electrochemically Driven Enzyme-Catalyzed Reduction of Carbon Dioxide to Formate using Formate Dehydrogenase and an Artificial Cofactor

Increasing levels of carbon dioxide in the atmosphere and the growing need for energy necessitate a shift toward reliance on renewable energy sources and the utilization of carbon dioxide. Thus, producing carbonaceous fuel by the electrochemical reduction of carbon dioxide has been very appealing. We have focused on addressing the principal challenges of poor selectivity and poor energy efficiency in the electrochemical reduction of carbon dioxide. We have demonstrated here a viable pathway for the efficient and continuous electrochemical reduction of CO₂ to formate using the metal-independent enzyme type of formate dehydrogenase (FDH) derived from C andida boidinii yeast. This type of FDH is attractive because it is commercially produced. In natural metabolic processes, this type of metal-independent FDH oxidizes formate to carbon dioxide using NAD⁺ as a cofactor. We show that FDH can catalyze the reverse process to generate formate when the natural cofactor NADH is replaced with an artificial cofactor, the methyl viologen radical cation. The methyl viologen radical cation is generated in situ, electrochemically. Our approach relies on the special properties of methyl viologen as a "unidirectional" redox cofactor for the conversion of CO₂ to formate. Methyl viologen (in the oxidized form) does not catalyze formate oxidation, while the methyl viologen radical cation is an effective cofactor for the reduction of carbon dioxide. Thus, although the thermodynamic driving force is favorable for the oxidized form of methyl viologen to oxidize formate to carbon dioxide, the kinetic factors are not favorable. Only the reverse reaction of carbon dioxide reduction to formate is kinetically viable with the cofactor, methyl viologen radical cation. Binding free energy calculated from atomistic molecular dynamics (MD) simulations consolidate our understanding of these special binding properties of the methyl viologen radical cation and its ability to facilitate the two-electron reduction of carbon dioxide to formate in metal-independent FDH. By carrying out the reactions in a novel three-compartment cell, we have demonstrated the continuous production of formate at high energy efficiency and yield. This cell configuration uses judiciously selected ion-exchange membranes to separate the reaction compartments to preserve the yields of the methyl viologen radical cation and formate. By the electroregeneration of the methyl viologen radical cation at -0.44 V versus the normal hydrogen electrode, we could produce formate at 20 mV negative to the reversible electrode potential for carbon dioxide reduction to formate. Our results are in sharp contrast to the large overpotentials of -800 to -1000 mV required on metal catalysts, vindicating the selectivity and kinetic facility provided by FDH. Formate yields as high as 97% ± 1% could be realized by avoiding the adventitious reoxidation of the methyl viologen radical cation by molecular oxygen. We anticipate that the insights from the electrochemical studies and the MD simulations to be useful in redesigning the metal-independent FDH and alternate artificial cofactors to achieve even higher rates of conversion.

Dynamic Role of the G Protein in Stabilizing the Active State of the Adenosine A2A Receptor

Agonist binding in the extracellular region of the G protein-coupled adenosine A2A receptor increases its affinity to the G proteins in the intracellular region, and vice versa. The structural basis for this effect is not evident from the crystal structures of A_2AR in various conformational states since it stems from the receptor dynamics. Using atomistic molecular dynamics simulations on four different conformational states of the adenosine A_2A receptor, we observed that the agonists show decreased ligand mobility, lower entropy of the extracellular loops in the active-intermediate state compared with the inactive state. In contrast, the entropy of the intracellular region increases to prime the receptor for coupling the G protein. Coupling of the G protein to A_2AR shrinks the agonist binding site, making tighter receptor agonist contacts with an increase in the strength of allosteric communication compared with the active-intermediate state. These insights provide a strong basis for structure-based ligand design studies.

•

GPCR conformation dynamics reveals the forward and backward allosteric mechanism

•

Agonist binding increases the entropy in the intracellular region of the GPCR

•

G protein binding shrinks the receptor-ligand contacts in the extracellular region

•

Increased allostery between G protein and agonist in the GPCR-G protein complex

NMR analysis of free and lipid nanodisc anchored CEACAM1 membrane proximal peptides with Ca2+/CaM

CEACAM1, a homotypic transmembrane receptor with 12 or 72 amino acid cytosolic domain isoforms, is converted from inactive cis-dimers to active trans-dimers by calcium-calmodulin (Ca²⁺/CaM). Previously, the weak binding of Ca²⁺/CaM to the human 12 AA cytosolic domain was studied using C-terminal anchored peptides. We now show the binding of ¹⁵N labeled Phe-454 cytosolic domain peptides in solution or membrane anchored using NMR demonstrates a significant role for the lipid bilayer. Although binding is increased by the mutation Phe454Ala, this mutation was previously shown to abrogate actin binding. On the other hand, Ca²⁺/CaM binding is abrogated by phosphorylation of nearby Thr-457, a post-translation modification required for actin binding and subsequent in vitro lumen formation. Binding of Ca²⁺/CaM to a membrane proximal peptide from the long 72 AA cytosolic domain anchored to lipid nanodiscs was very weak compared to lipid free conditions, suggesting membrane specific effects between the two isoforms. NMR analysis of ¹⁵N labeled Ca²⁺/CaM with unlabeled peptides showed the C-lobe of Ca²⁺/CaM is involved in peptide interactions, and hydrophobic residues such as Met-109, Val-142 and Met-144 play important roles in binding peptide. This information was incorporated into transmembrane models of CEACAM1 binding to Ca²⁺/CaM. The lack of Ca²⁺/CaM binding to phosphorylated Thr-457, a residue we have previously shown to be phosphorylated by CaMK2D, also dependent on Ca²⁺/CaM, suggests stepwise binding of the cytosolic domain first to Ca²⁺/CaM and then to actin.

Engineering Salt Bridge Networks between Transmembrane Helices Confers Thermostability in G-Protein-Coupled Receptors

Introduction of specific point mutations has been an effective strategy in enhancing the thermostability of G-protein-coupled receptors (GPCRs). Our previous work showed that a specific residue position on transmembrane helix 6 (TM6) in class A GPCRs consistently yields thermostable mutants. The crystal structure of human chemokine receptor CCR5 also showed increased thermostability upon mutation of two positions, A233D^6.33 and K303E^7.59. With the goal of testing the transferability of these two thermostabilizing mutations in other chemokine receptors, we tested the mutations A237D^6.33 and R307E^7.59 in human CCR3 for thermostability and aggregation properties in detergent solution. Interestingly, the double mutant exhibited a 6-10-fold decrease in the aggregation propensity of the wild-type protein. This is in stark contrast to the two single mutants whose aggregation properties resemble the wild type (WT). Moreover, unlike in CCR5, the two single mutants separately showed no increase in thermostability compared to the wild-type CCR3, while the double-mutant A237D^6.33/R307E^7.59 confers an increase of 2.6 °C in the melting temperature compared to the WT. Extensive all-atom molecular dynamics (MD) simulations in detergent micelles show that a salt bridge network between transmembrane helices TM3, TM6, and TM7 that is absent in the two single mutants confers stability in the double mutant. The free energy surface of the double mutant shows conformational homogeneity compared to the single mutants. An annular n-dodecyl maltoside detergent layer packs tighter to the hydrophobic surface of the double-mutant CCR3 compared to the single mutants providing additional stability. The purification of other C-C chemokine receptors lacking such stabilizing residues may benefit from the incorporation of these two point mutations.

Bitopic Inhibition of ATP and Substrate Binding in Ser/Thr Kinases through a Conserved Allosteric Mechanism

Protein kinases achieve substrate selective phosphorylation through their conformational flexibility and dynamic interaction with the substrate. Designing substrate selective or kinase selective small molecule inhibitors remains a challenge because of a lack of understanding of the dynamic mechanism by which substrates are selected by the kinase. Using a combination of all-atom molecular dynamics simulations and FRET sensors, we have delineated an allosteric mechanism that results in interaction among the DFG motif, G-loop, and activation loop and structurally links the nucleotide and substrate binding interfaces in protein kinase Cα and three other Ser/Thr kinases. ATP-competitive staurosporine analogues engage this allosteric switch region located just outside the ATP binding site to displace substrate binding to varying degrees. These inhibitors function as bitopic ligands by occupying the ATP binding site and interacting with the allosteric switch region. The conserved mechanism identified in this study can be exploited to select and design bitopic inhibitors for kinases.

Allosteric Activation of Striatal-Enriched Protein Tyrosine Phosphatase (STEP, PTPN5) by a Fragment-like Molecule

Protein tyrosine phosphatase non-receptor type 5 (PTPN5, STEP) is a brain specific phosphatase that regulates synaptic function and plasticity by modulation of N-methyl-d-aspartate receptor (NMDAR) and α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPAR) trafficking. Dysregulation of STEP has been linked to neurodegenerative and neuropsychiatric diseases, highlighting this enzyme as an attractive therapeutic target for drug discovery. Selective targeting of STEP with small molecules has been hampered by high conservation of the active site among protein tyrosine phosphatases. We report the discovery of the first small molecule allosteric activator for STEP that binds to the phosphatase domain. Allosteric binding is confirmed by both X-ray and ¹⁵N NMR experiments, and specificity has been demonstrated by an enzymatic test cascade. Molecular dynamics simulations indicate stimulation of enzymatic activity by a long-range allosteric mechanism. To allow the scientific community to make use of this tool, we offer to provide the compound in the course of an open innovation initiative.

C-NHEJ without indels is robust and requires synergistic function of distinct XLF domains

To investigate the fidelity of canonical non-homologous end joining (C-NHEJ), we developed an assay to detect EJ between distal ends of two Cas9-induced chromosomal breaks that are joined without causing insertion/deletion mutations (indels). Here we find that such EJ requires several core C-NHEJ factors, including XLF. Using variants of this assay, we find that C-NHEJ is required for EJ events that use 1-2, but not ≥3, nucleotides of terminal microhomology. We also investigated XLF residues required for EJ without indels, finding that one of two binding domains is essential (L115 or C-terminal lysines that bind XRCC4 and KU/DNA, respectively), and that disruption of one of these domains sensitizes XLF to mutations that affect its dimer interface, which we examined with molecular dynamic simulations. Thus, C-NHEJ, including synergistic function of distinct XLF domains, is required for EJ of chromosomal breaks without indels.