Thermodynamics of calmodulin binding to cardiac and skeletal muscle ryanodine receptor ion channels

Kinetics and mapping of Ca-driven calmodulin conformations on skeletal and cardiac muscle ryanodine receptors

Calmodulin transduces [Ca2+] information regulating the rhythmic Ca2+ cycling between the sarcoplasmic reticulum and cytoplasm during contraction and relaxation in cardiac and skeletal muscle. However, the structural dynamics by which calmodulin modulates the sarcoplasmic reticulum Ca2+ release channel, the ryanodine receptor, at physiologically relevant [Ca2+] is unknown. Using fluorescence lifetime FRET, we resolve different structural states of calmodulin and Ca2+-driven shifts in the conformation of calmodulin bound to ryanodine receptor. Skeletal and cardiac ryanodine receptor isoforms show different calmodulin-ryanodine receptor conformations, as well as binding and structural kinetics with 0.2-ms resolution, which reflect different functional roles of calmodulin. These FRET methods provide insight into the physiological calmodulin-ryanodine receptor structural states, revealing additional distinct structural states that complement cryo-EM models that are based on less physiological conditions. This technology will drive future studies on pathological calmodulin-ryanodine receptor interactions and dynamics with other important ryanodine receptor bound modulators.

Isolated Cardiac Ryanodine Receptor Function Varies Between Mammals

Concerted robust opening of cardiac ryanodine receptors' (RyR2) Ca2+ release 1oplasmic reticulum (SR) is fundamental for normal systolic cardiac function. During diastole, infrequent spontaneous RyR2 openings mediate the SR Ca2+ leak that normally constrains SR Ca2+ load. Abnormal large diastolic RyR2-mediated Ca2+ leak events can cause delayed after depolarizations (DADs) and arrhythmias. The RyR2-associated mechanisms underlying these processes are being extensively studied at multiple levels utilizing various model animals. Since there are well-described species-specific differences in cardiac intracellular Ca2+ handing in situ, we tested whether or not single RyR2 function in vitro retains this species specificity. We isolated RyR2-rich heavy SR microsomes from mouse, rat, rabbit, and human ventricular muscle and quantified RyR2 function using identical solutions and methods. The single RyR2 cytosolic Ca2+ sensitivity was similar across these species. However, there were significant species differences in single RyR2 mean open times in both systole and diastole-like solutions. In diastole-like solutions, single rat/mouse RyR2 open probability and frequency of long openings (> 6 ms) were similar, but these values were significantly greater than those of either single rabbit or human RyR2s. We propose these in vitro single RyR2 functional differences across species stem from the species-specific RyR2 regulatory environment present in the source tissue. Our results show the single rabbit RyR2 functional attributes, particularly in diastole-like conditions, replicate those of single human RyR2 best among the species tested.

Exploration of Tilmicosin Cardiotoxicity in Rats and the Protecting Role of the Rhodiola rosea Extract: Potential Roles of Cytokines, Antioxidant, Apoptotic, and Anti-Fibrotic Pathways

Tilmicosin (TIL) is a common macrolide antibiotic in veterinary medicine. High doses of TIL can have adverse cardiovascular effects. This study examined the effects of Rhodiola rosea (RHO) that have anti-inflammatory, antioxidant, and anti-fibrotic effects on tilmicosin (TIL)-induced cardiac injury targeting anti-inflammatory, antioxidant, apoptotic, and anti-apoptotic signaling pathways with anti-fibrotic outcomes. Thirty-six male Wistar albino rats were randomly divided into groups of six rats each. Rats received saline as a negative control, CARV 1 mL orally (10 mg/kg BW), and RHO 1 mL orally at 400 mg/kg BW daily for 12 consecutive days. The TIL group once received a single subcutaneous injection (SC) dose of TIL (75 mg/kg BW) on the sixth day of the experiment to induce cardiac damage. The standard group (CARV + TIL) received CARV daily for 12 consecutive days with a single TIL SC injection 1 h after CARV administration only on the sixth day of study and continued for another six successive days on CARV. The protective group (RHO + TIL) received RHO daily for the same period as in CARV + TIL-treated rats and with the dosage mentioned before. Serum was extracted at the time of the rat's scarification at 13 days of study and examined for biochemical assessments in serum lactate dehydrogenase (LDH), cardiac troponin I (cTI), and creatine phosphokinase (CK-MB). Protein carbonyl (PC) contents, malondialdehyde (MDA), and total antioxidant capacity (TAC) in cardiac homogenate were used to measure these oxidative stress markers. Quantitative RT-PCR was used to express interferon-gamma (INF-γ), cyclooxygenase-2 (COX-2), OGG1, BAX, caspase-3, B-cell lymphoma-2 (Bcl-2), and superoxide dismutase (SOD) genes in cardiac tissues, which are correlated with inflammation, antioxidants, and apoptosis. Alpha-smooth muscle actin (α-SMA), calmodulin (CaMKII), and other genes associated with Ca2+ hemostasis and fibrosis were examined using IHC analysis in cardiac cells (myocardium). TIL administration significantly increased the examined cardiac markers, LDH, cTI, and CK-MB. TIL administration also increased ROS, PC, and MDA while decreasing antioxidant activities (TAC and SOD mRNA) in cardiac tissues. Serum inflammatory cytokines and genes of inflammatory markers, DNA damage (INF-γ, COX-2), and apoptotic genes (caspase-3 and BAX) were upregulated with downregulation of the anti-apoptotic gene Bcl-2 as well as the DNA repair OGG1 in cardiac tissues. Furthermore, CaMKII and α-SMA genes were upregulated at cellular levels using cardiac tissue IHC analysis. On the contrary, pretreatment with RHO and CARV alone significantly decreased the cardiac injury markers induced by TIL, inflammatory and anti-inflammatory cytokines, and tissue oxidative-antioxidant parameters. INF-γ, COX-2, OGG1, BAX, and caspase-3 mRNA were downregulated, as observed by real-time PCR, while SOD and Bcl-2 mRNA were upregulated. Furthermore, the CaMKII and α-SMA genes' immune reactivities were significantly decreased in the RHO-pretreated rats.

Computational Analysis of Binding Interactions between the Ryanodine Receptor Type 2 and Calmodulin

Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to a variety of cardiac arrhythmias, such as catecholaminergic polymorphic ventricular tachycardia (CPVT). RyR2 is regulated by calmodulin (CaM), and mutations that disrupt their interaction can cause aberrant calcium release, leading to an arrhythmia. It was recently shown that increasing the RyR2-CaM binding affinity could rescue a defective CPVT-related RyR2 channel to near wild-type behavior. However, the interactions that determine the binding affinity at the RyR2-CaM binding interface are not well understood. In this study, we identify the key domains and interactions, including several new interactions, involved in the binding of CaM to RyR2. Also, our comparison between the wild-type and V3599K mutant suggests how the RyR2-CaM binding affinity can be increased via a change in the central and N-terminal lobe binding contacts for CaM. This computational approach provides new insights into the effect of a mutation at the RyR2-CaM binding interface, and it may find utility in drug design for the future treatment of cardiac arrhythmias.

Neuron-derived factors negatively modulate ryanodine receptor-mediated calcium release in cultured mouse astrocytes

Changes in intracellular Ca²⁺ concentration ([Ca²⁺]_i) produced by ryanodine receptor (RyR) agonist, caffeine (caf), and ionotropic agonists: N-methyl-d-aspartate (NMDA) receptor (NMDAR) agonist, NMDA and P2X7 receptor (P2X7R) agonist, 3'-O-(4-benzoyl)benzoyl adenosine 5'-triphosphate (BzATP) were measured in cultured mouse cortical astrocytes loaded with the fluorescent calcium indicator Fluo3-AM in a confocal laser scanning microscope. In mouse astrocytes cultured in standard medium (SM), treatment with caf increased [Ca²⁺]_i, with a peak response occurring about 10 min after stimulus application. Peak responses to NMDA or BzATP were observed about <1 min and 4.5 min post stimulus, respectively. Co-treatment with NMDA or BzATP did not alter the peak response to caf in astrocytes cultured in SM, the absence of the effects being most likely due to asynchrony between the response to caf, NMDA and BzATP. Incubation of astrocytes with neuron-condition medium (NCM) for 24 h totally abolished the caf-evoked [Ca²⁺]_i increase. In NCM-treated astrocytes, peak of [Ca²⁺]_i rise evoked by NMDA was delayed to about 3.5 min, and that induced by BzATP occurred about three minutes earlier than in SM. The results show that neurons secrete factors that negatively modulate RyR-mediated Ca²⁺-induced Ca²⁺ release (CICR) in astrocytes and alter the time course of Ca²⁺ responses to ionotropic stimuli.

The structural basis of ryanodine receptor ion channel function

Large-conductance Ca²⁺ release channels known as ryanodine receptors (RyRs) mediate the release of Ca²⁺ from an intracellular membrane compartment, the endo/sarcoplasmic reticulum. There are three mammalian RyR isoforms: RyR1 is present in skeletal muscle; RyR2 is in heart muscle; and RyR3 is expressed at low levels in many tissues including brain, smooth muscle, and slow-twitch skeletal muscle. RyRs form large protein complexes comprising four 560-kD RyR subunits, four ∼12-kD FK506-binding proteins, and various accessory proteins including calmodulin, protein kinases, and protein phosphatases. RyRs share ∼70% sequence identity, with the greatest sequence similarity in the C-terminal region that forms the transmembrane, ion-conducting domain comprising ∼500 amino acids. The remaining ∼4,500 amino acids form the large regulatory cytoplasmic "foot" structure. Experimental evidence for Ca²⁺, ATP, phosphorylation, and redox-sensitive sites in the cytoplasmic structure have been described. Exogenous effectors include the two Ca²⁺ releasing agents caffeine and ryanodine. Recent work describing the near atomic structures of mammalian skeletal and cardiac muscle RyRs provides a structural basis for the regulation of the RyRs by their multiple effectors.

High-Throughput Screens to Discover Small-Molecule Modulators of Ryanodine Receptor Calcium Release Channels

Using time-resolved fluorescence resonance energy transfer (FRET), we have developed and validated the first high-throughput screening (HTS) method to discover compounds that modulate an intracellular Ca²⁺ channel, the ryanodine receptor (RyR), for therapeutic applications. Intracellular Ca²⁺ regulation is critical for striated muscle function, and RyR is a central player. At resting [Ca²⁺], an increased propensity of channel opening due to RyR dysregulation is associated with severe cardiac and skeletal myopathies, diabetes, and neurological disorders. This leaky state of the RyR is an attractive target for pharmacological agents to treat such pathologies. Our FRET-based HTS detects RyR binding of accessory proteins calmodulin (CaM) or FKBP12.6. Under conditions that mimic a pathological state, we carried out a screen of the 727-compound NIH Clinical Collection, which yielded six compounds that reproducibly changed FRET by >3 SD. Dose-response of FRET and [³H]ryanodine binding readouts reveal that five hits reproducibly alter RyR1 structure and activity. One compound increased FRET and inhibited RyR1, which was only significant at nM [Ca²⁺], and accentuated without CaM present. These properties characterize a compound that could mitigate RyR1 leak. An excellent Z' factor and the tight correlation between structural and functional readouts validate this first HTS method to identify RyR modulators.

Lobe-specific calmodulin binding to different ryanodine receptor isoforms

Ryanodine receptors (RyRs) are large ion channels that are responsible for the release of Ca(2+) from the sarcoplasmic/endoplasmic reticulum. Calmodulin (CaM) is a Ca(2+) binding protein that can affect the channel open probability at both high and low Ca(2+) concentrations, shifting the Ca(2+) dependencies of channel opening in an isoform-specific manner. Here we analyze the binding of CaM and its individual domains to three different RyR regions using isothermal titration calorimetry. We compared binding to skeletal muscle (RyR1) and cardiac (RyR2) isoforms, under both Ca(2+)-loaded and Ca(2+)-free conditions. CaM can bind all three regions in both isoforms, but the binding modes differ appreciably in two segments. The results highlight a Ca(2+)/CaM and apoCaM binding site in the C-terminal fifth of the channel. This binding site is the target for malignant hyperthermia and central core disease mutations in RyR1, which affect the energetics and mode of CaM binding.

Temperature and transmural region influence functional measurements in unloaded left ventricular cardiomyocytes

Intact cardiomyocytes are increasingly being used to investigate the molecular mechanisms of contraction and to screen new therapeutic compounds. The function of the cardiomyocytes is often measured from the calcium transients and sarcomere length profiles. We studied the role of experimental temperature and transmural region on indices of function in freshly isolated, unloaded cardiomyocytes. Intact cardiomyocytes were isolated from the subendocardium, midmyocardium, and subepicardium of 3-month-old Sprague-Dawley rats. Myocytes from each region were studied at 25°C, 31°C, and 37°C. Cytosolic calcium transients were measured using Fura-2 fluorescence, whereas sarcomere length shortening and relengthening profiles were measured using high-speed video capture. For both the calcium transients and sarcomere length profiles, the time to peak and the time to half relaxation decreased significantly with increasing temperature. Increasing temperature also raised the minimum and maximum calcium levels of all cells. Of note, there was a reduced coefficient of variation (standard deviation divided by the mean) at higher temperatures for calcium fluorescence amplitudes, time to peak calcium, and rates of sarcomeric shortening and relengthening. The amplitudes and minimum of the calcium transients were significantly dependent on transmural region, and several sarcomere length parameters exhibited statistical interactions between temperature and transmural region. Together, these results show that biological variability can be reduced by performing experiments at 37°C rather than at room temperature, and by isolating cells from a specific transmural region. Adopting these procedures will improve the statistical power of subsequent analyses and increase the efficiency of future experiments.

FRET detection of calmodulin binding to the cardiac RyR2 calcium release channel

Calmodulin (CaM) binding to the type 2 ryanodine receptor (RyR2) regulates Ca release from the cardiac sarcoplasmic reticulum (SR). However, the structural basis of CaM regulation of the RyR2 is poorly defined, and the presence of other potential CaM binding partners in cardiac myocytes complicates resolution of CaM's regulatory interactions with RyR2. Here, we show that a fluorescence-resonance-energy-transfer (FRET)-based approach can effectively resolve RyR2 CaM binding, both in isolated SR membrane vesicles and in permeabilized ventricular myocytes. A small FRET donor was targeted to the RyR2 cytoplasmic assembly via fluorescent labeling of the FKBP12.6 subunit. Acceptor fluorophore was attached at discrete positions within either the N- or the C-lobe of CaM. FRET between FKBP12.6 and CaM bound to SR vesicles indicated CaM binding at a single high-affinity site within 60 Å of FKBP12.6. Micromolar Ca increased the apparent affinity of CaM binding and slowed CaM dissociation, but did not significantly affect maximal FRET efficiency at saturating CaM. FRET was strongest when the acceptor was attached at either of two positions within CaM's N-lobe versus sites in CaM's C-lobe, providing CaM orientation information. In permeabilized ventricular myocytes, FKBP12.6 and CaM colocalized to Z-lines, and the efficiency of energy transfer to both the N- and C-lobes of CaM was comparable to that observed in SR vesicle experiments. Results also indicate that both the location and orientation of CaM binding on the RyR2 are very similar to the skeletal muscle RyR1 isoform. Specific binding of CaM to functional RyR2 channels in the cardiac myocyte environment can be monitored using FKBP biosensors and FRET.

Energy variational analysis of ions in water and channels: Field theory for primitive models of complex ionic fluids

Ionic solutions are mixtures of interacting anions and cations. They hardly resemble dilute gases of uncharged noninteracting point particles described in elementary textbooks. Biological and electrochemical solutions have many components that interact strongly as they flow in concentrated environments near electrodes, ion channels, or active sites of enzymes. Interactions in concentrated environments help determine the characteristic properties of electrodes, enzymes, and ion channels. Flows are driven by a combination of electrical and chemical potentials that depend on the charges, concentrations, and sizes of all ions, not just the same type of ion. We use a variational method EnVarA (energy variational analysis) that combines Hamilton's least action and Rayleigh's dissipation principles to create a variational field theory that includes flow, friction, and complex structure with physical boundary conditions. EnVarA optimizes both the action integral functional of classical mechanics and the dissipation functional. These functionals can include entropy and dissipation as well as potential energy. The stationary point of the action is determined with respect to the trajectory of particles. The stationary point of the dissipation is determined with respect to rate functions (such as velocity). Both variations are written in one Eulerian (laboratory) framework. In variational analysis, an "extra layer" of mathematics is used to derive partial differential equations. Energies and dissipations of different components are combined in EnVarA and Euler-Lagrange equations are then derived. These partial differential equations are the unique consequence of the contributions of individual components. The form and parameters of the partial differential equations are determined by algebra without additional physical content or assumptions. The partial differential equations of mixtures automatically combine physical properties of individual (unmixed) components. If a new component is added to the energy or dissipation, the Euler-Lagrange equations change form and interaction terms appear without additional adjustable parameters. EnVarA has previously been used to compute properties of liquid crystals, polymer fluids, and electrorheological fluids containing solid balls and charged oil droplets that fission and fuse. Here we apply EnVarA to the primitive model of electrolytes in which ions are spheres in a frictional dielectric. The resulting Euler-Lagrange equations include electrostatics and diffusion and friction. They are a time dependent generalization of the Poisson-Nernst-Planck equations of semiconductors, electrochemistry, and molecular biophysics. They include the finite diameter of ions. The EnVarA treatment is applied to ions next to a charged wall, where layering is observed. Applied to an ion channel, EnVarA calculates a quick transient pile-up of electric charge, transient and steady flow through the channel, stationary "binding" in the channel, and the eventual accumulation of salts in "unstirred layers" near channels. EnVarA treats electrolytes in a unified way as complex rather than simple fluids. Ad hoc descriptions of interactions and flow have been used in many areas of science to deal with the nonideal properties of electrolytes. It seems likely that the variational treatment can simplify, unify, and perhaps derive and improve those descriptions.