©1996 by The American Society for Biochemistry and Molecular Biology, Inc.
Preferential Binding of MyoD-E12 versus Myogenin-E12 to the Murine Sarcoma Virus Enhancer in Vitro(*)

(Received for publication, December 26, 1995)

Piotr J. Czernik (1) Charlotte A. Peterson (1) (2) (3) Barry K. Hurlburt (1)(§)

From the  (1)Department of Biochemistry and Molecular Biology and the (2)Department of Internal Medicine, University of Arkansas for Medical Sciences, and the (3)John McClellan Veterans Administration Hospital, Little Rock, Arkansas 72205

ABSTRACT
INTRODUCTION
MATERIALS AND METHODS
RESULTS
DISCUSSION
FOOTNOTES
ACKNOWLEDGEMENTS
REFERENCES

ABSTRACT

The MyoD family of transcription factors regulates muscle-specific gene expression in vertebrates. In the adult rat, MyoD mRNA accumulates predominately in fast-twitch muscle, in particular type IIb and/or IIx fibers, whereas Myogenin mRNA is restricted to slow-twitch type I muscle fibers. Transgenic mice expressing the avian v-ski oncogene from the murine sarcoma virus (MSV) promoter-enhancer display preferential hypertrophy of type IIb fast-twitch muscle apparently because of the restricted expression of the transgene. We tested the hypothesis that preferential interactions of MyoD, as a heterodimer with E12, with the MSV enhancer, which has six E-box targets for MyoD family proteins, could contribute to v-ski gene expression in IIb muscle fibers. A series of quantitative binding studies was performed using an electrophoretic mobility shift assay to test MyoD-E12 versus Myogenin-E12 binding to the MSV enhancer. Our results indicate that MyoD-E12 binds the MSV enhancer with higher affinity and higher cooperativity than Myogenin-E12. Interestingly, MyoD-E12 bound all of the individual E-boxes tested with positive cooperativity indicating DNA-mediated dimerization of the protein subunits.


INTRODUCTION

Activation of muscle-specific gene expression in vertebrates is controlled by the MyoD family of transcription factors, MyoD, Myogenin, MRF4, and Myf-5, which are exclusively expressed in skeletal muscle (1, 2, 3, 4, 5) . Ectopic expression of any one member of this family in non-muscle cells, such as 10T1/2 or NIH3T3 fibroblasts, results in phenotypic conversion of that cell to muscle, indicating the trans-dominant nature of these proteins(3, 6) . The four family members share 80% homology in a 60-amino acid region known as the basic helix-loop-helix (bHLH) (^1)domain (7) . The basic region is responsible for DNA binding activity (8) and is involved in transactivation of regulated genes, whereas the helix-loop-helix motif constitutes an interface for dimerization with other bHLH proteins(9, 10, 11) . MyoD family proteins are not thought to form efficient homodimers; rather they readily heterodimerize with other HLH proteins, including the E12 and E47 products of the E2A gene (12) . Heterodimers are capable of activating transcription of muscle-specific genes as a result of binding to the E-box consensus sequence (CANNTG) in muscle-specific enhancers, presumably facilitating the formation of transcription-competent RNA polymerase-promoter complexes. The roles of these transcription factors in myogenesis has been reviewed recently(13) .

Distinctions between these seemingly redundant transcription factors have been elucidated by developmental expression profiles and null mutant studies from which a general model has been proposed(14) . Whereas myoD or myf-5 null mice have fairly normal muscle development, the double mutant is lethal, as are myogenin null mice(15, 16, 17, 18, 19) . There is a temporal pattern of expression for the MyoD family members in the mouse: Myf-5 is expressed first, followed shortly after by Myogenin, MyoD, and MRF4. In the neonate, all four proteins are expressed at a time when the fetal isoforms of myosin heavy chain (MyHC), rather than the adult isoforms predominate. Following birth and maturation of innervation, Myf-5 expression is reduced to an undetectable level, whereas MRF4 expression remains relatively high. Interestingly, both MyoD and Myogenin expression is also reduced, however, to different extents in different fiber types. MyoD expression predominates over that of Myogenin in fast-twitch type IIb/IIx muscle fibers, whereas Myogenin is preferentially expressed in slow-twitch type I fibers(20, 21) . During this time, the fetal isoform of MyHC is replaced by the adult isoforms. The fiber type-specific expression of MyoD and Myogenin is responsive to hormonal and neural signals. Denervation of either fiber type results in increased expression of both MyoD and Myogenin, as if the cells are reverting to the fetal form; adult MyHC is reduced, and fetal MyHC is reexpressed. Cross-reinnervation of a slow muscle (soleus) with a fast motor neuron or thyroid hormone treatment, both of which induce transformation to fast-twitch fiber types, increased the level of MyoD expression and lowered Myogenin expression in the affected myofibers (20) . These results suggest that the bHLH factors regulate different sets of genes in fetal and adult muscle. In the adult, they may be involved in controlling fiber type-specific gene expression in response to external signals.

The viral ski gene (v-ski) encodes a nuclear oncoprotein capable of causing morphological transformation, allowing cell growth in soft agar and inducing myogenesis in cultured avian cells(22, 23, 24) . The avian c-ski proto-oncogene encodes a nuclear protein of unknown function which can bind DNA in the presence of other nuclear proteins(25) . However, when expressed at high levels, similar to that of v-ski during infection, c-ski acquires the oncogenic and myogenic activities of v-ski(26) . This result led Colmenares et al. (26) to conclude that differences in the level of expression between v-ski and c-ski rather than differences in the encoded proteins per se was sufficient to account for the disparate characteristics of the two evolutionarily related proteins.

c-ski expression driven by the promoter-enhancer region of the murine sarcoma virus (MSV) long terminal repeat has unexpected effects in transgenic mice. In three independent lines of mice, hypertrophy of skeletal muscle was observed. One of these lines(8566) was studied in detail, and the hypertrophy of muscle was found to be the result of specific hypertrophy of type IIb fibers, the fastest of the fast-twitch fibers(27) . c-ski expression and the hypertrophic phenotype become elevated at approximately 12 days postpartum, the time at which innervation matures, and the adult isoforms of contractile proteins accumulate. c-ski mRNA accumulation in MSV-ski transgenic mice was reduced in denervated hindlimb muscle or in response to chemically induced hypothyroidism(27) .

The obvious correlation between the expression patterns of the endogenous myoD gene and the MSV-v-ski transgene suggests the possibility that the MSV promoter-enhancer may be responding to MyoD activity specifically in adult IIb muscle fibers. In this case, the MSV enhancer would be a useful tool, not only in studies of differences in MyoD and Myogenin activity, but also for analyses of changes in MyoD activity during development. The idea that MyoD in type IIb fibers may regulate the MSV enhancer is supported by the finding that the bHLH transcription factor E47 can bind and transactivate the murine leukemia virus promoter-enhancer in vivo(28) . The murine leukemia virus and MSV enhancers, as well as those from all type C retroviruses, are highly homologous(29) . Based on these facts, we postulated that the tissue-specific expression of v-ski in the transgenic mice may result from the action of muscle-specific regulators available in fast fibers but absent in slow fibers, of which MyoD is the most likely candidate. As a first test of this hypothesis, we asked if there were intrinsic differences in the ability of MyoD (fast-specific) and Myogenin (slow-specific) to bind the MSV enhancer. A quantitative, in vitro protein-DNA binding analysis of the MSV enhancer region with MyoD-E12 and Myogenin-E12 heterodimers was performed. Our results indicate that MyoD-E12 has higher affinity for the MSV enhancer and that it binds with higher cooperativity than Myogenin-E12. Taken together our results indicate that MyoD is likely to contribute to the fast fiber-specific expression of c-ski in transgenic mice.


MATERIALS AND METHODS

Preparation of DNA Fragments

The MSV enhancer region contains six E-boxes that we refer to as E1, E2, etc. from 5` to 3` (Fig. 1, A and B). Four overlapping MSV enhancer fragments were amplified using Vent DNA polymerase (New England Biolabs) with pMSVski plasmid as a template and synthetic primers. pMSVski was the generous gift of Stephen H. Hughes (NCI-FCRDC). The resulting DNA fragments contained sequences from 11-76, 11-152, 51-276, and 127-276 (see Fig. 1). These fragments were cloned into SmaI site of the pUC118 polylinker, generating pPC21, pPC31, pPC41, and pPC52, respectively. The nucleotide sequences of the cloned fragments were confirmed. For electrophoretic mobility shift assay (EMSA), DNA fragments were amplified with different combinations of P end-labeled ``reverse'' and ``-20'' universal primers complementary to sequences in the pUC vector and primers complementary to the MSV enhancer (Fig. 1B). DNA fragments E23 and E6, which contain MSV enhancer sequences from 86-124, and 212-241, respectively, were synthesized. Fragment E2 was created with synthetic oligonucleotides for E23 fragment where E3-box was mutated (described below). Other DNA fragments with mutated E-boxes were amplified using mutagenic primers and either the reverse or -20 primers. Mutagenic primers were designed to introduce T C and G A base transitions in the targeted E-box sequences (CANNTG CANNCA). These mutations were effective in abolishing MyoD-E12 and Myogenin-E12 binding to DNA fragments containing a single mutated E-box (data not shown). Fragments E2345 and E45 were generated by BfaI restriction enzyme digestion of E23456 and E456 fragments, respectively. Control DNA that contained the pUC118 polylinker region amplified with reverse and -20 primers was used as a negative control in EMSA. Also, a control DNA fragment with an E-box positioned near the end of the DNA fragment showed no differences in binding affinity for either protein compared with a similar DNA fragment with the E-box more centrally located (data not shown). All synthetic DNAs were made with a Milligen Expedite DNA Synthesizer by the UAMS Molecular Biology Core Facility and were end labeled with [-P]ATP and T4 polynucleotide kinase prior to polymerase chain reaction or annealing. Radioactive DNA substrates for EMSA were gel purified, and the concentrations were determined by absorbance at 260 nm.


Figure 1: Sequence of the MSV enhancer region and DNA fragments used for EMSA. Panel A, the sequence of the MSV enhancer with the six E-boxes enclosed. Binding sites for other transcription factors are shown below the horizontal bars and labeled. MCREF-1, LVa; GR, glucocorticoid receptor; NF-1, nuclear factor 1; CAT, CCAAT box-binding protein. Panel B, the DNA fragments tested in EMSA.



MyoD, Myogenin, and E12 Proteins

Full-length mouse MyoD, Myogenin, and E12 proteins were overexpressed in Escherichia coli strain BL21(DELysS) from T7 expression plasmids. Plasmids pT5-MyoD and pT5-E12 were the generous gifts of Charles P. Emerson, Jr. (University of Pennsylvania School of Medicine). Purified Myogenin was the generous gift of Soheila J. Maleki (University of Arkansas Medical Sciences). The expression and purification of MyoD, Myogenin, and E12 will be described elsewhere. (^2)The three proteins were judged to be approximately 90-95% pure following sodium dodecyl sulfate-polyacrylamide gel electrophoresis and Coomassie Brilliant Blue staining. Small aliquots of purified protein were stored at -20 °C and thawed freshly for each experiment. The concentrations of MyoD, Myogenin, and E12 were determined by absorbance at 280 nm using extinction coefficients of 17,750, 21,980 and 17,900 M, respectively.

EMSA

Proteins for heterodimer experiments were diluted from frozen stocks to 800 nM on ice in MCK buffer (10 mM Hepes/KOH, pH 7.9, 80 mM KCl, 1 mM EDTA, 1 mM dithiothreitol, 0.05% Triton X-100, 5% glycerol). Equal volumes of MyoD or Myogenin were mixed with E12, and the resulting solutions were equilibrated for 20 min at room temperature. Protein samples were prepared in 2-fold dilutions from 400 nM to 50 pM and mixed with equal volumes of DNA substrates diluted to 20 pM in MCK buffer. Thus, binding reactions contained 10 pM DNA and proteins in concentrations ranging from 25 pM to 200 nM, with respect to the heterodimer. Binding reactions were allowed to equilibrate for at least 30 min prior to electrophoresis. Native, 1-mm thick, 6% acrylamide gels (50:1 acrylamide:bisacrylamide) in 1 times TBE (89 mM Tris borate, 2 mM EDTA) were prepared in Protean IIxi (Bio-Rad) cooled apparatus at least 1 h prior to electrophoresis and prerun at 200 V constant voltage to reach a constant current of approximately 18 mA. Gel temperature was maintained at 17 °C with a water bath circulator at all times. Binding reactions were loaded on running gels, and then the voltage was raised to 250 V. Resolved gels were fixed in 10% acetic acid, 7% methanol for 15 min and dried. Gels were exposed to Reflection x-ray film (DuPont) for 12 h and 72 h at -70 °C with intensifying screens. During development of the EMSA we tested several different concentrations of acrylamide and various temperatures. Higher concentrations of acrylamide decreased the absolute value and compressed the range of the apparent binding constants presumably because of ``caging effects''(30) . We also tested temperatures ranging from 17 to 24 °C for EMSA. DNA complexes with MyoD-E12 were much more stable at the elevated temperatures than Myogenin-E12; however, both complexes were quite stable in the EMSA at 17 °C (data not shown).

Analysis of Data

Undersaturated films from the 12-h exposure were developed, scanned by an Image Densitometer (Bio-Rad), and the relative optical densities of free DNA bands were measured. Apparent amounts of bound DNA were calculated and plotted against respective protein concentrations. The resulting data were fitted with a sigmoidal (power logistic) function using Origin (MicroCal).

In this function, A1 is the initial y axis level, A2 is the final y axis level, x(0) is the center point of the curve, and p is the rate of increase or slope of the curve. Our analysis of theoretical binding curves for single and multiple protein/DNA interactions showed that the variable p exponent existing in the sigmoidal fitting function can be used as a relative measure of binding cooperativity. A p value of 1.0 indicates no cooperativity; p > 1 or p <1 indicates positive or negative cooperativity, respectively; and the magnitude of p is indicative of the degree of cooperativity. It is important to note that the p value is a relative value describing the slope of the fitted curve and is not directly related to the Hill coefficient.

For the purposes of this work, we define the binding constant (K) as the concentration of protein required to bind 50% of the total concentration of DNA. We use the term K rather than K(D) because we are considering multiple equilibria where either MyoD or Myogenin is interacting with another protein, E12, and binding to multiple, nonequivalent target sites. Considering a single binding site, there are at least two formal types of equilibrium which could be occurring. The MyoD family of proteins is commonly thought to bind as dimers, in which case, there is a monomer-dimer equilibrium, for example MyoX (MyoD or Myogenin) dimerizing with E12, described as follows.

where

There is also a protein-DNA equilibrium, for example MyoX-E12 interacting with DNA, described as follows.

where

Since there are no data currently to describe the protein/protein interactions, we are unable to assign accurate values to the dissociation constant for the protein/DNA interaction.

Alternatively, the proteins may bind the target DNA successively as monomers, as described.

where

In this case, the units of K(D) would be M^2 and cannot be derived directly from our data.


RESULTS

EMSA Analysis of MSV Enhancer

The MSV enhancer region under study contains six E-box sequences, which were numbered from 1 to 6 starting from the 5` end (E1, E2, etc; Fig. 1). DNA binding sites for six different nuclear factors were previously identified in the same region (Fig. 1A). To examine MyoD-E12 and Myogenin-E12 binding properties to DNA targets containing multiple E-box motifs, we created sets of nested deletion DNA fragments. The DNA fragment containing all six E-boxes (E123456) was too large for EMSA. Therefore, the first set contained DNA fragments where E-boxes were successively deleted from the 5` end of the E23456 fragment. A second set of nested deletions was generated from the 3` end of fragment E23456. Another 3` set of deletions was analyzed on parent DNA fragment E123 to test the contributions of E1. Using the EMSA described in the previous section, we measured apparent equilibrium binding constants (K) for protein-DNA complexes formed with MyoD-E12 or Myogenin-E12 heterodimers and MSV enhancer fragments.

Analysis of 5` Deletions of E23456

The results of our saturation experiments for 5` deletions of fragment E23456 in which a limiting amount of DNA (10 pM) was incubated with increasing concentrations of MyoD-E12 or Myogenin-E12 (25 pM to 200 nM) are shown in Fig. 2, C and D. A representative autoradiogram for MyoD-E12 interactions with the E23456 DNA fragment is presented in Fig. 2A and for Myogenin-E12 interactions with the E3456 DNA fragment in Fig. 2B. The unbound DNA signal was quantitated by densitometry, and the percentage of bound DNA was plotted versus the concentration of protein used (Fig. 2, C and D). For all saturation experiments, the concentration of protein required to shift 50% of the DNA was taken as the K.


Figure 2: EMSA for interaction of MyoD-E12 and Myogenin-E12 with 5` deletions of the E23456 DNA fragment. MyoD-E12 or Myogenin-E12 and P-labeled DNA fragments (10 pM) were equilibrated and electrophoresed through native acrylamide gels as described under ``Materials and Methods.'' Panel A, representative autoradiogram of MyoD-E12 and the E23456 DNA fragment. Panel B, representative autoradiogram of Myogenin-E12 and the E3456 DNA fragment. Panels C and D, the amount of unbound DNA in each lane was quantitated by densitometry and used to calculate the concentration of bound DNA. The percentage of bound DNA is plotted versus the concentration of protein (calculated as the heterodimer). Arrows indicate resolved protein-DNA complexes. Panel C, MyoD-E12; panel D, Myogenin-E12. , E23456; bullet, E3456; , E456; , E56; , E6.



The apparent binding constants for the parent E23456 fragment in this set were 0.2 nM and 0.3 nM for MyoD-E12 and Myogenin-E12, respectively. The K values for MyoD-E12 increased 16-fold from 0.2 nM for the E23456 fragment to 3.2 nM for the E6 fragment (Table 1). There were only very slight changes in K when E3 or E5 was deleted from the fragments, indicating that these E-boxes do not contribute significantly to binding affinity by MyoD-E12. The changes in the p value in Table 1for MyoD-E12 do not correlate with changes in affinity. Although removal of E2 from the target DNA increases the K 4-fold, there is no significant change in cooperativity. However, deletions of E3 and E4 reduce the positive cooperativity of binding.



The 15-fold increase in K for Myogenin-E12, as a result of deleting all but E2 from the parent E23456 fragment, is similar to that of MyoD-E12, and the absolute values for K are similar (Table 1). However, the large increase in K (10-fold) caused by deleting E2 indicates that this E-box is perhaps the most important for Myogenin-E12 binding. The presence of E4 in the fragment does not enhance Myogenin-E12 binding as it did for MyoD-E12. The p values determined from the curves in Fig. 2D indicate that there is positive cooperativity in binding for all DNA fragments tested except E6, but less than that observed for MyoD-E12 binding to the same fragments.

Analysis of 3` Deletions of E23456

The results obtained for MyoD-E12 and Myogenin-E12 complexes with 5` deletion DNA fragments showed that the presence of E2 increases the affinity of Myogenin-E12 binding, whereas both E2 and E4 enhance formation of MyoD-E12-DNA complexes. Our experiments also indicated that cooperative binding is important in the formation of complexes with these proteins, but it is most pronounced for MyoD-E12. To test these ideas further, we analyzed MyoD-E12 and Myogenin-E12 binding with 3` deletion fragments of the MSV enhancer using assays performed as described above for the 5` deletion fragments. A representative autoradiogram for MyoD-E12 interactions with the E234 DNA fragment is shown in Fig. 3A and for Myogenin-E12 interactions with the E2345 DNA fragment in Fig. 3B. The data from densitometric analysis are presented in Fig. 3, C and D. The values for K and cooperativity are summarized in Table 2. The parent E23456 fragment in this set was the same as in 5` deletion group. For MyoD-E12, deletion of E6 caused a 12-fold decrease in affinity, and binding was less cooperative. Removal of E5, generating the E234 fragment, almost completely restored binding to the level observed with E23456 (K = 0.3 nM), but without the high level of cooperativity observed in binding the E23456 fragment. Subsequent deletion of E4 resulted in a fragment (E23) that was bound by MyoD-E12 with nearly the same affinity and cooperativity as for fragment E2345. Finally, deletion of E3, creating fragment E2 with only one E-box, lowered binding affinity and cooperativity relative to E23.


Figure 3: EMSA for interaction of MyoD-E12 and Myogenin-E12 with 3` deletions of the E23456 DNA fragment. MyoD-E12 or Myogenin-E12 and P-labeled DNA fragments (10 pM) were equilibrated and electrophoresed through native acrylamide gels as described under ``Materials and Methods.'' Panel A, representative autoradiogram of MyoD-E12 and the E234 DNA fragment. Panel B, representative autoradiogram of Myogenin-E12 and the E2345 DNA fragment. Panels C and D, the amount of unbound DNA in each lane was quantitated by densitometry and used to calculate the concentration of bound DNA. The percentage of bound DNA is plotted versus the concentration of protein (calculated as the heterodimer). Arrows indicate resolved protein-DNA complexes. Panel C, MyoD-E12; panel D, Myogenin-E12. , E23456; bullet, E2345; , E234; , E23; , E2.





Deletion of E6 from the E23456 fragment reduced the affinity of Myogenin-E12 binding very much like it did MyoD-E12 binding with an increase in K of approximately 12-fold but without the concomitant reduction in cooperativity. There was a slight increase in affinity, approximately 2-fold, when E5 was removed, leaving E234. Removal of E4 had devastating effects on binding: the E23 and E2 fragments were bound by Myogenin-E12 with K values of 26 and 23 nM, respectively, and all positive cooperativity of binding was lost.

Analysis of 3` Deletions of E123

Since fragment E123456 was too large for EMSA we tested the contribution of E1 in enhancer binding using E123 as the parent fragment and removed E-boxes from the 3` end. A representative autoradiogram for MyoD-E12 interactions with the E123 DNA fragment is shown in Fig. 4A and for Myogenin-E12 interactions with the E123 DNA fragment in Fig. 4B. The data from the densitometric analysis are presented in Fig. 4, C and D, and the derived values for K and cooperativity are summarized in Table 3. MyoD-E12 bound the E123 with a K of 3.7 nM and modest positive cooperativity. Unexpectedly, the E12 fragment was bound with nearly 4-fold higher affinity (K = 1 nM), as if E3 had a negative effect on fragment binding. Removal of E2 from the E12 fragment had no significant difference in affinity. Myogenin-E12 bound the E123 fragment with nearly the same affinity as by MyoD-E12, whereas removal of E3 caused a nearly 5-fold increase in K. Fragment E1 was bound by Myogenin-E12 with the same affinity as E123. These data indicate that the presence of E3 has opposite effects on binding by MyoD-E12 and Myogenin-E12. The results from the deletion of E3 from fragment E23456 do not corroborate the effects of E3 deletion from fragment E123. Thus, there must be contributions to binding from the context of the various E-boxes.


Figure 4: EMSA for interaction of MyoD-E12 and Myogenin-E12 with 3` deletions of the E123 DNA fragment. MyoD-E12 or Myogenin-E12 and P-labeled DNA fragments (10 pM) were equilibrated and electrophoresed through native acrylamide gels as described under ``Materials and Methods.'' Panel A, representative autoradiogram of MyoD-E12 and the E123 DNA fragment. Panel B, representative autoradiogram of Myogenin-E12 and the E123 DNA fragment. Panels C and D, the amount of unbound DNA in each lane was quantitated by densitometry and used to calculate the concentration of bound DNA. The percentage of bound DNA is plotted versus the concentration of protein (calculated as the heterodimer). Arrows indicate resolved protein-DNA complexes. Panel C, MyoD-E12; panel D, Myogenin-E12. , E123; bullet, E12; , E1.





Analysis of Internal Mutations

To clarify the contributions of the individual E-boxes in the MSV enhancer to binding by MyoD-E12 and Myogenin-E12, we tested internal mutations of E2, E3, E4, and E5 using the EMSA. A representative autoradiogram of MyoD-E12 interactions with the E13 DNA fragment is shown in Fig. 5A and of Myogenin-E12 interactions with the E13 DNA fragment in Fig. 5B. The binding profiles are shown in Fig. 5, C and D, and the derived values for and K and p are presented in Table 4. Compared with the E123 fragment, the E13 fragment was bound with slightly higher affinity by MyoD-E12, but only 10% as well by Myogenin-E12. Thus, E2, in this context, has significant and opposite effects on MyoD-E12 and Myogenin-E12 interactions. Mutation of E3 within an E23456 parent fragment caused a small reduction in MyoD-E12 binding affinity and a much larger 15-fold loss in Myogenin-E12 binding. Compared with the E2345 fragment, the E235 fragment was bound with approximately 4-fold higher affinity by MyoD-E12 and with slightly increased affinity by Myogenin-E12. Removal of E5 from the E456 fragment had no significant effects on MyoD-E12 binding but increased Myogenin-E12 binding 4-fold. Taken together, the results from all of the EMSA experiments indicate that contributions of the various E-boxes within the MSV enhancer differ with respect to MyoD-E12 and Myogenin-E12 as well as with the context of the DNA fragment that encompasses them.


Figure 5: EMSA for interaction of MyoD-E12 and Myogenin-E12 with internal E-box mutations of the E23456 DNA fragment. MyoD-E12 or Myogenin-E12 and P-labeled DNA fragments (10 pM) were equilibrated and electrophoresed through native acrylamide gels as described under ``Materials and Methods.'' Panel A, representative autoradiogram of MyoD-E12 and the E13 DNA fragment. Panel B, representative autoradiogram of Myogenin-E12 and the E13 DNA fragment. Panels C and D, the amount of unbound DNA in each lane was quantitated by densitometry and used to calculate the concentration of bound DNA. The percentage of bound DNA is plotted versus the concentration of protein (calculated as the heterodimer). Arrows indicate resolved protein-DNA complexes. Panel C, MyoD-E12; panel D, Myogenin-E12. , E13; bullet, E2456; , E235; , E46.






DISCUSSION

We tested the hypothesis that MSV-v-ski transgene expression in fast-twitch muscle fibers was responsive to MyoD activity by examining the DNA binding characteristics of MyoD-E12 and Myogenin-E12 (slow fiber-specific) for the MSV enhancer in vitro. The results reported here are our initial analyses of this system and indicate complex interactions between these proteins bound to the six E-boxes in terms of affinity and cooperativity. Further analyses will be required to understand fully the molecular mechanisms involved in regulating expression of the MSV enhancer-promoter. Nevertheless, our results clearly indicate that MyoD-E12 bound nearly all of the DNA fragments tested with higher affinity and higher cooperativity than Myogenin-E12. Thus, we conclude that fiber type IIb-specific expression of MSV-ski in transgenic mice likely results in part from activation by MyoD, which preferentially accumulates in these fibers.

MyoD-E12 binding to all the DNA fragments containing single E-boxes (E1, E2, or E6) showed positive cooperativity, whereas only E1 was bound by Myogenin-E12 cooperatively (Table 1Table 2Table 3). In addition, MyoD-E12 bound all and Myogenin-E12 bound most of the DNA fragments containing multiple E-boxes with cooperativity. Thus, there are multiple contributions of cooperativity to DNA binding for these transcription factors. These results are strikingly similar to the interactions of repressor with O(R)1, O(R)2, and O(R)3(31) . repressor dimerization is enhanced by DNA binding, thus there is cooperativity observed in binding a single operator. Repressor also shows cooperativity in binding tandemly arranged target DNAs, via protein/protein interactions between repressor bound to O(R)1 and O(R)2. Binding to O(R)3 is not cooperative, however, because repressor at O(R)2 interacting with repressor at O(R)1 cannot interact with repressor at O(R)3. However, mutations in O(R)1, which disallow repressor binding, promote repressor binding to O(R)2 and O(R)3 with cooperativity. Such high levels of cooperativity dramatically increase the regulatory response due to small changes in protein concentration or activation (31) .

Cooperativity of dimeric protein interactions with a single DNA target site can arise from two distinct mechanisms. The stepwise binding of monomers can show cooperativity by increasing the affinity of binding for the second monomer in the presence of the first bound monomer. For example, this is what is observed for the glucocorticoid receptor DNA binding domain. (^3)Alternatively, cooperativity can also arise for binding of dimeric proteins to a single binding site if protein dimerization is significantly enhanced by DNA binding. It should also be noted that these two mechanisms are not mutually exclusive; both can be occurring. Although our results show cooperativity in binding a single site, they do not help distinguish between whether these proteins bind DNA as dimers, stepwise as monomers, or both. We are currently conducting experiments to clarify the mode of binding utilized by MyoD-E12.

Studies on DNA recognition by bHLH muscle transcription factors showed that the CANNTG E-box sequence establishes a minimal DNA binding site for these proteins. Consensus DNA recognition sequences for MyoD-E12 heterodimer (32) and Myogenin homodimer (and heterodimer with E12 as proposed in (33) ) derived in affinity selection experiments show a high level of similarity. The MSV enhancer contains six E-box sequences. Four of the E-boxes (E2, E4, E5, E6) have the sequence CAGATG, which is different from the consensus sequences determined in selection experiments for either MyoD or Myogenin, but is present in enhancers of muscle-specific promoters (34, 35, 36) . Two of the E-boxes (E1, E3) have the sequence CAGCTG, matching the consensus binding site from selection experiments with Myogenin. The crystal structure of the MyoD bHLH domain-DNA complex (8) revealed protein contacts with the bases of the primary E-box determinants, CA and TG, and some contacts to 5`- and 3`-flanking residues in the double-stranded TCAACAGCTGTTGA DNA, indicating that specific bases outside an E-box may be involved in recognition and binding. It has also been postulated that the central two bases of the E-box motif could contribute to distinct MyoD and Myogenin-DNA binding preferences. Binding site selection experiments showed that Myogenin has a preference for a symmetrical binding site, whereas MyoD binding sequences showed a nonpalindromic distribution of bases in central NN bases of an E-box. In general, our results show only slight differences in affinity for both proteins in combination with E12 binding to single E-box containing DNA fragments (except for E2) and that binding occurs with low to moderate affinity. Thus, MyoD and Myogenin prefer similar DNA targets and bind them with similar affinity.

On the other hand, numerous in vivo experiments have shown that MyoD and Myogenin have distinct functions in the regulation of muscle-specific gene expression. One proposed solution to this apparent contradiction is that heterodimerization, rather than homodimerization, with other bHLH proteins such as ubiquitous E12 and E47 or Id, leads to MyoD- or Myogenin-specific gene activation or inactivation. It has also been shown that other regulatory DNA elements in combination with E-box sequences can differentially contribute to muscle-specific gene activation(37) . These regulatory models, although attractive and well documented, share one weakness: namely, E-box sequences are widespread in genomes (average of one E-box in 256 bases) and are present in the promoters of genes not specifically expressed in muscle. Nevertheless, it has been shown that under certain experimental conditions a single E-box motif can be responsible for Myogenin-dependent gene activation(38) . In our studies, we showed that MyoD-E12 and Myogenin-E12 DNA binding affinity to targets containing multiple E-box sequences is much higher than interactions of these proteins with single DNA fragments containing a single E-box. This indicates that multiple E-box motifs, when present in promoter regions, could contribute to high affinity DNA binding for these proteins and simultaneously provide an additional level to MyoD and Myogenin DNA recognition specificity. Cooperative binding of MyoD-E47 to DNA targets containing multiple E-box motifs (39) and E-box cooperation in the promoter of the rat acetylcholine receptor subunit gene(40, 41) enhance transcriptional activity. We demonstrated that mutations of certain E-box sequences within MSV enhancer significantly reduce affinity of MyoD-E12 and Myogenin-E12 in cooperative DNA binding presumably through lack of communication between proteins. Conversely, the same situation could arise in vivo when two or more transcriptionally active proteins compete for overlapping DNA recognition sites. This competition could result in a delicate balance of repression and activation. Such a situation can be easily predicted in the case of MSV enhancer, where E2 and E4 constitute internal parts of the glucocorticoid response element. Binding of the glucocorticoid receptor in response to steroid hormone treatment could inhibit MyoD-E12 or Myogenin-E12 binding. In fact, clinical application of glucocorticoids in cancer therapy leads to repression of muscle-specific gene expression and a loss of muscle mass(42) . There are similar sequence relationships in other muscle-specific promoters where, in addition to glucocorticoid response element/E-box overlapping motifs, other transcription factor binding sites (e.g. LVa, Ets-1) consistently share the same stretch of bases with E-box sequences.

We have reported the results from an initial probe into the interactions of transcription factors with the MSV enhancer. It is clear from our results that complex protein/protein interactions and protein/DNA interactions are likely to contribute to enhancer activity. Specifically, cooperativity of binding plays a much larger role than previously expected. Currently, we are pursuing experiments that will further address the mechanisms contributing to cooperative binding, the biological relevance of MyoD-E12 and Myogenin-E12 DNA binding within MSV enhancer, and for possible modulation of MSV enhancer activity through glucocorticoid receptor/MyoD, Myogenin competition for their DNA binding sites.


FOOTNOTES

*
This work was supported by grants from the American Heart Association-Arkansas Affiliate (to B. K. H. and C. A. P.), the National Science Foundation Experimental Program to Support Competitive Research (to C. A. P.), and the American Cancer Society (to B. K. H.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore by hereby marked ``advertisement'' in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

§
To whom correspondence should be addressed: Dept. of Biochemistry and Molecular Biology, University of Arkansas for Medical Sciences, 4301 W. Markham St., Little Rock, AR 72205. Tel.: 501-686-5370; Fax: 501-686-8169.

(^1)
The abbreviations used are: bHLH, basic helix-loop-helix; MyHC, myosin heavy chain; MSV, murine sarcoma virus; EMSA, electrophoretic mobility shift assay.

(^2)
S. J. Maleki and B. K. Hurlburt, manuscript in preparation.

(^3)
J. Hill and C. A. Royer, personal communication.


ACKNOWLEDGEMENTS

We are grateful to Simon Hughes, who was involved in the conception of these experiments, and we thank Gary Bannon, Chris Baumann, Paul Gollnick, Kevin Raney and Cathy Royer for critiques on the manuscript. We also thank Cathy Royer for assistance in data analysis and for providing results prior to publication. We are especially grateful to Soheila J. Maleki for providing purified Myogenin, Charles P. Emerson Jr. for providing E12 and MyoD expression plasmids, and Stephen H. Hughes for pMSV-ski.


REFERENCES

  1. Davis, R. L., Weintraub, H., and Lassar, A. B. (1987) Cell 51, 987-1000 [Medline] [Order article via Infotrieve]
  2. Edmondson, D. G., and Olson, E. N. (1989) Genes & Dev. 3, 628-640
  3. Braun, T., Buschhausen-Decker, G., Bober, E., Tannich, E., and Arnold, H. H. (1989) EMBO J. 8, 701-709 [Abstract]
  4. Rhodes, S. J., and Konieczny, S. F. (1989) Genes & Dev. 3, 2050-2061
  5. Wright, W. E., Sassoon, D. A., and Lin, V. K. (1989) Cell 56, 607-617 [Medline] [Order article via Infotrieve]
  6. Yutzey, K. E., Rhodes, S. J., and Konieczny, S. F. (1990) Mol. Cell. Biol. 10, 3934-3944 [Medline] [Order article via Infotrieve]
  7. Murre, C., McCaw, P. S., and Baltimore, D. (1989) Cell 56, 1351-1359
  8. Ma, P. C. M., Rould, M. A., Weintraub, H., and Pabo, C. O. (1994) Cell 77, 451-459 [Medline] [Order article via Infotrieve]
  9. Murre, C., McCaw, P. S., Vaessin, H., Caudy, M., Jan, L. Y., Jan, Y. N., Cabrera, C. V., Buskin, J. N., Hauschka, S. D., Lassar, A. B., Weintraub, H., and Baltimore, D. (1989) Cell 58, 537-544 [Medline] [Order article via Infotrieve]
  10. Shirakata, M., Friedman, F. K., Wei, Q., and Paterson, B. M. (1993) Genes & Dev. 7, 2456-2470
  11. Voronova, A., and Baltimore, D. (1990) Proc. Natl. Acad. Sci. U. S. A. 87, 4722-4726 [Abstract]
  12. Tapscott, S. J., Davis, R. L., Thayer, M. J., Cheng, P. F., Weintraub, H., and Lassar, A. B. (1988) Science 242, 405-411 [Medline] [Order article via Infotrieve]
  13. Dias, P., Dilling, M., and Houghton, P. (1994) Semin. Diag. Pathol. 11, 3-14
  14. Weintraub, H. (1993) Cell 75, 1241-1244 [Medline] [Order article via Infotrieve]
  15. Braun, T., Rudnicki, M. A., Arnold, H. H., and Jaenisch, R. (1992) Cell 71, 369-382 [Medline] [Order article via Infotrieve]
  16. Rudnicki, M. A., Braun, T., Hinuma, S., and Jaenisch, R. (1992) Cell 71, 383-390 [Medline] [Order article via Infotrieve]
  17. Rudnicki, M. A., Schnegelsberg, P. N. J., Stead, R. H., Braun, T., Arnold, H. H., and Jaenisch, R. (1993) Cell 75, 1351-1359 [Medline] [Order article via Infotrieve]
  18. Hasty, P., Bradley, A., Morris, J. H., Edmondson, D. G., Venturi, J., Olson, E. N., and Klein, W. H. (1993) Nature 364, 532-535 [CrossRef][Medline] [Order article via Infotrieve]
  19. Nabeshima, Y., Hanaoka, K., Hayasaka, M., Esumi, E., Li, S., Noaka, I., and Nabeshima, Y. (1993) Nature 364, 501-506 [CrossRef][Medline] [Order article via Infotrieve]
  20. Voytik, S. L., Przyborski, M., Badylak, S. F., and Konieczny, S. F. (1993) Dev. Dynamics 198, 214-224 [Medline] [Order article via Infotrieve]
  21. Petropoulos, C. J., Payne, W., Salter, D. W., and Hughes, S. H. (1992) J. Virol. 66, 3391-3397 [Abstract]
  22. Li, Y., Turck, C. M., Teumer, J. K., and Stavnezer, E. (1986) J. Virol. 57, 1065-1072 [Medline] [Order article via Infotrieve]
  23. Colmenares, C., and Stavnezer, E. (1989) Cell 59, 293-303 [Medline] [Order article via Infotrieve]
  24. Stavnezer, E., Brodeur, D., and Brennan, L. A. (1989) Mol. Cell. Biol. 9, 4038-4045 [Medline] [Order article via Infotrieve]
  25. Nagase, T., Mizuguchi, G., Nomura, N., Ishizaki, R., Ueno, Y., and Ishii, S. (1990) Nucleic Acids Rec. 18, 337-343
  26. Colmenares, C., Sutrave, P., Hughes, S. H., and Stavnezer, E. (1991) J. Virol. 65, 4929-4935 [Medline] [Order article via Infotrieve]
  27. Leferovich, J. M., Lana, D. P., Sutrave, P., Hughes, S. H., and Kelly, A. M. (1995) J. Neurosci. 15, 596-603 [Abstract]
  28. Lawrenz-Smith, S. C., and Thomas, C. Y. (1995) J. Virol. 69, 4142-4148 [Abstract]
  29. Golemis, E. A., Speck, N. A., and Hopkins, N. (1989) J. Virol. 64, 534-542
  30. Lane, D., Prentki, P., and Chandler, M. (1992) Microbiol. Rev. 56, 509-528 [Abstract]
  31. Ptashne, M. (1986) A Genetic Switch , Cell Press and Blackwell Scientific Publications, San Francisco, CA
  32. Blackwell, K. T., and Weintraub, H. (1990) Science 23, 1104-1110
  33. Wright, W. E., Binder, M., and Funk, W. (1991) Mol. Cell. Biol. 11, 4104-4110 [Medline] [Order article via Infotrieve]
  34. Zingg, J. M., Pedraza, G. M., and Jost, J. P. (1991) Nucleic Acids Res. 19, 6433-6439 [Abstract]
  35. Vaidya, T. B., Rhodes, S. J., Moore, J. L., Sherman, D. A., Konieczny, S. F., and Taparowsky, E. J. (1992) Gene (Amst.) 116, 223-230
  36. Yee, S. P., and Rigby, P. W. (1993) Genes & Dev. 7, 1277-1289
  37. Amacher, S. L., Buskin, J. N., and Hauschka, S. D. (1993) Mol. Cell. Biol. 13, 2753-2764 [Abstract]
  38. Edmondson, D. G., Cheng, T. C., Cserjesi, P., Chakraborty, T., and Olson, E. N. (1992) Mol. Cell. Biol. 12, 3665-3677 [Abstract]
  39. Bengal, E., Flores, O., Rangarajan, N., Chen, A., and Weintraub, H. (1994) Proc. Natl. Acad. Sci. U. S. A. 91, 6221-6225 [Abstract]
  40. Berberich, C., Durr, I., Koenen, M., and Witzemann, V. (1993) Eur. J. Biochem. 216, 395-404 [Abstract]
  41. Simon, A. M., and Burden, S. J. (1993) Mol. Cell. Biol. 13, 5133-5140 [Abstract]
  42. Seene, T. (1994) J. Steroid Biochem. Mol. Biol. 50, 1-4 [CrossRef][Medline] [Order article via Infotrieve]

©1996 by The American Society for Biochemistry and Molecular Biology, Inc.