1
|
Zaki AM, Çınaroğlu SS, Rahman T, Patel S, Biggin PC. Plasticity of the selectivity filter is essential for permeation in lysosomal TPC2 channels. Proc Natl Acad Sci U S A 2024; 121:e2320153121. [PMID: 39074274 PMCID: PMC11317647 DOI: 10.1073/pnas.2320153121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2023] [Accepted: 06/24/2024] [Indexed: 07/31/2024] Open
Abstract
Two-pore channels are pathophysiologically important Na+- and Ca2+-permeable channels expressed in lysosomes and other acidic organelles. Unlike most other ion channels, their permeability is malleable and ligand-tuned such that when gated by the signaling lipid PI(3,5)P2, they are more Na+-selective than when gated by the Ca2+ mobilizing messenger nicotinic acid adenine dinucleotide phosphate. However, the structural basis that underlies such plasticity and single-channel behavior more generally remains poorly understood. A recent Cryo-electron microscopy (cryo-EM) structure of TPC2 bound to PI(3,5)P2 in a proposed open-channel conformation provided an opportunity to address this via molecular dynamics (MD) simulation. To our surprise, simulations designed to compute conductance through this structure revealed almost no Na+ permeation events even at very high transmembrane voltages. However further MD simulations identified a spontaneous transition to a dramatically different conformation of the selectivity filter that involved expansion and a flip in the orientation of two core asparagine residues. This alternative filter conformation was remarkably stable and allowed Na+ to flow through the channel leading to a conductance estimate that was in very good agreement with direct single-channel measurements. Furthermore, this conformation was more permeable for Na+ over Ca2+. Our results have important ramifications not just for understanding the control of ion selectivity in TPC2 channels but also more broadly in terms of how ion channels discriminate ions.
Collapse
Affiliation(s)
- Afroditi-Maria Zaki
- Department of Biochemistry, Structural Bioinformatics and Computational Biochemistry, University of Oxford, OxfordOX1 3QU, United Kingdom
| | - Süleyman Selim Çınaroğlu
- Department of Biochemistry, Structural Bioinformatics and Computational Biochemistry, University of Oxford, OxfordOX1 3QU, United Kingdom
| | - Taufiq Rahman
- Department of Pharmacology, University of Cambridge, CambridgeCB2 1PD, United Kingdom
| | - Sandip Patel
- Department of Cell and Developmental Biology, University College London, LondonWC1E, 6BT, United Kingdom
| | - Philip C. Biggin
- Department of Biochemistry, Structural Bioinformatics and Computational Biochemistry, University of Oxford, OxfordOX1 3QU, United Kingdom
| |
Collapse
|
2
|
Weckel-Dahman H, Carlsen R, Swanson JM. Multiscale Responsive Kinetic Modeling: Quantifying Biomolecular Reaction Flux under Varying Electrochemical Conditions. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.08.01.606205. [PMID: 39131358 PMCID: PMC11312519 DOI: 10.1101/2024.08.01.606205] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 08/13/2024]
Abstract
Attaining a complete thermodynamic and kinetic characterization for processes involving multiple interconnected rare-event transitions remains a central challenge in molecular biophysics. This challenge is amplified when the process must be understood under a range of reaction conditions. Herein, we present a condition-responsive kinetic modeling framework that can combine the strengths of bottom-up rate quantification from multiscale simulations with top-down solution refinement using experimental data. Although this framework can be applied to any process, we demonstrate its use for electrochemically driven transport through channels and transporters. Using the Cl- /H+ antiporter ClC-ec1 as a model system, we show how robust and predictive kinetic solutions can be obtained when the solution space is grounded by thermodynamic constraints, seeded through multiscale rate quantification, and further refined with experimental data, such as electrophysiology assays. Turning to the Shaker K+ channel, we demonstrate that robust solutions and biophysical insights can also be obtained with sufficient experimental data. This multi-pathway method proves capable of identifying single-pathway dominant mechanisms but also highlights that competing and off-pathway flux is still essential to replicate experimental findings and to describe concentration-dependent channel rectification.
Collapse
Affiliation(s)
- Hannah Weckel-Dahman
- Department of Chemistry, University of Utah, Salt Lake City, UT, 84112 – United States of America
| | - Ryan Carlsen
- Department of Chemistry, University of Utah, Salt Lake City, UT, 84112 – United States of America
| | - Jessica M.J. Swanson
- Department of Chemistry, University of Utah, Salt Lake City, UT, 84112 – United States of America
| |
Collapse
|
3
|
Chen C, Han P, Qing Y. Metabolic heterogeneity in tumor microenvironment - A novel landmark for immunotherapy. Autoimmun Rev 2024:103579. [PMID: 39004158 DOI: 10.1016/j.autrev.2024.103579] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Revised: 04/10/2024] [Accepted: 07/09/2024] [Indexed: 07/16/2024]
Abstract
The surrounding non-cancer cells and tumor cells that make up the tumor microenvironment (TME) have various metabolic rhythms. TME metabolic heterogeneity is influenced by the intricate network of metabolic control within and between cells. DNA, protein, transport, and microbial levels are important regulators of TME metabolic homeostasis. The effectiveness of immunotherapy is also closely correlated with alterations in TME metabolism. The response of a tumor patient to immunotherapy is influenced by a variety of variables, including intracellular metabolic reprogramming, metabolic interaction between cells, ecological changes within and between tumors, and general dietary preferences. Although immunotherapy and targeted therapy have made great strides, their use in the accurate identification and treatment of tumors still has several limitations. The function of TME metabolic heterogeneity in tumor immunotherapy is summarized in this article. It focuses on how metabolic heterogeneity develops and is regulated as a tumor progresses, the precise molecular mechanisms and potential clinical significance of imbalances in intracellular metabolic homeostasis and intercellular metabolic coupling and interaction, as well as the benefits and drawbacks of targeted metabolism used in conjunction with immunotherapy. This offers insightful knowledge and important implications for individualized tumor patient diagnosis and treatment plans in the future.
Collapse
Affiliation(s)
- Chen Chen
- The First Affiliated Hospital of Ningbo University, Ningbo 315211, Zhejiang, China
| | - Peng Han
- Harbin Medical University Cancer Hospital, Harbin 150081, Heilongjiang, China.
| | - Yanping Qing
- The First Affiliated Hospital of Ningbo University, Ningbo 315211, Zhejiang, China.
| |
Collapse
|
4
|
Takasawa T, Matsui T, Watanabe G, Kodera Y. Molecular dynamics simulations reveal differences in the conformational stability of FtsZs derived from Staphylococcus aureus and Bacillus subtilis. Sci Rep 2024; 14:16043. [PMID: 38992051 PMCID: PMC11239868 DOI: 10.1038/s41598-024-66763-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2024] [Accepted: 07/03/2024] [Indexed: 07/13/2024] Open
Abstract
FtsZ is highly conserved among bacteria and plays an essential role in bacterial cell division. The tense conformation of FtsZ bound to GTP assembles into a straight filament via head-to-tail associations, and then the upper subunit of FtsZ hydrolyzes GTP bound to the lower FtsZ subunit. The subunit with GDP bound disassembles accompanied by a conformational change in the subunit from the tense to relaxed conformation. Although crystal structures of FtsZ derived from several bacterial species have been determined, the conformational change from the relaxed to tense conformation has only been observed in Staphylococcus aureus FtsZ (SaFtsZ). Recent cryo-electron microscopy analyses revealed the three-dimensional reconstruction of the protofilament, in which tense molecules assemble via head-to-tail associations. However, the lower resolution of the protofilament suggested that the flexibility of the FtsZ protomers between the relaxed and tense conformations caused them to form in less-strict alignments. Furthermore, this flexibility may also prevent FtsZs other than SaFtsZ from crystalizing in the tense conformation, suggesting that the flexibility of bacterial FtsZs differs. In this study, molecular dynamics simulations were performed using SaFtsZ and Bacillus subtilis FtsZ in several situations, which suggested that different features of the FtsZs affect their conformational stability.
Collapse
Affiliation(s)
- Taichi Takasawa
- Department of Physics, School of Science, Kitasato University, 1-15-1 Kitasato, Minami-Ku, Sagamihara, Kanagawa, 252-0373, Japan
| | - Takashi Matsui
- Department of Physics, School of Science, Kitasato University, 1-15-1 Kitasato, Minami-Ku, Sagamihara, Kanagawa, 252-0373, Japan.
- Center for Disease Proteomics, School of Science, Kitasato University, 1-15-1 Kitasato, Minami-Ku, Sagamihara, Kanagawa, 252-0373, Japan.
| | - Go Watanabe
- Department of Data Science, School of Frontier Engineering, Kitasato University, 1-15-1 Kitasato, Minami-Ku, Sagamihara, Kanagawa, 252-0373, Japan.
- Kanagawa Institute of Industrial Science and Technology (KISTEC), 705-1 Shimoimaizumi, Ebina, Kanagawa, 243-0435, Japan.
| | - Yoshio Kodera
- Department of Physics, School of Science, Kitasato University, 1-15-1 Kitasato, Minami-Ku, Sagamihara, Kanagawa, 252-0373, Japan
- Center for Disease Proteomics, School of Science, Kitasato University, 1-15-1 Kitasato, Minami-Ku, Sagamihara, Kanagawa, 252-0373, Japan
| |
Collapse
|
5
|
Yang H, Zhou D, Zhou Z, Duan M, Yu H. Mechanistic Insight into the Mechanical Unfolding of the Integral Membrane Diacylglycerol Kinase. JACS AU 2024; 4:1422-1435. [PMID: 38665647 PMCID: PMC11040704 DOI: 10.1021/jacsau.3c00829] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/26/2023] [Revised: 02/29/2024] [Accepted: 03/01/2024] [Indexed: 04/28/2024]
Abstract
The essential forces stabilizing membrane proteins and governing their folding and unfolding are difficult to decipher. Single-molecule atomic force spectroscopy mechanically unfolds individual membrane proteins and quantifies their dynamics and energetics. However, it remains challenging to structurally assign unfolding intermediates precisely and to deduce dominant interactions between specific residues that facilitate either the localized stabilization of these intermediates or the global assembly of membrane proteins. Here, we performed force spectroscopy experiments and multiscale molecular dynamics simulations to study the unfolding pathway of diacylglycerol kinase (DGK), a small trimeric multispan transmembrane enzyme. The remarkable agreement between experiments and simulations allowed precise structural assignment and interaction analysis of unfolding intermediates, bypassing existing limitations on structural mapping, and thus provided mechanistic explanations for the formation of these states. DGK unfolding was found to proceed with structural segments varying in size that do not correlate with its secondary structure. We identified intermolecular side-chain packing interactions as one of the major contributions to the stability of unfolding intermediates. Mutagenesis creating packing defects induced a dramatic decrease in the mechano-stability of corresponding intermediates and also in the thermo-stability of DGK trimer, in good agreement with predictions from simulations. Hence, the molecular determinants of the mechano- and thermo-stability of a membrane protein can be identified at residue resolution. The accurate structural assignment established and microscopic mechanism revealed in this work may substantially expand the scope of single-molecule studies of membrane proteins.
Collapse
Affiliation(s)
- Huiying Yang
- School
of Physics, Huazhong University of Science
and Technology, Wuhan 430074, China
| | - Daihong Zhou
- School
of Physics, Huazhong University of Science
and Technology, Wuhan 430074, China
| | - Zhangyi Zhou
- School
of Physics, Huazhong University of Science
and Technology, Wuhan 430074, China
| | - Mojie Duan
- Innovation
Academy for Precision Measurement Science and Technology, Chinese
Academy of Sciences, Wuhan 430071, China
| | - Hao Yu
- School
of Physics, Huazhong University of Science
and Technology, Wuhan 430074, China
| |
Collapse
|
6
|
Mondal R, Vaissier Welborn V. Dynamics accelerate the kinetics of ion diffusion through channels: Continuous-time random walk models beyond the mean field approximation. J Chem Phys 2024; 160:144109. [PMID: 38597306 DOI: 10.1063/5.0188469] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Accepted: 03/14/2024] [Indexed: 04/11/2024] Open
Abstract
Ion channels are proteins that play a significant role in physiological processes, including neuronal excitability and signal transduction. However, the precise mechanisms by which these proteins facilitate ion diffusion through cell membranes are not well understood. This is because experimental techniques to characterize ion channel activity operate on a time scale too large to understand the role of the various protein conformations on diffusion. Meanwhile, computational approaches operate on a time scale too short to rationalize the observed behavior at the microscopic scale. In this paper, we present a continuous-time random walk model that aims to bridge the scales between the atomistic models of ion channels and the experimental measurement of their conductance. We show how diffusion slows down in complex systems by using 3D lattices that map out the pore geometry of two channels: Nav1.7 and gramicidin. We also introduce spatial and dynamic site disorder to account for system heterogeneity beyond the mean field approximation. Computed diffusion coefficients show that an increase in spatial disorder slows down diffusion kinetics, while dynamic disorder has the opposite effect. Our results imply that microscopic or phenomenological models based on the potential of mean force data overlook the functional importance of protein dynamics on ion diffusion through channels.
Collapse
Affiliation(s)
- Ronnie Mondal
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
- Macromolecules Innovation Institute, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Valerie Vaissier Welborn
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
- Macromolecules Innovation Institute, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
7
|
Fedida D, Sastre D, Dou Y, Westhoff M, Eldstrom J. Evaluating sequential and allosteric activation models in IKs channels with mutated voltage sensors. J Gen Physiol 2024; 156:e202313465. [PMID: 38294435 PMCID: PMC10829594 DOI: 10.1085/jgp.202313465] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 11/30/2023] [Accepted: 01/03/2024] [Indexed: 02/01/2024] Open
Abstract
The ion-conducting IKs channel complex, important in cardiac repolarization and arrhythmias, comprises tetramers of KCNQ1 α-subunits along with 1-4 KCNE1 accessory subunits and calmodulin regulatory molecules. The E160R mutation in individual KCNQ1 subunits was used to prevent activation of voltage sensors and allow direct determination of transition rate data from complexes opening with a fixed number of 1, 2, or 4 activatable voltage sensors. Markov models were used to test the suitability of sequential versus allosteric models of IKs activation by comparing simulations with experimental steady-state and transient activation kinetics, voltage-sensor fluorescence from channels with two or four activatable domains, and limiting slope currents at negative potentials. Sequential Hodgkin-Huxley-type models approximately describe IKs currents but cannot explain an activation delay in channels with only one activatable subunit or the hyperpolarizing shift in the conductance-voltage relationship with more activatable voltage sensors. Incorporating two voltage sensor activation steps in sequential models and a concerted step in opening via rates derived from fluorescence measurements improves models but does not resolve fundamental differences with experimental data. Limiting slope current data that show the opening of channels at negative potentials and very low open probability are better simulated using allosteric models of activation with one transition per voltage sensor, which implies that movement of all four sensors is not required for IKs conductance. Tiered allosteric models with two activating transitions per voltage sensor can fully account for IKs current and fluorescence activation kinetics in constructs with different numbers of activatable voltage sensors.
Collapse
Affiliation(s)
- David Fedida
- Department of Anesthesiology, Pharmacology and Therapeutics, University of British Columbia, Vancouver, Canada
| | - Daniel Sastre
- Department of Anesthesiology, Pharmacology and Therapeutics, University of British Columbia, Vancouver, Canada
| | - Ying Dou
- Department of Anesthesiology, Pharmacology and Therapeutics, University of British Columbia, Vancouver, Canada
| | - Maartje Westhoff
- Department of Anesthesiology, Pharmacology and Therapeutics, University of British Columbia, Vancouver, Canada
| | - Jodene Eldstrom
- Department of Anesthesiology, Pharmacology and Therapeutics, University of British Columbia, Vancouver, Canada
| |
Collapse
|
8
|
Stevens MJ, Rempe SLB. Insight into the K channel's selectivity from binding of K +, Na + and water to N-methylacetamide. Faraday Discuss 2024; 249:195-209. [PMID: 37846738 DOI: 10.1039/d3fd00110e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2023]
Abstract
In potassium channels that conduct K+ selectively over Na+, which sites are occupied by K+ or water and the mechanism of selectivity are unresolved questions. The combination of the energetics and the constraints imposed by the protein structure yield the selective permeation and occupancy. To gain insight into the combination of structure and energetics, we performed density functional theory (DFT) calculations of multiple N-methyl acetamide (NMA) ligands binding to K+ and Na+, relative to hydrated K+ and Na+. NMA is an analogue of the amino acid backbone and provides the carbonyl binding to the ions that occurs in most binding sites of the K+ channel. Unconstrained optimal structures are obtained through geometry optimization calculations of the NMA ligand binding. The complexes formed by 8 NMA binding to the cations have the O atoms positioned in nearly identical locations as the O atoms in the selectivity filter. The transfer free energies between bulk water and K+ or Na+ bound to 8 NMA are almost identical, implying there is no selectivity by a single site. For water optimized with 8 NMA, binding is weak and O atoms are not positioned as in the K+ channel selectivity filter, suggesting that the ions are much more favored than water. Optimal structures of 8 NMA binding with two cations (K+ or Na+) are stable and have lower binding free energy than the optimal structures with just one cation. However, in the Na+ case, the optimal structure deforms and does not match the K+ channel; that is, two bound Na+ are destabilizing. In contrast, the two K+ structure is stabilized and the selectivity free energy favors K+. Overall, this study shows that binding site occupancy and the mechanism for K+ selectivity involves multiple K+ binding in multiple neighboring layers or sites of the K+ channel selectivity filter.
Collapse
Affiliation(s)
- Mark J Stevens
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, NM 87185, USA.
| | - Susan L B Rempe
- Center for Integrated Nanotechnologies, Sandia National Laboratories, Albuquerque, NM 87185, USA.
| |
Collapse
|
9
|
Ngo K, Lopez Mateos D, Han Y, Rouen KC, Ahn SH, Wulff H, Clancy CE, Yarov-Yarovoy V, Vorobyov I. Elucidating molecular mechanisms of protoxin-II state-specific binding to the human NaV1.7 channel. J Gen Physiol 2024; 156:e202313368. [PMID: 38127314 PMCID: PMC10737443 DOI: 10.1085/jgp.202313368] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2023] [Revised: 09/08/2023] [Accepted: 12/04/2023] [Indexed: 12/23/2023] Open
Abstract
Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells, and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-II (PTx2) has high selectivity for hNaV1.7 and is a valuable scaffold for designing novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD and the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.
Collapse
Affiliation(s)
- Khoa Ngo
- Biophysics Graduate Group, University of California, Davis, Davis, CA, USA
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, USA
| | - Diego Lopez Mateos
- Biophysics Graduate Group, University of California, Davis, Davis, CA, USA
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, USA
| | - Yanxiao Han
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, USA
| | - Kyle C. Rouen
- Biophysics Graduate Group, University of California, Davis, Davis, CA, USA
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, USA
| | - Surl-Hee Ahn
- Department of Chemical Engineering, University of California, Davis, Davis, CA, USA
| | - Heike Wulff
- Department of Pharmacology, University of California, Davis, Davis, CA, USA
| | - Colleen E. Clancy
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, USA
- Department of Pharmacology, University of California, Davis, Davis, CA, USA
- Center for Precision Medicine and Data Science, University of California, Davis, Davis, CA, USA
| | - Vladimir Yarov-Yarovoy
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, USA
- Department of Anesthesiology and Pain Medicine, University of California, Davis, Davis, CA, USA
| | - Igor Vorobyov
- Department of Physiology and Membrane Biology, University of California, Davis, Davis, CA, USA
- Department of Pharmacology, University of California, Davis, Davis, CA, USA
| |
Collapse
|
10
|
Agrawal NR, Duan C, Wang R. Nature of Overcharging and Charge Inversion in Electrical Double Layers. J Phys Chem B 2024; 128:303-311. [PMID: 38150660 DOI: 10.1021/acs.jpcb.3c04739] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2023]
Abstract
Understanding overcharging and charge inversion is one of the long-standing challenges in soft matter and biophysics. To study these phenomena, we employ the modified Gaussian renormalized fluctuation theory, which allows for the self-consistent accounting of spatially varying ionic strength as well as the spatial variations in dielectric permittivity and excluded volume effects. The underlying dependence of overcharging on the electrostatic coupling is elucidated by varying the surface charge, counterion valency, and dielectric contrast. Consistent with simulations, three characteristic regimes corresponding to weak, moderate, and strong coupling are identified. Important features like the inversion of zeta potential, crowding, and ionic layering at the surface are successfully captured. For weak coupling, there is no overcharging. In the moderate coupling regime, overcharging increases with the surface charge. Finally, in the strong coupling regime, ionic crowding and saturation in overcharging are observed. Our theory predicts a nonmonotonic dependence of charge inversion on multivalent salt concentration as well as the addition of monovalent salt, in quantitative agreement with experiments.
Collapse
Affiliation(s)
- Nikhil R Agrawal
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720-1462, United States
| | - Chao Duan
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720-1462, United States
| | - Rui Wang
- Department of Chemical and Biomolecular Engineering, University of California, Berkeley, Berkeley, California 94720-1462, United States
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
11
|
Xia J, Gao H, Pan S, Huang T, Zhang L, Sui K, Gao J, Liu X, Jiang L. Light-Augmented Multi-ion Interaction in MXene Membrane for Simultaneous Water Treatment and Osmotic Power Generation. ACS NANO 2023; 17:25269-25278. [PMID: 38071658 DOI: 10.1021/acsnano.3c08487] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/27/2023]
Abstract
The mixing of wastewater and natural water releases abundant osmotic energy. Harvesting this energy could significantly reduce the energy and economic cost of water treatment, leading to sustainable wastewater treatment technology. Yet, such energy harvesting is highly challenging because it requires a material that is highly permeable to nontoxic ions while rejecting toxic ions in wastewater to reach high power density and prevent environmental pollution. In this work, we demonstrate that a light-augmented biomimetic multi-ion interaction in an MXene membrane can simultaneously realize high permeability of Na+ ions for enhanced osmotic power generation and high selectivity to heavy metal ions up to a ratio of 2050 for wastewater treatment. The Na+ permeability is enhanced by the photothermal effect of the MXene membrane. The transport of heavy metal ions, however, is suppressed because, under angstrom-confinement, heavy metal ions are strongly electrostatically repelled by the increased number of permeating Na+ ions. As a result, the membrane can stably generate osmotic power from simulated industrial wastewater, and the power density can be enhanced by 4 times under light illumination of approximate 1 sun intensity. This work highlights the importance of multi-ion interaction for the transport properties of ionic materials, which remains rarely investigated and poorly understood in previous studies.
Collapse
Affiliation(s)
- Jiaxiang Xia
- State Key Laboratory of Bio-Fibers and Eco-textiles, College of Materials Science and Engineering, Shandong Collaborative Innovation Center of Marine Biobased Fibers and Ecological Textiles, Qingdao University, Qingdao 266071, P. R. China
- Qingdao Institute of Bioenergy and Bioprocess Technology Chinese Academy of Sciences, Qingdao 266101, P. R. China
| | - Hongfei Gao
- Qingdao Institute of Bioenergy and Bioprocess Technology Chinese Academy of Sciences, Qingdao 266101, P. R. China
| | - Shangfa Pan
- Qingdao Institute of Bioenergy and Bioprocess Technology Chinese Academy of Sciences, Qingdao 266101, P. R. China
| | - Tao Huang
- State Key Laboratory of Bio-Fibers and Eco-textiles, College of Materials Science and Engineering, Shandong Collaborative Innovation Center of Marine Biobased Fibers and Ecological Textiles, Qingdao University, Qingdao 266071, P. R. China
- Qingdao Institute of Bioenergy and Bioprocess Technology Chinese Academy of Sciences, Qingdao 266101, P. R. China
| | - Li Zhang
- State Key Laboratory of Bio-Fibers and Eco-textiles, College of Materials Science and Engineering, Shandong Collaborative Innovation Center of Marine Biobased Fibers and Ecological Textiles, Qingdao University, Qingdao 266071, P. R. China
- Qingdao Institute of Bioenergy and Bioprocess Technology Chinese Academy of Sciences, Qingdao 266101, P. R. China
| | - Kunyan Sui
- State Key Laboratory of Bio-Fibers and Eco-textiles, College of Materials Science and Engineering, Shandong Collaborative Innovation Center of Marine Biobased Fibers and Ecological Textiles, Qingdao University, Qingdao 266071, P. R. China
| | - Jun Gao
- Qingdao Institute of Bioenergy and Bioprocess Technology Chinese Academy of Sciences, Qingdao 266101, P. R. China
- Shandong Energy Institute, Qingdao 266101, P. R. China
- Qingdao New Energy Shandong Laboratory, Qingdao 266101, P. R. China
| | - Xueli Liu
- State Key Laboratory of Bio-Fibers and Eco-textiles, College of Materials Science and Engineering, Shandong Collaborative Innovation Center of Marine Biobased Fibers and Ecological Textiles, Qingdao University, Qingdao 266071, P. R. China
| | - Lei Jiang
- CAS Key Laboratory of Bio-inspired Materials and Interfacial Science, Technical Institute of Physics and Chemistry Chinese Academy of Sciences, Beijing 100190, P. R. China
| |
Collapse
|
12
|
Papadourakis M, Sinenka H, Matricon P, Hénin J, Brannigan G, Pérez-Benito L, Pande V, van Vlijmen H, de Graaf C, Deflorian F, Tresadern G, Cecchini M, Cournia Z. Alchemical Free Energy Calculations on Membrane-Associated Proteins. J Chem Theory Comput 2023; 19:7437-7458. [PMID: 37902715 PMCID: PMC11017255 DOI: 10.1021/acs.jctc.3c00365] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2023] [Indexed: 10/31/2023]
Abstract
Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.
Collapse
Affiliation(s)
- Michail Papadourakis
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| | - Hryhory Sinenka
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Pierre Matricon
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Jérôme Hénin
- Laboratoire
de Biochimie Théorique UPR 9080, CNRS and Université Paris Cité, 75005 Paris, France
| | - Grace Brannigan
- Center
for Computational and Integrative Biology, Rutgers University−Camden, Camden, New Jersey 08103, United States of America
- Department
of Physics, Rutgers University−Camden, Camden, New Jersey 08102, United States
of America
| | - Laura Pérez-Benito
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Vineet Pande
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Herman van Vlijmen
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Chris de Graaf
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Francesca Deflorian
- Sosei
Heptares, Steinmetz Building,
Granta Park, Great Abington, Cambridge CB21 6DG, United
Kingdom
| | - Gary Tresadern
- CADD,
In Silico Discovery, Janssen Research &
Development, Turnhoutseweg 30, 2340 Beerse, Belgium
| | - Marco Cecchini
- Institut
de Chimie de Strasbourg, UMR7177, CNRS, Université de Strasbourg, F-67083 Strasbourg Cedex, France
| | - Zoe Cournia
- Biomedical
Research Foundation, Academy of Athens, 4 Soranou Ephessiou, 11527 Athens, Greece
| |
Collapse
|
13
|
Wadud MA, Karal MAS, Moniruzzaman M, Rashid MMO. Effects of membrane potentials on the electroporation of giant unilamellar vesicles. PLoS One 2023; 18:e0291496. [PMID: 37699026 PMCID: PMC10497157 DOI: 10.1371/journal.pone.0291496] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Accepted: 08/31/2023] [Indexed: 09/14/2023] Open
Abstract
Living organisms maintain a resting membrane potential, which plays an important role in various biophysical and biological processes. In the context of medical applications, irreversible electroporation (IRE) is a non-thermal and minimally invasive technique that utilizes precisely controlled electric field pulses of micro- to millisecond durations to effectively ablate cancer and tumor cells. Previous studies on IRE-induced rupture of cell-mimetic giant unilamellar vesicles (GUVs) have primarily been conducted in the absence of membrane potentials. In this study, we investigated the electroporation of GUVs, including parameters such as the rate constant of rupture and the probability of rupture, in the presence of various negative membrane potentials. The membranes of GUVs were prepared using lipids and channel forming proteins. As the membrane potential increased from 0 to -90 mV, the rate constant of rupture showed a significant increase from (7.5 ± 1.6)×10-3 to (35.6 ± 5.5)×10-3 s-1. The corresponding probability of rupture also exhibited a notable increase from 0.40 ± 0.05 to 0.68 ± 0.05. To estimate the pore edge tension, the electric tension-dependent logarithm of the rate constant was fitted with the Arrhenius equation for different membrane potentials. The presence of membrane potential did not lead to any significant changes in the pore edge tension. The increase in electroporation is reasonably explained by the decrease in the prepore free energy barrier. The choice of buffer used in GUVs can significantly influence the kinetics of electroporation. This study provides valuable insights that can contribute to the application of electroporation techniques in the biomedical field.
Collapse
Affiliation(s)
- Md. Abdul Wadud
- Department of Physics, Bangladesh University of Engineering and Technology, Dhaka, Bangladesh
| | | | - Md. Moniruzzaman
- Department of Physics, Bangladesh University of Engineering and Technology, Dhaka, Bangladesh
| | - Md. Mamun Or Rashid
- Department of Pharmacy, Noakhali Science and Technology University, Noakhali, Bangladesh
| |
Collapse
|
14
|
Wazawa T, Nagai T. Joule heating involving ion currents through channel proteins. Biophys Physicobiol 2023; 20:e200030. [PMID: 38124793 PMCID: PMC10728626 DOI: 10.2142/biophysico.bppb-v20.0030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Accepted: 06/26/2023] [Indexed: 12/23/2023] Open
Abstract
Ion currents associated with channel proteins in the presence of membrane potential are ubiquitous in cellular and organelle membranes. When an ion current occurs through a channel protein, Joule heating should occur. However, this Joule heating seems to have been largely overlooked in biology. Here we show theoretical investigation of Joule heating involving channel proteins in biological processes. We used electrochemical potential to derive the Joule's law for an ion current through an ion transport protein in the presence of membrane potential, and we suggest that heat production and absorption can occur. Simulation of temperature distribution around a single channel protein with the Joule heating revealed that the temperature increase was as small as <10-3 K, although an ensemble of channel proteins was suggested to exhibit a noticeable temperature increase. Thereby, we theoretically investigated the Joule heating of systems containing ensembles of channel proteins. Nerve is known to undergo rapid heat production followed by heat absorption during the action potential, and our simulation of Joule heating for a squid giant axon combined with the Hodgkin-Huxley model successfully reproduced the feature of the heat. Furthermore, we extended the theory of Joule heating to uncoupling protein 1 (UCP1), a solute carrier family transporter, which is important to the non-shivering thermogenesis in brown adipose tissue mitochondria (BATM). Our calculations showed that the Joule heat involving UCP1 was comparable to the literature calorimetry data of BATM. Joule heating of ion transport proteins is likely to be one of important mechanisms of cellular thermogenesis.
Collapse
Affiliation(s)
| | - Takeharu Nagai
- SANKEN, Osaka University, Ibaraki, Osaka 567-0047, Japan
| |
Collapse
|
15
|
Chen AY, Brooks BR, Damjanovic A. Ion channel selectivity through ion-modulated changes of selectivity filter p Ka values. Proc Natl Acad Sci U S A 2023; 120:e2220343120. [PMID: 37339196 PMCID: PMC10293820 DOI: 10.1073/pnas.2220343120] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Accepted: 04/26/2023] [Indexed: 06/22/2023] Open
Abstract
In bacterial voltage-gated sodium channels, the passage of ions through the pore is controlled by a selectivity filter (SF) composed of four glutamate residues. The mechanism of selectivity has been the subject of intense research, with suggested mechanisms based on steric effects, and ion-triggered conformational change. Here, we propose an alternative mechanism based on ion-triggered shifts in pKa values of SF glutamates. We study the NavMs channel for which the open channel structure is available. Our free-energy calculations based on molecular dynamics simulations suggest that pKa values of the four glutamates are higher in solution of K+ ions than in solution of Na+ ions. Higher pKa in the presence of K+ stems primarily from the higher population of dunked conformations of the protonated Glu sidechain, which exhibit a higher pKa shift. Since pKa values are close to the physiological pH, this results in predominant population of the fully deprotonated state of glutamates in Na+ solution, while protonated states are predominantly populated in K+ solution. Through molecular dynamics simulations we calculate that the deprotonated state is the most conductive, the singly protonated state is less conductive, and the doubly protonated state has significantly reduced conductance. Thus, we propose that a significant component of selectivity is achieved through ion-triggered shifts in the protonation state, which favors more conductive states for Na+ ions and less conductive states for K+ ions. This mechanism also suggests a strong pH dependence of selectivity, which has been experimentally observed in structurally similar NaChBac channels.
Collapse
Affiliation(s)
- Ada Y. Chen
- Department of Physics & Astronomy, Johns Hopkins University, Baltimore, MD21218
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, NIH, Bethesda, MD20892
| | - Bernard R. Brooks
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, NIH, Bethesda, MD20892
| | - Ana Damjanovic
- Laboratory of Computational Biology, National Heart, Lung and Blood Institute, NIH, Bethesda, MD20892
- Department of Biophysics, Johns Hopkins University, Baltimore, MD21218
| |
Collapse
|
16
|
Ngo K, Mateos DL, Han Y, Rouen KC, Ahn SH, Wulff H, Clancy CE, Yarov-Yarovoy V, Vorobyov I. Elucidating Molecular Mechanisms of Protoxin-2 State-specific Binding to the Human Na V1.7 Channel. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.27.530360. [PMID: 36909474 PMCID: PMC10002706 DOI: 10.1101/2023.02.27.530360] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/06/2023]
Abstract
Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-2 (PTx2) has high selectivity for hNaV1.7 and serves as a valuable scaffold to design novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD as well as the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.
Collapse
Affiliation(s)
- Khoa Ngo
- Biophysics Graduate Group, University of California, Davis
- Department of Physiology and Membrane Biology, University of California, Davis
| | - Diego Lopez Mateos
- Biophysics Graduate Group, University of California, Davis
- Department of Physiology and Membrane Biology, University of California, Davis
| | - Yanxiao Han
- Department of Physiology and Membrane Biology, University of California, Davis
| | - Kyle C. Rouen
- Biophysics Graduate Group, University of California, Davis
- Department of Physiology and Membrane Biology, University of California, Davis
| | - Surl-Hee Ahn
- Department of Chemical Engineering, University of California, Davis
| | - Heike Wulff
- Department of Pharmacology, University of California, Davis
| | - Colleen E. Clancy
- Department of Physiology and Membrane Biology, University of California, Davis
- Department of Pharmacology, University of California, Davis
- Center for Precision Medicine and Data Science, University of California, Davis
| | - Vladimir Yarov-Yarovoy
- Department of Physiology and Membrane Biology, University of California, Davis
- Department of Anesthesiology and Pain Medicine, University of California, Davis
| | - Igor Vorobyov
- Department of Physiology and Membrane Biology, University of California, Davis
- Department of Pharmacology, University of California, Davis
| |
Collapse
|
17
|
Abrahamyan A, Eldstrom J, Sahakyan H, Karagulyan N, Mkrtchyan L, Karapetyan T, Sargsyan E, Kneussel M, Nazaryan K, Schwarz JR, Fedida D, Vardanyan V. Mechanism of external K+ sensitivity of KCNQ1 channels. J Gen Physiol 2023; 155:213880. [PMID: 36809486 PMCID: PMC9960071 DOI: 10.1085/jgp.202213205] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2022] [Revised: 12/20/2022] [Accepted: 01/31/2023] [Indexed: 02/23/2023] Open
Abstract
KCNQ1 voltage-gated K+ channels are involved in a wide variety of fundamental physiological processes and exhibit the unique feature of being markedly inhibited by external K+. Despite the potential role of this regulatory mechanism in distinct physiological and pathological processes, its exact underpinnings are not well understood. In this study, using extensive mutagenesis, molecular dynamics simulations, and single-channel recordings, we delineate the molecular mechanism of KCNQ1 modulation by external K+. First, we demonstrate the involvement of the selectivity filter in the external K+ sensitivity of the channel. Then, we show that external K+ binds to the vacant outermost ion coordination site of the selectivity filter inducing a diminution in the unitary conductance of the channel. The larger reduction in the unitary conductance compared to whole-cell currents suggests an additional modulatory effect of external K+ on the channel. Further, we show that the external K+ sensitivity of the heteromeric KCNQ1/KCNE complexes depends on the type of associated KCNE subunits.
Collapse
Affiliation(s)
- Astghik Abrahamyan
- Molecular Neuroscience Group, Institute of Molecular Biology, National Academy of Sciences of the Republic of Armenia , Yerevan, Armenia
| | - Jodene Eldstrom
- Department of Anesthesiology, Pharmacology and Therapeutics, University of British Columbia , Vancouver, BC, Canada
| | - Harutyun Sahakyan
- Laboratory of Computational Modeling of Biological Processes, Institute of Molecular Biology of National Academy of Sciences of the Republic of Armenia , Yerevan, Armenia
| | - Nare Karagulyan
- Molecular Neuroscience Group, Institute of Molecular Biology, National Academy of Sciences of the Republic of Armenia , Yerevan, Armenia
| | - Liana Mkrtchyan
- Molecular Neuroscience Group, Institute of Molecular Biology, National Academy of Sciences of the Republic of Armenia , Yerevan, Armenia
| | - Tatev Karapetyan
- Molecular Neuroscience Group, Institute of Molecular Biology, National Academy of Sciences of the Republic of Armenia , Yerevan, Armenia
| | - Ernest Sargsyan
- Molecular Neuroscience Group, Institute of Molecular Biology, National Academy of Sciences of the Republic of Armenia , Yerevan, Armenia
| | - Matthias Kneussel
- Institute for Molecular Neurogenetics, Center for Molecular Neurobiology Hamburg , Hamburg, Germany
| | - Karen Nazaryan
- Laboratory of Computational Modeling of Biological Processes, Institute of Molecular Biology of National Academy of Sciences of the Republic of Armenia , Yerevan, Armenia
| | - Jürgen R Schwarz
- Institute for Molecular Neurogenetics, Center for Molecular Neurobiology Hamburg , Hamburg, Germany
| | - David Fedida
- Department of Anesthesiology, Pharmacology and Therapeutics, University of British Columbia , Vancouver, BC, Canada
| | - Vitya Vardanyan
- Molecular Neuroscience Group, Institute of Molecular Biology, National Academy of Sciences of the Republic of Armenia , Yerevan, Armenia
| |
Collapse
|
18
|
Alberini G, Alexis Paz S, Corradi B, Abrams CF, Benfenati F, Maragliano L. Molecular Dynamics Simulations of Ion Permeation in Human Voltage-Gated Sodium Channels. J Chem Theory Comput 2023; 19:2953-2972. [PMID: 37116214 DOI: 10.1021/acs.jctc.2c00990] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/30/2023]
Abstract
The recent determination of cryo-EM structures of voltage-gated sodium (Nav) channels has revealed many details of these proteins. However, knowledge of ionic permeation through the Nav pore remains limited. In this work, we performed atomistic molecular dynamics (MD) simulations to study the structural features of various neuronal Nav channels based on homology modeling of the cryo-EM structure of the human Nav1.4 channel and, in addition, on the recently resolved configuration for Nav1.2. In particular, single Na+ permeation events during standard MD runs suggest that the ion resides in the inner part of the Nav selectivity filter (SF). On-the-fly free energy parametrization (OTFP) temperature-accelerated molecular dynamics (TAMD) was also used to calculate two-dimensional free energy surfaces (FESs) related to single/double Na+ translocation through the SF of the homology-based Nav1.2 model and the cryo-EM Nav1.2 structure, with different realizations of the DEKA filter domain. These additional simulations revealed distinct mechanisms for single and double Na+ permeation through the wild-type SF, which has a charged lysine in the DEKA ring. Moreover, the configurations of the ions in the SF corresponding to the metastable states of the FESs are specific for each SF motif. Overall, the description of these mechanisms gives us new insights into ion conduction in human Nav cryo-EM-based and cryo-EM configurations that could advance understanding of these systems and how they differ from potassium and bacterial Nav channels.
Collapse
Affiliation(s)
- Giulio Alberini
- 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
| | - Sergio Alexis Paz
- Departamento de Química Teórica y Computacional, Facultad de Ciencias Químicas, Universidad Nacional de Córdoba, X5000HUA Córdoba, Argentina
- Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Instituto de Fisicoquímica de Córdoba (INFIQC), X5000HUA Córdoba, Argentina
| | - Beatrice Corradi
- 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
| | - Cameron F Abrams
- Department of Chemical and Biological Engineering, Drexel University, Philadelphia, Pennsylvania 19104, United States
| | - 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
|
19
|
Wang Y. Multidisciplinary Advances Address the Challenges in Developing Drugs against Transient Receptor Potential Channels to Treat Metabolic Disorders. ChemMedChem 2023; 18:e202200562. [PMID: 36530131 DOI: 10.1002/cmdc.202200562] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2022] [Revised: 12/01/2022] [Accepted: 12/16/2022] [Indexed: 12/23/2022]
Abstract
Transient receptor potential (TRP) channels are cation channels that regulate key physiological and pathological processes in response to a broad range of stimuli. Moreover, they systemically regulate the release of hormones, metabolic homeostasis, and complications of diabetes, which positions them as promising therapeutic targets to combat metabolic disorders. Nevertheless, there are significant challenges in the design of TRP ligands with high potency and durability. Herein we summarize the four challenges as hydrophobicity, selectivity, mono-target therapy, and interspecies discrepancy. We present 1134 TRP ligands with diversified modes of TRP-ligand interaction and provide a detailed discussion of the latest strategies, especially cryogenic electron microscopy (cryo-EM) and computational methods. We propose solutions to address the challenges with a critical analysis of advances in membrane partitioning, polypharmacology, biased agonism, and biochemical screening of transcriptional modulators. They are fueled by the breakthrough from cryo-EM, chemoinformatics and bioinformatics. The discussion is aimed to shed new light on designing next-generation drugs to treat obesity, diabetes and its complications, with optimal hydrophobicity, higher mode selectivity, multi-targeting and consistent activities between human and rodents.
Collapse
Affiliation(s)
- Yibing Wang
- School of Kinesiology, Shanghai University of Sport, Shanghai, 200438, P. R. China.,Shanghai Frontiers Science Research Base of Exercise and Metabolic Health, Shanghai, 200438, P. R. China
| |
Collapse
|
20
|
Schackert F, Biedermann J, Abdolvand S, Minniberger S, Song C, Plested AJR, Carloni P, Sun H. Mechanism of Calcium Permeation in a Glutamate Receptor Ion Channel. J Chem Inf Model 2023; 63:1293-1300. [PMID: 36758214 PMCID: PMC9976283 DOI: 10.1021/acs.jcim.2c01494] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/11/2023]
Abstract
The α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptors (AMPARs) are neurotransmitter-activated cation channels ubiquitously expressed in vertebrate brains. The regulation of calcium flux through the channel pore by RNA-editing is linked to synaptic plasticity while excessive calcium influx poses a risk for neurodegeneration. Unfortunately, the molecular mechanisms underlying this key process are mostly unknown. Here, we investigated calcium conduction in calcium-permeable AMPAR using Molecular Dynamics (MD) simulations with recently introduced multisite force-field parameters for Ca2+. Our calculations are consistent with experiment and explain the distinct calcium permeability in different RNA-edited forms of GluA2. For one of the identified metal binding sites, multiscale Quantum Mechanics/Molecular Mechanics (QM/MM) simulations further validated the results from MD and revealed small but reproducible charge transfer between the metal ion and its first solvation shell. In addition, the ion occupancy derived from MD simulations independently reproduced the Ca2+ binding profile in an X-ray structure of an NaK channel mimicking the AMPAR selectivity filter. This integrated study comprising X-ray crystallography, multisite MD, and multiscale QM/MM simulations provides unprecedented insights into Ca2+ permeation mechanisms in AMPARs, and paves the way for studying other biological processes in which Ca2+ plays a pivotal role.
Collapse
Affiliation(s)
- Florian
Karl Schackert
- Computational
Biomedicine (IAS-5/INM-9), Forschungszentrum
Jülich GmbH, 52428 Jülich, Germany,Department
of Physics, RWTH Aachen University, 52062 Aachen, Germany
| | - Johann Biedermann
- Institute
of Biology, Cellular Biophysics, Humboldt
Universität zu Berlin, 10115 Berlin, Germany,Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany
| | - Saeid Abdolvand
- Institute
of Biology, Cellular Biophysics, Humboldt
Universität zu Berlin, 10115 Berlin, Germany,Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany
| | - Sonja Minniberger
- Institute
of Biology, Cellular Biophysics, Humboldt
Universität zu Berlin, 10115 Berlin, Germany,Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany
| | - Chen Song
- Center
for Quantitative Biology, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, 100871, China,Peking-Tsinghua
Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, 100871, China
| | - Andrew J. R. Plested
- Institute
of Biology, Cellular Biophysics, Humboldt
Universität zu Berlin, 10115 Berlin, Germany,Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany
| | - Paolo Carloni
- Computational
Biomedicine (IAS-5/INM-9), Forschungszentrum
Jülich GmbH, 52428 Jülich, Germany,Department
of Physics, RWTH Aachen University, 52062 Aachen, Germany,
| | - Han Sun
- Leibniz
Forschungsinstitut für Molekulare Pharmakologie, 13125 Berlin, Germany,Institute
of Chemistry, TU Berlin, Straße des 17 Juni 135, 10623 Berlin, Germany,
| |
Collapse
|
21
|
Mateos DL, Yarov-Yarovoy V. Structural modeling of peptide toxin-ion channel interactions using RosettaDock. Proteins 2023. [PMID: 36729043 DOI: 10.1002/prot.26474] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2022] [Revised: 12/09/2022] [Accepted: 01/30/2023] [Indexed: 02/03/2023]
Abstract
Voltage-gated ion channels play essential physiological roles in action potential generation and propagation. Peptidic toxins from animal venoms target ion channels and provide useful scaffolds for the rational design of novel channel modulators with enhanced potency and subtype selectivity. Despite recent progress in obtaining experimental structures of peptide toxin-ion channel complexes, structural determination of peptide toxins bound to ion channels in physiologically important states remains challenging. Here we describe an application of RosettaDock approach to the structural modeling of peptide toxins interactions with ion channels. We tested this approach on 10 structures of peptide toxin-ion channel complexes and demonstrated that it can sample near-native structures in all tested cases. Our approach will be useful for improving the understanding of the molecular mechanism of natural peptide toxin modulation of ion channel gating and for the structural modeling of novel peptide-based ion channel modulators.
Collapse
Affiliation(s)
- Diego Lopez Mateos
- Department of Physiology and Membrane Biology, University of California Davis, Davis, California, USA.,Biophysics Graduate Group, University of California Davis, Davis, California, USA
| | - Vladimir Yarov-Yarovoy
- Department of Physiology and Membrane Biology, University of California Davis, Davis, California, USA.,Biophysics Graduate Group, University of California Davis, Davis, California, USA.,Department of Anesthesiology and Pain Medicine, University of California Davis, Davis, California, USA
| |
Collapse
|
22
|
He L, Zhang T, Zhu C, Yan T, Liu J. Crown Ether-Based Ion Transporters in Bilayer Membranes. Chemistry 2023; 29:e202300044. [PMID: 36723493 DOI: 10.1002/chem.202300044] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/07/2023] [Revised: 01/28/2023] [Accepted: 01/30/2023] [Indexed: 02/02/2023]
Abstract
Bilayer membranes that enhance the stability of the cell are essential for cell survival, separating and protecting the interior of the cell from its external environment. Membrane-based channel proteins are crucial for sustaining cellular activities. However, dysfunction of these proteins would induce serial channelopathies, which could be substituted by artificial ion channel analogs. Crown ethers (CEs) are widely studied in the area of artificial ion channels owing to their intrinsic host-guest interaction with different kinds of organic and inorganic ions. Other advantages such as lower price, chemical stability, and easier modification also make CE a research hotspot in the field of synthetic transmembrane nanopores. And numerous CEs-based membrane-active synthetic ion channels were designed and fabricated in the past decades. Herein, the recent progress of CEs-based synthetic ion transporters has been comprehensively summarized in this review, including their design principles, functional mechanisms, controllable properties, and biomedical applications. Furthermore, this review has been concluded by discussing the future opportunities and challenges facing this research field. It is anticipated that this review could offer some inspiration for the future fabrication of novel CEs-derived ion transporters with more advanced structures, properties, and practical applications.
Collapse
Affiliation(s)
- Lei He
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology, Ministry of Education, Key Laboratory of Organosilicon Material Technology of Zhejiang Province, Hangzhou Normal University, 311121, Hangzhou, P. R. China
| | - Tianlong Zhang
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology, Ministry of Education, Key Laboratory of Organosilicon Material Technology of Zhejiang Province, Hangzhou Normal University, 311121, Hangzhou, P. R. China
| | - Canhong Zhu
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology, Ministry of Education, Key Laboratory of Organosilicon Material Technology of Zhejiang Province, Hangzhou Normal University, 311121, Hangzhou, P. R. China
| | - Tengfei Yan
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology, Ministry of Education, Key Laboratory of Organosilicon Material Technology of Zhejiang Province, Hangzhou Normal University, 311121, Hangzhou, P. R. China
| | - Junqiu Liu
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology, Ministry of Education, Key Laboratory of Organosilicon Material Technology of Zhejiang Province, Hangzhou Normal University, 311121, Hangzhou, P. R. China
| |
Collapse
|
23
|
Jeong KJ, Jeong S, Lee S, Son CY. Predictive Molecular Models for Charged Materials Systems: From Energy Materials to Biomacromolecules. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2023; 35:e2204272. [PMID: 36373701 DOI: 10.1002/adma.202204272] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Revised: 07/05/2022] [Indexed: 06/16/2023]
Abstract
Electrostatic interactions play a dominant role in charged materials systems. Understanding the complex correlation between macroscopic properties with microscopic structures is of critical importance to develop rational design strategies for advanced materials. But the complexity of this challenging task is augmented by interfaces present in the charged materials systems, such as electrode-electrolyte interfaces or biological membranes. Over the last decades, predictive molecular simulations that are founded in fundamental physics and optimized for charged interfacial systems have proven their value in providing molecular understanding of physicochemical properties and functional mechanisms for diverse materials. Novel design strategies utilizing predictive models have been suggested as promising route for the rational design of materials with tailored properties. Here, an overview of recent advances in the understanding of charged interfacial systems aided by predictive molecular simulations is presented. Focusing on three types of charged interfaces found in energy materials and biomacromolecules, how the molecular models characterize ion structure, charge transport, morphology relation to the environment, and the thermodynamics/kinetics of molecular binding at the interfaces is discussed. The critical analysis brings two prominent field of energy materials and biological science under common perspective, to stimulate crossover in both research field that have been largely separated.
Collapse
Affiliation(s)
- Kyeong-Jun Jeong
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| | - Seungwon Jeong
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| | - Sangmin Lee
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| | - Chang Yun Son
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), Pohang, 790-784, South Korea
| |
Collapse
|
24
|
Kurki M, Poso A, Bartos P, Miettinen MS. Structure of POPC Lipid Bilayers in OPLS3e Force Field. J Chem Inf Model 2022; 62:6462-6474. [PMID: 36044537 PMCID: PMC9795559 DOI: 10.1021/acs.jcim.2c00395] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
It is crucial for molecular dynamics simulations of biomembranes that the force field parameters give a realistic model of the membrane behavior. In this study, we examined the OPLS3e force field for the carbon-hydrogen order parameters SCH of POPC (1-palmitoyl-2-oleoylphosphatidylcholine) lipid bilayers at varying hydration conditions and ion concentrations. The results show that OPLS3e behaves similarly to the CHARMM36 force field and relatively accurately follows the experimentally measured SCH for the lipid headgroup, the glycerol backbone, and the acyl tails. Thus, OPLS3e is a good choice for POPC bilayer simulations under many biologically relevant conditions. The exception are systems with an abundancy of ions, as similarly to most other force fields OPLS3e strongly overestimates the membrane-binding of cations, especially Ca2+. This leads to undesirable positive charge of the membrane surface and drastically lowers the concentration of Ca2+ in the surrounding solvent, which might cause issues in systems sensitive to correct charge distribution profiles across the membrane.
Collapse
Affiliation(s)
- Milla Kurki
- School
of Pharmacy, University of Eastern Finland, Kuopio Campus, Yliopistonranta 1
C, P.O. Box 1627, 70211 Kuopio, Finland
| | - Antti Poso
- School
of Pharmacy, University of Eastern Finland, Kuopio Campus, Yliopistonranta 1
C, P.O. Box 1627, 70211 Kuopio, Finland
| | - Piia Bartos
- School
of Pharmacy, University of Eastern Finland, Kuopio Campus, Yliopistonranta 1
C, P.O. Box 1627, 70211 Kuopio, Finland,
| | - Markus S. Miettinen
- Department
of Chemistry, University of Bergen, 5007 Bergen, Norway,Computational
Biology Unit, Department of Informatics, University of Bergen, 5007 Bergen, Norway
| |
Collapse
|
25
|
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
|
26
|
Effects of thermostats/barostats on physical properties of liquids by molecular dynamics simulations. J Mol Liq 2022. [DOI: 10.1016/j.molliq.2022.120116] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
27
|
Cao Y, Yang R, Wang W, Jiang S, Yang C, Wang Q, Liu N, Xue Y, Lee I, Meng X, Yuan Z. Effects of membrane lipids on phospholamban pentameric channel structure and ion transportation mechanisms. Int J Biol Macromol 2022; 224:766-775. [DOI: 10.1016/j.ijbiomac.2022.10.164] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2022] [Revised: 10/10/2022] [Accepted: 10/19/2022] [Indexed: 11/05/2022]
|
28
|
Zhu Z, Deng Z, Wang Q, Wang Y, Zhang D, Xu R, Guo L, Wen H. Simulation and Machine Learning Methods for Ion-Channel Structure Determination, Mechanistic Studies and Drug Design. Front Pharmacol 2022; 13:939555. [PMID: 35837274 PMCID: PMC9275593 DOI: 10.3389/fphar.2022.939555] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Accepted: 06/07/2022] [Indexed: 11/13/2022] Open
Abstract
Ion channels are expressed in almost all living cells, controlling the in-and-out communications, making them ideal drug targets, especially for central nervous system diseases. However, owing to their dynamic nature and the presence of a membrane environment, ion channels remain difficult targets for the past decades. Recent advancement in cryo-electron microscopy and computational methods has shed light on this issue. An explosion in high-resolution ion channel structures paved way for structure-based rational drug design and the state-of-the-art simulation and machine learning techniques dramatically improved the efficiency and effectiveness of computer-aided drug design. Here we present an overview of how simulation and machine learning-based methods fundamentally changed the ion channel-related drug design at different levels, as well as the emerging trends in the field.
Collapse
Affiliation(s)
- Zhengdan Zhu
- Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
- Beijing Institute of Big Data Research, Beijing, China
| | - Zhenfeng Deng
- DP Technology, Beijing, China
- School of Pharmaceutical Sciences, Peking University, Beijing, China
| | | | | | - Duo Zhang
- Academy for Advanced Interdisciplinary Studies, Peking University, Beijing, China
- DP Technology, Beijing, China
| | - Ruihan Xu
- DP Technology, Beijing, China
- National Engineering Research Center of Visual Technology, Peking University, Beijing, China
| | | | - Han Wen
- DP Technology, Beijing, China
| |
Collapse
|
29
|
Vallée C, Howlin BJ, Lewis R. Single ion free energy calculation in ASIC1: the importance of the HG loop. Phys Chem Chem Phys 2022; 24:13824-13830. [PMID: 35616206 DOI: 10.1039/d2cp01563c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Acid Sensing Ion Channels (ASICs) are one of the most studied channels of the Epithelial Sodium Channel/Degenerin (ENaC/DEG) superfamily. They are responsible for excitatory responses following acidification of the extracellular medium and are involved in several important physiological roles. The ASIC1 subunit can form a functional homotrimeric channel and its structure is currently the most characterised of the whole ENaC/DEG family. Here we computed the free energy profiles for single ion permeation in two different structures of ASIC1 using both Na+ and Cl- as permeating ions. The first structure is the open structure of the channel from the PDB entry 4NTW, and the second structure is the closed structure with the re-entrant loop which contains the highly conserved 'HG' motif form PDB entry 6VTK. Both structures show cation selective free energy profiles, however the profiles of the permeating Na+ differ significantly between the two structures. Indeed, whereas there is only a small energetically favorable (-0.5 kcal mol-1) location for Na+ in the open channel (4NTW) near the end of the pore, we observed a clear ion binding site (-7.8 kcal mol-1) located in between the 'GAS' belt and the 'HG' loop for the channel containing the re-entrant loop (6VTK). Knowing that the 'GAS' motif was determined as the selectivity filter, our results support previous observations while addressing the importance of the 'HG' motif for the interactions between the pore and the permeating cations.
Collapse
Affiliation(s)
- Cédric Vallée
- Leverhulme Quantum Biology Doctoral Training Centre, University of Surrey, Guildford GU2 5XH, UK. .,Department of Chemistry, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford GU2 7XH, UK.,School of Veterinary Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford GU2 7AL, UK
| | - Brendan J Howlin
- Leverhulme Quantum Biology Doctoral Training Centre, University of Surrey, Guildford GU2 5XH, UK. .,Department of Chemistry, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford GU2 7XH, UK
| | - Rebecca Lewis
- Leverhulme Quantum Biology Doctoral Training Centre, University of Surrey, Guildford GU2 5XH, UK. .,School of Veterinary Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford GU2 7AL, UK
| |
Collapse
|
30
|
Wong SL, Awatade NT, Astore MA, Allan KM, Carnell MJ, Slapetova I, Chen PC, Setiadi J, Pandzic E, Fawcett LK, Widger JR, Whan RM, Griffith R, Ooi CY, Kuyucak S, Jaffe A, Waters SA. Molecular Dynamics and Theratyping in Airway and Gut Organoids Reveal R352Q-CFTR Conductance Defect. Am J Respir Cell Mol Biol 2022; 67:99-111. [PMID: 35471184 PMCID: PMC9273222 DOI: 10.1165/rcmb.2021-0337oc] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
A significant challenge to making targeted cystic fibrosis transmembrane conductance regulator (CFTR) modulator therapies accessible to all individuals with cystic fibrosis (CF) are many mutations in the CFTR gene that can cause CF, most of which remain uncharacterized. Here, we characterized the structural and functional defects of the rare CFTR mutation R352Q, with a potential role contributing to intrapore chloride ion permeation, in patient-derived cell models of the airway and gut. CFTR function in differentiated nasal epithelial cultures and matched intestinal organoids was assessed using an ion transport assay and forskolin-induced swelling assay, respectively. CFTR potentiators (VX-770, GLPG1837, and VX-445) and correctors (VX-809, VX-445, with or without VX-661) were tested. Data from R352Q-CFTR were compared with data of 20 participants with mutations with known impact on CFTR function. R352Q-CFTR has residual CFTR function that was restored to functional CFTR activity by CFTR potentiators but not the corrector. Molecular dynamics simulations of R352Q-CFTR were carried out, which indicated the presence of a chloride conductance defect, with little evidence supporting a gating defect. The combination approach of in vitro patient-derived cell models and in silico molecular dynamics simulations to characterize rare CFTR mutations can improve the specificity and sensitivity of modulator response predictions and aid in their translational use for CF precision medicine.
Collapse
Affiliation(s)
- Sharon L Wong
- University of New South Wales, 7800, School of Women's and Children's Health, Faculty of Medicine, Sydney, New South Wales, Australia.,University of New South Wales, 7800, Molecular and Integrative Cystic Fibrosis Research Centre (miCF_RC), Sydney, New South Wales, Australia
| | - Nikhil T Awatade
- University of New South Wales, 7800, School of Women's and Children's Health, Faculty of Medicine, Sydney, New South Wales, Australia.,University of New South Wales, 7800, Molecular and Integrative Cystic Fibrosis Research Centre (miCF_RC), Sydney, New South Wales, Australia
| | - Miro A Astore
- The University of Sydney, 4334, School of Physics, Sydney, New South Wales, Australia
| | - Katelin M Allan
- University of New South Wales, 7800, School of Women's and Children's Health, Faculty of Medicine, Sydney, New South Wales, Australia.,University of New South Wales, 7800, Molecular and Integrative Cystic Fibrosis Research Centre (miCF_RC), Sydney, New South Wales, Australia
| | - Michael J Carnell
- University of New South Wales, 7800, Biomedical Imaging Facility, Mark Wainwright Analytical Centre, Sydney, New South Wales, Australia
| | - Iveta Slapetova
- University of New South Wales, 7800, Biomedical Imaging Facility, Mark Wainwright Analytical Centre, Sydney, New South Wales, Australia
| | - Po-Chia Chen
- The University of Sydney, 4334, School of Physics, Sydney, New South Wales, Australia
| | - Jeffry Setiadi
- The University of Sydney, 4334, School of Physics, Sydney, New South Wales, Australia
| | - Elvis Pandzic
- University of New South Wales, 7800, Biomedical Imaging Facility, Mark Wainwright Analytical Cen, Sydney, New South Wales, Australia
| | - Laura K Fawcett
- University of New South Wales, 7800, School of Women's and Children's Health, Faculty of Medicine, Sydney, New South Wales, Australia.,University of New South Wales, 7800, Molecular and Integrative Cystic Fibrosis Research Centre (miCF_RC), Sydney, New South Wales, Australia.,Sydney Children's Hospital Randwick, 63623, Department of Respiratory Medicine, Randwick, New South Wales, Australia
| | - John R Widger
- University of New South Wales, 7800, School of Women's and Children's Health, Faculty of Medicine, Sydney, New South Wales, Australia.,University of New South Wales, 7800, Molecular and Integrative Cystic Fibrosis Research Centre (miCF_RC), Sydney, New South Wales, Australia.,Sydney Children's Hospital Randwick, 63623, Department of Respiratory Medicine, Randwick, New South Wales, Australia
| | - Renee M Whan
- University of New South Wales, 7800, Biomedical Imaging Facility, Mark Wainwright Analytical Centre, Sydney, New South Wales, Australia
| | - Renate Griffith
- University of New South Wales, 7800, School of Chemistry, Sydney, New South Wales, Australia
| | - Chee Y Ooi
- Sydney Children's Hospital Randwick, Gastroenterology, Sydney, New South Wales, Australia
| | - Serdar Kuyucak
- The University of Sydney, 4334, School of Physics, Sydney, New South Wales, Australia
| | - Adam Jaffe
- Sydney Children`s Hospital, Respiratory Medicine, Sydney, New South Wales, Australia.,University of New South Wales, 7800, School of Women`s and Children`s Health, Sydney, New South Wales, Australia
| | - Shafagh A Waters
- Sydney Children's Hospital, Department of Respiratory Medicine, Sydney, New South Wales, Australia.,Univeristy of New South Wales, School of Women's and Children's Health, Sydney, New South Wales, Australia;
| |
Collapse
|
31
|
Salahub DR. Multiscale molecular modelling: from electronic structure to dynamics of nanosystems and beyond. Phys Chem Chem Phys 2022; 24:9051-9081. [PMID: 35389399 DOI: 10.1039/d1cp05928a] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Important contemporary biological and materials problems often depend on interactions that span orders of magnitude differences in spatial and temporal dimensions. This Tutorial Review attempts to provide an introduction to such fascinating problems through a series of case studies, aimed at beginning researchers, graduate students, postdocs and more senior colleagues who are changing direction to focus on multiscale aspects of their research. The choice of specific examples is highly personal, with examples either chosen from our own work or outstanding multiscale efforts from the literature. I start with various embedding schemes, as exemplified by polarizable continuum models, 3-D RISM, molecular DFT and frozen-density embedding. Next, QM/MM (quantum mechanical/molecular mechanical) techniques are the workhorse of pm-to-nm/ps-to-ns simulations; examples are drawn from enzymes and from nanocatalysis for oil-sands upgrading. Using polarizable force-fields in the QM/MM framework represents a burgeoning subfield; with examples from ion channels and electron dynamics in molecules subject to strong external fields, probing the atto-second dynamics of the electrons with RT-TDDFT (real-time - time-dependent density functional theory) eventually coupled with nuclear motion through the Ehrenfest approximation. This is followed by a section on coarse graining, bridging dimensions from atoms to cells. The penultimate chapter gives a quick overview of multiscale approaches that extend into the meso- and macro-scales, building on atomistic and coarse-grained techniques to enter the world of materials engineering, on the one hand, and cell biology, on the other. A final chapter gives just a glimpse of the burgeoning impact of machine learning on the structure-dynamics front. I aim to capture the excitement of contemporary leading-edge breakthroughs in the description of physico-chemical systems and processes in complex environments, with only enough historical content to provide context and aid the next generation of methodological development. While I aim also for a clear description of the essence of methodological breakthroughs, equations are kept to a minimum and detailed formalism and implementation details are left to the references. My approach is very selective (case studies) rather than exhaustive. I think that these case studies should provide fodder to build as complete a reference tree on multiscale modelling as the reader may wish, through forward and backward citation analysis. I hope that my choices of cases will excite interest in newcomers and help to fuel the growth of multiscale modelling in general.
Collapse
Affiliation(s)
- Dennis R Salahub
- Department of Chemistry, Department of Physics and Astronomy, CMS-Centre for Molecular Simulation, IQST-Institute for Quantum Science and Technology, Quantum Alberta, University of Calgary, Calgary, Alberta, T2N 1N4, Canada.
| |
Collapse
|
32
|
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
|
33
|
Ion permeation, selectivity, and electronic polarization in fluoride channels. Biophys J 2022; 121:1336-1347. [PMID: 35151630 PMCID: PMC9034187 DOI: 10.1016/j.bpj.2022.02.019] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2021] [Revised: 02/01/2022] [Accepted: 02/09/2022] [Indexed: 12/16/2022] Open
Abstract
Fluoride channels (Flucs) export toxic F- from the cytoplasm. Crystallography and mutagenesis have identified several conserved residues crucial for fluoride transport, but the permeation mechanism at the molecular level has remained elusive. Herein, we have applied constant-pH molecular dynamics and free-energy-sampling methods to investigate fluoride permeation through a Fluc protein from Escherichia coli. We find that fluoride is facile to permeate in its charged form, i.e., F-, by traversing through a non-bonded network. The extraordinary F- selectivity is gained by the hydrogen-bonding capability of the central binding site and the Coulombic filter at the channel entrance. The F- permeation rate calculated using an electronically polarizable force field is significantly more accurate compared with the experimental value than that calculated using a more standard additive force field, suggesting an essential role for electronic polarization in the F--Fluc interactions.
Collapse
|
34
|
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: 7] [Impact Index Per Article: 3.5] [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
|
35
|
Jäger M, Koslowski T, Wolf S. Predicting Ion Channel Conductance via Dissipation-Corrected Targeted Molecular Dynamics and Langevin Equation Simulations. J Chem Theory Comput 2021; 18:494-502. [PMID: 34928150 DOI: 10.1021/acs.jctc.1c00426] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Ion channels are important proteins for physiological information transfer and functional control. To predict the microscopic origins of their voltage-conductance characteristics, here we applied dissipation-corrected targeted molecular dynamics in combination with Langevin equation simulations to potassium diffusion through the gramicidin A channel as a test system. Performing a nonequilibrium principal component analysis on backbone dihedral angles, we find coupled protein-ion dynamics to occur during ion transfer. The dissipation-corrected free energy profiles correspond well to predictions from other biased simulation methods. The incorporation of an external electric field in Langevin simulations enables the prediction of macroscopic observables in the form of I-V characteristics.
Collapse
Affiliation(s)
- Miriam Jäger
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104 Freiburg, Germany
| | - Thorsten Koslowski
- Institute of Physical Chemistry, University of Freiburg, 79104 Freiburg, Germany
| | - Steffen Wolf
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104 Freiburg, Germany
| |
Collapse
|
36
|
Zhang S, Gao N, Liu K. Insights on the intercalation mechanism of the coal-bearing kaolinite intercalation based on experimental investigation and molecular dynamics simulations. CHEMICAL PAPERS 2021. [DOI: 10.1007/s11696-021-01803-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
37
|
Ding C, Zhao Y, Su J. Electropumping Phenomenon in Modified Carbon Nanotubes. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2021; 37:12318-12326. [PMID: 34644087 DOI: 10.1021/acs.langmuir.1c01793] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Controlling the water transport in a given direction is essential to the design of novel nanofluidic devices, which is still a challenge because of thermal fluctuations on the nanoscale. In this work, we find an interesting electropumping phenomenon for charge-modified carbon nanotubes (CNTs) through a series of molecular dynamics simulations. In electric fields, the flowing counterions on the CNT inner surface provide a direct driving force for water conduction. Specifically, the dynamics of cations and anions exhibit distinct behaviors that lead to thoroughly different water dynamics in positively and negatively charged CNTs. Because of the competition between the increased ion number and ion-CNT interaction, the cation flux displays an interesting maximum behavior with the increase in surface charge density; however, the anion flux rises further at higher charge density because it is less attractive to the surface. Thus, the anion flux is always several times larger than cation flux that induces a higher water flux in positive CNTs with nearly 100% pumping efficiency, which highly exceeds the efficiency of pristine CNTs. With the change in charge density, the translocation time, occupancy number, and radial density profiles for water and ions also demonstrate a nontrivial difference for positive and negative CNTs. Furthermore, the ion flux exhibits an excellent linear relationship with the field strength, leading to the same water flux behavior. For the change in salt concentration, the pumping efficiency for positive CNTs is also nearly 100%. Our results provide significant new insight into the ionic transport through modified CNTs and should be helpful for the design of nanometer water pumps.
Collapse
|
38
|
Vallée C, Howlin B, Lewis R. Ion Selectivity in the ENaC/DEG Family: A Systematic Review with Supporting Analysis. Int J Mol Sci 2021; 22:ijms222010998. [PMID: 34681656 PMCID: PMC8536179 DOI: 10.3390/ijms222010998] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/17/2021] [Revised: 10/01/2021] [Accepted: 10/05/2021] [Indexed: 12/16/2022] Open
Abstract
The Epithelial Sodium Channel/Degenerin (ENaC/DEG) family is a superfamily of sodium-selective channels that play diverse and important physiological roles in a wide variety of animal species. Despite their differences, they share a high homology in the pore region in which the ion discrimination takes place. Although ion selectivity has been studied for decades, the mechanisms underlying this selectivity for trimeric channels, and particularly for the ENaC/DEG family, are still poorly understood. This systematic review follows PRISMA guidelines and aims to determine the main components that govern ion selectivity in the ENaC/DEG family. In total, 27 papers from three online databases were included according to specific exclusion and inclusion criteria. It was found that the G/SxS selectivity filter (glycine/serine, non-conserved residue, serine) and other well conserved residues play a crucial role in ion selectivity. Depending on the ion type, residues with different properties are involved in ion permeability. For lithium against sodium, aromatic residues upstream of the selectivity filter seem to be important, whereas for sodium against potassium, negatively charged residues downstream of the selectivity filter seem to be important. This review provides new perspectives for further studies to unravel the mechanisms of ion selectivity.
Collapse
Affiliation(s)
- Cédric Vallée
- Leverhulme Quantum Biology Doctoral Training Centre, University of Surrey, Guildford GU2 5XH, UK; (C.V.); (B.H.)
- Department of Chemistry, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford GU2 7XH, UK
- School of Veterinary Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford GU2 7AL, UK
| | - Brendan Howlin
- Leverhulme Quantum Biology Doctoral Training Centre, University of Surrey, Guildford GU2 5XH, UK; (C.V.); (B.H.)
- Department of Chemistry, Faculty of Engineering and Physical Sciences, University of Surrey, Guildford GU2 7XH, UK
| | - Rebecca Lewis
- Leverhulme Quantum Biology Doctoral Training Centre, University of Surrey, Guildford GU2 5XH, UK; (C.V.); (B.H.)
- School of Veterinary Medicine, Faculty of Health and Medical Sciences, University of Surrey, Guildford GU2 7AL, UK
- Correspondence:
| |
Collapse
|
39
|
Şterbuleac D. Molecular dynamics: a powerful tool for studying the medicinal chemistry of ion channel modulators. RSC Med Chem 2021; 12:1503-1518. [PMID: 34671734 PMCID: PMC8459385 DOI: 10.1039/d1md00140j] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2021] [Accepted: 07/21/2021] [Indexed: 01/10/2023] Open
Abstract
Molecular dynamics (MD) simulations allow researchers to investigate the behavior of desired biological targets at ever-decreasing costs with ever-increasing precision. Among the biological macromolecules, ion channels are remarkable transmembrane proteins, capable of performing special biological processes and revealing a complex regulatory matrix, including modulation by small molecules, either endogenous or exogenous. Recently, given the developments in ion channel structure determination and accessibility of bio-computational techniques, MD and related tools are becoming increasingly popular in the intense research area regarding ligand-channel interactions. This review synthesizes and presents the most important fields of MD involvement in investigating channel-molecule interactions, including, but not limited to, deciphering the binding modes of ligands to their ion channel targets and the mechanisms through which chemical compounds exert their effect on channel function. Special attention is devoted to the importance of more elaborate methods, such as free energy calculations, while principles regarding drug design and discovery are highlighted. Several technical aspects involving the creation and simulation of channel-molecule MD systems (ligand parameterization, proper membrane setup, system building, etc.) are also presented.
Collapse
Affiliation(s)
- Daniel Şterbuleac
- Department of Health and Human Development, "Ştefan cel Mare" University of Suceava Str. Universităţii 13, 720229, E Building Suceava Romania
- Department of Forestry and Environmental Protection, "Ştefan cel Mare" University of Suceava Str. Universităţii 13, 720229, E Building Suceava Romania
- Integrated Center for Research, Development and Innovation in Advanced Materials, Nanotechnologies and Distributed Systems for Fabrication and Control (MANSiD), "Ştefan cel Mare" University of Suceava Str. Universităţii 13 720229 Suceava Romania
| |
Collapse
|
40
|
Conductance selectivity of Na + across the K + channel via Na + trapped in a tortuous trajectory. Proc Natl Acad Sci U S A 2021; 118:2017168118. [PMID: 33741736 DOI: 10.1073/pnas.2017168118] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Ion selectivity of the potassium channel is crucial for regulating electrical activity in living cells; however, the mechanism underlying the potassium channel selectivity that favors large K+ over small Na+ remains unclear. Generally, Na+ is not completely excluded from permeation through potassium channels. Herein, the distinct nature of Na+ conduction through the prototypical KcsA potassium channel was examined. Single-channel current recordings revealed that, at a high Na+ concentration (200 mM), the channel was blocked by Na+, and this blocking was relieved at high membrane potentials, suggesting the passage of Na+ across the channel. At a 2,000 mM Na+ concentration, single-channel Na+ conductance was measured as one-eightieth of the K+ conductance, indicating that the selectivity filter allows substantial conduit of Na+ Molecular dynamics simulations revealed unprecedented atomic trajectories of Na+ permeation. In the selectivity filter having a series of carbonyl oxygen rings, a smaller Na+ was distributed off-center in eight carbonyl oxygen-coordinated sites as well as on-center in four carbonyl oxygen-coordinated sites. This amphipathic nature of Na+ coordination yielded a continuous but tortuous path along the filter. Trapping of Na+ in many deep free energy wells in the filter caused slow elution. Conversely, K+ is conducted via a straight path, and as the number of occupied K+ ions increased to three, the concerted conduction was accelerated dramatically, generating the conductance selectivity ratio of up to 80. The selectivity filter allows accommodation of different ion species, but the ion coordination and interactions between ions render contrast conduction rates, constituting the potassium channel conductance selectivity.
Collapse
|
41
|
Hempel T, Del Razo MJ, Lee CT, Taylor BC, Amaro RE, Noé F. Independent Markov decomposition: Toward modeling kinetics of biomolecular complexes. Proc Natl Acad Sci U S A 2021; 118:e2105230118. [PMID: 34321356 PMCID: PMC8346863 DOI: 10.1073/pnas.2105230118] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
To advance the mission of in silico cell biology, modeling the interactions of large and complex biological systems becomes increasingly relevant. The combination of molecular dynamics (MD) simulations and Markov state models (MSMs) has enabled the construction of simplified models of molecular kinetics on long timescales. Despite its success, this approach is inherently limited by the size of the molecular system. With increasing size of macromolecular complexes, the number of independent or weakly coupled subsystems increases, and the number of global system states increases exponentially, making the sampling of all distinct global states unfeasible. In this work, we present a technique called independent Markov decomposition (IMD) that leverages weak coupling between subsystems to compute a global kinetic model without requiring the sampling of all combinatorial states of subsystems. We give a theoretical basis for IMD and propose an approach for finding and validating such a decomposition. Using empirical few-state MSMs of ion channel models that are well established in electrophysiology, we demonstrate that IMD models can reproduce experimental conductance measurements with a major reduction in sampling compared with a standard MSM approach. We further show how to find the optimal partition of all-atom protein simulations into weakly coupled subunits.
Collapse
Affiliation(s)
- Tim Hempel
- Department of Mathematics and Computer Science, Freie Universität Berlin, 14195 Berlin, Germany
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Mauricio J Del Razo
- Department of Mathematics and Computer Science, Freie Universität Berlin, 14195 Berlin, Germany
- Van't Hoff Institute for Molecular Sciences, University of Amsterdam, 1090 GD Amsterdam, The Netherlands
- Korteweg-de Vries Institute for Mathematics, University of Amsterdam, 1090 GE Amsterdam, The Netherlands
- Dutch Institute for Emergent Phenomena, 1090 GL Amsterdam, The Netherlands
| | - Christopher T Lee
- Department of Mechanical and Aerospace Engineering, University of California San Diego, La Jolla, CA 92093
| | - Bryn C Taylor
- Biomedical Sciences Graduate Program, University of California San Diego, La Jolla, CA 92093
| | - Rommie E Amaro
- Department of Chemistry & Biochemistry, University of California San Diego, La Jolla, CA 92093;
| | - Frank Noé
- Department of Mathematics and Computer Science, Freie Universität Berlin, 14195 Berlin, Germany;
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
- Department of Chemistry, Rice University, Houston, TX 77005
| |
Collapse
|
42
|
Zhang Y, Jiang Y, Peng J, Zhang H. Rational Design of Nonbonded Point Charge Models for Divalent Metal Cations with Lennard-Jones 12-6 Potential. J Chem Inf Model 2021; 61:4031-4044. [PMID: 34313132 DOI: 10.1021/acs.jcim.1c00580] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Exploring a metal-involved biochemical process at a molecular level often requires a reliable description of metal properties in aqueous solution by classical nonbonded models. An additional C4 term for considering ion-induced dipole interactions was previously proposed to supplement the widely used Lennard-Jones 12-6 potential (known as the 12-6-4 LJ-type model) with good accuracy. Here, we demonstrate an alternative to modeling divalent metal cations (M2+) with the traditional 12-6 LJ potential by developing nonbonded point charge models for use with 11 water models: TIP3P, SPC/E, SPC/Eb, TIP4P-Ew, TIP4P-D, and TIP4P/2005 and the more recent OPC3, TIP3P-FB, OPC, TIP4P-FB, and a99SB-disp. Our designed models simultaneously reproduce the experimental hydration free energy, ion-oxygen distance, and coordination number in the first hydration shell accurately for most of the metal cations, an accuracy equivalent to that of the complex 12-6-4 LJ-type and double exponential potential models. A systematic comparison with the existing M2+ models is presented as well in terms of effective ion radii, diffusion constants, water exchange rates, and ion-water interactions. Molecular dynamics simulations of metal substitution in Escherichia coli glyoxalase I variants show the great potential of our new models for metalloproteins.
Collapse
Affiliation(s)
- Yongguang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Yang Jiang
- Department of Chemistry, Pennsylvania State University, University Park 16802, Pennsylvania, United States
| | - Jiarong Peng
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| | - Haiyang Zhang
- Department of Biological Science and Engineering, School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing 100083, China
| |
Collapse
|
43
|
Jing Z, Rackers JA, Pratt LR, Liu C, Rempe SB, Ren P. Thermodynamics of ion binding and occupancy in potassium channels. Chem Sci 2021; 12:8920-8930. [PMID: 34257893 PMCID: PMC8246295 DOI: 10.1039/d1sc01887f] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2021] [Accepted: 06/01/2021] [Indexed: 12/15/2022] Open
Abstract
Potassium channels modulate various cellular functions through efficient and selective conduction of K+ ions. The mechanism of ion conduction in potassium channels has recently emerged as a topic of debate. Crystal structures of potassium channels show four K+ ions bound to adjacent binding sites in the selectivity filter, while chemical intuition and molecular modeling suggest that the direct ion contacts are unstable. Molecular dynamics (MD) simulations have been instrumental in the study of conduction and gating mechanisms of ion channels. Based on MD simulations, two hypotheses have been proposed, in which the four-ion configuration is an artifact due to either averaged structures or low temperature in crystallographic experiments. The two hypotheses have been supported or challenged by different experiments. Here, MD simulations with polarizable force fields validated by ab initio calculations were used to investigate the ion binding thermodynamics. Contrary to previous beliefs, the four-ion configuration was predicted to be thermodynamically stable after accounting for the complex electrostatic interactions and dielectric screening. Polarization plays a critical role in the thermodynamic stabilities. As a result, the ion conduction likely operates through a simple single-vacancy and water-free mechanism. The simulations explained crystal structures, ion binding experiments and recent controversial mutagenesis experiments. This work provides a clear view of the mechanism underlying the efficient ion conduction and demonstrates the importance of polarization in ion channel simulations. Polarization shapes the energy landscape of ion conduction in potassium channels.![]()
Collapse
Affiliation(s)
- Zhifeng Jing
- Department of Biomedical Engineering, The University of Texas at Austin Austin Texas 78712 USA
| | - Joshua A Rackers
- Center for Integrated Nanotechnologies, Sandia National Laboratories Albuquerque New Mexico 87185 USA
| | - Lawrence R Pratt
- Department of Chemical and Biomolecular Engineering, Tulane University New Orleans Louisiana 70118 USA
| | - Chengwen Liu
- Department of Biomedical Engineering, The University of Texas at Austin Austin Texas 78712 USA
| | - Susan B Rempe
- Center for Integrated Nanotechnologies, Sandia National Laboratories Albuquerque New Mexico 87185 USA
| | - Pengyu Ren
- Department of Biomedical Engineering, The University of Texas at Austin Austin Texas 78712 USA
| |
Collapse
|
44
|
Majumdar J, Moid M, Dasgupta C, Maiti PK. Dielectric Profile and Electromelting of a Monolayer of Water Confined in Graphene Slit Pore. J Phys Chem B 2021; 125:6670-6680. [PMID: 34107687 DOI: 10.1021/acs.jpcb.1c02266] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A monolayer of water confined between two parallel graphene sheets exists in many different phases and exhibits fascinating dielectric properties that have been studied in experiments. In this work, we use molecular dynamics simulations to study how the dielectric properties of a confined monolayer of water is affected by its structure. We consider six of the popular nonpolarizable water models-SPC/E, SPC/Fw, TIP3P, TIP3P_M (modified), TIP4P-2005, and TIP4P-2005f-and find that the in-plane structure of the water molecules at ambient temperature and pressure is strongly dependent on the water model: all the 3-point water models considered here show square ice formation, whereas no such structural ordering is observed for the 4-point water models. This allows us to investigate the role of the in-plane structure of the water monolayer on its dielectric profile. Our simulations show an anomalous perpendicular dielectric constant compared to the bulk, and the models that do not exhibit ice formation show very different dielectric response along the channel width compared to models that exhibit square ice formation. We also demonstrate the occurrence of electromelting of the in-plane ordered water under the application of a perpendicular electric field and find that the critical field for electromelting strongly depends on the water model. Together, we have shown the dependence of confined water properties on the different water structures that it may take when sandwiched between bilayer graphene. These remarkable properties of confined water can be exploited in various nanofluidic devices, artificial ion channels, and molecular sieving.
Collapse
Affiliation(s)
- Jeet Majumdar
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Mohd Moid
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| | - Chandan Dasgupta
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India.,International Centre for Theoretical Sciences, Bangalore 560089, India
| | - Prabal K Maiti
- Centre for Condensed Matter Theory, Department of Physics, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
45
|
Lynagh T, Flood E, Boiteux C, Sheikh ZP, Allen TW, Pless SA. Determinants of ion selectivity in ASIC1a- and ASIC2a-containing acid-sensing ion channels. J Gen Physiol 2021; 152:133617. [PMID: 31952079 PMCID: PMC7062507 DOI: 10.1085/jgp.201812297] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2019] [Accepted: 12/09/2019] [Indexed: 01/10/2023] Open
Abstract
Trimeric acid-sensing ion channels (ASICs) contribute to neuronal signaling by converting extracellular acidification into excitatory sodium currents. Previous work with homomeric ASIC1a implicates conserved leucine (L7') and consecutive glycine-alanine-serine (GAS belt) residues near the middle, and conserved negatively charged (E18') residues at the bottom of the pore in ion permeation and/or selectivity. However, a conserved mechanism of ion selectivity throughout the ASIC family has not been established. We therefore explored the molecular determinants of ion selectivity in heteromeric ASIC1a/ASIC2a and homomeric ASIC2a channels using site-directed mutagenesis, electrophysiology, and molecular dynamics free energy simulations. Similar to ASIC1a, E18' residues create an energetic preference for sodium ions at the lower end of the pore in ASIC2a-containing channels. However, and in contrast to ASIC1a homomers, ion permeation through ASIC2a-containing channels is not determined by L7' side chains in the upper part of the channel. This may be, in part, due to ASIC2a-specific negatively charged residues (E59 and E62) that lower the energy of ions in the upper pore, thus making the GAS belt more important for selectivity. This is confirmed by experiments showing that the L7'A mutation has no effect in ASIC2a, in contrast to ASIC1a, where it eliminated selectivity. ASIC2a triple mutants eliminating both L7' and upper charges did not lead to large changes in selectivity, suggesting a different role for L7' in ASIC2a compared with ASIC1a channels. In contrast, we observed measurable changes in ion selectivity in ASIC2a-containing channels with GAS belt mutations. Our results suggest that ion conduction and selectivity in the upper part of the ASIC pore may differ between subtypes, whereas the essential role of E18' in ion selectivity is conserved. Furthermore, we demonstrate that heteromeric channels containing mutations in only one of two ASIC subtypes provide a means of functionally testing mutations that render homomeric channels nonfunctional.
Collapse
Affiliation(s)
- Timothy Lynagh
- Center for Biopharmaceuticals, Department of Drug Design and Pharmacology, University of Copenhagen, Copenhagen, Denmark
| | - Emelie Flood
- School of Science, RMIT University, Melbourne, Australia
| | - Céline Boiteux
- School of Science, RMIT University, Melbourne, Australia
| | - Zeshan Pervez Sheikh
- Center for Biopharmaceuticals, Department of Drug Design and Pharmacology, University of Copenhagen, Copenhagen, Denmark
| | - Toby W Allen
- School of Science, RMIT University, Melbourne, Australia
| | - Stephan A Pless
- Center for Biopharmaceuticals, Department of Drug Design and Pharmacology, University of Copenhagen, Copenhagen, Denmark
| |
Collapse
|
46
|
Zequn Z, Jiangfang L. Molecular Insights Into the Gating Kinetics of the Cardiac hERG Channel, Illuminated by Structure and Molecular Dynamics. Front Pharmacol 2021; 12:687007. [PMID: 34168566 PMCID: PMC8217747 DOI: 10.3389/fphar.2021.687007] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2021] [Accepted: 05/20/2021] [Indexed: 11/13/2022] Open
Abstract
The rapidly activating delayed rectifier K+ current generated by the cardiac hERG potassium channel encoded by KCNH2 is the most important reserve current for cardiac repolarization. The unique inward rectification characteristics of the hERG channel depend on the gating regulation, which involves crucial structural domains and key single amino acid residues in the full-length hERG channel. Identifying critical molecules involved in the regulation of gating kinetics for the hERG channel requires high-resolution structures and molecular dynamics simulation models. Based on the latest progress in hERG structure and molecular dynamics simulation research, summarizing the molecules involved in the changes in the channel state helps to elucidate the unique gating characteristics of the channel and the reason for its high affinity to cardiotoxic drugs. In this review, we aim to summarize the significant advances in understanding the voltage gating regulation of the hERG channel based on its structure obtained from cryo-electron microscopy and computer simulations, which reveal the critical roles of several specific structural domains and amino acid residues.
Collapse
Affiliation(s)
- Zheng Zequn
- Department of Cardiovascular, Medical College, Ningbo University, Ningbo, China
| | - Lian Jiangfang
- Department of Cardiovascular, Lihuili Hospital Affiliated to Ningbo University, Ningbo, China
| |
Collapse
|
47
|
Negoro T, Hirata K, Lisy JM, Ishiuchi SI, Fujii M. Potassium and sodium ion complexes with a partial peptide of the selectivity filter in K + channels studied by cold ion trap infrared spectroscopy: the effect of hydration. Phys Chem Chem Phys 2021; 23:12045-12050. [PMID: 34075971 DOI: 10.1039/d1cp00936b] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022]
Abstract
Potassium channels allow K+ to rapidly diffuse, while the selectivity filter (SF) actively blocks Na+. The presence of water in the SF during ion translocation remains under debate due to the experimental and computational challenges in characterizing the interactions between water, ions, and the SF. Our bottom-up approach has been applied to a system composed of a partial peptide of the SF (Ac-tyrosine-NHMe) with a metal ion and a single water molecule to probe these interactions. The IR photodissociation spectra of M+Ac-tyrosine-NHMe(H2O) (M = Na, K) combined with quantum chemical calculations revealed that the water molecule binding sites are ion-dependent. In addition, the ion-peptide distances are elongated significantly for the K+ complex in comparison to the Na+ complex by the addition of a single water molecule. This striking structural difference with the water molecule is discussed in relation to ion selectivity and translocation within the K+ channel.
Collapse
Affiliation(s)
- Takumi Negoro
- Laboratory for Chemistry and Life Science, Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8503, Japan. and School of Life Science and Technology, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, Kanagawa 226-8503, Japan
| | - Keisuke Hirata
- Laboratory for Chemistry and Life Science, Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8503, Japan. and Tokyo Tech World Research Hub Initiative (WRHI), Institute of Innovation Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, 226-8503, Japan. and Department of Chemistry, School of Science, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo, 152-8550, Japan
| | - James M Lisy
- Tokyo Tech World Research Hub Initiative (WRHI), Institute of Innovation Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, 226-8503, Japan. and Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, IL 61801, USA
| | - Shun-Ichi Ishiuchi
- Laboratory for Chemistry and Life Science, Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8503, Japan. and Tokyo Tech World Research Hub Initiative (WRHI), Institute of Innovation Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, 226-8503, Japan. and Department of Chemistry, School of Science, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro-ku, Tokyo, 152-8550, Japan
| | - Masaaki Fujii
- Laboratory for Chemistry and Life Science, Institute of Innovative Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama 226-8503, Japan. and School of Life Science and Technology, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, Kanagawa 226-8503, Japan and Tokyo Tech World Research Hub Initiative (WRHI), Institute of Innovation Research, Tokyo Institute of Technology, 4259 Nagatsuta-cho, Midori-ku, Yokohama, 226-8503, Japan.
| |
Collapse
|
48
|
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
|
49
|
Carnevale V, Delemotte L, Howard RJ. Molecular Dynamics Simulations of Ion Channels. Trends Biochem Sci 2021; 46:621-622. [PMID: 33941431 DOI: 10.1016/j.tibs.2021.04.005] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2021] [Revised: 04/06/2021] [Accepted: 04/06/2021] [Indexed: 10/21/2022]
Affiliation(s)
- Vincenzo Carnevale
- Institute for Computational Molecular Science, Temple University, Philadelphia, PA 19122, USA
| | - Lucie Delemotte
- Department of Applied Physics, KTH Royal Institute of Technology, Science for Life Laboratory, Stockholm, Sweden.
| | - Rebecca J Howard
- Department of Biochemistry and Biophysics and Science for Life Laboratory, Stockholm University, 17165 Solna, Sweden
| |
Collapse
|
50
|
Nemchinova M, Melcr J, Wassenaar TA, Marrink SJ, Guskov A. Asymmetric CorA Gating Mechanism as Observed by Molecular Dynamics Simulations. J Chem Inf Model 2021; 61:2407-2417. [PMID: 33886304 PMCID: PMC8154316 DOI: 10.1021/acs.jcim.1c00261] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
![]()
The CorA family of
proteins plays a housekeeping role in the homeostasis
of divalent metal ions in many bacteria and archaea as well as in
mitochondria of eukaryotes, rendering it an important target to study
the mechanisms of divalent transport and regulation across different
life domains. Despite numerous studies, the mechanistic details of
the channel gating and the transport of the metal ions are still not
entirely understood. Here, we use all-atom and coarse-grained molecular
dynamics simulations combined with in vitro experiments
to investigate the influence of divalent cations on the function of
CorA. Simulations reveal pronounced asymmetric movements of monomers
that enable the rotation of the α7 helix and the cytoplasmic
subdomain with the subsequent formation of new interactions and the
opening of the channel. These computational results are functionally
validated using site-directed mutagenesis of the intracellular cytoplasmic
domain residues and biochemical assays. The obtained results infer
a complex network of interactions altering the structure of CorA to
allow gating. Furthermore, we attempt to reconcile the existing gating
hypotheses for CorA to conclude the mechanism of transport of divalent
cations via these proteins.
Collapse
Affiliation(s)
- Mariia Nemchinova
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Josef Melcr
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Tsjerk A Wassenaar
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Siewert J Marrink
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| | - Albert Guskov
- Groningen Institute for Biomolecular Sciences and Biotechnology, University of Groningen, 9747AG Groningen, The Netherlands
| |
Collapse
|