Novel Mutation in Spectrin-like Repeat 1 of Dystrophin Central Domain Causes Protein Misfolding and Mild Becker Muscular Dystrophy*

Background: Knowledge of dystrophin modifications at the protein molecular level is required to understand genotype/phenotype link in muscular dystrophies. Results: A new mutation in DMD gene is accompanied with dystrophin structural and functional modifications. Conclusion: This Becker muscular dystrophy is due to misfolding of the dystrophin rod domain. Significance: We provide new insights in understanding the molecular bases of myopathy at the protein level. Mutations in the dystrophin gene without disruption of the reading frame often lead to Becker muscular dystrophy, but a genotype/phenotype correlation is difficult to establish. Amino acid substitutions may disrupt binding capacities of dystrophin and have a major impact on the functionality of this protein. We have identified two brothers (ages 8 and 10 years) with very mild proximal weakness, recurrent abdominal pain, and moderately elevated serum creatine kinase levels. Gene sequencing revealed a novel mutation in exon 11 of the dystrophin gene (c.1280T>C) leading to a L427P amino acid substitution in repeat 1 of the central rod domain. Immunostaining of skeletal muscle showed weak staining of the dystrophin region encoded by exons 7 and 8 corresponding to the end of the actin-binding domain 1 and the N-terminal part of hinge 1. Spectrofluorescence and circular dichroism analysis of the domain repeat 1-2 (R1-2) revealed partial misfolding of the L427P mutated protein as well as a reduced refolding rate after denaturation. Based on computational homology models of the wild-type and mutated R1-2, a molecular dynamics study showed an alteration in the flexibility of the structure, which also strongly affects the conformational space available in the N-terminal region of the fragment. Our results suggest that this missense mutation hinders the dynamic properties of the entire N-terminal region of dystrophin.

Mutations in the dystrophin gene cause severe Duchenne (DMD) 2 or the milder Becker muscular dystrophy (BMD) (1)(2)(3). The diagnosis of BMD can be made by molecular testing of the dystrophin gene, a muscle biopsy evaluated by immunostaining, clinical presentation, and family history. Although BMD is milder than DMD, the phenotype varies greatly from person to person and can involve skeletal muscle, heart, and central nervous system (4,5). The dystrophin gene, which is 2.4 Mb and has 79 exons, is the largest known gene in humans. Dystrophin, the product of this gene, is a large cytoskeletal protein of 427 kDa situated along the inner side of the sarcolemma. It comprises four main domains (6 -8). The N-terminal domain constitutes the first actin-binding domain with two calponin homology domains (CH). The major part of the protein is a long central filamentous domain comprising 24 repeating units of about 110 residues each, similar to spectrin repeats. This central domain has several known partners, such as neuronal nitric oxide synthase, filamentous actin, and membrane lipids (9). A Cys-rich domain near the C terminus interacts with the intrinsic sarcolemmal protein ␤-dystroglycan, and the C-terminal domain interacts with dystrobrevins and syntrophins.
Various types of dystrophin gene mutations have been reported (10 -13). Large deletions and duplications are the most common, accounting for ϳ75% of the total. Most of the remaining 25% are small insertions, deletions, exonic point mutations, or splice site mutations. The mutations most often seen in children with BMD are in-frame deletions that allow the expression of an internally truncated dystrophin. In contrast, those causing DMD tend to be out of frame, and very little or no dystrophin is expressed (3). This reading frame rule holds true for 93-96% of DMD or BMD cases with the exceptions being * This work was supported by the Association Française contre les Myopathies and Iowa Wellstone Muscular Dystrophy Cooperative Research Center Grant NS053672. □ S This article contains supplemental Table S1, Figs. S1-S3, and additional references. 1  related to mutations affecting RNA splicing or alternative translation initiation (14 -16). A genotype to phenotype correlation in BMD is likely to be related to the residual functions of the mutated dystrophin. In fact, truncated "mini-dystrophin" constructs are proposed as candidate forms of gene or cell therapy for DMD. To take up this challenge, protein structural studies of mutated or truncated forms of dystrophin will help to predict the functionality of such constructs (17)(18)(19)(20).
In the present paper, we report a new point mutation in brothers presenting with a mild BMD phenotype. The mutation is located in exon 11, which encodes part of the first repeat of the dystrophin rod domain. We describe modified physical properties of this protein domain caused by a leucine to proline amino acid substitution and suggest a correlation between protein misfolding and dysfunction that leads to disease.

EXPERIMENTAL PROCEDURES
The study was exempted from human investigational research study by the Institutional Review Board. Genetic analysis of the DMD gene was performed by Athena Diagnostics (Worcester, MA) as a routine diagnostic test.

Muscle Histopathology
Immunofluorescence was carried out on unfixed cryosections of the patient's muscle biopsy as part of routine diagnostic testing at the University of Iowa. All of the anti-dystrophin protocols are 1-h incubation with primary antibodies and 30-min incubation with secondary antibodies. Detailed information about the primary antibodies is listed in supplemental Table 1. This panel of anti-dystrophin antibodies covers epitopes from the N terminus through the rod domain to the C terminus. Antibodies to the dystroglycans, sarcoglycans, utrophin, neuronal nitric-oxide synthase, spectrin, caveolin-3, and dysferlin were also used to evaluate the biopsy as described previously (21). The secondary antibodies are goat anti-mouse or goat anti-rabbit tagged with Alexa Fluor 488 or Alexa Fluor 594 (Molecular Probes, Invitrogen).

In Vitro Protein Biophysical Study
Materials-The pGEX-4T1 plasmid vector, GSTrap TM HP column were purchased from GE Healthcare. The E. coli BL21(DE3) bacteria were supplied by Invitrogen and restriction enzymes by New England Biolabs.
Cloning-The plasmid pTG11025, kindly provided by Dr. S. Braun (Transgene, Strasbourg, France), is an E. coli plasmid that carries the cDNA encoding the full-length Dp427m isoform of human dystrophin (NCBI Nucleotide Data Base NM-004006). Recombinant dystrophin rod domain wild-type tandem repeat R1-2 was cloned downstream of the sequence of glutathione S-transferase (GST) into pGEX-4T1 vector, using the BamHI restriction site of the preScission thrombin recognition sequence and XhoI restriction site. The L427P mutation was introduced by site-directed mutagenesis carried out by PCR using complementary oligonucleotide primers containing a CCG codon according to E. coli codon usage in place of the patient mutated CCA codon. In addition, this codon allows a new EcoRI restriction site for mutant selection. The cloned fragments were verified by sequencing using the dideoxy chain termination method.
Boundaries of Wild-type and Mutated Proteins-Boundaries of recombinant R1-2 were extended by assuming the original alignment of spectrin-type triple-helical repeats in dystrophin as proposed previously (22). The two recombinant proteins were extended by six additional residues at the C terminus to improve solubility and stability as reported previously (23,24). The sequence range is E-338 to R-563, and two GS residues were residual at the N terminus after thrombin cleavage.
Preparation of Wild-type and Mutated Proteins-Plasmid constructs were used to transform into the E. coli BL21(DE3) strain, and cultures were then grown to an A 600 nm of 0.5 at 37°C in LB broth supplemented with 100 g/ml ampicillin. Expression of GST fusion proteins was induced by addition of 0.5 mM isopropyl-␤-D-thiogalactopyranoside. After 4 h, bacteria were harvested by centrifugation at 4,000 ϫ g for 15 min at 4°C. As described previously, pellets from 500-ml cultures were resuspended in 20 ml of ice-cold lysis buffer (20 mM Tris-HCl, pH 7.5, 150 mM NaCl) and incubated with 0.5 mg⅐ml Ϫ1 lysozyme and 40 g⅐ml Ϫ1 DNase on ice for 30 min before being broken up by sonication. After centrifugation at 10,000 ϫ g for 20 min at 4°C, filter-clarified extracts were loaded onto a preequilibrated 5-ml GSTrap TM FF column at a flow rate of 0.25 ml⅐min Ϫ1 . The column was washed with 10 column volumes of lysis buffer at a flow rate of 1 ml⅐min Ϫ1 , and 50 NIH units of thrombin were then added. On-column cleavage took place for 48 h at 4°C, and fractions of 5 ml were collected in the washing buffer. Fractions were concentrated on 10K Centricon, and the final purification steps consisted of ion exchange chromatography on HiTrap Q FF (GE Healthcare). Protein concentration was assessed using absorbance at 280 nm and purity by SDS-PAGE.
Circular Dichroism (CD) Measurements-CD spectra were acquired at 20°C with a path length of 0.2 cm on a JASCO J-815 spectropolarimeter, equipped with a Peltier device for temperature control. From the molar ellipticity at 222 nm [ 222 ], the mean residue ellipticity [] MRW was calculated using the mean residue molar concentration of the proteins, and the percentage of ␣-helix was obtained using a 100% ␣-helix value of Ϫ36,000 degrees⅐cm 2 ⅐dmol Ϫ1 at 222 nm as described previously (25). The thermal unfolding of each construct was monitored by CD at 222 nm at concentrations of 2.5 M in TNE buffer (20 mM Tris-HCl, pH 7.5, 150 mM NaCl, 0,1 mM EDTA) with temperature increasing 1 degree/min from 15 to 85°C. As refolding after heating was observed, denaturation was demonstrated to be reversible (see "Results"), CD signals were fitted to a two-state transition with the following equation where y is the CD signal at 222 nm at temperature T, ␣ N and ␣ D are the intercepts of the native and denatured states, respectively, and ␤ N and ␤ D are the slopes of the native and denatured states, respectively; ⌬Gun is the free energy of unfolding and R is the gas constant in cal⅐mol Ϫ1 ⅐K Ϫ1 . This equation can be rewritten to facilitate nonlinear regression analysis by Sigma-Plot 10.0 in Windows as follows (25,26), where T m is the melting temperature and ⌬T is the width of the unfolding transition.
Steady-state Fluorescence Measurements-Tryptophan fluorescence spectra of native R1-2 WT and R1-2 L427P proteins at 1 M in TNE buffer were recorded at 295-nm excitation wavelength (bandwidth, 1 nm) on a Fluorolog 3 spectrofluorometer (Horiba Jobin-Yvon, France), using 10 ϫ 4-mm quartz cuvette at 20°C. After appropriate correction, we obtained the emission spectra of nondenatured or urea-denatured proteins. For urea denaturation equilibrium experiments, the proteins were previously incubated for 2 h in TNE with 0.25-8 M urea. The concentration of urea was determined by refractive index measurements. After fluorescence spectra acquisition, we obtained the maximum wavelength ( max ) and intensity of fluorescence from the emission spectra. The max , y, was then plotted as a function of urea concentration, U, and fitted to a onetransition, two-state reversible process as described previously (25) and given in the following equation, where ␣ is the intercept of the native (N) and denatured (D) states max , ␤ the slope of the native (N) and denatured (D) states, U 50% is the urea concentration at the mid-point of denaturation, m is the slope of the transition between the native and denatured states, R is the gas constant, 1.987 cal⅐mol Ϫ1 ⅐K Ϫ1 , and T the temperature in K. The curves were fitted by nonlinear regression with Sigma plot 10.0 in Windows. The term mU 50% in Equation 3 is equal to the free energy of denaturation in the absence of urea, ⌬Gun.
For 8-anilino-1-naphthalenesulfonic acid (ANS) fluorescence measurements, proteins were diluted at 1 M in TNE containing 20 M ANS. ANS was excited at 350 nm and emission scanned from 460 to 620 nm. All background spectra were subtracted.
Refolding and Unfolding Kinetics by Stopped-flow Fluorescence Spectroscopy-Stopped-flow data were recorded on a BioLogic SFM-3, MOS-250 instrument (Grenoble, France). The dead time of the stopped-flow device was 2.2 ms. For refolding measurements, stock solutions of 8 M urea-treated proteins were prepared for 2 h in TNE buffer. The refolding reactions were initiated at 25°C by 10-fold dilution of the ureatreated proteins in TNE buffer. The final protein concentrations were 0.5 M. The time-dependent fluorescence of tryptophan changes was monitored at excitation and emission wavelengths of 295 and 345 nm, respectively. Curves were derived from the averages of at least 10 individual kinetic data image from an H&E stained cryosection photographed using the 20ϫ objective. There are several atrophic fibers and a few mildly hypertrophic fibers. Scale bar, 50 m. Immunofluorescence images compare a pediatric, nondystrophic control biopsy with the BMD patient's biopsy. The antibodies are described in supplemental Table 1. Each image from the BMD patient illustrates the same region of the biopsy; the lowercase a identifies the same muscle fiber in each serial cryosection. The dystrophin regions recognized by each anti-dystrophin antibody are printed beside each side-by-side pair of control BMD images. All immunofluorescence images were obtained on a Zeiss 710 laser confocal microscope using a 20ϫ objective. Scale bar, 100 m. points after subtraction of the background urea buffer signal. Using Biokine software (BioLogic), the experimental results were fitted with exponential functions to obtain the rate constants.
Molecular Models and Dynamics Simulations-A homology model of the R1-2 tandem repeat was obtained previously via the I-Tasser server using a spectrin repeat structural pattern (27). The mutant form L427P of the protein was modeled by the leucine residue substitution with the Yasara program (28). Each protein was described by parameters of the Yamber3 force field, placed in an explicit solvent box using the TIP3P water model, and neutralized (0.9% NaCl). Both systems were simulated under periodic boundary and in the NVT ensemble (298 K). Structures were relaxed during a 1-ns molecular dynamics simulation, providing a stable root mean square deviation over the last 100 ps. Production trajectories of 12 ns were collected at 2-ps intervals for each of the two structural models (R1-2 WT and L427P mutant). Secondary structure measurements were performed using the DSSP program (29). Mean molecular surfaces and molecular hydrophobicity potential were provided by the Platinum webserver.

Clinical and Genetic Findings
The index case is a 10-year-old male who was evaluated initially by a pediatric gastroenterologist because of severe recurrent abdominal pain that started at 8 years of age. The investigations did not reveal any specific gastrointestinal disease, but transaminases were found to be elevated (alanine aminotransferase, 119 IU, normal 30 -65 IU; and aspartate aminotransferase, 81 IU, normal 15-37 IU) along with elevated creatine kinase (CK; 1,632 IU, normal 35-250 IU). The elevated CK triggered evaluation by a child neurologist who noted moderate enlargement of the gastrocnemius muscles and heel cord tightness. The patient had trouble stooping all the way to the ground as well as recovering from a lunge position without using a Gower's maneuver, particularly after multiple trials. Motor testing suggested a mild weakness in hip flexors with Medical Research Council scale of 4. Although the medical history did not suggest any perinatal or developmental problems, a school performance evaluation suggested a form of learning disability with dyslexia. A three-generation pedigree was unremarkable for neuromuscular disease. However, shortly following the patient's diagnosis, his 8-year-old brother presented with recurrent abdominal pain and greatly elevated CK (3,440 IU/liter).
A negative DMD deletion and duplication analysis was followed by complete sequencing of the gene, which identified a novel nucleotide transition TϾC at position 1280. The Leiden Muscular Dystrophy Pages mutation data base does not include this missense mutation. This mutation in codon 427 results in the amino acid substitution leucine to proline (L427P). The same genetic variant of the dystrophin gene was subsequently identified in the index patient's 8-year-old-brother as well as their mother.

Muscle Histology
Routine hematoxylin and eosin staining of a biopsy from the left vastus lateralis muscle showed only mild variation in fiber size due to scattered moderately atrophic and a few hypertrophic muscle fibers (Fig. 1). However, no necrotic and only rare regenerating muscle fibers were seen. Internal nuclei were mildly increased, and rare muscle fibers had clefts or splitting. No endomysial fibrosis and no lymphocytic inflammation were seen in the specimen. Immunofluorescence staining showed greatly reduced to absent immunoreactivity near the N terminus with antibodies recognizing epitopes encoded by exons 7 and 8 (see laser confocal fluorescence microscopy panels in Fig.  1). Dystrophin immunostaining abnormalities were absent or only mild using antibodies to other epitopes. Other proteins of the dystrophin-glycoprotein complex (dystroglycans and sarcoglycans), caveolin-3, and dysferlin stained normally; utrophin was not detected at the sarcolemma (data not shown).

Structural Analysis of Mutated Repeat 1-2 of Dystrophin
Protein Characterization- Fig. 2A shows a schematic representation of human dystrophin and an enlargement of the region from the CH2 of the actin-binding domain 1 to the start of repeat 3 together with the respective encoding exons. The exact sequence of this region together with the exon boundaries are given in supplemental Fig. S1. An alignment of this region for other species shows that the Leu residue in position 427 is highly conserved in dystrophin from mammals and is also present in chicken, zebra fish, and Drosophila. In the same position, a hydrophobic residue is also found in Caenorhabditis elegans (Fig. 2B). The L427P mutation found in our patients is located in helix C of the repeat R1 and is located at a d position known to stabilize hydrophobic helix-helix interacting clusters in the spectrin like triple coiled-coils (22,27,30). The mutation introduces a Pro residue that is known to be a helix breaker (31)(32)(33).
The wild-type and mutated proteins were obtained in comparable purity as assessed by Coomassie Blue-stained SDSpolyacrylamide gels (Fig. 3A). N-terminal sequencing confirmed the expected cleavage by thrombin for both proteins.  MAY 25, 2012 • VOLUME 287 • NUMBER 22

JOURNAL OF BIOLOGICAL CHEMISTRY 18157
CD spectra (Fig. 3B) display the typical features of proteins with a predominant ␣-helix folding as reported previously for dystrophin repeats, i.e. two minima troughs at 222 and 208 nm. However, the spectra from the wild-type and mutated proteins are different. First, the values at 222 nm are half as deep in the mutated compared with the wild type, indicating 73 and 45% of helicity for the wild-type and mutated proteins, respectively. This indicates that the mutated protein is in a different native state compared with the wild-type protein. Second, the ratio of the [222]/[208] values is about 1 for the wild-type protein whereas it is about 0.85 for the mutated protein. Taken together, this indicates that the coiled-coil would not be fully formed in the mutated protein compared with the wild type.
As shown in Fig. 3C, only a slight increase of tryptophan fluorescence intensity was observed for the mutated protein, with no change at the 338-nm maximal emission wavelength. The extent of exposure of hydrophobic surfaces in the R1-2 WT and R1-2 L427P was measured by their ability to bind to the fluorescent dye ANS. Fig. 3D shows that the fluorescence parameters of ANS, intensity and maximal emission wavelength, are in agreement with an extended hydrophobic surface area for the mutated protein compared with the wild type. Onedimensional proton NMR spectra corroborated these observations (supplemental Fig. S2). Several resonances at high field and the dispersion of the amide protons are indicative of a folded wild-type protein. In contrast, the spectrum of the mutated protein rather corresponds to a different folding state of the protein as indicated by the broad NH resonances and the absence of high field signals.
Thermal Denaturation-Thermal denaturation was followed by CD at 222-nm wavelength. The two proteins were denatured in a one-step mechanism (Fig. 4A). However, it should be noted that because the mutated protein has a less helical content than the wild type protein, the starting mean residue ellipticity values of CD at 222 nm are far from identical, in agreement with the spectra first obtained. Spectra acquired at the end of the heating procedure showed that the proteins are all denaturated (Fig. 4, B and C). As large refolding of about 75-80% of ellipticity after heating was observed for the two proteins (Fig. 4, B and C) the equations describing the thermodynamic process of denaturation can be used. Therefore, the temperatures at mid-denaturation were determined to be 53.4 and 50°C for the wild-type and the mutated protein, respectively (Fig. 4A). This indicates a lower thermal stability associated with the mutation.
Urea Unfolding and Refolding-The equilibrium denaturation was studied using tryptophan fluorescence. When denatured in 8 M urea, the intensity of tryptophan fluorescence increased, and the maximal emission wavelength was 348 nm for both proteins (Fig. 5A). The tryptophan fluorescence data follow- ing urea denaturation indicate that R1-2 wild-type and R1-2 L427P proteins are unfolded in a single-step process with identical urea concentrations of 4.10 M at mid-denaturation (Fig. 5B).
Refolding kinetics of the two proteins were then monitored by stopped-flow tryptophan fluorescence. These proteins were first unfolded in 8.0 M urea and then rapidly diluted to a final urea concentration of 0.8 M. The refolding data fitted well to a two-exponential process for both proteins. The fastest phase accounted for 92 and 68% of the amplitude with calculated rate constants of 23 s Ϫ1 and 11 s Ϫ1 for R1-2 WT and R1-2 L427P, respectively. The slowest phase accounted for 8 and 32% of the amplitude with rate constants of 0.02 s Ϫ1 and 0.06 s Ϫ1 for R1-2 WT and R1-2 L427P, respectively (Fig. 5C). This indicates a slower refolding for the mutant compared with the wild type.
Molecular Modeling-The models for the R1-2 WT and L427P were built using I-Tasser and Yasara, respectively. The two proteins exhibit the triple-helical coiled-coil three-dimensional structure typical of the spectrin-like repeats (Fig. 6). The two proteins contained six tryptophan residues, two of them being in repeat 1 and four of them in repeat 2. The Trp-431 situated in the helix C of repeat 1 is 4 residues from the mutation site (Fig. 6, zoom). Therefore, its environment can be modified by the mutation, and this can explain the increase in intrinsic tryptophan fluorescence observed in the mutated protein compared with the wild-type protein. In fact, the nonnative state of the mutated protein due to a modification of the folding around the mutation site likely exposes the Trp-431 to a more hydrophilic environment leading to the observed increase of its fluorescence. The environment of the other tryptophan residues is not modified, and the averaged fluorescence of the six tryptophan residues is only slightly modified.
Molecular Simulations-As expected from experimental structures and theoretical observations (31)(32)(33), the proline insertion caused by the L427P mutation induces R1 helix C destabilization, related to a short helix bend and a break in the H-bonding chain defining the helical structure. Atomic fluctuations shown by the root mean square fluctuation of the 12 ns of molecular dynamics trajectories are clearly increased at the mutation spot (Fig. 7A). The a and d hydrophobic positions of the heptad pattern specific to coiled-coil superstructures remain effective around the mutation. However, the total molecular surface increases by 4%, and the hydrophobic surface contribution moves from 31 to 33%, indicating that the coiledcoil arrangement in the R1 L427P is less stable than the wildtype R1. The main impact of the mutation on the R1-2 dynamics profile concerns the stiffening of the linker region supported by the important freezing in the mobility of the A/B loop. This modification in the flexibility of the coiled-coil structure also strongly affects the conformational space available to the N-terminal region of the R1-2 protein and the relative position of R1 repeat toward the R2 repeat (Fig. 7B). Principal component analysis showed that the first and the third most contributing components to be global motion are the most strongly modified by the mutation (supplemental Fig. S3).

DISCUSSION
Clinical variability in dystrophinopathy-related BMD ranges between a severe, Duchenne-like phenotype and a late onset mild disease. In-frame deletions in the mid-portion of dystrophin gene are responsible for most of the milder cases, but the size of the deletion does not appear to be a key factor in determining the severity of disease (4,34,35). Serum CK levels also show little correlation with the severity of clinical symptoms. Severe, even lethal cardiomyopathy may develop in patients with relatively minimal skeletal muscle weakness, although the genotype correlation has not been determined for these cases (36). Central nervous system involvement and gastrointestinal manifestations even in the early disease stages are more typical  27. The tandem repeat is constituted by the two repeats R1 and R2, each of them comprising three helices, HA, HB, and HC for R1 and HAЈ, HBЈ, and HCЈ for R2. The N and C termini are denoted N and C, respectively The tryptophan residues are denoted as red W, and the L427 site of the mutation is shown in Repeat 1. The region of the mutation is enlarged on the right, showing that the Trp-431 is very near the mutation pointing into the interior of the molecule in the wild-type model. Leucine 427 (dotted gray sphere) points to the interior of the molecule. Therefore, the mutation L427P leads to the proline residue (blue sticks) pointing also to the interior with putative helix breaking; this likely is able to destabilize the folding of these three helices.
of DMD than BMD (37,38). Our index case presented with subtle proximal skeletal muscle weakness and no cardiac involvement. It was not certain whether his recurrent abdominal pain was related to gastrointestinal hypomotility or other dysfunctions. The serum CK elevation and histopathology of the muscle biopsy were milder than is typical of BMD, but both dystrophin immunostaining and DMD gene sequencing confirmed this diagnosis.
A phenotype-genotype correlation is difficult to establish in BMD. In-frame deletions constitute the majority of the mutations, and it appears that mutations of the N terminus of dystrophin affecting the actin-binding domain 1 (39,40) or the Cys-rich domain (41) are more severe than the others. However, deletions around exons 45-55 are in a large number of cases accompanied by mild skeletal muscle disease but a possibly severe cardiomyopathy. Large deletions in the central rod domain (between exons 45 and 59) have been seen in late onset BMD (2). By combining genetic and protein information, it appears that there is a strong association between specific structural modifications of the mutated dystrophin and the age of onset of the cardiomyopathy (42).
Point mutations of the dystrophin gene are less common than large deletions and duplications. Mutated dystrophin is expressed in skeletal muscle of cases with missense point mutations, but to date there are very few publications that evaluate potential genotype-phenotype correlations. Previous studies of point mutations within the actin-binding domain of patients with either DMD or BMD phenotypes have concluded that the dystrophin produced is mostly misfolded, rather unstable and prone to aggregation (39,40). We recently showed that a double point mutation in the central rod domain repeat R23 of a DMD patient leads to a substantially lower thermal and chemical stability as well as slower refolding of the mutated protein compared with the wild-type protein (25).
In this work, the novel c.1280TϾC DMD mutation introduces a proline residue in the helix C of repeat 1. Proline is a residue known as a helix breaker, even though it does not always preclude the helical structure (31)(32)(33) and is found in some helices of wild-type dystrophin repeats (22). In the present case, the data demonstrate that the helical content of the mutated R1-2 fragment is dramatically decreased compared with the wild-type protein, showing that the mutated protein does not FIGURE 7. C␣ root mean square fluctuation over 12-ns molecular dynamics trajectories. A, ␣-helices delimitation (27) for the R1 (A, B, and C) and R2 (AЈ, BЈ, and CЈ) repeats and mutation location (*) are reported along the x axis in correspondence with the residue numbers in the dystrophin sequence. C␣ root mean square fluctuation (in Å) is plotted for each residue of R1-2 WT (black) and R1-2 L427P (green). B, flexibility variations are mapped on the mutant structure. The mutation spot is identified with a sphere. Higher rigidity (blue) or higher flexibility (red) spots of L427P mutant by comparison with the WT structure are shown on the molecular models, respectively. retain a folded state compared with the wild-type protein. This is consistent with the hydrophobic surface changes detected by ANS fluorescence. The mutation also affects thermal unfolding and urea refolding kinetics. The molecular modeling and dynamics presented here are in agreement with these experimental observations and clearly show that the single mutation L427P has striking local and global effects on the dynamic and mechanical behavior of the filamentous structure of the tandem repeat R1-2 of dystrophin. Taken together, these results demonstrate that the proline residue induces a break in the helix which in turn destabilizes the entire helix toward the C terminus of the molecule. The dystrophin protein as expressed in the patient's skeletal muscle is found to be subsarcolemmal just as in control muscle. However, the protein is poorly visualized by immunofluorescence using antibodies that recognize peptides encoded by exons 7 and 8. Given that the area bearing the L427P mutation is encoded by exon 11, the lack of staining of the full-length dystrophin in the patient's biopsy is likely due to conformational changes in hinge 1 consecutive to the rigidity changes we observed in the N-terminal part of R1-2.
It thus can be concluded that, even though full-length dystrophin is detected in the patient, the missense mutation modifies peptide stability and leads to misfolding of the N-terminal and proximal rod domains. This in turn has a negative effect on the mechanical or binding functions of dystrophin leading to BMD.