Regulation of the crossbridge cycle that drives muscle contraction involves a reconfiguration of the troponin–tropomyosin complex on actin filaments. By comparing atomic models of troponin–tropomyosin fitted to cryo-EM structures of inhibited and Ca2+-activated thin filaments, we find that tropomyosin pivots rather than rolls or slides across actin as generally thought. We propose that pivoting can account for the Ca2+ activation that initiates muscle contraction and then relaxation influenced by troponin-I (TnI). Tropomyosin is well-known to occupy either of three meta-stable configurations on actin, regulating access of myosin motorheads to their actin-binding sites and thus the crossbridge cycle. At low Ca2+ concentrations, tropomyosin is trapped by TnI in an inhibitory B-state that sterically blocks myosin binding to actin, leading to muscle relaxation. Ca2+ binding to TnC draws TnI away from tropomyosin, while tropomyosin moves to a C-state location over actin. This partially relieves the steric inhibition and allows weak binding of myosin heads to actin, which then transition to strong actin-bound configurations, fully activating the thin filament. Nevertheless, the reconfiguration that accompanies the initial Ca2+-sensitive B-state/C-state shift in troponin–tropomyosin on actin remains uncertain and at best is described by moderate-resolution cryo-EM reconstructions. Our recent computational studies indicate that intermolecular residue-to-residue salt-bridge linkage between actin and tropomyosin is indistinguishable in B- and C-state thin filament configurations. We show here that tropomyosin can pivot about relatively fixed points on actin to accompany B-state/C-state structural transitions. We argue that at low Ca2+ concentrations C-terminal TnI domains attract tropomyosin, causing it to bend and then pivot toward the TnI, thus blocking myosin binding and contraction.
Introduction
Head-to-tail-linked tropomyosin coiled coils form cables that span the length of actin-containing thin filaments in virtually all animal cells (Brown and Cohen, 2005; Hitchcock-DeGregori, 2008; Holmes and Lehman, 2008). The positioning of the tropomyosin cables across the surface of actin is not fixed, and in striated muscle fibers is defined by troponin, Ca2+ levels, and myosin-head interaction (Gordon et al., 2000; Holmes and Lehman, 2008; Lehman, 2016). By occupying one of three meta-stable configurations, tropomyosin is thought to influence the access of myosin motorheads to their actin-binding sites (Gordon et al., 2000; Poole et al., 2006; Geeves, 2012; Lehman, 2016). Accordingly, tropomyosin controls muscle contractility either by inhibiting or facilitating the myosin crossbridge cycle on actin. At low intracellular Ca2+ concentrations in resting, relaxed muscles, tropomyosin is thought to be trapped by the C-terminal domain of troponin-I (TnI) subunit of the troponin complex in an inhibitory B-state that blocks the myosin head binding to actin and therefore the initiation of the crossbridge cycle. Ca2+ binding to troponin-C (TnC) of the complex draws the inhibitory C-terminal domain of TnI away from actin–tropomyosin to relocate onto the TnC N-lobe domain (Marston and Zamora, 2020; Yamada et al., 2020; Lehman et al., 2021). Thus, the reversible binding and translocation of this relatively short domain between TnC and actin–tropomyosin is a lynchpin in the muscle regulatory process (Ramos, 1999; Blumenschein et al., 2006; Galińska-Rakoczy et al., 2008; Sheng and Jin, 2016; Wong et al., 2019; Hornos et al., 2021; Tobacman and Cammarato, 2021). The change in TnI association is accompanied by an ∼15 Å movement of tropomyosin to a C-state location over actin, relieving the steric inhibition and allowing weak binding of myosin heads to the now partially exposed actin-binding interface (Lehman et al., 1994; Gordon et al., 2000; Poole et al., 2006; Lehman, 2016; Risi et al., 2017; Risi et al., 2021a; Yamada et al., 2020). Unfettered by TnI, the tropomyosin cable is thought to oscillate about the C-state position on actin with patches of weakly bound myosin heads then transitioning to strongly actin-bound configurations (Maytum et al., 2008; Geeves, 2012). In turn, myosin binding causes a further tropomyosin movement another ∼25 Å away from the B-state to an M-state position. The open M-state configuration of tropomyosin then facilitates the binding of neighboring myosin heads, thereby fully activating the thin filament (Vibert et al., 1997; Behrmann et al., 2012; Lehman, 2016). This stepwise forward and reverse process is integral to the “Three-State Steric Model” that has been proposed to explain the on–off switching mechanism governing all striated muscles (McKillop and Geeves, 1993; Geeves, 2012; Lehman, 2017). At the same time, crossbridge cycling on thin filaments is additionally influenced by varying myosin head recruitment from thick filaments, possibly modulated by phosphorylation, strain, or the direct effects of Ca2+ (Lehman, 1978; Linari et al., 2015; Fusi et al., 2016; Kampourakis et al., 2016; Craig and Padrón, 2022; Ma et al., 2022). In this article, we concentrate on Ca2+-sensitive gating of myosin binding sites on actin by troponin–tropomyosin.
Despite the attractive simplicity of the steric model’s description of thin filament regulation, the exact structural changes that are integral to the Ca2+-sensitive B-state/C-state shifts in troponin and tropomyosin on actin remain uncertain. For example, precise delineation of B- and C-state tropomyosin positions on actin are difficult to define at atomic resolution because corresponding interactions between thin filament components are dynamic and occur over a huge structural scale (discussed in Rynkiewicz et al., 2022). At best, the transitions between tropomyosin, troponin, and actin are revealed as static moderate-resolution cryo-EM end-points (Yamada et al., 2020; Risi et al., 2021a). In contrast, the solution of the larger myosin-induced movement of the tropomyosin to the M-state position on actin appears to be less problematic (Doran et al., 2020; Doran et al., 2023a; Baldo et al., 2021). In this case, the reconfiguration is captured as distinct and well-defined M-state end-points (Milligan et al., 1990; Vibert et al., 1997; Poole et al., 2006; Behrmann et al., 2012; Doran et al., 2020; Doran et al., 2023a; Doran et al., 2023b; Risi et al., 2021b).
Tropomyosin is ideally adapted to interact end-to-end and continuously along actin filaments (Hitchcock-DeGregori, 2008; Holmes and Lehman, 2008; Orzechowski et al., 2014). The N- and C-terminal ends of coiled-coil tropomyosin are specialized to form short four-helix bundles with each other (Greenfield et al., 2006; Frey et al., 2010; Orzechowski et al., 2014). The resulting head-to-tail overlapping interactions thereby produce unbroken, longitudinally oriented tropomyosin cables traversing thin filaments. Moreover, tropomyosin is a modular molecule and each of its seven pseudorepeats interacts sequentially with a neighboring actin subunit that sits along the thin filament actin helix (McLachlan and Stewart, 1976; Brown and Cohen, 2005; Li et al., 2011; Orzechowski et al., 2014). However, as mentioned, the azimuthal position of the tropomyosin cables is not fixed. Indeed, the alternative azimuthal positional states assumed by the cable and spanning the length of the filament are thought to be governed by clustered electrostatic interactions between actin subunits and tropomyosin pseudorepeats (Brown and Cohen, 2005; Holmes and Lehman, 2008; Li et al., 2011; Pavadai et al., 2020a; Pavadai et al., 2020b), including residues Lys 326 and Lys 328 on actin subunits and glutamates and aspartates on each tropomyosin pseudorepeat (Brown et al., 2005; Li et al., 2011; Pavadai et al., 2020a; Pavadai et al., 2020b; Rynkiewicz et al., 2022). These positional states are separated from each other by low-energy barriers and thus are likely to easily alternate amongst themselves unless influenced further by troponin or myosin (Lehman et al., 2000; Pirani et al., 2005; Li et al., 2011; Risi et al., 2021a).
The relative alignment of actin and tropomyosin pseudorepeats and the azimuthal locations of the tropomyosin coiled-coil on thin filaments indicate that intermolecular residue-to-residue salt bridge formation between the two proteins is much the same in B- and C-state thin filament configurations (Pavadai et al., 2020a; Pavadai et al., 2020b; Rynkiewicz et al., 2022; Table 1). Considering this equivalence, how then is it possible for tropomyosin to reposition, let alone slide or roll over actin between B- and C-states, as has been variously proposed (McLachlan and Stewart, 1976; Holthauzen et al., 2004; Li et al., 2010; Behrmann et al., 2012; von der Ecken et al., 2015; Risi et al., 2017; Lehman et al., 2020)? We suggest here that neither sliding nor rolling motions describe B- to C-state tropomyosin transitions accurately; instead, we propose that central domains of tropomyosin molecules pivot back and forth about relatively fixed points near their head-to-tail junctions, and that this pivoting accounts for the apparent B- to C-state positional shift. In effect, pivoting changes the average center-of-mass of the tropomyosin on actin and thus the azimuthal position of tropomyosin’s superhelical axis in B- and C-states but does not alter the binding site linkages between tropomyosin and actin. We argue that at resting low-Ca2+ concentration, C-terminal domains of TnI attract the edge of the tropomyosin coiled-coil which then bends toward TnI. This motion produces the B-state filament configuration without disrupting contact between tropomyosin and canonical actin binding sites. We expect that muscle activation and ensuing C-terminal TnI detachment from actin–tropomyosin involves tropomyosin pivoting back to its C-state orientation. Our view is supported by new computational studies that complement existing cryo-EM reconstructions.
Materials and methods
Atomic models matched to cryo-EM volumes
B- and C-state atomic models of the human cardiac thin filament (PDB IDs 7UTL and 7UTI, respectively) were taken from Rynkiewicz et al. (2022) and used as reference starting structures for analysis. To date, these models most accurately and precisely fit to actin and tropomyosin contained in cryo-EM reconstructions of thin filament preserved under low- and high-Ca2+ conditions (Yamada et al., 2020; Risi et al., 2022; see Fig. S1). Generation of these models involved molecular dynamics flexible fitting routines (MDFF; Trabuco et al., 2009) and extensive energy minimization (Liebschner et al., 2019) of atomic models of tropomyosin, the troponin core-domain, TnT1, and the C-terminal TnI when bound to actin–tropomyosin and fitted to the cryo-EM structures (as described in Rynkiewicz et al. [2022] and below). The resulting 7UTI and 7UTL atomic models are representations only of structures that are resolved in the cryo-EM reconstructions as noted by Yamada et al. (2020) and Risi et al. (2021a).
MD simulation of 7UTI and 7UTL models
Since the 7UTI and 7UTL models do not contain full 284-residue tropomyosin chains, the tropomyosin was built out and the model extended to include 28 actin subunits, 8 tropomyosin coiled coils, and 4 troponin complexes in a periodic system, which was then used for MD simulations performed in an explicit solvent in the presence of 150 mM NaCl and 3 mM MgCl2. The periodic boundary system was chosen to simulate what in essence is an infinitely long filament. First, the PDB structures were translated to the origin and aligned along the z-axis. Next, a copy of the structure was moved up to extend the filament by aligning actin chain B onto chain P (7UTL) or chain D onto chain R (7UTI). For the periodic boundary, the helical rotation of the actin must be 167.1° rather than the 166.5° observed in the PDB structure. To relieve any strain on the filament caused by forcing the filament into the periodic system, starting outside the tropomyosin head-to-tail overlap junction and proceeding in each direction along the filament sequentially by pseudorepeat, the actins with their associated tropomyosin/troponin chains were rotated 1°. This results in the outer pseudorepeats being rotated a total of 6°. At this point, the tropomyosin chains from each copy were joined together at the closest backbone positions to create final full-length tropomyosin chains consisting of residues 1–36 of chains W and Z and residues 37–284 of chains a and b or residues 1–72 of chains g and h and residues 73–284 of chains i and j for 7UTL. For 7UTI, the full tropomyosin consisted of residues 1–42 of chains W and b and residues 43–284 of chains e and f or residues 1–53 of chains g and h and residues 54–284 of chains i and j. The full-length chains were then minimized locally using α-helical constraints obtained by running the program Coot (Emsley et al., 2010) at residues 34–39, 70–75, 40–45, 51–57, 61–66, 103–108, 140–145, 185–190, and 221–226 to fix any poor bonds created by the above manipulations of the structure.
This full tropomyosin and the associated actin and troponin chains were then copied using the actin helical symmetry to create the full filament. Before adding solvent molecules, an additional four actin subunits, as well as tropomyosin and troponin T1 chains were added to the ends of the model to serve as placeholders for the periodic boundary. Solvent atoms were then added using the solvate plugin in the graphics program VMD (Humphrey et al., 1996) with a final box size of 200 Å × 200 Å × 766.196 Å for 7UTL and 200 Å × 200 Å × 764.718 Å for 7UTI. After the removal of the placeholder chains, the system was neutralized and salt concentrations were adjusted to 0.15 M NaCl and 3 mM MgCl2 using the ionize plugin in VMD (Humphrey et al., 1996). Before the simulation, the system was minimized in the CHARMM36m force field (Huang et al., 2017) and a non-bonded cutoff of 14 Å with switching using NAMD version 2.13 tools (Phillips et al., 2005; Phillips et al., 2020) in two steps. First, the solvent molecules were minimized and then the whole system was minimized with a slow release of harmonic constraints to the starting coordinates. Next, the system was heated to 300 K at constant volume with constraints to the minimized coordinates followed by a slow release of the constraints at constant temperature and volume. At this point, the pressure in the system was brought up to about 1 bar using the Nosé-Hoover Langevin piston, which was maintained through the remainder of the simulation. The simulation was then allowed to proceed without constraints. For analysis, each of the four copies of tropomyosin/troponin in each model was treated as an independent replica.
MD simulation of thin filament segments containing C-terminal TnI and tropomyosin pseudorepeats 3–5
MD was also run on a short three-actin long thin filament model containing the 7UTL B-state configuration of the C-terminal domain of TnI (residues 155–210) on actin, but here superposed on its cognate, the three-actin long C-state actin–tropomyosin structure. The corresponding three pseudorepeat long tropomyosin included residues 90–192, thus the truncated tropomyosin lacked its head-to-tail overlap and interactions with TnT1. To prevent the ends of the tropomyosin helices from unwinding during the simulation, the backbone φ and ψ angles of the first three residues and the backbone O-NH n,n+4 hydrogen bond distances over the first five residues were constrained to known α-helical values. To prevent the tropomyosin segment from dissociating from actin during MD, a harmonic distance constraint was imposed to maintain side chains between actin residues K326 and K328 and eight of the closest neighboring glutamates on tropomyosin within a 0.8 Å standard deviation distribution of their 3.9 Å average tip-to-tip distances. This procedure kept tropomyosin at 37.1 ± 0.8 Å radius while preventing wobble about the tropomyosin superhelix along actin. However, tropomyosin was free to move azimuthally as a unit and to bend locally; thus, movement of the tropomyosin toward TnI could be assessed. The system was minimized and MD was carried out as above in explicit solvent.
Other methods
Interaction energy calculations were performed using the NAMD Energy plugin in VMD (Humphrey et al., 1996). All calculations and simulations used the CHARMM36m force field (Huang et al., 2017). Tropomyosin’s helical parameters were calculated using the program, TWISTER (Strelkov and Burkhard, 2002). Molecular graphics were carried out using UCSF-Chimera and VMD programs (Pettersen et al., 2004; Humphrey et al., 1996).
Online supplemental material
Fig. S1 illustrates the fitting of C-state and B-state atomic models of the human cardiac thin filament to matching high- and low-Ca2+ filament reconstructions. Fig. S2 compares the interaction energetics between residues of actin and tropomyosin along the above B-state and C-state atomic models as well as the corresponding energetics during MD simulations carried out on these models. Videos 1 and 2 show the movement of C-state tropomyosin toward TnI during MD.
Results
Approach
As mentioned, steric regulation of muscle contraction involves a sequential process in which thin filaments are activated in at least two steps. Thin filaments first transition from a blocked, low-Ca2+ B-state characteristic of relaxed muscles to an intermediate Ca2+-induced C-state and then to a final myosin-induced M-state characteristic of actively contracting muscle. We previously examined the final myosin-induced, M-state tropomyosin transitions on actin at atomic level resolution (Kiani et al, 2019; Doran et al, 2020, 2023a; Doran and Lehman, 2021). In the current analysis, we dissect the initial Ca2+-dependent tropomyosin transition by comparing the organization of low-Ca2+ B-state and high-Ca2+ C-state thin filaments at near-atomic resolution. Our objective is to better define how TnI mediates tropomyosin’s steric interference of myosin binding sites on actin in resting muscle and then how such steric hindrance is reversed in contracting muscle. Ever-improved atomic models defined by cryo-EM reconstructions of cardiac thin filaments now provide an opportunity to determine tropomyosin organization, interactions, energetics, and dynamics on actin in low- and high-Ca2+ states (Yamada et al., 2020; Pavadai et al., 2020a; Pavadai et al., 2020b; Risi et al., 2021a; Rynkiewicz et al., 2022). By examining tropomyosin configurations with a precision previously lacking, we can disambiguate prior confusion about tropomyosin rolling or sliding between B- and C-regulatory states on actin. As we report, tropomyosin appears to do neither, rather it reconfigures locally by TnI-induced pivoting.
Tropomyosin asymmetries
When viewed from the surface of each actin subunit, tropomyosin coiled coils display seven broad and narrow motifs that reflect the twisting of their two-chained coiled-coil helices (Fig. 1, A and B). This face-on appearance replicates that of the wide and tapered tropomyosin densities seen in surface views of F-actin–tropomyosin reconstructions (see Fig. S1). Each of the seven modular tropomyosin pseudorepeats is associated with a neighboring actin subunit along the thin filament, suggesting perfect symmetry. However, individual pseudorepeats are not identical copies of each other (Brown et al., 2005; Hitchcock-DeGregori, 2008; Li et al., 2010; Pavadai et al., 2020a; Pavadai et al., 2020b; hence the term pseudorepeat). Thus, the potential regional steric effects on myosin binding may vary from one actin subunit to the next. Moreover, the tropomyosin coiled coil is thought to behave as a dynamic semiflexible rod (Li et al., 2010) creating additional positional uncertainty. Any deviation from strict superhelical symmetry no doubt is also influenced by tropomyosin form-fitting to its actin substrate. For example, idealized axial projection of B- and C-state tropomyosin superhelices along actin may not account for subtle pseudorepeat non-uniformity or capture additional azimuthal variance associated with the effects of troponin on tropomyosin tracks along actin (Fig. 1). It is very likely that current comparison of B- and C-state actin–tropomyosin conformation provides only partial cues needed to characterize the roles of tropomyosin and troponin in the steric regulatory mechanism (Tobacman, 2021; Rynkiewicz et al, 2022). To fully elucidate this regulatory process, we suggest that the local and aggregate variance in tropomyosin position along with actin and the corresponding effect of TnI on B-state tropomyosin must be rigorously and systematically evaluated. We find that this can be best achieved by comparison of recent atomic-resolution B- and C-state filament models as described below.
Mapping tropomyosin in the thin filament models
Atomic models of B- and C-state filaments (PDB IDs 7UTI and 7UTL) were generated to accurately fit to the Yamada et al. (2020) cryo-EM reconstructions of cardiac thin filaments preserved under low- and high-Ca2+ conditions (EMDB 0728 and 0729, respectively; Fig. 1, A and B). These refined PDB coordinates (Rynkiewicz et al, 2022) were chosen over those previously reported by Yamada et al. (2020) and subsequently employed by Risi et al. (2022; PDB IDs 6KN7 and 6KN8) since the longitudinal alignment of tropomyosin on actin was incorrectly calibrated in the earlier reports (see Pavadai et al, 2020a, 2020b; Baldo et al, 2021; Rynkiewicz et al, 2022). The revised 7UTI and 7UTL coordinates and corresponding refined tropomyosin:actin geometry define residue-to-residue contacts most relevant to elucidating the mode of tropomyosin motions between B- and C-state configurations; they do not include hypothesized representations of unresolved densities in the EM maps.
Comparison of the updated B- and C-state models clearly confirms the well-recognized difference in the azimuthal positioning of tropomyosin on actin (Fig. 1, C and D), as originally reported (Lehman et al., 1994; Vibert et al., 1997; Poole et al., 2006; Yamada et al., 2020; Risi et al., 2021a). The difference is greatest between the two tropomyosin cables over pseudorepeats 3–5 (cf. Yamada et al., 2020; Pavadai et al., 2020a; Pavadai et al., 2020b). Here, B-state pseudorepeats approach the extended form of the C-terminal domain of TnI that is released from TnC at low Ca2+ levels and runs roughly parallel to tropomyosin, as described by Yamada et al. (2020) and Lehman et al. (2021). Thus, tropomyosin seemingly is drawn toward the TnI domain producing the well-recognized regional azimuthal shift (Fig. 1 C; and Fig. 2, C–G).
Superposition of the new models shows that the superhelical paths formed by B- and C-state tropomyosin are separated by ∼9 to 16 Å midway along the above-described region of the tropomyosin molecule (Fig. 1 E). When side chains are included in the calculations, the outer edge of B-state tropomyosin and the inner edge of C-state tropomyosin are separated by up to 30 Å. However, our calculations show that the distinction between B- and C-state tropomyosin paths is below 2.4 Å between residues 248 and 265 (Fig. 1 E). This proximity occurs close to the C-terminal end of tropomyosin as the chains approach the head-to-tail domain. It is in this region that tropomyosin colocalizes and binds strongly to helical domains of TnT (Pavadai et al., 2020a; Pavadai et al., 2020b; Risi et al., 2022), possibly serving to pin both B- and C-state tropomyosin chains on actin at a common site.
Tropomyosin pivoting
As described above, the superposition of the B- and C-state models clearly shows the component tropomyosin strands running in parallel as they wrap around the actin filament helix. Both B- and C-regulatory state maps show that both use the same set of 10 acidic residues distributed along tropomyosin to bind to actin. These residues interact closely with Lys 328 present on each successive actin subunit of the thin filaments (Fig. 1, A and B). The salt-bridge interactions formed between tropomyosin and actin are documented in Table 1 and Fig. S2).
Transverse sections of the superposed models further illustrate periodic interactions between tropomyosin and actin as well as confirm the distinct and varying B- and C-state configurations assumed by tropomyosin on actin (Figs. 2 and 3). When sectioning through the superposed maps, a conspicuous and complementary pattern of tropomyosin B/C-state profiles repeats itself roughly every 40 residues, i.e., occurs once for every tropomyosin pseudorepeat. The pattern is most obvious at points of superposition close to the above-mentioned actin and TnI binding sites where one of the two α-helical chains from the B-state tropomyosin and one of the two helices from the C-state coiled coil overlap and appear as concentric rings in the superposed sections (Fig. 2, C–G; and Fig. 3, A–F). In marked contrast to such juxtaposition, the other non-overlapping pair of chains diverge and orient themselves in opposite directions away from the well-matched pair (Figs. 2 and 3). Here, the latter B-state tropomyosin chain is angled toward the outer domain of actin, while the C-state chain is oriented toward the inner domain of actin (Figs. 2 and 3) as if the tropomyosin coiled-coil bends in the B-state by forming attractive electrostatic linkages and hydrophobic interactions with TnI (best visualized in Fig. 2, C–G). Not surprisingly, the magnitude of the bending of the B-state tropomyosin coiled-coil away from the C-state copy is greatest for tropomyosin residues contained in pseudorepeats 3 to 5 that lie next to the TnI C-terminal extension on B-state filaments (Figs. 1, 2, and 3). It therefore is likely that the B-state to C-state reorientation is caused locally by C-terminal TnI attachment and detachment to and from tropomyosin, as previously concluded (Yamada et al., 2020). Indeed, the bending appears to diminish toward pseudorepeats adjacent to the tropomyosin head-to-tail intermolecular junction, and, for example, is minimal over pseudorepeat seven as tropomyosin converges on the tropomyosin overlapping domain (Fig. 1), namely, where, as mentioned above, tropomyosin may also interact strongly with the TnT1 tail domain of TnT (Pavadai et al., 2020a; Pavadai et al., 2020b; Risi et al., 2022). With tropomyosin repositioning minimal near its head-to-tail overlap domain and greatest over more central domains (pseudorepeats 3–5), it follows that the B/C-state tropomyosin reconfiguration does not represent a global rolling or sliding motion of tropomyosin coiled-coil cables over actin but rather a pivoting between relatively fixed ends. Bear in mind that the tropomyosin pseudorepeat contact points on successive actin residues mentioned above to Lys 328 (Fig. 1, A–C; and Table 1) are unchanged in the B- and C-regulatory states. The observation is consistent with the view that these contacts are pivot point fulcrums for the tropomyosin motion (Table 1), as docking studies on actin–tropomyosin segments originally suggested (Pavadai et al., 2020a; Pavadai et al., 2020b; Doran et al., 2020).
Energetics of B- and C-state thin-filament interactions
We previously described an array of polar interactions between TnI and tropomyosin (Lehman et al., 2021) and now conclude that they likely facilitate tropomyosin’s pivoting between B- and C-state orientations. Indeed, measurement of thin-filament interaction energy that includes C-terminal TnI polar residues Gln 156, Arg 162, Arg 170, Lys 174, Glu 182, Arg 186, Arg 192, Asp 196, and Ser 199 increases the total interaction energy of tropomyosin and actin/troponin in B-state thin filaments and thereby favors the C-state to B-state transition (Table 2). In addition, non-polar contacts involving the TnI switch-peptide and the H4 helix also likely enhance B-state linkage with tropomyosin. In total these interactions, no doubt, are responsible for the B-state-directed tropomyosin pivoting initiated over tropomyosin pseudorepeats 3–5. Conversely, the energetic cost of pivoting tropomyosin toward the B-state position when the C-terminal TnI is sequestered by TnC at high Ca2+ concentration would diminish tropomyosin movement from its default C-state position.
MD simulations of B- and C-state thin filaments
In order to further evaluate tropomyosin organization and its motions on actin, MD simulations were performed on the 7UTI and 7UTL atomic models of C- and B-state thin filaments containing tropomyosin, the troponin core domain, TnT, and the C-terminal domain of TnI in the latter model. Tropomyosin was noted to oscillate locally about its two starting positions while maintaining its ∼37 Å radius and the canonical salt bridge pivot points between B/C-state helices along all actin subunits described above for the duration of the 120 ns MD simulation (Fig. 2, D–H; and Fig. S2). The stability of these interactions supports our contention that tropomyosin pivoting is a native process and not a cryo-EM artifact. We also noted that during MD simulations, side-chains of C-state tropomyosin Glu 114 and 139 appear to face and sample surrounding solvent, while in the B-state motions of E114 and E139 side-chains facing TnI are more limited due to contact with C-terminal TnI domain residues R192 and R170 (Fig. 2, D–H), supporting our notion that salt-bridge linkage, at least in part, drives complex formation between TnI and tropomyosin. In this regard, the interaction between tropomyosin and the C-terminal domain of TnI, for example via Arg 192, was also shown by NMR performed on thin-filament peptides (Wong et al., 2019; Hornos et al., 2021).
Examination of Fig. 2, D–H, suggests that during MD, a fraction of the side-chains extending from C-state tropomyosin residues 114 or 139 visit the region occupied by the B-state C-terminal TnI peptide as if primed to interact at the onset of relaxation (highlighted by a slanted arrow in Fig. 2, D–H). The size of the periodic boundary system used to mimic an infinitely long thin filament (∼3 million atoms) precluded capturing the global movement of C-state tropomyosin to its B-state position during MD. To visualize the potential for such a transition, we carried out MD simulations on a reduced size C/B-state chimera containing actin and tropomyosin pseudorepeats 3–5 in their C-state positions, as well as C-terminal residues 155–210 of TnI in its B-state position; thus, TnI was placed in a position set to draw C-state tropomyosin to its B-state position. The anticipated pivoting motion of tropomyosin occurred within 10–40 ns of the onset of MD, during which salt bridges formed between residues E114 and R192 as well as between E139 and R170 (Videos 1 and 2). The simulation suggests that tropomyosin can pivot as a global unit in the direction of the C-terminal domain of TnI. However, a full description of the interaction would require understanding the docking of C-terminal TnI to actin–tropomyosin once released from the N-lobe of TnC under low-Ca2+ muscle relaxing conditions. In this regard, several questions remain unanswered. For example, do the components of the TnI domain, namely its switch peptide, H4-helix, and unstructured tail-end bind to actin and tropomyosin sequentially, and in what order?
Discussion
X-ray fiber diffraction of intact tissue has previously been used to estimate the azimuthal shift of tropomyosin undergoing the B/C-state transition (Poole et al., 2006). The ∼23 Å (i.e., 15°) azimuthal movement that was calculated is considerably larger than has been observed by cryo-EM and validated here. However, contributions from the C-terminal extension of TnI (Yamada et al., 2020; Risi et al., 2021a) may not have been adequately delineated from those of tropomyosin in the fiber diffraction analysis. Thus, the inclusion of TnI contributions may have overstated the actual amount of tropomyosin translation. Similarly, unaccounted TnI contributions as well as the effects of sample flattening may also have exaggerated the amount of B/C-state tropomyosin relocation detected in reconstructions of negatively stained filaments (Lehman et al., 1994; Vibert et al., 1997; Pirani et al., 2005).
Our observations raise obvious questions about the origin and functional significance of the seemingly small B/C-state transitions, which we now can answer. The previously proposed, rather large 15° shift in tropomyosin position between low- and high-Ca2+ B- and C-state positions on actin provided a heuristic understanding of troponin–tropomyosin–based steric regulation of muscle contractility. However, we now posit that steric regulation primarily involves trapping the tropomyosin coiled-coil very close to the C-state over myosin binding sites on actin. As mentioned, the effect is brought about by interactions of tropomyosin with the C-terminal domain of TnI, and the resulting B/C-state pivoting is a manifestation of this. Such pivoting does not require or involve large-scale B/C-state translation of tropomyosin over the actin surface. TnI-induced pivoting does accentuate the magnitude of myosin binding site obstruction by drawing tropomyosin over the target site, but the action of TnI clasping onto tropomyosin and trapping the coiled-coil in an inhibitory state may be what is crucial.
It has long been known that tropomyosin can inhibit actomyosin ATPase in vitro (Lehrer and Morris, 1982, 1984) and by itself partially obstructs myosin binding sites on actin in the absence of troponin (Vibert et al., 1997; Lehman et al., 2000; Poole et al., 2006). However, troponin-free tropomyosin as well as C-state tropomyosin is thought to oscillate on actin (Pirani et al., 2005; Maytum et al., 2008). Thus, tropomyosin acting alone would be ineffective as an on–off switching mechanism, potentially allowing limited access for weakly binding myosin heads to thin filaments in B- or C-state configurations. This view is supported by observations made in vitro showing that the inhibitory effect of tropomyosin dissipates as myosin concentration is increased, enabling significant numbers of myosin heads to bind to actin–tropomyosin and transition to the M-state (Lehrer and Morris, 1982; Lehrer and Morris, 1984). Hence, TnI in the troponin-regulated thin filament is needed to augment tropomyosin-linked effects by stabilizing steric interference and preventing tropomyosin azimuthal drift. The combined troponin–tropomyosin action thus produces robust actomyosin inhibition and activation. Here, TnI is the choreographer directing the pivoting performed during the tropomyosin pirouette. Pivoting over a common set of residues is particularly advantageous for the B/C-state transition as it sidesteps the energetic cost involved in switching protein–protein interactions that tropomyosin sliding or rolling over actin would otherwise entail.
Changes observed in muscle fiber diffraction patterns 50 yr ago demonstrated tropomyosin’s regulatory movement over thin filament actin (Huxley, 1972; Haselgrove, 1972; Vibert et al., 1972; Parry and Squire, 1973). These seminal studies presaged the work reported here and no doubt are consistent with the troponin–tropomyosin regulatory mechanism that we propose.
Data availability
The data underlying Figs. 1, 2, and 3; Tables 1 and 2; and Fig. S1 and Fig. S2 are derived directly from PDB coordinates 7UTI and 7UTL deposited in the RCSB Protein Data Bank (https://www.rcsb.org/). Fig. S1 can be reproduced from the same PDB files and from coordinates deposited for cryo-EM volumes (EMDB IDs 0729, 0728, 22965, and 22964) in the Electron Microscopy Data Bank (https://www.ebi.ac.uk/emdb/) using the Fit-in-Map tool in the UCSF-Chimera suite of programs (Petterson et al., 2004).
Acknowledgments
Henk L. Granzier served as editor.
These studies were supported by National Institutes of Health grant R01HL036153 to W. Lehman. Computational work was carried out in-house by using resources provided by the Massachusetts Green High Performance Computing Center.
Author contributions: W. Lehman thought of the general approach taken, and along with M.J. Rynkiewicz developed the principal concepts presented. M.J. Rynkiewicz performed MD simulations and calculated interaction energies. W. Lehman and M.J. Rynkiewicz prepared figures and tables. W. Lehman wrote the manuscript.
References
Author notes
The authors declare no competing interests exist.