In silico assessment of the conduction mechanism of the Ryanodine Receptor 1 reveals previously unknown exit pathways

Eu3+ detects two functionally distinct luminal Ca2+ binding sites in ryanodine receptors

Ryanodine receptors (RyRs) are Ca2+ release channels, gated by Ca2+ in the cytosol and the sarcoplasmic reticulum lumen. Their regulation is impaired in certain cardiac and muscle diseases. Although a lot of data is available on the luminal Ca2+ regulation of RyR, its interpretation is complicated by the possibility that the divalent ions used to probe the luminal binding sites may contaminate the cytoplasmic sites by crossing the channel pore. In this study, we used Eu3+, an impermeable agonist of Ca2+ binding sites, as a probe to avoid this complication and to gain more specific information about the function of the luminal Ca2+ sensor. Single-channel currents were measured from skeletal muscle and cardiac RyRs (RyR1 and RyR2) using the lipid bilayer technique. We show that RyR2 is activated by the luminal addition of Ca2+, whereas RyR1 is inhibited. These results were qualitatively reproducible using Eu3+. The luminal regulation of RyR1 carrying a mutation associated with malignant hyperthermia was not different from that of the wild-type. RyR1 inhibition by Eu3+ was extremely voltage dependent, whereas RyR2 activation did not depend on the membrane potential. These results suggest that the RyR1 inhibition site is in the membrane's electric field (channel pore), whereas the RyR2 activation site is outside. Using in silico analysis and previous results, we predicted putative Ca2+ binding site sequences. We propose that RyR2 bears an activation site, which is missing in RyR1, but both isoforms share the same inhibitory Ca2+ binding site near the channel gate.

The significance of bioelectricity on all levels of organization of an organism. Part 1: From the subcellular level to cells

Bioelectricity plays an essential role in the structural and functional organization of biological organisms. In this first article of our three-part series, we summarize the importance of bioelectricity for the basic structural level of biological organization, i.e. from the subcellular level (charges, ion channels, molecules and cell organelles) to cells.

Molecular Dynamics Simulations of the Cardiac Ryanodine Receptor Type 2 (RyR2) Gating Mechanism

Mutations in the cardiac ryanodine receptor type 2 (RyR2) have been linked to fatal cardiac arrhythmias such as catecholaminergic polymorphic ventricular tachycardia (CPVT). While many CPVT mutations are associated with an increase in Ca²⁺ leak from the sarcoplasmic reticulum, the mechanistic details of RyR2 channel gating are not well understood, and this poses a barrier in the development of new pharmacological treatments. To address this, we explore the gating mechanism of the RyR2 using molecular dynamics (MD) simulations. We test the effect of changing the conformation of certain structural elements by constructing chimera RyR2 structures that are derived from the currently available closed and open cryo-electron microscopy (cryo-EM) structures, and we then use MD simulations to relax the system. Our key finding is that the position of the S4-S5 linker (S4S5L) on a single subunit can determine whether the channel as a whole is open or closed. Our analysis reveals that the position of the S4S5L is regulated by interactions with the U-motif on the same subunit and with the S6 helix on an adjacent subunit. We find that, in general, channel gating is crucially dependent on high percent occupancy interactions between adjacent subunits. We compare our interaction analysis to 49 CPVT1 mutations in the literature and find that 73% appear near a high percent occupancy interaction between adjacent subunits. This suggests that disruption of cooperative, high percent occupancy interactions between adjacent subunits is a primary cause of channel leak and CPVT in mutant RyR2 channels.

Building Biological Relevance Into Integrative Modelling of Macromolecular Assemblies

Recent advances in structural biophysics and integrative modelling methods now allow us to decipher the structures of large macromolecular assemblies. Understanding the dynamics and mechanisms involved in their biological function requires rigorous integration of all available data. We have developed a complete modelling pipeline that includes analyses to extract biologically significant information by consistently combining automated and interactive human-guided steps. We illustrate this idea with two examples. First, we describe the ryanodine receptor, an ion channel that controls ion flux across the cell membrane through transitions between open and closed states. The conformational changes associated with the transitions are small compared to the considerable system size of the receptor; it is challenging to consistently track these states with the available cryo-EM structures. The second example involves homologous recombination, in which long filaments of a recombinase protein and DNA catalyse the exchange of homologous DNA strands to reliably repair DNA double-strand breaks. The nucleoprotein filament reaction intermediates in this process are short-lived and heterogeneous, making their structures particularly elusive. The pipeline we describe, which incorporates experimental and theoretical knowledge combined with state-of-the-art interactive and immersive modelling tools, can help overcome these challenges. In both examples, we point to new insights into biological processes that arise from such interdisciplinary approaches.

Blocking effect of ferritin on the ryanodine receptor-isoform 2

Iron, an essential element for most living organism, participates in a wide variety of physiological processes. Disturbance in iron homeostasis has been associated with numerous pathologies, particularly in the heart and brain, which are the most susceptible organs. Under iron-overload conditions, the generation of reactive oxygen species leads to impairment in Ca2+ signaling, fundamentally implicated in cardiac and neuronal physiology. Since iron excess is accompanied by increased expression of iron-storage protein, ferritin, we examined whether ferritin has an effect on the ryanodine receptor - isoform 2 (RYR2), which is one of the major components of Ca2+ signaling. Using the method of planar lipid membranes, we show that ferritin induced an abrupt, permanent blockage of the RYR2 channel. The ferritin effect was strongly voltage dependent and competitively antagonized by cytosolic TEA+, an impermeant RYR2 blocker. Our results collectively indicate that monomeric ferritin highly likely blocks the RYR2 channel by a direct electrostatic interaction within the wider region of the channel permeation pathway.

Electrical polarity-dependent gating and a unique subconductance of RyR2 induced by S-adenosyl methionine via the ATP binding site

S-Adenosyl-l-methionine (SAM) was used to probe the functional effects exerted via the RyR2 adenine nucleotide binding site. Single channel experiments revealed that SAM applied to the cytoplasmic face of RyR2 had complex voltage dependent effects on channel gating and conductance. At positive transmembrane holding potentials, SAM caused a striking reduction in channel openings and a reduced channel conductance. In contrast, at negative potentials SAM promoted a clearly resolved subconductance state. At membrane potentials between -75 and -25 mV the open probability of the subconductance state was independent of voltage. ATP, but not the non-adenosine based RyR activator 4-chloro-m-cresol interfered with the effects of SAM at both negative and positive potentials. This suggests that ATP and SAM interact with a common binding site. Molecular docking showed SAM bound to the adenine nucleotide-binding site and formed a hydrogen bond to Glu4886 in the C-terminal end of the S6 alpha helix. In this configuration SAM may alter the conformation of the RyR2 ion conduction pathway. This work provides novel insight into potential functional outcomes of ligand binding to the RyR adenine nucleotide binding site.

Atomistic Details of Charge/Space Competition in the Ca2+ Selectivity of Ryanodine Receptors

Ryanodine receptors (RyRs) are ion channels responsible for the fast release of Ca²⁺ from the sarco/endoplasmic reticulum to the cytosol and show a selectivity of Ca²⁺ over monovalent cations. By utilizing a recently developed multisite Ca²⁺ model in molecular dynamic simulations, we show that multiple cations accumulate in the upper selectivity filter of RyRs, and the small size and high valence of Ca²⁺ make it preferable to K⁺ in competition for space in this confined region of negative electrostatic potential. The presence of Ca²⁺ in the upper selectivity filter significantly increases the energy barrier of K⁺ permeation, while the presence of K⁺ has little impact on the Ca²⁺ permeation. Our results provide the atomistic details of the charge/space competition mechanism for the ion selectivity of RyRs, which ensures the robustness of their Ca²⁺ release function. The mechanism could be utilized in protein- and nanoengineering for valence selectivity of ion species.

Interaction of a dirhamnolipid biosurfactant with sarcoplasmic reticulum calcium ATPase (SERCA1a)

The interaction of a dirhamnolipid biosurfactant secreted by Pseudomonas aeruginosa with calcium ATPase from sarcoplasmic reticulum (SR) was studied by means of different approaches, such as enzyme activity, fluorescence spectroscopy, Fourier transform infrared spectroscopy (FTIR), differential scanning calorimetry (DSC), and molecular docking simulations. The ATP hydrolysis activity was fully inhibited by incubation with dirhamnolipid (diRL) up to 0.1 mM concentration, corresponding to a surfactant concentration below membrane solubilization threshold. Surfactant-protein interaction induced conformational changes in the protein observed by an increase in the accessibility of tryptophan residues to the aqueous phase and by changes in the secondary structure of the protein as seen by fluorescence and FTIR spectroscopy. As a consequence, the protein become more unstable and denatured at lower temperatures, as seen by enzyme activity and DSC studies. Finally, these results were explained at molecular level throughout molecular docking simulations. It is concluded that there is a specific dirhamnolipid-protein interaction not related to the surface activity of the surfactant but to the particular physicochemical properties of the biosurfactant molecule.

The Ca2+ permeation mechanism of the ryanodine receptor revealed by a multi-site ion model

Ryanodine receptors (RyR) are ion channels responsible for the release of Ca²⁺ from the sarco/endoplasmic reticulum and play a crucial role in the precise control of Ca²⁺ concentration in the cytosol. The detailed permeation mechanism of Ca²⁺ through RyR is still elusive. By using molecular dynamics simulations with a specially designed Ca²⁺ model, we show that multiple Ca²⁺ ions accumulate in the upper selectivity filter of RyR1, but only one Ca²⁺ can occupy and translocate in the narrow pore at a time, assisted by electrostatic repulsion from the Ca²⁺ within the upper selectivity filter. The Ca²⁺ is nearly fully hydrated with the first solvation shell intact during the whole permeation process. These results suggest a remote knock-on permeation mechanism and one-at-a-time occupation pattern for the hydrated Ca²⁺ within the narrow pore, uncovering the basis underlying the high permeability and low selectivity of the RyR channels.

Molecular model of the ferroportin intracellular gate and implications for the human iron transport cycle and hemochromatosis type 4A

Ferroportin 1 (FPN1) is a major facilitator superfamily transporter that is essential for proper maintenance of human iron homeostasis at the systemic and cellular level. FPN1 dysfunction leads to the progressive accumulation of iron in reticuloendothelial cells, causing hemochromatosis type 4A (or ferroportin disease), an autosomal dominant disorder that displays large phenotypic heterogeneity. Although crystal structures have unveiled the outward- and inward-facing conformations of the bacterial homolog Bdellovibrio bacteriovorus Fpn (or Bd2019) and calcium has recently been identified as an essential cofactor, our molecular understanding of the iron transport mechanism remains incomplete. Here, we used a combination of molecular modeling, molecular dynamics simulations, and Ala site-directed mutagenesis, followed by complementary in vitro functional analyses, to explore the structural architecture of the human FPN1 intracellular gate. We reveal an interdomain network that involves 5 key amino acids and is likely very important for stability of the iron exporter facing the extracellular milieu. We also identify inter- and intradomain interactions that rely on the 2 Asp84 and Asn174 critical residues and do not exist in the bacterial homolog. These interactions are thought to play an important role in the modulation of conformational changes during the transport cycle. We interpret these results in the context of hemochromatosis type 4A, reinforcing the idea that different categories of loss-of-function mutations exist. Our findings provide an unprecedented view of the human FPN1 outward-facing structure and the particular function of the so-called "gating residues" in the mechanism of iron export.-Guellec, J., Elbahnsi, A., Le Tertre, M., Uguen, K., Gourlaouen, I., Férec, C., Ka, C., Callebaut, I., Le Gac, G. Molecular model of the ferroportin intracellular gate and implications for the human iron transport cycle and hemochromatosis type 4A.

Correlation of phenotype with genotype and protein structure in RYR1-related disorders

Variants in the skeletal muscle ryanodine receptor 1 gene (RYR1) result in a spectrum of RYR1-related disorders. Presentation during infancy is typical and ranges from delayed motor milestones and proximal muscle weakness to severe respiratory impairment and ophthalmoplegia. We aimed to elucidate correlations between genotype, protein structure and clinical phenotype in this rare disease population. Genetic and clinical data from 47 affected individuals were analyzed and variants mapped to the cryo-EM RyR1 structure. Comparisons of clinical severity, motor and respiratory function and symptomatology were made according to the mode of inheritance and affected RyR1 structural domain(s). Overall, 49 RYR1 variants were identified in 47 cases (dominant/de novo, n = 35; recessive, n = 12). Three variants were previously unreported. In recessive cases, facial weakness, neonatal hypotonia, ophthalmoplegia/paresis, ptosis, and scapular winging were more frequently observed than in dominant/de novo cases (all, p < 0.05). Both dominant/de novo and recessive cases exhibited core myopathy histopathology. Clinically severe cases were typically recessive or had variants localized to the RyR1 cytosolic shell domain. Motor deficits were most apparent in the MFM-32 standing and transfers dimension, [median (IQR) 85.4 (18.8)% of maximum score] and recessive cases exhibited significantly greater overall motor function impairment compared to dominant/de novo cases [79.7 (18.8)% vs. 87.5 (17.7)% of maximum score, p = 0.03]. Variant mapping revealed patterns of clinical severity across RyR1 domains, including a structural plane of interest within the RyR1 cytosolic shell, in which 84% of variants affected the bridging solenoid. We have corroborated genotype-phenotype correlations and identified RyR1 regions that may be especially sensitive to structural modification.

Molecular Mechanism of Conductance Enhancement in Narrow Cation-Selective Membrane Channels

Membrane proteins known as ryanodine receptors (RyRs) display large conductance of ∼1 nS and nearly ideal charge selectivity. Both properties are inversely correlated in other large-conductance but nonselective biological nanopores (i.e., α-hemolysin) used as industrial biosensors. Although recent cryo-electron microscopy structures of RyR2 show similarities to K⁺- and Na⁺-selective channels, it remains unclear whether similar ion conduction mechanisms occur in RyR2. Here, we combine microseconds of all-atom molecular dynamics (MD) simulations with mutagenesis and electrophysiology experiments to investigate large K⁺ conductance and charge selectivity (cation vs anion) in an open-state structure of RyR2. Our results show that a water-mediated knock-on mechanism enhances the cation permeation. The polar Q4863 ring may function as a confinement zone amplifying charge selectivity, while the cytoplasmic vestibule can contribute to the efficiency of the cation attraction. We also provide direct evidence that the rings of acidic residues at the channel vestibules are critical for both conductance and charge discrimination in RyRs.