Search

The ryanodine receptor type 2 (RyR2) is composed of four subunits that control calcium (Ca) release in cardiac cells. RyR2 serves primarily as a Ca sensor and can respond to rapid sub-millisecond pulses of Ca while remaining shut at resting concentrations. However, it is not known how the four subunits interact for the RyR2 to function as an effective Ca sensor. To address this question, and to understand the role of subunit cooperativity in Ca-mediated signal transduction, we have developed a computational model of the RyR2 composed of four interacting subunits. We first analyze the statistical properties of a single RyR2 tetramer, where each subunit can exist in a closed or open conformation. Our findings indicate that the number of subunits in the open state is a crucial parameter that dictates RyR2 kinetics. We find that three or four open subunits are required for the RyR2 to harness cooperative interactions to respond to sub-millisecond changes in Ca, while at the same time remaining shut at the resting Ca levels in the cardiac cell. If the required number of open subunits is lowered to one or two, the RyR2 cannot serve as a robust Ca sensor, as the large cooperativity required to stabilize the closed state prevents channel activation. Using this four-subunit model, we analyze the kinetics of Ca release from a RyR2 cluster. We show that the closure of a cluster of RyR2 channels is highly sensitive to the balance of cooperative interactions between closed and open subunits. Based on this result, we analyze how specific interactions between RyR2 subunits can induce persistent Ca leak from the sarcoplasmic reticulum (SR), which is believed to be arrhythmogenic. Thus, these results provide a framework to analyze how a pharmacologic or genetic modification of RyR2 subunit cooperativity can induce abnormal Ca cycling that can potentially lead to life-threatening arrhythmias., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2022 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The ionic atmosphere around nucleic acids remains only partially understood at atomic-level detail. Ion counting (IC) experiments provide a quantitative measure of the ionic atmosphere around nucleic acids and, as such, are a natural route for testing quantitative theoretical approaches. In this article, we replicate IC experiments involving duplex DNA in NaCl(aq) using molecular dynamics (MD) simulation, the three-dimensional reference interaction site model (3D-RISM), and nonlinear Poisson-Boltzmann (NLPB) calculations and test against recent buffer-equilibration atomic emission spectroscopy measurements. Further, we outline the statistical mechanical basis for interpreting IC experiments and clarify the use of specific concentration scales. Near physiological concentrations, MD simulation and 3D-RISM estimates are close to experimental results, but at higher concentrations (>0.7 M), both methods underestimate the number of condensed cations and overestimate the number of excluded anions. The effect of DNA charge on ion and water atmosphere extends 20-25 Å from its surface, yielding layered density profiles. Overall, ion distributions from 3D-RISMs are relatively close to those from corresponding MD simulations, but with less Na(+) binding in grooves and tighter binding to phosphates. NLPB calculations, on the other hand, systematically underestimate the number of condensed cations at almost all concentrations and yield nearly structureless ion distributions that are qualitatively distinct from those generated by both MD simulation and 3D-RISM. These results suggest that MD simulation and 3D-RISM may be further developed to provide quantitative insight into the characterization of the ion atmosphere around nucleic acids and their effect on structure and stability., (Copyright © 2014 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Several isotypes of the structural protein tubulin have been characterized. Their expression offers a plausible explanation for differences regarding microtubule function. Although sequence variation between tubulin isotypes occurs throughout the entire protein, it is the extreme carboxy-terminal tails (CTTs) that exhibit the greatest concentration of differences. In humans, the CTTs range in length from 9 to 25 residues and because of a considerable number of glutamic acid residues, contain over 1/3 of tubulin's total electrostatic charge. The CTTs are believed to be highly disordered and their precise function has yet to be determined. However, their absence has been shown to result in altered microtubule stability and a reduction in the interaction with several microtubule-associated proteins (MAPs). To characterize the role that CTTs play in microtubule function, we examined the global conformational differences within a set of nine human beta-tubulin isotypes using replica exchange molecular dynamics simulations. Through the analysis of the resulting configuration ensembles, we quantified differences such as the CTTs sequence influence on overall flexibility and average secondary structure. Although only minor variations between each CTT were observed, we suggest that these differences may be significant enough to affect interactions with MAPs, thereby influencing important properties such as microtubule assembly and stability.