Hypertrophic cardiomyopathy mutations in the pliant and light chain-binding regions of the lever arm of human β-cardiac myosin have divergent effects on myosin function

  1. Makenna M Morck
  2. Debanjan Bhowmik
  3. Divya Pathak
  4. Aminah Dawood
  5. James Spudich
  6. Kathleen M Ruppel  Is a corresponding author
  1. Stanford Cardiovascular Institute, Stanford University School of Medicine, United States
  2. Department of Biochemistry, Stanford University School of Medicine, United States
9 figures, 2 tables and 3 additional files

Figures

Location of hypertrophic cardiomyopathy-causing mutations along the lever arm.

(A) Homology model of the prestroke β-cardiac myosin subfragment 1 (S1; Homburger et al., 2016) highlighting the five mutations examined in the present study: D778V, L781P, S782N, A797T, and F834L and their locations along the lever arm. The essential light chain (ELC) and regulatory light chain (RLC) are shown in maroon and green, respectively, and the hook joint is highlighted in the red circle. (B) Alignment of β-cardiac myosin lever arms for several model organisms, demonstrating the high degree of conservation of the lever arm across species. Residues are colored by hydrophobicity, where red is most hydrophobic, and blue is most hydrophilic. Mutated residues studied here are indicated with colored stars.

Light chain loading of wild type (WT) vs mutant myosins.

(A) Representative gel used to quantitate light chain loading (L781P 2-hep is shown). For each sample, myosin was titrated in a denaturing SDS-PAGE gel, stained with Coomassie, and scanned for 700 nm fluorescence. See Materials and Methods. (B) Light chain loading across all WT 2-hep samples measured. Mean essential light chain/heavy chain (ELC/HC) = 0.19 ± 0.02, regulatory light chain/heavy chain (RLC/HC) = 0.11 ± 0.02 (mean ± SD). Expected ELC/HC and RLC/HC ratios, if Coomassie staining were unbiased, are 0.168 and 0.146. (C) Light chain loading of mutant 2-hep myosins normalized to WT controls. For each mutant, light chain loading was not statistically different from WT.

Figure 2—source data 1

Raw uncropped gel image.

Uncropped representative gel used to quantitate light chain loading (L781P 2-hep is shown) ± annotation.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig2-data1-v1.pptx
Figure 2—source data 2

Ratio of ELC and RLC to myosin 2-hep heavy chain.

(2B) Ratio of pixel density of ELC/RLC over WT 2-hep myosin heavy chain, based on fluorescent imaging of Coomassie stained gels. (2 C) Ratio of LC/HC for mutants over same-day WT LC/HC controls.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig2-data2-v1.xlsx
Figure 3 with 1 supplement
Actin-activated ATPase activity and in vitro motility velocity of 2-hep human β-cardiac myosin lever arm mutants vs wild-type (WT).

(A–E) Representative actin-activated ATPase curves of 2-hep constructs for each lever arm mutant, normalized to their WT controls. Each plot shows a single biological replicate (one of two, see Supplementary file 1), where error bars represent the SD of three technical replicates. Where error bars are not shown, error is smaller than the size of the data point. Mutant data are plotted against their prep-matched WT 2-hep control (gray) and normalized to the WT kcat value. Curves are fit to Michaelis-Menten kinetics, and the shaded areas represent 95% CI of the fits. Only D778V produced a significant difference: the kcat was increased by 16 ± 9% (average of two independent biological replicates compared to WT controls). (F) In vitro motility velocities of lever arm mutants. Each data point represents the average velocity of the mutant 2-hep on a single slide normalized to the WT 2-hep velocity from the same slide. Bars represent the average of the data points. D778V increased the velocity by 46 ± 8% over WT, L781P reduced the velocity by 30 ± 8% compared to WT, A797T reduced the velocity by 5 ± 6%, and S782N and F834L had no significant effects on motility velocity (mean ± SD of the data points shown). **** indicates p≤0.0001, * indicates p≤0.05.

Figure 3—source data 1

Rate vs (actin) for WT and mutant 2-hep myosins.

Turnover rate per second vs (actin) of the actin-activated ATPase activity of each mutant 2-hep myosin along with its same day WT control. Raw rates are depicted in the left-hand boxes and rates normalized to WT 2-hep control are shown in the right-hand boxes. Data shown is for all technical replicates of the depicted representative biological replicate for each mutant/WT pair.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig3-data1-v1.xlsx
Figure 3—figure supplement 1
Full actin-activated ATPase results.

(A–O) Each plot shows a single biological replicate of the ATPase results reported here (raw data, not normalized). Data points represent the average of three technical replicates, and error bars show SD (where error bars are not shown, error is smaller than the size of the data point). Curves are fit to Michaelis-Menten kinetics, and shaded areas represent the SE of the fit. Every curve is a separate biological replicate (i.e. each WT 2-hep, WT 25-hep, mutant 2-hep, and mutant 25-hep shown are from a separate batch of cells and protein preparation). Data are shown in the same order presented in Supplementary file 1, where fitting results for each curve are detailed.

Figure 3—figure supplement 1—source data 1

Turnover rate per second vs (actin) of the actin-activated ATPase activity of each mutant 2-hep myosin and each mutant 25-hep myosin along with its same day WT 2-hep control.

Raw rates are depicted in the left-hand boxes, rates of mutant 2-hep normalized to WT 2-hep control are shown in the middle boxes, and rates of mutant 25-hep normalized to mutant 2-hep are shown in the right-hand side boxes. For the five WT 2-hep/25-hep biological replicates, the raw rates are shown. Data shown is for all technical replicates of all biological replicates.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig3-figsupp1-data1-v1.xlsx
Figure 4 with 2 supplements
Harmonic force spectroscopy measurements of pliant region mutations in short subfragment 1 (sS1) human β-cardiac myosin.

(A) Detachment rate of sS1 myosin as a function of load force, average of all molecules (wild type [WT]: n=10 molecules, D778V: n=12 molecules, L781P: n=22 molecules, S782N: n=11 molecules). WT curve is shown in gray, D778V curve is shown in red, L781P curve is shown in purple, and S782N curve is shown in cyan. Dotted lines show error propagated from SEM of fitted parameters. (B) k0, the detachment rate at zero load force, for WT and each mutant sS1. D778V increased the detachment rate by 115 ± 9% (mean ± SEM), while L781P reduced the detachment rate by 31 ± 10%, and S782N had no significant effect on detachment rate. (C) δ, the measure of force sensitivity, for WT and each mutant sS1. All the pliant region mutations reduced the force sensitivity: D778V by 35 ± 11%, L781P by 22 ± 8%, and S782N by 21 ± 10% (mean ± SEM). (D) Step size for WT and each mutation. Only L781P affected step size, reducing it by 42 ± 22% (mean ± SEM). Each data point represents the average value for one molecule, and bars show the average of the data points. ** indicates p≤0.01, *** indicates p≤0.001, **** indicates p≤0.0001.

Figure 4—source data 1

Detachment rate, δ, and step size for all molecules of WT and pliant region mutant myosins.

Detachment rate k0, δ, and step size (nm) for all molecules of WT and mutant sS1 myosin, as measured by harmonic force spectroscopy. Each row represents an individual molecule. Molecules are show in the same order for each group of measurements.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig4-data1-v1.xlsx
Figure 4—figure supplement 1
Representative raw data and analysis for harmonic force spectroscopyHFS measurements.

(A–D) Detection of single actin-myosin interaction events (representative) from the change in the amplitude and the phase (region marked with cyan color for bead 1 [x1] and yellow color for bead 2 [x2]) of the oscillating trapped beads for wild type (WT) (A), D778V (B), L781P (C), and S782N (D). Stage oscillation is shown in the top panel of (A–H). Detachment rates (one point corresponds to one event) obtained from HFS measurements of one representative molecule of WT (E), D778V (F), L781P (G), and S782N (H) at different assistive (<0), zero and resistive (>0) load forces. (I–L) Representative fitting using Arrhenius equation (corrected for harmonic force, equation 1) of the detachment rates obtained at different load forces for one molecule each of WT (I), D778V (J), L781P (K), and S782N (L).(M–P) Detachment rate curves obtained for every molecule of WT (M), D778V (N), L781P (O), and S782N (P), where the thick line shows the average of all of the molecules (same as Figure 3A). (Q – T) Fitting of histogram of displacements (obtained from each event) of one representative molecule of WT (Q), D778V (R), L781P (S), and S782N (T).

Figure 4—figure supplement 1—source data 1

(I–L) Data fof representative fitting using Arrhenius equation (corrected for harmonic force, equation 1) of the detachment rates (Kdet) obtained at different load forces for one molecule each of WT (I), D778V (J), L781P (K), and S782N (L).

(Q–T) Data of displacement measurements binned in the histogram of displacements (obtained from each event) of one representative molecule of WT (Q), D778V (R), L781P (S), and S782N (T).

https://cdn.elifesciences.org/articles/76805/elife-76805-fig4-figsupp1-data1-v1.xlsx
Figure 4—figure supplement 2
Harmonic force spectroscopy measurements across unique protein preparations.

(A–C) Plots show the same data from Figure 3B–D with data points (representing unique molecules) from independent protein preparations shown in different colors and with different shapes. Wild type (WT), L781P, and S782N molecules come from two unique protein preparations, and D778V molecules come from three unique protein preparations.

Figure 4—figure supplement 2—source data 1

Detachment rate, delta, and step size for all molecules of WT and pliant region mutant myosins across separate myosin preps.

Data matches data from Figure 4—source data 1, with the addition of the column entitled ‘prep number; which denotes the myosin prep from which each molecule originated.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig4-figsupp2-data1-v1.xlsx
Effects of pliant region mutations on ensemble duty ratio, average force, and average power.

Calculation of (A) duty ratio, (B) average force, and (C) average power output at resistive forces based on measured detachment rate k0, actin-activated ATPase rate kcat, force-dependent detachment rate kdet(F), and step size (see equations 2-6). Wild type (WT) curves are shown in gray, D778V curves are shown in red, L781P curves are shown in purple, and S782N curves are shown in cyan, where dotted lines in the same color show error, propagated from the SEM of the individual molecules.

Actin-activated ATPase rates of 2-hep vs 25-hep β-cardiac myosin constructs.

(A–F) Representative actin-activated ATPase curves for each mutant, normalized to the mutant 2-hep control kcat. (A) is reproduced from Vander Roest et al., 2021. Each data point represents the average of three technical replicates of one biological replicate (one of two, see Figure 3—figure supplement 1 and Supplementary file 1), and error bars represent SDs. Where error bars are not shown, error is smaller than the size of the data point. Curves are fitted to Michaelis-Menten kinetics, and shaded areas represent the 95% CI of the fits. Arrows with percentages on each graph show the percent change from the 2-hep to the matched 25-hep (average of both biological replicates, not just the representative curve shown). See Supplementary file 1 for full results.

Figure 6—source data 1

Rate vs (actin) for wild type (WT) and mutant 2- vs 25-hep myosin.

Turnover rate per second vs (actin) of the actin-activated ATPase activity of each mutant 25-hep and 2-hep myosin pair. Raw rates are depicted in the left-hand boxes and 25-hep myosin rates normalized its mutant 2-hep myosin control are shown in the right-hand boxes. Data shown is for all technical replicates of the depicted representative biological replicate for each mutant 2-hep/25-hep pair.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig6-data1-v1.xlsx
Figure 7 with 2 supplements
Single ATP turnover kinetics for 25-hep myosin lever arm mutants compared to wild type (WT).

(A–F) Single mant-ATP turnover curves for WT 25-hep and each mutant. Thin curves show curves fitted to a biexponential decay and normalized to the fitted Y0=1.0 and plateau=0.0 from each replicate, while thick lines show average fitted curves across all replicates. Note that in several cases, the average curve obscures individual replicate curves underneath. The dotted black line represents a simulated single-exponential decay with the slow rate of the average curve, and the dotted gray line represents a simulated single-exponential decay with the fast rate of the average curve. (G) Percent slow phase of multiple replicates of each mutant. Pairwise comparisons show that only A797T and F834L significantly reduce the super relaxed state (SRX) state, resulting in a 31 ± 4% and 41 ± 4% (mean ± SEM) reduction in slow phase, respectively. S782N resulted in an 11 ± 4% (mean ± SEM) increase in SRX. * indicates p≤0.05, **** indicates p≤0.0001.

Figure 7—source data 1

Fast and slow rates and percent slow turnover for all 25-hep single turnover experiments.

Each row represents an individual technical replicate, where rows are matched by replicate across the three data groups.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig7-data1-v1.xlsx
Figure 7—figure supplement 1
Representative single ATP turnover results for 25-hep myosins.

(A–F) Representative single ATP turnover data (semi-transparent data points) for each wild type (WT) and mutant 25-hep, fitted to a biexponential decay (solid line) and normalized to the fitted Y0=1.0 and fitted plateau=0.0. The gray dotted line shows a simulated single-exponential decay with the fast rate of the shown fitted curve, and the black dotted line represents a simulated single-exponential decay with the slow rate of the shown fitted curve.

Figure 7—figure supplement 1—source data 1

Representative single turnover data for a single technical replicate of each WT and mutant 25-hep.

Time (seconds) vs normalized fluorescence, where Y0=1.0 of the fitted curve and 0.0 corresponds to the plateau value of the fitted curve. Raw (non-normalized) fluorescence prior to normalization is also shown.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig7-figsupp1-data1-v1.xlsx
Figure 7—figure supplement 2
Single ATP turnover kinetics for wild type (WT) vs mutant 2-hep myosins.

(A–F) Single mant-ATP turnover curves for WT 2-hep and each mutant. Thin curves show curves fitted to a biexponential decay and normalized to the fitted Y0=1.0 and plateau=0.0 from each replicate, while thick lines show average fitted curves across all replicates. Average curves obscure individual replicate curves underneath in some cases. The dotted black line represents a simulated single-exponential decay with the slow rate of the average curve, and the dotted gray line represents a simulated single-exponential decay with the fast rate of the average curve. (G) Percent slow phase of multiple replicates of each WT and mutant 2-hep. For all mutations, the slow phase is similarly low as compared to the WT control.

Figure 7—figure supplement 2—source data 1

Fast and slow rates and percent slow turnover for all 2-hep single turnover experiments.

Each row represents an individual technical replicate, where rows are matched by replicate across the three data groups.

https://cdn.elifesciences.org/articles/76805/elife-76805-fig7-figsupp2-data1-v1.xlsx
Contributions of lever arm position to the folded state structure.

(A–F) Homology modeled structures of human β-cardiac myosin in the interacting heads motif (IHM) structural state from (A) (Alamo et al., 2017), (B) (Nag et al., 2017), and (C) (Robert-Paganin et al., 2018) and cryo-EM solved structures of smooth muscle myosin in the IHM structural state from (D) (Scarff et al., 2020), (E) (Yang et al., 2020), and (F) (Heissler et al., 2021). In each structure, the pliant region is colored yellow, the essential light chain (ELC)-binding region is colored blue, and the regulatory light chain (RLC)-binding region is colored pink. The myosin heads are in gray, ELCs are in dark red and brown, RLCs are in green, and subfragment 2 tail regions are in cyan. Homology modeled structures (A–C) show a markedly different lever arm structure as compared to experimentally determined structures (D–F). (G) A797 (light green) is located in the region where the lever arm bends between the ELC and RLC. This region has been previously implicated in formation of the IHM (Houdusse and Cohen, 1996). (H) F834 (blue) is located very near to the hook or ankle joint at the end of the lever arm. This joint has been previously shown to be important in the formation of the IHM (Brown et al., 2011). (G) and (H) are from a homology modeled structure of human β-cardiac myosin in the prestroke state (Farman et al., 2014; see Figure 1A).

Technical details of the dual-beam optical trap experiment.

(A) The sample chamber (bottom) is shown in an orientation suitable for an inverted microscope. On top, the typical arrangement of the protein complexes on top of a platform bead is depicted. (B) A typical stroking event—as expected during actin-myosin interaction in a harmonic force spectroscopic (HFS) setup—is depicted. In a standard HFS experiment, the stage oscillates, resulting in a variety of assistive or resistive external load forces applied to the stroking myosin.

Tables

Key resources table
Reagent type (species) or resourceDesignationSource or referenceIdentifiersAdditional information
Gene (Homo sapiens)MYH7NCBI GeneGene ID: 4625
Gene (H. sapiens)MYL3NCBI GeneGene ID: 4634N-terminal FLAG-TEV tag added
Gene (H. sapiens)MYL2NCBI GeneGene ID: 4633N-terminal His-TEV tag added
Strain, strain background (Escherichia coli)BJ5183-AD-1Agilent200,157
Strain, strain background (E. coli)Rosetta (DE3) pLysSSigma-Aldrich70,956
Cell line (H. sapiens)HEK 293TATCCCRL-3216
Cell line (H. sapiens)C2C12ATCCCRL-1772
Transfected construct (H. sapiens)pAdEasy-1-MYH7 sS110.1073/pnas.1309493110
Transfected construct (H. sapiens)pAdEasy-1-MYH7 2hep10.1038/nsmb.3408
Transfected construct (H. sapiens)pAdEasy-1-MYH7 25hep10.1038/nsmb.3408
Biological sample (Bos taurus)Bovine cardiac acetone powderPelfreez57,195
Recombinant DNA reagentpShuttle-CMV- MYH7 sS110.1073/pnas.1309493110
Recombinant DNA reagentpShuttle-CMV- MYH7 2hep10.1038/nsmb.3408
Recombinant DNA reagentpShuttle-CMV- MYH7 25hep10.1038/nsmb.3408
Sequence-based reagentD778V SThis paperMutagenesis primerGGAGGAAATGAGGGTCGAGAGGCTGAGCC
Sequence-based reagentD778V ASThis paperMutagenesis primerGGCTCAGCCTCTCGACCCTCATTTCCTCC
Sequence-based reagentL781P SThis paperMutagenesis primerGATGATGCGGCTCGGCCTCTCGTCCCT
Sequence-based reagentL781P ASThis paperMutagenesis primerAAGGACGAGAGGCCGAGCCGCATCATC
Sequence-based reagentS782N SThis paperMutagenesis primerTGAGGGACGAGAGGCTGAACCGCATCATC
Sequence-based reagentS782N ASThis paperMutagenesis primerGATGATGCGGTTCAGCCTCTCGTCCCTCA
Sequence-based reagentA797T SThis paperMutagenesis primerGTACTCCATTCTGGTGAGCACACCTCGGG
Sequence-based reagentA797T ASThis paperMutagenesis primerCCCGAGGTGTGCTCACCAGAATGGAGTAC
Sequence-based reagentF834L SThis paperMutagenesis primerCTGGATGAAGCTCTACTTAAAGATCAAGCCGCTG
Sequence-based reagentF834L ASThis paperMutagenesis primerCAGCGGCTTGATCTTTAAGTAGAGCTTCATCCAG
Software, algorithmFAST10.1016 /j.celrep.2015.04.006Filament tracking and velocity measurement software
Author response table 1
Full numbers of biological and technical replicate.
AssayProteinBiological ReplicatesTechnical Replicates
Light-chain loading
assay
(Each biological replicate is from a single unique protein prep diluted across 5 gel lanes, implying 5 technical replicates per biological replicate)
WT 2-hep945
D778V 2-hep315
L781P 2-hep315
S782N 2-hep315
A797T 2-hep315
F834L 2-hep315
Actin-activated
ATPase assay
(Biological replicates are from separate protein preps, technical replicates are from separate wells on the same plate)
WT 2-hep1545
WT 25-hep515
D778V 2-hep26
D778V 25-hep26
L781P 2-hep26
L781P 25-hep26
S782N 2-hep26
S782N 25-hep26
A797T 2-hep26

Additional files

Supplementary file 1

Full actin-activated ATPase results.

For each mutation, two independent biological experiments were performed with freshly prepared myosins, where each biological replicate was measured with technical triplicates (6 total replicates for each protein). For WT 2- vs 25-hep, five biological replicates were measured with technical triplicates (15 total replicates for each protein). The 2-hep and 25-hep mutant myosins were prepared in tandem with a WT 2-hep control for comparison, given that actin-activated ATPase results show a slight drift from day-to-day. Results in the fourth and fifth columns, respectively, show the fitted values for kcat (s–1) and Kapp (μM) ± SE of the fit for each biological replicate. The kcat ratio of mutant 2-hep/WT 2-hep and mutant 25-hep/mutant 2-hep for each independent biological replicate is shown in the sixth and seventh columns, respectively, where the error is propagated from SE of the fit for each measurement. In the rightmost column (average mutant 25-hep/mutant 2-hep kcat ratio), statistically significant differences for mutant ratios vs WT ratio are shown, where * indicates p≤0.05, ** indicates p≤0.01, and *** indicates p≤0.001.

https://cdn.elifesciences.org/articles/76805/elife-76805-supp1-v1.pptx
Supplementary file 2

Full single ATP turnover results.

Single ATP turnover experiments were performed for both 2-hep and 25-hep myosins. Results show n, the number of independent experiments for each protein, average percent fast and slow phase ± SD, and average fast and slow rates ± SD. Stars denote statistically significant changes as compared to the WT control of the same protein construct, where * indicates p≤0.05 and **** indicates p≤0.0001.

https://cdn.elifesciences.org/articles/76805/elife-76805-supp2-v1.pptx
Transparent reporting form
https://cdn.elifesciences.org/articles/76805/elife-76805-transrepform1-v1.docx

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Makenna M Morck
  2. Debanjan Bhowmik
  3. Divya Pathak
  4. Aminah Dawood
  5. James Spudich
  6. Kathleen M Ruppel
(2022)
Hypertrophic cardiomyopathy mutations in the pliant and light chain-binding regions of the lever arm of human β-cardiac myosin have divergent effects on myosin function
eLife 11:e76805.
https://doi.org/10.7554/eLife.76805