1
|
Wilson MA, Pohorille A. Structure and Computational Electrophysiology of Ac-LS3, a Synthetic Ion Channel. J Phys Chem B 2022; 126:8985-8999. [PMID: 36306164 DOI: 10.1021/acs.jpcb.2c05965] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/11/2023]
Abstract
Computer simulations are reported on Ac-LS3, a synthetic ion channel, containing 21 residues with a Leu-Ser-Ser-Leu-Leu-Ser-Leu heptad repeat, which forms ions channels upon application of voltage. A hexameric, coiled-coil bundle initially positioned perpendicular to the membrane settled into a stable, tilted structure after 1.5 μs, most likely to improve contacts between the non-polar exterior of the channel and the hydrophobic core of the membrane. Once tilted, the bundle remained in this state during subsequent simulations of nearly 10 μs at voltages ranging from 200 to -100 mV. In contrast, attempts to identify a stable pentameric structure failed, thus supporting the hypothesis that the channel is a hexamer. Results at 100 mV were used to reconstruct the free energy profiles for K+ and Cl- in the channel. This was done by way of several methods in which results of molecular dynamics (MD) simulations were combined with the electrodiffusion model. Two of them developed recently do not require knowledge of the diffusivity. Instead, they utilize one-sided density profiles and committor probabilities. The consistency between different methods is very good, supporting the utility of the newly developed methods for reconstructing free energies of ions in channels. The flux of K+, which accounts for most of the current through the channel, calculated directly from MD matches well the total measured current. However, the current of Cl- is somewhat overestimated, possibly due to a slightly unbalanced force field involving chloride. The current-voltage dependence was also reconstructed by way of a recently developed, efficient method that requires simulations only at a single voltage, yielding good agreement with the experiment. Taken together, the results demonstrate that computational electrophysiology has become a reliable tool for studying how channels mediate ion transport through membranes.
Collapse
Affiliation(s)
- Michael A Wilson
- Exobiology Branch, MS239-4, NASA Ames Research Center, Moffett Field, California94035, United States.,SETI Institute, 189 Bernardo Avenue, Suite 200, Mountain View, California94043, United States
| | - Andrew Pohorille
- Exobiology Branch, MS239-4, NASA Ames Research Center, Moffett Field, California94033, United States.,Department of Pharmaceutical Chemistry, University of California San Francisco, San Francisco, California94132, United States
| |
Collapse
|
2
|
Lin YC, Luo YL. Unifying Single-Channel Permeability From Rare-Event Sampling and Steady-State Flux. Front Mol Biosci 2022; 9:860933. [PMID: 35495625 PMCID: PMC9043130 DOI: 10.3389/fmolb.2022.860933] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2022] [Accepted: 03/07/2022] [Indexed: 11/18/2022] Open
Abstract
Various all-atom molecular dynamics (MD) simulation methods have been developed to compute free energies and crossing rates of ions and small molecules through ion channels. However, a systemic comparison across different methods is scarce. Using a carbon nanotube as a model of small conductance ion channel, we computed the single-channel permeability for potassium ion using umbrella sampling, Markovian milestoning, and steady-state flux under applied voltage. We show that a slightly modified inhomogeneous solubility-diffusion equation yields a single-channel permeability consistent with the mean first passage time (MFPT) based method. For milestoning, applying cylindrical and spherical bulk boundary conditions yield consistent MFPT if factoring in the effective bulk concentration. The sensitivity of the MFPT to the output frequency of collective variables is highlighted using the convergence and symmetricity of the inward and outward MFPT profiles. The consistent transport kinetic results from all three methods demonstrated the robustness of MD-based methods in computing ion channel permeation. The advantages and disadvantages of each technique are discussed, focusing on the future applications of milestoning in more complex systems.
Collapse
|
3
|
Guardiani C, Cecconi F, Chiodo L, Cottone G, Malgaretti P, Maragliano L, Barabash ML, Camisasca G, Ceccarelli M, Corry B, Roth R, Giacomello A, Roux B. Computational methods and theory for ion channel research. ADVANCES IN PHYSICS: X 2022; 7:2080587. [PMID: 35874965 PMCID: PMC9302924 DOI: 10.1080/23746149.2022.2080587] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/28/2022] [Accepted: 05/15/2022] [Indexed: 06/15/2023] Open
Abstract
Ion channels are fundamental biological devices that act as gates in order to ensure selective ion transport across cellular membranes; their operation constitutes the molecular mechanism through which basic biological functions, such as nerve signal transmission and muscle contraction, are carried out. Here, we review recent results in the field of computational research on ion channels, covering theoretical advances, state-of-the-art simulation approaches, and frontline modeling techniques. We also report on few selected applications of continuum and atomistic methods to characterize the mechanisms of permeation, selectivity, and gating in biological and model channels.
Collapse
Affiliation(s)
- C. Guardiani
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Rome, Italy
| | - F. Cecconi
- CNR - Istituto dei Sistemi Complessi, Rome, Italy and Istituto Nazionale di Fisica Nucleare, INFN, Roma1 section. 00185, Roma, Italy
| | - L. Chiodo
- Department of Engineering, Campus Bio-Medico University, Rome, Italy
| | - G. Cottone
- Department of Physics and Chemistry-Emilio Segrè, University of Palermo, Palermo, Italy
| | - P. Malgaretti
- Helmholtz Institute Erlangen-Nürnberg for Renewable Energy (IEK-11), Forschungszentrum Jülich, Erlangen, Germany
| | - L. Maragliano
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Ancona, Italy, and Center for Synaptic Neuroscience and Technology, Istituto Italiano di Tecnologia, Genova, Italy
| | - M. L. Barabash
- Department of Materials Science and Nanoengineering, Rice University, Houston, TX 77005, USA
| | - G. Camisasca
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Rome, Italy
- Dipartimento di Fisica, Università Roma Tre, Rome, Italy
| | - M. Ceccarelli
- Department of Physics and CNR-IOM, University of Cagliari, Monserrato 09042-IT, Italy
| | - B. Corry
- Research School of Biology, The Australian National University, Canberra, ACT 2600, Australia
| | - R. Roth
- Institut Für Theoretische Physik, Eberhard Karls Universität Tübingen, Tübingen, Germany
| | - A. Giacomello
- Dipartimento di Ingegneria Meccanica e Aerospaziale, Sapienza Università di Roma, Rome, Italy
| | - B. Roux
- Department of Biochemistry & Molecular Biology, University of Chicago, Chicago IL, USA
| |
Collapse
|
4
|
Ray D, Stone SE, Andricioaei I. Markovian Weighted Ensemble Milestoning (M-WEM): Long-Time Kinetics from Short Trajectories. J Chem Theory Comput 2021; 18:79-95. [PMID: 34910499 DOI: 10.1021/acs.jctc.1c00803] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We introduce a rare-event sampling scheme, named Markovian Weighted Ensemble Milestoning (M-WEM), which inlays a weighted ensemble framework within a Markovian milestoning theory to efficiently calculate thermodynamic and kinetic properties of long-time-scale biomolecular processes from short atomistic molecular dynamics simulations. M-WEM is tested on the Müller-Brown potential model, the conformational switching in alanine dipeptide, and the millisecond time-scale protein-ligand unbinding in a trypsin-benzamidine complex. Not only can M-WEM predict the kinetics of these processes with quantitative accuracy but it also allows for a scheme to reconstruct a multidimensional free-energy landscape along additional degrees of freedom, which are not part of the milestoning progress coordinate. For the ligand-receptor system, the experimental residence time, association and dissociation kinetics, and binding free energy could be reproduced using M-WEM within a simulation time of a few hundreds of nanoseconds, which is a fraction of the computational cost of other currently available methods, and close to 4 orders of magnitude less than the experimental residence time. Due to the high accuracy and low computational cost, the M-WEM approach can find potential applications in kinetics and free-energy-based computational drug design.
Collapse
Affiliation(s)
- Dhiman Ray
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Sharon Emily Stone
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States
| | - Ioan Andricioaei
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States.,Department of Physics and Astronomy, University of California Irvine, Irvine, California 92697, United States
| |
Collapse
|
5
|
Jiang W, Lin YC, Botello-Smith W, Contreras JE, Harris AL, Maragliano L, Luo YL. Free energy and kinetics of cAMP permeation through connexin26 via applied voltage and milestoning. Biophys J 2021; 120:2969-2983. [PMID: 34214529 DOI: 10.1016/j.bpj.2021.06.024] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2021] [Revised: 04/08/2021] [Accepted: 06/17/2021] [Indexed: 11/18/2022] Open
Abstract
The connexin family is a diverse group of highly regulated wide-pore channels permeable to biological signaling molecules. Despite the critical roles of connexins in mediating selective molecular signaling in health and disease, the basis of molecular permeation through these pores remains unclear. Here, we report the thermodynamics and kinetics of binding and transport of a second messenger, adenosine-3',5'-cyclophosphate (cAMP), through a connexin26 hemichannel (Cx26). First, inward and outward fluxes of cAMP molecules solvated in KCl solution were obtained from 4 μs of ± 200 mV simulations. These fluxes data yielded a single-channel permeability of cAMP and cAMP/K+ permeability ratio consistent with experimentally measured values. The results from voltage simulations were then compared with the potential of mean force (PMF) and the mean first passage times (MFPTs) of a single cAMP without voltage, obtained from a total of 16.5 μs of Voronoi-tessellated Markovian milestoning simulations. Both the voltage simulations and the milestoning simulations revealed two cAMP-binding sites, for which the binding constants KD and dissociation rates koff were computed from PMF and MFPTs. The protein dipole inside the pore produces an asymmetric PMF, reflected in unequal cAMP MFPTs in each direction once within the pore. The free energy profiles under opposite voltages were derived from the milestoning PMF and revealed the interplay between voltage and channel polarity on the total free energy. In addition, we show how these factors influence the cAMP dipole vector during permeation, and how cAMP affects the local and nonlocal pore diameter in a position-dependent manner.
Collapse
Affiliation(s)
- Wenjuan Jiang
- Department of Pharmaceutical Sciences, College of Pharmacy, Western University of Health Sciences, Pomona, California
| | - Yi-Chun Lin
- Department of Pharmaceutical Sciences, College of Pharmacy, Western University of Health Sciences, Pomona, California
| | - Wesley Botello-Smith
- Department of Pharmaceutical Sciences, College of Pharmacy, Western University of Health Sciences, Pomona, California
| | - Jorge E Contreras
- Department of Physiology and Membrane Biology, School of Medicine, University of California, Davis, California.
| | - Andrew L Harris
- Department of Pharmacology, Physiology, and Neuroscience, New Jersey Medical School, Rutgers, The State University of New Jersey, Newark, New Jersey.
| | - Luca Maragliano
- Department of Life and Environmental Sciences, Polytechnic University of Marche, Ancona, Italy; Center for Synaptic Neuroscience and Technology, Italian Institute of Technology, Genoa, Italy.
| | - Yun Lyna Luo
- Department of Pharmaceutical Sciences, College of Pharmacy, Western University of Health Sciences, Pomona, California.
| |
Collapse
|
6
|
Wilson MA, Pohorille A. Electrophysiological Properties from Computations at a Single Voltage: Testing Theory with Stochastic Simulations. ENTROPY 2021; 23:e23050571. [PMID: 34066581 PMCID: PMC8148522 DOI: 10.3390/e23050571] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Revised: 04/24/2021] [Accepted: 04/28/2021] [Indexed: 12/13/2022]
Abstract
We use stochastic simulations to investigate the performance of two recently developed methods for calculating the free energy profiles of ion channels and their electrophysiological properties, such as current–voltage dependence and reversal potential, from molecular dynamics simulations at a single applied voltage. These methods require neither knowledge of the diffusivity nor simulations at multiple voltages, which greatly reduces the computational effort required to probe the electrophysiological properties of ion channels. They can be used to determine the free energy profiles from either forward or backward one-sided properties of ions in the channel, such as ion fluxes, density profiles, committor probabilities, or from their two-sided combination. By generating large sets of stochastic trajectories, which are individually designed to mimic the molecular dynamics crossing statistics of models of channels of trichotoxin, p7 from hepatitis C and a bacterial homolog of the pentameric ligand-gated ion channel, GLIC, we find that the free energy profiles obtained from stochastic simulations corresponding to molecular dynamics simulations of even a modest length are burdened with statistical errors of only 0.3 kcal/mol. Even with many crossing events, applying two-sided formulas substantially reduces statistical errors compared to one-sided formulas. With a properly chosen reference voltage, the current–voltage curves can be reproduced with good accuracy from simulations at a single voltage in a range extending for over 200 mV. If possible, the reference voltages should be chosen not simply to drive a large current in one direction, but to observe crossing events in both directions.
Collapse
Affiliation(s)
- Michael A. Wilson
- Exobiology Branch, MS 239-4, NASA Ames Research Center, Moffett Field, CA 94035, USA;
- SETI Institute, 189 Bernardo Ave, Suite 200, Mountain View, CA 94043, USA
| | - Andrew Pohorille
- Exobiology Branch, MS 239-4, NASA Ames Research Center, Moffett Field, CA 94035, USA;
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94132, USA
- Correspondence: ; Tel.: +1-650-604-5759
| |
Collapse
|
7
|
Alberini G, Benfenati F, Maragliano L. Structural Mechanism of ω-Currents in a Mutated Kv7.2 Voltage Sensor Domain from Molecular Dynamics Simulations. J Chem Inf Model 2021; 61:1354-1367. [PMID: 33570938 PMCID: PMC8023575 DOI: 10.1021/acs.jcim.0c01407] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
![]()
Activation of voltage-gated
ion channels is regulated by conformational
changes of the voltage sensor domains (VSDs), four water- and ion-impermeable
modules peripheral to the central, permeable pore domain. Anomalous
currents, defined as ω-currents, have been recorded in response
to mutations of residues on the VSD S4 helix and associated with ion
fluxes through the VSDs. In humans, gene defects in the potassium
channel Kv7.2 result in a broad range of epileptic disorders, from
benign neonatal seizures to severe epileptic encephalopathies. Experimental
evidence suggests that the R207Q mutation in S4, associated with peripheral
nerve hyperexcitability, induces ω-currents at depolarized potentials,
but the fine structural details are still elusive. In this work, we
use atom-detailed molecular dynamics simulations and a refined model
structure of the Kv7.2 VSD in the active conformation in a membrane/water
environment to study the effect of R207Q and four additional mutations
of proven clinical importance. Our results demonstrate that the R207Q
mutant shows the most pronounced increase of hydration in the internal
VSD cavity, a feature favoring the occurrence of ω-currents.
Free energy and kinetics calculations of sodium permeation through
the native and mutated VSD indicate as more favorable the formation
of a cationic current in the latter. Overall, our simulations establish
a mechanistic linkage between genetic variations and their physiological
outcome, by providing a computational description that includes both
thermodynamic and kinetic features of ion permeation associated with
ω-currents.
Collapse
Affiliation(s)
- Giulio Alberini
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy.,Department of Experimental Medicine, Università degli Studi di Genova, Viale Benedetto XV, 3, 16132 Genova, Italy
| | - Fabio Benfenati
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy.,IRCCS Ospedale Policlinico San Martino, Largo Rosanna Benzi, 10, 16132 Genova, Italy
| | - Luca Maragliano
- Center for Synaptic Neuroscience and Technology (NSYN@UniGe), Istituto Italiano di Tecnologia, Largo Rosanna Benzi, 10, 16132 Genova, Italy.,Department of Life and Environmental Sciences, Polytechnic University of Marche, Via Brecce Bianche, 60131 Ancona, Italy
| |
Collapse
|
8
|
Abstract
Mechanosensing is a key feature through which organisms can receive inputs from the environment and convert them into specific functional and behavioral outputs. Mechanosensation occurs in many cells and tissues, regulating a plethora of molecular processes based on the distribution of forces and stresses both at the cell membrane and at the intracellular organelles levels, through complex interactions between cells’ microstructures, cytoskeleton, and extracellular matrix. Although several primary and secondary mechanisms have been shown to contribute to mechanosensation, a fundamental pathway in simple organisms and mammals involves the presence of specialized sensory neurons and the presence of different types of mechanosensitive ion channels on the neuronal cell membrane. In this contribution, we present a review of the main ion channels which have been proven to be significantly involved in mechanotransduction in neurons. Further, we discuss recent studies focused on the biological mechanisms and modeling of mechanosensitive ion channels’ gating, and on mechanotransduction modeling at different scales and levels of details.
Collapse
|
9
|
Rao S, Klesse G, Lynch CI, Tucker SJ, Sansom MSP. Molecular Simulations of Hydrophobic Gating of Pentameric Ligand Gated Ion Channels: Insights into Water and Ions. J Phys Chem B 2021; 125:981-994. [PMID: 33439645 PMCID: PMC7869105 DOI: 10.1021/acs.jpcb.0c09285] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2020] [Revised: 12/13/2020] [Indexed: 12/30/2022]
Abstract
Ion channels are proteins which form gated nanopores in biological membranes. Many channels exhibit hydrophobic gating, whereby functional closure of a pore occurs by local dewetting. The pentameric ligand gated ion channels (pLGICs) provide a biologically important example of hydrophobic gating. Molecular simulation studies comparing additive vs polarizable models indicate predictions of hydrophobic gating are robust to the model employed. However, polarizable models suggest favorable interactions of hydrophobic pore-lining regions with chloride ions, of relevance to both synthetic carriers and channel proteins. Electrowetting of a closed pLGIC hydrophobic gate requires too high a voltage to occur physiologically but may inform designs for switchable nanopores. Global analysis of ∼200 channels yields a simple heuristic for structure-based prediction of (closed) hydrophobic gates. Simulation-based analysis is shown to provide an aid to interpretation of functional states of new channel structures. These studies indicate the importance of understanding the behavior of water and ions within the nanoconfined environment presented by ion channels.
Collapse
Affiliation(s)
- Shanlin Rao
- Department
of Biochemistry, University of Oxford, Oxford, U.K.
| | - Gianni Klesse
- Clarendon
Laboratory, Department of Physics, University
of Oxford, Oxford, U.K.
| | | | - Stephen J. Tucker
- Clarendon
Laboratory, Department of Physics, University
of Oxford, Oxford, U.K.
| | | |
Collapse
|