An orthogonalized PYR1-based CID module with reprogrammable ligand-binding specificity.

Plants sense abscisic acid (ABA) using chemical-induced dimerization (CID) modules, including the receptor PYR1 and HAB1, a phosphatase inhibited by ligand-activated PYR1. This system is unique because of the relative ease with which ligand recognition can be reprogrammed. To expand the PYR1 system, we designed an orthogonal '*' module, which harbors a dimer interface salt bridge; X-ray crystallographic, biochemical and in vivo analyses confirm its orthogonality. We used this module to create PYR1*MANDI/HAB1* and PYR1*AZIN/HAB1*, which possess nanomolar sensitivities to their activating ligands mandipropamid and azinphos-ethyl. Experiments in Arabidopsis thaliana and Saccharomyces cerevisiae demonstrate the sensitive detection of banned organophosphate contaminants using living biosensors and the construction of multi-input/output genetic circuits. Our new modules enable ligand-programmable multi-channel CID systems for plant and eukaryotic synthetic biology that can empower new plant-based and microbe-based sensing modalities.

Deep computational phenotyping of genomic variants impacting the SET domain of KMT2C reveal molecular mechanisms for their dysfunction.

Introduction: Kleefstra Syndrome type 2 (KLEFS-2) is a genetic, neurodevelopmental disorder characterized by intellectual disability, infantile hypotonia, severe expressive language delay, and characteristic facial appearance, with a spectrum of other distinct clinical manifestations. Pathogenic mutations in the epigenetic modifier type 2 lysine methyltransferase KMT2C have been identified to be causative in KLEFS-2 individuals. Methods: This work reports a translational genomic study that applies a multidimensional computational approach for deep variant phenotyping, combining conventional genomic analyses, advanced protein bioinformatics, computational biophysics, biochemistry, and biostatistics-based modeling. We use standard variant annotation, paralog annotation analyses, molecular mechanics, and molecular dynamics simulations to evaluate damaging scores and provide potential mechanisms underlying KMT2C variant dysfunction. Results: We integrated data derived from the structure and dynamics of KMT2C to classify variants into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). When compared with controls, these variants show values reflecting alterations in molecular fitness in both structure and dynamics. Discussion: We demonstrate that our 3D models for KMT2C variants suggest distinct mechanisms that lead to their imbalance and are not predictable from sequence alone. Thus, the missense variants studied here cause destabilizing effects on KMT2C function by different biophysical and biochemical mechanisms which we adeptly describe. This new knowledge extends our understanding of how variations in the KMT2C gene cause the dysfunction of its methyltransferase enzyme product, thereby bearing significant biomedical relevance for carriers of KLEFS2-associated genomic mutations.

A multi-layered computational structural genomics approach enhances domain-specific interpretation of Kleefstra syndrome variants in EHMT1.

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretations of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within the domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism for variant dysfunction.

A Multi-Layered Computational Structural Genomics Approach Enhances Domain-Specific Interpretation of Kleefstra Syndrome Variants in EHMT1.

This study investigates the functional significance of assorted variants of uncertain significance (VUS) in euchromatic histone lysine methyltransferase 1 (EHMT1), which is critical for early development and normal physiology. EHMT1 mutations cause Kleefstra syndrome and are linked to various human cancers. However, accurate functional interpretation of these variants are yet to be made, limiting diagnoses and future research. To overcome this, we integrate conventional tools for variant calling with computational biophysics and biochemistry to conduct multi-layered mechanistic analyses of the SET catalytic domain of EHMT1, which is critical for this protein function. We use molecular mechanics and molecular dynamics (MD)-based metrics to analyze the SET domain structure and functional motions resulting from 97 Kleefstra syndrome missense variants within this domain. Our approach allows us to classify the variants in a mechanistic manner into SV (Structural Variant), DV (Dynamic Variant), SDV (Structural and Dynamic Variant), and VUS (Variant of Uncertain Significance). Our findings reveal that the damaging variants are mostly mapped around the active site, substrate binding site, and pre-SET regions. Overall, we report an improvement for this method over conventional tools for variant interpretation and simultaneously provide a molecular mechanism of variant dysfunction.

Beyond structural bioinformatics for genomics with dynamics characterization of an expanded KRAS mutational landscape.

Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations, which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.

Fragment-based screening by protein-detected NMR spectroscopy.

Fragment-based drug discovery (FBDD) identifies low molecular weight compounds that can be developed into ligands with high affinity and selectivity for therapeutic targets. Screening fragment libraries (<10,000 molecules) with biophysical techniques against macromolecules provides information about novel chemical spaces that bind the macromolecule and scaffolds that can be modified to increase potency. A fragment-screening pipeline requires a standardized protocol for target selection, library assembly and maintenance, library screening, and hit validation to ensure hit integrity. Herein, the fundamental aspects of a fragment screening pipeline-focusing on protein-detected NMR data collection and analysis-are discussed in detail for researchers to use as a resource in their FBDD projects. Selected screening targets must undergo rigorous stability and buffer testing by NMR spectroscopy to ensure the protein structure is stable for the entire screen. Biophysical instrumentation that rapidly measures protein thermostability is helpful in buffer screening. Molecules in fragment libraries are analyzed computationally and physically, stored at appropriate temperatures, and multiplexed in well plates for library conservation. The screening protocol is streamlined using liquid handling robotics for sample preparation and customized Python scripts for protein-detected NMR data analysis. Molecules identified from the screen are titrated to determine their binding site(s) and Kd values and confirmed with an orthogonal biophysical assay. This detailed FBDD screening pipeline developed by the Program in Chemical Biology at the Medical College of Wisconsin has successfully screened many unrelated target proteins to identified novel molecules that selectively bind to these target proteins.

Design of stimulus-responsive two-state hinge proteins.

In nature, proteins that switch between two conformations in response to environmental stimuli structurally transduce biochemical information in a manner analogous to how transistors control information flow in computing devices. Designing proteins with two distinct but fully structured conformations is a challenge for protein design as it requires sculpting an energy landscape with two distinct minima. Here we describe the design of "hinge" proteins that populate one designed state in the absence of ligand and a second designed state in the presence of ligand. X-ray crystallography, electron microscopy, double electron-electron resonance spectroscopy, and binding measurements demonstrate that despite the significant structural differences the two states are designed with atomic level accuracy and that the conformational and binding equilibria are closely coupled.

Development and validation of a purification system for functional full-length human SR-B1 and CD36.

Scavenger receptor class B type 1 (SR-B1) and CD36 are both members of the class B scavenger receptor family that play important roles in lipoprotein metabolism and atherosclerotic disease. SR-B1 is the primary receptor for high-density lipoproteins, while CD36 is the receptor responsible for the internalization of oxidized low-density lipoproteins. Despite their importance, class B scavenger receptor structure has only been studied by functional domain or peptide fragments-there are currently no reports of utilizing purified full-length protein. Here we report the successful expression and purification of full-length human SR-B1 and CD36 using an Spodoptera frugiperda insect cell system. We demonstrate that both SR-B1 and CD36 retained their normal functions in Spodoptera frugiperda cells, including lipoprotein binding, lipid transport, and the formation of higher order oligomers in the plasma membrane. Purification schemes for both scavenger receptors were optimized and their purity was confirmed by SDS-PAGE. Both purified scavenger receptors were assessed for stability by thermal shift assay and shown to maintain stable melting temperatures up to 6 weeks post-purification. Microscale thermophoresis was used to demonstrate that purified SR-B1 and CD36 were able to bind their native lipoprotein ligands. Further, there was no difference in affinity of SR-B1 for high-density lipoprotein or CD36 for oxidized low-density lipoprotein, when comparing glycosylated and deglycosylated receptors. These studies mark a significant step forward in creating physiologically relevant tools to study scavenger receptor function and lay the groundwork for future functional studies and determination of receptor structure.

Metamorphic protein folding as evolutionary adaptation.

Metamorphic proteins switch reversibly between multiple distinct, stable structures, often with different functions. It was previously hypothesized that metamorphic proteins arose as intermediates in the evolution of a new fold - rare and transient exceptions to the 'one sequence, one fold' paradigm. However, as described herein, mounting evidence suggests that metamorphic folding is an adaptive feature, preserved and optimized over evolutionary time as exemplified by the NusG family and the chemokine XCL1. Analysis of extant protein families and resurrected protein ancestors demonstrates that large regions of sequence space are compatible with metamorphic folding. As a category that enhances biological fitness, metamorphic proteins are likely to employ fold switching to perform important biological functions and may be more common than previously thought.

Beyond Structural Bioinformatics for Genomics with Dynamics Characterization of an Expanded KRAS Mutational Landscape.

Current capabilities in genomic sequencing outpace functional interpretations. Our previous work showed that 3D protein structure calculations enhance mechanistic understanding of genetic variation in sequenced tumors and patients with rare diseases. The KRAS GTPase is among the critical genetic factors driving cancer and germline conditions. Because KRAS-altered tumors frequently harbor one of three classic hotspot mutations, nearly all studies have focused on these mutations, leaving significant functional ambiguity across the broader KRAS genomic landscape observed in cancer and non-cancer diseases. Herein, we extend structural bioinformatics with molecular simulations to study an expanded landscape of 86 KRAS mutations. We identify multiple coordinated changes strongly associated with experimentally established KRAS biophysical and biochemical properties. The patterns we observe span hotspot and non-hotspot alterations, which can all dysregulate Switch regions, producing mutation-restricted conformations with different effector binding propensities. We experimentally measured mutation thermostability and identified shared and distinct patterns with simulations. Our results indicate mutation-specific conformations which show potential for future research into how these alterations reverberate into different molecular and cellular functions. The data we present is not predictable using current genomic tools, demonstrating the added functional information derived from molecular simulations for interpreting human genetic variation.

Prospects for targeting ACKR1 in cancer and other diseases.

The chemokine network is comprised of a family of signal proteins that encode messages for cells displaying chemokine G-protein coupled receptors (GPCRs). The diversity of effects on cellular functions, particularly directed migration of different cell types to sites of inflammation, is enabled by different combinations of chemokines activating signal transduction cascades on cells displaying a combination of receptors. These signals can contribute to autoimmune disease or be hijacked in cancer to stimulate cancer progression and metastatic migration. Thus far, three chemokine receptor-targeting drugs have been approved for clinical use: Maraviroc for HIV, Plerixafor for hematopoietic stem cell mobilization, and Mogalizumab for cutaneous T-cell lymphoma. Numerous compounds have been developed to inhibit specific chemokine GPCRs, but the complexity of the chemokine network has precluded more widespread clinical implementation, particularly as anti-neoplastic and anti-metastatic agents. Drugs that block a single signaling axis may be rendered ineffective or cause adverse reactions because each chemokine and receptor often have multiple context-specific functions. The chemokine network is tightly regulated at multiple levels, including by atypical chemokine receptors (ACKRs) that control chemokine gradients independently of G-proteins. ACKRs have numerous functions linked to chemokine immobilization, movement through and within cells, and recruitment of alternate effectors like ß-arrestins. Atypical chemokine receptor 1 (ACKR1), previously known as the Duffy antigen receptor for chemokines (DARC), is a key regulator that binds chemokines involved in inflammatory responses and cancer proliferation, angiogenesis, and metastasis. Understanding more about ACKR1 in different diseases and populations may contribute to the development of therapeutic strategies targeting the chemokine network.

A modified ELISA assay differentiates CCL20 locked dimers from wild-type monomers.

A novel engineered CCL20 locked dimer (CCL20LD) is nearly identical to the naturally occurring chemokine CCL20 but blocks CCR6-mediated chemotaxis and offers a new approach to treat the diseases of psoriasis and psoriatic arthritis. Methods for quantifying CCL20LD serum levels are needed to assess pharmacokinetics parameters and evaluate drug delivery, metabolism, and toxicity. Existing ELISA kits fail to discriminate between CCL20LD and the natural chemokine, CCL20WT (the wild type monomer). Herein, we tested several available CCL20 monoclonal antibodies to be able to identify one clone that can be used both as a capture and a detection antibody (with biotin-labeling) to specifically detect CCL20LD with high specificity. After validation using recombinant proteins, the CCL20LD-selective ELISA was used to analyze blood samples from CCL20LD treated mice, demonstrating the utility of this novel assay for preclinical development of a biopharmaceutical lead compound for psoriatic disease.

De novo design of small beta barrel proteins.

Small beta barrel proteins are attractive targets for computational design because of their considerable functional diversity despite their very small size (<70 amino acids). However, there are considerable challenges to designing such structures, and there has been little success thus far. Because of the small size, the hydrophobic core stabilizing the fold is necessarily very small, and the conformational strain of barrel closure can oppose folding; also intermolecular aggregation through free beta strand edges can compete with proper monomer folding. Here, we explore the de novo design of small beta barrel topologies using both Rosetta energy-based methods and deep learning approaches to design four small beta barrel folds: Src homology 3 (SH3) and oligonucleotide/oligosaccharide-binding (OB) topologies found in nature and five and six up-and-down-stranded barrels rarely if ever seen in nature. Both approaches yielded successful designs with high thermal stability and experimentally determined structures with less than 2.4 Å rmsd from the designed models. Using deep learning for backbone generation and Rosetta for sequence design yielded higher design success rates and increased structural diversity than Rosetta alone. The ability to design a large and structurally diverse set of small beta barrel proteins greatly increases the protein shape space available for designing binders to protein targets of interest.

Mechanistic Insight into the Suppression of Polyglutamine Aggregation by SRCP1.

Protein aggregation is a hallmark of the polyglutamine diseases. One potential treatment for these diseases is suppression of polyglutamine aggregation. Previous work identified the cellular slime mold Dictyostelium discoideum as being naturally resistant to polyglutamine aggregation. Further work identified serine-rich chaperone protein 1 (SRCP1) as a protein that is both necessary in Dictyostelium and sufficient in human cells to suppress polyglutamine aggregation. Therefore, understanding how SRCP1 suppresses aggregation may be useful for developing therapeutics for the polyglutamine diseases. Here we utilized a de novo protein modeling approach to generate predictions of SRCP1's structure. Using our best-fit model, we generated mutants that were predicted to alter the stability of SRCP1 and tested these mutants' stability in cells. Using these data, we identified top models of SRCP1's structure that are consistent with the C-terminal region of SRCP1 forming a ß-hairpin with a highly dynamic N-terminal region. We next generated a series of peptides that mimic the predicted ß-hairpin and validated that they inhibit aggregation of a polyglutamine-expanded mutant huntingtin exon 1 fragment in vitro. To further assess mechanistic details of how SRCP1 inhibits polyglutamine aggregation, we utilized biochemical assays to determine that SRCP1 inhibits secondary nucleation in a manner dependent upon the regions flanking the polyglutamine tract. Finally, to determine if SRCP1 more could generally suppress protein aggregation, we confirmed that it was sufficient to inhibit aggregation of polyglutamine-expanded ataxin-3. Together these studies provide details into the structural and mechanistic basis of the inhibition of protein aggregation by SRCP1.

Fragment-based drug discovery of small molecule ligands for the human chemokine CCL28.

The mucosal chemokine CCL28 is a promising target for immunotherapy drug development due to its elevated expression level in epithelial cells and critical role in creating and maintaining an immunosuppressive tumor microenvironment. Using sulfotyrosine as a probe, NMR chemical shift mapping identified a potential receptor-binding hotspot on the human CCL28 surface. CCL28 was screened against 2,678 commercially available chemical fragments by 2D NMR, yielding thirteen verified hits. Computational docking predicted that two fragments could occupy adjoining subsites within the sulfotyrosine recognition cleft. Dual NMR titrations confirmed their ability to bind CCL28 simultaneously, thereby validating an initial fragment pair for linking and merging strategies to design high-potency CCL28 inhibitors.

DPP4-Truncated CXCL12 Alters CXCR4/ACKR3 Signaling, Osteogenic Cell Differentiation, Migration, and Senescence.

Bone marrow skeletal stem cells (SSCs) secrete many cytokines including stromal derived factor-1 or CXCL12, which influences cell proliferation, migration, and differentiation. All CXCL12 splice variants are rapidly truncated on their N-terminus by dipeptidyl peptidase 4 (DPP4). This includes the common variant CXCL12 alpha (1-68) releasing a much less studied metabolite CXCL12(3-68). Here, we found that CXCL12(3-68) significantly inhibited SSC osteogenic differentiation and RAW-264.7 cell osteoclastogenic differentiation and induced a senescent phenotype in SSCs. Importantly, pre-incubation of SSCs with CXCL12(3-68) significantly diminished their ability to migrate toward CXCL12(1-68) in transwell migration assays. Using a high-throughput G-protein-coupled receptor (GPCR) screen (GPCRome) and bioluminescent resonance energy transfer molecular interaction assays, we revealed that CXCL12(3-68) acts via the atypical cytokine receptor 3-mediated ß-arrestin recruitment and as a competitive antagonist to CXCR4-mediated signaling. Finally, a reverse phase protein array assay revealed that DPP4-cleaved CXCL12 possesses a different downstream signaling profile from that of intact CXCL12 or controls. The data presented herein provides insights into regulation of CXCL12 signaling. Importantly, it demonstrates that DPP4 proteolysis of CXCL12 generates a metabolite with significantly different and previously overlooked bioactivity that helps explain discrepancies in the literature. This also contributes to an understanding of the molecular mechanisms of osteoporosis and bone fracture repair and could potentially significantly affect the interpretation of experimental outcomes with clinical consequences in other fields where CXCL12 is vital, including cancer biology, immunology, cardiovascular biology, neurobiology, and associated pathologies.

Selective and Cell-Active PBRM1 Bromodomain Inhibitors Discovered through NMR Fragment Screening.

PBRM1 is a subunit of the PBAF chromatin remodeling complex that uniquely contains six bromodomains. PBRM1 can operate as a tumor suppressor or tumor promoter. PBRM1 is a tumor promoter in prostate cancer, contributing to migratory and immunosuppressive phenotypes. Selective chemical probes targeting PBRM1 bromodomains are desired to elucidate the association between aberrant PBRM1 chromatin binding and cancer pathogenesis and the contributions of PBRM1 to immunotherapy. Previous PBRM1 inhibitors unselectively bind SMARCA2 and SMARCA4 bromodomains with nanomolar potency. We used our protein-detected NMR screening pipeline to screen 1968 fragments against the second PBRM1 bromodomain, identifying 17 hits with Kd values from 45 µM to >2 mM. Structure-activity relationship studies on the tightest-binding hit resulted in nanomolar inhibitors with selectivity for PBRM1 over SMARCA2 and SMARCA4. These chemical probes inhibit the association of full-length PBRM1 to acetylated histone peptides and selectively inhibit growth of a PBRM1-dependent prostate cancer cell line.

Conformational selection guides ß-arrestin recruitment at a biased G protein-coupled receptor.

G protein-coupled receptors (GPCRs) recruit ß-arrestins to coordinate diverse cellular processes, but the structural dynamics driving this process are poorly understood. Atypical chemokine receptors (ACKRs) are intrinsically biased GPCRs that engage ß-arrestins but not G proteins, making them a model system for investigating the structural basis of ß-arrestin recruitment. Here, we performed nuclear magnetic resonance (NMR) experiments on 13CH3-ε-methionine-labeled ACKR3, revealing that ß-arrestin recruitment is associated with conformational exchange at key regions of the extracellular ligand-binding pocket and intracellular ß-arrestin-coupling region. NMR studies of ACKR3 mutants defective in ß-arrestin recruitment identified an allosteric hub in the receptor core that coordinates transitions among heterogeneously populated and selected conformational states. Our data suggest that conformational selection guides ß-arrestin recruitment by tuning receptor dynamics at intracellular and extracellular regions.

Rapid biosensor development using plant hormone receptors as reprogrammable scaffolds.

A general method to generate biosensors for user-defined molecules could provide detection tools for a wide range of biological applications. Here, we describe an approach for the rapid engineering of biosensors using PYR1 (Pyrabactin Resistance 1), a plant abscisic acid (ABA) receptor with a malleable ligand-binding pocket and a requirement for ligand-induced heterodimerization, which facilitates the construction of sense-response functions. We applied this platform to evolve 21 sensors with nanomolar to micromolar sensitivities for a range of small molecules, including structurally diverse natural and synthetic cannabinoids and several organophosphates. X-ray crystallography analysis revealed the mechanistic basis for new ligand recognition by an evolved cannabinoid receptor. We demonstrate that PYR1-derived receptors are readily ported to various ligand-responsive outputs, including enzyme-linked immunosorbent assay (ELISA)-like assays, luminescence by protein-fragment complementation and transcriptional circuits, all with picomolar to nanomolar sensitivity. PYR1 provides a scaffold for rapidly evolving new biosensors for diverse sense-response applications.

The Extended N-Terminal Domain Confers Atypical Chemokine Receptor Properties to CXCR3-B.

The chemokine receptor CXCR3 plays a critical role in immune cell recruitment and activation. CXCR3 exists as two main isoforms, CXCR3-A and CXCR3-B, resulting from alternative splicing. Although the two isoforms differ only by the presence of an N-terminal extension in CXCR3-B, they have been attributed divergent functional effects on cell migration and proliferation. CXCR3-B is the more enigmatic isoform and the mechanisms underlying its function and signaling remain elusive. We therefore undertook an in-depth cellular and molecular comparative study of CXCR3-A and CXCR3-B, investigating their activation at different levels of the signaling cascades, including G protein coupling, ß-arrestin recruitment and modulation of secondary messengers as well as their downstream gene response elements. We also compared the subcellular localization of the two isoforms and their trafficking under resting and stimulated conditions along with their ability to internalize CXCR3-related chemokines. Here, we show that the N-terminal extension of CXCR3-B drastically affects receptor features, modifying its cellular localization and preventing G protein coupling, while preserving ß-arrestin recruitment and chemokine uptake capacities. Moreover, we demonstrate that gradual truncation of the N terminus leads to progressive recovery of surface expression and G protein coupling. Our study clarifies the molecular basis underlying the divergent effects of CXCR3 isoforms, and emphasizes the ß-arrestin-bias and the atypical nature of CXCR3-B.