Destabilization of the TWIST1/E12 complex dimerization following the R154P point-mutation of TWIST1: an in silico approach.
Ontology highlight
ABSTRACT: The bHLH transcription factor TWIST1 plays a key role in the embryonic development and in tumorigenesis. Some loss-of-function mutations of the TWIST1 gene have been shown to cause an autosomal dominant craniosynostosis, known as the Saethre-Chotzen syndrome (SCS). Although the functional impacts of many TWIST1 mutations have been experimentally reported, little is known on the molecular mechanisms underlying their loss-of-function. In a previous study, we highlighted the predictive value of in silico molecular dynamics (MD) simulations in deciphering the molecular function of TWIST1 residues.Here, since the substitution of the arginine 154 amino acid by a glycine residue (R154G) is responsible for the SCS phenotype and the substitution of arginine 154 by a proline experimentally decreases the dimerizing ability of TWIST1, we investigated the molecular impact of this point mutation using MD approaches. Consistently, MD simulations highlighted a clear decrease in the stability of the ?-helix during the dimerization of the mutated R154P TWIST1/E12 dimer compared to the wild-type TE complex, which was further confirmed in vitro using immunoassays.Our study demonstrates that MD simulations provide a structural explanation for the loss-of-function associated with the SCS TWIST1 mutation and provides a proof of concept of the predictive value of these MD simulations. This in silico methodology could be used to determine reliable pharmacophore sites, leading to the application of docking approaches in order to identify specific inhibitors of TWIST1 complexes.
<h4>Background</h4>The bHLH transcription factor TWIST1 plays a key role in the embryonic development and in tumorigenesis. Some loss-of-function mutations of the TWIST1 gene have been shown to cause an autosomal dominant craniosynostosis, known as the Saethre-Chotzen syndrome (SCS). Although the functional impacts of many TWIST1 mutations have been experimentally reported, little is known on the molecular mechanisms underlying their loss-of-function. In a previous study, we highlighted the pred ...[more]
Project description:The TWIST1 embryonic transcription factor displays biphasic functions during the course of carcinogenesis. It facilitates the escape of cells from oncogene-induced fail-safe programs (senescence, apoptosis) and their consequent neoplastic transformation. Additionally, it promotes the epithelial-to-mesenchymal transition and the initiation of the metastatic spread of cancer cells. Interestingly, cancer cells recurrently remain dependent on TWIST1 for their survival and/or proliferation, making TWIST1 their Achilles' heel. TWIST1 has been reported to form either homodimeric or heterodimeric complexes mainly in association with the E bHLH class I proteins. These complexes display distinct, sometimes even antagonistic, functions during development and unequal prometastatic functions in prostate cancer cells. Using a tethered dimer strategy, we successively assessed the ability of TWIST1 dimers to cooperate with an activated version of RAS in human mammary epithelial cell transformation, to provide mice with the ability to spontaneously develop breast tumors, and lastly to maintain a senescence program at a latent state in several breast cancer cell lines. We demonstrate that the TWIST1-E12 complex, unlike the homodimer, is an oncogenic form of TWIST1 in mammary epithelial cells and that efficient binding of both partners is a prerequisite for its activity. The detection of the heterodimer in human premalignant lesions by a proximity ligation assay, at a stage preceding the initiation of the metastatic cascade, is coherent with such an oncogenic function. TWIST1-E protein heterodimeric complexes may thus constitute the main active forms of TWIST1 with regard to senescence inhibition over the time course of breast tumorigenesis.
Project description:The TWIST1 bHLH transcription factor controls embryonic development and cancer processes. Although molecular and genetic analyses have provided a wealth of data on the role of bHLH transcription factors, very little is known on the molecular mechanisms underlying their binding affinity to the E-box sequence of the promoter. Here, we used an in silico model of the TWIST1/E12 (TE) heterocomplex and performed molecular dynamics (MD) simulations of its binding to specific (TE-box) and modified E-box sequences. We focused on (i) active E-box and inactive E-box sequences, on (ii) modified active E-box sequences, as well as on (iii) two box sequences with modified adjacent bases the AT- and TA-boxes. Our in silico models were supported by functional in vitro binding assays. This exploration highlighted the predominant role of protein side-chain residues, close to the heart of the complex, at anchoring the dimer to DNA sequences, and unveiled a shift towards adjacent ((-1) and (-1*)) bases and conserved bases of modified E-box sequences. In conclusion, our study provides proof of the predictive value of these MD simulations, which may contribute to the characterization of specific inhibitors by docking approaches, and their use in pharmacological therapies by blocking the tumoral TWIST1/E12 function in cancers.
Project description:Diffuse invasion into adjacent brain matter by glioblastoma (GBM) is largely responsible for their dismal prognosis. Previously, we showed that the TWIST1 (TW) bHLH transcription factor and its regulated gene periostin (POSTN) promote invasive phenotypes of GBM cells. Since TW functional effects are regulated by phosphorylation and dimerization, we investigated how phosphorylation of serine 68 in TW regulates TW dimerization, POSTN expression, and invasion in glioma cells. Compared with wild-type TW, the hypophosphorylation mutant, TW(S68A), impaired TW heterodimerization with the E12 bHLH transcription factor and cell invasion in vitro but had no effect on TW homodimerization. Overexpression of TW:E12 forced dimerization constructs (FDCs) increased glioma cell invasion and upregulated pro-invasive proteins, including POSTN, in concert with cytoskeletal reorganization. By contrast, TW:TW homodimer FDCs inhibited POSTN expression and cell invasion in vitro. Further, phosphorylation of analogous PXSP phosphorylation sites in TW:E12 FDCs (TW S68 and E12 S139) coordinately regulated POSTN and PDGFRa mRNA expression. These results suggested that TW regulates pro-invasive phenotypes in part through coordinated phosphorylation events in TW and E12 that promote heterodimer formation and regulate downstream targets. This new mechanistic understanding provides potential therapeutic strategies to inhibit TW-POSTN signaling in GBM and other cancers.
Project description:Faustoviruses are the first giant viruses of amoebae isolated on Vermamoeba vermiformis. They are distantly related to African swine fever virus, the causative agent of lethal hemorrhagic fever in domestic pigs. Structural studies have shown the presence of a double protein layer encapsidating the double-stranded DNA genome of Faustovirus E12, the prototype strain. The major capsid protein (MCP) forming the external layer has been shown to be 645-amino acid-long. Unexpectedly, its encoding sequence has been found to be scattered along a 17 kbp-large genomic region. Using RNA-seq, we studied expression of Faustovirus E12 genes at nine time points over its entire replicative cycle. Paired-end 250 bp-long read sequencing on MiSeq instrument and double-round spliced alignment enabled the identification of 26 different splice-junctions. Reads corresponding to junctions represented 2% of mapped reads and mostly matched with the predicted MCP encoding sequences. Moreover, our study enabled describing a 1,939 bp-long transcript that corresponds to the MCP, delineating 13 exons. At least two types of introns coexist in the MCP gene: group I introns that can self-splice (n = 5) and spliceosome-like introns with non-canonical splice sites (n = 7). All splice-sites were non-canonical with five types of donor/acceptor splice-sites among which AA/TG was the most frequent association.
Project description:A homozygous mutation in growth hormone 1 (GH1) was recently identified in an individual with growth failure. This mutation, c.705G>C, causes replacement of cysteine at position 53 of the 191-amino-acid sequence of 22 kDa human GH (hGH) with serine (p.C53S). This hGH molecule (hereafter referred to as GH-C53S) lacks the disulfide bond between p.Cys-53 and p.Cys-165, which is highly conserved among species. It has been reported previously that monomeric GH-C53S has reduced bioactivity compared with WT GH (GH-WT) because of its decreased ability to bind and activate the GH receptor in vitro In this study, we discovered that substitution of p.Cys-53 in hGH significantly increased formation of hGH dimers in pituitary cells. We expressed His-tagged hGH variants in the cytoplasm of genetically modified Rosetta-gami B DE3 Escherichia coli cells, facilitating high-yield production. We observed that the bioactivity of monomeric GH-C53S is 25.2% of that of GH-WT and that dimeric GH-C53S-His has no significant bioactivity in cell proliferation assays. We also found that the expression of GH-C53S in pituitary cells deviates from that of GH-WT. GH-C53S was exclusively stained in the Golgi apparatus, and no secretory granules formed for this variant, impairing its stimulated release. In summary, the unpaired Cys-165 in GH-C53S forms a disulfide bond linking two hGH molecules in pituitary cells. We conclude that the GH-C53S dimer is inactive and responsible for the growth failure in the affected individual.
Project description:Leigh syndrome is a frequent, heterogeneous pediatric presentation of mitochondrial oxidative phosphorylation (OXPHOS) disease, manifesting with psychomotor retardation and necrotizing lesions in brain deep gray matter. OXPHOS occurs at the inner mitochondrial membrane through the integrated activity of five protein complexes, of which complex V (CV) functions in a dimeric form to directly generate adenosine triphosphate (ATP). Mutations in several different structural CV subunits cause Leigh syndrome; however, dimerization defects have not been associated with human disease. We report four Leigh syndrome subjects from three unrelated Ashkenazi Jewish families harboring a homozygous splice-site mutation (c.87 + 1G>C) in a novel CV subunit disease gene, USMG5. The Ashkenazi population allele frequency is 0.57%. This mutation produces two USMG5 transcripts, wild-type and lacking exon 3. Fibroblasts from two Leigh syndrome probands had reduced wild-type USMG5 mRNA expression and undetectable protein. The mutation did not alter monomeric CV expression, but reduced both CV dimer expression and ATP synthesis rate. Rescue with wild-type USMG5 cDNA in proband fibroblasts restored USMG5 protein, increased CV dimerization and enhanced ATP production rate. These data demonstrate that a recurrent USMG5 splice-site founder mutation in the Ashkenazi Jewish population causes autosomal recessive Leigh syndrome by reduction of CV dimerization and ATP synthesis.
Project description:Rapid improvements in the detection and tracking of early-stage tumor progression aim to guide decisions regarding cancer treatments as well as predict metastatic recurrence in patients following surgery. Mathematical models may have the potential to further assist in estimating metastatic risk, particularly when paired with in vivo tumor data that faithfully represent all stages of disease progression. Herein, we describe mathematical analysis that uses data from mouse models of spontaneous metastasis developing after surgical removal of orthotopically implanted primary tumors. Both presurgical (primary tumor) growth and postsurgical (metastatic) growth were quantified using bioluminescence and were then used to generate a mathematical formalism based on general laws of the disease (i.e., dissemination and growth). The model was able to fit and predict pre/postsurgical data at the level of the individual as well as the population. Our approach also enabled retrospective analysis of clinical data describing the probability of metastatic relapse as a function of primary tumor size. In these data-based models, interindividual variability was quantified by a key parameter of intrinsic metastatic potential. Critically, our analysis identified a highly nonlinear relationship between primary tumor size and postsurgical survival, suggesting possible threshold limits for the utility of tumor size as a predictor of metastatic recurrence. These findings represent a novel use of clinically relevant models to assess the impact of surgery on metastatic potential and may guide optimal timing of treatments in neoadjuvant (presurgical) and adjuvant (postsurgical) settings to maximize patient benefit.
Project description:Cannabinoid receptor 1 (CB1) is a promising drug target involved in many physiological processes. Using atomistic molecular dynamics (MD) simulations, we examined the structural effect of F237L mutation on CB1, a mutation that has qualitatively similar effects to allosteric ligand ORG27569 binding. This mutation showed a global effect on CB1 conformations. Among the observed effects, TM6 outward movement and the conformational change of the NPxxY motif upon receptor activation by CB1 agonist CP55940 were hindered compared to wt CB1. Within the orthosteric binding site, CP55940 interactions with CB1 were altered. Our results revealed that allosteric perturbations introduced by the mutation had a global impact on receptor conformations, suggesting that the mutation site is a key region for allosteric modulation in CB1.
Project description:Autosomal dominant mutations in the gene encoding the basic helix-loop-helix transcription factor Twist1 are associated with limb and craniofacial defects in humans with Saethre-Chotzen syndrome. The molecular mechanism underlying these phenotypes is poorly understood. We show that ectopic expression of the related basic helix-loop-helix factor Hand2 phenocopies Twist1 loss of function in the limb and that the two factors have a gene dosage-dependent antagonistic interaction. Dimerization partner choice by Twist1 and Hand2 can be modulated by protein kinase A- and protein phosphatase 2A-regulated phosphorylation of conserved helix I residues. Notably, multiple Twist1 mutations associated with Saethre-Chotzen syndrome alter protein kinase A-mediated phosphorylation of Twist1, suggesting that misregulation of Twist1 dimerization through either stoichiometric or post-translational mechanisms underlies phenotypes of individuals with Saethre-Chotzen syndrome.
Project description:Mutations in the juxtamembrane and tyrosine kinase domains of the KIT receptor have been implicated in several cancers and are known to promote tumorigenesis. However, the pathophysiological manifestations of mutations in the extracellular domain remain unknown. In this study, we examined the impact of a mutation in the extracellular domain of KIT on mast cell tumorigenesis. A KIT mutant with an Asn508Ile variation (N508I) in the extracellular domain derived from a canine mast cell tumor was introduced into IC-2 cells. The IC-2(N508I) cells proliferated in a cytokine-independent manner and showed KIT auto-phosphorylation. Subcutaneous injection of IC-2(N508I) cells into the dorsal area of immunodeficient BALB/c-nu/nu mice resulted in the formation of solid tumors, but tumor progression was abrogated by treatment with a tyrosine kinase inhibitor (STI571). In addition, the N508I mutant KIT protein dimerized in the absence of the natural ligand, stem cell factor. Structure modeling indicates that the increased hydrophobicity of the mutant led to the stabilization of KIT dimers. These results suggest that this extracellular domain mutation confers a ligand-independent tumorigenic phenotype to mast cells by KIT auto-dimerization that is STI571-sensitive. This is the first report demonstrating the tumorigenic potential of a mutation in the extracellular domain of KIT.