Structural basis for DNA recognition by FOXC2

Unravelling the therapeutic potential of forkhead box proteins in breast cancer: An update (Review)

Breast cancer, a prominent cause of mortality among women, develops from abnormal growth of breast tissue, thereby rendering it one of the most commonly detected cancers in the female population. Although numerous treatment strategies are available for breast cancer, discordance in terms of effective treatment and response still exists. Recently, the potential of signaling pathways and transcription factors has gained substantial attention in the cancer community; therefore, understanding their role will assist researchers in comprehending the onset and advancement of breast cancer. Forkhead box (FOX) proteins, which are important transcription factors, are considered crucial regulators of various cellular activities, including cell division and proliferation. The present study explored several subclasses of FOX proteins and their possible role in breast carcinogenesis, followed by the interaction between microRNA (miRNA) and FOX proteins. This interaction is implicated in promoting cell infiltration into surrounding tissues, ultimately leading to metastasis. The various roles that FOX proteins play in breast cancer development, their intricate relationships with miRNA, and their involvement in therapeutic resistance highlight the complexity of breast cancer dynamics. Therefore, recognizing the progress and challenges in current treatments is crucial because, despite advancements, persistent disparities in treatment effectiveness underscore the need for ongoing research, with future studies emphasizing the necessity for targeted strategies that account for the multifaceted aspects of breast cancer.

Proximal telomeric decompaction due to telomere shortening drives FOXC1-dependent myocardial senescence

Telomeres, TTAGGGn DNA repeat sequences located at the ends of eukaryotic chromosomes, play a pivotal role in aging and are targets of DNA damage response. Although we and others have demonstrated presence of short telomeres in genetic cardiomyopathic and heart failure cardiomyocytes, little is known about the role of telomere lengths in cardiomyocyte. Here, we demonstrate that in heart failure patient cardiomyocytes, telomeres are shortened compared to healthy controls. We generated isogenic human induced pluripotent stem cell derived cardiomyocytes (hiPSC-CMs) with short telomeres (sTL-CMs) and normal telomeres (nTL-CMs) as model. Compared to nTL-CMs, short telomeres result in cardiac dysfunction and expression of senescent markers. Using Hi-C and RNASeq, we observe that short telomeres induced TAD insulation decrease near telomeric ends and this correlated with a transcription upregulation in sTL-CMs. FOXC1, a key transcription factor involved in early cardiogenesis, was upregulated in sTL-CMs and its protein levels were negatively correlated with telomere lengths in heart failure patients. Overexpression of FOXC1 induced hiPSC-CM aging, mitochondrial and contractile dysfunction; knockdown of FOXC1 rescued these phenotypes. Overall, the work presented demonstrate that increased chromatin accessibility due to telomere shortening resulted in the induction of FOXC1-dependent expression network responsible for contractile dysfunction and myocardial senescence.

FOXC1 regulates endothelial CD98 (LAT1/4F2hc) expression in retinal angiogenesis and blood-retina barrier formation

Angiogenesis, the growth of new blood vessels from pre-existing vasculature, is essential for the development of new organ systems, but transcriptional control of angiogenesis remains incompletely understood. Here we show that FOXC1 is essential for retinal angiogenesis. Endothelial cell (EC)-specific loss of Foxc1 impairs retinal vascular growth and expression of Slc3a2 and Slc7a5, which encode the heterodimeric CD98 (LAT1/4F2hc) amino acid transporter and regulate the intracellular transport of essential amino acids and activation of the mammalian target of rapamycin (mTOR). EC-Foxc1 deficiency diminishes mTOR activity, while administration of the mTOR agonist MHY-1485 rescues perturbed retinal angiogenesis. EC-Foxc1 expression is required for retinal revascularization and resolution of neovascular tufts in a model of oxygen-induced retinopathy. Foxc1 is also indispensable for pericytes, a critical component of the blood-retina barrier during retinal angiogenesis. Our findings establish FOXC1 as a crucial regulator of retinal vessels and identify therapeutic targets for treating retinal vascular disease.

Structural characterization of the DNA binding mechanism of retinoic acid-related orphan receptor gamma

Retinoic acid-related orphan receptor gamma (RORγ) plays critical roles in regulating various biological processes and has been linked to immunodeficiency disorders and cancers. DNA recognition is essential for RORγ to exert its functions. However, the underlying mechanism of the DNA binding by RORγ remains unclear. In this study, we present the crystal structure of the complex of RORγ1 DNA-binding domain (RORγ1-DBD)/direct repeat DNA element DR2 at 2.3 Å resolution. We demonstrate that RORγ1-DBD binds the DR2 motif as a homodimer, with the C-terminal extension (CTE) region of RORγ1-DBD contributing to the DNA recognition and the formation of dimeric interface. Further studies reveal that REV-ERB-DBD and RXR-DBD, also bind the DR2 site as a homodimer, while NR4A2-DBD binds DR2 as a monomer. Our research uncovers a binding mechanism of RORγ1 to the DR2 site and provides insights into the biological functions of RORγ1 and the broader RORs subfamily.

Structural insights into the HDAC4-MEF2A-DNA complex and its implication in long-range transcriptional regulation

Class IIa Histone deacetylases (HDACs), including HDAC4, 5, 7 and 9, play key roles in multiple important developmental and differentiation processes. Recent studies have shown that class IIa HDACs exert their transcriptional repressive function by interacting with tissue-specific transcription factors, such as members of the myocyte enhancer factor 2 (MEF2) family of transcription factors. However, the molecular mechanism is not well understood. In this study, we determined the crystal structure of an HDAC4-MEF2A-DNA complex. This complex adopts a dumbbell-shaped overall architecture, with a 2:4:2 stoichiometry of HDAC4, MEF2A and DNA molecules. In the complex, two HDAC4 molecules form a dimer through the interaction of their glutamine-rich domain (GRD) to form the stem of the 'dumbbell'; while two MEF2A dimers and their cognate DNA molecules are bridged by the HDAC4 dimer. Our structural observations were then validated using biochemical and mutagenesis assays. Further cell-based luciferase reporter gene assays revealed that the dimerization of HDAC4 is crucial in its ability to repress the transcriptional activities of MEF2 proteins. Taken together, our findings not only provide the structural basis for the assembly of the HDAC4-MEF2A-DNA complex but also shed light on the molecular mechanism of HDAC4-mediated long-range gene regulation.

Insights into the molecular roles of FOXR2 in the pathology of primary pediatric brain tumors

Forkhead box gene R2 (FOXR2) belongs to the family of FOX genes which codes for highly conserved transcription factors (TFs) with critical roles in biological processes ranging from development to organogenesis to metabolic and immune regulation to cellular homeostasis. A number of FOX genes are associated with cancer development and progression and poor prognosis. A growing body of evidence suggests that FOXR2 is an oncogene. Studies suggested important roles for FOXR2 in cancer cell growth, metastasis, and drug resistance. Recent studies showed that FOXR2 is overexpressed by a subset of newly identified entities of embryonal tumors. This review discusses the role(s) FOXR2 plays in the pathology of pediatric brain cancers and its potential as a therapeutic target.

Homology modeling of Forkhead box protein C2: identification of potential inhibitors using ligand and structure-based virtual screening

Overexpression of Forkhead box protein C2 (FOXC2) has been associated with different types of carcinomas. FOXC2 plays an important role in the initiation and maintenance of the epithelial-mesenchymal transition (EMT) process, which is essential for the development of higher-grade tumors with an enhanced ability for metastasis. Thus, FOXC2 has become a therapeutic target for the development of anticancer drugs. MC-1-F2, the only identified experimental inhibitor of FOXC2, interacts with the full length of FOXC2. However, only the DNA-binding domain (DBD) of FOXC2 has resolved crystal structure. In this work, a three-dimensional (3D) structure of the full-length FOXC2 using homology modeling was developed and used for structure-based drug design (SBDD). The quality of this 3D model of the full-length FOXC2 was evaluated using MolProbity, ERRAT, and ProSA modules. Molecular dynamics (MD) simulation was also carried out to verify its stability. Ligand-based drug design (LBDD) was carried out to identify similar analogues for MC-1-F2 against 15 million compounds from ChEMBL and ZINC databases. 792 molecules were retrieved from this similarity search. De novo SBDD was performed against the full-length 3D structure of FOXC2 through homology modeling to identify novel inhibitors. The combination of LBDD and SBDD helped in gaining a better insight into the binding of MC-1-F2 and its analogues against the full length of the FOXC2. The binding free energy of the top hits was further investigated using MD simulations and MM/GBSA calculations to result in eight promising hits as lead compounds targeting FOXC2.

Structures of p53/BCL-2 complex suggest a mechanism for p53 to antagonize BCL-2 activity

Mitochondrial apoptosis is strictly controlled by BCL-2 family proteins through a subtle network of protein interactions. The tumor suppressor protein p53 triggers transcription-independent apoptosis through direct interactions with BCL-2 family proteins, but the molecular mechanism is not well understood. In this study, we present three crystal structures of p53-DBD in complex with the anti-apoptotic protein BCL-2 at resolutions of 2.3-2.7 Å. The structures show that two loops of p53-DBD penetrate directly into the BH3-binding pocket of BCL-2. Structure-based mutations at the interface impair the p53/BCL-2 interaction. Specifically, the binding sites for p53 and the pro-apoptotic protein Bax in the BCL-2 pocket are mostly identical. In addition, formation of the p53/BCL-2 complex is negatively correlated with the formation of BCL-2 complexes with pro-apoptotic BCL-2 family members. Defects in the p53/BCL-2 interaction attenuate p53-mediated cell apoptosis. Overall, our study provides a structural basis for the interaction between p53 and BCL-2, and suggests a molecular mechanism by which p53 regulates transcription-independent apoptosis by antagonizing the interaction of BCL-2 with pro-apoptotic BCL-2 family members.

Endothelial FOXC1 and FOXC2 promote intestinal regeneration after ischemia-reperfusion injury

Intestinal ischemia underlies several clinical conditions and can result in the loss of the intestinal mucosal barrier. Ischemia-induced damage to the intestinal epithelium is repaired by stimulation of intestinal stem cells (ISCs), and paracrine signaling from the vascular niche regulates intestinal regeneration. Here, we identify FOXC1 and FOXC2 as essential regulators of paracrine signaling in intestinal regeneration after ischemia-reperfusion (I/R) injury. Vascular endothelial cell (EC)- and lymphatic EC (LEC)-specific deletions of Foxc1, Foxc2, or both in mice worsen I/R-induced intestinal damage by causing defects in vascular regrowth, expression of chemokine CXCL12 and Wnt activator R-spondin 3 (RSPO3) in blood ECs (BECs) and LECs, respectively, and activation of Wnt signaling in ISCs. Both FOXC1 and FOXC2 directly bind to regulatory elements of the CXCL12 and RSPO3 loci in BECs and LECs, respectively. Treatment with CXCL12 and RSPO3 rescues the I/R-induced intestinal damage in EC- and LEC-Foxc mutant mice, respectively. This study provides evidence that FOXC1 and FOXC2 are required for intestinal regeneration by stimulating paracrine CXCL12 and Wnt signaling.

Structural basis of the farnesoid X receptor/retinoid X receptor heterodimer on inverted repeat DNA

Farnesoid X receptor (FXR) is a ligand-activated transcription factor known as bile acid receptor (BAR). FXR plays critical roles in various biological processes, including metabolism, immune inflammation, liver regeneration and liver carcinogenesis. FXR forms a heterodimer with the retinoid X receptor (RXR) and binds to diverse FXR response elements (FXREs) to exert its various biological functions. However, the mechanism by which the FXR/RXR heterodimer binds the DNA elements remains unclear. In this study, we aimed to use structural, biochemical and bioinformatics analyses to study the mechanism of FXR binding to the typical FXRE, such as the IR1 site, and the heterodimer interactions in the FXR-DBD/RXR-DBD complex. Further biochemical assays showed that RAR, THR and NR4A2 do not form heterodimers with RXR when bound to the IR1 sites, which indicates that IR1 may be a unique binding site for the FXR/RXR heterodimer. Our studies may provide a further understanding of the dimerization specificity of nuclear receptors.

FOXI3 pathogenic variants cause one form of craniofacial microsomia

Craniofacial microsomia (CFM; also known as Goldenhar syndrome), is a craniofacial developmental disorder of variable expressivity and severity with a recognizable set of abnormalities. These birth defects are associated with structures derived from the first and second pharyngeal arches, can occur unilaterally and include ear dysplasia, microtia, preauricular tags and pits, facial asymmetry and other malformations. The inheritance pattern is controversial, and the molecular etiology of this syndrome is largely unknown. A total of 670 patients belonging to unrelated pedigrees with European and Chinese ancestry with CFM, are investigated. We identify 18 likely pathogenic variants in 21 probands (3.1%) in FOXI3. Biochemical experiments on transcriptional activity and subcellular localization of the likely pathogenic FOXI3 variants, and knock-in mouse studies strongly support the involvement of FOXI3 in CFM. Our findings indicate autosomal dominant inheritance with reduced penetrance, and/or autosomal recessive inheritance. The phenotypic expression of the FOXI3 variants is variable. The penetrance of the likely pathogenic variants in the seemingly dominant form is reduced, since a considerable number of such variants in affected individuals were inherited from non-affected parents. Here we provide suggestive evidence that common variation in the FOXI3 allele in trans with the pathogenic variant could modify the phenotypic severity and accounts for the incomplete penetrance.

Molecular basis for DNA recognition by the maternal pioneer transcription factor FoxH1

Forkhead box H1 (FoxH1) is an essential maternal pioneer factor during embryonic development that binds to specific GG/GT-containing DNA target sequences. Here we have determined high-resolution structures of three FoxH1 proteins (from human, frog and fish species) and four DNAs to clarify the way in which FoxH1 binds to these sites. We found that the protein-DNA interactions extend to both the minor and major DNA grooves and are thus almost twice as extensive as those of other FOX family members. Moreover, we identified two specific amino acid changes in FoxH1 that allowed the recognition of GG/GT motifs. Consistent with the pioneer factor activity of FoxH1, we found that its affinity for nucleosomal DNA is even higher than for linear DNA fragments. The structures reported herein illustrate how FoxH1 binding to distinct DNA sites provides specificity and avoids cross-regulation by other FOX proteins that also operate during the maternal-zygotic transition and select canonical forkhead sites.

Structure-based virtual screening identified novel FOXM1 inhibitors as the lead compounds for ovarian cancer

Ovarian cancer (OC) is a gynecological tumor with possibly the worst prognosis, its 5-year survival rate being only 47.4%. The first line of therapy prescribed is chemotherapy consisting of platinum and paclitaxel. The primary reason for treatment failure is drug resistance. FOXM1 protein has been found to be closely associated with drug resistance, and inhibition of FOXM1 expression sensitizes cisplatin-resistant ovarian cancer cells. Combining existing first-line chemotherapy drugs with FOXM1 prolongs the overall survival of patients, therefore, FOXM1 is considered a potential therapeutic target in ovarian cancer. Previous research conducted by our team revealed a highly credible conformation of FOXM1 which enables binding by small molecules. Based on this conformation, the current study conducted virtual screening to determine a new structural skeleton for FOXM1 inhibitors which would enhance their medicinal properties. DZY-4 showed the highest affinity towards FOXM1, and its inhibitory effect on proliferation and migration of ovarian cancer at the cellular level was better than or equal to that of cisplatin, while its efficacy was equivalent to that of cisplatin in a nude mouse model. In this study, the anti-tumor effect of DZY-4 is reported for the first time. DZY-4 shows potential as a drug that can be used for ovarian cancer treatment, as well as a drug lead for future research.

FOXL2 and FOXA1 cooperatively assemble on the TP53 promoter in alternative dimer configurations

Although both the p53 and forkhead box (FOX) family proteins are key transcription factors associated with cancer progression, their direct relationship is unknown. Here, we found that FOX family proteins bind to the non-canonical homotypic cluster of the p53 promoter region (TP53). Analysis of crystal structures of FOX proteins (FOXL2 and FOXA1) bound to the p53 homotypic cluster indicated that they interact with a 2:1 stoichiometry accommodated by FOX-induced DNA allostery. In particular, FOX proteins exhibited distinct dimerization patterns in recognition of the same p53-DNA; dimer formation of FOXA1 involved protein–protein interaction, but FOXL2 did not. Biochemical and biological functional analyses confirmed the cooperative binding of FOX proteins to the TP53 promoter for the transcriptional activation of TP53. In addition, up-regulation of TP53 was necessary for FOX proteins to exhibit anti-proliferative activity in cancer cells. These analyses reveal the presence of a discrete characteristic within FOX family proteins in which FOX proteins regulate the transcription activity of the p53 tumor suppressor via cooperative binding to the TP53 promoter in alternative dimer configurations.

Genetic landscape of FOXC2 mutations in lymphedema-distichiasis syndrome: Different mechanism of pathogenicity for mutations in different domains

Lymphedema-dissociated syndrome (LDS), of which the pathogenesis is not fully understood, afflicts many patients. In this study, we investigated the effect of FOXC2 gene loss-of-function on the development of LDS disease.Two Han Chinese families with LDS were recruited in this study, pathogenic mutations were identified by Sanger sequencing. Reverse-transcription PCR, subcellular localization, dual fluorescein enzymes, and other in vitro experiments were used to study the functional effects of eight FOXC2 mutations. Two pathogenic FOXC2 duplication mutations (c.930_936dup and c.931-937dup) were identified in the two families. Both mutations caused uneven distribution in the nucleus and a chromatin contraction phenotype, weakening the DNA binding activity and transcription activity. We then performed functional analysis on six additional mutations in different domains of FOXC2 that were reported to cause LDS. We found mutations located in the forkhead domain and central region dramatically reduced the transactivation ability, while mutations in activation domain-2 enhanced this ability. All 8 mutations down-regulated the transcription of ANGPT2 and affected the activity of the ERK-RAS pathway, which may cause abnormal formation of lymphatic vessels. Our findings also showed that all 8 mutations decreased the ability of interaction between FOXC2 and the Wnt4 promoter, suggesting mutations in FOXC2 may also affect the Wnt4-Frizzled-RYK signaling pathway, leading the abnormal differentiation of the meibomian glands into hair follicle cells during the embryonic period and causing distichiasis. This study expanded and revealed the potential pathogenesis mechanism.

EMDLP: Ensemble multiscale deep learning model for RNA methylation site prediction

BACKGROUND

Recent research recommends that epi-transcriptome regulation through post-transcriptional RNA modifications is essential for all sorts of RNA. Exact identification of RNA modification is vital for understanding their purposes and regulatory mechanisms. However, traditional experimental methods of identifying RNA modification sites are relatively complicated, time-consuming, and laborious. Machine learning approaches have been applied in the procedures of RNA sequence features extraction and classification in a computational way, which may supplement experimental approaches more efficiently. Recently, convolutional neural network (CNN) and long short-term memory (LSTM) have been demonstrated achievements in modification site prediction on account of their powerful functions in representation learning. However, CNN can learn the local response from the spatial data but cannot learn sequential correlations. And LSTM is specialized for sequential modeling and can access both the contextual representation but lacks spatial data extraction compared with CNN. There is strong motivation to construct a prediction framework using natural language processing (NLP), deep learning (DL) for these reasons.

RESULTS

This study presents an ensemble multiscale deep learning predictor (EMDLP) to identify RNA methylation sites in an NLP and DL way. It organically combines the dilated convolution and Bidirectional LSTM (BiLSTM), which helps to take better advantage of the local and global information for site prediction. The first step of EMDLP is to represent the RNA sequences in an NLP way. Thus, three encodings, e.g., RNA word embedding, One-hot encoding, and RGloVe, which is an improved learning method of word vector representation based on GloVe, are adopted to decipher sites from the viewpoints of the local and global information. Then, a dilated convolutional Bidirectional LSTM network (DCB) model is constructed with the dilated convolutional neural network (DCNN) followed by BiLSTM to extract potential contributing features for methylation site prediction. Finally, these three encoding methods are integrated by a soft vote to obtain better predictive performance. Experiment results on m¹A and m⁶A reveal that the area under the receiver operating characteristic(AUROC) of EMDLP obtains respectively 95.56%, 85.24%, and outperforms the state-of-the-art models. To maximize user convenience, a user-friendly webserver for EMDLP was publicly available at http://www.labiip.net/EMDLP/index.php ( http://47.104.130.81/EMDLP/index.php ).

CONCLUSIONS

We developed a predictor for m¹A and m⁶A methylation sites.

Metastasis and MAPK Pathways

Cancer is a leading cause of death worldwide. In many cases, the treatment of the disease is limited due to the metastasis of cells to distant locations of the body through the blood and lymphatic drainage. Most of the anticancer therapeutic options focus mainly on the inhibition of tumor cell growth or the induction of cell death, and do not consider the molecular basis of metastasis. The aim of this work is to provide a comprehensive review focusing on cancer metastasis and the mitogen-activated protein kinase (MAPK) pathway (ERK/JNK/P38 signaling) as a crucial modulator of this process.

Motif models proposing independent and interdependent impacts of nucleotides are related to high and low affinity transcription factor binding sites in Arabidopsis

Position weight matrix (PWM) is the traditional motif model representing the transcription factor (TF) binding sites. It proposes that the positions contribute independently to TFs binding affinity, although this hypothesis does not fit the data perfectly. This explains why PWM hits are missing in a substantial fraction of ChIP-seq peaks. To study various modes of the direct binding of plant TFs, we compiled the benchmark collection of 111 ChIP-seq datasets for Arabidopsis thaliana, and applied the traditional PWM, and two alternative motif models BaMM and SiteGA, proposing the dependencies of the positions. The variation in the stringency of the recognition thresholds for the models proposed that the hits of PWM, BaMM, and SiteGA models are associated with the sites of high/medium, any, and low affinity, respectively. At the medium recognition threshold, about 60% of ChIP-seq peaks contain PWM hits consisting of conserved core consensuses, while BaMM and SiteGA provide hits for an additional 15% of peaks in which a weaker core consensus is compensated through intra-motif dependencies. The presence/absence of these dependencies in the motifs of alternative/traditional models was confirmed by the dependency logo DepLogo visualizing the position-wise partitioning of the alignments of predicted sites. We exemplify the detailed analysis of ChIP-seq profiles for plant TFs CCA1, MYC2, and SEP3. Gene ontology (GO) enrichment analysis revealed that among the three motif models, the SiteGA had the highest portions of genes with the significantly enriched GO terms among all predicted genes. We showed that both alternative motif models provide for traditional PWM greater extensions in predicted sites for TFs MYC2/SEP3 with condition/tissue specific functions, compared to those for TF CCA1 with housekeeping functions. Overall, the combined application of standard and alternative motif models is beneficial to detect various modes of the direct TF-DNA interactions in the maximal portion of ChIP-seq loci.

Mechanistic Insights into the Preference for Tandem Binding Sites in DNA Recognition by FOXM1

FOXM1 is an essential proliferation-associated transcription factor that controls the activation of a number of cell cycle regulatory genes. Unlike other forkhead box (FOX) transcription factors, FOXM1 has been shown to prefer binding tandem regulatory DNA sites. However, the underlying reason for such preference is not clear. Here, we showed that the tandem DNA motif, named DIV2, is widely distributed in the promoter region of FOXM1 target genes. The binding of FOXM1 on the DIV2 site differs dramatically from other sites, which is in a highly cooperative fashion, with a much enhanced thermal stability and can be clearly detected by EMSA. The crystal structure of FOXM1 in complex with the DIV2 DNA reveals that the cooperative binding is likely to be driven by intermolecular protein-protein interactions (PPIs). Further half-site spacer insertion assays showed that FOXM1 can bind another site, DIV0, in a similar manner to the DIV2 site. Given the high occurrence of the DIV2 and DIV0 sites in FOXM1 target genes, our results suggest that FOXM1 prefers tandem DNA sites to enable cooperative DNA recognition, and such binding characteristics may further confer its specificity during transcriptional regulation.

Toward a mechanistic understanding of DNA binding by forkhead transcription factors and its perturbation by pathogenic mutations

Forkhead box (FOX) proteins are an evolutionarily conserved family of transcription factors that play numerous regulatory roles in eukaryotes during developmental and adult life. Dysfunction of FOX proteins has been implicated in a variety of human diseases, including cancer, neurodevelopment disorders and genetic diseases. The FOX family members share a highly conserved DNA-binding domain (DBD), which is essential for DNA recognition, binding and function. Since the first FOX structure was resolved in 1993, >30 FOX structures have been reported to date. It is clear now that the structure and DNA recognition mechanisms vary among FOX members; however, a systematic review on this aspect is lacking. In this manuscript, we present an overview of the mechanisms by which FOX transcription factors bind DNA, including protein structures, DNA binding properties and disease-causing mutations. This review should enable a better understanding of FOX family transcription factors for basic researchers and clinicians.

Human FoxP Transcription Factors as Tractable Models of the Evolution and Functional Outcomes of Three-Dimensional Domain Swapping

The association of two or more proteins to adopt a quaternary complex is one of the most widespread mechanisms by which protein function is modulated. In this scenario, three-dimensional domain swapping (3D-DS) constitutes one plausible pathway for the evolution of protein oligomerization that exploits readily available intramolecular contacts to be established in an intermolecular fashion. However, analysis of the oligomerization kinetics and thermodynamics of most extant 3D-DS proteins shows its dependence on protein unfolding, obscuring the elucidation of the emergence of 3D-DS during evolution, its occurrence under physiological conditions, and its biological relevance. Here, we describe the human FoxP subfamily of transcription factors as a feasible model to study the evolution of 3D-DS, due to their significantly faster dissociation and dimerization kinetics and lower dissociation constants in comparison to most 3D-DS models. Through the biophysical and functional characterization of FoxP proteins, relevant structural aspects highlighting the evolutionary adaptations of these proteins to enable efficient 3D-DS have been ascertained. Most biophysical studies on FoxP suggest that the dynamics of the polypeptide chain are crucial to decrease the energy barrier of 3D-DS, enabling its fast oligomerization under physiological conditions. Moreover, comparison of biophysical parameters between human FoxP proteins in the context of their minute sequence differences suggests differential evolutionary strategies to favor homoassociation and presages the possibility of heteroassociations, with direct impacts in their gene regulation function.

Application of alternative de novo motif recognition models for analysis of structural heterogeneity of transcription factor binding sites: a case study of FOXA2 binding sites

The most popular model for the search of ChIP-seq data for transcription factor binding sites (TFBS) is the positional weight matrix (PWM). However, this model does not take into account dependencies between nucleotide occurrences in different site positions. Currently, two recently proposed models, BaMM and InMoDe, can do as much. However, application of these models was usually limited only to comparing their recognition accuracies with that of PWMs, while none of the analyses of the co-prediction and relative positioning of hits of different models in peaks has yet been performed. To close this gap, we propose the pipeline called MultiDeNA. This pipeline includes stages of model training, assessing their recognition accuracy, scanning ChIP-seq peaks and their classification based on scan results. We applied our pipeline to 22 ChIP-seq datasets of TF FOXA2 and considered PWM, dinucleotide PWM (diPWM), BaMM and InMoDe models. The combination of these four models allowed a significant increase in the fraction of recognized peaks compared to that for the sole PWM model: the increase was 26.3 %. The BaMM model provided the main contribution to the recognition of sites. Although the major fraction of predicted peaks contained TFBS of different models with coincided positions, the medians of the fraction of peaks containing the predictions of sole models were 1.08, 0.49, 4.15 and 1.73 % for PWM, diPWM, BaMM and InMoDe, respectively. Thus, FOXA2 BSs were not fully described by only a sole model, which indicates theirs heterogeneity. We assume that the BaMM model is the most successful in describing the structure of the FOXA2 BS in ChIP-seq datasets under study.

Comprehensive phenotypic and functional analysis of dominant and recessive FOXE3 alleles in ocular developmental disorders

The forkhead transcription factor FOXE3 is critical for vertebrate eye development. Recessive and dominant variants cause human ocular disease but the full range of phenotypes and mechanisms of action for the two classes of variants are unknown. We identified FOXE3 variants in individuals with congenital eye malformations and carried out in vitro functional analysis on selected alleles. Sixteen new recessive and dominant families, including six novel variants, were identified. Analysis of new and previously reported genetic and clinical data demonstrated a broad phenotypic range with an overlap between recessive and dominant disease. Most families with recessive alleles, composed of truncating and forkhead-domain missense variants, had severe corneal opacity (90%; sclerocornea in 47%), aphakia (83%) and microphthalmia (80%), but some had milder features including isolated cataract. The phenotype was most variable for recessive missense variants, suggesting that the functional consequences may be highly dependent on the type of amino acid substitution and its position. When assessed, aniridia or iris hypoplasia were noted in 89% and optic nerve anomalies in 60% of recessive cases, indicating that these defects are also common and may be underrecognized. In dominant pedigrees, caused by extension variants, normal eye size (96%), cataracts (99%) and variable anterior segment anomalies were seen in most, but some individuals had microphthalmia, aphakia or sclerocornea, more typical of recessive disease. Functional studies identified variable effects on the protein stability, DNA binding, nuclear localization and transcriptional activity for recessive FOXE3 variants, whereas dominant alleles showed severe impairment in all areas and dominant-negative characteristics.

FOXC2 Autoregulates Its Expression in the Pulmonary Endothelium After Endotoxin Stimulation in a Histone Acetylation-Dependent Manner

The innate immune response of pulmonary endothelial cells (EC) to lipopolysaccharide (LPS) induces Forkhead box protein C2 (FOXC2) activation through Toll Like Receptor 4 (TLR4). The mechanisms by which FOXC2 expression is regulated in lung EC under LPS stimulation remain unclear. We postulated that FOXC2 regulates its own expression in sepsis, and its transcriptional autoregulation directs lymphatic EC cell-fate decision. Bioinformatic analysis identified potential FOXC2 binding sites in the FOXC2 promoter. In human lung EC, we verified using chromatin immunoprecipitation (ChIP) and luciferase assays that FOXC2 bound to its own promoter and stimulated its expression after LPS stimulation. Chemical inhibition of histone acetylation by garcinol repressed LPS-induced histone acetylation in the FOXC2 promoter region, and disrupted LPS-mediated FOXC2 binding and transcriptional activation. CRISPR/dCas9/gRNA directed against FOXC2-binding-element (FBE) suppressed LPS-stimulated FOXC2 binding and autoregulation by blocking FBEs in the FOXC2 promoter, and repressed expression of lymphatic EC markers. In a neonatal mouse model of sterile sepsis, LPS-induced FOXC2 binding to FBE and FOXC2 expression in lung EC was attenuated with garcinol treatment. These data reveal a new mechanism of LPS-induced histone acetylation-dependent FOXC2 autoregulation.

Imbalance between Expression of FOXC2 and Its lncRNA in Lymphedema-Distichiasis Caused by Frameshift Mutations

Forkhead-box C2 (FOXC2) is a transcription factor involved in lymphatic system development. FOXC2 mutations cause Lymphedema-distichiasis syndrome (LD). Recently, a natural antisense was identified, called lncRNA FOXC2-AS1, which increases FOXC2 mRNA stability. No studies have evaluated FOXC2 and FOXC2-AS1 blood expression in LD and healthy subjects. Here, we show that FOXC2 and FOXC-AS1 expression levels were similar in both controls and patients, and a significantly higher amount of both RNAs was observed in females. A positive correlation between FOXC2 and FOXC2-AS1 expression was found in both controls and patients, excluding those with frameshift mutations. In these patients, the FOXC2-AS1/FOXC2 ratio was about 1:1, while it was higher in controls and patients carrying other types of mutations. The overexpression or silencing of FOXC2-AS1 determined a significant increase or reduction in FOXC2 wild-type and frameshift mutant proteins, respectively. Moreover, confocal and bioinformatic analysis revealed that these variations caused the formation of nuclear proteins aggregates also involving DNA. In conclusion, patients with frameshift mutations presented lower values of the FOXC2-AS1/FOXC2 ratio, due to a decrease in FOXC2-AS1 expression. The imbalance between FOXC2 mRNA and its lncRNA could represent a molecular mechanism to reduce the amount of FOXC2 misfolded proteins, protecting cells from damage.

Structural insight into the molecular mechanism of p53-mediated mitochondrial apoptosis

The tumor suppressor p53 is mutated in approximately half of all human cancers. p53 can induce apoptosis through mitochondrial membrane permeabilization by interacting with and antagonizing the anti-apoptotic proteins BCL-xL and BCL-2. However, the mechanisms by which p53 induces mitochondrial apoptosis remain elusive. Here, we report a 2.5 Å crystal structure of human p53/BCL-xL complex. In this structure, two p53 molecules interact as a homodimer, and bind one BCL-xL molecule to form a ternary complex with a 2:1 stoichiometry. Mutations at the p53 dimer interface or p53/BCL-xL interface disrupt p53/BCL-xL interaction and p53-mediated apoptosis. Overall, our current findings of the bona fide structure of p53/BCL-xL complex reveal the molecular basis of the interaction between p53 and BCL-xL, and provide insight into p53-mediated mitochondrial apoptosis.

The structure of human tumor suppressor p53 in complex with the antiapoptotic protein BCL-xL reveals the basis of the p53–BCL-xL interaction and provides insight into the mechanisms of p53-mediated mitochondrial apoptosis.

Mechanism of forkhead transcription factors binding to a novel palindromic DNA site

Forkhead transcription factors bind a canonical consensus DNA motif, RYAAAYA (R = A/G, Y = C/T), as a monomer. However, the molecular mechanisms by which forkhead transcription factors bind DNA as a dimer are not well understood. In this study, we show that FOXO1 recognizes a palindromic DNA element DIV2, and mediates transcriptional regulation. The crystal structure of FOXO1/DIV2 reveals that the FOXO1 DNA binding domain (DBD) binds the DIV2 site as a homodimer. The wing1 region of FOXO1 mediates the dimerization, which enhances FOXO1 DNA binding affinity and complex stability. Further biochemical assays show that FOXO3, FOXM1 and FOXI1 also bind the DIV2 site as homodimer, while FOXC2 can only bind this site as a monomer. Our structural, biochemical and bioinformatics analyses not only provide a novel mechanism by which FOXO1 binds DNA as a homodimer, but also shed light on the target selection of forkhead transcription factors.

The Role of Forkhead Box Family in Bone Metabolism and Diseases

Forkhead box (Fox) family, an evolutionarily conserved family of transcription factors carrying the "Forkhead" motif, plays an indispensable role in human health and disease. Fox family genes are involved in cell differentiation, proliferation and apoptosis, embryonic development, aging, glucose and lipid metabolism, and immune regulation. The regulatory role of the Fox family in the context of bone metabolism and orthopedic diseases is an emerging research hotspot. In this review, we highlight the major molecular mechanisms underlying the regulatory role of Fox factors in bone metabolism, bone development, bone homeostasis, and bone diseases associated with inhibition or upregulation of Fox factors. In addition, we discuss the emerging evidence in the realm of Fox factor-based therapeutics.

Identification of Long Noncoding RNAs Involved in Differentiation and Survival of Vascular Smooth Muscle Cells

Long noncoding RNAs (lncRNAs) have recently been implicated in many pathophysiological cardiovascular processes, including vascular remodeling and atherosclerosis. However, the functional role of lncRNAs in the differentiation, proliferation, and apoptosis of vascular smooth muscle cells (VSMCs) is largely unknown. In this study, differentially expressed lncRNAs in synthetic and contractile human VSMCs were screened using RNA sequencing. Among the seven selected lncRNAs, the expression of MSC-AS1, MBNL1-AS1, and GAS6-AS2 was upregulated, whereas the expression of NR2F1-AS1, FUT8-AS1, FOXC2-AS1, and CTD-2207P18.2 was reduced upon VSMC differentiation. We focused on the NR2F1-AS1 and FOXC2-AS1 lncRNAs and showed that their knockdown significantly reduced the expression of smooth muscle contractile marker genes (ACTA2, CNN1, and TAGLN). Furthermore, FOXC2-AS1 was found to regulate cell proliferation and apoptosis through Akt/mTOR signaling, and affect Notch signaling, which is a key regulator of the contractile phenotype of VSMCs. Taken together, we identified novel lncRNAs involved in VSMC proliferation and differentiation and FOXC2-AS1 as a multifunctional regulator for vascular homeostasis and associated diseases.

Structural basis of tropifexor as a potent and selective agonist of farnesoid X receptor

Farnesoid X receptor (FXR) is considered as a potential target for the treatment of several liver disorders such as primary biliary cholangitis (PBC) and primary sclerosing cholangitis (PSC). Tropifexor is a highly potent and non-steroidal FXR agonist that has progressed into phase II clinical trials in patients with PBC. The clinical trials demonstrate that tropifexor improved serum markers of patients with liver diseases and lower side effect such as pruritus that might be implicated with TGR5 activation. However, the molecular mechanism of the potency and selectivity of tropifexor remains unclear. In this study, the binding affinity of FXR and tropifexor is measured by isothermal titration calorimetry (ITC) assays. The crystal structure of the FXR/tropifexor complex is determined at 2.7 Å resolution to explain the molecular mechanism of tropifexor bound to FXR-LBD. Structural comparison with other FXR/agonists structures reveals the conformational change in the FXR/tropifexor structure. Moreover the structural superposition of TGR5/tropifexor indicates that the steric hindrance between tropifexor and TGR5 might be a possible explanation for the impotency arises of tropifexor to TGR5. Overall, our analyses might provide an insight into the molecular mechanism of tropifexor binding to FXR-LBD and account for the high selectivity of tropifexor for FXR versus TGR5.

Structural Basis for DNA Recognition by FOXG1 and the Characterization of Disease-causing FOXG1 Mutations

Forkhead box G1 (FOXG1) is a transcription factor mainly expressed in the brain that plays a critical role in the development and regionalization of the forebrain. Aberrant expression of FOXG1 has implications in FOXG1 syndrome, a serious neurodevelopmental disorder. Here, we report the crystal structure of the FOXG1 DNA-binding domain (DBD) in complex with the forkhead consensus DNA site DBE2 at the resolution of 1.6 Å. FOXG1-DBD adopts a typical winged helix fold. Compared to those of other FOX-DBD/DBE2 structures, the N terminus, H3 helix and wing2 region of FOXG1-DBD exhibit differences in DNA recognition. The FOXG1-DBD wing2 region adopts a unique architecture composed of two β-strands that differs from all other known FOX-DBD wing2 folds. Mutation assays revealed that the disease-causing mutations within the FOXG1-DBD affect DNA binding, protein thermal stability, or both. Our report provides initial insight into how FOXG1 binds DNA and sheds light on how disease-causing mutations in FOXG1-DBD affect its DNA-binding ability.

Axenfeld-Rieger syndrome-associated mutants of the transcription factor FOXC1 abnormally regulate NKX2-5 in model zebrafish embryos

FOXC1 is a member of the forkhead family of transcription factors, and whose function is poorly understood. A variety of FOXC1 mutants have been identified in patients diagnosed with the autosomal dominant disease Axenfeld-Rieger syndrome, which is mainly characterized by abnormal development of the eyes, particularly those who also have accompanying congenital heart defects (CHD). However, the role of FOXC1 in CHD, and how these mutations might impact FOXC1 function, remains elusive. Our previous work provided one clue to possible function, demonstrating that zebrafish foxc1a, an orthologue of human FOXC1 essential for heart development, directly regulates the expression of nkx2.5, encoding a transcriptional regulator of cardiac progenitor cells. Abnormal expression of Nkx2-5 leads to CHD in mice and is also associated with CHD patients. Whether this link extends to the human system, however, requires investigation. In this study, we demonstrate that FOXC1 does regulate human NKX2-5 expression in a dose-dependent manner via direct binding to its proximal promoter. A comparison of FOXC1 mutant function in the rat cardiac cell line H9c2 and zebrafish embryos suggested that the zebrafish embryos might serve as a more representative model system than the H9c2 cells. Finally, we noted that three of the Axenfeld-Rieger syndrome FOXC1 mutations tested increased, whereas a fourth repressed the expression of NKX2-5 These results imply that mutant FOXC1s might play etiological roles in CHD by abnormally regulating NKX2-5 in the patients. And zebrafish embryos can serve as a useful in vivo platform for rapidly evaluating disease-causing roles of mutated genes.

Structural basis of binding of homodimers of the nuclear receptor NR4A2 to selective Nur-responsive DNA elements

Proteins of nuclear receptor subfamily 4 group A (NR4A), including NR4A1/NGFI-B, NR4A2/Nurr1, and NR4A3/NOR-1, are nuclear transcription factors that play important roles in metabolism, apoptosis, and proliferation. NR4A proteins recognize DNA response elements as monomers or dimers to regulate the transcription of a variety of genes involved in multiple biological processes. In this study, we determined two crystal structures of the NR4A2 DNA-binding domain (NR4A2-DBD) bound to two Nur-responsive elements: an inverted repeat and an everted repeat at 2.6-2.8 Å resolution. The structures revealed that two NR4A2-DBD molecules bind independently to the everted repeat, whereas two other NR4A2-DBD molecules form a novel dimer interface on the inverted repeat. Moreover, substitution of the interfacial residue valine 298 to lysine as well as mutation of DNA bases involved in the interactions abolished the dimerization. Overall, our structural, biochemical, and bioinformatics analyses provide a molecular basis for the binding of the NR4A2 protein dimers to NurREs and advance our understanding of the dimerization specificity of nuclear receptors.

Crystal Structure of FOXC2 in Complex with DNA Target

Forkhead transcription factor C2 (FOXC2) is a transcription factor regulating vascular and lymphatic development, and its mutations are linked to lymphedema-distichiasis syndrome. FOXC2 is also a crucial regulator of the epithelial-mesenchymal transition processes essential for tumor metastasis. Here, we report the crystal structure of the FOXC2-DNA-binding domain in complex with its cognate DNA. The crystal structure provides the basis of DNA sequence recognition by FOXC2 for the T/CAAAC motif. Helix 3 makes the majority of the DNA-protein interactions and confers the DNA sequence specificity. The computational energy calculation results also validate the structural observations. The FOXC2 and DNA complex structure provides a detailed picture of protein and DNA interactions, which allows us to predict its DNA recognition specificity and impaired functions in mutants identified in human patients.