Gating Charge Calculations by Computational Electrophysiology Simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Machtens JP
PROVIDER: S-EPMC5389965 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:A series of ab initio (density functional) calculations were carried out on side chains of a set of amino acids, plus water, from the (intracellular) gating region of the KcsA K(+) channel. Their atomic coordinates, except hydrogen, are known from X-ray structures [D.A. Doyle, J.M. Cabral, R.A. Pfuetzner, A. Kuo, J.M. Gulbis, S.L. Cohen, B.T. Chait, R. MacKinnon, The structure of the potassium channel: molecular basis of K(+) conduction and selectivity, Science 280 (1998) 69-77; R. MacKinnon, S.L. Cohen, A. Kuo, A. Lee, B.T. Chait, Structural conservation in prokaryotic and eukaryotic potassium channels, Science 280 (1998) 106-109; Y. Jiang, A. Lee, J. Chen, M. Cadene, B.T. Chait, R. MacKinnon, The open pore conformation of potassium channels. Nature 417 (2001) 523-526], as are the coordinates of some water oxygen atoms. The 1k4c structure is used for the starting coordinates. Quantum mechanical optimization, in spite of the starting configuration, places the atoms in positions much closer to the 1j95, more tightly closed, configuration. This state shows four water molecules forming a "basket" under the Q119 side chains, blocking the channel. When a hydrated K(+) approaches this "basket", the optimized system shows a strong set of hydrogen bonds with the K(+) at defined positions, preventing further approach of the K(+) to the basket. This optimized structure with hydrated K(+) added shows an ice-like 12 molecule nanocrystal of water. If the water molecules exchange, unless they do it as a group, the channel will remain blocked. The "basket" itself appears to be very stable, although it is possible that the K(+) with its hydrating water molecules may be more mobile, capable of withdrawing from the gate. It is also not surprising that water essentially freezes, or forms a kind of glue, in a nanometer space; this agrees with experimental results on a rather different, but similarly sized (nm dimensions) system [K.B. Jinesh, J.W.M. Frenken, Capillary condensation in atomic scale friction: how water acts like a glue, Phys. Rev. Lett. 96 (2006) 166103/1-4]. It also agrees qualitatively with simulations on channels [A. Anishkin, S. Sukharev, Water dynamics and dewetting transitions in the small mechanosensitive channel MscS, Biophys. J. 86 (2004) 2883-2895; O. Beckstein, M.S.P. Sansom, Liquid-vapor oscillations of water in hydrophobic nanopores, Proc. Natl Acad. Sci. U. S. A. 100 (2003) 7063-7068] and on featureless channel-like systems [J. Lu, M.E. Green, Simulation of water in a pore with charges: application to a gating mechanism for ion channels, Prog. Colloid Polym. Sci. 103 (1997) 121-129], in that it forms a boundary on water that is not obvious from the liquid state. The idea that a structure is stable, even if individual molecules exchange, is well known, for example from the hydration shell of ions. We show that when charges are added in the form of protons to the domains (one proton per domain), the optimized structure is open. No stable water hydrogen bonds hold it together; an opening of 11.0 A appears, measured diagonally between non-neighboring domains as glutamine 119 carbonyl O-O distance. This is comparable to the opening in the MthK potassium channel structure that is generally agreed to be open. The appearance of the opening is in rather good agreement with that found by Perozo and coworkers. In contrast, in the uncharged structure this diagonal distance is 6.5 A, and the water "basket" constricts the uncharged opening still further, with the ice-like structure that couples the K(+) ion to the gating region freezing the entrance to the channel. Comparison with our earlier model for voltage gated channels suggests that a similar mechanism may apply in those channels.
| S-EPMC1989151 | biostudies-literature
Project description:Dynorphins are endogenous neuropeptides that function as ligands for the κ-opioid receptor. In addition to opioid activity, dynorphins can induce several pathological effects such as neurological dysfunctions and cell death. Previous studies have suggested that Dynorphin A (DynA) mediates some pathogenic actions through formation of transient pores in lipid domains of the plasma membrane. Here, we use planar bilayer electrophysiology to show that DynA induces pore formation in negatively charged membranes. We find a large variability in pore conformations showing equilibrium conductance fluctuations, what disregards electroporation as the dominant mechanism of pore formation. Ion selectivity measurements showing cationic selectivity indicate that positive protein charges of DynA are stabilized by phosphatidyl serine negative charges in the formation of combined structures. We complement our study with computational simulations that assess the stability of diverse peptide arrangements in the hydrophobic core of the bilayer. We show that DynA is capable of assembling in charged membranes to form water-filled pores that conduct ions.
| S-EPMC8718563 | biostudies-literature
Project description:Developing an understanding of the mechanism of voltage-gated ion channels in molecular terms requires knowledge of the structure of the active and resting conformations. Although the active-state conformation is known from x-ray structures, an atomic resolution structure of a voltage-dependent ion channel in the resting state is not currently available. This has motivated various efforts at using computational modeling methods and molecular dynamics (MD) simulations to provide the missing information. A comparison of recent computational results reveals an emerging consensus on voltage-dependent gating from computational modeling and MD simulations. This progress is highlighted in the broad context of preexisting work about voltage-gated channels.
| S-EPMC3514734 | biostudies-literature
Project description:Ionotropic glutamate receptors (iGluRs) are tetrameric ligand-gated ion channels essential to all aspects of brain function, including higher order processes such as learning and memory. For decades, electrophysiology was the primary means for characterizing the function of iGluRs and gaining mechanistic insight. Since the turn of the century, structures of isolated water-soluble domains and transmembrane-domain-containing constructs have provided the basis for formulating mechanistic hypotheses. Because these structures only represent sparse, often incomplete snapshots during iGluR activation, significant gaps in knowledge remain regarding structures, energetics, and dynamics of key substates along the functional processes. Some of these gaps have recently been filled by molecular dynamics simulations and theoretical modeling. In this Account, I describe our work in the latter arena toward characterizing iGluR gating motions and developing a formalism for calculating thermodynamic and kinetic properties of stationary gating. The structures of iGluR subunits have a highly modular architecture, in which the ligand-binding domain and the transmembrane domain are well separated and connected by flexible linkers. The ligand-binding domain in turn is composed of two subdomains. During activation, agonist binding induces the closure of the intersubdomain cleft. The cleft closure leads to the outward pulling of a linker tethered to the extracellular terminus of the major pore-lining helix of the transmembrane domain, thereby opening the channel. This activation model based on molecular dynamics simulations was validated by residue-specific information from electrophysiological data on cysteine mutants. A further critical test was made through introducing glycine insertions in the linker. Molecular dynamics simulations showed that, with lengthening by glycine insertions, the linker became less effective in pulling the pore-lining helix, leading to weaker stabilization of the channel-open state. In full agreement, single-channel recordings showed that the channel open probability decreased progressively as the linker was lengthened by glycine insertions. Crystal structures of ligand-binding domains showing different degrees of cleft closure between full and partial agonists suggested a simple mechanism for one subtype of iGluRs, but mysteries surrounded a second subtype, where the ligand-binding domains open to similar degrees when bound with either full or partial agonists. Our free energy simulations now suggest that broadening of the free energy basin for cleft closure is a plausible solution. A theoretical basis for these mechanistic hypotheses on partial agonisms was provided by a model for the free energy surface of a full receptor, where the stabilization by cleft closure is transmitted via the linker to the channel-open state. This model can be implemented by molecular dynamics simulations to predict thermodynamic and kinetics properties of stationary gating that are amenable to direct test by single-channel recordings. Close integration between computation and electrophysiology holds great promises in revealing the conformations of key substates in functional processes and the mechanisms of disease-associated mutations.
| S-EPMC5398286 | biostudies-literature
Project description:We here investigate the Electronic Circular Dichroism (ECD) Spectra of two representative Guanine-rich sequences folded in a Quadruple helix (GQ), by using a recently developed fragment diabatisation based excitonic model (FrDEx). FrDEx can include charge transfer (CT) excited states and consider the effect of the surrounding monomers on the local excitations (LEs). When applied to different structures generated by molecular dynamics simulations on a fragment of the human telomeric sequence (Tel21/22), FrDEx provides spectra fully consistent with the experimental one and in good agreement with that provided by quantum mechanical (QM) method used for its parametrization, i.e., TD-M05-2X. We show that the ECD spectrum is moderately sensitive to the conformation adopted by the bases of the loops and more significantly to the thermal fluctuations of the Guanine tetrads. In particular, we show how changes in the overlap of the tetrads modulate the intensity of the ECD signal. We illustrate how this correlates with changes in the character of the excitonic states at the bottom of the La and Lb bands, with larger LE and CT involvement of bases that are more closely stacked. As an additional test, we utilised FrDEx to compute the ECD spectrum of the monomeric and dimeric forms of a GQ forming sequence T30695 (5'TGGGTGGGTGGGTGGG3'), i.e., a system containing up to 24 Guanine bases, and demonstrated the satisfactory reproduction of the experimental and QM reference results. This study provides new insights on the effects modulating the ECD spectra of GQs and, more generally, further validates FrDEx as an effective tool to predict and assign the spectra of closely stacked multichromophore systems.
| S-EPMC8398971 | biostudies-literature
Project description:Voltage sensors regulate the conformations of voltage-dependent ion channels and enzymes. Their nearly switchlike response as a function of membrane voltage comes from the movement of positively charged amino acids, arginine or lysine, across the membrane field. We used mutations with natural and unnatural amino acids, electrophysiological recordings, and x-ray crystallography to identify a charge transfer center in voltage sensors that facilitates this movement. This center consists of a rigid cyclic "cap" and two negatively charged amino acids to interact with a positive charge. Specific mutations induce a preference for lysine relative to arginine. By placing lysine at specific locations, the voltage sensor can be stabilized in different conformations, which enables a dissection of voltage sensor movements and their relation to ion channel opening.
| S-EPMC2869078 | biostudies-literature
Project description:KV11.1 voltage-gated K(+) channels are noted for unusually slow activation, fast inactivation, and slow deactivation kinetics, which tune channel activity to provide vital repolarizing current during later stages of the cardiac action potential. The bulk of charge movement in human ether-a-go-go-related gene (hERG) is slow, as is return of charge upon repolarization, suggesting that the rates of hERG channel opening and, critically, that of deactivation might be determined by slow voltage sensor movement, and also by a mode-shift after activation. To test these ideas, we compared the kinetics and voltage dependence of ionic activation and deactivation with gating charge movement. At 0 mV, gating charge moved ∼threefold faster than ionic current, which suggests the presence of additional slow transitions downstream of charge movement in the physiological activation pathway. A significant voltage sensor mode-shift was apparent by 24 ms at +60 mV in gating currents, and return of charge closely tracked pore closure after pulses of 100 and 300 ms duration. A deletion of the N-terminus PAS domain, mutation R4AR5A or the LQT2-causing mutation R56Q gave faster-deactivating channels that displayed an attenuated mode-shift of charge. This indicates that charge movement is perturbed by N- and C-terminus interactions, and that these domain interactions stabilize the open state and limit the rate of charge return. We conclude that slow on-gating charge movement can only partly account for slow hERG ionic activation, and that the rate of pore closure has a limiting role in the slow return of gating charges.
| S-EPMC4375626 | biostudies-literature
Project description:Thermal stability of proteins is crucial for both biotechnological and therapeutic applications. Rational protein engineering therefore frequently aims at increasing thermal stability by introducing stabilizing mutations. The accurate prediction of the thermodynamic consequences caused by mutations, however, is highly challenging as thermal stability changes are caused by alterations in the free energy of folding. Growing computational power, however, increasingly allows us to use alchemical free energy simulations, such as free energy perturbation or thermodynamic integration, to calculate free energy differences with relatively high accuracy. In this article, we present an automated protocol for setting up alchemical free energy calculations for mutations of naturally occurring amino acids (except for proline) that allows an unprecedented, automated screening of large mutant libraries. To validate the developed protocol, we calculated thermodynamic stability differences for 109 mutations in the microbial Ribonuclease Barnase. The obtained quantitative agreement with experimental data illustrates the potential of the approach in protein engineering and design.
| S-EPMC2872215 | biostudies-literature
Project description:The pig is commonly used as an experimental model of human heart disease, including for the study of mechanisms of arrhythmia. However, there exist differences between human and porcine cellular electrophysiology: The pig action potential (AP) has a deeper phase-1 notch, a longer duration at 50% repolarization, and higher plateau potentials than human. Ionic differences underlying the AP include larger rapid delayed-rectifier and smaller inward-rectifier K+-currents (IKr and IK1 respectively) in humans. AP steady-state rate-dependence and restitution is steeper in pigs. Porcine Ca2+ transients can have two components, unlike human. Although a reliable computational model for human ventricular myocytes exists, one for pigs is lacking. This hampers translation from results obtained in pigs to human myocardium. Here, we developed a computational model of the pig ventricular cardiomyocyte AP using experimental datasets of the relevant ionic currents, Ca2+-handling, AP shape, AP duration restitution, and inducibility of triggered activity and alternans. To properly capture porcine Ca2+ transients, we introduced a two-step process with a faster release in the t-tubular region, followed by a slower diffusion-induced release from a non t-tubular subcellular region. The pig model behavior was compared with that of a human ventricular cardiomyocyte (O'Hara-Rudy) model. The pig, but not the human model, developed early afterdepolarizations (EADs) under block of IK1, while IKr block led to EADs in the human but not in the pig model. At fast rates (pacing cycle length = 400 ms), the human cell model was more susceptible to spontaneous Ca2+ release-mediated delayed afterdepolarizations (DADs) and triggered activity than pig. Fast pacing led to alternans in human but not pig. Developing species-specific models incorporating electrophysiology and Ca2+-handling provides a tool to aid translating antiarrhythmic and arrhythmogenic assessment from the bench to the clinic.
| S-EPMC8277015 | biostudies-literature
Project description:Fluorescent molecular rotors (FMRs) belong to an important class of environment-sensitive dyes capable of acting as nanoprobes in the measurement of viscosity and polarity of their micro-environment. FMRs have found widespread applications in various research fields, ranging from analytical to biochemical sciences, for example in intracellular imaging studies or in volatile organic compound detection. Here, a computational investigation of a recently proposed FMR, namely 4-(diphenylamino)phthalonitrile (DPAP), in various chemical environments is presented. A purposely developed molecular mechanics force field is proposed and then applied to simulate the rotor in a high- and low-polar solvent (i.e., acetonitrile, tetrahydrofuran, o-xylene and cyclohexane), a polymer matrix and a lipid membrane. Subtle effects of the molecular interactions with the embedding medium, the structural fluctuations of the rotor and its rotational dynamics are analyzed in some detail. The results correlate with a previous work, thus supporting the reliability of the model, and provide further insights into the environment-specific properties of the dye. In particular, it is shown how molecular diffusion and rotational correlation times of the FMR are affected by the surrounding medium and how the molecular orientation of the dye becomes anisotropic once immersed in the lipid bilayer. Moreover, a qualitative correlation between the FMR rotational dynamics and the fluorescence lifetime is detected, a result in line with the observed viscosity dependence of its emission. Finally, optical absorption spectra are computed and successfully compared with their experimental counterparts.
| S-EPMC5886372 | biostudies-literature