Basic Region Variants of the MAX b-HLH-LZ preferentially form heterodimers with the MYC b-HLH-LZ to bind the E-box, rather than binding as homodimers.

The MYC associated factor X (MAX) is the heterodimeric partner of the MYC paralogs (MYC, MYCN and MYCL). When deregulated, high level of the MYC paralogs contribute to all aspects of tumorigenesis and tumor growth. MAX can also heterodimerize with the MXD proteins, MNT and MGA. Heterodimerization and sequence specific DNA binding to the E-Box sequences at gene promoters is controlled by their heterodimerization with the MAX b-HLH-LZ. As a heterodimer with MAX, MYC proteins activate genes involved in cell metabolism, growth and proliferation whereas MXD proteins, MNT and MGA repress them. MAX can also bind to the E-Bos sequence as a homodimer. Being devoid of a transactivation domain it can act as an antagonist of the MYC/MAX heterodimers. Variants of MAX have been reported to be linked to cancer. These variants are either not expressed, inactivated or lead to missense mutations. This has led to the notion that MAX may have a tumor suppressor role. Here, we characterize three of those variants with missense mutations in the basic region, i.e. E32K, R35P and R35C. We analyzed their heterodimerization with the b-HLH-LZ of MYC and their DNA binding properties as homo-and heterodimers. The R35C variant b-HLH-LZ was found to have a markedly increased affinity for the b-HLH-LZ of MYC. We also observed that all three b-HLH-LZ variants have a lower affinity as homodimers for the E-Box than the WT. This was shown to lead to a preferential binding of all the heterodimeric b-LHLH-LZ to the E-Box. This effect is exacerbated in the case of the R35C variant. We argue that this preferential binding of MYC as heterodimers with these variants to E-Box sequences could contribute to tumorigenesis. Hence, our results suggest that, mechanistically, the MAX homodimer bound to the E-Box could act as a tumor suppressor. MATERIALS AND METHODSO_ST_ABSMolecular modelingC_ST_ABSThe open source version 1.7.6.0 of Pymol was used for modeling and molecular rendering [1]. The crystal structure of the MAX homodimer bound to the E-Box (1HLO [2]) was used as a template for the generation of the models. The variants were generated using the mutagenesis function in the wizard. The conformation of the K32 side chain was manually set in order to avoid introducing steric clashes with DNA. Protein expression and purificationThe cDNA, coding for the MAX b-HLH-LZ (Max* hereafter, residues 22-103, UniProt entry P61244-1) to which are added the GSGC residues in c-terminal, inserted in the pET3a vector was already available in the laboratory [3] and was used as a template to generate the plasmids with inserts coding for each of the mutants (E32K, R35C and R35P) through quick-change PCR with Q5 DNA polymerase and DpnI from New England Biolabs. The primers used were purchased from IDT DNA, their sequences are listed in Table S1. Sequence for each construct was confirmed by Sanger sequencing at the Plateforme de sequencage SANGER - Centre de recherche du CHU de Quebec - Universite Laval. The primary structure for the basic region of each construct is given in Fig. 2A. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=137 SRC="FIGDIR/small/715400v1_fig2.gif" ALT="Figure 2"> View larger version (41K): org.highwire.dtl.DTLVardef@1b05d5eorg.highwire.dtl.DTLVardef@1c1d692org.highwire.dtl.DTLVardef@ee469dorg.highwire.dtl.DTLVardef@15e0ba4_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOFigure 2.C_FLOATNO Structure schematics, specific and non-specific interactions dictating specificity and stability of binding of the basic region of MAX to the canonical (CACGTG) E-Box. A. Primary structure for the basic region of MAX and each of the variants. Positions making the most important contacts with the E-box are indicated by black arrows. Positions for the variants studied here are colored according to the Zappo colour scheme, following their physico-chemical properties: red for negative, blue for positive, magenta for proline and yellow for cysteine. B. The side chain (carboxylate) of E32 receives H-Bonds from the CA nucleobases in the leading strand (white carbon atoms). R35 and R36 make a salt bridges with phosphate groups while and the guanidino moiety of R36 makes a specific H-Bond with the nucleobase of the G in the strand of the reverse complement (cyan carbon atoms). C. The R35C mutation removes one non-specific salt-bridge at the interface of the complex. D. The aliphatic portion of the K side chain in the E32K variant is unable to accept the H-Bonds from the CA nucleobases and leads to the stabilisation of the complex and the helical structure of the basic region. E. In addition to removing a salt-bride, the Pro residue in the R35P kinks the path of the basic region, prevents the establishment of the specific H-Bonds mandatory for recognition of the E-Box and leads to unfolding of the helical state. C_FIG The MYC b-HLH-LZ (Myc*), the Max*WT b-HLH-LZ and its variants were expressed and purified as previously described [3,4] After lyophilisation, the b-HLH-LZs were kept at -20{degrees}C and solubilised in Myc buffer (50 mM NaCl, 50 mM NaH2PO4 pH 5.5) for Myc* or PBS for Max* at a final concentration of 1 mM before use. Circular dichroismAll circular dichroism (CD) measurements were performed on a Jasco J-810 spectropolarimeter equipped with a Peltier-type thermostat. The instrument was routinely calibrated using an aqueous solution of d-10-(+)-camphorsulfonic acid at 290.5 nm. Samples were prepared as follows: Max* (either WT or a variant) was diluted in 100 {micro}l 2X CD buffer (40 mM KCl, 11.4 mM K2HPO4, 28.6 mM KH2PO4, pH 6.8) and the volume adjusted to 106 {micro}l with PBS. 10 {micro}l TCEP 16 mM were added, and the volume further adjusted to 192 {micro}l with ddH2O before samples were incubated overnight at room temperature. After reduction, Myc* was added and the volume adjusted to 198 {micro}l with Myc buffer (Na2HPO4 0.95 mM, NaH2PO4 49.05 mM, 50 mM NaCl, pH 5.5). The DNA complexes were prepared as follows. After a 10 minutes incubation of the protein samples at room temperature, 0, 1 or 2 {micro}l of 2 mM of specific or non-specific DNA duplexes in 10 mM Tris pH 8.0 were added and the volume adjusted to 200 {micro}l with 10 mM Tris pH 8.0. The strands of the specific probe were: 5-ATT ACC CAC GTG TCC T*AC-3 and 5-GTA GGA CAC GTG GGT* AAT-3 (with the E-box sequence underlined) and the non-specific probe: 5-ATT ACC TCC GGA TCC T*AC-3 and 5-GTA GGA TCC GGA GGT* AAT-3 (Integrated DNA Technologies). Samples were further incubated for 10 minutes at room temperature and transferred to a 1 mm path length quartz cuvette. All spectra were recorded from 250 to 195 nm at 0.1 nm intervals by accumulating 10 spectra at 25 {degrees}C. Thermal denaturations were recorded at 222 nm from 5 to 95 {degrees}C at a heating rate of 1 {degrees}C/min. CD signal for spectra and thermal denaturations was corrected by substracting the signal from corresponding spectra or thermal denaturation either for buffer alone or the appropriate DNA duplex. CD signal was then converted to mean residue ellipticity using the following formula [5]: [{theta}] = {delta} {middle dot} MRW/(10{middle dot}c l) where [{theta}] is the mean residue ellipticity in deg {middle dot} cm2 dmol-1, {delta} is the CD signal in millidegrees, MRW is the mean residue weight, c is the concentration in mg/ml and l is the pathlength in mm. For the heterodimers, the concentration used was the sum of Max* and Myc* and the MRW was determined using a weighted average.