Crystal structure of the human NKX2.5 homeodomain in complex with DNA target

DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

Prioritizing cardiovascular disease-associated variants altering NKX2-5 and TBX5 binding through an integrative computational approach

Cardiovascular diseases (CVDs) are the leading cause of death worldwide and are heavily influenced by genetic factors. Genome-wide association studies have mapped >90% of CVD-associated variants within the noncoding genome, which can alter the function of regulatory proteins, such as transcription factors (TFs). However, due to the overwhelming number of single-nucleotide polymorphisms (SNPs) (>500,000) in genome-wide association studies, prioritizing variants for in vitro analysis remains challenging. In this work, we implemented a computational approach that considers support vector machine (SVM)-based TF binding site classification and cardiac expression quantitative trait loci (eQTL) analysis to identify and prioritize potential CVD-causing SNPs. We identified 1535 CVD-associated SNPs within TF footprints and putative cardiac enhancers plus 14,218 variants in linkage disequilibrium with genotype-dependent gene expression in cardiac tissues. Using ChIP-seq data from two cardiac TFs (NKX2-5 and TBX5) in human-induced pluripotent stem cell-derived cardiomyocytes, we trained a large-scale gapped k-mer SVM model to identify CVD-associated SNPs that altered NKX2-5 and TBX5 binding. The model was tested by scoring human heart TF genomic footprints within putative enhancers and measuring in vitro binding through electrophoretic mobility shift assay. Five variants predicted to alter NKX2-5 (rs59310144, rs6715570, and rs61872084) and TBX5 (rs7612445 and rs7790964) binding were prioritized for in vitro validation based on the magnitude of the predicted change in binding and are in cardiac tissue eQTLs. All five variants altered NKX2-5 and TBX5 DNA binding. We present a bioinformatic approach that considers tissue-specific eQTL analysis and SVM-based TF binding site classification to prioritize CVD-associated variants for in vitro analysis.

Prioritizing Cardiovascular Disease-Associated Variants Altering NKX2-5 Binding through an Integrative Computational Approach

Cardiovascular diseases (CVDs) are the leading cause of death worldwide and are heavily influenced by genetic factors. Genome-wide association studies (GWAS) have mapped > 90% of CVD-associated variants within the non-coding genome, which can alter the function of regulatory proteins, like transcription factors (TFs). However, due to the overwhelming number of GWAS single nucleotide polymorphisms (SNPs) (>500,000), prioritizing variants for in vitro analysis remains challenging. In this work, we implemented a computational approach that considers support vector machine (SVM)-based TF binding site classification and cardiac expression quantitative trait loci (eQTL) analysis to identify and prioritize potential CVD-causing SNPs. We identified 1,535 CVD-associated SNPs that occur within human heart footprints/enhancers and 9,309 variants in linkage disequilibrium (LD) with differential gene expression profiles in cardiac tissue. Using hiPSC-CM ChIP-seq data from NKX2-5 and TBX5, two cardiac TFs essential for proper heart development, we trained a large-scale gapped k-mer SVM (LS-GKM-SVM) predictive model that can identify binding sites altered by CVD-associated SNPs. The computational predictive model was tested by scoring human heart footprints and enhancers in vitro through electrophoretic mobility shift assay (EMSA). Three variants (rs59310144, rs6715570, and rs61872084) were prioritized for in vitro validation based on their eQTL in cardiac tissue and LS-GKM-SVM prediction to alter NKX2-5 DNA binding. All three variants altered NKX2-5 DNA binding. In summary, we present a bioinformatic approach that considers tissue-specific eQTL analysis and SVM-based TF binding site classification to prioritize CVD-associated variants for in vitro experimental analysis.

Modeling a variant of unknown significance in the Drosophila ortholog of the human cardiogenic gene NKX2.5

Sequencing of human genome samples has unearthed genetic variants for which functional testing is necessary to validate their clinical significance. We used the Drosophila system to analyze a variant of unknown significance in the human congenital heart disease gene NKX2.5 (also known as NKX2-5). We generated an R321N allele of the NKX2.5 ortholog tinman (tin) to model a human K158N variant and tested its function in vitro and in vivo. The R321N Tin isoform bound poorly to DNA in vitro and was deficient in activating a Tin-dependent enhancer in tissue culture. Mutant Tin also showed a significantly reduced interaction with a Drosophila T-box cardiac factor named Dorsocross1. We generated a tinR321N allele using CRISPR/Cas9, for which homozygotes were viable and had normal heart specification, but showed defects in the differentiation of the adult heart that were exacerbated by further loss of tin function. We propose that the human K158N variant is pathogenic through causing a deficiency in DNA binding and a reduced ability to interact with a cardiac co-factor, and that cardiac defects might arise later in development or adult life.

Using Drosophila to model a variant of unknown significance in the human cardiogenic gene Nkx2.5

Sequencing of human genome samples has unearthed genetic variants for which functional testing is necessary to validate their clinical significance. We used the Drosophila system to analyze a variant of unknown significance in the human congenital heart disease gene, Nkx2 . 5 . We generated an R321N allele of the Nkx2 . 5 ortholog tinman ( tin ) to model a human K158N variant and tested its function in vitro and in vivo. The R321N Tin isoform bound poorly to DNA in vitro and was deficient in activating a Tin-dependent enhancer in tissue culture. Mutant Tin also showed a significantly reduced interaction with a Drosophila Tbox cardiac factor named Dorsocross1. We generated a tin R321N allele using CRISPR/Cas9, for which homozygotes were viable and had normal heart specification, but showed defects in the differentiation of the adult heart that were exacerbated by further loss of tin function. We conclude that the human K158N mutation is likely pathogenic through causing both a deficiency in DNA binding and a reduced ability to interact with a cardiac cofactor, and that cardiac defects might arise later in development or adult life.

Disease-associated non-coding variants alter NKX2-5 DNA-binding affinity

Genome-wide association studies (GWAS) have mapped over 90 % of disease- or trait-associated variants within the non-coding genome, like cis-regulatory elements (CREs). Non-coding single nucleotide polymorphisms (SNPs) are genomic variants that can change how DNA-binding regulatory proteins, like transcription factors (TFs), interact with the genome and regulate gene expression. NKX2-5 is a TF essential for proper heart development, and mutations affecting its function have been associated with congenital heart diseases (CHDs). However, establishing a causal mechanism between non-coding genomic variants and human disease remains challenging. To address this challenge, we identified 8475 SNPs predicted to alter NKX2-5 DNA-binding using a position weight matrix (PWM)-based predictive model. Five variants were prioritized for in vitro validation; four of them are associated with traits and diseases that impact cardiovascular health. The impact of these variants on NKX2-5 binding was evaluated with electrophoretic mobility shift assay (EMSA) using purified recombinant NKX2-5 homeodomain. Binding curves were constructed to determine changes in binding between variant and reference alleles. Variants rs7350789, rs7719885, rs747334, and rs3892630 increased binding affinity, whereas rs61216514 decreased binding by NKX2-5 when compared to the reference genome. Our findings suggest that differential TF-DNA binding affinity can be key in establishing a causal mechanism of pathogenic variants.

A de novo heterozygous HOXA11 variant in a patient with mesomelic dysplasia with urogenital abnormalities

Mesomelic dysplasias are a genetically and clinically heterogeneous group of diseases with more than 10 types defined. This article presents an 18-year-old female patient with normal intelligence and a multisystem phenotype including disproportionate short stature, scoliosis, mesomelic limb shortening, radial bowing, short fourth to fifth metacarpals and metatarsals, fusions in the carpal/tarsal bones, operated pes equinovarus, primary amenorrhea, uterine hypoplasia, vesicoureteral reflux, and chronic kidney disease. Whole-exome sequencing revealed a de novo heterozygous c.881T>G (p.Met294Arg) variant in HOXA11 (NM_005523.6) gene. The variant was located in the homeodomain of HOXA11 and predicted to alter DNA-binding ability of the protein. In silico analyses indicated that the variant could promote the alterations in the protein-protein interaction. The possible functional effect of the variant was supposed as dominant-negative. Hoxa11-mutant mice have been reported to exhibit homeotic transformations in the thoracic and sacral vertebrae, zeugopodal phenotype in forelimb and hindlimb, and urogenital abnormalities. Although mice models were reported as mesomelic dysplasia and urogenital abnormalities (MDUGA), this phenotype has not yet been reported in humans. This was the first case with MDUGA putatively related to a de novo variant in HOXA11.

The transcriptional regulator HDP1 controls expansion of the inner membrane complex during early sexual differentiation of malaria parasites

Transmission of Plasmodium falciparum and other malaria parasites requires their differentiation from asexual blood stages into gametocytes, the non-replicative sexual stage necessary to infect the mosquito vector. This transition involves changes in gene expression and chromatin reorganization that result in the activation and silencing of stage-specific genes. However, the genomes of malaria parasites have been noted for their limited number of transcriptional and chromatin regulators, and the molecular mediators of these changes remain largely unknown. We recently identified homeodomain protein 1 (HDP1) as a DNA-binding protein, first expressed in gametocytes, that enhances the expression of key genes critical for early sexual differentiation. The discovery of HDP1 marks a new class of transcriptional regulator in malaria parasites outside of the better-characterized ApiAP2 family. Here, using molecular biology, biochemistry and microscopy techniques, we show that HDP1 is essential for gametocyte maturation, facilitating the necessary upregulation of inner membrane complex components during early gametocytogenesis that gives P. falciparum gametocytes their characteristic shape.

Screening of NKX2.5 gene in Moroccan Tetralogy of Fallot (TOF) patients: worldwide mutation rate comparisons show a significant association between R25C variant and TOF phenotype

Background

Tetralogy of Fallot is the most prevalent cyanotic congenital heart disease, occurring in 1/3 600 live births. This disorder comprises ventricular septal defect, right ventricular outflow obstruction, over-riding aorta, and right ventricular hypertrophy. The present study aims to reveal the spectrum of Nk2 homeobox 5 (NKX2-5) variants identified in a Moroccan non-syndromic tetralogy of Fallot cohort and to compare mutation rate with different studies from all over the world. Thirty-one patients with non-syndromic tetralogy of Fallot were recruited in this cross-sectional study. DNAs were extracted, and coding regions of NKX2.5 were PCR-amplified and sequenced. The obtained sequences were analyzed using different bioinformatics tools. Statistical comparisons were carried out using the R software.

Results

R25C mutation was found in two patients, in association with the E21E variant. The latter variant was frequently observed in the population and seems to have a potential altering effect on the splicing process. The NKX2.5 mutation rate in our tetralogy of Fallot population is around 6.4%, and no significant difference was noticed in comparison with previous studies. At the same time, a comparison of R25C mutation rate between atrial septal defect and tetralogy of Fallot worldwide populations shows a particular association of R25C mutation with tetralogy of Fallot phenotype.

Conclusions

This study reveals a consistency between our NKX2.5 mutation rate and those of different tetralogy of Fallot populations around the world. Our findings suggest a possible combined effect of R25C mutation and E21E variant on the carriers and emphasize particularly the significant association of R25C mutation with tetralogy of Fallot, which highlights the importance of an anticipative screening for TOF phenotype among the carriers’ offspring at the perinatal period.

Case Report: A Novel NKX2-5 Mutation in a Family With Congenital Heart Defects, Left Ventricular Non-compaction, Conduction Disease, and Sudden Cardiac Death

The NKX2-5 gene encodes for a transcription factor crucial for cardiac cell differentiation and proliferation. It was the first gene associated with congenital heart disease (CHD) in humans and has been linked to conduction disorders or cardiomyopathies. However, an overlapping phenotype is not frequent in the literature. We describe a family with a novel missense mutation in the NKX2-5 gene (p.Gln181Pro) with numerous antecedents with atrial septal defect (ASD), left ventricular non-compaction (LVNC), conduction disease, and sudden cardiac death (SCD).

Crystal Structures of Ternary Complexes of MEF2 and NKX2-5 Bound to DNA Reveal a Disease Related Protein-Protein Interaction Interface

MEF2 and NKX2-5 transcription factors interact with each other in cardiogenesis and are necessary for normal heart formation. Despite evidence suggesting that these two transcription factors function synergistically and possibly through direct physical interactions, molecular mechanisms by which they interact are not clear. Here we determined the crystal structures of ternary complexes of MEF2 and NKX2-5 bound to myocardin enhancer DNA in two crystal forms. These crystal structures are the first example of human MADS-box/homeobox ternary complex structures involved in cardiogenesis. Our structures reveal two possible modes of interactions between MEF2 and NKX2-5: MEF2 and NKX bind to adjacent DNA sites to recognize DNA in cis; and MEF2 and NKX bind to different DNA strands to interact with each other in trans via a conserved protein-protein interface observed in both crystal forms. Disease-related mutations are mapped to the observed protein-protein interface. Our structural studies provide a starting point to understand and further study the molecular mechanisms of the interactions between MEF2 and NKX2.5 and their roles in cardiogenesis.

An update on genetic variants of the NKX2-5

NKX2-5 is a homeodomain transcription factor that plays a crucial role in heart development. It is the first gene where a single genetic variant (GV) was found to be associated with congenital heart diseases in humans. In this study, we carried out a comprehensive survey of NKX2-5 GVs to build a unified, curated, and updated compilation of all available GVs. We retrieved a total of 1,380 unique GVs. From these, 970 had information on their frequency in the general population and 143 have been linked to pathogenic phenotypes in humans. In vitro effect was ascertained for 38 GVs. The homeodomain had the biggest cluster of pathogenic variants in the protein: 49 GVs in 60 residues, 23 in its third α-helix, where 11 missense variants may affect protein-DNA interaction or the hydrophobic core. We also pinpointed the likely location of pathogenic GVs in four linear motifs. These analyses allowed us to assign a putative explanation for the effect of 90 GVs. This study pointed to reliable pathogenicity for GVs in helix 3 of the homeodomain and may broaden the scope of functional and structural studies that can be done to better understand the effect of GVs in NKX2-5 function.

Variants in NKX2-5 and FLNC Cause Dilated Cardiomyopathy and Sudden Cardiac Death

BACKGROUND

Dilated cardiomyopathy (DCM) is an important cause of heart failure. Variants in >50 genes have been reported to cause DCM, but causative variants have been found in less than half of familial cases. Variants causing DCM in Iceland have not been reported before.

METHODS

We performed a genome-wide association study on DCM based on whole genome sequencing. We tested the association of 32.5 million sequence variants in 424 cases and 337 689 population controls in Iceland.

RESULTS

We identified 2 DCM variants in established cardiomyopathy genes, a missense variant p.Phe145Leu in NKX2-5 carried by 1 in 7100 Icelanders ( P=7.0×10^-12) and a frameshift variant p.Phe1626Serfs*40 in FLNC carried by 1 in 3600 Icelanders ( P=2.1×10^-10). Both variants associate with heart failure and sudden cardiac death. Additionally, p.Phe145Leu in NKX2-5 associates with high degree atrioventricular block and atrial septal defect ( P<1.4×10^-4). The penetrance of serious heart disease among carriers of the NKX2-5 variant is high and higher than that of the FLNC variant.

CONCLUSIONS

Two rare variants in NKX2-5 and FLNC, carried by 1 in 2400 Icelanders, cause familial DCM in Iceland. These genes have recently been associated with DCM. Given the serious consequences of these variants, we suggest screening for them in individuals with DCM and their family members, with subsequent monitoring of carriers, offering early intervention.

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.

Can All-Atom Molecular Dynamics Simulations Quantitatively Describe Homeodomain-DNA Binding Equilibria?

We systematically investigate the applicability of a molecular dynamics-based setup for the calculations of standard binding free energies of biologically relevant protein-DNA complexes. The free energies are extracted from a potential of mean force calculated using umbrella sampling simulations. Two protein-DNA systems derived from a homeodomain transcription factor complex are studied in order to investigate the binding of both disordered and globular proteins. Free energies and trajectories obtained using two modern molecular mechanical force fields are compared to each other and to experimental data. The temperature dependence of the calculated standard binding free energies is investigated by performing all simulations over a range of temperatures. We show that the values of standard binding free energies obtained from these simulations are overestimated compared to experimental results. Significant differences are observed between the two protein-DNA systems and between the two force fields, which are explained by different propensities to form inter- and intramolecular contacts. The number of protein-DNA contacts increases with increasing temperature, in agreement with the experimentally known temperature dependence of enthalpies of binding. However, conclusions about the temperature dependence of the standard binding free energies cannot be made with confidence, as the differences among the values are on the order of statistical uncertainty.

Computational modeling suggests impaired interactions between NKX2.5 and GATA4 in individuals carrying a novel pathogenic D16N NKX2.5 mutation

NKX2.5, a homeobox containing gene, plays an important role in embryonic heart development and associated mutations are linked with various cardiac abnormalities. We sequenced the NKX2.5 gene in 100 congenital heart disease (CHD) patients and 200 controls. Our analysis revealed a total of 7 mutations, 3 in intronic region, 3 in coding region and 1 in 3’ UTR. Of the above mutations, one mutation was found to be associated with tetralogy of fallot (TOF) and two (rs2277923 and a novel mutation, D16N) were strongly associated with VSD. A novel missense mutation, D16N (p-value =0.009744), located in the tinman (TN) region and associated with ventricular septal defect (VSD), is the most significant findings of this study. Computational analysis revealed that D16N mutation is pathogenic in nature. Through the molecular modeling, docking and molecular dynamics simulation studies, we have identified the location of mutant D16N in NKX2.5 and its interaction map with other partners at the atomic level. We found NKX2.5-GATA4 complex is stable, however, in case of mutant we observed significant conformational changes and loss of key polar interactions, which might be a cause of the pathogenic behavior. This study underscores the structural basis of D16N pathogenic mutation in the regulation of NKX2.5 and how this mutation renders the structural-functional divergence that possibly leading towards the diseased state.

Biallelic mutations in the homeodomain of NKX6-2 underlie a severe hypomyelinating leukodystrophy

Hypomyelinating leukodystrophies are genetically heterogeneous disorders with overlapping clinical and neuroimaging features reflecting variable abnormalities in myelin formation. We report on the identification of biallelic inactivating mutations in NKX6-2, a gene encoding a transcription factor regulating multiple developmental processes with a main role in oligodendrocyte differentiation and regulation of myelin-specific gene expression, as the cause underlying a previously unrecognized severe variant of hypomyelinating leukodystrophy. Five affected subjects (three unrelated families) were documented to share biallelic inactivating mutations affecting the NKX6-2 homeobox domain. A trio-based whole exome sequencing analysis in the first family detected a homozygous frameshift change [c.606delinsTA; p.(Lys202Asnfs*?)]. In the second family, homozygosity mapping coupled to whole exome sequencing identified a homozygous nucleotide substitution (c.565G>T) introducing a premature stop codon (p.Glu189*). In the third family, whole exome sequencing established compound heterozygosity for a non-conservative missense change affecting a key residue participating in DNA binding (c.599G>A; p.Arg200Gln) and a nonsense substitution (c.589C>T; p.Gln197*), in both affected siblings. The clinical presentation was homogeneous, with four subjects having severe motor delays, nystagmus and absent head control, and one individual showing gross motor delay at the age of 6 months. All exhibited neuroimaging that was consistent with hypomyelination. These findings define a novel, severe form of leukodystrophy caused by impaired NKX6-2 function.

New Insights into Cooperative Binding of Homeodomain Transcription Factors PREP1 and PBX1 to DNA

PREP1 and PBX1 are homeodomain (HD) transcription factors that play crucial roles in embryonic development. Here, we present the first biophysical characterization of a PREP1 HD, and the NMR spectroscopic study of its DNA binding pocket. The data show that residues flanking the HD participate in DNA binding. The kinetic parameters for DNA binding of individual PREP1 and PBX1 HDs, and of their combination, show that isolated PREP1 and PBX1 HDs bind to DNA in a cooperative manner. A novel PREP1 motif, flanking the HD at the C-terminus, is required for cooperativity.

Complex Interdependence Regulates Heterotypic Transcription Factor Distribution and Coordinates Cardiogenesis

Transcription factors (TFs) are thought to function with partners to achieve specificity and precise quantitative outputs. In the developing heart, heterotypic TF interactions, such as between the T-box TF TBX5 and the homeodomain TF NKX2-5, have been proposed as a mechanism for human congenital heart defects. We report extensive and complex interdependent genomic occupancy of TBX5, NKX2-5, and the zinc finger TF GATA4 coordinately controlling cardiac gene expression, differentiation, and morphogenesis. Interdependent binding serves not only to co-regulate gene expression but also to prevent TFs from distributing to ectopic loci and activate lineage-inappropriate genes. We define preferential motif arrangements for TBX5 and NKX2-5 cooperative binding sites, supported at the atomic level by their co-crystal structure bound to DNA, revealing a direct interaction between the two factors and induced DNA bending. Complex interdependent binding mechanisms reveal tightly regulated TF genomic distribution and define a combinatorial logic for heterotypic TF regulation of differentiation.

A direct fate exclusion mechanism by Sonic hedgehog-regulated transcriptional repressors

Sonic hedgehog (Shh) signaling patterns the vertebrate spinal cord by activating a group of transcriptional repressors in distinct neural progenitors of somatic motor neuron and interneuron subtypes. To identify the action of this network, we performed a genome-wide analysis of the regulatory actions of three key ventral determinants in mammalian neural tube patterning: Nkx2.2, Nkx6.1 and Olig2. Previous studies have demonstrated that each factor acts predominantly as a transcriptional repressor, at least in part, to inhibit alternative progenitor fate choices. Here, we reveal broad and direct repression of multiple alternative fates as a general mechanism of repressor action. Additionally, the repressor network targets multiple Shh signaling components providing negative feedback to ongoing Shh signaling. Analysis of chromatin organization around Nkx2.2-, Nkx6.1- and Olig2-bound regions, together with co-analysis of engagement of the transcriptional activator Sox2, indicate that repressors bind to, and probably modulate the action of, neural enhancers. Together, the data suggest a model for neural progenitor specification downstream of Shh signaling, in which Nkx2.2 and Olig2 direct repression of alternative neural progenitor fate determinants, an action augmented by the overlapping activity of Nkx6.1 in each cell type. Integration of repressor and activator inputs, notably activator inputs mediated by Sox2, is probably a key mechanism in achieving cell type-specific transcriptional outcomes in mammalian neural progenitor fate specification.

A novel NKX2.6 mutation associated with congenital ventricular septal defect

Congenital heart disease (CHD) is the most common birth defect and is the most prevalent non-infectious cause of infant death. Aggregating evidence demonstrates that genetic defects are involved in the pathogenesis of CHD. However, CHD is genetically heterogeneous and the genetic determinants for CHD in an overwhelming majority of patients remain unknown. In this study, the coding regions and splice junctions of the NKX2.6 gene, which encodes a homeodomain transcription factor crucial for cardiovascular development, were sequenced in 210 unrelated CHD patients. As a result, a novel heterozygous NKX2.6 mutation, p.K152Q, was identified in an index patient with ventricular septal defect (VSD). Genetic analysis of the proband's available family members showed that the mutation cosegregated with VSD transmitted as an autosomal dominant trait with complete penetrance. The missense mutation was absent in 400 control chromosomes and the altered amino acid was completely conserved evolutionarily across species. Due to unknown transcriptional targets of NKX2.6, the functional characteristics of the identified mutation at transcriptional activity were analyzed by using NKX2.5 as a surrogate. Alignment between human NKX2.6 and NKX2.5 proteins displayed that K152Q-mutant NKX2.6 was equivalent to K158Q-mutant NKX2.5, and introduction of K158Q into NKX2.5 significantly reduced its transcriptional activating function when compared with its wild-type counterpart. This study firstly links NKX2.6 loss-of-function mutation with increased susceptibility to isolated VSD, providing novel insight into the molecular mechanism underpinning VSD and contributing to the development of new preventive and therapeutic strategies for this common form of CHD.

NuProPlot: nucleic acid and protein interaction analysis and plotting program

Growing numbers of protein and nucleic acid complex structures are being determined and deposited in the Protein Data Bank and the Nucleic Acid Database. With the increasing complexity of these structures, it is challenging to analyse and visualize the three-dimensional interactions. The currently available programs for such analysis and visualization are limited in their applications. They can only analyse a subset of protein–nucleic acid complexes and require multiple iterations before obtaining plots that are suitable for presentation. An interactive web-based program,NuProPlot(http://www.nuproplot.com), has been developed which can automatically identify hydrogen, electrostatic and van der Waals interactions between proteins and nucleic acids and generate a plot showing all of the interactions. Protein–DNA and protein–RNA interactions can be visualized in simple two-dimensional schematics. Interactive schematic drawing options allow selection of the plotted area and repositioning of the individual interactions for better legibility.NuProPlotis a fully automated and user-friendly program providing various custom options.NuProPlotrepresents a greatly improved option for analysis and presentation of protein–nucleic acid interactions.

Role of DNA conformation & energetic insights in Msx-1-DNA recognition as revealed by molecular dynamics studies on specific and nonspecific complexes

In most of homeodomain-DNA complexes, glutamine or lysine is present at 50th position and interacts with 5th and 6th nucleotide of core recognition region. Molecular dynamics simulations of Msx-1-DNA complex (Q50-TG) and its variant complexes, that is specific (Q50K-CC), nonspecific (Q50-CC) having mutation in DNA and (Q50K-TG) in protein, have been carried out. Analysis of protein-DNA interactions and structure of DNA in specific and nonspecific complexes show that amino acid residues use sequence-dependent shape of DNA to interact. The binding free energies of all four complexes were analysed to define role of amino acid residue at 50th position in terms of binding strength considering the variation in DNA on stability of protein-DNA complexes. The order of stability of protein-DNA complexes shows that specific complexes are more stable than nonspecific ones. Decomposition analysis shows that N-terminal amino acid residues have been found to contribute maximally in binding free energy of protein-DNA complexes. Among specific protein-DNA complexes, K50 contributes more as compared to Q50 towards binding free energy in respective complexes. The sequence dependence of local conformation of DNA enables Q50/Q50K to make hydrogen bond with nucleotide(s) of DNA. The changes in amino acid sequence of protein are accommodated and stabilized around TAAT core region of DNA having variation in nucleotides.

A novel NKX2.5 loss-of-function mutation associated with congenital bicuspid aortic valve

Bicuspid aortic valve (BAV) is the most common form of congenital cardiovascular defect in humans and is associated with substantial morbidity and mortality. Emerging evidence demonstrates that genetic risk factors play an important role in the pathogenesis of BAV. However, BAV is a genetically heterogenous disorder, and the genetic defects underpinning BAV in most patients remain to be identified. In the present study, the coding exons and flanking introns of the NKX2.5 gene, which encodes a homeodomain-containing transcription factor essential for the normal development of the aortic valve, were sequenced in 142 unrelated patients with BAV. The available relatives of the mutation carrier and 200 unrelated healthy subjects used as controls were also genotyped for NKX2.5. The functional characteristics of the mutation were delineated by using a dual-luciferase reporter assay system. As a result, a novel heterozygous NKX2.5 mutation, p.K192X, was identified in a family with BAV transmitted in an autosomal dominant pattern. The nonsense mutation was absent in 400 control chromosomes. Functional analyses revealed that the mutant NKX2.5 had no transcriptional activity compared with its wild-type counterpart. Furthermore, the mutation abolished the synergistic transcriptional activation between NKX2.5 and GATA5, another transcription factor crucial for the aortic valvular morphogenesis. In conclusion, this study is the first to link an NKX2.5 loss-of-function mutation to enhanced susceptibility to human BAV, providing novel insight into the molecular mechanism of BAV and suggesting potential implications for genetic counseling and clinical care of families presenting with BAV.

A novel LMX1B mutation in a family with end-stage renal disease of 'unknown cause'

End-stage renal disease (ESRD) presenting in a familial autosomal dominant pattern points to an underlying monogenic cause. Nail-patella syndrome (NPS) is an autosomal dominant disorder that may lead to ESRD caused by mutations in the transcription factor LMX1B. Renal-limited forms of this disease, termed nail-patella-like renal disease (NPLRD), and LMX1B nephropathy have recently been described. We report a large family, from the North East of England, with seven affected members with varying phenotypes of renal disease, ranging from ESRD at 28 years of age to microscopic haematuria and proteinuria and relatively preserved renal function. In this family, there were no extra-renal manifestations to suggest NPS. Genome-wide linkage studies and inheritance by descent (IBD) suggested disease loci on Chromosome 1 and 9. Whole exome sequencing (WES) analysis identified a novel sequence variant (p.R249Q) in the LMX1B gene in each of the three samples submitted, which was confirmed using Sanger sequencing. The variant segregated with the disease in all affected individuals. In silico modelling revealed that R249 is putatively located in close proximity to the DNA phosphoskeleton, supporting a role for this residue in the interaction between the LMX1B homeodomain and its target DNA. WES and analysis of potential target genes, including CD2AP, NPHS2, COL4A3, COL4A4 and COL4A5, did not reveal any co-inherited pathogenic variants. In conclusion, we confirm a novel LMX1B mutation in a large family with an autosomal dominant pattern of nephropathy. This report confirms that LMX1B mutations may cause a glomerulopathy without extra-renal manifestations. A molecular genetic diagnosis of LMX1B nephropathy thus provides a definitive diagnosis, prevents the need for renal biopsies and allows at risk family members to be screened.

NKX2-6 mutation predisposes to familial atrial fibrillation

Atrial fibrillation (AF) is the most common form of sustained cardiac arrhythmia and is associated with substantially increased morbidity and mortality rates. Aggregating evidence demonstrates that genetic defects are involved in the pathogenesis of AF and a number of AF-associated genes have been identified. Nevertheless, AF is a genetically heterogeneous disorder and the genetic components underpinning AF in an overwhelming majority of patients remain unclear. In this study, the entire coding exons and splice junction sites of the NK2 homeobox 6 (NKX2-6) gene, which encodes a homeodomain transcription factor important for cardiovascular development, were sequenced in 150 unrelated patients with lone AF, and a novel heterozygous NKX2-6 mutation, p.Q175H, was identified in an index patient. Genetic analysis of the available family members of the mutation carrier revealed that the mutation co-segregated with AF transmitted in an autosomal dominant pattern. The missense mutation was absent in the 200 unrelated ethnically matched healthy individuals used as controls and the altered amino acid was completely conserved evolutionarily among species. Due to unknown transcriptional targets of NKX2-6, the functional characteristics of the mutation as regards transcriptional activity were analyzed using NKX2-5 as a surrogate. Alignment between human NKX2-6 and NKX2-5 proteins displayed that the Q175H-mutant NKX2-6 was equivalent to the Q181H-mutant NKX2-5, and the introduction of Q181H into NKX2-5 significantly decreased its transcriptional activity at the atrial natriuretic factor promoter. The present study firstly associates genetically defective NKX2-6 with enhanced susceptibility to AF, providing novel insight into the molecular mechanisms underlying AF and suggesting potential strategies for the antenatal prophylaxis and personalized treatment of AF.

A mouse model of human congenital heart disease: high incidence of diverse cardiac anomalies and ventricular noncompaction produced by heterozygous Nkx2-5 homeodomain missense mutation

BACKGROUND

Heterozygous human mutations of NKX2-5 are highly penetrant and associated with varied congenital heart defects. The heterozygous knockout of murine Nkx2-5, in contrast, manifests less profound cardiac malformations, with low disease penetrance. We sought to study this apparent discrepancy between human and mouse genetics. Because missense mutations in the NKX2-5 homeodomain (DNA-binding domain) are the most frequently reported type of human mutation, we replicated this genetic defect in a murine knockin model.

METHODS AND RESULTS

We generated a murine model in a 129/Sv genetic background by knocking-in an Nkx2-5 homeodomain missense mutation previously identified in humans. The mutation was located at homeodomain position 52Arg→Gly (R52G). All the heterozygous neonatal Nkx2-5(+/R52G) mice demonstrated a prominent trabecular layer in the ventricular wall, so called noncompaction, along with diverse cardiac anomalies, including atrioventricular septal defects, Ebstein malformation of the tricuspid valve, and perimembranous and muscular ventricular septal defects. In addition, P10 Nkx2-5(+/R52G) mice demonstrated atrial sepal anomalies, with significant increase in the size of the interatrial communication and fossa ovalis, and decrease in the length of the flap valve compared with control Nkx2-5(+/+) or Nkx2-5(+/-) mice.

CONCLUSIONS

The results of our study demonstrate that heterozygous missense mutation in the murine Nkx2-5 homeodomain (R52G) is highly penetrant and result in pleiotropic cardiac effects. Thus, in contrast to heterozygous Nkx2-5 knockout mice, the effects of the heterozygous knockin mimic findings in humans with heterozygous missense mutation in NKX2-5 homeodomain.

CEH-28 activates dbl-1 expression and TGF-β signaling in the C. elegans M4 neuron

M4 is a multifunctional neuron in the Caenorhabditis elegans pharynx that can both stimulate peristaltic contractions of the muscles in the pharyngeal isthmus and function systemically to regulate an enhanced sensory response under hypoxic conditions. Here we identify a third function for M4 that depends on activation of the TGF-β family gene dbl-1 by the homeodomain transcription factor CEH-28. dbl-1 is expressed in M4 and a subset of other neurons, and we show CEH-28 specifically activates dbl-1 expression in M4. Characterization of the dbl-1 promoter indicates that CEH-28 targets an M4-specific enhancer within the dbl-1 promoter region, while expression in other neurons is mediated by separate regulatory sequences. Unlike ceh-28 mutants, dbl-1 mutants do not exhibit M4 synaptic and signaling defects. Instead, both ceh-28 and dbl-1 mutants exhibit morphological defects in the g1 gland cells located adjacent to M4 in the pharynx, and these defects can be partially rescued by M4-specific expression of dbl-1 in these mutants. Identical gland cell defects are observed in sma-6 and daf-4 mutants defective in the receptor for DBL-1, but they are not observed in sma-2 and sma-3 mutants lacking the R-Smads functioning downstream of this receptor. Together these results identify a novel neuroendocrine function for M4 and provide evidence for an R-Smad-independent mechanism for DBL-1 signaling in C. elegans.

Structure solution of DNA-binding proteins and complexes with ARCIMBOLDO libraries

Protein-DNA interactions play a major role in all aspects of genetic activity within an organism, such as transcription, packaging, rearrangement, replication and repair. The molecular detail of protein-DNA interactions can be best visualized through crystallography, and structures emphasizing insight into the principles of binding and base-sequence recognition are essential to understanding the subtleties of the underlying mechanisms. An increasing number of high-quality DNA-binding protein structure determinations have been witnessed despite the fact that the crystallographic particularities of nucleic acids tend to pose specific challenges to methods primarily developed for proteins. Crystallographic structure solution of protein-DNA complexes therefore remains a challenging area that is in need of optimized experimental and computational methods. The potential of the structure-solution program ARCIMBOLDO for the solution of protein-DNA complexes has therefore been assessed. The method is based on the combination of locating small, very accurate fragments using the program Phaser and density modification with the program SHELXE. Whereas for typical proteins main-chain α-helices provide the ideal, almost ubiquitous, small fragments to start searches, in the case of DNA complexes the binding motifs and DNA double helix constitute suitable search fragments. The aim of this work is to provide an effective library of search fragments as well as to determine the optimal ARCIMBOLDO strategy for the solution of this class of structures.

Crystallization and preliminary X-ray analysis of the cardiac transcription factor complex of NKX2.5 and TBX5 with DNA

Heart development depends on timely expression of genes regulated by combinatorial interactions of master cardiac transcription factors. To elucidate the molecular basis of their interactions, a ternary complex of cardiac transcription factors, NKX2.5 and TBX5, and their target DNA was studied using X-ray crystallography. Here, the purification, crystallization and preliminary X-ray crystallographic analyses of the NKX2.5 homeodomain and TBX5 DNA-binding domain complex with a DNA element from the -252 promoter region of the atrial natriuretic factor are reported. The crystal diffracted to 2.88 Å resolution and belonged to space group P21, with unit-cell parameters a = 69.30, b = 77.78, c = 77.60 Å, β = 108.31°. Two sets of ternary complexes are present in an asymmetric unit with a solvent content of 54%.

Mutational spectrum of the NKX2-5 gene in patients with lone atrial fibrillation

Atrial fibrillation (AF) is the most common form of sustained cardiac arrhythmia in humans and is responsible for substantial morbidity and mortality worldwide. Emerging evidence indicates that abnormal cardiovascular development is involved in the pathogenesis of AF. In this study, the coding exons and splice sites of the NKX2-5 gene, which encodes a homeodomain-containing transcription factor essential for cardiovascular genesis, were sequenced in 146 unrelated patients with lone AF as well as the available relatives of the mutation carriers. A total of 700 unrelated ethnically matched healthy individuals used as controls were genotyped. The disease-causing potential of the identified NKX2-5 variations was predicted by MutationTaster and PolyPhen-2. The functional characteristics of the mutant NKX2-5 proteins were analyzed using a dual-luciferase reporter assay system. As a result, two heterozygous NKX2-5 mutations, including a previously reported p.E21Q and a novel p.T180A mutation, were identified in two families with AF transmitted in an autosomal dominant pattern. The mutations co-segregated with AF in the families with complete penetrance. The detected substitutions, which altered the amino acids highly conserved evolutionarily across species, were absent in 700 control individuals and were both predicted to be causative. Functional analyses demonstrated that the NKX2-5 mutants were associated with significantly decreased transcriptional activity compared with their wild-type counterpart. The findings expand the spectrum of NKX2-5 mutations linked to AF and provide additional evidence that dysfunctional NKX2-5 may confer vulnerability to AF, suggesting the potential benefit for the early prophylaxis and personalized treatment of AF.

Prevalence and spectrum of Nkx2.5 mutations associated with idiopathic atrial fibrillation

OBJECTIVE

The aim of this study was to evaluate the prevalence and spectrum of Nkx2.5 mutations associated with idiopathic atrial fibrillation (AF).

METHODS

A cohort of 136 unrelated patients with idiopathic atrial fibrillation and 200 unrelated, ethnically matched healthy controls were enrolled. The coding exons and splice junctions of the Nkx2.5 gene were sequenced in 136 atrial fibrillation patients, and the available relatives of mutation carriers and 200 controls were subsequently genotyped for the identified mutations. The functional characteristics of the mutated Nkx2.5 gene were analyzed using a dual-luciferase reporter assay system.

RESULTS

Two novel heterozygous Nkx2.5 mutations (p.N19D and p.F186S) were identified in 2 of the 136 unrelated atrial fibrillation cases, with a mutational prevalence of approximately 1.47%. These missense mutations co-segregated with atrial fibrillation in the families and were absent in the 400 control chromosomes. Notably, 2 mutation carriers also had congenital atrial septal defects and atrioventricular block. Multiple alignments of the Nkx2.5 protein sequences across various species revealed that the altered amino acids were completely conserved evolutionarily. Functional analysis demonstrated that the mutant Nkx2.5 proteins were associated with significantly reduced transcriptional activity compared to their wild-type counterpart.

CONCLUSION

These findings associate the Nkx2.5 loss-of-function mutation with atrial fibrillation and atrioventricular block and provide novel insights into the molecular mechanism involved in the pathogenesis of atrial fibrillation. These results also have potential implications for early prophylaxis and allele-specific therapy of this common arrhythmia.

Functional characterization of a novel mutation in NKX2-5 associated with congenital heart disease and adult-onset cardiomyopathy

BACKGROUND

The transcription factor NKX2-5 is crucial for heart development, and mutations in this gene have been implicated in diverse congenital heart diseases and conduction defects in mouse models and humans. Whether NKX2-5 mutations have a role in adult-onset heart disease is unknown.

METHODS AND RESULTS

Mutation screening was performed in 220 probands with adult-onset dilated cardiomyopathy. Six NKX2-5 coding sequence variants were identified, including 3 nonsynonymous variants. A novel heterozygous mutation, I184M, located within the NKX2-5 homeodomain, was identified in 1 family. A subset of family members had congenital heart disease, but there was an unexpectedly high prevalence of dilated cardiomyopathy. Functional analysis of I184M in vitro demonstrated a striking increase in protein expression when transfected into COS-7 cells or HL-1 cardiomyocytes because of reduced degradation by the Ubiquitin-proteasome system. In functional assays, DNA-binding activity of I184M was reduced, resulting in impaired activation of target genes despite increased expression levels of mutant protein.

CONCLUSIONS

Certain NKX2-5 homeodomain mutations show abnormal protein degradation via the Ubiquitin-proteasome system and partially impaired transcriptional activity. We propose that this class of mutation can impair heart development and mature heart function and contribute to NKX2-5-related cardiomyopathies with graded severity.

A novel NKX2.5 loss-of-function mutation responsible for familial atrial fibrillation

Atrial fibrillation (AF) represents the most common form of sustained cardiac arrhythmia and accounts for substantial morbidity and mortality. Increasing evidence demonstrates that abnormal cardiovascular development is involved in the pathogenesis of AF. In this study, the coding exons and splice sites of the NKX2.5 gene, which encodes a homeodomain-containing transcription factor pivotal for normal cardiovascular morphogenesis, were sequenced in 110 unrelated index patients with familial AF. The available relatives of the mutation carrier and 200 unrelated ethnically-matched healthy individuals serving as controls were subsequently genotyped. The disease-causing potential of the identified NKX2.5 variation was predicted by MutationTaster. The functional characteristics of the mutant NKX2.5 protein were analyzed using a dual-luciferase reporter assay system. As a result, a novel heterozygous NKX2.5 mutation, p.F145S, was identified in a family with AF transmitted as an autosomal dominant trait, which co-segregated with AF in the family with complete penetrance. The detected substitution, which altered the amino acid completely conserved evolutionarily across species, was absent in 400 control chromosomes and was automatically predicted to be causative. Functional analysis demonstrated that the NKX2.5 mutant was associated with significantly decreased transcriptional activity compared with its wild-type counterpart. To the best of our knowledge, this is the first report on the association of the NKX2.5 loss-of-function mutation with increased susceptibility to familial AF. The findings of the present study provide novel insights into the molecular mechanism underlying AF, suggesting the potential implications for the early prophylaxis and allele-specific therapy of AF.