Myofilament-associated proteins with intrinsic disorder (MAPIDs) and their resolution by computational modeling.

The cardiac sarcomere is a cellular structure in the heart that enables muscle cells to contract. Dozens of proteins belong to the cardiac sarcomere, which work in tandem to generate force and adapt to demands on cardiac output. Intriguingly, the majority of these proteins have significant intrinsic disorder that contributes to their functions, yet the biophysics of these intrinsically disordered regions (IDRs) have been characterized in limited detail. In this review, we first enumerate these myofilament-associated proteins with intrinsic disorder (MAPIDs) and recent biophysical studies to characterize their IDRs. We secondly summarize the biophysics governing IDR properties and the state-of-the-art in computational tools toward MAPID identification and characterization of their conformation ensembles. We conclude with an overview of future computational approaches toward broadening the understanding of intrinsic disorder in the cardiac sarcomere.

Prediction of Kv11.1 potassium channel PAS-domain variants trafficking via machine learning.

Congenital long QT syndrome (LQTS) is characterized by a prolonged QT-interval on an electrocardiogram (ECG). An abnormal prolongation in the QT-interval increases the risk for fatal arrhythmias. Genetic variants in several different cardiac ion channel genes, including KCNH2, are known to cause LQTS. Here, we evaluated whether structure-based molecular dynamics (MD) simulations and machine learning (ML) could improve the identification of missense variants in LQTS-linked genes. To do this, we investigated KCNH2 missense variants in the Kv11.1 channel protein shown to have wild type (WT) like or class II (trafficking-deficient) phenotypes in vitro. We focused on KCNH2 missense variants that disrupt normal Kv11.1 channel protein trafficking, as it is the most common phenotype for LQTS-associated variants. Specifically, we used computational techniques to correlate structural and dynamic changes in the Kv11.1 channel protein PAS domain (PASD) with Kv11.1 channel protein trafficking phenotypes. These simulations unveiled several molecular features, including the numbers of hydrating waters and hydrogen bonding pairs, as well as folding free energy scores, that are predictive of trafficking. We then used statistical and machine learning (ML) (Decision tree (DT), Random forest (RF), and Support vector machine (SVM)) techniques to classify variants using these simulation-derived features. Together with bioinformatics data, such as sequence conservation and folding energies, we were able to predict with reasonable accuracy (≈75%) which KCNH2 variants do not traffic normally. We conclude that structure-based simulations of KCNH2 variants localized to the Kv11.1 channel PASD led to an improvement in classification accuracy. Therefore, this approach should be considered to complement the classification of variant of unknown significance (VUS) in the Kv11.1 channel PASD.

Inhibitory and stimulatory micropeptides preferentially bind to different conformations of the cardiac calcium pump.

The ATP-dependent ion pump sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) sequesters Ca2+ in the endoplasmic reticulum to establish a reservoir for cell signaling. Because of its central importance in physiology, the activity of this transporter is tightly controlled via direct interactions with tissue-specific regulatory micropeptides that tune SERCA function to match changing physiological conditions. In the heart, the micropeptide phospholamban (PLB) inhibits SERCA, while dwarf open reading frame (DWORF) stimulates SERCA. These competing interactions determine cardiac performance by modulating the amplitude of Ca2+ signals that drive the contraction/relaxation cycle. We hypothesized that the functions of these peptides may relate to their reciprocal preferences for SERCA binding; SERCA binds PLB more avidly at low cytoplasmic [Ca2+] but binds DWORF better when [Ca2+] is high. In the present study, we demonstrated this opposing Ca2+ sensitivity is due to preferential binding of DWORF and PLB to different intermediate states that SERCA samples during the Ca2+ transport cycle. We show PLB binds best to the SERCA E1-ATP state, which prevails at low [Ca2+]. In contrast, DWORF binds most avidly to E1P and E2P states that are more populated when Ca2+ is elevated. Moreover, FRET microscopy revealed dynamic shifts in SERCA-micropeptide binding equilibria during cellular Ca2+ elevations. A computational model showed that DWORF exaggerates changes in PLB-SERCA binding during the cardiac cycle. These results suggest a mechanistic basis for inhibitory versus stimulatory micropeptide function, as well as a new role for DWORF as a modulator of dynamic oscillations of PLB-SERCA regulatory interactions.

Fluorescence lifetime imaging microscopy reveals sodium pump dimers in live cells.

The sodium-potassium ATPase (Na/K-ATPase, NKA) establishes ion gradients that facilitate many physiological functions including action potentials and secondary transport processes. NKA comprises a catalytic subunit (alpha) that interacts closely with an essential subunit (beta) and regulatory transmembrane micropeptides called FXYD proteins. In the heart, a key modulatory partner is the FXYD protein phospholemman (PLM, FXYD1), but the stoichiometry of the alpha-beta-PLM regulatory complex is unknown. Here, we used fluorescence lifetime imaging and spectroscopy to investigate the structure, stoichiometry, and affinity of the NKA-regulatory complex. We observed a concentration-dependent binding of the subunits of NKA-PLM regulatory complex, with avid association of the alpha subunit with the essential beta subunit as well as lower affinity alpha-alpha and alpha-PLM interactions. These data provide the first evidence that, in intact live cells, the regulatory complex is composed of two alpha subunits associated with two beta subunits, decorated with two PLM regulatory subunits. Docking and molecular dynamics (MD) simulations generated a structural model of the complex that is consistent with our experimental observations. We propose that alpha-alpha subunit interactions support conformational coupling of the catalytic subunits, which may enhance NKA turnover rate. These observations provide insight into the pathophysiology of heart failure, wherein low NKA expression may be insufficient to support formation of the complete regulatory complex with the stoichiometry (alpha-beta-PLM)2.

Molecular Insights into the MLCK Activation by CaM.

Calmodulin (CaM) is a universal regulatory protein that modulates numerous cellular processes by using calcium (Ca2+) as the signal. In smooth muscle cells (SMC), one major target of CaM is myosin light chain kinase (MLCK), a kinase that phosphorylates the myosin regulatory light chain and thereby regulates cell contraction. In the absence of CaM, MLCK remains inhibited by its autoinhibitory domain (AID). While it is well established that CaM activates MLCK, the molecular interactions between these two proteins remain elusive due to the lack of structural data. In this work, we constructed a molecular model of mammalian CaM (mCaM) in complex with MLCK leveraging AlphaFold, published biochemical data, and protein-protein docking. The model, along with a strategic set of CaM mutants including a inhibitory variant soybean CaM isoform 4 (sCaM-4), was subject to molecular dynamics (MD) simulations. Using principal component analysis (PCA), we mapped out the transition path for the removal of the AID from the MLCK kinase domain to provide molecular basis of MLCK activation. Additionally, we established MLCK conformations that correspond to the active and inactive states of the kinase. We showed that mCaM and sCaM-4 cause MLCK to undergo the transition to the active and inactive states, respectively. Using two structural metrics, we computed the probabilities of MLCK activation by different CaM variants, which were in good agreement with the experimental data. Distributions along these metrics revealed that different inhibitory CaM variants impair MLCK activation through unique mechanisms. We finally identified molecular contacts that contribute to the MLCK activation by CaM. Overall, we report a de novo molecular model of CaM-MLCK that provides insights into the molecular mechanism of MLCK activation by CaM. The mechanism requires effective removal of the AID while preserving an active configuration of the kinase domain. This mechanism may be shared by other MLCK isoforms and potentially other structurally similar kinases with CaM-mediated regulatory domains.

Disruption of Z-Disc Function Promotes Mechanical Dysfunction in Human Myocardium: Evidence for a Dual Myofilament Modulatory Role by Alpha-Actinin 2.

The ACTN2 gene encodes α-actinin 2, located in the Z-disc of the sarcomeres in striated muscle. In this study, we sought to investigate the effects of an ACTN2 missense variant of unknown significance (p.A868T) on cardiac muscle structure and function. Left ventricular free wall samples were obtained at the time of cardiac transplantation from a heart failure patient with the ACTN2 A868T heterozygous variant. This variant is in the EF 3-4 domain known to interact with titin and α-actinin. At the ultrastructural level, ACTN2 A868T cardiac samples presented small structural changes in cardiomyocytes when compared to healthy donor samples. However, contractile mechanics of permeabilized ACTN2 A868T variant cardiac tissue displayed higher myofilament Ca2+ sensitivity of isometric force, reduced sinusoidal stiffness, and faster rates of tension redevelopment at all Ca2+ levels. Small-angle X-ray diffraction indicated increased separation between thick and thin filaments, possibly contributing to changes in muscle kinetics. Molecular dynamics simulations indicated that while the mutation does not significantly impact the structure of α-actinin on its own, it likely alters the conformation associated with titin binding. Our results can be explained by two Z-disc mediated communication pathways: one pathway that involves α-actinin's interaction with actin, affecting thin filament regulation, and the other pathway that involves α-actinin's interaction with titin, affecting thick filament activation. This work establishes the role of α-actinin 2 in modulating cross-bridge kinetics and force development in the human myocardium as well as how it can be involved in the development of cardiac disease.

Modulation of P2X4 pore closure by magnesium, potassium, and ATP.

The P2X4 receptor plays a prominent role in cellular responses to extracellular ATP. Through classical all-atom molecular dynamics (MD) simulations totaling 24 µs we have investigated how metal-complexed ATP stabilizes the channel's open state and prevents its closing. We have identified two metal-binding sites, Mg2+ and potassium K+, one at the intersection of the three subunits in the ectodomain (MBS1) and the second one near the ATP-binding site (MBS2), similar to those characterized in Gulf Coast P2X. Our data indicate that when Mg2+ and K+ ions are complexed with ATP, the channel is locked into an open state. Interestingly, irrespective of the number of bound ATP molecules, Mg2+ ions bound to the MBS2 impeded the collapse of the open-state protein to a closed state by stabilizing the ATP-protein interactions. However, when Mg2+ in the MBS2 was replaced with K+ ions, as might be expected when in equilibrium with an extracellular solution, the interactions between the subunits were weakened and the pore collapsed. This collapse was apparent when fewer than two ATPs were bound to MBS2 in the presence of K+. Therefore, the different capacities of common cations to stabilize the channel may underlie a mechanism governing P2X4 channel gating in physiological systems. This study therefore provides structural insights into the differential modulation of ATP activation of P2X4 by Mg2+ and K+.

Crowding within synaptic junctions influences the degradation of nucleotides by CD39 and CD73 ectonucleotidases.

Synapsed cells can communicate using exocytosed nucleotides like adenosine triphosphate (ATP). Ectonucleotidases localized to synaptic junctions degrade nucleotides into metabolites like adenosine monophosphate (AMP) or adenosine. Oftentimes nucleotide degradation occurs in a sequential manner, of which ATP degradation by CD39 and CD73 is a representative example. Here, CD39 first converts ATP and adenosine diphosphate (ADP) into AMP, after which AMP is dephosphorylated into adenosine by CD73. Hence, the concerted activity of CD39 and CD73 can help shape cellular responses to extracellular ATP. In a previous study, we demonstrated that coupled CD39 and CD73 activity within synapse-like junctions is strongly controlled by the enzymes' co-localization, their surface charge densities, and the electrostatic potential of the surrounding cell membranes. In this study, we demonstrate that crowders within synaptic junctions, which can include globular proteins like cytokines and membrane-bound proteins, impact coupled CD39 and CD73 ectonucleotidase activity and, in turn, the availability of intrasynapse ATP. Specifically, we developed a spatially explicit, reaction-diffusion model for the coupled conversion of ATP → AMP and AMP → adenosine in a model synaptic junction with crowders that is solved via the finite element method. Our modeling results suggest that the association rate for ATP to CD39 is strongly influenced by the density of intrasynaptic protein crowders, as increasing crowder density generally suppressed ATP association kinetics. Much of this suppression can be rationalized based on a loss of configurational entropy. The surface charges of crowders can further influence the association rate, with the surprising result that favorable crowder-nucleotide electrostatic interactions can yield CD39 association rates that are faster than crowder-free configurations. However, attractive crowder-nucleotide interactions decrease the rate and efficiency of adenosine production, which in turn increases the availability of ATP and AMP within the synapse relative to crowder-free configurations. These findings highlight how CD39 and CD73 ectonucleotidase activity, electrostatics, and crowding within synapses influence the availability of nucleotides for intercellular communication.

Calmodulin's Interdomain Linker Is Optimized for Dynamics Signal Transmission and Calcium Binding.

Linkers are ubiquitous in multidomain proteins. These linkers are integral to protein functions, and accumulating evidence suggests that the linkers' versatile roles are encoded in their sequences. However, a molecular picture of how amino acid differences in the linker influence protein function is still lacking. By using extensive Gaussian-accelerated MD coupled with dynamic network analysis, we reveal the molecular bases underlying the linker's role in Calmodulin (CaM), a highly conserved Ca2+-signaling hub in eukaryotes. Three CaM constructs comprising a wild-type linker, a flexible linker (four glycines at position D78-S81), and a rigid linker (four prolines at position D78-S81) were simulated. We show that the flexible linker resembles the wild type in allowing CaM to sample a large ensemble of conformations while the rigid linker confines the sampling. Our simulations recapture experimental observations that target binding enhances the Ca2+ affinity to CaM's EF-hand sites at the N-domain. However, only the wild-type linker can both correctly capture the Ca2+ binding order and maintain the α-helical structure of the domain. The other two constructs either bind Ca2+ in an incorrect order or exhibit unfolding of an N-domain helix. We demonstrate that the wild-type linker achieves these outcomes by transmitting interdomain dynamics efficiently. This was evidenced by stronger (anti)correlations among the linker residues, decoupling of the hydrogen bonds between A1-A15 and V35-E45, and structuring of the N-domain for Ca2+ binding. This decoupling was not evident for the other two constructs. Lastly, we show that the wild-type linker's optimal transmission stems from its thermodynamically favorable strain and solvation relative to the other two constructs. Our results show how the linker sequence tunes CaM function, suggesting possible mechanisms for changes in linker properties such as mutations or post-translational modifications to modulate protein/substrate binding.

Mutation-Specific Differences in Kv7.1 (KCNQ1) and Kv11.1 (KCNH2) Channel Dysfunction and Long QT Syndrome Phenotypes.

The electrocardiogram (ECG) empowered clinician scientists to measure the electrical activity of the heart noninvasively to identify arrhythmias and heart disease. Shortly after the standardization of the 12-lead ECG for the diagnosis of heart disease, several families with autosomal recessive (Jervell and Lange-Nielsen Syndrome) and dominant (Romano-Ward Syndrome) forms of long QT syndrome (LQTS) were identified. An abnormally long heart rate-corrected QT-interval was established as a biomarker for the risk of sudden cardiac death. Since then, the International LQTS Registry was established; a phenotypic scoring system to identify LQTS patients was developed; the major genes that associate with typical forms of LQTS were identified; and guidelines for the successful management of patients advanced. In this review, we discuss the molecular and cellular mechanisms for LQTS associated with missense variants in KCNQ1 (LQT1) and KCNH2 (LQT2). We move beyond the "benign" to a "pathogenic" binary classification scheme for different KCNQ1 and KCNH2 missense variants and discuss gene- and mutation-specific differences in K+ channel dysfunction, which can predispose people to distinct clinical phenotypes (e.g., concealed, pleiotropic, severe, etc.). We conclude by discussing the emerging computational structural modeling strategies that will distinguish between dysfunctional subtypes of KCNQ1 and KCNH2 variants, with the goal of realizing a layered precision medicine approach focused on individuals.

Non-Canonical Interaction between Calmodulin and Calcineurin Contributes to the Differential Regulation of Plant-Derived Calmodulins on Calcineurin.

Calmodulin (CaM) serves as an important Ca2+ signaling hub that regulates many protein signaling pathways. Recently, it was demonstrated that plant CaM homologues can regulate mammalian targets, often in a manner that opposes the impact of the mammalian CaM (mCaM). However, the molecular basis of how CaM homologue mutations differentially impact target activation is unclear. To understand these mechanisms, we examined two CaM isoforms found in soybean plants that differentially regulate a mammalian target, calcineurin (CaN). These CaM isoforms, sCaM-1 and sCaM-4, share >90 and ∼78% identity with the mCaM, respectively, and activate CaN with comparable or reduced activity relative to mCaM. We used molecular dynamics (MD) simulations and fluorometric assays of CaN-dependent dephosphorylation of MUF-P to probe whether calcium and protein-protein binding interactions are altered by plant CaMs relative to mCaM as a basis for differential CaN regulation. In the presence of CaN, we found that the two sCaMs' Ca2+ binding properties, such as their predicted coordination of Ca2+ and experimentally measured EC50 [Ca2+] values are comparable to mCaM. Furthermore, the binding of CaM to the CaM binding region (CaMBR) in CaN is comparable among the three CaMs, as evidenced by MD-predicted binding energies and experimentally measured EC50 [CaM] values. However, mCaM and sCaM-1 exhibited binding with a secondary region of CaN's regulatory domain that is weakened for sCaM-4. We speculate that this secondary interaction affects the turnover rate (kcat) of CaN based on our modeling of enzyme activity, which is consistent with our experimental data. Together, our data describe how plant-derived CaM variants alter CaN activity through enlisting interactions other than those directly influencing Ca2+ binding and canonical CaMBR binding, which may additionally play a role in the differential regulation of other mammalian targets.

Assessing the Role of Calmodulin's Linker Flexibility in Target Binding.

Calmodulin (CaM) is a highly-expressed Ca2+ binding protein known to bind hundreds of protein targets. Its binding selectivity to many of these targets is partially attributed to the protein's flexible alpha helical linker that connects its N- and C-domains. It is not well established how its linker mediates CaM's binding to regulatory targets yet. Insights into this would be invaluable to understanding its regulation of diverse cellular signaling pathways. Therefore, we utilized Martini coarse-grained (CG) molecular dynamics simulations to probe CaM/target assembly for a model system: CaM binding to the calcineurin (CaN) regulatory domain. The simulations were conducted assuming a 'wild-type' calmodulin with normal flexibility of its linker, as well as a labile, highly-flexible linker variant to emulate structural changes that could be induced, for instance, by post-translational modifications. For the wild-type model, 98% of the 600 simulations across three ionic strengths adopted a bound complex within 2 µs of simulation time; of these, 1.7% sampled the fully-bound state observed in the experimentally-determined crystallographic structure. By calculating the mean-first-passage-time for these simulations, we estimated the association rate to be ka= 8.7 × 108 M-1 s-1, which is similar to the diffusion-limited, experimentally-determined rate of 2.2 × 108 M-1 s-1. Furthermore, our simulations recapitulated its well-known inverse relationship between the association rate and the solution ionic strength. In contrast, although over 97% of the labile linker simulations formed tightly-bound complexes, only 0.3% achieved the fully-bound configuration. This effect appears to stem from a difference in the ensembles of extended and collapsed states which are controlled by the linker flexibility. Therefore, our simulations suggest that variations in the CaM linker's propensity for alpha helical secondary structure can modulate the kinetics of target binding.

A Matched-Filter-Based Algorithm for Subcellular Classification of T-System in Cardiac Tissues.

In mammalian ventricular cardiomyocytes, invaginations of the surface membrane form the transverse tubular system (T-system), which consists of transverse tubules (TTs) that align with sarcomeres and Z-lines as well as longitudinal tubules (LTs) that are present between Z-lines in some species. In many cardiac disease etiologies, the T-system is perturbed, which is believed to promote spatially heterogeneous, dyssynchronous Ca2+ release and inefficient contraction. In general, T-system characterization approaches have been directed primarily at isolated cells and do not detect subcellular T-system heterogeneity. Here, we present MatchedMyo, a matched-filter-based algorithm for subcellular T-system characterization in isolated cardiomyocytes and millimeter-scale myocardial sections. The algorithm utilizes "filters" representative of TTs, LTs, and T-system absence. Application of the algorithm to cardiomyocytes isolated from rat disease models of myocardial infarction (MI), dilated cardiomyopathy induced via aortic banding, and sham surgery confirmed and quantified heterogeneous T-system structure and remodeling. Cardiomyocytes from post-MI hearts exhibited increasing T-system disarray as proximity to the infarct increased. We found significant (p < 0.05, Welch's t-test) increases in LT density within cardiomyocytes proximal to the infarct (12 ± 3%, data reported as mean ± SD, n = 3) versus sham (4 ± 2%, n = 5), but not distal to the infarct (7 ± 1%, n = 3). The algorithm also detected decreases in TTs within 5° of the myocyte minor axis for isolated aortic banding (36 ± 9%, n = 3) and MI cardiomyocytes located intermediate (37 ± 4%, n = 3) and proximal (34 ± 4%, n = 3) to the infarct versus sham (57 ± 12%, n = 5). Application of bootstrapping to rabbit MI tissue revealed distal sections comprised 18.9 ± 1.0% TTs, whereas proximal sections comprised 10.1 ± 0.8% TTs (p < 0.05), a 46.6% decrease. The matched-filter approach therefore provides a robust and scalable technique for T-system characterization from isolated cells through millimeter-scale myocardial sections.

Simulation of P2X-mediated calcium signalling in microglia.

KEY POINTS: A computational model of P2X channel activation in microglia was developed that includes downfield Ca2+ -dependent signalling pathways. This model provides quantitative insights into how diverse signalling pathways in microglia converge to control microglial function. ABSTRACT: Microglia function is orchestrated through highly coupled signalling pathways that depend on calcium (Ca2+ ). In response to extracellular ATP, transient increases in intracellular Ca2+ driven through the activation of purinergic receptors, P2X and P2Y, are sufficient to promote cytokine synthesis. Although the steps comprising the pathways bridging purinergic receptor activation with transcriptional responses have been probed in great detail, a quantitative model for how these steps collectively control cytokine production has not been established. Here we developed a minimal computational model that quantitatively links extracellular stimulation of two prominent ionotropic purinergic receptors, P2X4 and P2X7, with the graded production of a gene product, namely the tumour necrosis factor α (TNFα) cytokine. In addition to Ca2+ handling mechanisms common to eukaryotic cells, our model includes microglia-specific processes including ATP-dependent P2X4 and P2X7 activation, activation of nuclear factor of activated T-cells (NFAT) transcription factors, and TNFα production. Parameters for this model were optimized to reproduce published data for these processes, where available. With this model, we determined the propensity for TNFα production in microglia, subject to a wide range of ATP exposure amplitudes, frequencies and durations that the cells could encounter in vivo. Furthermore, we have investigated the extent to which modulation of the signal transduction pathways influence TNFα production. Our results suggest that pulsatile stimulation of P2X4 via micromolar ATP may be sufficient to promote TNFα production, whereas high-amplitude ATP exposure is necessary for production via P2X7. Furthermore, under conditions that increase P2X4 expression, for instance, following activation by pathogen-associated molecular factors, P2X4-associated TNFα production is greatly enhanced. Given that Ca2+ homeostasis in microglia is profoundly important to its function, this computational model provides a quantitative framework to explore hypotheses pertaining to microglial physiology.

3D dSTORM imaging reveals novel detail of ryanodine receptor localization in rat cardiac myocytes.

KEY POINTS: Using 3D direct stochastic optical reconstruction microscopy (dSTORM), we developed novel approaches to quantitatively describe the nanoscale, 3D organization of ryanodine receptors (RyRs) in cardiomyocytes. Complex arrangements of RyR clusters were observed in 3D space, both at the cell surface and within the cell interior, with allocation to dyadic and non-dyadic pools. 3D imaging importantly allowed discernment of clusters overlapping in the z-axis, for which detection was obscured by conventional 2D imaging techniques. Thus, RyR clusters were found to be significantly smaller than previous 2D estimates. Ca2+ release units (CRUs), i.e. functional groupings of neighbouring RyR clusters, were similarly observed to be smaller than earlier reports. Internal CRUs contained more RyRs in more clusters than CRUs on the cell surface, and yielded longer duration Ca2+ sparks. ABSTRACT: Cardiomyocyte contraction is dependent on Ca2+ release from ryanodine receptors (RyRs). However, the precise localization of RyRs remains unknown, due to shortcomings of imaging techniques which are diffraction limited or restricted to 2D. We aimed to determine the 3D nanoscale organization of RyRs in rat cardiomyocytes by employing direct stochastic optical reconstruction microscopy (dSTORM) with phase ramp technology. Initial observations at the cell surface showed an undulating organization of RyR clusters, resulting in their frequent overlap in the z-axis and obscured detection by 2D techniques. Non-overlapping clusters were imaged to create a calibration curve for estimating RyR number based on recorded fluorescence blinks. Employing this method at the cell surface and interior revealed smaller RyR clusters than 2D estimates, as erroneous merging of axially aligned RyRs was circumvented. Functional groupings of RyR clusters (Ca2+ release units, CRUs), contained an average of 18 and 23 RyRs at the surface and interior, respectively, although half of all CRUs contained only a single 'rogue' RyR. Internal CRUs were more tightly packed along z-lines than surface CRUs, contained larger and more numerous RyR clusters, and constituted ∼75% of the roughly 1 million RyRs present in an average cardiomyocyte. This complex internal 3D geometry was underscored by correlative imaging of RyRs and t-tubules, which enabled quantification of dyadic and non-dyadic RyR populations. Mirroring differences in CRU size and complexity, Ca2+ sparks originating from internal CRUs were of longer duration than those at the surface. These data provide novel, nanoscale insight into RyR organization and function across cardiomyocytes.

Molecular Basis of S100A1 Activation at Saturating and Subsaturating Calcium Concentrations.

The S100A1 protein mediates a wide variety of physiological processes through its binding of calcium (Ca(2+)) and endogenous target proteins. S100A1 presents two Ca(2+)-binding domains: a high-affinity "canonical" EF (cEF) hand and a low-affinity "pseudo" EF (pEF) hand. Accumulating evidence suggests that both Ca(2+)-binding sites must be saturated to stabilize an open state conducive to peptide recognition, yet the pEF hand's low affinity limits Ca(2+) binding at normal physiological concentrations. To understand the molecular basis of Ca(2+) binding and open-state stabilization, we performed 100 ns molecular dynamics simulations of S100A1 in the apo/holo (Ca(2+)-free/bound) states and a half-saturated state, for which only the cEF sites are Ca(2+)-bound. Our simulations indicate that the pattern of oxygen coordination about Ca(2+) in the cEF relative to the pEF site contributes to the former's higher affinity, whereas Ca(2+) binding strongly reshapes the protein's conformational dynamics by disrupting ß-sheet coupling between EF hands. Moreover, modeling of the half-saturated configuration suggests that the open state is unstable and reverts toward a closed state in the absence of the pEF Ca(2+) ion. These findings indicate that Ca(2+) binding at the cEF site alone is insufficient to stabilize opening; thus, posttranslational modification of the protein may be required for target peptide binding at subsaturating intracellular Ca(2+) levels.

Enzyme localization, crowding, and buffers collectively modulate diffusion-influenced signal transduction: Insights from continuum diffusion modeling.

Biochemical reaction networks consisting of coupled enzymes connect substrate signaling events with biological function. Substrates involved in these reactions can be strongly influenced by diffusion "barriers" arising from impenetrable cellular structures and macromolecules, as well as interactions with biomolecules, especially within crowded environments. For diffusion-influenced reactions, the spatial organization of diffusion barriers arising from intracellular structures, non-specific crowders, and specific-binders (buffers) strongly controls the temporal and spatial reaction kinetics. In this study, we use two prototypical biochemical reactions, a Goodwin oscillator, and a reaction with a periodic source/sink term to examine how a diffusion barrier that partitions substrates controls reaction behavior. Namely, we examine how conditions representative of a densely packed cytosol, including reduced accessible volume fraction, non-specific interactions, and buffers, impede diffusion over nanometer length-scales. We find that diffusion barriers can modulate the frequencies and amplitudes of coupled diffusion-influenced reaction networks, as well as give rise to "compartments" of decoupled reactant populations. These effects appear to be intensified in the presence of buffers localized to the diffusion barrier. These findings have strong implications for the role of the cellular environment in tuning the dynamics of signaling pathways.

Electrostatic channeling in P. falciparum DHFR-TS: Brownian dynamics and Smoluchowski modeling.

We perform Brownian dynamics simulations and Smoluchowski continuum modeling of the bifunctional Plasmodium falciparum dihydrofolate reductase-thymidylate synthase (P. falciparum DHFR-TS) with the objective of understanding the electrostatic channeling of dihydrofolate generated at the TS active site to the DHFR active site. The results of Brownian dynamics simulations and Smoluchowski continuum modeling suggest that compared to Leishmania major DHFR-TS, P. falciparum DHFR-TS has a lower but significant electrostatic-mediated channeling efficiency (?15-25%) at physiological pH (7.0) and ionic strength (150 mM). We also find that removing the electric charges from key basic residues located between the DHFR and TS active sites significantly reduces the channeling efficiency of P. falciparum DHFR-TS. Although several protozoan DHFR-TS enzymes are known to have similar tertiary and quaternary structure, subtle differences in structure, active-site geometry, and charge distribution appear to influence both electrostatic-mediated and proximity-based substrate channeling.

A model study of sequential enzyme reactions and electrostatic channeling.

We study models of two sequential enzyme-catalyzed reactions as a basic functional building block for coupled biochemical networks. We investigate the influence of enzyme distributions and long-range molecular interactions on reaction kinetics, which have been exploited in biological systems to maximize metabolic efficiency and signaling effects. Specifically, we examine how the maximal rate of product generation in a series of sequential reactions is dependent on the enzyme distribution and the electrostatic composition of its participant enzymes and substrates. We find that close proximity between enzymes does not guarantee optimal reaction rates, as the benefit of decreasing enzyme separation is countered by the volume excluded by adjacent enzymes. We further quantify the extent to which the electrostatic potential increases the efficiency of transferring substrate between enzymes, which supports the existence of electrostatic channeling in nature. Here, a major finding is that the role of attractive electrostatic interactions in confining intermediate substrates in the vicinity of the enzymes can contribute more to net reactive throughput than the directional properties of the electrostatic fields. These findings shed light on the interplay of long-range interactions and enzyme distributions in coupled enzyme-catalyzed reactions, and their influence on signaling in biological systems.

Pathogenic TNNI1 variants disrupt sarcomere contractility resulting in hypo- and hypercontractile muscle disease.

Troponin I (TnI) regulates thin filament activation and muscle contraction. Two isoforms, TnI-fast (TNNI2) and TnI-slow (TNNI1), are predominantly expressed in fast- and slow-twitch myofibers, respectively. TNNI2 variants are a rare cause of arthrogryposis, whereas TNNI1 variants have not been conclusively established to cause skeletal myopathy. We identified recessive loss-of-function TNNI1 variants as well as dominant gain-of-function TNNI1 variants as a cause of muscle disease, each with distinct physiological consequences and disease mechanisms. We identified three families with biallelic TNNI1 variants (F1: p.R14H/c.190-9G>A, F2 and F3: homozygous p.R14C), resulting in loss of function, manifesting with early-onset progressive muscle weakness and rod formation on histology. We also identified two families with a dominantly acting heterozygous TNNI1 variant (F4: p.R174Q and F5: p.K176del), resulting in gain of function, manifesting with muscle cramping, myalgias, and rod formation in F5. In zebrafish, TnI proteins with either of the missense variants (p.R14H; p.R174Q) incorporated into thin filaments. Molecular dynamics simulations suggested that the loss-of-function p.R14H variant decouples TnI from TnC, which was supported by functional studies showing a reduced force response of sarcomeres to submaximal [Ca2+] in patient myofibers. This contractile deficit could be reversed by a slow skeletal muscle troponin activator. In contrast, patient myofibers with the gain-of-function p.R174Q variant showed an increased force to submaximal [Ca2+], which was reversed by the small-molecule drug mavacamten. Our findings demonstrated that TNNI1 variants can cause muscle disease with variant-specific pathomechanisms, manifesting as either a hypo- or a hypercontractile phenotype, suggesting rational therapeutic strategies for each mechanism.