Dilated cardiomyopathy-associated skeletal muscle actin (ACTA1) mutation R256H disrupts actin structure and function and causes cardiomyocyte hypocontractility.

Skeletal muscle actin (ACTA1) mutations are a prevalent cause of skeletal myopathies consistent with ACTA1's high expression in skeletal muscle. Rare de novo mutations in ACTA1 associated with combined cardiac and skeletal myopathies have been reported, but ACTA1 represents only ~20% of the total actin pool in cardiomyocytes, making its role in cardiomyopathy controversial. Here we demonstrate how a mutation in an actin isoform expressed at low levels in cardiomyocytes can cause cardiomyopathy by focusing on a unique ACTA1 variant, R256H. We previously identified this variant in a family with dilated cardiomyopathy, who had reduced systolic function without clinical skeletal myopathy. Using a battery of multiscale biophysical tools, we show that R256H has potent effects on ACTA1 function at the molecular scale and in human cardiomyocytes. Importantly, we demonstrate that R256H acts in a dominant manner, where the incorporation of small amounts of mutant protein into thin filaments is sufficient to disrupt molecular contractility, and that this effect is dependent on the presence of troponin and tropomyosin. To understand the structural basis of this change in regulation, we resolved a structure of R256H filaments using cryoelectron microscopy, and we see alterations in actin's structure that have the potential to disrupt interactions with tropomyosin. Finally, we show that ACTA1R256H/+ human-induced pluripotent stem cell cardiomyocytes demonstrate reduced contractility and sarcomeric organization. Taken together, we demonstrate that R256H has multiple effects on ACTA1 function that are sufficient to cause reduced contractility and establish a likely causative relationship between ACTA1 R256H and clinical cardiomyopathy.

Novel coenzyme Q6 genetic variant increases susceptibility to pneumococcal disease.

Acute lower respiratory tract infection (ALRI) remains a major worldwide cause of childhood mortality, compelling innovation in prevention and treatment. Children in Papua New Guinea (PNG) experience profound morbidity from ALRI caused by Streptococcus pneumoniae. As a result of evolutionary divergence, the human PNG population exhibits profound genetic variation and diversity. To address unmet health needs of children in PNG, we tested whether genetic variants increased ALRI morbidity. Whole-exome sequencing of a pilot child cohort identified homozygosity for a novel single-nucleotide variant (SNV) in coenzyme Q6 (COQ6) in cases with ALRI. COQ6 encodes a mitochondrial enzyme essential for biosynthesis of ubiquinone, an electron acceptor in the electron transport chain. A significant association of SNV homozygosity with ALRI was replicated in an independent ALRI cohort (P = 0.036). Mice homozygous for homologous mouse variant Coq6 exhibited increased mortality after pneumococcal lung infection, confirming causality. Bone marrow chimeric mice further revealed that expression of variant Coq6 in recipient (that is, nonhematopoietic) tissues conferred increased mortality. Variant Coq6 maintained ubiquinone biosynthesis, while accelerating metabolic remodeling after pneumococcal challenge. Identification of this COQ6 variant provides a genetic basis for increased pneumonia susceptibility in PNG and establishes a previously unrecognized role for the enzyme COQ6 in regulating inflammatory-mediated metabolic remodeling.

Danicamtiv reduces myosin's working stroke but enhances contraction by activating the thin filament.

Heart failure is a leading cause of death worldwide, and even with current treatments, the 5-year transplant-free survival rate is only ∼50-70%. As such, there is a need to develop new treatments for patients that improve survival and quality of life. Recently, there have been efforts to develop small molecules for heart failure that directly target components of the sarcomere, including cardiac myosin. One such molecule, danicamtiv, recently entered phase II clinical trials; however, its mechanism of action and direct effects on myosin's mechanics and kinetics are not well understood. Using optical trapping techniques, stopped flow transient kinetics, and in vitro reconstitution assays, we found that danicamtiv reduces the size of cardiac myosin's working stroke, and in contrast to studies in muscle fibers, we found that it does not affect actomyosin detachment kinetics at the level of individual crossbridges. We demonstrate that danicamtiv accelerates actomyosin association kinetics, leading to increased recruitment of myosin crossbridges and subsequent thin filament activation at physiologically-relevant calcium concentrations. Finally, we computationally model how the observed changes in mechanics and kinetics at the level of single crossbridges contribute to increased cardiac contraction and improved diastolic function compared to the related myotrope, omecamtiv mecarbil. Taken together, our results have important implications for the design of new sarcomeric-targeting compounds for heart failure. SIGNIFICANCE STATEMENT: Heart failure is a leading cause of death worldwide, and there is a need to develop new treatments that improve outcomes for patients. Recently, the myosin-binding small molecule danicamtiv entered clinical trials for heart failure; however, its mechanism at the level of single myosin crossbridges is not well understood. We determined the molecular mechanism of danicamtiv and showed how drug-induced molecular changes can mechanistically increase heart contraction. Moreover, we demonstrate fundamental differences between danicamtiv and the related myosin-binding small molecule omecamtiv mecarbil that explain the improved diastolic function seen with danicamtiv. Our results have important implications for the design of new therapeutics for heart failure.

An Unbiased Proteomic Platform for Activity-based Arginylation Profiling.

Protein arginylation is an essential posttranslational modification (PTM) catalyzed by arginyl-tRNA-protein transferase 1 (ATE1) in mammalian systems. Arginylation features a post-translational conjugation of an arginyl to a protein, making it extremely challenging to differentiate from translational arginine residues with the same mass in a protein sequence. Here we present a general activity-based arginylation profiling (ABAP) platform for the unbiased discovery of arginylation substrates and their precise modification sites. This method integrates isotopic arginine labeling into an ATE1 assay utilizing biological lysates (ex vivo) rather than live cells, thus eliminating translational bias derived from the ribosomal activity and enabling bona fide arginylation identification using isotopic features. ABAP has been successfully applied to an array of peptide, protein, cell, patient, and animal tissue samples using 20 µg sample input, with 229 unique arginylation sites revealed from human proteomes. Representative sites were validated and followed up for their biological functions. The developed platform is globally applicable to the aforementioned sample types and therefore paves the way for functional studies of this difficult-to-characterize protein modification.

Structural dynamics of the intrinsically disordered linker region of cardiac troponin T.

The cardiac troponin complex, composed of troponins I, T, and C, plays a central role in regulating the calcium-dependent interactions between myosin and the thin filament. Mutations in troponin can cause cardiomyopathies; however, it is still a major challenge to connect how changes in sequence affect troponin's function. Recent high-resolution structures of the thin filament revealed critical insights into the structure-function relationship of troponin, but there remain large, unresolved segments of troponin, including the troponin-T linker region that is a hotspot for cardiomyopathy mutations. This linker region is predicted to be intrinsically disordered, with behaviors that are not well described by traditional structural approaches; however, this proposal has not been experimentally verified. Here, we used a combination of single-molecule Förster resonance energy transfer (FRET), molecular dynamics simulations, and functional reconstitution assays to investigate the troponin-T linker region. We show that in the context of both isolated troponin and the fully regulated troponin complex, the linker behaves as a dynamic, intrinsically disordered region. This region undergoes polyampholyte expansion in the presence of high salt and distinct conformational changes during the assembly of the troponin complex. We also examine the ΔE160 hypertrophic cardiomyopathy mutation in the linker and demonstrate that it does not affect the conformational dynamics of the linker, rather it allosterically affects interactions with other troponin complex subunits, leading to increased molecular contractility. Taken together, our data clearly demonstrate the importance of disorder within the troponin-T linker and provide new insights into the molecular mechanisms driving the pathogenesis of cardiomyopathies.

Insights into posttranslational regulation of skeletal muscle contractile function by the acetyltransferases, p300 and CBP.

Mice with skeletal muscle-specific and inducible double knockout of the lysine acetyltransferases, p300 (E1A binding protein p300) and CBP (cAMP-response element-binding protein binding protein), referred to as i-mPCKO, demonstrate a dramatic loss of contractile function in skeletal muscle and ultimately die within 7 days. Given that many proteins involved in ATP generation and cross-bridge cycling are acetylated, we investigated whether these processes are dysregulated in skeletal muscle from i-mPCKO mice and, thus, whether they could underlie the rapid loss of muscle contractile function. Just 4-5 days after inducing knockout of p300 and CBP in skeletal muscle from adult i-mPCKO mice, there was ∼90% reduction in ex vivo contractile function in the extensor digitorum longus (EDL) and a ∼65% reduction in in vivo ankle dorsiflexion torque, as compared with wild type (WT; i.e., Cre negative) littermates. Despite this profound loss of contractile force in i-mPCKO mice, there were no genotype-driven differences in fatigability during repeated contractions, nor were there genotype differences in mitochondrial-specific pathway enrichment of the proteome, intermyofibrillar mitochondrial volume, or mitochondrial respiratory function. As it relates to cross-bridge cycling, remarkably, the overt loss of contractile function in i-mPCKO muscle was reversed in permeabilized fibers supplied with exogenous Ca2+ and ATP, with active tension being similar between i-mPCKO and WT mice, regardless of Ca2+ concentration. Actin-myosin motility was also similar in skeletal muscle from i-mPCKO and WT mice. In conclusion, neither mitochondrial abundance/function, nor actomyosin cross-bridge cycling, are the underlying driver of contractile dysfunction in i-mPCKO mice.NEW & NOTEWORTHY The mechanism underlying dramatic loss of muscle contractile function with inducible deletion of both E1A binding protein p300 (p300) and cAMP-response element-binding protein binding protein (CBP) in skeletal muscle remains unknown. Here, we find that impairments in mitochondrial function or cross-bridge cycling are not the underlying mechanism of action. Future work will investigate other aspects of excitation-contraction coupling, such as Ca2+ handling and membrane excitability, as contractile function could be rescued by permeabilizing skeletal muscle, which provides exogenous Ca2+ and bypasses membrane depolarization.

Dilated cardiomyopathy-associated skeletal muscle actin (ACTA1) mutation R256H disrupts actin structure and function and causes cardiomyocyte hypocontractility.

Skeletal muscle actin (ACTA1) mutations are a prevalent cause of skeletal myopathies consistent with ACTA1's high expression in skeletal muscle. Rare de novo mutations in ACTA1 associated with combined cardiac and skeletal myopathies have been reported, but ACTA1 represents only ~20% of the total actin pool in cardiomyocytes, making its role in cardiomyopathy controversial. Here we demonstrate how a mutation in an actin isoform expressed at low levels in cardiomyocytes can cause cardiomyopathy by focusing on a unique ACTA1 mutation, R256H. We previously identified this mutation in multiple family members with dilated cardiomyopathy (DCM), who had reduced systolic function without clinical skeletal myopathy. Using a battery of multiscale biophysical tools, we show that R256H has potent functional effects on ACTA1 function at the molecular scale and in human cardiomyocytes. Importantly, we demonstrate that R256H acts in a dominant manner, where the incorporation of small amounts of mutant protein into thin filaments is sufficient to disrupt molecular contractility, and that this effect is dependent on the presence of troponin and tropomyosin. To understand the structural basis of this change in regulation, we resolved a structure of R256H filaments using Cryo-EM, and we see alterations in actin's structure that have the potential to disrupt interactions with tropomyosin. Finally, we show that ACTA1R256H/+ human induced pluripotent stem cell cardiomyocytes demonstrate reduced contractility and sarcomeric disorganization. Taken together, we demonstrate that R256H has multiple effects on ACTA1 function that are sufficient to cause reduced contractility and establish a likely causative relationship between ACTA1 R256H and clinical cardiomyopathy.

Assessing Cardiac Contractility From Single Molecules to Whole Hearts.

Fundamentally, the heart needs to generate sufficient force and power output to dynamically meet the needs of the body. Cardiomyocytes contain specialized structures referred to as sarcomeres that power and regulate contraction. Disruption of sarcomeric function or regulation impairs contractility and leads to cardiomyopathies and heart failure. Basic, translational, and clinical studies have adapted numerous methods to assess cardiac contraction in a variety of pathophysiological contexts. These tools measure aspects of cardiac contraction at different scales ranging from single molecules to whole organisms. Moreover, these studies have revealed new pathogenic mechanisms of heart disease leading to the development of novel therapies targeting contractility. In this review, the authors explore the breadth of tools available for studying cardiac contractile function across scales, discuss their strengths and limitations, highlight new insights into cardiac physiology and pathophysiology, and describe how these insights can be harnessed for therapeutic candidate development and translational.

Harnessing molecular mechanism for precision medicine in dilated cardiomyopathy caused by a mutation in troponin T.

Familial dilated cardiomyopathy (DCM) is frequently caused by autosomal dominant point mutations in genes involved in diverse cellular processes, including sarcomeric contraction. While patient studies have defined the genetic landscape of DCM, genetics are not currently used in patient care, and patients receive similar treatments regardless of the underlying mutation. It has been suggested that a precision medicine approach based on the molecular mechanism of the underlying mutation could improve outcomes; however, realizing this approach has been challenging due to difficulties linking genotype and phenotype and then leveraging this information to identify therapeutic approaches. Here, we used multiscale experimental and computational approaches to test whether knowledge of molecular mechanism could be harnessed to connect genotype, phenotype, and drug response for a DCM mutation in troponin T, deletion of K210. Previously, we showed that at the molecular scale, the mutation reduces thin filament activation. Here, we used computational modeling of this molecular defect to predict that the mutant will reduce cellular and tissue contractility, and we validated this prediction in human cardiomyocytes and engineered heart tissues. We then used our knowledge of molecular mechanism to computationally model the effects of a small molecule that can activate the thin filament. We demonstrate experimentally that the modeling correctly predicts that the small molecule can partially rescue systolic dysfunction at the expense of diastolic function. Taken together, our results demonstrate how molecular mechanism can be harnessed to connect genotype and phenotype and inspire strategies to optimize mechanism-based therapeutics for DCM.

Insights into post-translational regulation of skeletal muscle contractile function by the acetyltransferases, p300 and CBP.

Mice with skeletal muscle-specific inducible double knockout of the lysine acetyltransferases, p300 (E1A binding protein p300) and CBP (cAMP-response element-binding protein binding protein), referred to as i-mPCKO, demonstrate a dramatic loss of contractile function in skeletal muscle and ultimately die within 7 days. Given that many proteins involved in ATP generation and cross-bridge cycling are acetylated, we investigated whether these processes are dysregulated in skeletal muscle from i-mPCKO mice and thus could underlie the rapid loss of muscle contractile function. Just 4-5 days after inducing knockout of p300 and CBP in skeletal muscle from adult i-mPCKO mice, there was ∼90% reduction in ex vivo contractile function in the extensor digitorum longus (EDL) and a ∼65% reduction in in vivo ankle dorsiflexion torque, as compared to wildtype (WT; i.e. Cre negative) littermates. Despite the profound loss of contractile force in i-mPCKO mice, there were no genotype-driven differences in fatigability during repeated contractions, nor were there genotype differences in mitochondrial specific pathway enrichment of the proteome, intermyofibrillar mitochondrial volume or mitochondrial respiratory function. As it relates to cross-bridge cycling, remarkably, the overt loss of contractile function in i-mPCKO muscle was reversed in permeabilized fibers supplied with exogenous Ca 2+ and ATP, with active tension being similar between i-mPCKO and WT mice, regardless of Ca 2+ concentration. Actin-myosin motility was also similar in skeletal muscle from i-mPCKO and WT mice. In conclusion, neither mitochondrial abundance/function, nor actomyosin cross-bridge cycling, are the underlying driver of contractile dysfunction in i-mPCKO mice. New & Noteworthy: The mechanism underlying dramatic loss of muscle contractile function with inducible deletion of both p300 and CBP in skeletal muscle remains unknown. Here we find that impairments in mitochondrial function or cross-bridge cycling are not the underlying mechanism of action. Future work will investigate other aspects of excitation-contraction coupling, such as Ca 2+ handling and membrane excitability, as contractile function could be rescued by permeabilizing skeletal muscle, which provides exogenous Ca 2+ and bypasses membrane depolarization.

CRISPR Activation Reverses Haploinsufficiency and Functional Deficits Caused by TTN Truncation Variants.

BACKGROUND: TTN truncation variants (TTNtvs) are the most common genetic lesion identified in individuals with dilated cardiomyopathy, a disease with high morbidity and mortality rates. TTNtvs reduce normal TTN (titin) protein levels, produce truncated proteins, and impair sarcomere content and function. Therapeutics targeting TTNtvs have been elusive because of the immense size of TTN, the rarity of specific TTNtvs, and incomplete knowledge of TTNtv pathogenicity. METHODS: We adapted CRISPR activation using dCas9-VPR to functionally interrogate TTNtv pathogenicity and develop a therapeutic in human cardiomyocytes and 3-dimensional cardiac microtissues engineered from induced pluripotent stem cell models harboring a dilated cardiomyopathy-associated TTNtv. We performed guide RNA screening with custom TTN reporter assays, agarose gel electrophoresis to quantify TTN protein levels and isoforms, and RNA sequencing to identify molecular consequences of TTN activation. Cardiomyocyte epigenetic assays were also used to nominate DNA regulatory elements to enable cardiomyocyte-specific TTN activation. RESULTS: CRISPR activation of TTN using single guide RNAs targeting either the TTN promoter or regulatory elements in spatial proximity to the TTN promoter through 3-dimensional chromatin interactions rescued TTN protein deficits disturbed by TTNtvs. Increasing TTN protein levels normalized sarcomere content and contractile function despite increasing truncated TTN protein. In addition to TTN transcripts, CRISPR activation also increased levels of myofibril assembly-related and sarcomere-related transcripts. CONCLUSIONS: TTN CRISPR activation rescued TTNtv-related functional deficits despite increasing truncated TTN levels, which provides evidence to support haploinsufficiency as a relevant genetic mechanism underlying heterozygous TTNtvs. CRISPR activation could be developed as a therapeutic to treat a large proportion of TTNtvs.

Multiscale biophysical models of cardiomyopathies reveal complexities challenging existing dogmas.

Mutations in sarcomeric proteins, including myosin, cause a variety of cardiomyopathies. A prominent hypothesis has been that myosin mutations causing hypercontractility of the motor lead to hypertrophic cardiomyopathy, while those causing hypocontractility lead to dilated cardiomyopathy; however, recent biophysical studies using multiscale computational and experimental models have revealed complexities not captured by this hypothesis. We summarize recent publications in Biophysical Journal challenging this dogma and highlighting the need for multiscale modeling of these complex diseases.

Divergent Molecular Phenotypes in Point Mutations at the Same Residue in Beta-Myosin Heavy Chain Lead to Distinct Cardiomyopathies.

In genetic cardiomyopathies, a frequently described phenomenon is how similar mutations in one protein can lead to discrete clinical phenotypes. One example is illustrated by two mutations in beta myosin heavy chain (ß-MHC) that are linked to hypertrophic cardiomyopathy (HCM) (Ile467Val, I467V) and left ventricular non-compaction (LVNC) (Ile467Thr, I467T). To investigate how these missense mutations lead to independent diseases, we studied the molecular effects of each mutation using recombinant human ß-MHC Subfragment 1 (S1) in in vitro assays. Both HCM-I467V and LVNC-I467T S1 mutations exhibited similar mechanochemical function, including unchanged ATPase and enhanced actin velocity but had opposing effects on the super-relaxed (SRX) state of myosin. HCM-I467V S1 showed a small reduction in the SRX state, shifting myosin to a more actin-available state that may lead to the "gain-of-function" phenotype commonly described in HCM. In contrast, LVNC-I467T significantly increased the population of myosin in the ultra-slow SRX state. Interestingly, molecular dynamics simulations reveal that I467T allosterically disrupts interactions between ADP and the nucleotide-binding pocket, which may result in an increased ADP release rate. This predicted change in ADP release rate may define the enhanced actin velocity measured in LVNC-I467T, but also describe the uncoupled mechanochemical function for this mutation where the enhanced ADP release rate may be sufficient to offset the increased SRX population of myosin. These contrasting molecular effects may lead to contractile dysregulation that initiates LVNC-associated signaling pathways that progress the phenotype. Together, analysis of these mutations provides evidence that phenotypic complexity originates at the molecular level and is critical to understanding disease progression and developing therapies.

Single-molecule mechanics and kinetics of cardiac myosin interacting with regulated thin filaments.

The cardiac cycle is a tightly regulated process wherein the heart generates force to pump blood to the body during systole and then relaxes during diastole. Disruption of this finely tuned cycle can lead to a range of diseases including cardiomyopathies and heart failure. Cardiac contraction is driven by the molecular motor myosin, which pulls regulated thin filaments in a calcium-dependent manner. In some muscle and nonmuscle myosins, regulatory proteins on actin tune the kinetics, mechanics, and load dependence of the myosin working stroke; however, it is not well understood whether or how thin-filament regulatory proteins tune the mechanics of the cardiac myosin motor. To address this critical gap in knowledge, we used single-molecule techniques to measure the kinetics and mechanics of the substeps of the cardiac myosin working stroke in the presence and absence of thin filament regulatory proteins. We found that regulatory proteins gate the calcium-dependent interactions between myosin and the thin filament. At physiologically relevant ATP concentrations, cardiac myosin's mechanics and unloaded kinetics are not affected by thin-filament regulatory proteins. We also measured the load-dependent kinetics of cardiac myosin at physiologically relevant ATP concentrations using an isometric optical clamp, and we found that thin-filament regulatory proteins do not affect either the identity or magnitude of myosin's primary load-dependent transition. Interestingly, at low ATP concentrations at both saturating and physiologically relevant subsaturating calcium concentrations, thin-filament regulatory proteins have a small effect on actomyosin dissociation kinetics, suggesting a mechanism beyond simple steric blocking. These results have important implications for the modeling of cardiac physiology and diseases.

Distinct effects of two hearing loss-associated mutations in the sarcomeric myosin MYH7b.

For decades, sarcomeric myosin heavy chain proteins were assumed to be restricted to striated muscle where they function as molecular motors that contract muscle. However, MYH7b, an evolutionarily ancient member of this myosin family, has been detected in mammalian nonmuscle tissues, and mutations in MYH7b are linked to hereditary hearing loss in compound heterozygous patients. These mutations are the first associated with hearing loss rather than a muscle pathology, and because there are no homologous mutations in other myosin isoforms, their functional effects were unknown. We generated recombinant human MYH7b harboring the D515N or R1651Q hearing loss-associated mutation and studied their effects on motor activity and structural and assembly properties, respectively. The D515N mutation had no effect on steady-state actin-activated ATPase rate or load-dependent detachment kinetics but increased actin sliding velocity because of an increased displacement during the myosin working stroke. Furthermore, we found that the D515N mutation caused an increase in the proportion of myosin heads that occupy the disordered-relaxed state, meaning more myosin heads are available to interact with actin. Although we found no impact of the R1651Q mutation on myosin rod secondary structure or solubility, we observed a striking aggregation phenotype when this mutation was introduced into nonmuscle cells. Our results suggest that each mutation independently affects MYH7b function and structure. Together, these results provide the foundation for further study of a role for MYH7b outside the sarcomere.

Functional assays reveal the pathogenic mechanism of a de novo tropomyosin variant identified in patient with dilated cardiomyopathy.

Dilated cardiomyopathy (DCM) is a leading cause of heart failure and a major indicator for heart transplant. Human genetic studies have identified over a thousand causal mutations for DCM in genes involved in a variety of cellular processes, including sarcomeric contraction. A substantial clinical challenge is determining the pathogenicity of novel variants in disease-associated genes. This challenge of connecting genotype and phenotype has frustrated attempts to develop effective, mechanism-based treatments for patients. Here, we identified a de novo mutation (T237S) in TPM1, the gene that encodes the thin filament protein tropomyosin, in a patient with DCM and conducted in vitro experiments to characterize the pathogenicity of this novel variant. We expressed recombinant mutant protein, reconstituted it into thin filaments, and examined the effects of the mutation on thin filament function. We show that the mutation reduces the calcium sensitivity of thin filament activation, as previously seen for known pathogenic mutations. Mechanistically, this shift is due to mutation-induced changes in tropomyosin positioning along the thin filament. We demonstrate that the thin filament activator omecamtiv mecarbil restores the calcium sensitivity of thin filaments regulated by the mutant tropomyosin, which lays the foundation for additional experiments to explore the therapeutic potential of this drug for patients harboring the T237S mutation. Taken together, our results suggest that the TPM1 T237S mutation is likely pathogenic and demonstrate how functional in vitro characterization of pathogenic protein variants in the lab might guide precision medicine in the clinic.

Drug specificity and affinity are encoded in the probability of cryptic pocket opening in myosin motor domains.

The design of compounds that can discriminate between closely related target proteins remains a central challenge in drug discovery. Specific therapeutics targeting the highly conserved myosin motor family are urgently needed as mutations in at least six of its members cause numerous diseases. Allosteric modulators, like the myosin-II inhibitor blebbistatin, are a promising means to achieve specificity. However, it remains unclear why blebbistatin inhibits myosin-II motors with different potencies given that it binds at a highly conserved pocket that is always closed in blebbistatin-free experimental structures. We hypothesized that the probability of pocket opening is an important determinant of the potency of compounds like blebbistatin. To test this hypothesis, we used Markov state models (MSMs) built from over 2 ms of aggregate molecular dynamics simulations with explicit solvent. We find that blebbistatin's binding pocket readily opens in simulations of blebbistatin-sensitive myosin isoforms. Comparing these conformational ensembles reveals that the probability of pocket opening correctly identifies which isoforms are most sensitive to blebbistatin inhibition and that docking against MSMs quantitatively predicts blebbistatin binding affinities (R2=0.82). In a blind prediction for an isoform (Myh7b) whose blebbistatin sensitivity was unknown, we find good agreement between predicted and measured IC50s (0.67 µM vs. 0.36 µM). Therefore, we expect this framework to be useful for the development of novel specific drugs across numerous protein targets.

Single Molecule Mechanics and Kinetics of Cardiac Myosin Interacting with Regulated Thin Filaments.

The cardiac cycle is a tightly regulated process wherein the heart generates force to pump blood to the body during systole and then relaxes during diastole. Disruption of this finely tuned cycle can lead to a range of diseases including cardiomyopathies and heart failure. Cardiac contraction is driven by the molecular motor myosin, which pulls regulated thin filaments in a calcium-dependent manner. In some muscle and non-muscle myosins, regulatory proteins on actin tune the kinetics, mechanics, and load dependence of the myosin working stroke; however, it is not well understood whether or how thin filament regulatory proteins tune the mechanics of the cardiac myosin motor. To address this critical gap in knowledge, we used single-molecule techniques to measure the kinetics and mechanics of the substeps of the cardiac myosin working stroke in the presence and absence of thin filament regulatory proteins. We found that regulatory proteins gate the calcium-dependent interactions between myosin and the thin filament. At physiologically relevant ATP concentrations, cardiac myosin's mechanics and unloaded kinetics are not affected by thin filament regulatory proteins. We also measured the load-dependent kinetics of cardiac myosin at physiologically relevant ATP concentrations using an isometric optical clamp, and we found that thin filament regulatory proteins do not affect either the identity or magnitude of myosin's primary load-dependent transition. Interestingly, at low ATP concentrations, thin filament regulatory proteins have a small effect on actomyosin dissociation kinetics, suggesting a mechanism beyond simple steric blocking. These results have important implications for both disease modeling and computational models of muscle contraction. Significance Statement: Human heart contraction is powered by the molecular motor ß-cardiac myosin, which pulls on thin filaments consisting of actin and the regulatory proteins troponin and tropomyosin. In some muscle and non-muscle systems, these regulatory proteins tune the kinetics, mechanics, and load dependence of the myosin working stroke. Despite having a central role in health and disease, it is not well understood whether the mechanics or kinetics of ß-cardiac myosin are affected by regulatory proteins. We show that regulatory proteins do not affect the mechanics or load-dependent kinetics of the working stroke at physiologically relevant ATP concentrations; however, they can affect the kinetics at low ATP concentrations, suggesting a mechanism beyond simple steric blocking. This has important implications for modeling of cardiac physiology and diseases.

Functional divergence of the sarcomeric myosin, MYH7b, supports species-specific biological roles.

Myosin heavy chain 7b (MYH7b) is an evolutionarily ancient member of the sarcomeric myosin family, which typically supports striated muscle function. However, in mammals, alternative splicing prevents MYH7b protein production in cardiac and most skeletal muscles and limits expression to a subset of specialized muscles and certain nonmuscle environments. In contrast, MYH7b protein is abundant in python cardiac and skeletal muscles. Although the MYH7b expression pattern diverges in mammals versus reptiles, MYH7b shares high sequence identity across species. So, it remains unclear how mammalian MYH7b function may differ from that of other sarcomeric myosins and whether human and python MYH7b motor functions diverge as their expression patterns suggest. Thus, we generated recombinant human and python MYH7b protein and measured their motor properties to investigate any species-specific differences in activity. Our results reveal that despite having similar working strokes, the MYH7b isoforms have slower actin-activated ATPase cycles and actin sliding velocities than human cardiac ß-MyHC. Furthermore, python MYH7b is tuned to have slower motor activity than human MYH7b because of slower kinetics of the chemomechanical cycle. We found that the MYH7b isoforms adopt a higher proportion of myosin heads in the ultraslow, super-relaxed state compared with human cardiac ß-MyHC. These findings are supported by molecular dynamics simulations that predict MYH7b preferentially occupies myosin active site conformations similar to those observed in the structurally inactive state. Together, these results suggest that MYH7b is specialized for slow and energy-conserving motor activity and that differential tuning of MYH7b orthologs contributes to species-specific biological roles.