A multiscale model linking ion-channel molecular dynamics and electrostatics to the cardiac action potential.
Ontology highlight
ABSTRACT:
SUBMITTER: Silva JR
PROVIDER: S-EPMC2700153 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Mathematical models for the action potential (AP) generation of the electrically excitable cells including the heart are involved different mechanisms including the voltage-dependent currents with nonlinear time- and voltage-gating properties. From the shape of the AP waveforms to the duration of the refractory periods or heart rhythms are greatly affected by the functions describing the features or the quantities of these ion channels. In this work, a mathematical measure to analyze the regional contributions of voltage-gated channels is defined by dividing the AP into phases, epochs, and intervals of interest. The contribution of each time-dependent current for the newly defined cardiomyocyte model is successfully calculated and it is found that the contribution of dominant ion channels changes substantially not only for each phase but also for different regions of the cardiac AP. Besides, the defined method can also be applied in all Hodgkin-Huxley types of electrically excitable cell models to be able to understand the underlying dynamics better.
| S-EPMC8288042 | biostudies-literature
Project description:Molecular-dynamics simulations were carried out to ascertain which of the potential multimeric forms of the transmembrane peptaibol channel, antiamoebin, is consistent with its measured conductance. Estimates of the conductance obtained through counting ions that cross the channel and by solving the Nernst-Planck equation yield consistent results, indicating that the motion of ions inside the channel can be satisfactorily described as diffusive. The calculated conductance of octameric channels is markedly higher than the conductance measured in single channel recordings, whereas the tetramer appears to be nonconducting. The conductance of the hexamer was estimated to be 115 ± 34 pS and 74 ± 20 pS, at 150 mV and 75 mV, respectively, in satisfactory agreement with the value of 90 pS measured at 75 mV. On this basis, we propose that the antiamoebin channel consists of six monomers. Its pore is large enough to accommodate K? and Cl? with their first solvation shells intact. The free energy barrier encountered by K? is only 2.2 kcal/mol whereas Cl? encounters a substantially higher barrier of nearly 5 kcal/mol. This difference makes the channel selective for cations. Ion crossing events are shown to be uncorrelated and follow Poisson statistics.
| S-EPMC3093558 | biostudies-literature
Project description:Many processes important to chemistry, materials science, and biology cannot be described without considering electronic and nuclear-level dynamics and their coupling to slower, cooperative motions of the system. These inherently multiscale problems require computationally efficient and accurate methods to converge statistical properties. In this paper, a method is presented that uses data directly from condensed phase ab initio simulations to develop reactive molecular dynamics models that do not require predefined empirical functions. Instead, the interactions used in the reactive model are expressed as linear combinations of interpolating functions that are optimized by using a linear least-squares algorithm. One notable benefit of the procedure outlined here is the capability to minimize the number of parameters requiring nonlinear optimization. The method presented can be generally applied to multiscale problems and is demonstrated by generating reactive models for the hydrated excess proton and hydroxide ion based directly on condensed phase ab initio molecular dynamics simulations. The resulting models faithfully reproduce the water-ion structural properties and diffusion constants from the ab initio simulations. Additionally, the free energy profiles for proton transfer, which is sensitive to the structural diffusion of both ions in water, are reproduced. The high fidelity of these models to ab initio simulations will permit accurate modeling of general chemical reactions in condensed phase systems with computational efficiency orders of magnitudes greater than currently possible with ab initio simulation methods, thus facilitating a proper statistical sampling of the coupling to slow, large-scale motions of the system.
| S-EPMC3432097 | biostudies-other
Project description:Many studies have revealed changes in specific protein channels due to physiological causes such as mutation and their effects on action potential duration changes. However, no studies have been conducted to predict the type of protein channel abnormalities that occur through an action potential (AP) shape. Therefore, in this study, we aim to predict the ion channel conductance that is altered from various AP shapes using a machine learning algorithm. We perform electrophysiological simulations using a single-cell model to obtain AP shapes based on variations in the ion channel conductance. In the AP simulation, we increase and decrease the conductance of each ion channel at a constant rate, resulting in 1,980 AP shapes and one standard AP shape without any changes in the ion channel conductance. Subsequently, we calculate the AP difference shapes between them and use them as the input of the machine learning model to predict the changed ion channel conductance. In this study, we demonstrate that the changed ion channel conductance can be predicted with high prediction accuracy, as reflected by an F1 score of 0.985, using only AP shapes and simple machine learning.
| S-EPMC8035260 | biostudies-literature
Project description:Tau is a microtubule-associated protein that promotes microtubule assembly and stability. This protein is implicated in several neurodegenerative diseases, including Alzheimer's. To date, the three-dimensional (3D) structure of tau has not been fully solved, experimentally. Even the most recent information is sometimes controversial in regard to how this protein folds, interacts, and behaves. Predicting the tau structure and its profile sheds light on the knowledge about its properties and biological function, such as the binding to microtubules (MT) and, for instance, the effect on ionic conductivity. Our findings on the tau structure suggest a disordered protein, with discrete portions of well-defined secondary structure, mostly at the microtubule binding region. In addition, the first molecular dynamics simulation of full-length tau along with an MT section was performed, unveiling tau structure when associated with MT and interaction sites. Electrostatics and conductivity were also examined to understand how tau affects the ions in the intracellular fluid environment. Our results bring a new insight into tau and tubulin MT proteins, their characteristics, and the structure?function relationship.
| S-EPMC6468555 | biostudies-literature
Project description:Channelrhodopsin-2 (ChR2) is a light-activated and non-selective cationic channel protein that can be easily expressed in specific neurons to control neuronal activity by light. Although ChR2 has been extensively used as an optogenetic tool in neuroscience research, the molecular mechanism of cation channel formation following retinal photoisomerization in ChR2 is not well understood. In this paper, studies of the closed and opened state ChR2 structures are presented. The formation of the cationic channel is elucidated in atomic detail using molecular dynamics simulations on the all-trans-retinal (ChR2-trans) configuration of ChR2 and its isomerization products, 13-cis-retinal (ChR2-cis) configuration, respectively. Photoisomerization of the retinal-chromophore causes the destruction of interactions among the crucial residues (e.g., E90, E82, N258, and R268) around the channel and the extended H-bond network mediated by numerous water molecules, which opens the pore. Steering molecular dynamics (SMD) simulations show that the electrostatic interactions at the binding sites in intracellular gate (ICG) and central gate (CG) can influence the transmembrane transport of Na+ in ChR2-cis obviously. Potential of mean force (PMF) constructed by SMD and umbrella sampling also found the existing energy wells at these two binding sites during the transportation of Na+. These wells partly hinder the penetration of Na+ into cytoplasm through the ion channel. This investigation provides a theoretical insight on the formation mechanism of ion channels and the mechanism of ion permeation.
| S-EPMC6695816 | biostudies-literature
Project description:We investigate the validity of continuum electrostatics in the gramicidin A channel using a recently determined high-resolution structure. The potential and electric field acting on ions in and around the channel are computed by solving Poisson's equation. These are then used in Brownian dynamics simulations to obtain concentration profiles and the current passing through the channel. We show that regardless of the effective dielectric constant used for water in the channel or the channel protein, it is not possible to reproduce all the experimental data on gramicidin A; thus, continuum electrostatics cannot provide a valid framework for the description of ion dynamics in gramicidin channels. Using experimental data and molecular dynamics simulations as guides, we have constructed potential energy profiles that can satisfactorily describe the available physiological data. These profiles provide useful benchmarks for future potential of mean force calculations of permeating ions from molecular dynamics simulations of gramicidin A. They also offer a convenient starting point for studying structure-function relationships in modified gramicidin channels.
| S-EPMC1302233 | biostudies-other
Project description:Chaotic spiral or scroll wave dynamics can be found in diverse systems. In cardiac dynamics, spiral or scroll waves of electrical excitation determine the dynamics during life-threatening arrhythmias like ventricular fibrillation. In numerical studies it was found that chaotic episodes of spiral and scroll waves can be transient, thus they terminate spontaneously. We show in this study that this behavior can also be observed using models which describe the ion channel dynamics of human cardiomyocytes (Bueno-Orovio-Cherry-Fenton model and the Ten Tusscher-Noble-Noble-Panfilov model). For both models we find that the average lifetime of the chaotic transients grows exponentially with the system size. With this behavior, we classify the systems into the group of type-II supertransients. We observe a significant difference of the breakup behavior between the models, which results in a distinct dynamics during the final phase just before the termination. The observation of a (temporally) stable single-spiral state affects the prevailing description of the dynamics of type-II supertransients as being "quasi-stationary" and also the feasibility of predicting the spontaneous termination of the spiral wave dynamics. In the long term, the relation between the breakup behavior of spiral waves and properties of chaotic transients like predictability or average transient lifetime may contribute to an improved understanding and classification of cardiac arrhythmias.
| S-EPMC6713330 | biostudies-literature
Project description:Transient receptor potential ankyrin 1 (TRPA1) is a transmembrane protein channeling the influx of calcium ions. As a polymodal nocisensor, TRPA1 can be activated by thermal, mechanical stimuli and a wide range of chemically damaging molecules including small volatile environmental toxicants and endogenous algogenic lipids. After activation by such compounds, the ion channel opens up, its central pore widens allowing calcium influx into the cytosol inducing signal transduction pathways. Afterwards, the calcium influx desensitizes irritant evoked responses and results in an inactive state of the ion channel. Recent experimental determination of structures of apo and holo forms of TRPA1 opened the way towards the design of new agonists, which can activate the ion channel. The present study is aimed at the elucidation of binding dynamics of agonists using experimental structures of TRPA1-agonist complexes at the atomic level applying molecular docking and dynamics methods accounting for covalent and non-covalent interactions. Following a test of docking methods focused on the final, holo structures, prerequisite binding modes were detected involving the apo forms. It was shown how reversible interactions with prerequisite binding sites contribute to structural changes of TRPA1 leading to covalent bonding of agonists. The proposed dynamics of action allowed a mechanism-based forecast of new, druggable binding sites of potent agonists.
| S-EPMC8540651 | biostudies-literature
Project description:Bacterial Gloeobacter violaceus pentameric ligand-gated ion channel (GLIC) is activated to cation permeation upon lowering the solution pH. Its function can be modulated by anesthetic halothane. In the present work, we integrate molecular dynamics (MD) and Brownian dynamics (BD) simulations to elucidate the ion conduction, charge selectivity, and halothane modulation mechanisms in GLIC, based on recently resolved X-ray crystal structures of the open-channel GLIC. MD calculations of the potential of mean force (PMF) for a Na(+) revealed two energy barriers in the extracellular domain (R109 and K38) and at the hydrophobic gate of transmembrane domain (I233), respectively. An energy well for Na(+) was near the intracellular entrance: the depth of this energy well was modulated strongly by the protonation state of E222. The energy barrier for Cl(-) was found to be 3-4 times higher than that for Na(+). Ion permeation characteristics were determined through BD simulations using a hybrid MD/continuum electrostatics approach to evaluate the energy profiles governing the ion movement. The resultant channel conductance and a near-zero permeability ratio (P(Cl)/P(Na)) were comparable to experimental data. On the basis of these calculations, we suggest that a ring of five E222 residues may act as an electrostatic gate. In addition, the hydrophobic gate region may play a role in charge selectivity due to a higher dehydration energy barrier for Cl(-) ions. The effect of halothane on the Na(+) PMF was also evaluated. Halothane was found to perturb salt bridges in GLIC that may be crucial for channel gating and open-channel stability, but had no significant impact on the single ion PMF profiles.
| S-EPMC3071019 | biostudies-literature