Computational modeling reveals key molecular properties and dynamic behavior of disruptor of telomeric silencing 1-like (DOT1L) and partnering complexes involved in leukemogenesis.

Disruptor of telomeric silencing 1-like (DOT1L) is the only non-SET domain histone lysine methyltransferase (KMT) and writer of H3K79 methylation on nucleosomes marked by H2B ubiquitination. DOT1L has elicited significant attention because of its interaction or fusion with members of the AF protein family in blood cell biology and leukemogenic transformation. Here, our goal was to extend previous structural information by performing a robust molecular dynamic study of DOT1L and its leukemogenic partners combined with mutational analysis. We show that statically and dynamically, D161, G163, E186, and F223 make frequent time-dependent interactions with SAM, while additional residues T139, K187, and N241 interact with SAM only under dynamics. Dynamics models reveal DOT1L, SAM, and H4 moving as one and show that more than twice the number of DOT1L residues interacts with these partners, relative to the static structure. Mutational analyses indicate that six of these residues are intolerant to substitution. We describe the dynamic behavior of DOT1L interacting with AF10 and AF9. Studies on the dynamics of a heterotrimeric complex of DOT1L1-AF10 illuminated describe coordinated motions that impact the relative position of the DOT1L HMT domain to the nucleosome. The molecular motions of the DOT1L-AF9 complex are less extensive and highly dynamic, resembling a swivel-like mechanics. Through molecular dynamics and mutational analysis, we extend the knowledge previous provided by static measurements. These results are important to consider when describing the biochemical properties of DOT1L, under normal and in disease conditions, as well as for the development of novel therapeutic agents.

NADPH oxidase 4 contributes to TRPV4-mediated endothelium-dependent vasodilation in human arterioles by regulating protein phosphorylation of TRPV4 channels.

Impaired endothelium-dependent vasodilation has been suggested to be a key component of coronary microvascular dysfunction (CMD). A better understanding of endothelial pathways involved in vasodilation in human arterioles may provide new insight into the mechanisms of CMD. The goal of this study is to investigate the role of TRPV4, NOX4, and their interaction in human arterioles and examine the underlying mechanisms. Arterioles were freshly isolated from adipose and heart tissues obtained from 71 patients without coronary artery disease, and vascular reactivity was studied by videomicroscopy. In human adipose arterioles (HAA), ACh-induced dilation was significantly reduced by TRPV4 inhibitor HC067047 and by NOX 1/4 inhibitor GKT137831, but GKT137831 did not further affect the dilation in the presence of TRPV4 inhibitors. GKT137831 also inhibited TRPV4 agonist GSK1016790A-induced dilation in HAA and human coronary arterioles (HCA). NOX4 transcripts and proteins were detected in endothelial cells of HAA and HCA. Using fura-2 imaging, GKT137831 significantly reduced GSK1016790A-induced Ca2+ influx in the primary culture of endothelial cells and TRPV4-WT-overexpressing human coronary artery endothelial cells (HCAEC). However, GKT137831 did not affect TRPV4-mediated Ca2+ influx in non-phosphorylatable TRPV4-S823A/S824A-overexpressing HCAEC. In addition, treatment of HCAEC with GKT137831 decreased the phosphorylation level of Ser824 in TRPV4. Finally, proximity ligation assay (PLA) revealed co-localization of NOX4 and TRPV4 proteins. In conclusion, both TRPV4 and NOX4 contribute to ACh-induced dilation in human arterioles from patients without coronary artery disease. NOX4 increases TRPV4 phosphorylation in endothelial cells, which in turn enhances TRPV4-mediated Ca2+ entry and subsequent endothelium-dependent dilation in human arterioles.

SOD2 acetylation on lysine 68 promotes stem cell reprogramming in breast cancer.

Mitochondrial superoxide dismutase (SOD2) suppresses tumor initiation but promotes invasion and dissemination of tumor cells at later stages of the disease. The mechanism of this functional switch remains poorly defined. Our results indicate that as SOD2 expression increases acetylation of lysine 68 ensues. Acetylated SOD2 promotes hypoxic signaling via increased mitochondrial reactive oxygen species (mtROS). mtROS, in turn, stabilize hypoxia-induced factor 2α (HIF2α), a transcription factor upstream of "stemness" genes such as Oct4, Sox2, and Nanog. In this sense, our findings indicate that SOD2K68Ac and mtROS are linked to stemness reprogramming in breast cancer cells via HIF2α signaling. Based on these findings we propose that, as tumors evolve, the accumulation of SOD2K68Ac turns on a mitochondrial pathway to stemness that depends on HIF2α and may be relevant for the progression of breast cancer toward poor outcomes.

Inorganic arsenic promotes luminal to basal transition and metastasis of breast cancer.

Inorganic arsenic (iAs/As2 O32- ) is an environmental toxicant found in watersheds around the world including in densely populated areas. iAs is a class I carcinogen known to target the skin, lungs, bladder, and digestive organs, but its role as a primary breast carcinogen remains controversial. Here, we examined a different possibility: that exposure to iAs promotes the transition of well-differentiated epithelial breast cancer cells characterized by estrogen and progesterone receptor expression (ER+/PR+), to more basal phenotypes characterized by active proliferation, and propensity to metastasis in vivo. Our results indicate two clear phenotypic responses to low-level iAs that depend on the duration of the exposure. Short-term pulses of iAs activate ER signaling, consistent with its reported pseudo-estrogen activity, but longer-term, chronic treatments for over 6 months suppresses both ER and PR expression and signaling. In fact, washout of these chronically exposed cells for up to 1 month failed to fully reverse the transcriptional and phenotypic effects of prolonged treatments, indicating durable changes in cellular physiologic identity. RNA-seq studies found that chronic iAs drives the transition toward more basal phenotypes characterized by impaired hormone receptor signaling despite the conservation of estrogen receptor expression. Because treatments for breast cancer patients are largely designed based on the detection of hormone receptor expression, our results suggest greater scrutiny of ER+ cancers in patients exposed to iAs, because these tumors may spawn more aggressive phenotypes than unexposed ER+ tumors, in particular, basal subtypes that tend to develop therapy resistance and metastasis.

The Role of the Histone Methyltransferase Enhancer of Zeste Homolog 2 (EZH2) in the Pathobiological Mechanisms Underlying Inflammatory Bowel Disease (IBD).

Regulatory T (Treg) cells expressing the transcription factor FOXP3 play a pivotal role in maintaining immunologic self-tolerance. We and others have shown previously that EZH2 is recruited to the FOXP3 promoter and its targets in Treg cells. To further address the role for EZH2 in Treg cellular function, we have now generated mice that lack EZH2 specifically in Treg cells (EZH2Δ/ΔFOXP3+). We find that EZH2 deficiency in FOXP3+ T cells results in lethal multiorgan autoimmunity. We further demonstrate that EZH2Δ/ΔFOXP3+ T cells lack a regulatory phenotype in vitro and secrete proinflammatory cytokines. Of special interest, EZH2Δ/ΔFOXP3+ mice develop spontaneous inflammatory bowel disease. Guided by these results, we assessed the FOXP3 and EZH2 gene networks by RNA sequencing in isolated intestinal CD4+ T cells from patients with Crohn's disease. Gene network analysis demonstrates that these CD4+ T cells display a Th1/Th17-like phenotype with an enrichment of gene targets shared by FOXP3 and EZH2. Combined, these results suggest that the inflammatory milieu found in Crohn's disease could lead to or result from deregulation of FOXP3/EZH2-enforced T cell gene networks contributing to the underlying intestinal inflammation.

Transcriptional Profiling Underscores the Role of Preprocurement Allograft Metabolism and Innate Immune Status on Outcomes in Human Liver Transplantation.

Objective: The adverse effects of ischemia-reperfusion injury (IRI) remain a principal barrier to a successful outcome after lifesaving orthotopic liver transplantation (OLT). Gene expression during different phases of IRI is dynamic and modified by individual exposures, making it attractive for identifying potential therapeutic targets for improving the number of suitable organs for transplantation and patient outcomes. However, data remain limited on the functional landscape of gene expression during liver graft IRI, spanning procurement to reperfusion and recovery. Therefore, we sought to characterize transcriptomic profiles of IRI during multiple phases in human OLT. Methods: We conducted clinical data analyses, histologic evaluation, and RNA sequencing of 17 consecutive human primary OLT. We performed liver allograft biopsies at 4 time points: baseline (B, before donor cross-clamp), at the end of cold ischemia (CI), during early reperfusion (ER, after revascularization), and during late reperfusion (LR). Data were generated and then recipients grouped by post-OLT outcomes categories: immediate allograft function (IAF; n = 11) versus early allograft dysfunction (EAD; n = 6) groups. Results: We observed that CI (vs B) modified a transcriptomic landscape enriched for a metabolic and immune process. Expression levels of hallmark inflammatory response genes were higher transitioning from CI to ER and decreased from ER to LR. IAF group predominantly showed higher bile and fatty acid metabolism activity during LR compared with EAD group, while EAD group maintained more immunomodulatory activities. Throughout all time points, EAD specimens exhibited decreased metabolic activity in both bile and fatty acid pathways. Conclusions: We report transcriptomic profiles of human liver allograft IRI from prepreservation in the donor to posttransplantation in the recipient. Immunomodulatory and metabolic landscapes across ER and LR phases were different between IAF and EAD allografts. Our study also highlights marker genes for these biological processes that we plan to explore as novel therapeutic targets or surrogate markers for severe allograft injury in clinical OLT.

Verteporfin induces lipid peroxidation and ferroptosis in pancreatic cancer cells.

Pancreatic ductal adenocarcinoma (PDAC) has extremely poor prognosis, with a 5-year survival rate of approximately 11 %. Yes-associated protein (YAP) is a major downstream effector of the Hippo-YAP pathway and plays a pivotal role in regulation of cell proliferation and organ regeneration and tumorigenesis. Activation of YAP signaling has been associated with PDAC progression and drug resistance. Verteporfin (VP) is a photosensitizer used for photodynamic therapy and previous work showed that it can function as a YAP inhibitor. The efficacy of VP on human cancer are being tested in several trials. In this study, we examined the effect of VP on reactive oxygen species (ROS) and lipid peroxidation in pancreatic cancer cells, by using fluorescent molecular probes and by measuring the levels of malondialdehyde, a metabolic byproduct and marker of lipid peroxidation. We found that VP causes rapid increase of both overall ROS and lipid peroxide levels, independent of light activation. These effects were not dependent on YAP, as knockdown of YAP did not cause ROS or lipid peroxidation or enhance VP-induced ROS production. Temoporfin, another photodynamic drug, did not show similar activities. In addition, VP treatment led to loss of cell membrane integrity and reduction of viability. Notably, the activity of VP to induce lipid peroxidation was neutralized by ferroptosis inhibitors ferrostatin-1 or liproxstatin-1. VP treatment also reduced the levels of glutathione peroxidase 4 (GPX4), an enzyme that protects against lipid peroxidation. These results indicate that VP can induce lipid peroxidation and ferroptosis in the absence of light activation. Our findings reveal a novel mechanism by which VP inhibits tumor growth and provide insights into development of new therapeutic strategies for the treatment of pancreatic cancer.

Type I interferon signaling promotes radioresistance in head and neck cancer.

Despite the promise of concurrent radiotherapy (RT) and immunotherapy in head and neck cancer (HNC), multiple randomized trials of this combination have had disappointing results. To evaluate potential immunologic mechanisms of RT resistance, we compared pre-treatment HNCs that developed RT resistance to a matched cohort that achieved curative status. Gene set enrichment analysis demonstrated that a pre-treatment pro-immunogenic tumor microenvironment (TME), including type II interferon [interferon gamma (IFNγ)] and tumor necrosis factor alpha (TNFα) signaling, predicted cure while type I interferon [interferon alpha (IFNα)] enrichment was associated with an immunosuppressive TME found in tumors that went on to recur. We then used immune deconvolution of RNA sequencing datasets to evaluate immunologic cell subset enrichment. This identified M2 macrophage signaling associated with type I IFN pathway expression in RT-recurrent disease. To further dissect mechanism, we then evaluated differential gene expression between pre-treatment and RT-resistant HNCs from sampled from the same patients at the same anatomical location in the oral cavity. Here, recurrent samples exhibited upregulation of type I IFN-stimulated genes (ISGs) including members of the IFN-induced protein with tetratricopeptide repeats (IFIT) and IFN-induced transmembrane (IFITM) gene families. While several ISGs were upregulated in each recurrent cancer, IFIT2 was significantly upregulated in all recurrent tumors when compared with the matched pre-RT specimens. Based on these observations, we hypothesized sustained type I IFN signaling through ISGs, such as IFIT2, may suppress the intra-tumoral immune response thereby promoting radiation resistance.

EHMT2 Inactivation in Pancreatic Epithelial Cells Shapes the Transcriptional Landscape and Inflammation Response of the Whole Pancreas.

The Euchromatic Histone Methyl Transferase Protein 2 (EHMT2), also known as G9a, deposits transcriptionally repressive chromatin marks that play pivotal roles in the maturation and homeostasis of multiple organs. Recently, we have shown that EHMT2 inactivation alters growth and immune gene expression networks, antagonizing KRAS-mediated pancreatic cancer initiation and promotion. Here, we elucidate the essential role of EHMT2 in maintaining a transcriptional landscape that protects organs from inflammation. Comparative RNA-seq studies between normal postnatal and young adult pancreatic tissue from EHMT2 conditional knockout animals ( EHMT2 fl/fl ) targeted to the exocrine pancreatic epithelial cells ( Pdx1-Cre and P48 Cre/+ ), reveal alterations in gene expression networks in the whole organ related to injury-inflammation-repair, suggesting an increased predisposition to damage. Thus, we induced an inflammation repair response in the EHMT2 fl/fl pancreas and used a data science-based approach to integrate RNA-seq-derived pathways and networks, deconvolution digital cytology, and spatial transcriptomics. We also analyzed the tissue response to damage at the morphological, biochemical, and molecular pathology levels. The EHMT2 fl/fl pancreas displays an enhanced injury-inflammation-repair response, offering insights into fundamental molecular and cellular mechanisms involved in this process. More importantly, these data show that conditional EHMT2 inactivation in exocrine cells reprograms the local environment to recruit mesenchymal and immunological cells needed to mount an increased inflammatory response. Mechanistically, this response is an enhanced injury-inflammation-repair reaction with a small contribution of specific EHMT2-regulated transcripts. Thus, this new knowledge extends the mechanisms underlying the role of the EHMT2-mediated pathway in suppressing pancreatic cancer initiation and modulating inflammatory pancreatic diseases.

Ehmt2 inactivation in pancreatic epithelial cells shapes the transcriptional landscape and inflammation response of the whole pancreas.

Introduction: The Euchromatic Histone Methyl Transferase Protein 2 (EHMT2), also known as G9a, deposits transcriptionally repressive chromatin marks that play pivotal roles in the maturation and homeostasis of multiple organs. Recently, we have shown that Ehmt2 inactivation in the mouse pancreas alters growth and immune gene expression networks, antagonizing Kras-mediated pancreatic cancer initiation and promotion. Here, we elucidate the essential role of Ehmt2 in maintaining a transcriptional landscape that protects organs from inflammation. Methods: Comparative RNA-seq studies between normal postnatal and young adult pancreatic tissue from Ehmt2 conditional knockout animals (Ehmt2 fl/fl ) targeted to the exocrine pancreatic epithelial cells (Pdx1-Cre and P48 Cre/+ ), reveal alterations in gene expression networks in the whole organ related to injury-inflammation-repair, suggesting an increased predisposition to damage. Thus, we induced an inflammation repair response in the Ehmt2 fl/fl pancreas and used a data science-based approach to integrate RNA-seq-derived pathways and networks, deconvolution digital cytology, and spatial transcriptomics. We also analyzed the tissue response to damage at the morphological, biochemical, and molecular pathology levels. Results and discussion: The Ehmt2 fl/fl pancreas displays an enhanced injury-inflammation-repair response, offering insights into fundamental molecular and cellular mechanisms involved in this process. More importantly, these data show that conditional Ehmt2 inactivation in exocrine cells reprograms the local environment to recruit mesenchymal and immunological cells needed to mount an increased inflammatory response. Mechanistically, this response is an enhanced injury-inflammation-repair reaction with a small contribution of specific Ehmt2-regulated transcripts. Thus, this new knowledge extends the mechanisms underlying the role of the Ehmt2-mediated pathway in suppressing pancreatic cancer initiation and modulating inflammatory pancreatic diseases.

Changes in Daily Apparent Diffusion Coefficient on Fully Quantitative Magnetic Resonance Imaging Correlate With Established Genomic Pathways of Radiation Sensitivity and Reveal Novel Biologic Associations.

PURPOSE: Changes in quantitative magnetic resonance imaging (qMRI) are frequently observed during chemotherapy or radiation therapy (RT). It is hypothesized that qMRI features are reflective of underlying tissue responses. It's unknown what underlying genomic characteristics underly qMRI changes. We hypothesized that qMRI changes may correlate with DNA damage response (DDR) capacity within human tumors. Therefore, we designed the current study to correlate qMRI changes from daily RT treatment with underlying tumor transcriptomic profiles. METHODS AND MATERIALS: Study participants were prospectively enrolled (National Clinical Trial 03500081). RNA expression levels for 757 genes from pretreatment biopsies were obtained using a custom panel that included signatures of radiation sensitivity and DDR. Daily qMRI data were obtained from a 1.5 Tesla MR linear accelerator. Using these images, d-slow, d-star, perfusion, and apparent diffusion coefficient-mean values in tumors were plotted per-fraction, over time, and associated with genomic pathways. RESULTS: A total of 1022 qMRIs were obtained from 39 patients and both genomic data and qMRI data from 27 total patients. For 20 of those patients, we also generated normal tissue transcriptomic data. Radio sensitivity index values most closely associated with tissue of origin. Multiple genomic pathways including DNA repair, peroxisome, late estrogen receptor responses, KRAS signaling, and UV response were significantly associated with qMRI feature changes (P < .001). CONCLUSIONS: Genomic pathway associations across metabolic, RT sensitivity, and DDR pathways indicate common tumor biology that may correlate with qMRI changes during a course of treatment. Such data provide hypothesis-generating novel mechanistic insight into the biologic meaning of qMRI changes during treatment and enable optimal selection of imaging biomarkers for biologically MR-guided RT.

Detailed structural-functional analysis of the Krüppel-like factor 16 (KLF16) transcription factor reveals novel mechanisms for silencing Sp/KLF sites involved in metabolism and endocrinology.

Krüppel-like factor (KLF) proteins have elicited significant attention due to their emerging key role in metabolic and endocrine diseases. Here, we extend this knowledge through the biochemical characterization of KLF16, unveiling novel mechanisms regulating expression of genes involved in reproductive endocrinology. We found that KLF16 selectively binds three distinct KLF-binding sites (GC, CA, and BTE boxes). KLF16 also regulated the expression of several genes essential for metabolic and endocrine processes in sex steroid-sensitive uterine cells. Mechanistically, we determined that KLF16 possesses an activation domain that couples to histone acetyltransferase-mediated pathways, as well as a repression domain that interacts with the histone deacetylase chromatin-remodeling system via all three Sin3 isoforms, suggesting a higher level of plasticity in chromatin cofactor selection. Molecular modeling combined with molecular dynamic simulations of the Sin3a-KLF16 complex revealed important insights into how this interaction occurs at an atomic resolution level, predicting that phosphorylation of Tyr-10 may modulate KLF16 function. Phosphorylation of KLF16 was confirmed by in vivo (32)P incorporation and controlled by a Y10F site-directed mutant. Inhibition of Src-type tyrosine kinase signaling as well as the nonphosphorylatable Y10F mutation disrupted KLF16-mediated gene silencing, demonstrating that its function is regulatable rather than constitutive. Subcellular localization studies revealed that signal-induced nuclear translocation and euchromatic compartmentalization constitute an additional mechanism for regulating KLF16 function. Thus, this study lends insights on key biochemical mechanisms for regulating KLF sites involved in reproductive biology. These data also contribute to the new functional information that is applicable to understanding KLF16 and other highly related KLF proteins.

Sequence-specific recruitment of heterochromatin protein 1 via interaction with Krüppel-like factor 11, a human transcription factor involved in tumor suppression and metabolic diseases.

Heterochromatin protein 1 (HP1) proteins are "gatekeepers" of epigenetic gene silencing that is mediated by lysine 9 of histone H3 methylation (H3K9me). Current knowledge supports a paradigm whereby HP1 proteins achieve repression by binding to H3K9me marks and interacting to H3K9 histone methyltransferases (HMTs), such as SUV39H1, which methylate this residue on adjacent nucleosomes thereby compacting chromatin and silencing gene expression. However, the mechanism underlying the recruitment of this epigenetic regulator to target gene promoters remains poorly characterized. In the current study, we reveal for the first time a mechanism whereby HP1 is recruited to promoters by a well characterized Krüppel-like transcription factor (KLF), in a sequence-specific manner, to mediate complex biological phenomena. A PXVXL HP1-interacting domain identified at position 487-491 of KLF11 mediates the binding of HP1α and KLF11 in vitro and in cultured cells. KLF11 also recruits HP1α and its histone methyltransferase, SUV39H1, to promoters to limit KLF11-mediated gene activation. Indeed, a KLF11ΔHP1 mutant derepresses KLF11-regulated cancer genes, by inhibiting HP1-SUV39H1 recruitment, decreasing H3K9me3, while increasing activation-associated marks. Biologically, impairment of KLF11-mediated HP1-HMT recruitment abolishes tumor suppression, providing direct evidence that HP1-HMTs act in a sequence-specific manner to achieve this function rather than its well characterized binding to methylated chromatin without intermediary. Collectively, these studies reveal a novel role for HP1 as a cofactor in tumor suppression, expand our mechanistic understanding of a KLF associated to human disease, and outline cellular and biochemical mechanisms underlying this phenomenon, increasing the specificity of targeting HP1-HMT complexes to gene promoters.

Pepsinogen/Proton Pump Co-Expression in Barrett's Esophageal Cells Induces Cancer-Associated Changes.

EDUCATIONAL OBJECTIVE: At the conclusion of this presentation, participants should better understand the carcinogenic potential of pepsin and proton pump expression in Barrett's esophagus. OBJECTIVE: Barrett's esophagus (BE) is a well-known risk factor for esophageal adenocarcinoma (EAC). Gastric H+ /K+ ATPase proton pump and pepsin expression has been demonstrated in some cases of BE; however, the contribution of local pepsin and proton pump expression to carcinogenesis is unknown. In this study, RNA sequencing was used to examine global transcriptomic changes in a BE cell line ectopically expressing pepsinogen and/or gastric H+ /K+ ATPase proton pumps. STUDY DESIGN: In vitro translational. METHODS: BAR-T, a human BE cell line devoid of expression of pepsinogen or proton pumps, was transduced by lentivirus-encoding pepsinogen (PGA5) and/or gastric proton pump subunits (ATP4A, ATP4B). Changes relative to the parental line were assessed by RNA sequencing. RESULTS: Top canonical pathways associated with protein-coding genes differentially expressed in pepsinogen and/or proton pump expressing BAR-T cells included those involved in the tumor microenvironment and epithelial-mesenchymal transition. Top upstream regulators of coding transcripts included TGFB1 and ERBB2, which are associated with the pathogenesis and prognosis of BE and EAC. Top upstream regulators of noncoding transcripts included p300-CBP, I-BET-151, and CD93, which have previously described associations with EAC or carcinogenesis. The top associated disease of both coding and noncoding transcripts was cancer. CONCLUSIONS: These data support the carcinogenic potential of pepsin and proton pump expression in BE and reveal molecular pathways affected by their expression. Further study is warranted to investigate the role of these pathways in carcinogenesis associated with BE. LEVEL OF EVIDENCE: NA Laryngoscope, 133:59-69, 2023.

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.

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.

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.

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.

The splanchnic mesenchyme is the tissue of origin for pancreatic fibroblasts during homeostasis and tumorigenesis.

Pancreatic cancer is characterized by abundant desmoplasia, a dense stroma composed of extra-cellular and cellular components, with cancer associated fibroblasts (CAFs) being the major cellular component. However, the tissue(s) of origin for CAFs remains controversial. Here we determine the tissue origin of pancreatic CAFs through comprehensive lineage tracing studies in mice. We find that the splanchnic mesenchyme, the fetal cell layer surrounding the endoderm from which the pancreatic epithelium originates, gives rise to the majority of resident fibroblasts in the normal pancreas. In a genetic mouse model of pancreatic cancer, resident fibroblasts expand and constitute the bulk of CAFs. Single cell RNA profiling identifies gene expression signatures that are shared among the fetal splanchnic mesenchyme, adult fibroblasts and CAFs, suggesting a persistent transcriptional program underlies splanchnic lineage differentiation. Together, this study defines the phylogeny of the mesenchymal component of the pancreas and provides insights into pancreatic morphogenesis and tumorigenesis.

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.