|
|
||||||||
1From the Departments of Ophthalmology and 2Medical Genetics, University of Alberta, Edmonton, Alberta, Canada; 3The New York Eye and Ear Infirmary, New York, New York; 4New York Medical College, Valhalla, New York; and the 5Department of Medical Genetics, University of British Columbia, Vancouver, British Columbia, Canada.
| Abstract |
|---|
|
|
|---|
METHODS. Sequencing DNA from patients with Axenfeld-Rieger malformation resulted in the identification of two novel missense mutations (G165R and R169P) in wing 2 of FOXC1. Site-directed mutagenesis was used to introduce these mutations, as well as previously reported mutation (M161K), into the FOXC1 cDNA. These FOXC1 mutants were evaluated to determine their ability to localize to the nucleus, bind DNA and activate gene expression.
RESULTS. Two novel missense mutations were identified in unrelated patients, in wing 2 of the FOXC1 forkhead domain. Because there had been no previous biochemical analysis, the mutation M161K was also investigated. All three mutant proteins localized correctly to the nucleus. The G165R mutation maintained wild-type levels of DNA binding; however, both the M161K and R169P mutations displayed reduced DNA binding ability. Biochemical analysis showed that all three mutations disrupt FOXC1s transactivation ability.
CONCLUSIONS. Biochemical analysis of mutations G165R and R169P and of a previously reported mutation, M161K, demonstrate the functional significance of wing 2. M161K and R169P disrupt DNA binding of FOXC1, consistent with the hypothesis that wing 2 is necessary for DNA binding. The results also suggest that wing 2 plays a role in gene activation. These results provide the first insights into how mutations in wing 2 disrupt FOXC1 function.
FOX (forkhead box) proteins are a unique subset of the wHTH proteins characterized by the presence of a loop C-terminal to ß-strands called wing 2, that is essential for FOX/winged helix proteins to bind DNA as a monomer.9 ,10 The three-dimensional solution structure of FOXA3/DNA has been determined, and the structure depicts wing 2 crossing the three-helix bundle in the vicinity of the N terminus of
-helix 1 and the C terminus of
-helix 3.1 The region spanning wings 1 and 2 appears flexible and may provide a way of connecting subdomains that produce contacts necessary for DNA binding.3 11
Stevens et al.11 proposed an evolutionary pattern for the formation of wing 2 from an initial unstructured distal region in a basic wHTH protein. They further suggested that the acquisition of bulky aromatic and/or hydrophobic residues in this distal region allows for association with the hydrophobic and/or aromatic residues in the three-helix bundle. The ability to bind DNA as a monomer is a functional advantage that may have strengthened associations between residues in wing 2 and
-helix 1 that resulted in a functional dependence on wing 2 for stable DNA binding.11 However, the consequences of alteration of wing 2 on FOX protein function has not been studied.
Mutations in FOXC1 lead to Axenfeld-Rieger (AR) malformations, a clinically heterogeneous condition with variable penetrance.12 13 14 15 16 17 18 19 20 Patients manifest posterior embryotoxon, adhesions of the iris and the cornea, iris hypoplasia, corectopia, and/or an underdeveloped or aberrantly developed angle between the cornea and the iris. Nonocular features such as dental anomalies and redundant periumbilical skin may also be present. On rare occasions, patients have congenital cardiac anomalies.12 15 21 22 23 Glaucoma is a progressive blinding condition that occurs in approximately half of patients with AR malformations.
A molecular dissection of FOXC1 has shown that disease-causing missense mutations in the FHD disrupt FOXC1 function.5 7 Mutations in
-helix 1, between
-helix 2 and 4, and those that reside in the recognition helix
-helix 3, have been characterized.5 7 These mutations have been found to impede the ability of FOXC1 to bind DNA and activate transcription as well as to reduce levels of protein produced.5 7 Molecular analysis of disease-causing missense mutations has enabled possible functional domains to be assigned to different regions in the FOXC1 FHD.5 7
We report herein, the identification, patient phenotype, and molecular analysis of two novel and one previously reported3 4 disease-causing missense mutations in wing 2 of the FHD of FOXC1. Molecular analyses reveal that all three mutations reduce the ability of FOXC1 to transactivate gene expression and that the R169P and M161K mutations also disrupt the ability of FOXC1 to bind DNA. Analysis of these mutations demonstrates the importance of wing 2, a poorly characterized region, to the functional integrity of FOXC1 DNA binding and FOX/wHTH proteins in general.
| Methods |
|---|
|
|
|---|
|
Plasmid Construction
FOXC1 pcDNA4 His/Max B TM (Invitrogen, Carlsbad, CA) has been described previously.5 Mutagenesis was performed with a kit (Quickchange; Stratagene, La Jolla, CA) with the addition of 10% dimethylsulfoxide and the appropriate primers. The mutagenic primer sequences were as follows: M161K, forward 5'-ccg gac tcc tac aac aag ttc gag aac ggc agc-3' and reverse 5'gct gcc gtt ctc gaa ctt gtt gta gga gtc cgg-3'; G165R, forward 5'-aac atg ttc gag aac cgc agc ttc ctg cgg cgg-3' and reverse 5'-ccg ccg cag gaa gct gcg gtt ctc gaa cat gtt-3'; and R169P, forward 5'-gag aac ggc agc ttc ctg ccg cgg cgg cgg cgc ttc aag-3' and reverse 5'-ctt gaa gcg ccg ccg ccg cgg cag gaa gct gcc gtt ctc-3'. Potential mutant constructs were sequenced with an automated sequencer (LI-COR, Lincoln, NE) and confirmed mutants were subcloned into the FOXC1 pcDNA4 His/Max vector and resequenced.
Cell Culture
COS-7 and HeLa cells were maintained in Dulbeccos modified Eagles medium supplemented with 10% fetal bovine serum. Transfections were performed with a reagent (FuGene6; Roche Diagnostics, Indianapolis, IN), according to the manufacturers protocol, with 2 µg DNA per 100-mm plate and 0.5 µg DNA per 25-mm well.
Immunofluorescence
COS-7 cells were grown and transfected on coverslips with 1 x 105 cell/mL placed into six-well tissue culture plates. After 24 hours, immunofluorescence was performed as described previously.5
Electrophoretic Mobility Shift Assays
COS-7 extracts containing recombinant FOXC1 were equalized for amounts of recombinant FOXC1 protein by Western analysis using the anti-Xpress antibody (Stratagene) against the vector encoded N-terminal Xpress tag. Protein extracts were incubated with 1.25 mM dithiothreitol (DTT); 0.3 µg sheared salmon sperm DNA; 0.125 µg poly dI.dC (Sigma-Aldrich); and 80,000 cpm 32P-dCTPlabeled, double-stranded DNA containing the in vitro preferred FOXC1-binding site (shown in italic; forward 5'-gat cca aag taa ata aa caa cag a-3'; reverse 5'-gat ctc tgt tgt tta ttt act ttg-3'). Reactions were incubated at room temperature for 30 minutes and resolved on a 6% polyacrylamide Tris/glycine/EDTA gel.
Dual Luciferase Reporter Assays
Transactivation assays were performed as described previously (Dual Luciferase Reporter Assay System; Promega, Madison, WI).5
Modeling
The homology model of FOXC1 was created by aligning residues 76-177 of FOXC1 (GenBank Accession Q12948) with residues 117-218 of the crystal structure of FOXA31 in Swiss-PdbViewer.24 These sequences align with no gaps (Fig. 2B ). Although FOXC2, whose abbreviated FHD was solved by nuclear magnetic resonance (NMR) spectroscopy,25 shares very high sequence identity with FOXC1 (95.4%; Fig. 2B ), only FOXA3 (68.8% identity) was used as a template for homology modeling of FOXC1, because that purified peptide contained the most complete wing 2. The backbone atoms of FOXC1 were "fitted" against the template, and the modeling project was submitted to the SWISS-MODEL server (http://www.expasy.org/spdbv/; provided in the public domain by the Swiss Institute of Bioinformatics, Geneva, Switzerland). The resultant FOXC1 layer was subjected to in silico mutagenesis using the "mutate" tool of the Swiss-PdbViewer. For each amino acid, all pair-wise interatomic contacts of less than 0.7 nm (between residues that are greater than 11 positions apart), including a term for solvent accessibility, are converted to pseudoenergy terms (analogous to threading analysis). Low-scoring regions indicate structurally favorable environments and roughly correspond to the predicted
-helices and ß-strands (not shown). The program ANOLEA26 was used to calculate pseudoenergy terms for individual residues in our wild-type and M161K models of FOXC1 through its knowledge-based atomic-level mean force potential. When M161 was changed to K161, the default (optimal) rotamer choice was selected to replace the wild-type side chain.
|
| Results |
|---|
|
|
|---|
C; G165R) was detected in patient 1. The transversion results in a glycine to arginine change. Another G to C transversion at codon position 169 (506G>C; R169P) was detected in patient 2. This transversion results in an arginine-to-proline change. Sequencing of the patients DNA and 100 normal control chromosomes confirmed that these mutations are not present in the normal population and are de novo mutations. A previously reported M161K mutation, resulting from a T-to-A mutation at codon position 161 (482T
A; M161K) in two different AR patients3 4 but not previously biochemically analyzed, was also included in our analysis.
|
|
|
BOX) activated luciferase expression at only 1% of wild-type levels. FOXC1 carrying G165R transactivated expression of the luciferase reporter at 36% of wild-type levels. M161K and R169P levels of transactivation were 20% and 23% of wild-type, respectively. These results indicate that although these mutant proteins are able to bind to the FOXC1-binding site in vitro and maintain nuclear localization, all three mutant proteins exhibit reduced levels of gene activation. Therefore, the overriding defect resulting from all three mutations is likely to be an impaired capability to activate transcription of target genes.
Molecular Modeling
A theoretical structure of the FHD of FOXC1 was generated by homology modeling using FOXA3 (Fig. 2A ) as the template. Of the closely related FHDs with known structure, the FOXA3 model has the most complete wing 2 region and also contains identical residues at the positions equivalent to M161, G165, and R169 in FOXC1 (Fig. 2B ). The model structures were submitted to the ANOLEA server which computes a nonlocal energy profile for assessing protein structure.26 Methionine 161 in FOXC1 was mutated to lysine in silico to predict structural defects in the M161K molecule. Fig. 6 illustrates that the M161K mutation may disrupt structurally favorable interactions in an interior region of FOXC1 where T88, C135, F136, and M161 are predicted to converge. Assuming there are no changes in the structure of the protein backbone after mutagenesis, the optimal rotamer for K161s side chain noticeably increased the ANOLEA scores of residues 88, 135, and 161 (Fig. 6C ). This pattern suggests that these residues are normally involved in a highly organized region of side chain packing (Fig. 6A ) such as would occur in a globular proteins hydrophobic core. Disruption of this core-like module to accommodate the M161K amino acid substitution may lead to localized instability or deformation of wing 2. Similar results were obtained with in silico mutagenesis analysis of FOXA3 (Fig. 6C ) that has identical amino acids at the equivalent positions (Fig. 2B ). In silico mutagenesis to create G165R and R169P models did not result in any significant changes to ANOLEA scores of nonlocal residues (not shown) as expected for solvent-exposed positions.
|
Nuclear Localization Disruption by Mutations G165R, M161K, or R169P
The nuclear localization signals (NLSs) of FOXC1 reside in the N and C terminus of the FHD.6 Although G165, M161, and R169 all reside in the C-terminal region required for nuclear localization, substitutions at these sites are not sufficient to disrupt the signal and perturb nuclear localization (Fig. 4 ). Therefore, the wing 2 mutations are tolerated in the mutant FOXC1 proteins and do not disturb the NLSs.
Effect of G165R on Wild-Type Levels of DNA Binding and Transactivation Ability
In the FOXA3/DNA crystal structure, position G165 faces outward and is not a predicted base, phosphate backbone, or water contact.1 Therefore, mutation at position 165 is not predicted to disrupt DNA binding. The fact that FOXC1 carrying G165R can maintain wild-type levels of DNA binding (Fig. 5B ) agrees with this prediction and indicates that the winged-helix structure necessary for DNA binding is not disrupted. The severe reduction in transactivation due to the G165R mutation implies a functional role for this position in transactivation. Previous characterization of two FOXC1 disease-causing missense mutations, F112S and I126M in the FHD of FOXC1 show similar disruptions in transactivation while maintaining normal DNA binding.6 Because of the disruptions of proteinprotein interactions, the G165R mutation may also disrupt the intramolecular interactions necessary for transcription activation. Previous reports suggest mutations F112S and I126M do not disrupt DNA binding because they are not involved in protein/DNA interactions.5 Rather, they may disrupt transactivation due to disturbances in proteinprotein interactions and/or intramolecular interactions.5 Previous characterization of the FOXC1 FHD has implicated the N terminus of
-helix 1 in transactivation and DNA binding.5 Wing 2 of FOXA3 meanders across the surface of the three-helix bundle in the vicinity of the N terminus of
-helix 1,1 establishing an environment receptive to interaction between regions of
-helix 1 and wing 2. Therefore we hypothesize that potential intramolecular interactions between wing 2 and
-helix 1 are disrupted by the G165R mutation, indirectly disrupting the transactivation ability while maintaining DNA-binding ability.
Effect of M161 on DNA-Binding and Transactivation Abilities
Stable formation of wing 2 is essential for DNA binding of FOX proteins.1 In our predicted structure of FOXC1, M161 is spatially arranged such that its side chain resides in close proximity to the side chains of T88, C135, and F136 in a region of reduced solvent accessibility (Fig. 6A ), where the cusp of wing 1 converges with
-helix 1 and the middle of wing 2. The substitution of polar lysine for nonpolar methionine at position 161, as occurs in both patients with the reported mutation, places a larger positively charged residue in the vicinity of a hydrophobic-corelike module. In silico mutagenesis of M161 to M161K, in which only the atomic constituents and rotation of the amino acid side chain at position 161 were altered, suggests that the new side chain would no longer closely associate with the T88/C135/F136 cluster (Fig. 6B ). We hypothesize that M161K perturbs the stability of wing 2 enough to reduce DNA-binding ability by affecting the interatomic nonlocal interactions between residues T88, C135, and M161.
Unlike G165R the loss of transactivation in M161K is most likely caused by the approximate 20% to 30% reduction in DNA binding of this mutation (Fig. 5B ). M161K is similar to a previous mutation, L86F, in
-helix 1 of the FHD of FOXC1.7 L86F shows a minor reduction of DNA binding ability but shows marked reduction in transactivation from a luciferase reporter,7 implicating L86F as having an important role in intramolecular interactions.22 In FOXA3, position M202 is one of 35 known residues that contribute to the hydrophobic core of the FHD.1 This position is conserved in FOXC1 and corresponds to position M161. We predict therefore that M161 may also have a role in intramolecular interactions that the M161K mutation disrupts, impairing transactivation.
Effect of Mutation at Position R169 on DNA Binding and Transactivation Ability
Based on the crystal structure of FOXA3,1 the R169P mutation occurs in a position that is predicted to be responsible for the sole side chainbase minor groove contact in wing.1 Our analysis predicts that R169P disrupts this base contact and potentially destabilizes the formation of wing 2, which is critical for DNA binding activity.1
The severe reduction in DNA binding from mutation R169P is translated into a marked reduction in transactivation of a luciferase reporter (Fig. 5C ). The severe reduction in DNA binding seen in FOXC1 R169P, and the knowledge that R169 makes a DNA contact in the related FOX protein FOXA3,1 supports the hypothesis that this reduced transactivation is a result of reduced DNA-binding capability.
The Role of Wing 2 in DNA-Binding and Transactivation of FOXC1
FOX proteins exhibit similar binding specificity to a core sequence (TAAAYA) due to the conservation of residues in the recognition helix1 ; however, it has been suggested that the DNA binding specificity may be determined by residues outside of
-helix 3.28 In sequence comparisons of FOX genes, residues of wings 1 and 2 have the most divergent sequences, suggesting their possible involvement in sequence recognition.10 29 Although mutations M161K and R169P reduced the DNA binding ability of FOXC1, none of the three wing 2 mutations discussed herein affected the binding specificity of FOXC1 (data not shown). It cannot be ruled out, however, that other residues in wing 2 play a role in sequence specificity by interacting with other residues in the helix bundle. Nevertheless, mutations G165R, M161K, and R169P, although apparently disrupting the stable formation of wing 2 and subsequent DNA binding, do not appear to disturb internal interactions between wing 2 and other regions involved in DNA binding specificity.
We propose that wing 2 plays an important role in DNA binding and the transactivation capability of FOXC1. Previous structural analyses of FOX proteins strongly support a role for wing 2 in DNA binding.1 11 30 Disrupting R169, the sole base contact made in the minor groove by wing 2, is sufficient to reduce FOXC1 DNA binding to 20% to 30% of wild-type levels and reduce the capability of the mutated protein to transactivate genes nearly fourfold (Fig. 5C ). Disrupting residue M161, which is not predicted to be involved in any direct base contacts,1 is also sufficient to reduce the DNA binding and transactivation ability of FOXC1 (Fig. 5B ). Clearly, wing 2 of the FOXC1 FHD plays an important role in DNA binding, whether it is through a direct base contact or through intramolecular interactions that may be required for correct conformation of the FOXC1 protein. Our analysis of the G165R mutation indicates that wing 2 also plays a key role in transactivation, a role not initially predicted outside of wing 2 involvement in DNA binding. Future studies of the functional role of wing 2 will not only enhance our understanding of the FOX/wHTH family, but will also expand our knowledge of what our results indicate is a key region necessary for FOX protein DNA binding and transactivation.
| Acknowledgements |
|---|
| Footnotes |
|---|
Submitted for publication February 17, 2004; revised March 19, 2004; accepted March 29, 2004.
Disclosure: T.C. Murphy, None; R.A. Saleem, None; T. Footz, None; R. Ritch, None; B. McGillivray, None; M.A. Walter, None
The publication costs of this article were defrayed in part by page charge payment. This article must therefore be marked "advertisement" in accordance with 18 U.S.C.
1734 solely to indicate this fact.
Corresponding author: Michael A. Walter, 832 Medical Sciences Building, University of Alberta, Edmonton AB, Canada T6G 2H7; mwalter{at}ualberta.ca.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. H. Strungaru, I. Dinu, and M. A. Walter Genotype-Phenotype Correlations in Axenfeld-Rieger Malformation and Glaucoma Patients with FOXC1 and PITX2 Mutations Invest. Ophthalmol. Vis. Sci., January 1, 2007; 48(1): 228 - 237. [Abstract] [Full Text] [PDF] |
||||
![]() |
Y. A. Ito, T. K. Footz, T. C. Murphy, W. Courtens, and M. A. Walter Analyses of a Novel L130F Missense Mutation in FOXC1 Arch Ophthalmol, January 1, 2007; 125(1): 128 - 135. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. Weisschuh, P. Dressler, F. Schuettauf, C. Wolf, B. Wissinger, and E. Gramer Novel Mutations of FOXC1 and PITX2 in Patients with Axenfeld-Rieger Malformations. Invest. Ophthalmol. Vis. Sci., September 1, 2006; 47(9): 3846 - 3852. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. Cella, J. P. Cabral de Vasconcellos, M. Barbosa de Melo, B. Kneipp, F. F. Costa, C. A. Longui, and V. P. Costa Structural Assessment of PITX2, FOXC1, CYP1B1, and GJA1 Genes in Patients with Axenfeld-Rieger Syndrome with Developmental Glaucoma Invest. Ophthalmol. Vis. Sci., May 1, 2006; 47(5): 1803 - 1809. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. B. Berry, M. A. Lines, J. M. Oas, T. Footz, D. A. Underhill, P. J. Gage, and M. A. Walter Functional interactions between FOXC1 and PITX2 underlie the sensitivity to FOXC1 gene dose in Axenfeld-Rieger syndrome and anterior segment dysgenesis Hum. Mol. Genet., March 15, 2006; 15(6): 905 - 919. [Abstract] [Full Text] [PDF] |
||||
![]() |
F. B. Berry, Y. Tamimi, M. V. Carle, O. J. Lehmann, and M. A. Walter The establishment of a predictive mutational model of the forkhead domain through the analyses of FOXC2 missense mutations identified in patients with hereditary lymphedema with distichiasis Hum. Mol. Genet., September 15, 2005; 14(18): 2619 - 2627. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. A. Saleem, S. Banerjee-Basu, T. C. Murphy, A. Baxevanis, and M. A. Walter Essential structural and functional determinants within the forkhead domain of FOXC1 Nucleic Acids Res., August 6, 2004; 32(14): 4182 - 4193. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |