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.

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 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.

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.

The dimeric form of CXCL12 binds to atypical chemokine receptor 1.

The pleiotropic chemokine CXCL12 is involved in diverse physiological and pathophysiological processes, including embryogenesis, hematopoiesis, leukocyte migration, and tumor metastasis. It is known to engage the classical receptor CXCR4 and the atypical receptor ACKR3. Differential receptor engagement can transduce distinct cellular signals and effects as well as alter the amount of free, extracellular chemokine. CXCR4 binds both monomeric and the more commonly found dimeric forms of CXCL12, whereas ACKR3 binds monomeric forms. Here, we found that CXCL12 also bound to the atypical receptor ACKR1 (previously known as Duffy antigen/receptor for chemokines or DARC). In vitro nuclear magnetic resonance spectroscopy and isothermal titration calorimetry revealed that dimeric CXCL12 bound to the extracellular N terminus of ACKR1 with low nanomolar affinity, whereas the binding affinity of monomeric CXCL12 was orders of magnitude lower. In transfected MDCK cells and primary human Duffy-positive erythrocytes, a dimeric, but not a monomeric, construct of CXCL12 efficiently bound to and internalized with ACKR1. This interaction between CXCL12 and ACKR1 provides another layer of regulation of the multiple biological functions of CXCL12. The findings also raise the possibility that ACKR1 can bind other dimeric chemokines, thus potentially further expanding the role of ACKR1 in chemokine retention and presentation.

Targeted biologic inhibition of both tumor cell-intrinsic and intercellular CLPTM1L/CRR9-mediated chemotherapeutic drug resistance.

Recurrence of therapy-resistant tumors is a principal problem in solid tumor oncology, particularly in ovarian cancer. Despite common complete responses to first line, platinum-based therapies, most women with ovarian cancer recur, and eventually, nearly all with recurrent disease develop platinum resistance. Likewise, both intrinsic and acquired resistance contribute to the dismal prognosis of pancreatic cancer. Our previous work and that of others has established CLPTM1L (cleft lip and palate transmembrane protein 1-like)/CRR9 (cisplatin resistance related protein 9) as a cytoprotective oncofetal protein that is present on the tumor cell surface. We show that CLPTM1L is broadly overexpressed and accumulated on the plasma membrane of ovarian tumor cells, while weakly or not expressed in normal tissues. High expression of CLPTM1L is associated with poor outcome in ovarian serous adenocarcinoma. Robust re-sensitization of resistant ovarian cancer cells to platinum-based therapy was achieved using human monoclonal biologics inhibiting CLPTM1L in both orthotopic isografts and patient-derived cisplatin resistant xenograft models. Furthermore, we demonstrate that in addition to cell-autonomous cytoprotection by CLPTM1L, extracellular CLPTM1L confers resistance to chemotherapeutic killing in an ectodomain-dependent fashion, and that this intercellular resistance mechanism is inhibited by anti-CLPTM1L biologics. Specifically, exosomal CLPTM1L from cisplatin-resistant ovarian carcinoma cell lines conferred resistance to cisplatin in drug-sensitive parental cell lines. CLPTM1L is present in extracellular vesicle fractions of tumor culture supernatants and in patients' serum with increasing abundance upon chemotherapy treatment. These findings have encouraging implications for the use of anti-CLPTM1L targeted biologics in the treatment of therapy-resistant tumors.

Bacterial expression of the phosphodiester-binding site of the cation-independent mannose 6-phosphate receptor for crystallographic and NMR studies.

The cation-independent mannose 6-phosphate receptor (CI-MPR) is a multifunctional protein that interacts with diverse ligands and plays central roles in autophagy, development, and tumor suppression. By delivering newly synthesized phosphomannosyl-containing acid hydrolases from the Golgi to endosomal compartments, CI-MPR is an essential component in the generation of lysosomes that are critical for the maintenance of cellular homeostasis. The ability of CI-MPR to interact with ∼60 different acid hydrolases is facilitated by its large extracellular region, with four out of its 15 domains binding phosphomannosyl residues. Although the glycan specificity of CI-MPR has been elucidated, the molecular basis of carbohydrate binding has not been determined for two out of these four carbohydrate recognition domains (CRD). Here we report expression of CI-MPR's CRD located in domain 5 that preferentially binds phosphodiester-containing glycans. Domain 5 of CI-MPR was expressed in Escherichia coli BL21 (DE3) cells as a fusion protein containing an N-terminal histidine tag and the small ubiquitin-like modifier (SUMO) protein. The His6-SUMO-CRD construct was recovered from inclusion bodies, refolded in buffer to facilitate disulfide bond formation, and subjected to Ni-NTA affinity chromatography and size exclusion chromatography. Surface plasmon resonance analyses demonstrated that the purified protein was active and bound phosphorylated glycans. Characterization by NMR spectroscopy revealed high quality (1)H-(15)N HSQC spectra. Additionally, crystallization conditions were identified and a crystallographic data set of the CRD was collected to 1.8Šresolution. Together, these studies demonstrate the feasibility of producing CI-MPR's CRD suitable for three-dimensional structure determination by NMR spectroscopic and X-ray crystallographic approaches.

Recombinant leukemia inhibitory factor suppresses human medullary thyroid carcinoma cell line xenografts in mice.

Medullary thyroid carcinoma (MTC) is a neoplasm of the endocrine system, which originates from parafollicular C-cells of the thyroid gland. For MTC therapy, the Food and Drug Administration recently approved vandetanib and cabozantinib, multi-kinase inhibitors targeting RET and other tyrosine kinase receptors of vascular endothelial growth factor, epidermal growth factor, or hepatocyte growth factor. Nevertheless, not all patients with the progressive MTC respond to these drugs, requiring the development of additional therapeutic modalities that have distinct activity. Previously, we reported that expression of activated Ras or Raf in the human MTC cell lines, TT and MZ-CRC-1, can induce growth arrest and RET downregulation via a leukemia inhibitory factor (LIF)-mediated autocrine/paracrine loop. In this study, we aimed to evaluate bacterially-produced recombinant human LIF for its efficacy to suppress human MTC xenografts in mice. Here, we report that, consistent with its effects in vitro, locally or systemically administered recombinant LIF effectively suppressed growth of TT and MZ-CRC-1 xenografts in mice. Further, as predicted from its effects in TT and MZ-CRC-1 cell cultures in vitro, recombinant LIF activated the JAK/STAT pathway and downregulated RET and E2F1 expression in tumors in mice. These results suggest that LIF is a potent cytostatic agent for MTC cells, which regulates unique mechanisms that are not targeted by currently available therapeutic agents.

Solution structure of CCL21 and identification of a putative CCR7 binding site.

CCL21 is a human chemokine that recruits normal immune cells and metastasizing tumor cells to lymph nodes through activation of the G protein-coupled receptor CCR7. The CCL21 structure solved by NMR contains a conserved chemokine domain followed by an extended, unstructured C-terminus that is not typical of most other chemokines. A sedimentation equilibrium study showed CCL21 to be monomeric. Chemical shift mapping indicates that the CCR7 N-terminus binds to the N-loop and third ß-strand of CCL21's chemokine domain. Details of CCL21-receptor recognition may enable structure-based drug discovery of novel antimetastatic agents.

Structure of the SCAN domain from the tumor suppressor protein MZF1.

The SCAN domain mediates interactions between members of a subfamily of zinc-finger transcription factors and is found in more than 60 C2H2 zinc finger genes in the human genome, including the tumor suppressor gene myeloid zinc finger 1 (MZF1). Glutathione-S-transferase pull-down assays showed that the MZF1 SCAN domain self-associates, and a Kd value of 600 nM was measured by intrinsic tryptophan fluorescence polarization. The MZF1 structure determined by NMR spectroscopy revealed a domain-swapped dimer. Each monomer consists of five alpha helices in two subdomains connected by the alpha2-alpha3 loop. Residues from helix 3 of each monomer compose the core of the dimer interface, while the alpha1-alpha2 loop and helix 2 pack against helices 3 and 5 from the opposing monomer. Comprehensive sequence analysis is coupled with the first high-resolution structure of a SCAN dimer to provide an initial view of the recognition elements that govern dimerization for this large family of transcription factors.