g_elpot: A Tool for Quantifying Biomolecular Electrostatics from Molecular Dynamics TrajectoriesClick to copy article linkArticle link copied!
- Andrei Y. Kostritskii*Andrei Y. Kostritskii*E-mail: [email protected]Institute of Biological Information Processing (IBI-1), Molekular- und Zellphysiologie, and JARA-HPC, Forschungszentrum Jülich, 52425 Jülich, GermanyDepartment of Physics, RWTH Aachen University, 52062 Aachen, GermanyInstitute of Clinical Pharmacology, RWTH Aachen University, 52062 Aachen, GermanyMore by Andrei Y. Kostritskii
- Claudia AllevaClaudia AllevaInstitute of Biological Information Processing (IBI-1), Molekular- und Zellphysiologie, and JARA-HPC, Forschungszentrum Jülich, 52425 Jülich, GermanyMore by Claudia Alleva
- Saskia CönenSaskia CönenInstitute of Biological Information Processing (IBI-1), Molekular- und Zellphysiologie, and JARA-HPC, Forschungszentrum Jülich, 52425 Jülich, GermanyInstitute of Clinical Pharmacology, RWTH Aachen University, 52062 Aachen, GermanyMore by Saskia Cönen
- Jan-Philipp Machtens*Jan-Philipp Machtens*E-mail: [email protected]. Phone: ++49 2461 614043.Institute of Biological Information Processing (IBI-1), Molekular- und Zellphysiologie, and JARA-HPC, Forschungszentrum Jülich, 52425 Jülich, GermanyInstitute of Clinical Pharmacology, RWTH Aachen University, 52062 Aachen, GermanyMore by Jan-Philipp Machtens
Abstract
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_0007.gif)
Electrostatic forces drive a wide variety of biomolecular processes by defining the energetics of the interaction between biomolecules and charged substances. Molecular dynamics (MD) simulations provide trajectories that contain ensembles of structural configurations sampled by biomolecules and their environment. Although this information can be used for high-resolution characterization of biomolecular electrostatics, it has not yet been possible to calculate electrostatic potentials from MD trajectories in a way allowing for quantitative connection to energetics. Here, we present g_elpot, a GROMACS-based tool that utilizes the smooth particle mesh Ewald method to quantify the electrostatics of biomolecules by calculating potential within water molecules that are explicitly present in biomolecular MD simulations. g_elpot can extract the global distribution of the electrostatic potential from MD trajectories and measure its time course in functionally important regions of a biomolecule. To demonstrate that g_elpot can be used to gain biophysical insights into various biomolecular processes, we applied the tool to MD trajectories of the P2X3 receptor, TMEM16 lipid scramblases, the secondary-active transporter GltPh, and DNA complexed with cationic polymers. Our results indicate that g_elpot is well suited for quantifying electrostatics in biomolecular systems to provide a deeper understanding of its role in biomolecular processes.
This publication is licensed under
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
License Summary*
You are free to share(copy and redistribute) this article in any medium or format within the parameters below:
Creative Commons (CC): This is a Creative Commons license.
Attribution (BY): Credit must be given to the creator.
Non-Commercial (NC): Only non-commercial uses of the work are permitted.
No Derivatives (ND): Derivative works may be created for non-commercial purposes, but sharing is prohibited.
*Disclaimer
This summary highlights only some of the key features and terms of the actual license. It is not a license and has no legal value. Carefully review the actual license before using these materials.
Introduction
Methods
Quantification of Electrostatics Based on Water-Molecule Potentials
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m001.gif)
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m002.gif)
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m003.gif)
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m004.gif)
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m005.gif)
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m006.gif)
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m007.gif)
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m008.gif)
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m009.gif)
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m010.gif)
Implementation
Figure 1
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_0002.gif)
Figure 1. Water-based scheme to calculate the electrostatic potential. First, the long-range part of the potential is calculated by solving the Poisson equation on the grid using the SPME method. Second, each water molecule is assigned a potential averaged over a sphere around its center of geometry, resulting in a distribution of water-molecule potentials. Third, the frame is fitted onto the reference structure and the grid of the final electrostatic-potential map is built around the biomolecule. Fourth, the water-molecule potentials are distributed over the nodes of the electrostatic-potential map, followed by normalization of the node potential based on the number of contributing water molecules. The final map is used for the quantitative analysis of electrostatics, e.g., to generate a potential profile.
Availability and Usage
Simulation Details
P2X3 System
TMEM16 Systems
GltPh Systems
DNA Systems
Calculating the Electrostatic-Potential Profiles
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_m011.gif)
Results and Discussion
Electrostatics within the Cation-Selective Pore of the P2X3 Receptor
Figure 2
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_0003.gif)
Figure 2. Calculation of the electrostatic potential within the pore of the P2X3 cation channel. (a) (Left) Structure of the P2X3 receptor transmembrane domain without the ATP-binding ectodomain and embedded in a lipid membrane. The oxygen and hydrogen atoms of water molecules and the phosphorus atoms around the protein are shown as red, white, and orange spheres, respectively. (Middle) The water-occupancy map, contoured at an occupancy level of 0.1. (Right) Electrostatic-potential map, obtained with β = 20 nm–1 and contoured at −400 mV. The maps were calculated from the first 100 ns block of the 1 μs trajectory. (b) Electrostatic-potential profiles along the P2X3 pore, calculated from the first 100 ns block of the 1 μs trajectory using different values for the β parameter. (c) Average electrostatic-potential profile along the P2X3 pore based on the whole 1 μs trajectory (β = 20 nm–1). Whiskers indicate the standard deviation of the potential, as calculated for 10 separate 100 ns blocks.
Electrostatics within the Proteolipidic Pores of TMEM16 Lipid Scramblases
Figure 3
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_0004.gif)
Figure 3. Electrostatic-potential profiles along the proteolipidic pore of TMEM16 lipid scramblases. (a) Dimeric structure of nhTMEM16 in a lipid membrane. Protomers are colored white and gray; phosphorus atoms are shown in orange. (b,c) Side views of nhTMEM16 in the mixed POPC:POPS membrane (b) and of TMEM16K in the POPC membrane (c), together with the examples of the masked water-occupancy maps used to calculate the electrostatic-potential profiles. The maps were contoured at an occupancy level of 0.1. Phosphorus atoms of POPC and POPS are shown in orange and cyan, respectively. (d) Average electrostatic-potential profiles along the proteolipidic pore of TMEM16K and nhTMEM16 in pure POPC and mixed POPC:POPS membranes. The profiles were averaged over 1 μs trajectories. Whiskers indicate the standard deviation of the potential, as calculated for 10 separate 100 ns blocks.
Electrostatics in the Substrate-Binding Site of the GltPh Transporter
Figure 4
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_0005.gif)
Figure 4. Sodium binding increases the electrostatic potential at the aspartate-binding site of GltPh. (a) Structure of the aspartate-bound GltPh protomer in the outward-facing conformation shown together with different occupancy states of the binding sites. The aspartate molecule is shown in violet, and Na+ ions are shown as blue spheres. The violet cross in the apo panel indicates the position where the potential was measured. Membrane borders are shown schematically. (b) Time courses of electrostatic potential in the aspartate-binding site (violet) and Na+ occupancy of the Na1 site (blue) in one of the protomers upon spontaneous Na+ binding in an MD simulation initiated from the apo state. The time courses are shown as running averages with 50 ns windows. (c) Electrostatic potential in the aspartate-binding site in different states. Each data point is an average potential from a 50 ns block of a single protomer trajectory, whiskers indicate the 5th and 95th percentiles, and a median value is shown for each state.
Local Effects of Polycations on DNA Electrostatics
Figure 5
![](/cms/10.1021/acs.jctc.0c01246/asset/images/medium/ct0c01246_0006.gif)
Figure 5. Local effects of polycations on DNA electrostatics. (a) Structure of DNA in complex with cationic polymers PVA (left) and PAA (right). The backbone and grooves of DNA are shown in red and white, respectively. Polymers are shown as spheres with nitrogen atoms in blue and carbon atoms in teal (PVA) and yellow (PAA); hydrogen atoms are omitted. (b) Radial distribution of electrostatic potential around the central axis of DNA either in free form or in complex with PVA or PAA, calculated either globally (top), or separately for major (middle) and minor (bottom) grooves. Whiskers indicate the standard deviation of the potential, calculated separately for 100 ns blocks.
Conclusions and Perspectives
Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jctc.0c01246.
Dependence of the short-range part of the potential on β; convergence of P2X3 profiles calculated from the raw SPME potential; time convergence of the g_elpot calculations; raw data on the electrostatic potential and Na1 occupancy in the GltPh simulations; chemical structure of PVA and PAA; groove-specific division of space in DNA systems; and simulation times for DNA and GltPh systems (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.
Acknowledgments
The authors thank Dr. Andrey Gurtovenko for providing trajectories of DNA with polycations and Dr. Christoph Fahlke, Dr. Ralf Hausmann, and Piersilvio Longo for critically reading the manuscript. A.Y.K. thanks Diana Kondinskaia for fruitful discussions. This work was supported by the Deutsche Forschungsgemeinschaft (German Research Foundation) to J.-P.M. (MA 7525/1-2, as part of the Research Unit FOR 2518, DynIon, project P4; and MA 7525/2-1, as part of the Research Unit FOR 6046, project P2) and to Dr. Ralf Hausmann (HA 6095/1-2), the Jülich-Aachen Research Alliance Center for Simulation and Data Science (JARA-CSD) School for Simulation and Data Science (SSD), and by a grant from the Interdisciplinary Centre for Clinical Research within the faculty of Medicine at the RWTH Aachen University (IZKF TN1-3/IA532003). The authors gratefully acknowledge the computing time granted through JARA on the supercomputer JURECA at Forschungszentrum Jülich and the supercomputer CLAIX at RWTH Aachen University.
References
This article references 79 other publications.
- 1Honig, B.; Nicholls, A. Classical electrostatics in biology and chemistry. Science 1995, 268, 1144– 1149, DOI: 10.1126/science.7761829Google Scholar1https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXlvVyqsrY%253D&md5=19b5365db27b7d14c861b67fc79b38cfClassical electrostatics in biology and chemistryHonig, Barry; Nicholls, AnthonyScience (Washington, D. C.) (1995), 268 (5214), 1144-9CODEN: SCIEAS; ISSN:0036-8075. (American Association for the Advancement of Science)A review with 75 refs. A major revival in the use of classical electrostatics as an approach to the study of charged and polar mols. in aq. soln. has been made possible through the development of fast numerical and computational methods to solve the Poisson-Boltzmann equation for solute mols. that have complex shapes and charge distributions. Graphical visualization of the calcd. electrostatic potentials generated by proteins and nucleic acids has revealed insights into the role of electrostatic interactions in a wide range of biol. phenomena. Classical electrostatics has also proved to be a successful quant. tool yielding accurate descriptions of elec. potentials, diffusion-limited processes, pH-dependent properties of proteins, ionic strength-dependent phenonema, and the solvation free energies of org. mols.
- 2Li, C.; Jia, Z.; Chakravorty, A.; Pahari, S.; Peng, Y.; Basu, S.; Koirala, M.; Panday, S. K.; Petukh, M.; Li, L.; Alexov, E. DelPhi Suite: New Developments and Review of Functionalities. J. Comput. Chem. 2019, 40, 2502– 2508, DOI: 10.1002/jcc.26006Google Scholar2https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXht1Cntb7J&md5=7b8a34361ca244d25c3e19b2a20bf86eDelPhi Suite: New Developments and Review of FunctionalitiesLi, Chuan; Jia, Zhe; Chakravorty, Arghya; Pahari, Swagata; Peng, Yunhui; Basu, Sankar; Koirala, Mahesh; Panday, Shailesh Kumar; Petukh, Marharyta; Li, Lin; Alexov, EmilJournal of Computational Chemistry (2019), 40 (28), 2502-2508CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)Electrostatic potential, energies, and forces affect virtually any process in mol. biol., however, computing these quantities is a difficult task due to irregularly shaped macromols. and the presence of water. Here, the authors report a new edition of the popular software package DelPhi along with describing its functionalities. The new DelPhi is a C++ object-oriented package supporting various levels of multiprocessing and memory distribution. It is demonstrated that multiprocessing results in significant improvement of computational time. Furthermore, for computations requiring large grid size (large macromol. assemblages), the approach of memory distribution is shown to reduce the requirement of RAM and thus permitting large-scale modeling to be done on Linux clusters with moderate architecture. The new release comes with new features, whose functionalities and applications are described as well.
- 3Baker, N. A.; Sept, D.; Joseph, S.; Holst, M. J.; McCammon, J. A. Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 10037– 10041, DOI: 10.1073/pnas.181342398Google Scholar3https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXmvFWisbc%253D&md5=1b861999ef12c6972e82e8ada0f387cbElectrostatics of nanosystems: application to microtubules and the ribosomeBaker, Nathan A.; Sept, David; Joseph, Simpson; Holst, Michael J.; McCammon, J. AndrewProceedings of the National Academy of Sciences of the United States of America (2001), 98 (18), 10037-10041CODEN: PNASA6; ISSN:0027-8424. (National Academy of Sciences)Evaluation of the electrostatic properties of biomols. has become a std. practice in mol. biophysics. Foremost among the models used to elucidate the electrostatic potential is the Poisson-Boltzmann equation; however, existing methods for solving this equation have limited the scope of accurate electrostatic calcns. to relatively small biomol. systems. Here we present the application of numerical methods to enable the trivially parallel soln. of the Poisson-Boltzmann equation for supramol. structures that are orders of magnitude larger in size. As a demonstration of this methodol., electrostatic potentials have been calcd. for large microtubule and ribosome structures. The results point to the likely role of electrostatics in a variety of activities of these structures.
- 4Callenberg, K. M.; Choudhary, O. P.; de Forest, G. L.; Gohara, D. W.; Baker, N. A.; Grabe, M. APBSmem: a graphical interface for electrostatic calculations at the membrane. PLoS One 2010, 5, e12722 DOI: 10.1371/journal.pone.0012722Google ScholarThere is no corresponding record for this reference.
- 5Dong, F.; Olsen, B.; Baker, N. A. Computational methods for biomolecular electrostatics. Methods Cell Biol. 2008, 84, 843– 870Google Scholar5https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXht1amurw%253D&md5=865ba17d85d1ef7e60d5e222f28bf107Computational methods for biomolecular electrostaticsDong, Feng; Olsen, Brett; Baker, Nathan A.Methods in Cell Biology (2008), 84 (Biophysical Tools for Biologists, Volume One), 843-870CODEN: MCBLAG; ISSN:0091-679X. (Elsevier)A review. An understanding of intermol. interactions is essential for insight into how cells develop, operate, communicate, and control their activities. Such interactions include several components: contributions from linear, angular, and torsional forces in covalent bonds, van der waals forces, as well as electrostatics. Among the various components of mol. interactions, electrostatics are of special importance because of their long range and their influence on polar or charged mols., including water, aq. ions, and amino or nucleic acids, which are some of the primary components of living systems. Electrostatics, therefore, play important roles in detg. the structure, motion, and function of a wide range of biol. mols. This chapter presents a brief overview of electrostatic interactions in cellular systems, with a particular focus on how computational tools can be used to investigate these types of interactions.
- 6Wang, L. Measurements and implications of the membrane dipole potential. Annu. Rev. Biochem. 2012, 81, 615– 635, DOI: 10.1146/annurev-biochem-070110-123033Google Scholar6https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XhtVGlsL7M&md5=7eeb5e3fc8ad92e8b59ef4a039138b9cMeasurements and implications of the membrane dipole potentialWang, LiguoAnnual Review of Biochemistry (2012), 81 (), 615-635CODEN: ARBOAW; ISSN:0066-4154. (Annual Reviews Inc.)A review. There are three kinds of membrane potentials: the surface potentials, resulting from the accumulation of charges at the membrane surfaces; the transmembrane potential, detd. by imbalance of charge in the aq. solns.; and the dipole potential, a membrane-internal potential from the dipolar components of the phospholipids and interface water. The abs. value of the dipole potential has been very difficult to measure, although its value has been estd. to be at 200-1000 mV from ion translocation rates (detd. by the planar lipid bilayer method), the surface potential of lipid monolayers (detd. by the lipid monolayer method), mol.-dynamics calcns., and electron scattering using cryoelectron microscopy (cryo-EM). Spectroscopy methods have also been used to monitor the dipole potential changes on the basis of the obsd. fluorescence changes of voltage-sensitive probes. The dipole potential accounts for the much larger permeability of a bare phospholipid membrane to anions than cations and affects the conformation and function of membrane proteins.
- 7Warshel, A.; Sharma, P. K.; Kato, M.; Parson, W. W. Modeling electrostatic effects in proteins. Biochim. Biophys. Acta, Proteins Proteomics 2006, 1764, 1647– 1676, DOI: 10.1016/j.bbapap.2006.08.007Google Scholar7https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xht1Sjs7%252FO&md5=78eb047d080cc02e209aa71d2ad278beModeling electrostatic effects in proteinsWarshel, Arieh; Sharma, Pankaz K.; Kato, Mitsunori; Parson, William W.Biochimica et Biophysica Acta, Proteins and Proteomics (2006), 1764 (11), 1647-1676CODEN: BBAPBW; ISSN:1570-9639. (Elsevier Ltd.)A review. Electrostatic energies provide what is perhaps the most effective tool for structure-function correlation of biol. mols. This review considers the current state of simulations of electrostatic energies in macromols. as well as the early developments of this field. We focus on the relationship between microscopic and macroscopic models, considering the convergence problems of the microscopic models and the fact that the dielec. consts. in semimacroscopic models depend on the definition and the specific treatment. The advances and the challenges in the field are illustrated considering a wide range of functional properties including pKa's, redox potentials, ion and proton channels, enzyme catalysis, ligand binding and protein stability. We conclude by pointing out that, despite the current problems and the significant misunderstandings in the field, there is an overall progress that should lead eventually to quant. descriptions of electrostatic effects in proteins and thus to quant. descriptions of the function of proteins.
- 8Cisneros, G. A.; Karttunen, M.; Ren, P.; Sagui, C. Classical electrostatics for biomolecular simulations. Chem. Rev. 2014, 114, 779– 814, DOI: 10.1021/cr300461dGoogle Scholar8https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXhtlGnurrF&md5=86c888e948cf518f01f2a93d66e1a7f7Classical Electrostatics for Biomolecular SimulationsCisneros, G. Andres; Karttunen, Mikko; Ren, Pengyu; Sagui, CelesteChemical Reviews (Washington, DC, United States) (2014), 114 (1), 779-814CODEN: CHREAY; ISSN:0009-2665. (American Chemical Society)In this review, the authors primarily treat electrostatic methods for classical biomol. simulations in explicit solvent, with special emphasis on the accuracy of the electrostatic representation. The first part of the article is dedicated to reviewing computational methods for dealing with the long-range nature of the electrostatic interactions, including traditional methods as well as recent extensions and developments. The second part of the article is dedicated to reviewing current methods for representing the mol. electronic charge cloud, that go well beyond the point-charge representation. In addn., the authors discuss various multiscale approaches at the quantum mechanics/mol. mechanics level and at the mol. mechanics/coarse-grain level. The general trends in both cases are toward more accuracy and proper treatment of electrostatics. The authors finish the article by describing other challenging developments and giving the perspective on the current state of the field.
- 9Hollingsworth, S. A.; Dror, R. O. Molecular Dynamics Simulation for All. Neuron 2018, 99, 1129– 1143, DOI: 10.1016/j.neuron.2018.08.011Google Scholar9https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhslGltLzL&md5=35e14022763288ed45c2a605cc900f61Molecular Dynamics Simulation for AllHollingsworth, Scott A.; Dror, Ron O.Neuron (2018), 99 (6), 1129-1143CODEN: NERNET; ISSN:0896-6273. (Cell Press)A review. The impact of mol. dynamics (MD) simulations in mol. biol. and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other biomols. in full at. detail and at very fine temporal resoln. Major improvements in simulation speed, accuracy, and accessibility, together with the proliferation of exptl. structural data, have increased the appeal of biomol. simulation to experimentalists-a trend particularly noticeable in, although certainly not limited to, neuroscience. Simulations have proven valuable in deciphering functional mechanisms of proteins and other biomols., in uncovering the structural basis for disease, and in the design and optimization of small mols., peptides, and proteins. Here we describe, in practical terms, the types of information MD simulations can provide and the ways in which they typically motivate further exptl. work.
- 10Kutzner, C.; Köpfer, D. A.; Machtens, J.-P.; de Groot, B. L.; Song, C.; Zachariae, U. Insights into the function of ion channels by computational electrophysiology simulations. Biochim. Biophys. Acta, Biomembr. 2016, 1858, 1741– 1752, DOI: 10.1016/j.bbamem.2016.02.006Google Scholar10https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XisFOrsb4%253D&md5=20860bb3a64b0140a3957639bf7da039Insights into the function of ion channels by computational electrophysiology simulationsKutzner, Carsten; Koepfer, David A.; Machtens, Jan-Philipp; de Groot, Bert L.; Song, Chen; Zachariae, UlrichBiochimica et Biophysica Acta, Biomembranes (2016), 1858 (7_Part_B), 1741-1752CODEN: BBBMBS; ISSN:0005-2736. (Elsevier B.V.)A review. Ion channels are of universal importance for all cell types and play key roles in cellular physiol. and pathol. Increased insight into their functional mechanisms is crucial to enable drug design on this important class of membrane proteins, and to enhance our understanding of some of the fundamental features of cells. This review presents the concepts behind the recently developed simulation protocol Computational Electrophysiol. (CompEL), which facilitates the atomistic simulation of ion channels in action. In addn., the review provides guidelines for its application in conjunction with the mol. dynamics software package GROMACS. We first lay out the rationale for designing CompEL as a method that models the driving force for ion permeation through channels the way it is established in cells, i.e., by electrochem. ion gradients across the membrane. This is followed by an outline of its implementation and a description of key settings and parameters helpful to users wishing to set up and conduct such simulations. In recent years, key mechanistic and biophys. insights have been obtained by employing the CompEL protocol to address a wide range of questions on ion channels and permeation. We summarize these recent findings on membrane proteins, which span a spectrum from highly ion-selective, narrow channels to wide diffusion pores. Finally we discuss the future potential of CompEL in light of its limitations and strengths. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.
- 11Lin, Y.-L.; Aleksandrov, A.; Simonson, T.; Roux, B. An Overview of Electrostatic Free Energy Computations for Solutions and Proteins. J. Chem. Theory Comput. 2014, 10, 2690– 2709, DOI: 10.1021/ct500195pGoogle Scholar11https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXnslaksLc%253D&md5=6aa50d77dc3326ae518f682336b9f4a8An Overview of Electrostatic Free Energy Computations for Solutions and ProteinsLin, Yen-Lin; Aleksandrov, Alexey; Simonson, Thomas; Roux, BenoitJournal of Chemical Theory and Computation (2014), 10 (7), 2690-2709CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A review. Free energy simulations for electrostatic and charging processes in complex mol. systems encounter specific difficulties owing to the long-range, 1/r Coulomb interaction. To calc. the solvation free energy of a simple ion, it is essential to take into account the polarization of nearby solvent but also the electrostatic potential drop across the liq.-gas boundary, however distant. The latter does not exist in a simulation model based on periodic boundary conditions because there is no phys. boundary to the system. An important consequence is that the ref. value of the electrostatic potential is not an ion in a vacuum. Also, in an infinite system, the electrostatic potential felt by a perturbing charge is conditionally convergent and dependent on the choice of computational conventions. Furthermore, with Ewald lattice summation and tinfoil conducting boundary conditions, the charges experience a spurious shift in the potential that depends on the details of the simulation system such as the vol. fraction occupied by the solvent. All these issues can be handled with established computational protocols, as reviewed here and illustrated for several small ions and three solvated proteins.
- 12Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N·log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089– 10092, DOI: 10.1063/1.464397Google Scholar12https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3sXks1Ohsr0%253D&md5=3c9f230bd01b7b714fd096d4d2e755f6Particle mesh Ewald: an N·log(N) method for Ewald sums in large systemsDarden, Tom; York, Darrin; Pedersen, LeeJournal of Chemical Physics (1993), 98 (12), 10089-92CODEN: JCPSA6; ISSN:0021-9606.An N·log(N) method for evaluating electrostatic energies and forces of large periodic systems is presented. The method is based on interpolation of the reciprocal space Ewald sums and evaluation of the resulting convolution using fast Fourier transforms. Timings and accuracies are presented for three large cryst. ionic systems.
- 13Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577– 8593, DOI: 10.1063/1.470117Google Scholar13https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXptlehtrw%253D&md5=092a679dd3bee08da28df41e302383a7A smooth particle mesh Ewald methodEssmann, Ulrich; Perera, Lalith; Berkowitz, Max L.; Darden, Tom; Lee, Hsing; Pedersen, Lee G.Journal of Chemical Physics (1995), 103 (19), 8577-93CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The previously developed particle mesh Ewald method is reformulated in terms of efficient B-spline interpolation of the structure factors. This reformulation allows a natural extension of the method to potentials of the form 1/rp with p ≥ 1. Furthermore, efficient calcn. of the virial tensor follows. Use of B-splines in the place of Lagrange interpolation leads to analytic gradients as well as a significant improvement in the accuracy. The authors demonstrate that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N). For biomol. systems with many thousands of atoms and this method permits the use of Ewald summation at a computational cost comparable to that of a simple truncation method of 10 Å or less.
- 14Aksimentiev, A.; Schulten, K. Imaging alpha-hemolysin with molecular dynamics: ionic conductance, osmotic permeability, and the electrostatic potential map. Biophys. J. 2005, 88, 3745– 3761, DOI: 10.1529/biophysj.104.058727Google Scholar14https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXksl2jsL4%253D&md5=23548fa7b86cb4b542040ef5906d8d2bImaging α-hemolysin with molecular dynamics: Ionic conductance, osmotic permeability, and the electrostatic potential mapAksimentiev, Aleksij; Schulten, KlausBiophysical Journal (2005), 88 (6), 3745-3761CODEN: BIOJAU; ISSN:0006-3495. (Biophysical Society)α-Hemolysin of Staphylococcus aureus is a self-assembling toxin that forms a water-filled transmembrane channel upon oligomerization in a lipid membrane. Apart from being one of the best-studied toxins of bacterial origin, α-hemolysin is the principal component in several biotechnol. applications, including systems for controlled delivery of small solutes across lipid membranes, stochastic sensors for small solutes, and an alternative to conventional technol. for DNA sequencing. Through large-scale mol. dynamics simulations, we studied the permeability of the α-hemolysin/lipid bilayer complex for water and ions. The studied system, composed of ∼300,000 atoms, included one copy of the protein, a patch of a DPPC lipid bilayer, and a 1 M water soln. of KCl. Monitoring the fluctuations of the pore structure revealed an asym., on av., cross section of the α-hemolysin stem. Applying external electrostatic fields produced a transmembrane ionic current; repeating simulations at several voltage biases yielded a current/voltage curve of α-hemolysin and a set of electrostatic potential maps. The selectivity of α-hemolysin to Cl- was found to depend on the direction and the magnitude of the applied voltage bias. The results of our simulations are in excellent quant. agreement with available exptl. data. Analyzing trajectories of all water mol., we computed the α-hemolysin's osmotic permeability for water as well as its electroosmotic effect, and characterized the permeability of its seven side channels. The side channels were found to connect seven His-144 residues surrounding the stem of the protein to the bulk soln.; the protonation of these residues was obsd. to affect the ion conductance, suggesting the seven His-144 to comprise the pH sensor that gates conductance of the α-hemolysin channel.
- 15Humphrey, W.; Dalke, A.; Schulten, K. VMD: visual molecular dynamics. J. Mol. Graphics 1996, 14, 33– 38, DOI: 10.1016/0263-7855(96)00018-5Google Scholar15https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28Xis12nsrg%253D&md5=1e3094ec3151fb85c5ff05f8505c78d5VDM: visual molecular dynamicsHumphrey, William; Dalke, Andrew; Schulten, KlausJournal of Molecular Graphics (1996), 14 (1), 33-8, plates, 27-28CODEN: JMGRDV; ISSN:0263-7855. (Elsevier)VMD is a mol. graphics program designed for the display and anal. of mol. assemblies, in particular, biopolymers such as proteins and nucleic acids. VMD can simultaneously display any no. of structures using a wide variety of rendering styles and coloring methods. Mols. are displayed as one or more "representations," in which each representation embodies a particular rendering method and coloring scheme for a selected subset of atoms. The atoms displayed in each representation are chosen using an extensive atom selection syntax, which includes Boolean operators and regular expressions. VMD provides a complete graphical user interface for program control, as well as a text interface using the Tcl embeddable parser to allow for complex scripts with variable substitution, control loops, and function calls. Full session logging is supported, which produces a VMD command script for later playback. High-resoln. raster images of displayed mols. may be produced by generating input scripts for use by a no. of photorealistic image-rendering applications. VMD has also been expressly designed with the ability to animate mol. dynamics (MD) simulation trajectories, imported either from files or from a direct connection to a running MD simulation. VMD is the visualization component of MDScope, a set of tools for interactive problem solving in structural biol., which also includes the parallel MD program NAMD, and the MDCOMM software used to connect the visualization and simulation programs, VMD is written in C++, using an object-oriented design; the program, including source code and extensive documentation, is freely available via anonymous ftp and through the World Wide Web.
- 16Enkavi, G.; Tajkhorshid, E. Simulation of spontaneous substrate binding revealing the binding pathway and mechanism and initial conformational response of GlpT. Biochemistry 2010, 49, 1105– 1114, DOI: 10.1021/bi901412aGoogle Scholar16https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXovVaktg%253D%253D&md5=be122d6c0290699b0453d120fc26f267Simulation of Spontaneous Substrate Binding Revealing the Binding Pathway and Mechanism and Initial Conformational Response of GlpTEnkavi, Giray; Tajkhorshid, EmadBiochemistry (2010), 49 (6), 1105-1114CODEN: BICHAW; ISSN:0006-2960. (American Chemical Society)Glycerol 3-phosphate transporter (GlpT) mediates the import of glycerol 3-phosphate (G3P) using the gradient of inorg. phosphate (Pi). To study the process and mechanism of substrate binding and to investigate the protein's initial response, we performed equil. simulations of wild-type GlpT and several of its mutant forms in membranes in the presence of all physiol. relevant substrates (Pi-, Pi2-, G3P-, and G3P2-). The simulations capture spontaneous substrate binding of GlpT, driven by the pos. electrostatic potential of the lumen. K80 is found to act as a "hook" making the first encounter with the substrate and guiding it toward the binding site, where it binds tightly to R45, a key binding site residue that acts as a "fork" holding the substrate. R269 establishes no direct contact with the substrate during the simulations, a surprising behavior given its structural pseudosymmetry to R45. In all substrate-bound systems, partial closing of the cytoplasmic half of GlpT was obsd. The substrate appears to stabilize the partially occluded state, as in the two apo simulations either no closing was obsd. or the protein reverted to its open form toward the end of the simulation, whereas in all substrate-bound systems, a stable partially closed state was produced. Along with the modulation of the periplasmic salt bridge network, these substrate-induced events destabilize the periplasmic half while inducing a closure in the cytoplasmic half, thus capturing the early stages of the proposed rocker-switch mechanism in GlpT.
- 17Arcario, M. J.; Tajkhorshid, E. Membrane-induced structural rearrangement and identification of a novel membrane anchor in talin F2F3. Biophys. J. 2014, 107, 2059– 2069, DOI: 10.1016/j.bpj.2014.09.022Google Scholar17https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhslGrsLzO&md5=c6de1b54606c2757ec9a32d21ca1e489Membrane-Induced Structural Rearrangement and Identification of a Novel Membrane Anchor in Talin F2F3Arcario, Mark J.; Tajkhorshid, EmadBiophysical Journal (2014), 107 (9), 2059-2069CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)Exptl. challenges assocd. with characterization of the membrane-bound form of talin have prevented us from understanding the mol. mechanism of its membrane-dependent integrin activation. Here, utilizing what we believe to be a novel membrane mimetic model, we present a reproducible model of membrane-bound talin obsd. across multiple independent simulations. We characterize both local and global membrane-induced structural transitions that successfully reconcile discrepancies between biochem. and structural studies and provide insight into how talin might modulate integrin function. Membrane binding of talin, captured in unbiased simulations, proceeds through three distinct steps: initial electrostatic recruitment of the F2 subdomain to anionic lipids via several basic residues; insertion of an initially buried, conserved hydrophobic anchor into the membrane; and assocn. of the F3 subdomain with the membrane surface through a large, interdomain conformational change. These latter two steps, to our knowledge, have not been obsd. or described previously. Electrostatic anal. shows talin F2F3 to be highly polarized, with a highly pos. underside, which we attribute to the initial electrostatic recruitment, and a neg. top face, which can help orient the protein optimally with respect to the membrane, thereby reducing the no. of unproductive membrane collision events.
- 18Frigo, M.; Johnson, S. G. The Design and Implementation of FFTW3. Proc. IEEE 2005, 93, 216– 231, DOI: 10.1109/JPROC.2004.840301Google ScholarThere is no corresponding record for this reference.
- 19Phillips, J. C.; Hardy, D. J.; Maia, J. D. C.; Stone, J. E.; Ribeiro, J. V.; Bernardi, R. C.; Buch, R.; Fiorin, G.; Hénin, J.; Jiang, W.; McGreevy, R.; Melo, M. C. R.; Radak, B. K.; Skeel, R. D.; Singharoy, A.; Wang, Y.; Roux, B.; Aksimentiev, A.; Luthey-Schulten, Z.; Kalé, L. V.; Schulten, K.; Chipot, C.; Tajkhorshid, E. Scalable molecular dynamics on CPU and GPU architectures with NAMD. J. Chem. Phys. 2020, 153, 044130 DOI: 10.1063/5.0014475Google Scholar19https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhsFajsL3J&md5=aa6378c15e48addc4b6b02523e55427fScalable molecular dynamics on CPU and GPU architectures with NAMDPhillips, James C.; Hardy, David J.; Maia, Julio D. C.; Stone, John E.; Ribeiro, Joao V.; Bernardi, Rafael C.; Buch, Ronak; Fiorin, Giacomo; Henin, Jerome; Jiang, Wei; McGreevy, Ryan; Melo, Marcelo C. R.; Radak, Brian K.; Skeel, Robert D.; Singharoy, Abhishek; Wang, Yi; Roux, Benoit; Aksimentiev, Aleksei; Luthey-Schulten, Zaida; Kale, Laxmikant V.; Schulten, Klaus; Chipot, Christophe; Tajkhorshid, EmadJournal of Chemical Physics (2020), 153 (4), 044130CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A review. NAMD is a mol. dynamics program designed for high-performance simulations of very large biol. objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodn. ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomol. force fields. Here, the authors review the main features of NAMD that allow both equil. and enhanced-sampling mol. dynamics simulations with numerical efficiency. The authors describe the underlying concepts used by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temp., pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mech./mol.-mech. descriptions. The authors detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at detg. free-energy differences of either alchem. or geometrical transformations and outline their applicability to specific problems. Last, the roadmap for the development of NAMD and the authors' current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biol. realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze are discussed. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu. (c) 2020 American Institute of Physics.
- 20Mansoor, S. E.; Lü, W.; Oosterheert, W.; Shekhar, M.; Tajkhorshid, E.; Gouaux, E. X-ray structures define human P2X(3) receptor gating cycle and antagonist action. Nature 2016, 538, 66– 71, DOI: 10.1038/nature19367Google Scholar20https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhsFWiurnJ&md5=db6de173d819f23ba4304d0290ba0000X-ray structures define human P2X3 receptor gating cycle and antagonist actionMansoor, Steven E.; Lu, Wei; Oosterheert, Wout; Shekhar, Mrinal; Tajkhorshid, Emad; Gouaux, EricNature (London, United Kingdom) (2016), 538 (7623), 66-71CODEN: NATUAS; ISSN:0028-0836. (Nature Publishing Group)P2X receptors are trimeric, non-selective cation channels activated by ATP that have important roles in the cardiovascular, neuronal and immune systems. Despite their central function in human physiol. and although they are potential targets of therapeutic agents, there are no structures of human P2X receptors. The mechanisms of receptor desensitization and ion permeation, principles of antagonism, and complete structures of the pore-forming transmembrane domains of these receptors remain unclear. Here we report X-ray crystal structures of the human P2X3 receptor in apo/resting, agonist-bound/open-pore, agonist-bound/closed-pore/desensitized and antagonist-bound/closed states. The open state structure harbours an intracellular motif we term the 'cytoplasmic cap', which stabilizes the open state of the ion channel pore and creates lateral, phospholipid-lined cytoplasmic fenestrations for water and ion egress. The competitive antagonists TNP-ATP and A-317491 stabilize the apo/resting state and reveal the interactions responsible for competitive inhibition. These structures illuminate the conformational rearrangements that underlie P2X receptor gating and provide a foundation for the development of new pharmacol. agents.
- 21Wolf, M. G.; Hoefling, M.; Aponte-Santamaría, C.; Grubmüller, H.; Groenhof, G. g_membed: Efficient insertion of a membrane protein into an equilibrated lipid bilayer with minimal perturbation. J. Comput. Chem. 2010, 31, 2169– 2174, DOI: 10.1002/jcc.21507Google Scholar21https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXnsFWjs7c%253D&md5=1f0353d68fb867fb684283b61d52eadbg_membed: Efficient insertion of a membrane protein into an equilibrated lipid bilayer with minimal perturbationWolf, Maarten G.; Hoefling, Martin; Aponte-Santamaria, Camilo; Grubmueller, Helmut; Groenhof, GerritJournal of Computational Chemistry (2010), 31 (11), 2169-2174CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)To efficiently insert a protein into an equilibrated and fully hydrated membrane with minimal membrane perturbation the authors present a computational tool, called g_membed, which is part of the Gromacs suite of programs. The input consists of an equilibrated membrane system, either flat or curved, and a protein structure in the right position and orientation with respect to the lipid bilayer. G_membed first decreases the width of the protein in the xy-plane and removes all mols. (generally lipids and waters) that overlap with the narrowed protein. Then the protein is grown back to its full size in a short mol. dynamics simulation (typically 1000 steps), thereby pushing the lipids away to optimally accommodate the protein in the membrane. After embedding the protein in the membrane, both the lipid properties and the hydration layer are still close to equil. Thus, only a short equilibration run (less then 1 ns in the cases tested) is required to re-equilibrate the membrane. Its simplicity makes g_membed very practical for use in scripting and high-throughput mol. dynamics simulations. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010.
- 22Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19– 25, DOI: 10.1016/j.softx.2015.06.001Google ScholarThere is no corresponding record for this reference.
- 23Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684– 3690, DOI: 10.1063/1.448118Google Scholar23https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL2cXmtlGksbY%253D&md5=5510dc00297d63b91ee3a7a4ae5aacb1Molecular dynamics with coupling to an external bathBerendsen, H. J. C.; Postma, J. P. M.; Van Gunsteren, W. F.; DiNola, A.; Haak, J. R.Journal of Chemical Physics (1984), 81 (8), 3684-90CODEN: JCPSA6; ISSN:0021-9606.In mol. dynamics (MD) simulations, the need often arises to maintain such parameters as temp. or pressure rather than energy and vol., or to impose gradients for studying transport properties in nonequil. MD. A method is described to realize coupling to an external bath with const. temp. or pressure with adjustable time consts. for the coupling. The method is easily extendable to other variables and to gradients, and can be applied also to polyat. mols. involving internal constraints. The influence of coupling time consts. on dynamical variables is evaluated. A leap-frog algorithm is presented for the general case involving constraints with coupling to both a const. temp. and a const. pressure bath.
- 24Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101 DOI: 10.1063/1.2408420Google Scholar24https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXosVCltg%253D%253D&md5=9c182b57bfc65bca6be23c8c76b4be77Canonical sampling through velocity rescalingBussi, Giovanni; Donadio, Davide; Parrinello, MicheleJournal of Chemical Physics (2007), 126 (1), 014101/1-014101/7CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The authors present a new mol. dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains const. during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. The authors illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liq. phases. Its performance is excellent and largely independent of the thermostat parameter also with regard to the dynamic properties.
- 25Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B. L.; Grubmüller, H.; MacKerell, A. D. CHARMM36m: an improved force field for folded and intrinsically disordered proteins. Nat. Methods 2017, 14, 71– 73, DOI: 10.1038/nmeth.4067Google Scholar25https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhvVSiu77I&md5=0aa151fbef2ee0b5e2cfb593c54330c2CHARMM36m: an improved force field for folded and intrinsically disordered proteinsHuang, Jing; Rauscher, Sarah; Nawrocki, Grzegorz; Ran, Ting; Feig, Michael; de Groot, Bert L.; Grubmuller, Helmut; MacKerell, Alexander D. JrNature Methods (2017), 14 (1), 71-73CODEN: NMAEA3; ISSN:1548-7091. (Nature Publishing Group)The all-atom additive CHARMM36 protein force field is widely used in mol. modeling and simulations. We present its refinement, CHARMM36m (http://mackerell.umaryland.edu/charmm_ff.shtml), with improved accuracy in generating polypeptide backbone conformational ensembles for intrinsically disordered peptides and proteins.
- 26Klauda, J. B.; Venable, R. M.; Freites, J. A.; O’Connor, J. W.; Tobias, D. J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A. D.; Pastor, R. W. Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types. J. Phys. Chem. B 2010, 114, 7830– 7843, DOI: 10.1021/jp101759qGoogle Scholar26https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXmsVGjtrk%253D&md5=04c4d3aa1a59740190994d7df80cbb71Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid TypesKlauda, Jeffery B.; Venable, Richard M.; Freites, J. Alfredo; O'Connor, Joseph W.; Tobias, Douglas J.; Mondragon-Ramirez, Carlos; Vorobyov, Igor; MacKerell, Alexander D., Jr.; Pastor, Richard W.Journal of Physical Chemistry B (2010), 114 (23), 7830-7843CODEN: JPCBFK; ISSN:1520-6106. (American Chemical Society)A significant modification to the additive all-atom CHARMM lipid force field (FF) is developed and applied to phospholipid bilayers with both choline and ethanolamine contg. head groups and with both satd. and unsatd. aliph. chains. Motivated by the current CHARMM lipid FF (C27 and C27r) systematically yielding values of the surface area per lipid that are smaller than exptl. ests. and gel-like structures of bilayers well above the gel transition temp., selected torsional, Lennard-Jones and partial at. charge parameters were modified by targeting both quantum mech. (QM) and exptl. data. QM calcns. ranging from high-level ab initio calcns. on small mols. to semiempirical QM studies on a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer in combination with exptl. thermodn. data were used as target data for parameter optimization. These changes were tested with simulations of pure bilayers at high hydration of the following six lipids: DPPC, 1,2-dimyristoyl-sn-phosphatidylcholine (DMPC), 1,2-dilauroyl-sn-phosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-phosphatidylcholine (POPC), 1,2-dioleoyl-sn-phosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-sn-phosphatidylethanolamine (POPE); simulations of a low hydration DOPC bilayer were also performed. Agreement with exptl. surface area is on av. within 2%, and the d. profiles agree well with neutron and x-ray diffraction expts. NMR deuterium order parameters (SCD) are well predicted with the new FF, including proper splitting of the SCD for the aliph. carbon adjacent to the carbonyl for DPPC, POPE, and POPC bilayers. The area compressibility modulus and frequency dependence of 13C NMR relaxation rates of DPPC and the water distribution of low hydration DOPC bilayers also agree well with expt. Accordingly, the presented lipid FF, referred to as C36, allows for mol. dynamics simulations to be run in the tensionless ensemble (NPT), and is anticipated to be of utility for simulations of pure lipid systems as well as heterogeneous systems including membrane proteins.
- 27Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926– 935, DOI: 10.1063/1.445869Google Scholar27https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL3sXksF2htL4%253D&md5=a1161334e381746be8c9b15a5e56f704Comparison of simple potential functions for simulating liquid waterJorgensen, William L.; Chandrasekhar, Jayaraman; Madura, Jeffry D.; Impey, Roger W.; Klein, Michael L.Journal of Chemical Physics (1983), 79 (2), 926-35CODEN: JCPSA6; ISSN:0021-9606.Classical Monte Carlo simulations were carried out for liq. H2O in the NPT ensemble at 25° and 1 atm using 6 of the simpler intermol. potential functions for the dimer. Comparisons were made with exptl. thermodn. and structural data including the neutron diffraction results of Thiessen and Narten (1982). The computed densities and potential energies agree with expt. except for the original Bernal-Fowler model, which yields an 18% overest. of the d. and poor structural results. The discrepancy may be due to the correction terms needed in processing the neutron data or to an effect uniformly neglected in the computations. Comparisons were made for the self-diffusion coeffs. obtained from mol. dynamics simulations.
- 28Brunner, J. D.; Lim, N. K.; Schenck, S.; Duerst, A.; Dutzler, R. X-ray structure of a calcium-activated TMEM16 lipid scramblase. Nature 2014, 516, 207– 212, DOI: 10.1038/nature13984Google Scholar28https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhvFGltLrO&md5=828d54af1de0c0bd05a31a1eed8a6fdfX-ray structure of a calcium-activated TMEM16 lipid scramblaseBrunner, Janine D.; Lim, Novandy K.; Schenck, Stephan; Duerst, Alessia; Dutzler, RaimundNature (London, United Kingdom) (2014), 516 (7530), 207-212CODEN: NATUAS; ISSN:0028-0836. (Nature Publishing Group)The TMEM16 family of proteins, also known as anoctamins, features a remarkable functional diversity. This family contains the long sought-after Ca2+-activated chloride channels as well as lipid scramblases and cation channels. Here we present the crystal structure of a TMEM16 family member from the fungus Nectria hematococca that operates as a Ca2+-activated lipid scramblase. Each subunit of the homodimeric protein contains ten transmembrane helixes and a hydrophilic membrane-traversing cavity that is exposed to the lipid bilayer as a potential site of catalysis. This cavity harbors a conserved Ca2+-binding site located within the hydrophobic core of the membrane. Mutations of residues involved in Ca2+ coordination affect both lipid scrambling in N. hematococca TMEM16 and ion conduction in the Cl- channel TMEM16A. The structure reveals the general architecture of the family and its mode of Ca2+ activation. It also provides insight into potential scrambling mechanisms and serves as a framework to unravel the conduction of ions in certain TMEM16 proteins.
- 29Bushell, S. R.; Pike, A. C.; Falzone, M. E.; Rorsman, N. J.; Ta, C. M.; Corey, R. A.; Newport, T. D.; Christianson, J. C.; Scofano, L. F.; Shintre, C. A.; Tessitore, A.; Chu, A.; Wang, Q.; Shrestha, L.; Mukhopadhyay, S. M. M.; Love, J. D.; Burgess-Brown, N. A.; Sitsapesan, R.; Stansfeld, P. J.; Huiskonen, J. T.; Tammaro, P.; Accardi, A.; Carpenter, E. P. The structural basis of lipid scrambling and inactivation in the endoplasmic reticulum scramblase TMEM16K. Nat. Commun. 2019, 10, 3956 DOI: 10.1038/s41467-019-11753-1Google Scholar29https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3MrksVCrtw%253D%253D&md5=f146bbc9a998091e2825f45ff419a859The structural basis of lipid scrambling and inactivation in the endoplasmic reticulum scramblase TMEM16KBushell Simon R; Pike Ashley C W; Shintre Chitra A; Tessitore Annamaria; Chu Amy; Wang Qinrui; Shrestha Leela; Mukhopadhyay Shubhashish M M; Burgess-Brown Nicola A; Carpenter Elisabeth P; Falzone Maria E; Accardi Alessio; Rorsman Nils J G; Ta Chau M; Scofano Lara F; Sitsapesan Rebecca; Tammaro Paolo; Rorsman Nils J G; Ta Chau M; Corey Robin A; Newport Thomas D; Wang Qinrui; Stansfeld Phillip J; Newport Thomas D; Christianson John C; Shintre Chitra A; Tessitore Annamaria; Chu Amy; Love James D; Love James D; Huiskonen Juha T; Accardi Alessio; Accardi AlessioNature communications (2019), 10 (1), 3956 ISSN:.Membranes in cells have defined distributions of lipids in each leaflet, controlled by lipid scramblases and flip/floppases. However, for some intracellular membranes such as the endoplasmic reticulum (ER) the scramblases have not been identified. Members of the TMEM16 family have either lipid scramblase or chloride channel activity. Although TMEM16K is widely distributed and associated with the neurological disorder autosomal recessive spinocerebellar ataxia type 10 (SCAR10), its location in cells, function and structure are largely uncharacterised. Here we show that TMEM16K is an ER-resident lipid scramblase with a requirement for short chain lipids and calcium for robust activity. Crystal structures of TMEM16K show a scramblase fold, with an open lipid transporting groove. Additional cryo-EM structures reveal extensive conformational changes from the cytoplasmic to the ER side of the membrane, giving a state with a closed lipid permeation pathway. Molecular dynamics simulations showed that the open-groove conformation is necessary for scramblase activity.
- 30Kostritskii, A. Y.; Machtens, J.-P. Molecular mechanisms of ion conduction and ion selectivity in TMEM16 lipid scramblases. Nat. Commun. 2021, in press, DOI: 10.1038/s41467-021-22724-w .Google ScholarThere is no corresponding record for this reference.
- 31Alleva, C.; Kovalev, K.; Astashkin, R.; Berndt, M. I.; Baeken, C.; Balandin, T.; Gordeliy, V.; Fahlke, C.; Machtens, J.-P. Na+-dependent gate dynamics and electrostatic attraction ensure substrate coupling in glutamate transporters. Sci. Adv. 2020, 6, eaba9854 DOI: 10.1126/sciadv.aba9854Google ScholarThere is no corresponding record for this reference.
- 32Lomize, M. A.; Pogozheva, I. D.; Joo, H.; Mosberg, H. I.; Lomize, A. L. OPM database and PPM web server: resources for positioning of proteins in membranes. Nucleic Acids Res. 2012, 40, D370– D376, DOI: 10.1093/nar/gkr703Google Scholar32https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhs12hurzJ&md5=b931126c56d227ab0eb6c5762c1dae6bOPM database and PPM web server: resources for positioning of proteins in membranesLomize, Mikhail A.; Pogozheva, Irina D.; Joo, Hyeon; Mosberg, Henry I.; Lomize, Andrei L.Nucleic Acids Research (2012), 40 (D1), D370-D376CODEN: NARHAD; ISSN:0305-1048. (Oxford University Press)The Orientations of Proteins in Membranes (OPM) database is a curated web resource that provides spatial positions of membrane-bound peptides and proteins of known three-dimensional structure in the lipid bilayer, together with their structural classification, topol. and intracellular localization. OPM currently contains more than 1200 transmembrane and peripheral proteins and peptides from approx. 350 organisms that represent approx. 3800 Protein Data Bank entries. Proteins are classified into classes, superfamilies and families and assigned to 21 distinct membrane types. Spatial positions of proteins with respect to the lipid bilayer are optimized by the PPM 2.0 method that accounts for the hydrophobic, hydrogen bonding and electrostatic interactions of the proteins with the anisotropic water-lipid environment described by the dielec. const. and hydrogen-bonding profiles. The OPM database is freely accessible at http://opm.phar.umich.edu. Data can be sorted, searched or retrieved using the hierarchical classification, source organism, localization in different types of membranes. The database offers downloadable coordinates of proteins and peptides with membrane boundaries. A gallery of protein images and several visualization tools are provided. The database is supplemented by the PPM server (http://opm.phar.umich.edu/server.php) which can be used for calcg. spatial positions in membranes of newly detd. proteins structures or theor. models.
- 33Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J. L.; Dror, R. O.; Shaw, D. E. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins: Struct., Funct., Bioinf. 2010, 78, 1950– 1958, DOI: 10.1002/prot.22711Google Scholar33https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXkvFegtLo%253D&md5=447a9004026e2b93f0f7beff165daa09Improved side-chain torsion potentials for the Amber ff99SB protein force fieldLindorff-Larsen, Kresten; Piana, Stefano; Palmo, Kim; Maragakis, Paul; Klepeis, John L.; Dror, Ron O.; Shaw, David E.Proteins: Structure, Function, and Bioinformatics (2010), 78 (8), 1950-1958CODEN: PSFBAF ISSN:. (Wiley-Liss, Inc.)Recent advances in hardware and software have enabled increasingly long mol. dynamics (MD) simulations of biomols., exposing certain limitations in the accuracy of the force fields used for such simulations and spurring efforts to refine these force fields. Recent modifications to the Amber and CHARMM protein force fields, for example, have improved the backbone torsion potentials, remedying deficiencies in earlier versions. Here, the authors further advance simulation accuracy by improving the amino acid side-chain torsion potentials of the Amber ff99SB force field. First, the authors used simulations of model alpha-helical systems to identify the four residue types whose rotamer distribution differed the most from expectations based on Protein Data Bank statistics. Second, the authors optimized the side-chain torsion potentials of these residues to match new, high-level quantum-mech. calcns. Finally, the authors used microsecond-timescale MD simulations in explicit solvent to validate the resulting force field against a large set of exptl. NMR measurements that directly probe side-chain conformations. The new force field, which the authors have termed Amber ff99SB-ILDN, exhibits considerably better agreement with the NMR data. Proteins 2010. © 2010 Wiley-Liss, Inc.
- 34Berger, O.; Edholm, O.; Jähnig, F. Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. Biophys. J. 1997, 72, 2002– 2013, DOI: 10.1016/S0006-3495(97)78845-3Google Scholar34https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXivVOjsL4%253D&md5=984ea49ced9d0e2ce912bf99220fd228Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperatureBerger, Oliver; Edholm, Olle; Jahnig, FritzBiophysical Journal (1997), 72 (5), 2002-2013CODEN: BIOJAU; ISSN:0006-3495. (Biophysical Society)Mol. dynamics simulations of 500 ps were performed on a system consisting of a bilayer of 64 mols. of the lipid dipalmitoylphosphatidylcholine and 23 water mols. per lipid at an isotropic pressure of 1 atm and 50°. Special attention was devoted to reproduce the correct d. of the lipid, because this quantity is known exptl. with a precision better than 1%. For this purpose, the Lennard-Jones parameters of the hydrocarbon chains were adjusted by simulating a system consisting of 128 pentadecane mols. and varying the Lennard-Jones parameters until the exptl. d. and heat of vaporization were obtained. With these parameters the lipid d. resulted in perfect agreement with the exptl. d. The orientational order parameter of the hydrocarbon chains agreed perfectly well with the exptl. values, which, because of its correlation with the area per lipid, makes it possible to give a proper est. of the area per lipid of 0.61±0.1 nm2.
- 35Joung, I. S.; Cheatham, T. E. Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations. J. Phys. Chem. B 2008, 112, 9020– 9041, DOI: 10.1021/jp8001614Google Scholar35https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXnvFGqtL4%253D&md5=aa489470ae1c7479bf0911710217bd28Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular SimulationsJoung, In Suk; Cheatham, Thomas E.Journal of Physical Chemistry B (2008), 112 (30), 9020-9041CODEN: JPCBFK; ISSN:1520-6106. (American Chemical Society)Alkali (Li+, Na+, K+, Rb+, and Cs+) and halide (F-, Cl-, Br-, and I-) ions play an important role in many biol. phenomena, roles that range from stabilization of biomol. structure, to influence on biomol. dynamics, to key physiol. influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomol. structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in soln. and the interactions of ions with other mols. At present, the best force fields for biomols. employ a simple additive, nonpolarizable, and pairwise potential for at. interaction. In this work, the authors describe their efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and soln. properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mech. treatment, the authors' goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodol. is general and can be extended to other ions and to polarizable force-field models. The authors' starting point centered on observations from long simulations of biomols. in salt soln. with the AMBER force fields where salt crystals formed well below their soly. limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Aqvist cation parameters. To provide a more appropriate balance, the authors reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, the authors calcd. hydration free energies of the solvated ions and also lattice energies (LE) and lattice consts. (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4PEW, and SPC/E. In addn. to well reproducing the soln. and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells.
- 36Drew, H. R.; Dickerson, R. E. Structure of a B-DNA dodecamer. III. Geometry of hydration. J. Mol. Biol. 1981, 151, 535– 556, DOI: 10.1016/0022-2836(81)90009-7Google Scholar36https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL38XmsVentg%253D%253D&md5=4378664a4da9ec83155a3f7e046027d9Structure of a B-DNA dodecamer. III. Geometry of hydrationDrew, Horace R.; Dickerson, Richard E.Journal of Molecular Biology (1981), 151 (3), 535-56CODEN: JMOBAK; ISSN:0022-2836.The B-DNA dodecamer, C-G-C-G-A-A-T-T-C-G-C-G, crystd. as slightly more than 1 full turn of right-handed B-DNA in space group P212121 with cell dimensions a = 24.87 Å, b = 40.39 Å, and c = 66.20 Å. X-ray anal. showed that it was surrounded by 72 ordered water mols., mostly assocd. with polar N and O atoms at exposed edges of base-pairs. Hydration in the major groove was mainly in the form of a monodentate monolayer. Hydration of backbone phosphate O atoms were not ordered, except when immobilized by the 5-Me groups of adjacent thymines. The minor groove was extensively and regularly hydrated, with a zigzag spine of 1st- and 2nd-shell hydration along the floor of the groove serving as a foundation for less-regular outer shells extending beyond the radius of the phosphate backbone. The spine network bridged purine N-3 and pyrimidine O-2 atoms in adjacent base pairs; it was regular in the A-A-T-T center, but was disrupted at the C-G-C-G ends, partly by the guanine N-2 amino groups. The minor groove hydration spine may be responsible for the stability of the B form of polymers contg. only A-T and I-C base pairs, and its disruption may explain the ease of transition to the A form of polymers with G-C pairs.
- 37Pérez, A.; Marchán, I.; Svozil, D.; Sponer, J.; Cheatham, T. E.; Laughton, C. A.; Orozco, M. Refinement of the AMBER force field for nucleic acids: improving the description of alpha/gamma conformers. Biophys. J. 2007, 92, 3817– 3829, DOI: 10.1529/biophysj.106.097782Google Scholar37https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXls1ygu7Y%253D&md5=e0f8d87e5fe7710e59cc41f5e8c25f7bRefinement of the AMBER force field for nucleic acids: improving the description of α/γ conformersPerez, Alberto; Marchan, Ivan; Svozil, Daniel; Sponer, Jiri; Cheatham, Thomas E., III; Laughton, Charles A.; Orozco, ModestoBiophysical Journal (2007), 92 (11), 3817-3829CODEN: BIOJAU; ISSN:0006-3495. (Biophysical Society)We present here the parmbsc0 force field, a refinement of the AMBER parm99 force field, where emphasis has been made on the correct representation of the α/γ concerted rotation in nucleic acids (NAs). The modified force field corrects overpopulations of the α/γ = (g+,t) backbone that were seen in long (more than 10 ns) simulations with previous AMBER parameter sets (parm94-99). The force field has been derived by fitting to high-level quantum mech. data and verified by comparison with very high-level quantum mech. calcns. and by a very extensive comparison between simulations and exptl. data. The set of validation simulations includes two of the longest trajectories published to date for the DNA duplex (200 ns each) and the largest variety of NA structures studied to date (15 different NA families and 97 individual structures). The total simulation time used to validate the force field includes near 1 μs of state-of-the-art mol. dynamics simulations in aq. soln.
- 38Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182– 7190, DOI: 10.1063/1.328693Google Scholar38https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL38XislSnuw%253D%253D&md5=a0a5617389f6cabbf2a405c649aadf03Polymorphic transitions in single crystals: A new molecular dynamics methodParrinello, M.; Rahman, A.Journal of Applied Physics (1981), 52 (12), 7182-90CODEN: JAPIAU; ISSN:0021-8979.A Lagrangian formulation is introduced; it can be used to make mol. dynamics (MD) calcns. on systems under the most general, externally applied, conditions of stress. In this formulation the MD cell shape and size can change according to dynamic equations given by this Lagrangian. This MD technique was used to the study of structural transitions of a Ni single crystal under uniform uniaxial compressive and tensile loads. Some results regarding the stress-strain relation obtained by static calcns. are invalid at finite temp. Under compressive loading, the model of Ni shows a bifurcation in its stress-strain relation; this bifurcation provides a link in configuration space between cubic and hexagonal close packing. Such a transition could perhaps be obsd. exptl. under extreme conditions of shock.
- 39Kondinskaia, D. A.; Kostritskii, A. Y.; Nesterenko, A. M.; Antipina, A. Y.; Gurtovenko, A. A. Atomic-Scale Molecular Dynamics Simulations of DNA-Polycation Complexes: Two Distinct Binding Patterns. J. Phys. Chem. B 2016, 120, 6546– 6554, DOI: 10.1021/acs.jpcb.6b03779Google Scholar39https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XpsFGlsLg%253D&md5=ab05858cf4c832df5a5027e197842427Atomic-Scale Molecular Dynamics Simulations of DNA-Polycation Complexes: Two Distinct Binding PatternsKondinskaia, Diana A.; Kostritskii, Andrei Yu.; Nesterenko, Alexey M.; Antipina, Alexandra Yu.; Gurtovenko, Andrey A.Journal of Physical Chemistry B (2016), 120 (27), 6546-6554CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Synthetic cationic polymers represent a promising class of delivery vectors for gene therapy. Here, we employ atomistic mol. dynamics simulations to gain insight into the structure and properties of complexes of DNA with four linear polycations: polyethylenimine (PEI), poly-L-lysine (PLL), polyvinylamine (PVA), and polyallylamine (PAA). These polycations differ in their polymer geometries, protonation states, and hydrophobicities of their backbone chains. Overall, our results demonstrate for the first time the existence of two distinct patterns of binding of DNA with polycations. For PEI, PLL, and PAA, the complex is stabilized by the electrostatic attraction between protonated amine groups of the polycation and phosphate groups of DNA. In contrast, PVA demonstrates an alternative binding pattern as it gets embedded into the DNA major groove. It is likely that both the polymer topol. and affinity of the backbone chain of PVA to the DNA groove are responsible for such behavior. The differences in binding patterns can have important biomedical implications: embedding PVA into a DNA groove makes it less sensitive to changes in the aq. environment (pH level, ionic strength, etc.) and could therefore hinder the intracellular release of genetic material from a delivery vector, leading to lower transfection activity.
- 40Michaud-Agrawal, N.; Denning, E. J.; Woolf, T. B.; Beckstein, O. MDAnalysis: a toolkit for the analysis of molecular dynamics simulations. J. Comput. Chem. 2011, 32, 2319– 2327, DOI: 10.1002/jcc.21787Google Scholar40https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXnvFalsr8%253D&md5=d567042c65cfdc1c81336a29137654bfMDAnalysis: A toolkit for the analysis of molecular dynamics simulationsMichaud-Agrawal, Naveen; Denning, Elizabeth J.; Woolf, Thomas B.; Beckstein, OliverJournal of Computational Chemistry (2011), 32 (10), 2319-2327CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)MDAnal. is an object-oriented library for structural and temporal anal. of mol. dynamics (MD) simulation trajectories and individual protein structures. It is written in the Python language with some performance-crit. code in C. It uses the powerful NumPy package to expose trajectory data as fast and efficient NumPy arrays. It has been tested on systems of millions of particles. Many common file formats of simulation packages including CHARMM, Gromacs, Amber, and NAMD and the Protein Data Bank format can be read and written. Atoms can be selected with a syntax similar to CHARMM's powerful selection commands. MDAnal. enables both novice and experienced programmers to rapidly write their own anal. tools and access data stored in trajectories in an easily accessible manner that facilitates interactive explorative anal. MDAnal. has been tested on and works for most Unix-based platforms such as Linux and Mac OS X. It is freely available under the GNU General Public License from http://mdanal.googlecode.com. © 2011 Wiley Periodicals, Inc. J Comput Chem 2011.
- 41North, R. A. Molecular physiology of P2X receptors. Physiol. Rev. 2002, 82, 1013– 1067, DOI: 10.1152/physrev.00015.2002Google Scholar41https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38Xot1Ggu78%253D&md5=d621c6ce761f5e2db7ec553ed0f5ac10Molecular physiology of P2X receptorsNorth, R. AlanPhysiological Reviews (2002), 82 (4), 1013-1067CODEN: PHREA7; ISSN:0031-9333. (American Physiological Society)A review. P2X receptors are membrane ion channels that open in response to the binding of extracellular ATP. Seven genes in vertebrates encode P2X receptor subunits, which are 40-50% identical in amino acid sequence. Each subunit has two transmembrane domains, sepd. by an extracellular domain (∼280 amino acids). Channels form as multimers of several subunits. Homomeric P2X1, P2X2, P2X3, P2X4, P2X5, and P2X7 channels and heteromeric P2X2/3 and P2X1/5 channels have been most fully characterized following heterologous expression. Some agonists (e.g., αβ-methylene ATP) and antagonists [e.g., 2',3'-O-(2,4,6-trinitrophenyl)-ATP] are strongly selective for receptors contg. P2X1 and P2X3 subunits. All P2X receptors are permeable to small monovalent cations; some have significant calcium or anion permeability. In many cells, activation of homomeric P2X7 receptors induces a permeability increase to larger org. cations including some fluorescent dyes and also signals to the cytoskeleton; these changes probably involve addnl. interacting proteins. P2X receptors are abundantly distributed, and functional responses are seen in neurons, glia, epithelia, endothelia, bone, muscle, and hemopoietic tissues. The mol. compn. of native receptors is becoming understood, and some cells express more than one type of P2X receptor. On smooth muscles, P2X receptors respond to ATP released from sympathetic motor nerves (e.g., in ejaculation). On sensory nerves, they are involved in the initiation of afferent signals in several viscera (e.g., bladder, intestine) and play a key role in sensing tissue-damaging and inflammatory stimuli. Paracrine roles for ATP signaling through P2X receptors are likely in neurohypophysis, ducted glands, airway epithelia, kidney, bone, and hemopoietic tissues. In the last case, P2X7 receptor activation stimulates cytokine release by engaging intracellular signaling pathways.
- 42Surprenant, A.; North, R. A. Signaling at purinergic P2X receptors. Annu. Rev. Physiol. 2009, 71, 333– 359, DOI: 10.1146/annurev.physiol.70.113006.100630Google Scholar42https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXjsValsbk%253D&md5=c9f3772f8f746ea0591ff30a3c7beb36Signaling at purinergic P2X receptorsSurprenant, Annmarie; North, R. AlanAnnual Review of Physiology (2009), 71 (), 333-359CODEN: ARPHAD; ISSN:0066-4278. (Annual Reviews Inc.)A review. P2X receptors are membrane cation channels gated by extracellular ATP. Seven P2X receptor subunits (P2X1-7) are widely distributed in excitable and nonexcitable cells of vertebrates. They play key roles in inter alia afferent signaling (including pain), regulation of renal blood flow, vascular endothelium, and inflammatory responses. We summarize the evidence for these and other roles, emphasizing exptl. work with selective receptor antagonists or with knockout mice. The receptors are trimeric membrane proteins: studies of the biophys. properties of mutated subunits expressed in heterologous cells have indicated parts of the subunits involved in ATP binding, ion permeation (including calcium permeability), and membrane trafficking. We review our current understanding of the mol. properties of P2X receptors, including how this understanding is informed by the identification of distantly related P2X receptors in simple eukaryotes.
- 43Illes, P.; Müller, C. E.; Jacobson, K. A.; Grutter, T.; Nicke, A.; Fountain, S. J.; Kennedy, C.; Schmalzing, G.; Jarvis, M. F.; Stojilkovic, S. S. Update of P2X receptor properties and their pharmacology: IUPHAR Review 30. Br. J. Pharmacol. 2021, 178, 489– 514, DOI: 10.1111/bph.15299Google Scholar43https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXktFehuw%253D%253D&md5=6cb4be20f4b37e5d18e523bdf7f1aef3Update of P2X receptor properties and their pharmacology: IUPHAR Review 30Illes, Peter; Mueller, Christa E.; Jacobson, Kenneth A.; Grutter, Thomas; Nicke, Annette; Fountain, Samuel J.; Kennedy, Charles; Schmalzing, Guenther; Jarvis, Michael F.; Stojilkovic, Stanko S.; King, Brian F.; Di Virgilio, FrancescoBritish Journal of Pharmacology (2021), 178 (3), 489-514CODEN: BJPCBM; ISSN:1476-5381. (Wiley-Blackwell)A review. The known seven mammalian receptor subunits (P2X1-7) form cationic channels gated by ATP. Three subunits compose a receptor channel. Each subunit is a polypeptide consisting of two transmembrane regions (TM1 and TM2), intracellular N- and C-termini, and a bulky extracellular loop. Crystn. allowed the identification of the 3D structure and gating cycle of P2X receptors. The agonist-binding pocket is located at the intersection of two neighboring subunits. In addn. to the mammalian P2X receptors, their primitive ligand-gated counterparts with little structural similarity have also been cloned. Selective agonists for P2X receptor subtypes are not available, but medicinal chem. supplied a range of subtype-selective antagonists, as well as pos. and neg. allosteric modulators. Knockout mice and selective antagonists helped to identify pathol. functions due to defective P2X receptors, such as male infertility (P2X1), hearing loss (P2X2), pain/cough (P2X3), neuropathic pain (P2X4), inflammatory bone loss (P2X5), and faulty immune reactions (P2X7).
- 44Briones, R.; Blau, C.; Kutzner, C.; de Groot, B. L.; Aponte-Santamaría, C. Gromaρs: A Gromacs-Based Toolset to Analyse Density Maps Derived from Molecular Dynamics Simulations. Biophys. J. 2019, 116, 4– 11, DOI: 10.1016/j.bpj.2018.11.3126Google Scholar44https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXisVKitrbK&md5=7e63b96203c0e8aa9a182c4098310eb5GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics SimulationsBriones, Rodolfo; Blau, Christian; Kutzner, Carsten; de Groot, Bert L.; Aponte-Santamaria, CamiloBiophysical Journal (2019), 116 (1), 4-11CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)The authors introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged d. maps from mol. dynamics simulations. GROmaρs efficiently computes d. maps by fast multi-Gaussian spreading of at. densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of at. av. localization during the simulations. Most importantly, it allows the comparison between computed and ref. maps (e.g., exptl.) through calcn. of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quant. contrast perturbed and control simulation data sets and to examine how much biomol. systems resemble both synthetic and exptl. d. maps. This was esp. advantageous for multimol. systems in which std. comparisons like RMSDs are difficult to compute. In addn., GROmaρs incorporates abs. and relative spatial free-energy ests. to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the d. calcn. Furthermore, masking of regions was implemented to speed up calcns. and to facilitate the comparison with exptl. maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and av. charge distribution in biomol. systems. The authors employed all these functionalities to inspect the localization of lipid and water mols. in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, the authors anticipate a high applicability of GROmaρs for the anal. of mol. dynamics simulations and their comparison with exptl. detd. densities.
- 45Chen, X.; Bründl, M.; Friesacher, T.; Stary-Weinzinger, A. Computational Insights Into Voltage Dependence of Polyamine Block in a Strong Inwardly Rectifying K+ Channel. Front. Pharmacol. 2020, 11, 721 DOI: 10.3389/fphar.2020.00721Google Scholar45https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhvFKktb7J&md5=e735c85eef42803dfb799d1736d99774Computational insights into voltage dependence of polyamine block in a strong inwardly rectifying K+ channelChen, Xingyu; Bruendl, Michael; Friesacher, Theres; Stary-Weinzinger, AnnaFrontiers in Pharmacology (2020), 11 (), 00721CODEN: FPRHAU; ISSN:1663-9812. (Frontiers Media S.A.)Inwardly rectifying potassium (KIR) channels play important roles in controlling cellular excitability and K+ ion homeostasis. Under physiol. conditions, KIR channels allow large K+ influx at potentials neg. to the equil. potential of K+ but permit little outward current at potentials pos. to the equil. potential of K+, due to voltage dependent block of outward K+ flux by cytoplasmic polyamines. These polycationic mols. enter the KIR channel pore from the intracellular side. They block K+ ion movement through the channel at depolarized potentials, thereby ensuring, for instance, the long plateau phase of the cardiac action potential. Key questions concerning how deeply these charged mols. migrate into the pore and how the steep voltage dependence arises remain unclear. Recent MD simulations on GIRK2 (=Kir3.2) crystal structures have provided unprecedented details concerning the conduction mechanism of a KIR channel. Here, we use MD simulations with applied field to provide detailed insights into voltage dependent block of putrescine, using the conductive state of the strong inwardly rectifying K+ channel GIRK2 as starting point. Our μs long simulations elucidate details about binding sites of putrescine in the pore and suggest that voltage dependent rectification arises from a dual mechanism.
- 46Nyblom, M.; Frick, A.; Wang, Y.; Ekvall, M.; Hallgren, K.; Hedfalk, K.; Neutze, R.; Tajkhorshid, E.; Törnroth-Horsefield, S. Structural and functional analysis of SoPIP2;1 mutants adds insight into plant aquaporin gating. J. Mol. Biol. 2009, 387, 653– 668, DOI: 10.1016/j.jmb.2009.01.065Google Scholar46https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXjsVOnu7c%253D&md5=7176199fe71199e1ea750db551d14c46Structural and Functional Analysis of SoPIP2;1 Mutants Adds Insight into Plant Aquaporin GatingNyblom, Maria; Frick, Anna; Wang, Yi; Ekvall, Mikael; Hallgren, Karin; Hedfalk, Kristina; Neutze, Richard; Tajkhorshid, Emad; Toernroth-Horsefield, SusannaJournal of Molecular Biology (2009), 387 (3), 653-668CODEN: JMOBAK; ISSN:0022-2836. (Elsevier Ltd.)Plant plasma membrane aquaporins facilitate water flux into and out of plant cells, thus coupling their cellular function to basic aspects of plant physiol. Posttranslational modifications of conserved phosphorylation sites, changes in cytoplasmic pH, and the binding of Ca2+ can regulate water transport activity by gating the plasma membrane aquaporins. A structural mechanism unifying these diverse biochem. signals has emerged for the spinach aquaporin SoPIP2;1, although several questions concerning the opening mechanism remain. Here, we describe the X-ray structures of the S115E and S274E single SoPIP2;1 mutants and the corresponding double mutant. Phosphorylation of these serines is believed to increase water transport activity of SoPIP2;1 by opening the channel. However, all mutants crystd. in a closed conformation, as confirmed by water transport assays, implying that neither substitution fully mimics the phosphorylated state. Nevertheless, a half-turn extension of transmembrane helix 1 occurs upon the substitution of Ser115, which draws the Cα atom of Glu31 10 Å away from its wild-type conformation, thereby disrupting the divalent cation binding site involved in the gating mechanism. Mutation of Ser274 disorders the C-terminus but no other significant conformational changes are obsd. Inspection of the hydrogen-bond interactions within loop D suggested that the phosphorylation of Ser188 may also produce an open channel, and this was supported by an increased water transport activity for the S188E mutant and mol. dynamics simulations. These findings add addnl. insight into the general mechanism of plant aquaporin gating.
- 47Ruiz, L.; Wu, Y.; Keten, S. Tailoring the water structure and transport in nanotubes with tunable interiors. Nanoscale 2015, 7, 121– 132, DOI: 10.1039/C4NR05407EGoogle Scholar47https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhslGisr%252FI&md5=9224107648d21b753e8d61e94101d943Tailoring the water structure and transport in nanotubes with tunable interiorsRuiz, Luis; Wu, Yuanqiao; Keten, SinanNanoscale (2015), 7 (1), 121-132CODEN: NANOHL; ISSN:2040-3372. (Royal Society of Chemistry)Self-assembly of cyclic peptide nanotubes (CPNs) in polymer thin films has opened up the possibility of creating sepn. membranes with tunable nanopores that can differentiate mols. at the sub-nanometer level. While it has been demonstrated that the interior chem. of the CPNs can be tailored by inserting functional groups in the nanopore lumen (mCPNs), a design strategy for picking the chem. modifications that lead to particular transport properties has not been established. Drawing from the knowledgebase of functional groups in natural amino acids, here we use mol. dynamics simulations to elucidate how bioinspired mutations influence the transport of water through mCPNs. We show that, at the nanoscale, factors besides the pore size, such as electrostatic interactions and steric effects, can dramatically change the transport properties. We recognize a novel asym. structure of water under nanoconfinement inside the chem. functionalized nanotubes and identify that the small non-polar glycine-mimic groups that minimize the steric constraints and confer a hydrophobic character to the nanotube interior are the fastest transporters of water. Our computationally developed expts. on a realistic material system circumvent synthetic challenges, and lay the foundation for bioinspired principles to tailor artificial nanochannels for sepn. applications such as desalination, ion-exchange and carbon capture.
- 48Kasimova, M. A.; Lindahl, E.; Delemotte, L. Determining the molecular basis of voltage sensitivity in membrane proteins. J. Gen. Physiol. 2018, 150, 1444– 1458, DOI: 10.1085/jgp.201812086Google Scholar48https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXmtlSrsL4%253D&md5=467ecade0184d6cc82276971a54ce988Determining the molecular basis of voltage sensitivity in membrane proteinsKasimova, Marina A.; Lindahl, Erik; Delemotte, LucieJournal of General Physiology (2018), 150 (10), 1444-1458CODEN: JGPLAD; ISSN:1540-7748. (Rockefeller University Press)Voltage-sensitive membrane proteins are united by their ability to transform changes in membrane potential into mech. work. They are responsible for a spectrum of physiol. processes in living organisms, including elec. signaling and cell-cycle progression. Moreover, the detection of these elements by using exptl. techniques is challenging because of the diversity of membrane proteins. Here, we provide a computational approach to predict voltage-sensing elements in any membrane protein, independent of its structure or function. It relies on an estn. of the propensity of a protein to respond to changes in membrane potential. We first show that this property correlates well with voltage sensitivity by applying our approach to a set of voltage-sensitive and voltage-insensitive membrane proteins. We further show that it correctly identifies authentic voltage-sensitive residues in the voltage-sensor domain of voltage-gated ion channels. Finally, we investigate six membrane proteins for which the voltage-sensing elements have not yet been characterized and identify residues and ions that might be involved in the response to voltage. The suggested approach is fast and simple and enables a characterization of voltage sensitivity that goes beyond mere identification of charges. We anticipate that its application before mutagenesis expts. will significantly reduce the no. of potential voltage-sensitive elements to be tested.
- 49Scheurer, M.; Rodenkirch, P.; Siggel, M.; Bernardi, R. C.; Schulten, K.; Tajkhorshid, E.; Rudack, T. PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations. Biophys. J. 2018, 114, 577– 583, DOI: 10.1016/j.bpj.2017.12.003Google Scholar49https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXis1Oguw%253D%253D&md5=cdcf9bda51716dae11aa7575c13abb61PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD SimulationsScheurer, Maximilian; Rodenkirch, Peter; Siggel, Marc; Bernardi, Rafael C.; Schulten, Klaus; Tajkhorshid, Emad; Rudack, TillBiophysical Journal (2018), 114 (3), 577-583CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)Mol. dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth gave a large amt. of simulation data that need to be filtered for relevant information to address specific biomedical and biochem. questions. One of the most relevant mol. properties that can be studied by all-atom MD simulations is the time-dependent evolution of the complex noncovalent interaction networks governing such fundamental aspects as mol. recognition, binding strength, and mech. and structural stability. Extg., evaluating, and visualizing noncovalent interactions is a key task in the daily work of structural biologists. The authors have developed PyContact, an easy-to-use, highly flexible, and intuitive graphical user interface-based application, designed to provide a toolkit to study biomol. interactions in MD trajectories. PyContact is designed to facilitate this task by enabling identification of relevant noncovalent interactions in a comprehensible manner. The implementation of PyContact as a standalone application enables rapid anal. and data visualization without any addnl. programming requirements, and also preserves full in-program customization and extension capabilities for advanced users. The statistical anal. representation is interactively combined with full mapping of the results on the mol. system through the synergistic connection between PyContact and VMD. The authors showcase the capabilities and scientific significance of PyContact by analyzing and visualizing in great detail the noncovalent interactions underlying the ion permeation pathway of the human P2X3 receptor. As a second application, the authors examine the protein-protein interaction network of the mech. ultrastable cohesin-dockering complex.
- 50Melcr, J.; Martinez-Seara, H.; Nencini, R.; Kolafa, J.; Jungwirth, P.; Ollila, O. H. S. Accurate Binding of Sodium and Calcium to a POPC Bilayer by Effective Inclusion of Electronic Polarization. J. Phys. Chem. B 2018, 122, 4546– 4557, DOI: 10.1021/acs.jpcb.7b12510Google Scholar50https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXmvVaqtro%253D&md5=a49868ee725f486dcc49522bca8038d0Accurate Binding of Sodium and Calcium to a POPC Bilayer by Effective Inclusion of Electronic PolarizationMelcr, Josef; Martinez-Seara, Hector; Nencini, Ricky; Kolafa, Jiri; Jungwirth, Pavel; Ollila, O. H. SamuliJournal of Physical Chemistry B (2018), 122 (16), 4546-4557CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Binding affinities and stoichiometries of Na+ and Ca2+ ions to phospholipid bilayers are of paramount significance in the properties and functionality of cellular membranes. Current ests. of binding affinities and stoichiometries of cations are, however, inconsistent due to limitations in the available exptl. and computational methods. In this work, the authors improve the description of the binding details of Na+ and Ca2+ ions to a 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) bilayer by implicitly including electronic polarization as a mean field correction, known as the electronic continuum correction (ECC). This is applied by scaling the partial charges of a selected state-of-the-art POPC lipid model for mol. dynamics simulations. The improved ECC-POPC model reproduces not only the exptl. measured structural parameters for the ion-free membrane, but also the response of lipid headgroup to a strongly bound cationic amphiphile, as well as the binding affinities of Na+ and Ca2+ ions. With the new model, the authors observe on the one side negligible binding of Na+ ions to POPC bilayer, while on the other side stronger interactions of Ca2+ primarily with phosphate oxygens, which is in agreement with the previous interpretations of the exptl. spectroscopic data. The present model results in Ca2+ ions forming complexes with one to three POPC mols. with almost equal probabilities, suggesting more complex binding stoichiometries than those from simple models used to interpret the NMR data previously. The results of this work pave the way to quant. mol. simulations with realistic electrostatic interactions of complex biochem. systems at cellular membranes.
- 51Melcrová, A.; Pokorna, S.; Pullanchery, S.; Kohagen, M.; Jurkiewicz, P.; Hof, M.; Jungwirth, P.; Cremer, P. S.; Cwiklik, L. The complex nature of calcium cation interactions with phospholipid bilayers. Sci. Rep. 2016, 6, 38035 DOI: 10.1038/srep38035Google Scholar51https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XitFWnsr7F&md5=a4be8d96f3217f88b2ab791199726e08The complex nature of calcium cation interactions with phospholipid bilayersMelcrova, Adela; Pokorna, Sarka; Pullanchery, Saranya; Kohagen, Miriam; Jurkiewicz, Piotr; Hof, Martin; Jungwirth, Pavel; Cremer, Paul S.; Cwiklik, LukaszScientific Reports (2016), 6 (), 38035CODEN: SRCEC3; ISSN:2045-2322. (Nature Publishing Group)Understanding interactions of calcium with lipid membranes at the mol. level is of great importance in light of their involvement in calcium signaling, assocn. of proteins with cellular membranes, and membrane fusion. We quantify these interactions in detail by employing a combination of spectroscopic methods with atomistic mol. dynamics simulations. Namely, time-resolved fluorescent spectroscopy of lipid vesicles and vibrational sum frequency spectroscopy of lipid monolayers are used to characterize local binding sites of calcium in zwitterionic and anionic model lipid assemblies, while dynamic light scattering and zeta potential measurements are employed for macroscopic characterization of lipid vesicles in calcium-contg. environments. To gain addnl. at.-level information, the expts. are complemented by mol. simulations that utilize an accurate force field for calcium ions with scaled charges effectively accounting for electronic polarization effects. We demonstrate that lipid membranes have substantial calcium-binding capacity, with several types of binding sites present. Significantly, the binding mode depends on calcium concn. with important implications for calcium buffering, synaptic plasticity, and protein-membrane assocn.
- 52Burnstock, G. Purinergic signalling and disorders of the central nervous system. Nat. Rev. Drug Discovery 2008, 7, 575– 590, DOI: 10.1038/nrd2605Google Scholar52https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXnvFShtrs%253D&md5=95a3c406c113eec093fd5437e98ee49ePurinergic signaling and disorders of the central nervous systemBurnstock, GeoffreyNature Reviews Drug Discovery (2008), 7 (7), 575-590CODEN: NRDDAG; ISSN:1474-1776. (Nature Publishing Group)A review. There is increasing interest in the role of purinergic signaling in CNS disorders, but few modulators have reached the clinic. Here, Burnstock provides an overview of the role of purinergic signaling in specific disorders, including brain trauma, ischemia, neurodegenerative, neuroinflammatory and neuropsychiatric diseases, and highlights specific purinergic receptor subtypes that represent promising therapeutic targets. Purines have key roles in neurotransmission and neuromodulation, with their effects being mediated by the purine and pyrimidine receptor subfamilies, P1, P2X and P2Y. Recently, purinergic mechanisms and specific receptor subtypes have been shown to be involved in various pathol. conditions including brain trauma and ischemia, neurodegenerative diseases involving neuroimmune and neuroinflammatory reactions, as well as in neuropsychiatric diseases, including depression and schizophrenia. This article reviews the role of purinergic signaling in CNS disorders, highlighting specific purinergic receptor subtypes, most notably A2A, P2X4 and P2X7, that might be therapeutically targeted for the treatment of these conditions.
- 53Pomorski, T.; Menon, A. K. Lipid flippases and their biological functions. Cell. Mol. Life Sci. 2006, 63, 2908– 2921, DOI: 10.1007/s00018-006-6167-7Google Scholar53https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXht1egtb0%253D&md5=10e9fa64e238c152ee43de250d62a07eLipid flippases and their biological functionsPomorski, T.; Menon, A. K.Cellular and Molecular Life Sciences (2006), 63 (24), 2908-2921CODEN: CMLSFI; ISSN:1420-682X. (Birkhaeuser Verlag)A review. The typically distinct phospholipid compn. of the 2 leaflets of a membrane bilayer is generated and maintained by bi-directional transport (flip-flop) of lipids between the leaflets. Specific membrane proteins, termed lipid flippases, play an essential role in this transport process. Energy-independent flippases allow common phospholipids to equilibrate rapidly between the two monolayers and also play a role in the biosynthesis of a variety of glycoconjugates such as glycosphingolipids, N-glycoproteins, and glycosylphosphatidylinositol (GPI)-anchored proteins. ATP-dependent flippases, including members of a conserved subfamily of P-type ATPases and ATP-binding cassette transporters, mediate the net transfer of specific phospholipids to one leaflet of a membrane and are involved in the creation and maintenance of transbilayer lipid asymmetry of membranes such as the plasma membrane of eukaryotes. Energy-dependent flippases also play a role in the biosynthesis of glycoconjugates such as bacterial lipopolysaccharide. Here, the authors summarize recent progress on the identification and characterization of the various flippases and the demonstration of their biol. functions.
- 54Jiang, T.; Yu, K.; Hartzell, H. C.; Tajkhorshid, E. Lipids and ions traverse the membrane by the same physical pathway in the nhTMEM16 scramblase. eLife 2017, 6, e28671 DOI: 10.7554/eLife.28671Google Scholar54https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXitlOit7nJ&md5=e124dd1a49738ae073b1abca3c95ef2eLipids and ions traverse the membrane by the same physical pathway in the nhTMEM16 scramblaseJiang, Tao; Yu, Kuai; Hartzell, H. Criss; Tajkhorshid, EmadeLife (2017), 6 (), e28671/1-e28671/26CODEN: ELIFA8; ISSN:2050-084X. (eLife Sciences Publications Ltd.)From bacteria to mammals, different phospholipid species are segregated between the inner and outer leaflets of the plasma membrane by ATP-dependent lipid transporters. Disruption of this asymmetry by ATP-independent phospholipid scrambling is important in cellular signaling, but its mechanism remains incompletely understood. Using MD simulations coupled with exptl. assays, we show that the surface hydrophilic transmembrane cavity exposed to the lipid bilayer on the fungal scramblase nhTMEM16 serves as the pathway for both lipid translocation and ion conduction across the membrane. Ca2+ binding stimulates its open conformation by altering the structure of transmembrane helixes that line the cavity. We have identified key amino acids necessary for phospholipid scrambling and validated the idea that ions permeate TMEM16 Cl- channels via a structurally homologous pathway by showing that mutation of two residues in the pore region of the TMEM16A Ca2+-activated Cl- channel convert it into a robust scramblase.
- 55Bethel, N. P.; Grabe, M. Atomistic insight into lipid translocation by a TMEM16 scramblase. Proc. Natl. Acad. Sci. U.S.A. 2016, 113, 14049– 14054, DOI: 10.1073/pnas.1607574113Google Scholar55https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhvV2gu7rI&md5=c73f9572ab589dba5374c160eca1b9b6Atomistic insight into lipid translocation by a TMEM16 scramblaseBethel, Neville P.; Grabe, MichaelProceedings of the National Academy of Sciences of the United States of America (2016), 113 (49), 14049-14054CODEN: PNASA6; ISSN:0027-8424. (National Academy of Sciences)The transmembrane protein 16 (TMEM16) family of membrane proteins includes both lipid scramblases and ion channels involved in olfaction, nociception, and blood coagulation. The crystal structure of the fungal Nectria haematococca TMEM16 (nhTMEM16) scramblase suggested a putative mechanism of lipid transport, whereby polar and charged lipid headgroups move through the low-dielec. environment of the membrane by traversing a hydrophilic groove on the membrane-spanning surface of the protein. Here, we use computational methods to explore the membrane-protein interactions involved in lipid scrambling. Fast, continuum membrane-bending calcns. reveal a global pattern of charged and hydrophobic surface residues that bends the membrane in a large-amplitude sinusoidal wave, resulting in bilayer thinning across the hydrophilic groove. Atomic simulations uncover two lipid headgroup-interaction sites flanking the groove. The cytoplasmic site nucleates headgroup-dipole stacking interactions that form a chain of lipid mols. that penetrate into the groove. In two instances, a cytoplasmic lipid interdigitates into this chain, crosses the bilayer, and enters the extracellular leaflet, and the reverse process happens twice as well. Continuum membrane-bending anal. carried out on homol. models of mammalian homologs shows that these family members also bend the membrane-even those that lack scramblase activity. Sequence alignments show that the lipid-interaction sites are conserved in many family members but less so in those with reduced scrambling ability. Our anal. provides insight into how large-scale membrane bending and protein chem. facilitate lipid permeation in the TMEM16 family, and we hypothesize that membrane interactions also affect ion permeation.
- 56Pedemonte, N.; Galietta, L. J. V. Structure and function of TMEM16 proteins (anoctamins). Physiol. Rev. 2014, 94, 419– 459, DOI: 10.1152/physrev.00039.2011Google Scholar56https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXptFygt70%253D&md5=9b0cf5bb181dc1ce259d3cab7e7de629Structure and function of TMEM16 proteins (anoctamins)Pedemonte, Nicoletta; Galietta, Luis J. V.Physiological Reviews (2014), 94 (2), 419-459CODEN: PHREA7; ISSN:0031-9333. (American Physiological Society)A review. TMEM16 proteins, also known as anoctamins, are involved in a variety of functions that include ion transport, phospholipid scrambling, and regulation of other membrane proteins. The first two members of the family, TMEM16A (anoctamin-1, AN01) and TMEM16B (anoctamin-2, AN02), function as Ca2+-activated Cl- channels (CaCCs), a type of ion channel that plays important functions such as transepithelial ion transport, smooth muscle contraction, olfaction, phototransduction, nociception, and control of neuronal excitability. Genetic ablation of TMEM16A in mice causes impairment of epithelial Cl- secretion, tracheal abnormalities, and block of gastrointestinal peristalsis. TMEM16A is directly regulated by cytosolic Ca2+ as well as indirectly by its interaction with calmodulin. Other members of the anoctamin family, such as TMEM16C, TMEM16D, TMEM16F, TMEM16G, and TMEM16J, may work as phospholipid scramblases and/or ion channels. In particular, TMEM16F (AN06) is a major contributor to the process of phosphatidylserine translocation from the inner to the outer leaflet of the plasma membrane. Intriguingly, TMEM16F is also assocd. with the appearance of anion/cation channels activated by very high Ca2+ concns. Furthermore, a TMEM16 protein expressed in Aspergillus fumigatus displays both ion channel and lipid scramblase activity. This finding suggests that dual function is an ancestral characteristic of TMEM16 proteins and that some members, such as TMEM16A and TMEM16B, have evolved to a pure channel function. Mutations in anoctamin genes (AN03, AN05, AN06, and AN010) cause various genetic diseases. These diseases suggest the involvement of anoctamins in a variety of cell functions whose link with ion transport and/or lipid scrambling needs to be clarified.
- 57Kalienkova, V.; Clerico Mosina, V.; Bryner, L.; Oostergetel, G. T.; Dutzler, R.; Paulino, C. Stepwise activation mechanism of the scramblase nhTMEM16 revealed by cryo-EM. eLife 2019, 8, e44364 DOI: 10.7554/eLife.44364Google Scholar57https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXisVegtb3K&md5=fd6a7da5a78a87f0439578ee967f76dbStepwise activation mechanism of the scramblase nhTMEM16 revealed by cryo- EMKalienkova, Valeria; Mosina, Vanessa Clerico; Bryner, Laura; Oostergetel, Gert T.; Dutzler, Raimund; Paulino, CristinaeLife (2019), 8 (), 44364CODEN: ELIFA8; ISSN:2050-084X. (eLife Sciences Publications Ltd.)Scramblases catalyze the movement of lipids between both leaflets of a bilayer. Whereas the X-ray structure of the protein nhTMEM16 has previously revealed the architecture of a Ca2+-dependent lipid scramblase, its regulation mechanism has remained elusive. Here, we have used cryo-electron microscopy and functional assays to address this question. Ca2+-bound and Ca2+-free conformations of nhTMEM16 in detergent and lipid nanodiscs illustrate the interactions with its environment and they reveal the conformational changes underlying its activation. In this process, Ca2+ binding induces a stepwise transition of the catalytic subunit cavity, converting a closed cavity that is shielded from the membrane in the absence of ligand, into a polar furrow that becomes accessible to lipid headgroups in the Ca2+-bound state. Addnl., our structures demonstrate how nhTMEM16 distorts the membrane at both entrances of the subunit cavity, thereby decreasing the energy barrier for lipid movement.
- 58Falzone, M. E.; Rheinberger, J.; Lee, B.-C.; Peyear, T.; Sasset, L.; Raczkowski, A. M.; Eng, E. T.; Di Lorenzo, A.; Andersen, O. S.; Nimigean, C. M.; Accardi, A. Structural basis of Ca2+-dependent activation and lipid transport by a TMEM16 scramblase. eLife 2019, 8, e43229 DOI: 10.7554/eLife.43229Google Scholar58https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXisVegtLvK&md5=bf21156b76417f26aa25c95a6203cfa2Structural basis of Ca2+-dependent activation and lipid transport by a TMEM16 scramblaseFalzone, Maria E.; Rheinberger, Jan; Lee, Byoung-Cheol; Peyear, Thasin; Sasset, Linda; Raczkowski, Ashleigh M.; Eng, Edward T.; Lorenzo, Annarita Di; Andersen, Olaf S.; Nimigean, Crina M.; Accardi, AlessioeLife (2019), 8 (), 43229CODEN: ELIFA8; ISSN:2050-084X. (eLife Sciences Publications Ltd.)The lipid distribution of plasma membranes of eukaryotic cells is asym. and phospholipid scramblases disrupt this asymmetry by mediating the rapid, nonselective transport of lipids down their concn. gradients. As a result, phosphatidylserine is exposed to the outer leaflet of membrane, an important step in extracellular signaling networks controlling processes such as apoptosis, blood coagulation, membrane fusion and repair. Several TMEM16 family members have been identified as Ca2+-activated scramblases, but the mechanisms underlying their Ca2+-dependent gating and their effects on the surrounding lipid bilayer remain poorly understood. Here, we describe three high-resoln. cryo-electron microscopy structures of a fungal scramblase from Aspergillus fumigatus, afTMEM16, reconstituted in lipid nanodiscs. These structures reveal that Ca2+-dependent activation of the scramblase entails global rearrangement of the transmembrane and cytosolic domains. These structures, together with functional expts., suggest that activation of the protein thins the membrane near the transport pathway to facilitate rapid transbilayer lipid movement.
- 59Feng, S.; Dang, S.; Han, T. W.; Ye, W.; Jin, P.; Cheng, T.; Li, J.; Jan, Y. N.; Jan, L. Y.; Cheng, Y. Cryo-EM studies of TMEM16F Calcium-Activated ion channel suggest features important for lipid scrambling. Cell Rep. 2019, 28, 567– 579, DOI: 10.1016/j.celrep.2019.06.023Google Scholar59https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtlGjtrjI&md5=be585b9b50fc27815a33b0aae7de9bdbCryo-EM Studies of TMEM16F Calcium-Activated Ion Channel Suggest Features Important for Lipid ScramblingFeng, Shengjie; Dang, Shangyu; Han, Tina Wei; Ye, Wenlei; Jin, Peng; Cheng, Tong; Li, Junrui; Jan, Yuh Nung; Jan, Lily Yeh; Cheng, YifanCell Reports (2019), 28 (2), 567-579.e4CODEN: CREED8; ISSN:2211-1247. (Cell Press)As a Ca2+-activated lipid scramblase and ion channel that mediates Ca2+ influx, TMEM16F relies on both functions to facilitate extracellular vesicle generation, blood coagulation, and bone formation. How a bona fide ion channel scrambles lipids remains elusive. Our structural analyses revealed the coexistence of an intact channel pore and PIP2-dependent protein conformation changes leading to membrane distortion. Correlated to the extent of membrane distortion, many tightly bound lipids are slanted. Structure-based mutagenesis studies further reveal that neutralization of some lipid-binding residues or those near membrane distortion specifically alters the onset of lipid scrambling, but not Ca2+ influx, thus identifying features outside of channel pore that are important for lipid scrambling. Together, our studies demonstrate that membrane distortion does not require open hydrophilic grooves facing the membrane interior and provide further evidence to suggest sep. pathways for lipid scrambling and ion permeation.
- 60Alvadia, C.; Lim, N. K.; Clerico Mosina, V.; Oostergetel, G. T.; Dutzler, R.; Paulino, C. Cryo-EM structures and functional characterization of the murine lipid scramblase TMEM16F. eLife 2019, 8, e44365 DOI: 10.7554/eLife.44365Google Scholar60https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXisVegt7rL&md5=0cd86dbb69ebf0ad0da6c047b9efc032Cryo-EM structures and functional characterization of the murine lipid scramblase TMEM16FAlvadia, Carolina; Lim, Novandy K.; Mosina, Vanessa Clerico; Oostergetel, Gert T.; Dutzler, Raimund; Paulino, CristinaeLife (2019), 8 (), 44365CODEN: ELIFA8; ISSN:2050-084X. (eLife Sciences Publications Ltd.)The lipid scramblase TMEM16F initiates blood coagulation by catalyzing the exposure of phosphatidylserine in platelets. The protein is part of a family of membrane proteins, which encompasses calcium-activated channels for ions and lipids. Here, we reveal features of murine TMEM16F (mTMEM16F) that underlie its function as a lipid scramblase and an ion channel. The cryo-EM data of mTMEM16F in absence and presence of Ca2+ define the ligand-free closed conformation of the protein and the structure of a Ca2+-bound intermediate. Both conformations resemble their counterparts of the scrambling-incompetent anion channel mTMEM16A, yet with distinct differences in the region of ion and lipid permeation. In conjunction with functional data, we demonstrate the relationship between ion conduction and lipid scrambling. Although activated by a common mechanism, both functions appear to be mediated by alternate protein conformations that are at equil. in the ligand-bound state.
- 61Lim, N. K.; Lam, A. K.; Dutzler, R. Independent activation of ion conduction pores in the double-barreled calcium-activated chloride channel TMEM16A. J. Gen. Physiol. 2016, 148, 375– 392, DOI: 10.1085/jgp.201611650Google Scholar61https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXkvF2qtLk%253D&md5=4cab74ef8af397f48127717903dfcf78Independent activation of ion conduction pores in the double-barreled calcium-activated chloride channel TMEM16ALim, Novandy K.; Lam, Andy K. M.; Dutzler, RaimundJournal of General Physiology (2016), 148 (5), 375-392, S19-S20CODEN: JGPLAD; ISSN:1540-7748. (Rockefeller University Press)The TMEM16 proteins constitute a family of membrane proteins with unusual functional breadth, including lipid scramblases and Cl- channels. Members of both these branches are activated by Ca2+, acting from the intracellular side, and probably share a common architecture, which was defined in the recent structure of the lipid scramblase nhTMEM16. The structural features of subunits and the arrangement of Ca2+-binding sites in nhTMEM16 suggest that the dimeric protein harbors two locations for catalysis that are independent with respect to both activation and lipid conduction. Here, we ask whether a similar independence is obsd. in the Ca2+-activated Cl- channel TMEM16A. For this purpose, we generated concatenated constructs contg. subunits with distinct activation and permeation properties. Our biochem. investigations demonstrate the integrity of concatemers after solubilization and purifn. During investigation by patch-clamp electrophysiol., the functional behavior of constructs contg. either two wild-type (WT) subunits or one WT subunit paired with a second subunit with compromised activation closely resembles TMEM16A. This resemblance extends to ion selectivity, conductance, and the concn. and voltage dependence of channel activation by Ca2+. Constructs combining subunits with different potencies for Ca2+ show a biphasic activation curve that can be described as a linear combination of the properties of its constituents. The functional independence is further supported by mutation of a putative pore-lining residue that changes the conduction properties of the mutated subunit. Our results strongly suggest that TMEM16A contains two ion conduction pores that are independently activated by Ca2+ binding to sites that are embedded within the transmembrane part of each subunit.
- 62Jeng, G.; Aggarwal, M.; Yu, W.-P.; Chen, T.-Y. Independent activation of distinct pores in dimeric TMEM16A channels. J. Gen. Physiol. 2016, 148, 393– 404, DOI: 10.1085/jgp.201611651Google Scholar62https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXkvF2qtrk%253D&md5=72fcb349ae446a73a14d16604751294aIndependent activation of distinct pores in dimeric TMEM16A channelsJeng, Grace; Aggarwal, Muskaan; Yu, Wei-Ping; Chen, Tsung-YuJournal of General Physiology (2016), 148 (5), 393-404, S13-S14CODEN: JGPLAD; ISSN:1540-7748. (Rockefeller University Press)The TMEM16 family encompasses Ca2+-activated Cl- channels (CaCCs) and lipid scramblases. These proteins are formed by two identical subunits, as confirmed by the recently solved crystal structure of a TMEM16 lipid scram-blase. However, the high-resoln. structure did not provide definitive information regarding the pore architecture of the TMEM16 channels. In this study, we express TMEM16A channels constituting two covalently linked subunits with different Ca2+ affinities. The dose-response curve of the heterodimer appears to be a weighted sum of two dose-response curves-one corresponding to the high-affinity subunit and the other to the low-affinity subunit. However, fluorescence resonance energy transfer expts. suggest that the covalently linked heterodimeric proteins fold and assemble as one mol. Together these results suggest that activation of the two TMEM16A subunits likely activate independently of each other. The Ca2+ activation curve for the heterodi-mer at a low Ca2+ concn. range ([Ca2+] < 5μM) is similar to that of the wild-type channel-the Hill coeffs. in both cases are significantly greater than one. This suggests that Ca2+ binding to one subunit of TMEM16A is sufficient to activate the channel and that each subunit contains more than one Ca2+-binding site. We also take advantage of the I-V curve rectification that results from mutation of a pore residue to address the pore architec-ture of the channel. By introducing the pore mutation and the mutation that alters Ca2+ affinity in the same or different subunits, we demonstrate that activation of different subunits appears to be assocd. with the opening of different pores. These results suggest that the TMEM16A CaCC may also adopt a"double-barrel"pore architecture, similar to that found in CLC channels and transporters.
- 63Kanner, B. I.; Bendahan, A. Binding order of substrates to the (sodium + potassium ion)-coupled L-glutamic acid transporter from rat brain. Biochemistry 1982, 21, 6327– 6330, DOI: 10.1021/bi00267a044Google Scholar63https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL38XlvFelsLo%253D&md5=2ff713407683ab142203cc789e4a3b61Binding order of substrates to the (sodium + potassium ion)-coupled L-glutamic acid transporter from rat brainKanner, Baruch I.; Bendahan, AnnieBiochemistry (1982), 21 (24), 6327-30CODEN: BICHAW; ISSN:0006-2960.Efflux of L-glutamic acid (I) [7440-23-5] from rat brain synaptic plasm membrane vesicles requires external K+. This requirement is satd. by concns. of ∼15 mequiv/L K+. In the absence of K+, I can be released from the vesicles in the presence of external I. This stimulation does not require external Na+ but is dependent on the external concn. of I. Half-maximal effects are obtained by concns. of ∼1 μM which are very similar to the apparent Km for I influx. Efflux of labeled I driven by external Na+ plus I requires internal Na+. The transporter apparently displays an asym. behavior toward Na+. This ion dissocs. much more slowly than does I on the external surface of the membrane but not on the internal surface. Furthermore, the transporter apparently translocates K+ in a step distinct from the I translocation step. The simplest explanation is that upon translocation of Na+ and I and their release to the inside, K+ binds to the transporter, enabling it to return to the outside to allow initiation of a new transport cycle.
- 64Zerangue, N.; Kavanaugh, M. P. Flux coupling in a neuronal glutamate transporter. Nature 1996, 383, 634– 637, DOI: 10.1038/383634a0Google Scholar64https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XmtlSqsbg%253D&md5=4d4328ffafe70486746bf0b8d3929a4fFlux coupling in a neuronal glutamate transporterZerangue, Noa; Kavanaugh, Michael P.Nature (London) (1996), 383 (6601), 634-637CODEN: NATUAS; ISSN:0028-0836. (Macmillan Magazines)Synaptic transmission is commonly terminated by diffusion and reuptake of neurotransmitter from the synaptic cleft. Glutamate reuptake prevents neurotoxicity and sets the lower limit for the concn. of extracellular glutamate, so it is important to understand the thermodn. of this process. Here we use voltage clamping with a pH-sensitive fluorescent dye to monitor elec. currents and pH changes assocd. with flux of glutamate mediated by the human neuronal glutamate transporter EAAT3. In contrast to a previous model, we find that three sodium ions and one proton are cotransported with each glutamate ion into the cell, while one potassium ion is transported out of the cell. This coupling can support a transmembrane glutamate concn. gradient ([Glu]in/[Glu]out) exceeding 106 under equil. conditions, and would allow the transporter to continue removing glutamate over a wide range of ionic conditions.
- 65Vandenberg, R. J.; Ryan, R. M. Mechanisms of Glutamate Transport. Physiol. Rev. 2013, 93, 1621– 1657, DOI: 10.1152/physrev.00007.2013Google Scholar65https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXhvVShtrrP&md5=b573b96bbaf6f03c835c471cf2aaf7b8Mechanisms of glutamate transportVandenberg, Robert J.; Ryan, Renae M.Physiological Reviews (2013), 93 (4), 1621-1657CODEN: PHREA7; ISSN:0031-9333. (American Physiological Society)A review. L-Glutamate is the predominant excitatory neurotransmitter in the mammalian central nervous system and plays important roles in a wide variety of brain functions, but it is also a key player in the pathogenesis of many neurol. disorders. The control of glutamate concns. is crit. to the normal functioning of the central nervous system, and in this review we discuss how glutamate transporters regulate glutamate concns. to maintain dynamic signaling mechanisms between neurons. In 2004, the crystal structure of a prokaryotic homolog of the mammalian glutamate transporter family of proteins was crystd. and its structure detd. This has paved the way for a better understanding of the structural basis for glutamate transporter function. In this review we provide a broad perspective of this field of research, but focus primarily on the more recent studies with a particular emphasis on how our understanding of the structure of glutamate transporters has generated new insights.
- 66Machtens, J.-P.; Kortzak, D.; Lansche, C.; Leinenweber, A.; Kilian, P.; Begemann, B.; Zachariae, U.; Ewers, D.; de Groot, B. L.; Briones, R.; Fahlke, C. Mechanisms of anion conduction by coupled glutamate transporters. Cell 2015, 160, 542– 553, DOI: 10.1016/j.cell.2014.12.035Google Scholar66https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXitFyls7c%253D&md5=36f8d3049997b4b1877e4dbc13b6630fMechanisms of anion conduction by coupled glutamate transportersMachtens, Jan-Philipp; Kortzak, Daniel; Lansche, Christine; Leinenweber, Ariane; Kilian, Petra; Begemann, Birgit; Zachariae, Ulrich; Ewers, David; de Groot, Bert L.; Briones, Rodolfo; Fahlke, ChristophCell (Cambridge, MA, United States) (2015), 160 (3), 542-553CODEN: CELLB5; ISSN:0092-8674. (Cell Press)Excitatory amino acid transporters (EAATs) are essential for terminating glutamatergic synaptic transmission. They are not only coupled glutamate/Na+/H+/K+ transporters, but they also function as anion-selective channels. EAAT anion channels regulate neuronal excitability, and gain-of-function mutations in these proteins result in ataxia and epilepsy. Here, the authors combined mol. dynamics simulations with fluorescence spectroscopy of prokaryotic homolog glutamate transporter GltPh and patch-clamp recordings of mammalian EAATs to det. how these transporters conduct anions. Whereas outward- and inward-facing GltPh conformations were nonconductive, lateral movement of the glutamate transport domain from intermediate transporter conformations resulted in the formation of an anion-selective conduction pathway. Fluorescence quenching of inserted Trp residues indicated the entry of anions into this pathway, and mutations of homologous pore-forming residues had analogous effects on GltPh simulations and EAAT2/EAAT4 measurements of single-channel currents and anion/cation selectivities. These findings provided a mechanistic framework of how neurotransmitter transporters can operate as anion-selective and ligand-gated ion channels.
- 67Bastug, T.; Heinzelmann, G.; Kuyucak, S.; Salim, M.; Vandenberg, R. J.; Ryan, R. M. Position of the Third Na+ Site in the Aspartate Transporter GltPh and the Human Glutamate Transporter, EAAT1. PLoS One 2012, 7, e33058 DOI: 10.1371/journal.pone.0033058Google Scholar67https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38Xks1Grs7Y%253D&md5=5b2a68466812e3ca7c2730b78dc4e0c0Position of the third Na+ site in the Aspartate transporter GltPh and the human glutamate transporter, EAAT1Bastug, Turgut; Heinzelmann, Germano; Kuyucak, Serdar; Salim, Marietta; Vandenberg, Robert J.; Ryan, Renae M.PLoS One (2012), 7 (3), e33058CODEN: POLNCL; ISSN:1932-6203. (Public Library of Science)Glutamate transport via the human excitatory amino acid transporters is coupled to the co-transport of three Na+ ions, one H+ and the counter-transport of one K+ ion. Transport by an archaeal homolog of the human glutamate transporters, GltPh, whose three dimensional structure is known is also coupled to three Na+ ions but only two Na+ ion binding sites have been obsd. in the crystal structure of GltPh. In order to fully utilize the GltPh structure in functional studies of the human glutamate transporters, it is essential to understand the transport mechanism of GltPh and accurately det. the no. and location of Na+ ions coupled to transport. Several sites have been proposed for the binding of a third Na+ ion from electrostatic calcns. and mol. dynamics simulations. In this study, we have performed detailed free energy simulations for GltPh and reveal a new site for the third Na+ ion involving the side chains of Threonine 92, Serine 93, Asparagine 310, Aspartate 312, and the backbone of Tyrosine 89. We have also studied the transport properties of alanine mutants of the coordinating residues Threonine 92 and Serine 93 in GltPh, and the corresponding residues in a human glutamate transporter, EAAT1. The mutant transporters have reduced affinity for Na+ compared to their wild type counterparts. These results confirm that Threonine 92 and Serine 93 are involved in the coordination of the third Na+ ion in GltPh and EAAT1.
- 68Kortzak, D.; Alleva, C.; Weyand, I.; Ewers, D.; Zimmermann, M. I.; Franzen, A.; Machtens, J.-P.; Fahlke, C. Allosteric gate modulation confers K+ coupling in glutamate transporters. EMBO J. 2019, 38, e101468 DOI: 10.15252/embj.2019101468Google ScholarThere is no corresponding record for this reference.
- 69Setiadi, J.; Kuyucak, S. Free-Energy Simulations Resolve the Low-Affinity Na+-High-Affinity Asp Binding Paradox in GltPh. Biophys. J. 2019, 117, 780– 789, DOI: 10.1016/j.bpj.2019.07.016Google Scholar69https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhsFSrs7nO&md5=730d8e0bc640faf763c5ebc5e004bcffFree-Energy Simulations Resolve the Low-Affinity Na+-High-Affinity Asp Binding Paradox in GltPhSetiadi, Jeffry; Kuyucak, SerdarBiophysical Journal (2019), 117 (4), 780-789CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)Glutamate transporters clear up excess extracellular glutamate by cotransporting three Na+ and one H+ with the countertransport of one K+. The archaeal homologs are selective to aspartate and only cotransport three Na+. The crystal structures of GltPh from archaea have been used in computational studies to understand the transport mechanism. Although some progress has been made with regard to the ligand-binding sites, a consistent picture of transport still eludes us. A major concern is the discrepancy between the computed binding free energies, which predict high-affinity Na+-low-affinity aspartate binding, and the exptl. results in which the opposite is obsd. Here, we show that the binding of the first two Na+ ions involves an intermediate state near the Na1 site, where two Na+ ions coexist and couple to aspartate with similar strengths, boosting its affinity. Binding free energies for Na+ and aspartate obtained using this intermediate state are in good agreement with the exptl. values. Thus, the paradox in binding affinities arises from the assumption that the ligands bind to the sites obsd. in the crystal structure following the order dictated by their binding free energies with no intermediate states. In fact, the presence of an intermediate state eliminates such a correlation between the binding free energies and the binding order. The intermediate state also facilitates transition of the first Na+ ion to its final binding site via a knock-on mechanism, which induces substantial conformational changes in the protein consistent with exptl. observations.
- 70Machtens, J.-P.; Briones, R.; Alleva, C.; de Groot, B. L.; Fahlke, C. Gating Charge Calculations by Computational Electrophysiology Simulations. Biophys. J. 2017, 112, 1396– 1405, DOI: 10.1016/j.bpj.2017.02.016Google Scholar70https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXkt1Sqs7c%253D&md5=416dc846d784f5ce514486d29d376372Gating Charge Calculations by Computational Electrophysiology SimulationsMachtens, Jan-Philipp; Briones, Rodolfo; Alleva, Claudia; de Groot, Bert L.; Fahlke, ChristophBiophysical Journal (2017), 112 (7), 1396-1405CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)Elec. cell signaling requires adjustment of ion channel, receptor, or transporter function in response to changes in membrane potential. For the majority of such membrane proteins, the mol. details of voltage sensing remain insufficiently understood. Here, we present a mol. dynamics simulation-based method to det. the underlying charge movement across the membrane-the gating charge-by measuring elec. capacitor properties of membrane-embedded proteins. We illustrate the approach by calcg. the charge transfer upon membrane insertion of the HIV gp41 fusion peptide, and validate the method on two prototypical voltage-dependent proteins, the Kv1.2 K+ channel and the voltage sensor of the Ciona intestinalis voltage-sensitive phosphatase, against exptl. data. We then use the gating charge anal. to study how the T1 domain modifies voltage sensing in Kv1.2 channels and to investigate the voltage dependence of the initial binding of two Na+ ions in Na+-coupled glutamate transporters. Our simulation approach quantifies various mechanisms of voltage sensing, enables direct comparison with expts., and supports mechanistic interpretation of voltage sensitivity by fractional amino acid contributions.
- 71Groeneveld, M.; Slotboom, D.-J. Na+:Aspartate Coupling Stoichiometry in the Glutamate Transporter Homologue GltPh. Biochemistry 2010, 49, 3511– 3513, DOI: 10.1021/bi100430sGoogle Scholar71https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXksVantLY%253D&md5=45ca89490d62bfb11bfb2fde29c0db46Na+:Aspartate Coupling Stoichiometry in the Glutamate Transporter Homologue GltPhGroeneveld, Maarten; Slotboom, Dirk-JanBiochemistry (2010), 49 (17), 3511-3513CODEN: BICHAW; ISSN:0006-2960. (American Chemical Society)The Na+ aspartate symporter GltPh from Pyrococcus horikoshii is the only member of the glutamate transporter family for which crystal structures have been detd. The cation:aspartate coupling stoichiometry is unknown, thus hampering the elucidation of the ion coupling mechanism. Here we measure transport of 22Na+ and [14C]aspartate in proteoliposomes contg. purified GltPh and demonstrate that three Na+ ions are symported with aspartate.
- 72Samal, S. K.; Dash, M.; Van Vlierberghe, S.; Kaplan, D. L.; Chiellini, E.; van Blitterswijk, C.; Moroni, L.; Dubruel, P. Cationic polymers and their therapeutic potential. Chem. Soc. Rev. 2012, 41, 7147– 7194, DOI: 10.1039/c2cs35094gGoogle Scholar72https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XhsVOnurfE&md5=418983b4d598ad31e8010d09cc1b60c5Cationic polymers and their therapeutic potentialSamal, Sangram Keshari; Dash, Mamoni; Van Vlierberghe, Sandra; Kaplan, David L.; Chiellini, Emo; van Blitterswijk, Clemens; Moroni, Lorenzo; Dubruel, PeterChemical Society Reviews (2012), 41 (21), 7147-7194CODEN: CSRVBR; ISSN:0306-0012. (Royal Society of Chemistry)A review. The last decade has witnessed enormous research focused on cationic polymers. Cationic polymers are the subject of intense research as non-viral gene delivery systems, due to their flexible properties, facile synthesis, robustness and proven gene delivery efficiency. Here, the authors review the most recent scientific advances in cationic polymers and their derivs. not only for gene delivery purposes but also for various alternative therapeutic applications. An overview of the synthesis and prepn. of cationic polymers is provided along with their inherent bioactive and intrinsic therapeutic potential. In addn., cationic polymer based biomedical materials are covered. Major progress in the fields of drug and gene delivery as well as tissue engineering applications is summarized in the present review.
- 73Kostritskii, A. Y.; Kondinskaia, D. A.; Nesterenko, A. M.; Gurtovenko, A. A. Adsorption of Synthetic Cationic Polymers on Model Phospholipid Membranes: Insight from Atomic-Scale Molecular Dynamics Simulations. Langmuir 2016, 32, 10402– 10414, DOI: 10.1021/acs.langmuir.6b02593Google Scholar73https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhsFens7bP&md5=1d3d471cd08d1aaa4c10582314cbe459Adsorption of Synthetic Cationic Polymers on Model Phospholipid Membranes: Insight from Atomic-Scale Molecular Dynamics SimulationsKostritskii, Andrei Yu.; Kondinskaia, Diana A.; Nesterenko, Alexey M.; Gurtovenko, Andrey A.Langmuir (2016), 32 (40), 10402-10414CODEN: LANGD5; ISSN:0743-7463. (American Chemical Society)Although synthetic cationic polymers represent a promising class of effective antibacterial agents, the mol. mechanisms behind their antimicrobial activity remain poorly understood. To this end, we employ at.-scale mol. dynamics simulations to explore adsorption of several linear cationic polymers of different chem. structure and protonation (polyallylamine (PAA), polyethylenimine (PEI), polyvinylamine (PVA) and poly-L-lysine (PLL)) on model bacterial membranes (4:1 mixt. of zwitterionic phosphatidylethanolamine (PE) and anionic phosphatidylglycerol (PG) lipids). Overall, our findings show that binding of polycations to the anionic membrane surface effectively neutralizes its charge, leading to the reorientation of water mols. close to the lipid/water interface and to the partial release of counterions to the water phase. In certain cases one has even an overcharging of the membrane, which was shown to be a cooperative effect of polymer charges and lipid counterions. Protonated amine groups of polycations are found to interact preferably with head groups of anionic lipids, giving rise to formation of hydrogen bonds and to a noticeable lateral immobilization of the lipids. While all the above findings are mostly defined by the overall charge of a polymer, we found that the polymer architecture also matters. In particular, PVA and PEI are able to accumulate anionic PG lipids on the membrane surface, leading to lipid segregation. In turn, PLL whose charge twice exceeds charges of PVA/PEI does not induce such lipid segregation due to its considerably less compact architecture and relatively long side chains. We also show that partitioning of a polycation into the lipid/water interface is an interplay between its protonation level (the overall charge) and hydrophobicity of the backbone. Therefore, a possible strategy in creating highly efficient antimicrobial polymeric agents could be in tuning these polycation's properties through proper combination of protonated and hydrophobic blocks.
- 74Kondinskaia, D. A.; Gurtovenko, A. A. Supramolecular complexes of DNA with cationic polymers: The effect of polymer concentration. Polymer 2018, 142, 277– 284, DOI: 10.1016/j.polymer.2018.03.048Google Scholar74https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXmslKlt7o%253D&md5=b803292ea83e9ebd5c210446b17035b6Supramolecular complexes of DNA with cationic polymers: The effect of polymer concentrationKondinskaia, Diana A.; Gurtovenko, Andrey A.Polymer (2018), 142 (), 277-284CODEN: POLMAG; ISSN:0032-3861. (Elsevier Ltd.)Supramol. complexes of DNA with cationic polymers are of tremendous importance as these polymers can successfully be used in gene therapy as delivery vectors of genetic material. In this work we employ atomistic mol. dynamics simulations to get an insight into the impact of polymer concn. on the structural and electrostatic properties of the complexes of DNA and cationic polymers. Four linear cationic polymers of different chem. structure and protonation state were studied: polyethylenimine (PEI), poly-L-lysine (PLL), polyvinylamine (PVA), and polyallylamine (PAA). For all considered polymers our computational findings clearly demonstrate that increasing the polymer concn. leads to the overcharging of the DNA mol. This is relevant as the overall pos. charge of the DNA-polycation complex is considered to be a prerequisite for efficient binding of the polyplexes to neg. charged cell membranes. However, the concn. effects themselves are found to be sensitive to the chem. structure of a polymer. In particular, flexible PEI chains, being characterized by the high affinity to DNA's phosphate groups, are able to bind to DNA in an independent manner. As a result, DNA and PEIs form compact complexes with the largest cumulative pos. charge among all four types of cationic polymers under study. In contrast, PVA chains show the affinity to the major groove of DNA due to the hydrophobic nature of their backbones. This leads to stronger interactions with the DNA mols. and to the competition between PVA chains for DNA binding sites. This competition is responsible for the obsd. satn. in the PVA-induced neutralizing of the DNA charges when the PVA concn. increases. As for PLL and PAA, both these polycations are found to be less effective in neutralizing the charge of the polyanionic DNA mols. as compared to PEI and PVA. PLL has relatively long side chains, so that some of them cannot access DNA phosphates and reside in aq. soln. In turn, PAA has a low protonation level, leading to a weak electrostatic binding to DNA. Overall, our findings can relate the properties of supramol. DNA-polycation complexes at elevated polymer concns. with the chem. structure of a polycation and therefore can be useful for the development of novel, highly efficient polycation-based delivery vectors of DNA.
- 75Rohs, R.; West, S. M.; Sosinsky, A.; Liu, P.; Mann, R. S.; Honig, B. The role of DNA shape in protein-DNA recognition. Nature 2009, 461, 1248– 1253, DOI: 10.1038/nature08473Google Scholar75https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhtlOjtrvL&md5=2c314c36cc47209d7869e3fbda04f5e5The role of DNA shape in protein-DNA recognitionRohs, Remo; West, Sean M.; Sosinsky, Alona; Liu, Peng; Mann, Richard S.; Honig, BarryNature (London, United Kingdom) (2009), 461 (7268), 1248-1253CODEN: NATUAS; ISSN:0028-0836. (Nature Publishing Group)The recognition of specific DNA sequences by proteins is thought to depend on 2 types of mechanism: one that involves the formation of H-bonds with specific bases, primarily in the major groove, and one involving sequence-dependent deformations of the DNA helix. Here, by comprehensively analyzing the 3-dimensional structures of protein-DNA complexes, the authors show that the binding of Arg residues to narrow minor grooves is a widely used mode for protein-DNA recognition. This readout mechanism exploits the phenomenon that narrow minor grooves strongly enhance the neg. electrostatic potential of the DNA. The nucleosome core particle offers a prominent example of this effect. Minor-groove narrowing is often assocd. with the presence of A-tracts, AT-rich sequences that exclude the flexible TpA step. These findings indicate that the ability to detect local variations in DNA shape and electrostatic potential is a general mechanism that enables proteins to use information in the minor groove, which otherwise offers few opportunities for the formation of base-specific H-bonds, to achieve DNA-binding specificity.
- 76Andrews, A. J.; Luger, K. Nucleosome structure(s) and stability: variations on a theme. Annu. Rev. Biophys. 2011, 40, 99– 117, DOI: 10.1146/annurev-biophys-042910-155329Google Scholar76https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXnvFaitbo%253D&md5=916cbdc89d17871e816f8e1b24926b90Nucleosome structure(s) and stability: variations on a themeAndrews, Andrew J.; Luger, KarolinAnnual Review of Biophysics (2011), 40 (), 99-117CODEN: ARBNCV; ISSN:1936-122X. (Annual Reviews Inc.)A review. Chromatin is a highly regulated, modular nucleoprotein complex that is central to many processes in eukaryotes. The organization of DNA into nucleosomes and higher-order structures has profound implications for DNA accessibility. Alternative structural states of the nucleosome, and the thermodn. parameters governing its assembly and disassembly, need to be considered in order to understand how access to nucleosomal DNA is regulated. Here, the authors provide a brief historical account of how the overriding perception regarding aspects of nucleosome structure has changed over the past 30 yr. The authors discuss recent tech. advances regarding nucleosome structure and its phys. characterization and review the evidence for alternative nucleosome conformations and their implications for nucleosome and chromatin dynamics.
- 77Hub, J. S.; de Groot, B. L. Detection of functional modes in protein dynamics. PLoS Comput. Biol. 2009, 5, e1000480, DOI: 10.1371/journal.pcbi.1000480Google Scholar77https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD1Mros1agug%253D%253D&md5=65120d75ba73b0d5fba27e45233f5ab2Detection of functional modes in protein dynamicsHub Jochen S; de Groot Bert LPLoS computational biology (2009), 5 (8), e1000480 ISSN:.Proteins frequently accomplish their biological function by collective atomic motions. Yet the identification of collective motions related to a specific protein function from, e.g., a molecular dynamics trajectory is often non-trivial. Here, we propose a novel technique termed "functional mode analysis" that aims to detect the collective motion that is directly related to a particular protein function. Based on an ensemble of structures, together with an arbitrary "functional quantity" that quantifies the functional state of the protein, the technique detects the collective motion that is maximally correlated to the functional quantity. The functional quantity could, e.g., correspond to a geometric, electrostatic, or chemical observable, or any other variable that is relevant to the function of the protein. In addition, the motion that displays the largest likelihood to induce a substantial change in the functional quantity is estimated from the given protein ensemble. Two different correlation measures are applied: first, the Pearson correlation coefficient that measures linear correlation only; and second, the mutual information that can assess any kind of interdependence. Detecting the maximally correlated motion allows one to derive a model for the functional state in terms of a single collective coordinate. The new approach is illustrated using a number of biomolecules, including a polyalanine-helix, T4 lysozyme, Trp-cage, and leucine-binding protein.
- 78Krivobokova, T.; Briones, R.; Hub, J. S.; Munk, A.; de Groot, B. L. Partial least-squares functional mode analysis: application to the membrane proteins AQP1, Aqy1, and CLC-ec1. Biophys. J. 2012, 103, 786– 796, DOI: 10.1016/j.bpj.2012.07.022Google Scholar78https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38Xht1ektLbN&md5=b40f9fbb32b434d8a8b71f0956e4dd9bPartial Least-Squares Functional Mode Analysis: Application to the Membrane Proteins AQP1, Aqy1, and CLC-ec1Krivobokova, Tatyana; Briones, Rodolfo; Hub, Jochen S.; Munk, Axel; de Groot, Bert L.Biophysical Journal (2012), 103 (4), 786-796CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)We introduce an approach based on the recently introduced functional mode anal. to identify collective modes of internal dynamics that maximally correlate to an external order parameter of functional interest. Input structural data can be either exptl. detd. structure ensembles or simulated ensembles, such as mol. dynamics trajectories. Partial least-squares regression is shown to yield a robust soln. to the multidimensional optimization problem, with a minimal and controllable risk of overfitting, as shown by extensive cross-validation. Several examples illustrate that the partial least-squares-based functional mode anal. successfully reveals the collective dynamics underlying the fluctuations in selected functional order parameters. Applications to T4 lysozyme, the Trp-cage, the aquaporin channels Aqy1 and hAQP1, and the CLC-ec1 chloride antiporter are presented in which the active site geometry, the hydrophobic solvent-accessible surface, channel gating dynamics, water permeability (pf), and a dihedral angle are defined as functional order parameters. The Aqy1 case reveals a gating mechanism that connects the inner channel gating residues with the protein surface, thereby providing an explanation of how the membrane may affect the channel. hAQP1 shows how the pf correlates with structural changes around the arom./arginine region of the pore. The CLC-ec1 application shows how local motions of the gating Glu148 couple to a collective motion that affects ion affinity in the pore.
- 79Bratanov, D.; Kovalev, K.; Machtens, J.-P.; Astashkin, R.; Chizhov, I.; Soloviov, D.; Volkov, D.; Polovinkin, V.; Zabelskii, D.; Mager, T.; Gushchin, I.; Rokitskaya, T.; Antonenko, Y.; Alekseev, A.; Shevchenko, V.; Yutin, N.; Rosselli, R.; Baeken, C.; Borshchevskiy, V.; Bourenkov, G.; Popov, A.; Balandin, T.; Büldt, G.; Manstein, D. J.; Rodriguez-Valera, F.; Fahlke, C.; Bamberg, E.; Koonin, E.; Gordeliy, V. Unique structure and function of viral rhodopsins. Nat. Commun. 2019, 10, 4939 DOI: 10.1038/s41467-019-12718-0Google Scholar79https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3MjjtFChsg%253D%253D&md5=f68e6b3fcfb1157bdb7cd373eaf7d8aeUnique structure and function of viral rhodopsinsBratanov Dmitry; Kovalev Kirill; Volkov Dmytro; Polovinkin Vitaly; Zabelskii Dmitrii; Alekseev Alexey; Shevchenko Vitaly; Baeken Christian; Borshchevskiy Valentin; Balandin Taras; Gordeliy Valentin; Bratanov Dmitry; Kovalev Kirill; Volkov Dmytro; Zabelskii Dmitrii; Alekseev Alexey; Shevchenko Vitaly; Baeken Christian; Borshchevskiy Valentin; Balandin Taras; Gordeliy Valentin; Bratanov Dmitry; Kovalev Kirill; Astashkin Roman; Polovinkin Vitaly; Gordeliy Valentin; Kovalev Kirill; Astashkin Roman; Soloviov Dmytro; Volkov Dmytro; Zabelskii Dmitrii; Gushchin Ivan; Alekseev Alexey; Shevchenko Vitaly; Borshchevskiy Valentin; Buldt Georg; Fahlke Christoph; Bamberg Ernst; Gordeliy Valentin; Kovalev Kirill; Alekseev Alexey; Shevchenko Vitaly; Kovalev Kirill; Machtens Jan-Philipp; Fahlke Christoph; Machtens Jan-Philipp; Chizhov Igor; Manstein Dietmar J; Soloviov Dmytro; Soloviov Dmytro; Soloviov Dmytro; Mager Thomas; Bamberg Ernst; Rokitskaya Tatyana; Antonenko Yuri; Yutin Natalya; Koonin Eugene; Rosselli Riccardo; Rodriguez-Valera Francisco; Bourenkov Gleb; Popov AlexanderNature communications (2019), 10 (1), 4939 ISSN:.Recently, two groups of rhodopsin genes were identified in large double-stranded DNA viruses. The structure and function of viral rhodopsins are unknown. We present functional characterization and high-resolution structure of an Organic Lake Phycodnavirus rhodopsin II (OLPVRII) of group 2. It forms a pentamer, with a symmetrical, bottle-like central channel with the narrow vestibule in the cytoplasmic part covered by a ring of 5 arginines, whereas 5 phenylalanines form a hydrophobic barrier in its exit. The proton donor E42 is placed in the helix B. The structure is unique among the known rhodopsins. Structural and functional data and molecular dynamics suggest that OLPVRII might be a light-gated pentameric ion channel analogous to pentameric ligand-gated ion channels, however, future patch clamp experiments should prove this directly. The data shed light on a fundamentally distinct branch of rhodopsins and may contribute to the understanding of virus-host interactions in ecologically important marine protists.
Cited By
This article is cited by 4 publications.
- Piya Patra, Yi Qin Gao. Sequence-Specific Structural Features and Solvation Properties of Transcription Factor Binding DNA Motifs: Insights from Molecular Dynamics Simulation. The Journal of Physical Chemistry B 2022, 126
(45)
, 9187-9206. https://doi.org/10.1021/acs.jpcb.2c05749
- Andrei Y. Kostritskii, Jan-Philipp Machtens. Domain- and state-specific shape of the electric field tunes voltage sensing in voltage-gated sodium channels. Biophysical Journal 2023, 122
(10)
, 1807-1821. https://doi.org/10.1016/j.bpj.2023.04.013
- Yuxuan Zhuang, Colleen M. Noviello, Ryan E. Hibbs, Rebecca J. Howard, Erik Lindahl. Differential interactions of resting, activated, and desensitized states of the α7 nicotinic acetylcholine receptor with lipidic modulators. Proceedings of the National Academy of Sciences 2022, 119
(43)
https://doi.org/10.1073/pnas.2208081119
- Annika Etges, Nicole Hellmig, Gudrun Walenda, Bassam G. Haddad, Jan-Philipp Machtens, Thomas Morosan, Lars Christian Rump, Ute I. Scholl. A Novel Homozygous KLHL3 Mutation as a Cause of Autosomal Recessive Pseudohypoaldosteronism Type II Diagnosed Late in Life. Nephron 2022, 146
(4)
, 418-428. https://doi.org/10.1159/000521626
Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.
Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.
The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.
Recommended Articles
Abstract
Figure 1
Figure 1. Water-based scheme to calculate the electrostatic potential. First, the long-range part of the potential is calculated by solving the Poisson equation on the grid using the SPME method. Second, each water molecule is assigned a potential averaged over a sphere around its center of geometry, resulting in a distribution of water-molecule potentials. Third, the frame is fitted onto the reference structure and the grid of the final electrostatic-potential map is built around the biomolecule. Fourth, the water-molecule potentials are distributed over the nodes of the electrostatic-potential map, followed by normalization of the node potential based on the number of contributing water molecules. The final map is used for the quantitative analysis of electrostatics, e.g., to generate a potential profile.
Figure 2
Figure 2. Calculation of the electrostatic potential within the pore of the P2X3 cation channel. (a) (Left) Structure of the P2X3 receptor transmembrane domain without the ATP-binding ectodomain and embedded in a lipid membrane. The oxygen and hydrogen atoms of water molecules and the phosphorus atoms around the protein are shown as red, white, and orange spheres, respectively. (Middle) The water-occupancy map, contoured at an occupancy level of 0.1. (Right) Electrostatic-potential map, obtained with β = 20 nm–1 and contoured at −400 mV. The maps were calculated from the first 100 ns block of the 1 μs trajectory. (b) Electrostatic-potential profiles along the P2X3 pore, calculated from the first 100 ns block of the 1 μs trajectory using different values for the β parameter. (c) Average electrostatic-potential profile along the P2X3 pore based on the whole 1 μs trajectory (β = 20 nm–1). Whiskers indicate the standard deviation of the potential, as calculated for 10 separate 100 ns blocks.
Figure 3
Figure 3. Electrostatic-potential profiles along the proteolipidic pore of TMEM16 lipid scramblases. (a) Dimeric structure of nhTMEM16 in a lipid membrane. Protomers are colored white and gray; phosphorus atoms are shown in orange. (b,c) Side views of nhTMEM16 in the mixed POPC:POPS membrane (b) and of TMEM16K in the POPC membrane (c), together with the examples of the masked water-occupancy maps used to calculate the electrostatic-potential profiles. The maps were contoured at an occupancy level of 0.1. Phosphorus atoms of POPC and POPS are shown in orange and cyan, respectively. (d) Average electrostatic-potential profiles along the proteolipidic pore of TMEM16K and nhTMEM16 in pure POPC and mixed POPC:POPS membranes. The profiles were averaged over 1 μs trajectories. Whiskers indicate the standard deviation of the potential, as calculated for 10 separate 100 ns blocks.
Figure 4
Figure 4. Sodium binding increases the electrostatic potential at the aspartate-binding site of GltPh. (a) Structure of the aspartate-bound GltPh protomer in the outward-facing conformation shown together with different occupancy states of the binding sites. The aspartate molecule is shown in violet, and Na+ ions are shown as blue spheres. The violet cross in the apo panel indicates the position where the potential was measured. Membrane borders are shown schematically. (b) Time courses of electrostatic potential in the aspartate-binding site (violet) and Na+ occupancy of the Na1 site (blue) in one of the protomers upon spontaneous Na+ binding in an MD simulation initiated from the apo state. The time courses are shown as running averages with 50 ns windows. (c) Electrostatic potential in the aspartate-binding site in different states. Each data point is an average potential from a 50 ns block of a single protomer trajectory, whiskers indicate the 5th and 95th percentiles, and a median value is shown for each state.
Figure 5
Figure 5. Local effects of polycations on DNA electrostatics. (a) Structure of DNA in complex with cationic polymers PVA (left) and PAA (right). The backbone and grooves of DNA are shown in red and white, respectively. Polymers are shown as spheres with nitrogen atoms in blue and carbon atoms in teal (PVA) and yellow (PAA); hydrogen atoms are omitted. (b) Radial distribution of electrostatic potential around the central axis of DNA either in free form or in complex with PVA or PAA, calculated either globally (top), or separately for major (middle) and minor (bottom) grooves. Whiskers indicate the standard deviation of the potential, calculated separately for 100 ns blocks.
References
This article references 79 other publications.
- 1Honig, B.; Nicholls, A. Classical electrostatics in biology and chemistry. Science 1995, 268, 1144– 1149, DOI: 10.1126/science.77618291https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXlvVyqsrY%253D&md5=19b5365db27b7d14c861b67fc79b38cfClassical electrostatics in biology and chemistryHonig, Barry; Nicholls, AnthonyScience (Washington, D. C.) (1995), 268 (5214), 1144-9CODEN: SCIEAS; ISSN:0036-8075. (American Association for the Advancement of Science)A review with 75 refs. A major revival in the use of classical electrostatics as an approach to the study of charged and polar mols. in aq. soln. has been made possible through the development of fast numerical and computational methods to solve the Poisson-Boltzmann equation for solute mols. that have complex shapes and charge distributions. Graphical visualization of the calcd. electrostatic potentials generated by proteins and nucleic acids has revealed insights into the role of electrostatic interactions in a wide range of biol. phenomena. Classical electrostatics has also proved to be a successful quant. tool yielding accurate descriptions of elec. potentials, diffusion-limited processes, pH-dependent properties of proteins, ionic strength-dependent phenonema, and the solvation free energies of org. mols.
- 2Li, C.; Jia, Z.; Chakravorty, A.; Pahari, S.; Peng, Y.; Basu, S.; Koirala, M.; Panday, S. K.; Petukh, M.; Li, L.; Alexov, E. DelPhi Suite: New Developments and Review of Functionalities. J. Comput. Chem. 2019, 40, 2502– 2508, DOI: 10.1002/jcc.260062https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXht1Cntb7J&md5=7b8a34361ca244d25c3e19b2a20bf86eDelPhi Suite: New Developments and Review of FunctionalitiesLi, Chuan; Jia, Zhe; Chakravorty, Arghya; Pahari, Swagata; Peng, Yunhui; Basu, Sankar; Koirala, Mahesh; Panday, Shailesh Kumar; Petukh, Marharyta; Li, Lin; Alexov, EmilJournal of Computational Chemistry (2019), 40 (28), 2502-2508CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)Electrostatic potential, energies, and forces affect virtually any process in mol. biol., however, computing these quantities is a difficult task due to irregularly shaped macromols. and the presence of water. Here, the authors report a new edition of the popular software package DelPhi along with describing its functionalities. The new DelPhi is a C++ object-oriented package supporting various levels of multiprocessing and memory distribution. It is demonstrated that multiprocessing results in significant improvement of computational time. Furthermore, for computations requiring large grid size (large macromol. assemblages), the approach of memory distribution is shown to reduce the requirement of RAM and thus permitting large-scale modeling to be done on Linux clusters with moderate architecture. The new release comes with new features, whose functionalities and applications are described as well.
- 3Baker, N. A.; Sept, D.; Joseph, S.; Holst, M. J.; McCammon, J. A. Electrostatics of nanosystems: application to microtubules and the ribosome. Proc. Natl. Acad. Sci. U.S.A. 2001, 98, 10037– 10041, DOI: 10.1073/pnas.1813423983https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD3MXmvFWisbc%253D&md5=1b861999ef12c6972e82e8ada0f387cbElectrostatics of nanosystems: application to microtubules and the ribosomeBaker, Nathan A.; Sept, David; Joseph, Simpson; Holst, Michael J.; McCammon, J. AndrewProceedings of the National Academy of Sciences of the United States of America (2001), 98 (18), 10037-10041CODEN: PNASA6; ISSN:0027-8424. (National Academy of Sciences)Evaluation of the electrostatic properties of biomols. has become a std. practice in mol. biophysics. Foremost among the models used to elucidate the electrostatic potential is the Poisson-Boltzmann equation; however, existing methods for solving this equation have limited the scope of accurate electrostatic calcns. to relatively small biomol. systems. Here we present the application of numerical methods to enable the trivially parallel soln. of the Poisson-Boltzmann equation for supramol. structures that are orders of magnitude larger in size. As a demonstration of this methodol., electrostatic potentials have been calcd. for large microtubule and ribosome structures. The results point to the likely role of electrostatics in a variety of activities of these structures.
- 4Callenberg, K. M.; Choudhary, O. P.; de Forest, G. L.; Gohara, D. W.; Baker, N. A.; Grabe, M. APBSmem: a graphical interface for electrostatic calculations at the membrane. PLoS One 2010, 5, e12722 DOI: 10.1371/journal.pone.0012722There is no corresponding record for this reference.
- 5Dong, F.; Olsen, B.; Baker, N. A. Computational methods for biomolecular electrostatics. Methods Cell Biol. 2008, 84, 843– 8705https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXht1amurw%253D&md5=865ba17d85d1ef7e60d5e222f28bf107Computational methods for biomolecular electrostaticsDong, Feng; Olsen, Brett; Baker, Nathan A.Methods in Cell Biology (2008), 84 (Biophysical Tools for Biologists, Volume One), 843-870CODEN: MCBLAG; ISSN:0091-679X. (Elsevier)A review. An understanding of intermol. interactions is essential for insight into how cells develop, operate, communicate, and control their activities. Such interactions include several components: contributions from linear, angular, and torsional forces in covalent bonds, van der waals forces, as well as electrostatics. Among the various components of mol. interactions, electrostatics are of special importance because of their long range and their influence on polar or charged mols., including water, aq. ions, and amino or nucleic acids, which are some of the primary components of living systems. Electrostatics, therefore, play important roles in detg. the structure, motion, and function of a wide range of biol. mols. This chapter presents a brief overview of electrostatic interactions in cellular systems, with a particular focus on how computational tools can be used to investigate these types of interactions.
- 6Wang, L. Measurements and implications of the membrane dipole potential. Annu. Rev. Biochem. 2012, 81, 615– 635, DOI: 10.1146/annurev-biochem-070110-1230336https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XhtVGlsL7M&md5=7eeb5e3fc8ad92e8b59ef4a039138b9cMeasurements and implications of the membrane dipole potentialWang, LiguoAnnual Review of Biochemistry (2012), 81 (), 615-635CODEN: ARBOAW; ISSN:0066-4154. (Annual Reviews Inc.)A review. There are three kinds of membrane potentials: the surface potentials, resulting from the accumulation of charges at the membrane surfaces; the transmembrane potential, detd. by imbalance of charge in the aq. solns.; and the dipole potential, a membrane-internal potential from the dipolar components of the phospholipids and interface water. The abs. value of the dipole potential has been very difficult to measure, although its value has been estd. to be at 200-1000 mV from ion translocation rates (detd. by the planar lipid bilayer method), the surface potential of lipid monolayers (detd. by the lipid monolayer method), mol.-dynamics calcns., and electron scattering using cryoelectron microscopy (cryo-EM). Spectroscopy methods have also been used to monitor the dipole potential changes on the basis of the obsd. fluorescence changes of voltage-sensitive probes. The dipole potential accounts for the much larger permeability of a bare phospholipid membrane to anions than cations and affects the conformation and function of membrane proteins.
- 7Warshel, A.; Sharma, P. K.; Kato, M.; Parson, W. W. Modeling electrostatic effects in proteins. Biochim. Biophys. Acta, Proteins Proteomics 2006, 1764, 1647– 1676, DOI: 10.1016/j.bbapap.2006.08.0077https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD28Xht1Sjs7%252FO&md5=78eb047d080cc02e209aa71d2ad278beModeling electrostatic effects in proteinsWarshel, Arieh; Sharma, Pankaz K.; Kato, Mitsunori; Parson, William W.Biochimica et Biophysica Acta, Proteins and Proteomics (2006), 1764 (11), 1647-1676CODEN: BBAPBW; ISSN:1570-9639. (Elsevier Ltd.)A review. Electrostatic energies provide what is perhaps the most effective tool for structure-function correlation of biol. mols. This review considers the current state of simulations of electrostatic energies in macromols. as well as the early developments of this field. We focus on the relationship between microscopic and macroscopic models, considering the convergence problems of the microscopic models and the fact that the dielec. consts. in semimacroscopic models depend on the definition and the specific treatment. The advances and the challenges in the field are illustrated considering a wide range of functional properties including pKa's, redox potentials, ion and proton channels, enzyme catalysis, ligand binding and protein stability. We conclude by pointing out that, despite the current problems and the significant misunderstandings in the field, there is an overall progress that should lead eventually to quant. descriptions of electrostatic effects in proteins and thus to quant. descriptions of the function of proteins.
- 8Cisneros, G. A.; Karttunen, M.; Ren, P.; Sagui, C. Classical electrostatics for biomolecular simulations. Chem. Rev. 2014, 114, 779– 814, DOI: 10.1021/cr300461d8https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXhtlGnurrF&md5=86c888e948cf518f01f2a93d66e1a7f7Classical Electrostatics for Biomolecular SimulationsCisneros, G. Andres; Karttunen, Mikko; Ren, Pengyu; Sagui, CelesteChemical Reviews (Washington, DC, United States) (2014), 114 (1), 779-814CODEN: CHREAY; ISSN:0009-2665. (American Chemical Society)In this review, the authors primarily treat electrostatic methods for classical biomol. simulations in explicit solvent, with special emphasis on the accuracy of the electrostatic representation. The first part of the article is dedicated to reviewing computational methods for dealing with the long-range nature of the electrostatic interactions, including traditional methods as well as recent extensions and developments. The second part of the article is dedicated to reviewing current methods for representing the mol. electronic charge cloud, that go well beyond the point-charge representation. In addn., the authors discuss various multiscale approaches at the quantum mechanics/mol. mechanics level and at the mol. mechanics/coarse-grain level. The general trends in both cases are toward more accuracy and proper treatment of electrostatics. The authors finish the article by describing other challenging developments and giving the perspective on the current state of the field.
- 9Hollingsworth, S. A.; Dror, R. O. Molecular Dynamics Simulation for All. Neuron 2018, 99, 1129– 1143, DOI: 10.1016/j.neuron.2018.08.0119https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXhslGltLzL&md5=35e14022763288ed45c2a605cc900f61Molecular Dynamics Simulation for AllHollingsworth, Scott A.; Dror, Ron O.Neuron (2018), 99 (6), 1129-1143CODEN: NERNET; ISSN:0896-6273. (Cell Press)A review. The impact of mol. dynamics (MD) simulations in mol. biol. and drug discovery has expanded dramatically in recent years. These simulations capture the behavior of proteins and other biomols. in full at. detail and at very fine temporal resoln. Major improvements in simulation speed, accuracy, and accessibility, together with the proliferation of exptl. structural data, have increased the appeal of biomol. simulation to experimentalists-a trend particularly noticeable in, although certainly not limited to, neuroscience. Simulations have proven valuable in deciphering functional mechanisms of proteins and other biomols., in uncovering the structural basis for disease, and in the design and optimization of small mols., peptides, and proteins. Here we describe, in practical terms, the types of information MD simulations can provide and the ways in which they typically motivate further exptl. work.
- 10Kutzner, C.; Köpfer, D. A.; Machtens, J.-P.; de Groot, B. L.; Song, C.; Zachariae, U. Insights into the function of ion channels by computational electrophysiology simulations. Biochim. Biophys. Acta, Biomembr. 2016, 1858, 1741– 1752, DOI: 10.1016/j.bbamem.2016.02.00610https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XisFOrsb4%253D&md5=20860bb3a64b0140a3957639bf7da039Insights into the function of ion channels by computational electrophysiology simulationsKutzner, Carsten; Koepfer, David A.; Machtens, Jan-Philipp; de Groot, Bert L.; Song, Chen; Zachariae, UlrichBiochimica et Biophysica Acta, Biomembranes (2016), 1858 (7_Part_B), 1741-1752CODEN: BBBMBS; ISSN:0005-2736. (Elsevier B.V.)A review. Ion channels are of universal importance for all cell types and play key roles in cellular physiol. and pathol. Increased insight into their functional mechanisms is crucial to enable drug design on this important class of membrane proteins, and to enhance our understanding of some of the fundamental features of cells. This review presents the concepts behind the recently developed simulation protocol Computational Electrophysiol. (CompEL), which facilitates the atomistic simulation of ion channels in action. In addn., the review provides guidelines for its application in conjunction with the mol. dynamics software package GROMACS. We first lay out the rationale for designing CompEL as a method that models the driving force for ion permeation through channels the way it is established in cells, i.e., by electrochem. ion gradients across the membrane. This is followed by an outline of its implementation and a description of key settings and parameters helpful to users wishing to set up and conduct such simulations. In recent years, key mechanistic and biophys. insights have been obtained by employing the CompEL protocol to address a wide range of questions on ion channels and permeation. We summarize these recent findings on membrane proteins, which span a spectrum from highly ion-selective, narrow channels to wide diffusion pores. Finally we discuss the future potential of CompEL in light of its limitations and strengths. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.
- 11Lin, Y.-L.; Aleksandrov, A.; Simonson, T.; Roux, B. An Overview of Electrostatic Free Energy Computations for Solutions and Proteins. J. Chem. Theory Comput. 2014, 10, 2690– 2709, DOI: 10.1021/ct500195p11https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXnslaksLc%253D&md5=6aa50d77dc3326ae518f682336b9f4a8An Overview of Electrostatic Free Energy Computations for Solutions and ProteinsLin, Yen-Lin; Aleksandrov, Alexey; Simonson, Thomas; Roux, BenoitJournal of Chemical Theory and Computation (2014), 10 (7), 2690-2709CODEN: JCTCCE; ISSN:1549-9618. (American Chemical Society)A review. Free energy simulations for electrostatic and charging processes in complex mol. systems encounter specific difficulties owing to the long-range, 1/r Coulomb interaction. To calc. the solvation free energy of a simple ion, it is essential to take into account the polarization of nearby solvent but also the electrostatic potential drop across the liq.-gas boundary, however distant. The latter does not exist in a simulation model based on periodic boundary conditions because there is no phys. boundary to the system. An important consequence is that the ref. value of the electrostatic potential is not an ion in a vacuum. Also, in an infinite system, the electrostatic potential felt by a perturbing charge is conditionally convergent and dependent on the choice of computational conventions. Furthermore, with Ewald lattice summation and tinfoil conducting boundary conditions, the charges experience a spurious shift in the potential that depends on the details of the simulation system such as the vol. fraction occupied by the solvent. All these issues can be handled with established computational protocols, as reviewed here and illustrated for several small ions and three solvated proteins.
- 12Darden, T.; York, D.; Pedersen, L. Particle mesh Ewald: An N·log(N) method for Ewald sums in large systems. J. Chem. Phys. 1993, 98, 10089– 10092, DOI: 10.1063/1.46439712https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK3sXks1Ohsr0%253D&md5=3c9f230bd01b7b714fd096d4d2e755f6Particle mesh Ewald: an N·log(N) method for Ewald sums in large systemsDarden, Tom; York, Darrin; Pedersen, LeeJournal of Chemical Physics (1993), 98 (12), 10089-92CODEN: JCPSA6; ISSN:0021-9606.An N·log(N) method for evaluating electrostatic energies and forces of large periodic systems is presented. The method is based on interpolation of the reciprocal space Ewald sums and evaluation of the resulting convolution using fast Fourier transforms. Timings and accuracies are presented for three large cryst. ionic systems.
- 13Essmann, U.; Perera, L.; Berkowitz, M. L.; Darden, T.; Lee, H.; Pedersen, L. G. A smooth particle mesh Ewald method. J. Chem. Phys. 1995, 103, 8577– 8593, DOI: 10.1063/1.47011713https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2MXptlehtrw%253D&md5=092a679dd3bee08da28df41e302383a7A smooth particle mesh Ewald methodEssmann, Ulrich; Perera, Lalith; Berkowitz, Max L.; Darden, Tom; Lee, Hsing; Pedersen, Lee G.Journal of Chemical Physics (1995), 103 (19), 8577-93CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The previously developed particle mesh Ewald method is reformulated in terms of efficient B-spline interpolation of the structure factors. This reformulation allows a natural extension of the method to potentials of the form 1/rp with p ≥ 1. Furthermore, efficient calcn. of the virial tensor follows. Use of B-splines in the place of Lagrange interpolation leads to analytic gradients as well as a significant improvement in the accuracy. The authors demonstrate that arbitrary accuracy can be achieved, independent of system size N, at a cost that scales as N log(N). For biomol. systems with many thousands of atoms and this method permits the use of Ewald summation at a computational cost comparable to that of a simple truncation method of 10 Å or less.
- 14Aksimentiev, A.; Schulten, K. Imaging alpha-hemolysin with molecular dynamics: ionic conductance, osmotic permeability, and the electrostatic potential map. Biophys. J. 2005, 88, 3745– 3761, DOI: 10.1529/biophysj.104.05872714https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2MXksl2jsL4%253D&md5=23548fa7b86cb4b542040ef5906d8d2bImaging α-hemolysin with molecular dynamics: Ionic conductance, osmotic permeability, and the electrostatic potential mapAksimentiev, Aleksij; Schulten, KlausBiophysical Journal (2005), 88 (6), 3745-3761CODEN: BIOJAU; ISSN:0006-3495. (Biophysical Society)α-Hemolysin of Staphylococcus aureus is a self-assembling toxin that forms a water-filled transmembrane channel upon oligomerization in a lipid membrane. Apart from being one of the best-studied toxins of bacterial origin, α-hemolysin is the principal component in several biotechnol. applications, including systems for controlled delivery of small solutes across lipid membranes, stochastic sensors for small solutes, and an alternative to conventional technol. for DNA sequencing. Through large-scale mol. dynamics simulations, we studied the permeability of the α-hemolysin/lipid bilayer complex for water and ions. The studied system, composed of ∼300,000 atoms, included one copy of the protein, a patch of a DPPC lipid bilayer, and a 1 M water soln. of KCl. Monitoring the fluctuations of the pore structure revealed an asym., on av., cross section of the α-hemolysin stem. Applying external electrostatic fields produced a transmembrane ionic current; repeating simulations at several voltage biases yielded a current/voltage curve of α-hemolysin and a set of electrostatic potential maps. The selectivity of α-hemolysin to Cl- was found to depend on the direction and the magnitude of the applied voltage bias. The results of our simulations are in excellent quant. agreement with available exptl. data. Analyzing trajectories of all water mol., we computed the α-hemolysin's osmotic permeability for water as well as its electroosmotic effect, and characterized the permeability of its seven side channels. The side channels were found to connect seven His-144 residues surrounding the stem of the protein to the bulk soln.; the protonation of these residues was obsd. to affect the ion conductance, suggesting the seven His-144 to comprise the pH sensor that gates conductance of the α-hemolysin channel.
- 15Humphrey, W.; Dalke, A.; Schulten, K. VMD: visual molecular dynamics. J. Mol. Graphics 1996, 14, 33– 38, DOI: 10.1016/0263-7855(96)00018-515https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28Xis12nsrg%253D&md5=1e3094ec3151fb85c5ff05f8505c78d5VDM: visual molecular dynamicsHumphrey, William; Dalke, Andrew; Schulten, KlausJournal of Molecular Graphics (1996), 14 (1), 33-8, plates, 27-28CODEN: JMGRDV; ISSN:0263-7855. (Elsevier)VMD is a mol. graphics program designed for the display and anal. of mol. assemblies, in particular, biopolymers such as proteins and nucleic acids. VMD can simultaneously display any no. of structures using a wide variety of rendering styles and coloring methods. Mols. are displayed as one or more "representations," in which each representation embodies a particular rendering method and coloring scheme for a selected subset of atoms. The atoms displayed in each representation are chosen using an extensive atom selection syntax, which includes Boolean operators and regular expressions. VMD provides a complete graphical user interface for program control, as well as a text interface using the Tcl embeddable parser to allow for complex scripts with variable substitution, control loops, and function calls. Full session logging is supported, which produces a VMD command script for later playback. High-resoln. raster images of displayed mols. may be produced by generating input scripts for use by a no. of photorealistic image-rendering applications. VMD has also been expressly designed with the ability to animate mol. dynamics (MD) simulation trajectories, imported either from files or from a direct connection to a running MD simulation. VMD is the visualization component of MDScope, a set of tools for interactive problem solving in structural biol., which also includes the parallel MD program NAMD, and the MDCOMM software used to connect the visualization and simulation programs, VMD is written in C++, using an object-oriented design; the program, including source code and extensive documentation, is freely available via anonymous ftp and through the World Wide Web.
- 16Enkavi, G.; Tajkhorshid, E. Simulation of spontaneous substrate binding revealing the binding pathway and mechanism and initial conformational response of GlpT. Biochemistry 2010, 49, 1105– 1114, DOI: 10.1021/bi901412a16https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXovVaktg%253D%253D&md5=be122d6c0290699b0453d120fc26f267Simulation of Spontaneous Substrate Binding Revealing the Binding Pathway and Mechanism and Initial Conformational Response of GlpTEnkavi, Giray; Tajkhorshid, EmadBiochemistry (2010), 49 (6), 1105-1114CODEN: BICHAW; ISSN:0006-2960. (American Chemical Society)Glycerol 3-phosphate transporter (GlpT) mediates the import of glycerol 3-phosphate (G3P) using the gradient of inorg. phosphate (Pi). To study the process and mechanism of substrate binding and to investigate the protein's initial response, we performed equil. simulations of wild-type GlpT and several of its mutant forms in membranes in the presence of all physiol. relevant substrates (Pi-, Pi2-, G3P-, and G3P2-). The simulations capture spontaneous substrate binding of GlpT, driven by the pos. electrostatic potential of the lumen. K80 is found to act as a "hook" making the first encounter with the substrate and guiding it toward the binding site, where it binds tightly to R45, a key binding site residue that acts as a "fork" holding the substrate. R269 establishes no direct contact with the substrate during the simulations, a surprising behavior given its structural pseudosymmetry to R45. In all substrate-bound systems, partial closing of the cytoplasmic half of GlpT was obsd. The substrate appears to stabilize the partially occluded state, as in the two apo simulations either no closing was obsd. or the protein reverted to its open form toward the end of the simulation, whereas in all substrate-bound systems, a stable partially closed state was produced. Along with the modulation of the periplasmic salt bridge network, these substrate-induced events destabilize the periplasmic half while inducing a closure in the cytoplasmic half, thus capturing the early stages of the proposed rocker-switch mechanism in GlpT.
- 17Arcario, M. J.; Tajkhorshid, E. Membrane-induced structural rearrangement and identification of a novel membrane anchor in talin F2F3. Biophys. J. 2014, 107, 2059– 2069, DOI: 10.1016/j.bpj.2014.09.02217https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhslGrsLzO&md5=c6de1b54606c2757ec9a32d21ca1e489Membrane-Induced Structural Rearrangement and Identification of a Novel Membrane Anchor in Talin F2F3Arcario, Mark J.; Tajkhorshid, EmadBiophysical Journal (2014), 107 (9), 2059-2069CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)Exptl. challenges assocd. with characterization of the membrane-bound form of talin have prevented us from understanding the mol. mechanism of its membrane-dependent integrin activation. Here, utilizing what we believe to be a novel membrane mimetic model, we present a reproducible model of membrane-bound talin obsd. across multiple independent simulations. We characterize both local and global membrane-induced structural transitions that successfully reconcile discrepancies between biochem. and structural studies and provide insight into how talin might modulate integrin function. Membrane binding of talin, captured in unbiased simulations, proceeds through three distinct steps: initial electrostatic recruitment of the F2 subdomain to anionic lipids via several basic residues; insertion of an initially buried, conserved hydrophobic anchor into the membrane; and assocn. of the F3 subdomain with the membrane surface through a large, interdomain conformational change. These latter two steps, to our knowledge, have not been obsd. or described previously. Electrostatic anal. shows talin F2F3 to be highly polarized, with a highly pos. underside, which we attribute to the initial electrostatic recruitment, and a neg. top face, which can help orient the protein optimally with respect to the membrane, thereby reducing the no. of unproductive membrane collision events.
- 18Frigo, M.; Johnson, S. G. The Design and Implementation of FFTW3. Proc. IEEE 2005, 93, 216– 231, DOI: 10.1109/JPROC.2004.840301There is no corresponding record for this reference.
- 19Phillips, J. C.; Hardy, D. J.; Maia, J. D. C.; Stone, J. E.; Ribeiro, J. V.; Bernardi, R. C.; Buch, R.; Fiorin, G.; Hénin, J.; Jiang, W.; McGreevy, R.; Melo, M. C. R.; Radak, B. K.; Skeel, R. D.; Singharoy, A.; Wang, Y.; Roux, B.; Aksimentiev, A.; Luthey-Schulten, Z.; Kalé, L. V.; Schulten, K.; Chipot, C.; Tajkhorshid, E. Scalable molecular dynamics on CPU and GPU architectures with NAMD. J. Chem. Phys. 2020, 153, 044130 DOI: 10.1063/5.001447519https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhsFajsL3J&md5=aa6378c15e48addc4b6b02523e55427fScalable molecular dynamics on CPU and GPU architectures with NAMDPhillips, James C.; Hardy, David J.; Maia, Julio D. C.; Stone, John E.; Ribeiro, Joao V.; Bernardi, Rafael C.; Buch, Ronak; Fiorin, Giacomo; Henin, Jerome; Jiang, Wei; McGreevy, Ryan; Melo, Marcelo C. R.; Radak, Brian K.; Skeel, Robert D.; Singharoy, Abhishek; Wang, Yi; Roux, Benoit; Aksimentiev, Aleksei; Luthey-Schulten, Zaida; Kale, Laxmikant V.; Schulten, Klaus; Chipot, Christophe; Tajkhorshid, EmadJournal of Chemical Physics (2020), 153 (4), 044130CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)A review. NAMD is a mol. dynamics program designed for high-performance simulations of very large biol. objects on CPU- and GPU-based architectures. NAMD offers scalable performance on petascale parallel supercomputers consisting of hundreds of thousands of cores, as well as on inexpensive commodity clusters commonly found in academic environments. It is written in C++ and leans on Charm++ parallel objects for optimal performance on low-latency architectures. NAMD is a versatile, multipurpose code that gathers state-of-the-art algorithms to carry out simulations in apt thermodn. ensembles, using the widely popular CHARMM, AMBER, OPLS, and GROMOS biomol. force fields. Here, the authors review the main features of NAMD that allow both equil. and enhanced-sampling mol. dynamics simulations with numerical efficiency. The authors describe the underlying concepts used by NAMD and their implementation, most notably for handling long-range electrostatics; controlling the temp., pressure, and pH; applying external potentials on tailored grids; leveraging massively parallel resources in multiple-copy simulations; and hybrid quantum-mech./mol.-mech. descriptions. The authors detail the variety of options offered by NAMD for enhanced-sampling simulations aimed at detg. free-energy differences of either alchem. or geometrical transformations and outline their applicability to specific problems. Last, the roadmap for the development of NAMD and the authors' current efforts toward achieving optimal performance on GPU-based architectures, for pushing back the limitations that have prevented biol. realistic billion-atom objects to be fruitfully simulated, and for making large-scale simulations less expensive and easier to set up, run, and analyze are discussed. NAMD is distributed free of charge with its source code at www.ks.uiuc.edu. (c) 2020 American Institute of Physics.
- 20Mansoor, S. E.; Lü, W.; Oosterheert, W.; Shekhar, M.; Tajkhorshid, E.; Gouaux, E. X-ray structures define human P2X(3) receptor gating cycle and antagonist action. Nature 2016, 538, 66– 71, DOI: 10.1038/nature1936720https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhsFWiurnJ&md5=db6de173d819f23ba4304d0290ba0000X-ray structures define human P2X3 receptor gating cycle and antagonist actionMansoor, Steven E.; Lu, Wei; Oosterheert, Wout; Shekhar, Mrinal; Tajkhorshid, Emad; Gouaux, EricNature (London, United Kingdom) (2016), 538 (7623), 66-71CODEN: NATUAS; ISSN:0028-0836. (Nature Publishing Group)P2X receptors are trimeric, non-selective cation channels activated by ATP that have important roles in the cardiovascular, neuronal and immune systems. Despite their central function in human physiol. and although they are potential targets of therapeutic agents, there are no structures of human P2X receptors. The mechanisms of receptor desensitization and ion permeation, principles of antagonism, and complete structures of the pore-forming transmembrane domains of these receptors remain unclear. Here we report X-ray crystal structures of the human P2X3 receptor in apo/resting, agonist-bound/open-pore, agonist-bound/closed-pore/desensitized and antagonist-bound/closed states. The open state structure harbours an intracellular motif we term the 'cytoplasmic cap', which stabilizes the open state of the ion channel pore and creates lateral, phospholipid-lined cytoplasmic fenestrations for water and ion egress. The competitive antagonists TNP-ATP and A-317491 stabilize the apo/resting state and reveal the interactions responsible for competitive inhibition. These structures illuminate the conformational rearrangements that underlie P2X receptor gating and provide a foundation for the development of new pharmacol. agents.
- 21Wolf, M. G.; Hoefling, M.; Aponte-Santamaría, C.; Grubmüller, H.; Groenhof, G. g_membed: Efficient insertion of a membrane protein into an equilibrated lipid bilayer with minimal perturbation. J. Comput. Chem. 2010, 31, 2169– 2174, DOI: 10.1002/jcc.2150721https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXnsFWjs7c%253D&md5=1f0353d68fb867fb684283b61d52eadbg_membed: Efficient insertion of a membrane protein into an equilibrated lipid bilayer with minimal perturbationWolf, Maarten G.; Hoefling, Martin; Aponte-Santamaria, Camilo; Grubmueller, Helmut; Groenhof, GerritJournal of Computational Chemistry (2010), 31 (11), 2169-2174CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)To efficiently insert a protein into an equilibrated and fully hydrated membrane with minimal membrane perturbation the authors present a computational tool, called g_membed, which is part of the Gromacs suite of programs. The input consists of an equilibrated membrane system, either flat or curved, and a protein structure in the right position and orientation with respect to the lipid bilayer. G_membed first decreases the width of the protein in the xy-plane and removes all mols. (generally lipids and waters) that overlap with the narrowed protein. Then the protein is grown back to its full size in a short mol. dynamics simulation (typically 1000 steps), thereby pushing the lipids away to optimally accommodate the protein in the membrane. After embedding the protein in the membrane, both the lipid properties and the hydration layer are still close to equil. Thus, only a short equilibration run (less then 1 ns in the cases tested) is required to re-equilibrate the membrane. Its simplicity makes g_membed very practical for use in scripting and high-throughput mol. dynamics simulations. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010.
- 22Abraham, M. J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J. C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1–2, 19– 25, DOI: 10.1016/j.softx.2015.06.001There is no corresponding record for this reference.
- 23Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. R. Molecular dynamics with coupling to an external bath. J. Chem. Phys. 1984, 81, 3684– 3690, DOI: 10.1063/1.44811823https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL2cXmtlGksbY%253D&md5=5510dc00297d63b91ee3a7a4ae5aacb1Molecular dynamics with coupling to an external bathBerendsen, H. J. C.; Postma, J. P. M.; Van Gunsteren, W. F.; DiNola, A.; Haak, J. R.Journal of Chemical Physics (1984), 81 (8), 3684-90CODEN: JCPSA6; ISSN:0021-9606.In mol. dynamics (MD) simulations, the need often arises to maintain such parameters as temp. or pressure rather than energy and vol., or to impose gradients for studying transport properties in nonequil. MD. A method is described to realize coupling to an external bath with const. temp. or pressure with adjustable time consts. for the coupling. The method is easily extendable to other variables and to gradients, and can be applied also to polyat. mols. involving internal constraints. The influence of coupling time consts. on dynamical variables is evaluated. A leap-frog algorithm is presented for the general case involving constraints with coupling to both a const. temp. and a const. pressure bath.
- 24Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. J. Chem. Phys. 2007, 126, 014101 DOI: 10.1063/1.240842024https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXosVCltg%253D%253D&md5=9c182b57bfc65bca6be23c8c76b4be77Canonical sampling through velocity rescalingBussi, Giovanni; Donadio, Davide; Parrinello, MicheleJournal of Chemical Physics (2007), 126 (1), 014101/1-014101/7CODEN: JCPSA6; ISSN:0021-9606. (American Institute of Physics)The authors present a new mol. dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains const. during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. The authors illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liq. phases. Its performance is excellent and largely independent of the thermostat parameter also with regard to the dynamic properties.
- 25Huang, J.; Rauscher, S.; Nawrocki, G.; Ran, T.; Feig, M.; de Groot, B. L.; Grubmüller, H.; MacKerell, A. D. CHARMM36m: an improved force field for folded and intrinsically disordered proteins. Nat. Methods 2017, 14, 71– 73, DOI: 10.1038/nmeth.406725https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhvVSiu77I&md5=0aa151fbef2ee0b5e2cfb593c54330c2CHARMM36m: an improved force field for folded and intrinsically disordered proteinsHuang, Jing; Rauscher, Sarah; Nawrocki, Grzegorz; Ran, Ting; Feig, Michael; de Groot, Bert L.; Grubmuller, Helmut; MacKerell, Alexander D. JrNature Methods (2017), 14 (1), 71-73CODEN: NMAEA3; ISSN:1548-7091. (Nature Publishing Group)The all-atom additive CHARMM36 protein force field is widely used in mol. modeling and simulations. We present its refinement, CHARMM36m (http://mackerell.umaryland.edu/charmm_ff.shtml), with improved accuracy in generating polypeptide backbone conformational ensembles for intrinsically disordered peptides and proteins.
- 26Klauda, J. B.; Venable, R. M.; Freites, J. A.; O’Connor, J. W.; Tobias, D. J.; Mondragon-Ramirez, C.; Vorobyov, I.; MacKerell, A. D.; Pastor, R. W. Update of the CHARMM all-atom additive force field for lipids: validation on six lipid types. J. Phys. Chem. B 2010, 114, 7830– 7843, DOI: 10.1021/jp101759q26https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXmsVGjtrk%253D&md5=04c4d3aa1a59740190994d7df80cbb71Update of the CHARMM All-Atom Additive Force Field for Lipids: Validation on Six Lipid TypesKlauda, Jeffery B.; Venable, Richard M.; Freites, J. Alfredo; O'Connor, Joseph W.; Tobias, Douglas J.; Mondragon-Ramirez, Carlos; Vorobyov, Igor; MacKerell, Alexander D., Jr.; Pastor, Richard W.Journal of Physical Chemistry B (2010), 114 (23), 7830-7843CODEN: JPCBFK; ISSN:1520-6106. (American Chemical Society)A significant modification to the additive all-atom CHARMM lipid force field (FF) is developed and applied to phospholipid bilayers with both choline and ethanolamine contg. head groups and with both satd. and unsatd. aliph. chains. Motivated by the current CHARMM lipid FF (C27 and C27r) systematically yielding values of the surface area per lipid that are smaller than exptl. ests. and gel-like structures of bilayers well above the gel transition temp., selected torsional, Lennard-Jones and partial at. charge parameters were modified by targeting both quantum mech. (QM) and exptl. data. QM calcns. ranging from high-level ab initio calcns. on small mols. to semiempirical QM studies on a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer in combination with exptl. thermodn. data were used as target data for parameter optimization. These changes were tested with simulations of pure bilayers at high hydration of the following six lipids: DPPC, 1,2-dimyristoyl-sn-phosphatidylcholine (DMPC), 1,2-dilauroyl-sn-phosphatidylcholine (DLPC), 1-palmitoyl-2-oleoyl-sn-phosphatidylcholine (POPC), 1,2-dioleoyl-sn-phosphatidylcholine (DOPC), and 1-palmitoyl-2-oleoyl-sn-phosphatidylethanolamine (POPE); simulations of a low hydration DOPC bilayer were also performed. Agreement with exptl. surface area is on av. within 2%, and the d. profiles agree well with neutron and x-ray diffraction expts. NMR deuterium order parameters (SCD) are well predicted with the new FF, including proper splitting of the SCD for the aliph. carbon adjacent to the carbonyl for DPPC, POPE, and POPC bilayers. The area compressibility modulus and frequency dependence of 13C NMR relaxation rates of DPPC and the water distribution of low hydration DOPC bilayers also agree well with expt. Accordingly, the presented lipid FF, referred to as C36, allows for mol. dynamics simulations to be run in the tensionless ensemble (NPT), and is anticipated to be of utility for simulations of pure lipid systems as well as heterogeneous systems including membrane proteins.
- 27Jorgensen, W. L.; Chandrasekhar, J.; Madura, J. D.; Impey, R. W.; Klein, M. L. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926– 935, DOI: 10.1063/1.44586927https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL3sXksF2htL4%253D&md5=a1161334e381746be8c9b15a5e56f704Comparison of simple potential functions for simulating liquid waterJorgensen, William L.; Chandrasekhar, Jayaraman; Madura, Jeffry D.; Impey, Roger W.; Klein, Michael L.Journal of Chemical Physics (1983), 79 (2), 926-35CODEN: JCPSA6; ISSN:0021-9606.Classical Monte Carlo simulations were carried out for liq. H2O in the NPT ensemble at 25° and 1 atm using 6 of the simpler intermol. potential functions for the dimer. Comparisons were made with exptl. thermodn. and structural data including the neutron diffraction results of Thiessen and Narten (1982). The computed densities and potential energies agree with expt. except for the original Bernal-Fowler model, which yields an 18% overest. of the d. and poor structural results. The discrepancy may be due to the correction terms needed in processing the neutron data or to an effect uniformly neglected in the computations. Comparisons were made for the self-diffusion coeffs. obtained from mol. dynamics simulations.
- 28Brunner, J. D.; Lim, N. K.; Schenck, S.; Duerst, A.; Dutzler, R. X-ray structure of a calcium-activated TMEM16 lipid scramblase. Nature 2014, 516, 207– 212, DOI: 10.1038/nature1398428https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhvFGltLrO&md5=828d54af1de0c0bd05a31a1eed8a6fdfX-ray structure of a calcium-activated TMEM16 lipid scramblaseBrunner, Janine D.; Lim, Novandy K.; Schenck, Stephan; Duerst, Alessia; Dutzler, RaimundNature (London, United Kingdom) (2014), 516 (7530), 207-212CODEN: NATUAS; ISSN:0028-0836. (Nature Publishing Group)The TMEM16 family of proteins, also known as anoctamins, features a remarkable functional diversity. This family contains the long sought-after Ca2+-activated chloride channels as well as lipid scramblases and cation channels. Here we present the crystal structure of a TMEM16 family member from the fungus Nectria hematococca that operates as a Ca2+-activated lipid scramblase. Each subunit of the homodimeric protein contains ten transmembrane helixes and a hydrophilic membrane-traversing cavity that is exposed to the lipid bilayer as a potential site of catalysis. This cavity harbors a conserved Ca2+-binding site located within the hydrophobic core of the membrane. Mutations of residues involved in Ca2+ coordination affect both lipid scrambling in N. hematococca TMEM16 and ion conduction in the Cl- channel TMEM16A. The structure reveals the general architecture of the family and its mode of Ca2+ activation. It also provides insight into potential scrambling mechanisms and serves as a framework to unravel the conduction of ions in certain TMEM16 proteins.
- 29Bushell, S. R.; Pike, A. C.; Falzone, M. E.; Rorsman, N. J.; Ta, C. M.; Corey, R. A.; Newport, T. D.; Christianson, J. C.; Scofano, L. F.; Shintre, C. A.; Tessitore, A.; Chu, A.; Wang, Q.; Shrestha, L.; Mukhopadhyay, S. M. M.; Love, J. D.; Burgess-Brown, N. A.; Sitsapesan, R.; Stansfeld, P. J.; Huiskonen, J. T.; Tammaro, P.; Accardi, A.; Carpenter, E. P. The structural basis of lipid scrambling and inactivation in the endoplasmic reticulum scramblase TMEM16K. Nat. Commun. 2019, 10, 3956 DOI: 10.1038/s41467-019-11753-129https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3MrksVCrtw%253D%253D&md5=f146bbc9a998091e2825f45ff419a859The structural basis of lipid scrambling and inactivation in the endoplasmic reticulum scramblase TMEM16KBushell Simon R; Pike Ashley C W; Shintre Chitra A; Tessitore Annamaria; Chu Amy; Wang Qinrui; Shrestha Leela; Mukhopadhyay Shubhashish M M; Burgess-Brown Nicola A; Carpenter Elisabeth P; Falzone Maria E; Accardi Alessio; Rorsman Nils J G; Ta Chau M; Scofano Lara F; Sitsapesan Rebecca; Tammaro Paolo; Rorsman Nils J G; Ta Chau M; Corey Robin A; Newport Thomas D; Wang Qinrui; Stansfeld Phillip J; Newport Thomas D; Christianson John C; Shintre Chitra A; Tessitore Annamaria; Chu Amy; Love James D; Love James D; Huiskonen Juha T; Accardi Alessio; Accardi AlessioNature communications (2019), 10 (1), 3956 ISSN:.Membranes in cells have defined distributions of lipids in each leaflet, controlled by lipid scramblases and flip/floppases. However, for some intracellular membranes such as the endoplasmic reticulum (ER) the scramblases have not been identified. Members of the TMEM16 family have either lipid scramblase or chloride channel activity. Although TMEM16K is widely distributed and associated with the neurological disorder autosomal recessive spinocerebellar ataxia type 10 (SCAR10), its location in cells, function and structure are largely uncharacterised. Here we show that TMEM16K is an ER-resident lipid scramblase with a requirement for short chain lipids and calcium for robust activity. Crystal structures of TMEM16K show a scramblase fold, with an open lipid transporting groove. Additional cryo-EM structures reveal extensive conformational changes from the cytoplasmic to the ER side of the membrane, giving a state with a closed lipid permeation pathway. Molecular dynamics simulations showed that the open-groove conformation is necessary for scramblase activity.
- 30Kostritskii, A. Y.; Machtens, J.-P. Molecular mechanisms of ion conduction and ion selectivity in TMEM16 lipid scramblases. Nat. Commun. 2021, in press, DOI: 10.1038/s41467-021-22724-w .There is no corresponding record for this reference.
- 31Alleva, C.; Kovalev, K.; Astashkin, R.; Berndt, M. I.; Baeken, C.; Balandin, T.; Gordeliy, V.; Fahlke, C.; Machtens, J.-P. Na+-dependent gate dynamics and electrostatic attraction ensure substrate coupling in glutamate transporters. Sci. Adv. 2020, 6, eaba9854 DOI: 10.1126/sciadv.aba9854There is no corresponding record for this reference.
- 32Lomize, M. A.; Pogozheva, I. D.; Joo, H.; Mosberg, H. I.; Lomize, A. L. OPM database and PPM web server: resources for positioning of proteins in membranes. Nucleic Acids Res. 2012, 40, D370– D376, DOI: 10.1093/nar/gkr70332https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXhs12hurzJ&md5=b931126c56d227ab0eb6c5762c1dae6bOPM database and PPM web server: resources for positioning of proteins in membranesLomize, Mikhail A.; Pogozheva, Irina D.; Joo, Hyeon; Mosberg, Henry I.; Lomize, Andrei L.Nucleic Acids Research (2012), 40 (D1), D370-D376CODEN: NARHAD; ISSN:0305-1048. (Oxford University Press)The Orientations of Proteins in Membranes (OPM) database is a curated web resource that provides spatial positions of membrane-bound peptides and proteins of known three-dimensional structure in the lipid bilayer, together with their structural classification, topol. and intracellular localization. OPM currently contains more than 1200 transmembrane and peripheral proteins and peptides from approx. 350 organisms that represent approx. 3800 Protein Data Bank entries. Proteins are classified into classes, superfamilies and families and assigned to 21 distinct membrane types. Spatial positions of proteins with respect to the lipid bilayer are optimized by the PPM 2.0 method that accounts for the hydrophobic, hydrogen bonding and electrostatic interactions of the proteins with the anisotropic water-lipid environment described by the dielec. const. and hydrogen-bonding profiles. The OPM database is freely accessible at http://opm.phar.umich.edu. Data can be sorted, searched or retrieved using the hierarchical classification, source organism, localization in different types of membranes. The database offers downloadable coordinates of proteins and peptides with membrane boundaries. A gallery of protein images and several visualization tools are provided. The database is supplemented by the PPM server (http://opm.phar.umich.edu/server.php) which can be used for calcg. spatial positions in membranes of newly detd. proteins structures or theor. models.
- 33Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J. L.; Dror, R. O.; Shaw, D. E. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins: Struct., Funct., Bioinf. 2010, 78, 1950– 1958, DOI: 10.1002/prot.2271133https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXkvFegtLo%253D&md5=447a9004026e2b93f0f7beff165daa09Improved side-chain torsion potentials for the Amber ff99SB protein force fieldLindorff-Larsen, Kresten; Piana, Stefano; Palmo, Kim; Maragakis, Paul; Klepeis, John L.; Dror, Ron O.; Shaw, David E.Proteins: Structure, Function, and Bioinformatics (2010), 78 (8), 1950-1958CODEN: PSFBAF ISSN:. (Wiley-Liss, Inc.)Recent advances in hardware and software have enabled increasingly long mol. dynamics (MD) simulations of biomols., exposing certain limitations in the accuracy of the force fields used for such simulations and spurring efforts to refine these force fields. Recent modifications to the Amber and CHARMM protein force fields, for example, have improved the backbone torsion potentials, remedying deficiencies in earlier versions. Here, the authors further advance simulation accuracy by improving the amino acid side-chain torsion potentials of the Amber ff99SB force field. First, the authors used simulations of model alpha-helical systems to identify the four residue types whose rotamer distribution differed the most from expectations based on Protein Data Bank statistics. Second, the authors optimized the side-chain torsion potentials of these residues to match new, high-level quantum-mech. calcns. Finally, the authors used microsecond-timescale MD simulations in explicit solvent to validate the resulting force field against a large set of exptl. NMR measurements that directly probe side-chain conformations. The new force field, which the authors have termed Amber ff99SB-ILDN, exhibits considerably better agreement with the NMR data. Proteins 2010. © 2010 Wiley-Liss, Inc.
- 34Berger, O.; Edholm, O.; Jähnig, F. Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperature. Biophys. J. 1997, 72, 2002– 2013, DOI: 10.1016/S0006-3495(97)78845-334https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK2sXivVOjsL4%253D&md5=984ea49ced9d0e2ce912bf99220fd228Molecular dynamics simulations of a fluid bilayer of dipalmitoylphosphatidylcholine at full hydration, constant pressure, and constant temperatureBerger, Oliver; Edholm, Olle; Jahnig, FritzBiophysical Journal (1997), 72 (5), 2002-2013CODEN: BIOJAU; ISSN:0006-3495. (Biophysical Society)Mol. dynamics simulations of 500 ps were performed on a system consisting of a bilayer of 64 mols. of the lipid dipalmitoylphosphatidylcholine and 23 water mols. per lipid at an isotropic pressure of 1 atm and 50°. Special attention was devoted to reproduce the correct d. of the lipid, because this quantity is known exptl. with a precision better than 1%. For this purpose, the Lennard-Jones parameters of the hydrocarbon chains were adjusted by simulating a system consisting of 128 pentadecane mols. and varying the Lennard-Jones parameters until the exptl. d. and heat of vaporization were obtained. With these parameters the lipid d. resulted in perfect agreement with the exptl. d. The orientational order parameter of the hydrocarbon chains agreed perfectly well with the exptl. values, which, because of its correlation with the area per lipid, makes it possible to give a proper est. of the area per lipid of 0.61±0.1 nm2.
- 35Joung, I. S.; Cheatham, T. E. Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular Simulations. J. Phys. Chem. B 2008, 112, 9020– 9041, DOI: 10.1021/jp800161435https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXnvFGqtL4%253D&md5=aa489470ae1c7479bf0911710217bd28Determination of Alkali and Halide Monovalent Ion Parameters for Use in Explicitly Solvated Biomolecular SimulationsJoung, In Suk; Cheatham, Thomas E.Journal of Physical Chemistry B (2008), 112 (30), 9020-9041CODEN: JPCBFK; ISSN:1520-6106. (American Chemical Society)Alkali (Li+, Na+, K+, Rb+, and Cs+) and halide (F-, Cl-, Br-, and I-) ions play an important role in many biol. phenomena, roles that range from stabilization of biomol. structure, to influence on biomol. dynamics, to key physiol. influence on homeostasis and signaling. To properly model ionic interaction and stability in atomistic simulations of biomol. structure, dynamics, folding, catalysis, and function, an accurate model or representation of the monovalent ions is critically necessary. A good model needs to simultaneously reproduce many properties of ions, including their structure, dynamics, solvation, and moreover both the interactions of these ions with each other in the crystal and in soln. and the interactions of ions with other mols. At present, the best force fields for biomols. employ a simple additive, nonpolarizable, and pairwise potential for at. interaction. In this work, the authors describe their efforts to build better models of the monovalent ions within the pairwise Coulombic and 6-12 Lennard-Jones framework, where the models are tuned to balance crystal and soln. properties in Ewald simulations with specific choices of well-known water models. Although it has been clearly demonstrated that truly accurate treatments of ions will require inclusion of nonadditivity and polarizability (particularly with the anions) and ultimately even a quantum mech. treatment, the authors' goal was to simply push the limits of the additive treatments to see if a balanced model could be created. The applied methodol. is general and can be extended to other ions and to polarizable force-field models. The authors' starting point centered on observations from long simulations of biomols. in salt soln. with the AMBER force fields where salt crystals formed well below their soly. limit. The likely cause of the artifact in the AMBER parameters relates to the naive mixing of the Smith and Dang chloride parameters with AMBER-adapted Aqvist cation parameters. To provide a more appropriate balance, the authors reoptimized the parameters of the Lennard-Jones potential for the ions and specific choices of water models. To validate and optimize the parameters, the authors calcd. hydration free energies of the solvated ions and also lattice energies (LE) and lattice consts. (LC) of alkali halide salt crystals. This is the first effort that systematically scans across the Lennard-Jones space (well depth and radius) while balancing ion properties like LE and LC across all pair combinations of the alkali ions and halide ions. The optimization across the entire monovalent series avoids systematic deviations. The ion parameters developed, optimized, and characterized were targeted for use with some of the most commonly used rigid and nonpolarizable water models, specifically TIP3P, TIP4PEW, and SPC/E. In addn. to well reproducing the soln. and crystal properties, the new ion parameters well reproduce binding energies of the ions to water and the radii of the first hydration shells.
- 36Drew, H. R.; Dickerson, R. E. Structure of a B-DNA dodecamer. III. Geometry of hydration. J. Mol. Biol. 1981, 151, 535– 556, DOI: 10.1016/0022-2836(81)90009-736https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL38XmsVentg%253D%253D&md5=4378664a4da9ec83155a3f7e046027d9Structure of a B-DNA dodecamer. III. Geometry of hydrationDrew, Horace R.; Dickerson, Richard E.Journal of Molecular Biology (1981), 151 (3), 535-56CODEN: JMOBAK; ISSN:0022-2836.The B-DNA dodecamer, C-G-C-G-A-A-T-T-C-G-C-G, crystd. as slightly more than 1 full turn of right-handed B-DNA in space group P212121 with cell dimensions a = 24.87 Å, b = 40.39 Å, and c = 66.20 Å. X-ray anal. showed that it was surrounded by 72 ordered water mols., mostly assocd. with polar N and O atoms at exposed edges of base-pairs. Hydration in the major groove was mainly in the form of a monodentate monolayer. Hydration of backbone phosphate O atoms were not ordered, except when immobilized by the 5-Me groups of adjacent thymines. The minor groove was extensively and regularly hydrated, with a zigzag spine of 1st- and 2nd-shell hydration along the floor of the groove serving as a foundation for less-regular outer shells extending beyond the radius of the phosphate backbone. The spine network bridged purine N-3 and pyrimidine O-2 atoms in adjacent base pairs; it was regular in the A-A-T-T center, but was disrupted at the C-G-C-G ends, partly by the guanine N-2 amino groups. The minor groove hydration spine may be responsible for the stability of the B form of polymers contg. only A-T and I-C base pairs, and its disruption may explain the ease of transition to the A form of polymers with G-C pairs.
- 37Pérez, A.; Marchán, I.; Svozil, D.; Sponer, J.; Cheatham, T. E.; Laughton, C. A.; Orozco, M. Refinement of the AMBER force field for nucleic acids: improving the description of alpha/gamma conformers. Biophys. J. 2007, 92, 3817– 3829, DOI: 10.1529/biophysj.106.09778237https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXls1ygu7Y%253D&md5=e0f8d87e5fe7710e59cc41f5e8c25f7bRefinement of the AMBER force field for nucleic acids: improving the description of α/γ conformersPerez, Alberto; Marchan, Ivan; Svozil, Daniel; Sponer, Jiri; Cheatham, Thomas E., III; Laughton, Charles A.; Orozco, ModestoBiophysical Journal (2007), 92 (11), 3817-3829CODEN: BIOJAU; ISSN:0006-3495. (Biophysical Society)We present here the parmbsc0 force field, a refinement of the AMBER parm99 force field, where emphasis has been made on the correct representation of the α/γ concerted rotation in nucleic acids (NAs). The modified force field corrects overpopulations of the α/γ = (g+,t) backbone that were seen in long (more than 10 ns) simulations with previous AMBER parameter sets (parm94-99). The force field has been derived by fitting to high-level quantum mech. data and verified by comparison with very high-level quantum mech. calcns. and by a very extensive comparison between simulations and exptl. data. The set of validation simulations includes two of the longest trajectories published to date for the DNA duplex (200 ns each) and the largest variety of NA structures studied to date (15 different NA families and 97 individual structures). The total simulation time used to validate the force field includes near 1 μs of state-of-the-art mol. dynamics simulations in aq. soln.
- 38Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. J. Appl. Phys. 1981, 52, 7182– 7190, DOI: 10.1063/1.32869338https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL38XislSnuw%253D%253D&md5=a0a5617389f6cabbf2a405c649aadf03Polymorphic transitions in single crystals: A new molecular dynamics methodParrinello, M.; Rahman, A.Journal of Applied Physics (1981), 52 (12), 7182-90CODEN: JAPIAU; ISSN:0021-8979.A Lagrangian formulation is introduced; it can be used to make mol. dynamics (MD) calcns. on systems under the most general, externally applied, conditions of stress. In this formulation the MD cell shape and size can change according to dynamic equations given by this Lagrangian. This MD technique was used to the study of structural transitions of a Ni single crystal under uniform uniaxial compressive and tensile loads. Some results regarding the stress-strain relation obtained by static calcns. are invalid at finite temp. Under compressive loading, the model of Ni shows a bifurcation in its stress-strain relation; this bifurcation provides a link in configuration space between cubic and hexagonal close packing. Such a transition could perhaps be obsd. exptl. under extreme conditions of shock.
- 39Kondinskaia, D. A.; Kostritskii, A. Y.; Nesterenko, A. M.; Antipina, A. Y.; Gurtovenko, A. A. Atomic-Scale Molecular Dynamics Simulations of DNA-Polycation Complexes: Two Distinct Binding Patterns. J. Phys. Chem. B 2016, 120, 6546– 6554, DOI: 10.1021/acs.jpcb.6b0377939https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XpsFGlsLg%253D&md5=ab05858cf4c832df5a5027e197842427Atomic-Scale Molecular Dynamics Simulations of DNA-Polycation Complexes: Two Distinct Binding PatternsKondinskaia, Diana A.; Kostritskii, Andrei Yu.; Nesterenko, Alexey M.; Antipina, Alexandra Yu.; Gurtovenko, Andrey A.Journal of Physical Chemistry B (2016), 120 (27), 6546-6554CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Synthetic cationic polymers represent a promising class of delivery vectors for gene therapy. Here, we employ atomistic mol. dynamics simulations to gain insight into the structure and properties of complexes of DNA with four linear polycations: polyethylenimine (PEI), poly-L-lysine (PLL), polyvinylamine (PVA), and polyallylamine (PAA). These polycations differ in their polymer geometries, protonation states, and hydrophobicities of their backbone chains. Overall, our results demonstrate for the first time the existence of two distinct patterns of binding of DNA with polycations. For PEI, PLL, and PAA, the complex is stabilized by the electrostatic attraction between protonated amine groups of the polycation and phosphate groups of DNA. In contrast, PVA demonstrates an alternative binding pattern as it gets embedded into the DNA major groove. It is likely that both the polymer topol. and affinity of the backbone chain of PVA to the DNA groove are responsible for such behavior. The differences in binding patterns can have important biomedical implications: embedding PVA into a DNA groove makes it less sensitive to changes in the aq. environment (pH level, ionic strength, etc.) and could therefore hinder the intracellular release of genetic material from a delivery vector, leading to lower transfection activity.
- 40Michaud-Agrawal, N.; Denning, E. J.; Woolf, T. B.; Beckstein, O. MDAnalysis: a toolkit for the analysis of molecular dynamics simulations. J. Comput. Chem. 2011, 32, 2319– 2327, DOI: 10.1002/jcc.2178740https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXnvFalsr8%253D&md5=d567042c65cfdc1c81336a29137654bfMDAnalysis: A toolkit for the analysis of molecular dynamics simulationsMichaud-Agrawal, Naveen; Denning, Elizabeth J.; Woolf, Thomas B.; Beckstein, OliverJournal of Computational Chemistry (2011), 32 (10), 2319-2327CODEN: JCCHDD; ISSN:0192-8651. (John Wiley & Sons, Inc.)MDAnal. is an object-oriented library for structural and temporal anal. of mol. dynamics (MD) simulation trajectories and individual protein structures. It is written in the Python language with some performance-crit. code in C. It uses the powerful NumPy package to expose trajectory data as fast and efficient NumPy arrays. It has been tested on systems of millions of particles. Many common file formats of simulation packages including CHARMM, Gromacs, Amber, and NAMD and the Protein Data Bank format can be read and written. Atoms can be selected with a syntax similar to CHARMM's powerful selection commands. MDAnal. enables both novice and experienced programmers to rapidly write their own anal. tools and access data stored in trajectories in an easily accessible manner that facilitates interactive explorative anal. MDAnal. has been tested on and works for most Unix-based platforms such as Linux and Mac OS X. It is freely available under the GNU General Public License from http://mdanal.googlecode.com. © 2011 Wiley Periodicals, Inc. J Comput Chem 2011.
- 41North, R. A. Molecular physiology of P2X receptors. Physiol. Rev. 2002, 82, 1013– 1067, DOI: 10.1152/physrev.00015.200241https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD38Xot1Ggu78%253D&md5=d621c6ce761f5e2db7ec553ed0f5ac10Molecular physiology of P2X receptorsNorth, R. AlanPhysiological Reviews (2002), 82 (4), 1013-1067CODEN: PHREA7; ISSN:0031-9333. (American Physiological Society)A review. P2X receptors are membrane ion channels that open in response to the binding of extracellular ATP. Seven genes in vertebrates encode P2X receptor subunits, which are 40-50% identical in amino acid sequence. Each subunit has two transmembrane domains, sepd. by an extracellular domain (∼280 amino acids). Channels form as multimers of several subunits. Homomeric P2X1, P2X2, P2X3, P2X4, P2X5, and P2X7 channels and heteromeric P2X2/3 and P2X1/5 channels have been most fully characterized following heterologous expression. Some agonists (e.g., αβ-methylene ATP) and antagonists [e.g., 2',3'-O-(2,4,6-trinitrophenyl)-ATP] are strongly selective for receptors contg. P2X1 and P2X3 subunits. All P2X receptors are permeable to small monovalent cations; some have significant calcium or anion permeability. In many cells, activation of homomeric P2X7 receptors induces a permeability increase to larger org. cations including some fluorescent dyes and also signals to the cytoskeleton; these changes probably involve addnl. interacting proteins. P2X receptors are abundantly distributed, and functional responses are seen in neurons, glia, epithelia, endothelia, bone, muscle, and hemopoietic tissues. The mol. compn. of native receptors is becoming understood, and some cells express more than one type of P2X receptor. On smooth muscles, P2X receptors respond to ATP released from sympathetic motor nerves (e.g., in ejaculation). On sensory nerves, they are involved in the initiation of afferent signals in several viscera (e.g., bladder, intestine) and play a key role in sensing tissue-damaging and inflammatory stimuli. Paracrine roles for ATP signaling through P2X receptors are likely in neurohypophysis, ducted glands, airway epithelia, kidney, bone, and hemopoietic tissues. In the last case, P2X7 receptor activation stimulates cytokine release by engaging intracellular signaling pathways.
- 42Surprenant, A.; North, R. A. Signaling at purinergic P2X receptors. Annu. Rev. Physiol. 2009, 71, 333– 359, DOI: 10.1146/annurev.physiol.70.113006.10063042https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXjsValsbk%253D&md5=c9f3772f8f746ea0591ff30a3c7beb36Signaling at purinergic P2X receptorsSurprenant, Annmarie; North, R. AlanAnnual Review of Physiology (2009), 71 (), 333-359CODEN: ARPHAD; ISSN:0066-4278. (Annual Reviews Inc.)A review. P2X receptors are membrane cation channels gated by extracellular ATP. Seven P2X receptor subunits (P2X1-7) are widely distributed in excitable and nonexcitable cells of vertebrates. They play key roles in inter alia afferent signaling (including pain), regulation of renal blood flow, vascular endothelium, and inflammatory responses. We summarize the evidence for these and other roles, emphasizing exptl. work with selective receptor antagonists or with knockout mice. The receptors are trimeric membrane proteins: studies of the biophys. properties of mutated subunits expressed in heterologous cells have indicated parts of the subunits involved in ATP binding, ion permeation (including calcium permeability), and membrane trafficking. We review our current understanding of the mol. properties of P2X receptors, including how this understanding is informed by the identification of distantly related P2X receptors in simple eukaryotes.
- 43Illes, P.; Müller, C. E.; Jacobson, K. A.; Grutter, T.; Nicke, A.; Fountain, S. J.; Kennedy, C.; Schmalzing, G.; Jarvis, M. F.; Stojilkovic, S. S. Update of P2X receptor properties and their pharmacology: IUPHAR Review 30. Br. J. Pharmacol. 2021, 178, 489– 514, DOI: 10.1111/bph.1529943https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3MXktFehuw%253D%253D&md5=6cb4be20f4b37e5d18e523bdf7f1aef3Update of P2X receptor properties and their pharmacology: IUPHAR Review 30Illes, Peter; Mueller, Christa E.; Jacobson, Kenneth A.; Grutter, Thomas; Nicke, Annette; Fountain, Samuel J.; Kennedy, Charles; Schmalzing, Guenther; Jarvis, Michael F.; Stojilkovic, Stanko S.; King, Brian F.; Di Virgilio, FrancescoBritish Journal of Pharmacology (2021), 178 (3), 489-514CODEN: BJPCBM; ISSN:1476-5381. (Wiley-Blackwell)A review. The known seven mammalian receptor subunits (P2X1-7) form cationic channels gated by ATP. Three subunits compose a receptor channel. Each subunit is a polypeptide consisting of two transmembrane regions (TM1 and TM2), intracellular N- and C-termini, and a bulky extracellular loop. Crystn. allowed the identification of the 3D structure and gating cycle of P2X receptors. The agonist-binding pocket is located at the intersection of two neighboring subunits. In addn. to the mammalian P2X receptors, their primitive ligand-gated counterparts with little structural similarity have also been cloned. Selective agonists for P2X receptor subtypes are not available, but medicinal chem. supplied a range of subtype-selective antagonists, as well as pos. and neg. allosteric modulators. Knockout mice and selective antagonists helped to identify pathol. functions due to defective P2X receptors, such as male infertility (P2X1), hearing loss (P2X2), pain/cough (P2X3), neuropathic pain (P2X4), inflammatory bone loss (P2X5), and faulty immune reactions (P2X7).
- 44Briones, R.; Blau, C.; Kutzner, C.; de Groot, B. L.; Aponte-Santamaría, C. Gromaρs: A Gromacs-Based Toolset to Analyse Density Maps Derived from Molecular Dynamics Simulations. Biophys. J. 2019, 116, 4– 11, DOI: 10.1016/j.bpj.2018.11.312644https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXisVKitrbK&md5=7e63b96203c0e8aa9a182c4098310eb5GROmaρs: A GROMACS-Based Toolset to Analyze Density Maps Derived from Molecular Dynamics SimulationsBriones, Rodolfo; Blau, Christian; Kutzner, Carsten; de Groot, Bert L.; Aponte-Santamaria, CamiloBiophysical Journal (2019), 116 (1), 4-11CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)The authors introduce a computational toolset, named GROmaρs, to obtain and compare time-averaged d. maps from mol. dynamics simulations. GROmaρs efficiently computes d. maps by fast multi-Gaussian spreading of at. densities onto a three-dimensional grid. It complements existing map-based tools by enabling spatial inspection of at. av. localization during the simulations. Most importantly, it allows the comparison between computed and ref. maps (e.g., exptl.) through calcn. of difference maps and local and time-resolved global correlation. These comparison operations proved useful to quant. contrast perturbed and control simulation data sets and to examine how much biomol. systems resemble both synthetic and exptl. d. maps. This was esp. advantageous for multimol. systems in which std. comparisons like RMSDs are difficult to compute. In addn., GROmaρs incorporates abs. and relative spatial free-energy ests. to provide an energetic picture of atomistic localization. This is an open-source GROMACS-based toolset, thus allowing for static or dynamic selection of atoms or even coarse-grained beads for the d. calcn. Furthermore, masking of regions was implemented to speed up calcns. and to facilitate the comparison with exptl. maps. Beyond map comparison, GROmaρs provides a straightforward method to detect solvent cavities and av. charge distribution in biomol. systems. The authors employed all these functionalities to inspect the localization of lipid and water mols. in aquaporin systems, the binding of cholesterol to the G protein coupled chemokine receptor type 4, and the identification of permeation pathways through the dermicidin antimicrobial channel. Based on these examples, the authors anticipate a high applicability of GROmaρs for the anal. of mol. dynamics simulations and their comparison with exptl. detd. densities.
- 45Chen, X.; Bründl, M.; Friesacher, T.; Stary-Weinzinger, A. Computational Insights Into Voltage Dependence of Polyamine Block in a Strong Inwardly Rectifying K+ Channel. Front. Pharmacol. 2020, 11, 721 DOI: 10.3389/fphar.2020.0072145https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BB3cXhvFKktb7J&md5=e735c85eef42803dfb799d1736d99774Computational insights into voltage dependence of polyamine block in a strong inwardly rectifying K+ channelChen, Xingyu; Bruendl, Michael; Friesacher, Theres; Stary-Weinzinger, AnnaFrontiers in Pharmacology (2020), 11 (), 00721CODEN: FPRHAU; ISSN:1663-9812. (Frontiers Media S.A.)Inwardly rectifying potassium (KIR) channels play important roles in controlling cellular excitability and K+ ion homeostasis. Under physiol. conditions, KIR channels allow large K+ influx at potentials neg. to the equil. potential of K+ but permit little outward current at potentials pos. to the equil. potential of K+, due to voltage dependent block of outward K+ flux by cytoplasmic polyamines. These polycationic mols. enter the KIR channel pore from the intracellular side. They block K+ ion movement through the channel at depolarized potentials, thereby ensuring, for instance, the long plateau phase of the cardiac action potential. Key questions concerning how deeply these charged mols. migrate into the pore and how the steep voltage dependence arises remain unclear. Recent MD simulations on GIRK2 (=Kir3.2) crystal structures have provided unprecedented details concerning the conduction mechanism of a KIR channel. Here, we use MD simulations with applied field to provide detailed insights into voltage dependent block of putrescine, using the conductive state of the strong inwardly rectifying K+ channel GIRK2 as starting point. Our μs long simulations elucidate details about binding sites of putrescine in the pore and suggest that voltage dependent rectification arises from a dual mechanism.
- 46Nyblom, M.; Frick, A.; Wang, Y.; Ekvall, M.; Hallgren, K.; Hedfalk, K.; Neutze, R.; Tajkhorshid, E.; Törnroth-Horsefield, S. Structural and functional analysis of SoPIP2;1 mutants adds insight into plant aquaporin gating. J. Mol. Biol. 2009, 387, 653– 668, DOI: 10.1016/j.jmb.2009.01.06546https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXjsVOnu7c%253D&md5=7176199fe71199e1ea750db551d14c46Structural and Functional Analysis of SoPIP2;1 Mutants Adds Insight into Plant Aquaporin GatingNyblom, Maria; Frick, Anna; Wang, Yi; Ekvall, Mikael; Hallgren, Karin; Hedfalk, Kristina; Neutze, Richard; Tajkhorshid, Emad; Toernroth-Horsefield, SusannaJournal of Molecular Biology (2009), 387 (3), 653-668CODEN: JMOBAK; ISSN:0022-2836. (Elsevier Ltd.)Plant plasma membrane aquaporins facilitate water flux into and out of plant cells, thus coupling their cellular function to basic aspects of plant physiol. Posttranslational modifications of conserved phosphorylation sites, changes in cytoplasmic pH, and the binding of Ca2+ can regulate water transport activity by gating the plasma membrane aquaporins. A structural mechanism unifying these diverse biochem. signals has emerged for the spinach aquaporin SoPIP2;1, although several questions concerning the opening mechanism remain. Here, we describe the X-ray structures of the S115E and S274E single SoPIP2;1 mutants and the corresponding double mutant. Phosphorylation of these serines is believed to increase water transport activity of SoPIP2;1 by opening the channel. However, all mutants crystd. in a closed conformation, as confirmed by water transport assays, implying that neither substitution fully mimics the phosphorylated state. Nevertheless, a half-turn extension of transmembrane helix 1 occurs upon the substitution of Ser115, which draws the Cα atom of Glu31 10 Å away from its wild-type conformation, thereby disrupting the divalent cation binding site involved in the gating mechanism. Mutation of Ser274 disorders the C-terminus but no other significant conformational changes are obsd. Inspection of the hydrogen-bond interactions within loop D suggested that the phosphorylation of Ser188 may also produce an open channel, and this was supported by an increased water transport activity for the S188E mutant and mol. dynamics simulations. These findings add addnl. insight into the general mechanism of plant aquaporin gating.
- 47Ruiz, L.; Wu, Y.; Keten, S. Tailoring the water structure and transport in nanotubes with tunable interiors. Nanoscale 2015, 7, 121– 132, DOI: 10.1039/C4NR05407E47https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXhslGisr%252FI&md5=9224107648d21b753e8d61e94101d943Tailoring the water structure and transport in nanotubes with tunable interiorsRuiz, Luis; Wu, Yuanqiao; Keten, SinanNanoscale (2015), 7 (1), 121-132CODEN: NANOHL; ISSN:2040-3372. (Royal Society of Chemistry)Self-assembly of cyclic peptide nanotubes (CPNs) in polymer thin films has opened up the possibility of creating sepn. membranes with tunable nanopores that can differentiate mols. at the sub-nanometer level. While it has been demonstrated that the interior chem. of the CPNs can be tailored by inserting functional groups in the nanopore lumen (mCPNs), a design strategy for picking the chem. modifications that lead to particular transport properties has not been established. Drawing from the knowledgebase of functional groups in natural amino acids, here we use mol. dynamics simulations to elucidate how bioinspired mutations influence the transport of water through mCPNs. We show that, at the nanoscale, factors besides the pore size, such as electrostatic interactions and steric effects, can dramatically change the transport properties. We recognize a novel asym. structure of water under nanoconfinement inside the chem. functionalized nanotubes and identify that the small non-polar glycine-mimic groups that minimize the steric constraints and confer a hydrophobic character to the nanotube interior are the fastest transporters of water. Our computationally developed expts. on a realistic material system circumvent synthetic challenges, and lay the foundation for bioinspired principles to tailor artificial nanochannels for sepn. applications such as desalination, ion-exchange and carbon capture.
- 48Kasimova, M. A.; Lindahl, E.; Delemotte, L. Determining the molecular basis of voltage sensitivity in membrane proteins. J. Gen. Physiol. 2018, 150, 1444– 1458, DOI: 10.1085/jgp.20181208648https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXmtlSrsL4%253D&md5=467ecade0184d6cc82276971a54ce988Determining the molecular basis of voltage sensitivity in membrane proteinsKasimova, Marina A.; Lindahl, Erik; Delemotte, LucieJournal of General Physiology (2018), 150 (10), 1444-1458CODEN: JGPLAD; ISSN:1540-7748. (Rockefeller University Press)Voltage-sensitive membrane proteins are united by their ability to transform changes in membrane potential into mech. work. They are responsible for a spectrum of physiol. processes in living organisms, including elec. signaling and cell-cycle progression. Moreover, the detection of these elements by using exptl. techniques is challenging because of the diversity of membrane proteins. Here, we provide a computational approach to predict voltage-sensing elements in any membrane protein, independent of its structure or function. It relies on an estn. of the propensity of a protein to respond to changes in membrane potential. We first show that this property correlates well with voltage sensitivity by applying our approach to a set of voltage-sensitive and voltage-insensitive membrane proteins. We further show that it correctly identifies authentic voltage-sensitive residues in the voltage-sensor domain of voltage-gated ion channels. Finally, we investigate six membrane proteins for which the voltage-sensing elements have not yet been characterized and identify residues and ions that might be involved in the response to voltage. The suggested approach is fast and simple and enables a characterization of voltage sensitivity that goes beyond mere identification of charges. We anticipate that its application before mutagenesis expts. will significantly reduce the no. of potential voltage-sensitive elements to be tested.
- 49Scheurer, M.; Rodenkirch, P.; Siggel, M.; Bernardi, R. C.; Schulten, K.; Tajkhorshid, E.; Rudack, T. PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations. Biophys. J. 2018, 114, 577– 583, DOI: 10.1016/j.bpj.2017.12.00349https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXis1Oguw%253D%253D&md5=cdcf9bda51716dae11aa7575c13abb61PyContact: Rapid, Customizable, and Visual Analysis of Noncovalent Interactions in MD SimulationsScheurer, Maximilian; Rodenkirch, Peter; Siggel, Marc; Bernardi, Rafael C.; Schulten, Klaus; Tajkhorshid, Emad; Rudack, TillBiophysical Journal (2018), 114 (3), 577-583CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)Mol. dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth gave a large amt. of simulation data that need to be filtered for relevant information to address specific biomedical and biochem. questions. One of the most relevant mol. properties that can be studied by all-atom MD simulations is the time-dependent evolution of the complex noncovalent interaction networks governing such fundamental aspects as mol. recognition, binding strength, and mech. and structural stability. Extg., evaluating, and visualizing noncovalent interactions is a key task in the daily work of structural biologists. The authors have developed PyContact, an easy-to-use, highly flexible, and intuitive graphical user interface-based application, designed to provide a toolkit to study biomol. interactions in MD trajectories. PyContact is designed to facilitate this task by enabling identification of relevant noncovalent interactions in a comprehensible manner. The implementation of PyContact as a standalone application enables rapid anal. and data visualization without any addnl. programming requirements, and also preserves full in-program customization and extension capabilities for advanced users. The statistical anal. representation is interactively combined with full mapping of the results on the mol. system through the synergistic connection between PyContact and VMD. The authors showcase the capabilities and scientific significance of PyContact by analyzing and visualizing in great detail the noncovalent interactions underlying the ion permeation pathway of the human P2X3 receptor. As a second application, the authors examine the protein-protein interaction network of the mech. ultrastable cohesin-dockering complex.
- 50Melcr, J.; Martinez-Seara, H.; Nencini, R.; Kolafa, J.; Jungwirth, P.; Ollila, O. H. S. Accurate Binding of Sodium and Calcium to a POPC Bilayer by Effective Inclusion of Electronic Polarization. J. Phys. Chem. B 2018, 122, 4546– 4557, DOI: 10.1021/acs.jpcb.7b1251050https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXmvVaqtro%253D&md5=a49868ee725f486dcc49522bca8038d0Accurate Binding of Sodium and Calcium to a POPC Bilayer by Effective Inclusion of Electronic PolarizationMelcr, Josef; Martinez-Seara, Hector; Nencini, Ricky; Kolafa, Jiri; Jungwirth, Pavel; Ollila, O. H. SamuliJournal of Physical Chemistry B (2018), 122 (16), 4546-4557CODEN: JPCBFK; ISSN:1520-5207. (American Chemical Society)Binding affinities and stoichiometries of Na+ and Ca2+ ions to phospholipid bilayers are of paramount significance in the properties and functionality of cellular membranes. Current ests. of binding affinities and stoichiometries of cations are, however, inconsistent due to limitations in the available exptl. and computational methods. In this work, the authors improve the description of the binding details of Na+ and Ca2+ ions to a 1-palmitoyl-2-oleoylphosphatidylcholine (POPC) bilayer by implicitly including electronic polarization as a mean field correction, known as the electronic continuum correction (ECC). This is applied by scaling the partial charges of a selected state-of-the-art POPC lipid model for mol. dynamics simulations. The improved ECC-POPC model reproduces not only the exptl. measured structural parameters for the ion-free membrane, but also the response of lipid headgroup to a strongly bound cationic amphiphile, as well as the binding affinities of Na+ and Ca2+ ions. With the new model, the authors observe on the one side negligible binding of Na+ ions to POPC bilayer, while on the other side stronger interactions of Ca2+ primarily with phosphate oxygens, which is in agreement with the previous interpretations of the exptl. spectroscopic data. The present model results in Ca2+ ions forming complexes with one to three POPC mols. with almost equal probabilities, suggesting more complex binding stoichiometries than those from simple models used to interpret the NMR data previously. The results of this work pave the way to quant. mol. simulations with realistic electrostatic interactions of complex biochem. systems at cellular membranes.
- 51Melcrová, A.; Pokorna, S.; Pullanchery, S.; Kohagen, M.; Jurkiewicz, P.; Hof, M.; Jungwirth, P.; Cremer, P. S.; Cwiklik, L. The complex nature of calcium cation interactions with phospholipid bilayers. Sci. Rep. 2016, 6, 38035 DOI: 10.1038/srep3803551https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XitFWnsr7F&md5=a4be8d96f3217f88b2ab791199726e08The complex nature of calcium cation interactions with phospholipid bilayersMelcrova, Adela; Pokorna, Sarka; Pullanchery, Saranya; Kohagen, Miriam; Jurkiewicz, Piotr; Hof, Martin; Jungwirth, Pavel; Cremer, Paul S.; Cwiklik, LukaszScientific Reports (2016), 6 (), 38035CODEN: SRCEC3; ISSN:2045-2322. (Nature Publishing Group)Understanding interactions of calcium with lipid membranes at the mol. level is of great importance in light of their involvement in calcium signaling, assocn. of proteins with cellular membranes, and membrane fusion. We quantify these interactions in detail by employing a combination of spectroscopic methods with atomistic mol. dynamics simulations. Namely, time-resolved fluorescent spectroscopy of lipid vesicles and vibrational sum frequency spectroscopy of lipid monolayers are used to characterize local binding sites of calcium in zwitterionic and anionic model lipid assemblies, while dynamic light scattering and zeta potential measurements are employed for macroscopic characterization of lipid vesicles in calcium-contg. environments. To gain addnl. at.-level information, the expts. are complemented by mol. simulations that utilize an accurate force field for calcium ions with scaled charges effectively accounting for electronic polarization effects. We demonstrate that lipid membranes have substantial calcium-binding capacity, with several types of binding sites present. Significantly, the binding mode depends on calcium concn. with important implications for calcium buffering, synaptic plasticity, and protein-membrane assocn.
- 52Burnstock, G. Purinergic signalling and disorders of the central nervous system. Nat. Rev. Drug Discovery 2008, 7, 575– 590, DOI: 10.1038/nrd260552https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1cXnvFShtrs%253D&md5=95a3c406c113eec093fd5437e98ee49ePurinergic signaling and disorders of the central nervous systemBurnstock, GeoffreyNature Reviews Drug Discovery (2008), 7 (7), 575-590CODEN: NRDDAG; ISSN:1474-1776. (Nature Publishing Group)A review. There is increasing interest in the role of purinergic signaling in CNS disorders, but few modulators have reached the clinic. Here, Burnstock provides an overview of the role of purinergic signaling in specific disorders, including brain trauma, ischemia, neurodegenerative, neuroinflammatory and neuropsychiatric diseases, and highlights specific purinergic receptor subtypes that represent promising therapeutic targets. Purines have key roles in neurotransmission and neuromodulation, with their effects being mediated by the purine and pyrimidine receptor subfamilies, P1, P2X and P2Y. Recently, purinergic mechanisms and specific receptor subtypes have been shown to be involved in various pathol. conditions including brain trauma and ischemia, neurodegenerative diseases involving neuroimmune and neuroinflammatory reactions, as well as in neuropsychiatric diseases, including depression and schizophrenia. This article reviews the role of purinergic signaling in CNS disorders, highlighting specific purinergic receptor subtypes, most notably A2A, P2X4 and P2X7, that might be therapeutically targeted for the treatment of these conditions.
- 53Pomorski, T.; Menon, A. K. Lipid flippases and their biological functions. Cell. Mol. Life Sci. 2006, 63, 2908– 2921, DOI: 10.1007/s00018-006-6167-753https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD2sXht1egtb0%253D&md5=10e9fa64e238c152ee43de250d62a07eLipid flippases and their biological functionsPomorski, T.; Menon, A. K.Cellular and Molecular Life Sciences (2006), 63 (24), 2908-2921CODEN: CMLSFI; ISSN:1420-682X. (Birkhaeuser Verlag)A review. The typically distinct phospholipid compn. of the 2 leaflets of a membrane bilayer is generated and maintained by bi-directional transport (flip-flop) of lipids between the leaflets. Specific membrane proteins, termed lipid flippases, play an essential role in this transport process. Energy-independent flippases allow common phospholipids to equilibrate rapidly between the two monolayers and also play a role in the biosynthesis of a variety of glycoconjugates such as glycosphingolipids, N-glycoproteins, and glycosylphosphatidylinositol (GPI)-anchored proteins. ATP-dependent flippases, including members of a conserved subfamily of P-type ATPases and ATP-binding cassette transporters, mediate the net transfer of specific phospholipids to one leaflet of a membrane and are involved in the creation and maintenance of transbilayer lipid asymmetry of membranes such as the plasma membrane of eukaryotes. Energy-dependent flippases also play a role in the biosynthesis of glycoconjugates such as bacterial lipopolysaccharide. Here, the authors summarize recent progress on the identification and characterization of the various flippases and the demonstration of their biol. functions.
- 54Jiang, T.; Yu, K.; Hartzell, H. C.; Tajkhorshid, E. Lipids and ions traverse the membrane by the same physical pathway in the nhTMEM16 scramblase. eLife 2017, 6, e28671 DOI: 10.7554/eLife.2867154https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXitlOit7nJ&md5=e124dd1a49738ae073b1abca3c95ef2eLipids and ions traverse the membrane by the same physical pathway in the nhTMEM16 scramblaseJiang, Tao; Yu, Kuai; Hartzell, H. Criss; Tajkhorshid, EmadeLife (2017), 6 (), e28671/1-e28671/26CODEN: ELIFA8; ISSN:2050-084X. (eLife Sciences Publications Ltd.)From bacteria to mammals, different phospholipid species are segregated between the inner and outer leaflets of the plasma membrane by ATP-dependent lipid transporters. Disruption of this asymmetry by ATP-independent phospholipid scrambling is important in cellular signaling, but its mechanism remains incompletely understood. Using MD simulations coupled with exptl. assays, we show that the surface hydrophilic transmembrane cavity exposed to the lipid bilayer on the fungal scramblase nhTMEM16 serves as the pathway for both lipid translocation and ion conduction across the membrane. Ca2+ binding stimulates its open conformation by altering the structure of transmembrane helixes that line the cavity. We have identified key amino acids necessary for phospholipid scrambling and validated the idea that ions permeate TMEM16 Cl- channels via a structurally homologous pathway by showing that mutation of two residues in the pore region of the TMEM16A Ca2+-activated Cl- channel convert it into a robust scramblase.
- 55Bethel, N. P.; Grabe, M. Atomistic insight into lipid translocation by a TMEM16 scramblase. Proc. Natl. Acad. Sci. U.S.A. 2016, 113, 14049– 14054, DOI: 10.1073/pnas.160757411355https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhvV2gu7rI&md5=c73f9572ab589dba5374c160eca1b9b6Atomistic insight into lipid translocation by a TMEM16 scramblaseBethel, Neville P.; Grabe, MichaelProceedings of the National Academy of Sciences of the United States of America (2016), 113 (49), 14049-14054CODEN: PNASA6; ISSN:0027-8424. (National Academy of Sciences)The transmembrane protein 16 (TMEM16) family of membrane proteins includes both lipid scramblases and ion channels involved in olfaction, nociception, and blood coagulation. The crystal structure of the fungal Nectria haematococca TMEM16 (nhTMEM16) scramblase suggested a putative mechanism of lipid transport, whereby polar and charged lipid headgroups move through the low-dielec. environment of the membrane by traversing a hydrophilic groove on the membrane-spanning surface of the protein. Here, we use computational methods to explore the membrane-protein interactions involved in lipid scrambling. Fast, continuum membrane-bending calcns. reveal a global pattern of charged and hydrophobic surface residues that bends the membrane in a large-amplitude sinusoidal wave, resulting in bilayer thinning across the hydrophilic groove. Atomic simulations uncover two lipid headgroup-interaction sites flanking the groove. The cytoplasmic site nucleates headgroup-dipole stacking interactions that form a chain of lipid mols. that penetrate into the groove. In two instances, a cytoplasmic lipid interdigitates into this chain, crosses the bilayer, and enters the extracellular leaflet, and the reverse process happens twice as well. Continuum membrane-bending anal. carried out on homol. models of mammalian homologs shows that these family members also bend the membrane-even those that lack scramblase activity. Sequence alignments show that the lipid-interaction sites are conserved in many family members but less so in those with reduced scrambling ability. Our anal. provides insight into how large-scale membrane bending and protein chem. facilitate lipid permeation in the TMEM16 family, and we hypothesize that membrane interactions also affect ion permeation.
- 56Pedemonte, N.; Galietta, L. J. V. Structure and function of TMEM16 proteins (anoctamins). Physiol. Rev. 2014, 94, 419– 459, DOI: 10.1152/physrev.00039.201156https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2cXptFygt70%253D&md5=9b0cf5bb181dc1ce259d3cab7e7de629Structure and function of TMEM16 proteins (anoctamins)Pedemonte, Nicoletta; Galietta, Luis J. V.Physiological Reviews (2014), 94 (2), 419-459CODEN: PHREA7; ISSN:0031-9333. (American Physiological Society)A review. TMEM16 proteins, also known as anoctamins, are involved in a variety of functions that include ion transport, phospholipid scrambling, and regulation of other membrane proteins. The first two members of the family, TMEM16A (anoctamin-1, AN01) and TMEM16B (anoctamin-2, AN02), function as Ca2+-activated Cl- channels (CaCCs), a type of ion channel that plays important functions such as transepithelial ion transport, smooth muscle contraction, olfaction, phototransduction, nociception, and control of neuronal excitability. Genetic ablation of TMEM16A in mice causes impairment of epithelial Cl- secretion, tracheal abnormalities, and block of gastrointestinal peristalsis. TMEM16A is directly regulated by cytosolic Ca2+ as well as indirectly by its interaction with calmodulin. Other members of the anoctamin family, such as TMEM16C, TMEM16D, TMEM16F, TMEM16G, and TMEM16J, may work as phospholipid scramblases and/or ion channels. In particular, TMEM16F (AN06) is a major contributor to the process of phosphatidylserine translocation from the inner to the outer leaflet of the plasma membrane. Intriguingly, TMEM16F is also assocd. with the appearance of anion/cation channels activated by very high Ca2+ concns. Furthermore, a TMEM16 protein expressed in Aspergillus fumigatus displays both ion channel and lipid scramblase activity. This finding suggests that dual function is an ancestral characteristic of TMEM16 proteins and that some members, such as TMEM16A and TMEM16B, have evolved to a pure channel function. Mutations in anoctamin genes (AN03, AN05, AN06, and AN010) cause various genetic diseases. These diseases suggest the involvement of anoctamins in a variety of cell functions whose link with ion transport and/or lipid scrambling needs to be clarified.
- 57Kalienkova, V.; Clerico Mosina, V.; Bryner, L.; Oostergetel, G. T.; Dutzler, R.; Paulino, C. Stepwise activation mechanism of the scramblase nhTMEM16 revealed by cryo-EM. eLife 2019, 8, e44364 DOI: 10.7554/eLife.4436457https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXisVegtb3K&md5=fd6a7da5a78a87f0439578ee967f76dbStepwise activation mechanism of the scramblase nhTMEM16 revealed by cryo- EMKalienkova, Valeria; Mosina, Vanessa Clerico; Bryner, Laura; Oostergetel, Gert T.; Dutzler, Raimund; Paulino, CristinaeLife (2019), 8 (), 44364CODEN: ELIFA8; ISSN:2050-084X. (eLife Sciences Publications Ltd.)Scramblases catalyze the movement of lipids between both leaflets of a bilayer. Whereas the X-ray structure of the protein nhTMEM16 has previously revealed the architecture of a Ca2+-dependent lipid scramblase, its regulation mechanism has remained elusive. Here, we have used cryo-electron microscopy and functional assays to address this question. Ca2+-bound and Ca2+-free conformations of nhTMEM16 in detergent and lipid nanodiscs illustrate the interactions with its environment and they reveal the conformational changes underlying its activation. In this process, Ca2+ binding induces a stepwise transition of the catalytic subunit cavity, converting a closed cavity that is shielded from the membrane in the absence of ligand, into a polar furrow that becomes accessible to lipid headgroups in the Ca2+-bound state. Addnl., our structures demonstrate how nhTMEM16 distorts the membrane at both entrances of the subunit cavity, thereby decreasing the energy barrier for lipid movement.
- 58Falzone, M. E.; Rheinberger, J.; Lee, B.-C.; Peyear, T.; Sasset, L.; Raczkowski, A. M.; Eng, E. T.; Di Lorenzo, A.; Andersen, O. S.; Nimigean, C. M.; Accardi, A. Structural basis of Ca2+-dependent activation and lipid transport by a TMEM16 scramblase. eLife 2019, 8, e43229 DOI: 10.7554/eLife.4322958https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXisVegtLvK&md5=bf21156b76417f26aa25c95a6203cfa2Structural basis of Ca2+-dependent activation and lipid transport by a TMEM16 scramblaseFalzone, Maria E.; Rheinberger, Jan; Lee, Byoung-Cheol; Peyear, Thasin; Sasset, Linda; Raczkowski, Ashleigh M.; Eng, Edward T.; Lorenzo, Annarita Di; Andersen, Olaf S.; Nimigean, Crina M.; Accardi, AlessioeLife (2019), 8 (), 43229CODEN: ELIFA8; ISSN:2050-084X. (eLife Sciences Publications Ltd.)The lipid distribution of plasma membranes of eukaryotic cells is asym. and phospholipid scramblases disrupt this asymmetry by mediating the rapid, nonselective transport of lipids down their concn. gradients. As a result, phosphatidylserine is exposed to the outer leaflet of membrane, an important step in extracellular signaling networks controlling processes such as apoptosis, blood coagulation, membrane fusion and repair. Several TMEM16 family members have been identified as Ca2+-activated scramblases, but the mechanisms underlying their Ca2+-dependent gating and their effects on the surrounding lipid bilayer remain poorly understood. Here, we describe three high-resoln. cryo-electron microscopy structures of a fungal scramblase from Aspergillus fumigatus, afTMEM16, reconstituted in lipid nanodiscs. These structures reveal that Ca2+-dependent activation of the scramblase entails global rearrangement of the transmembrane and cytosolic domains. These structures, together with functional expts., suggest that activation of the protein thins the membrane near the transport pathway to facilitate rapid transbilayer lipid movement.
- 59Feng, S.; Dang, S.; Han, T. W.; Ye, W.; Jin, P.; Cheng, T.; Li, J.; Jan, Y. N.; Jan, L. Y.; Cheng, Y. Cryo-EM studies of TMEM16F Calcium-Activated ion channel suggest features important for lipid scrambling. Cell Rep. 2019, 28, 567– 579, DOI: 10.1016/j.celrep.2019.06.02359https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhtlGjtrjI&md5=be585b9b50fc27815a33b0aae7de9bdbCryo-EM Studies of TMEM16F Calcium-Activated Ion Channel Suggest Features Important for Lipid ScramblingFeng, Shengjie; Dang, Shangyu; Han, Tina Wei; Ye, Wenlei; Jin, Peng; Cheng, Tong; Li, Junrui; Jan, Yuh Nung; Jan, Lily Yeh; Cheng, YifanCell Reports (2019), 28 (2), 567-579.e4CODEN: CREED8; ISSN:2211-1247. (Cell Press)As a Ca2+-activated lipid scramblase and ion channel that mediates Ca2+ influx, TMEM16F relies on both functions to facilitate extracellular vesicle generation, blood coagulation, and bone formation. How a bona fide ion channel scrambles lipids remains elusive. Our structural analyses revealed the coexistence of an intact channel pore and PIP2-dependent protein conformation changes leading to membrane distortion. Correlated to the extent of membrane distortion, many tightly bound lipids are slanted. Structure-based mutagenesis studies further reveal that neutralization of some lipid-binding residues or those near membrane distortion specifically alters the onset of lipid scrambling, but not Ca2+ influx, thus identifying features outside of channel pore that are important for lipid scrambling. Together, our studies demonstrate that membrane distortion does not require open hydrophilic grooves facing the membrane interior and provide further evidence to suggest sep. pathways for lipid scrambling and ion permeation.
- 60Alvadia, C.; Lim, N. K.; Clerico Mosina, V.; Oostergetel, G. T.; Dutzler, R.; Paulino, C. Cryo-EM structures and functional characterization of the murine lipid scramblase TMEM16F. eLife 2019, 8, e44365 DOI: 10.7554/eLife.4436560https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXisVegt7rL&md5=0cd86dbb69ebf0ad0da6c047b9efc032Cryo-EM structures and functional characterization of the murine lipid scramblase TMEM16FAlvadia, Carolina; Lim, Novandy K.; Mosina, Vanessa Clerico; Oostergetel, Gert T.; Dutzler, Raimund; Paulino, CristinaeLife (2019), 8 (), 44365CODEN: ELIFA8; ISSN:2050-084X. (eLife Sciences Publications Ltd.)The lipid scramblase TMEM16F initiates blood coagulation by catalyzing the exposure of phosphatidylserine in platelets. The protein is part of a family of membrane proteins, which encompasses calcium-activated channels for ions and lipids. Here, we reveal features of murine TMEM16F (mTMEM16F) that underlie its function as a lipid scramblase and an ion channel. The cryo-EM data of mTMEM16F in absence and presence of Ca2+ define the ligand-free closed conformation of the protein and the structure of a Ca2+-bound intermediate. Both conformations resemble their counterparts of the scrambling-incompetent anion channel mTMEM16A, yet with distinct differences in the region of ion and lipid permeation. In conjunction with functional data, we demonstrate the relationship between ion conduction and lipid scrambling. Although activated by a common mechanism, both functions appear to be mediated by alternate protein conformations that are at equil. in the ligand-bound state.
- 61Lim, N. K.; Lam, A. K.; Dutzler, R. Independent activation of ion conduction pores in the double-barreled calcium-activated chloride channel TMEM16A. J. Gen. Physiol. 2016, 148, 375– 392, DOI: 10.1085/jgp.20161165061https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXkvF2qtLk%253D&md5=4cab74ef8af397f48127717903dfcf78Independent activation of ion conduction pores in the double-barreled calcium-activated chloride channel TMEM16ALim, Novandy K.; Lam, Andy K. M.; Dutzler, RaimundJournal of General Physiology (2016), 148 (5), 375-392, S19-S20CODEN: JGPLAD; ISSN:1540-7748. (Rockefeller University Press)The TMEM16 proteins constitute a family of membrane proteins with unusual functional breadth, including lipid scramblases and Cl- channels. Members of both these branches are activated by Ca2+, acting from the intracellular side, and probably share a common architecture, which was defined in the recent structure of the lipid scramblase nhTMEM16. The structural features of subunits and the arrangement of Ca2+-binding sites in nhTMEM16 suggest that the dimeric protein harbors two locations for catalysis that are independent with respect to both activation and lipid conduction. Here, we ask whether a similar independence is obsd. in the Ca2+-activated Cl- channel TMEM16A. For this purpose, we generated concatenated constructs contg. subunits with distinct activation and permeation properties. Our biochem. investigations demonstrate the integrity of concatemers after solubilization and purifn. During investigation by patch-clamp electrophysiol., the functional behavior of constructs contg. either two wild-type (WT) subunits or one WT subunit paired with a second subunit with compromised activation closely resembles TMEM16A. This resemblance extends to ion selectivity, conductance, and the concn. and voltage dependence of channel activation by Ca2+. Constructs combining subunits with different potencies for Ca2+ show a biphasic activation curve that can be described as a linear combination of the properties of its constituents. The functional independence is further supported by mutation of a putative pore-lining residue that changes the conduction properties of the mutated subunit. Our results strongly suggest that TMEM16A contains two ion conduction pores that are independently activated by Ca2+ binding to sites that are embedded within the transmembrane part of each subunit.
- 62Jeng, G.; Aggarwal, M.; Yu, W.-P.; Chen, T.-Y. Independent activation of distinct pores in dimeric TMEM16A channels. J. Gen. Physiol. 2016, 148, 393– 404, DOI: 10.1085/jgp.20161165162https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXkvF2qtrk%253D&md5=72fcb349ae446a73a14d16604751294aIndependent activation of distinct pores in dimeric TMEM16A channelsJeng, Grace; Aggarwal, Muskaan; Yu, Wei-Ping; Chen, Tsung-YuJournal of General Physiology (2016), 148 (5), 393-404, S13-S14CODEN: JGPLAD; ISSN:1540-7748. (Rockefeller University Press)The TMEM16 family encompasses Ca2+-activated Cl- channels (CaCCs) and lipid scramblases. These proteins are formed by two identical subunits, as confirmed by the recently solved crystal structure of a TMEM16 lipid scram-blase. However, the high-resoln. structure did not provide definitive information regarding the pore architecture of the TMEM16 channels. In this study, we express TMEM16A channels constituting two covalently linked subunits with different Ca2+ affinities. The dose-response curve of the heterodimer appears to be a weighted sum of two dose-response curves-one corresponding to the high-affinity subunit and the other to the low-affinity subunit. However, fluorescence resonance energy transfer expts. suggest that the covalently linked heterodimeric proteins fold and assemble as one mol. Together these results suggest that activation of the two TMEM16A subunits likely activate independently of each other. The Ca2+ activation curve for the heterodi-mer at a low Ca2+ concn. range ([Ca2+] < 5μM) is similar to that of the wild-type channel-the Hill coeffs. in both cases are significantly greater than one. This suggests that Ca2+ binding to one subunit of TMEM16A is sufficient to activate the channel and that each subunit contains more than one Ca2+-binding site. We also take advantage of the I-V curve rectification that results from mutation of a pore residue to address the pore architec-ture of the channel. By introducing the pore mutation and the mutation that alters Ca2+ affinity in the same or different subunits, we demonstrate that activation of different subunits appears to be assocd. with the opening of different pores. These results suggest that the TMEM16A CaCC may also adopt a"double-barrel"pore architecture, similar to that found in CLC channels and transporters.
- 63Kanner, B. I.; Bendahan, A. Binding order of substrates to the (sodium + potassium ion)-coupled L-glutamic acid transporter from rat brain. Biochemistry 1982, 21, 6327– 6330, DOI: 10.1021/bi00267a04463https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaL38XlvFelsLo%253D&md5=2ff713407683ab142203cc789e4a3b61Binding order of substrates to the (sodium + potassium ion)-coupled L-glutamic acid transporter from rat brainKanner, Baruch I.; Bendahan, AnnieBiochemistry (1982), 21 (24), 6327-30CODEN: BICHAW; ISSN:0006-2960.Efflux of L-glutamic acid (I) [7440-23-5] from rat brain synaptic plasm membrane vesicles requires external K+. This requirement is satd. by concns. of ∼15 mequiv/L K+. In the absence of K+, I can be released from the vesicles in the presence of external I. This stimulation does not require external Na+ but is dependent on the external concn. of I. Half-maximal effects are obtained by concns. of ∼1 μM which are very similar to the apparent Km for I influx. Efflux of labeled I driven by external Na+ plus I requires internal Na+. The transporter apparently displays an asym. behavior toward Na+. This ion dissocs. much more slowly than does I on the external surface of the membrane but not on the internal surface. Furthermore, the transporter apparently translocates K+ in a step distinct from the I translocation step. The simplest explanation is that upon translocation of Na+ and I and their release to the inside, K+ binds to the transporter, enabling it to return to the outside to allow initiation of a new transport cycle.
- 64Zerangue, N.; Kavanaugh, M. P. Flux coupling in a neuronal glutamate transporter. Nature 1996, 383, 634– 637, DOI: 10.1038/383634a064https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADyaK28XmtlSqsbg%253D&md5=4d4328ffafe70486746bf0b8d3929a4fFlux coupling in a neuronal glutamate transporterZerangue, Noa; Kavanaugh, Michael P.Nature (London) (1996), 383 (6601), 634-637CODEN: NATUAS; ISSN:0028-0836. (Macmillan Magazines)Synaptic transmission is commonly terminated by diffusion and reuptake of neurotransmitter from the synaptic cleft. Glutamate reuptake prevents neurotoxicity and sets the lower limit for the concn. of extracellular glutamate, so it is important to understand the thermodn. of this process. Here we use voltage clamping with a pH-sensitive fluorescent dye to monitor elec. currents and pH changes assocd. with flux of glutamate mediated by the human neuronal glutamate transporter EAAT3. In contrast to a previous model, we find that three sodium ions and one proton are cotransported with each glutamate ion into the cell, while one potassium ion is transported out of the cell. This coupling can support a transmembrane glutamate concn. gradient ([Glu]in/[Glu]out) exceeding 106 under equil. conditions, and would allow the transporter to continue removing glutamate over a wide range of ionic conditions.
- 65Vandenberg, R. J.; Ryan, R. M. Mechanisms of Glutamate Transport. Physiol. Rev. 2013, 93, 1621– 1657, DOI: 10.1152/physrev.00007.201365https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3sXhvVShtrrP&md5=b573b96bbaf6f03c835c471cf2aaf7b8Mechanisms of glutamate transportVandenberg, Robert J.; Ryan, Renae M.Physiological Reviews (2013), 93 (4), 1621-1657CODEN: PHREA7; ISSN:0031-9333. (American Physiological Society)A review. L-Glutamate is the predominant excitatory neurotransmitter in the mammalian central nervous system and plays important roles in a wide variety of brain functions, but it is also a key player in the pathogenesis of many neurol. disorders. The control of glutamate concns. is crit. to the normal functioning of the central nervous system, and in this review we discuss how glutamate transporters regulate glutamate concns. to maintain dynamic signaling mechanisms between neurons. In 2004, the crystal structure of a prokaryotic homolog of the mammalian glutamate transporter family of proteins was crystd. and its structure detd. This has paved the way for a better understanding of the structural basis for glutamate transporter function. In this review we provide a broad perspective of this field of research, but focus primarily on the more recent studies with a particular emphasis on how our understanding of the structure of glutamate transporters has generated new insights.
- 66Machtens, J.-P.; Kortzak, D.; Lansche, C.; Leinenweber, A.; Kilian, P.; Begemann, B.; Zachariae, U.; Ewers, D.; de Groot, B. L.; Briones, R.; Fahlke, C. Mechanisms of anion conduction by coupled glutamate transporters. Cell 2015, 160, 542– 553, DOI: 10.1016/j.cell.2014.12.03566https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2MXitFyls7c%253D&md5=36f8d3049997b4b1877e4dbc13b6630fMechanisms of anion conduction by coupled glutamate transportersMachtens, Jan-Philipp; Kortzak, Daniel; Lansche, Christine; Leinenweber, Ariane; Kilian, Petra; Begemann, Birgit; Zachariae, Ulrich; Ewers, David; de Groot, Bert L.; Briones, Rodolfo; Fahlke, ChristophCell (Cambridge, MA, United States) (2015), 160 (3), 542-553CODEN: CELLB5; ISSN:0092-8674. (Cell Press)Excitatory amino acid transporters (EAATs) are essential for terminating glutamatergic synaptic transmission. They are not only coupled glutamate/Na+/H+/K+ transporters, but they also function as anion-selective channels. EAAT anion channels regulate neuronal excitability, and gain-of-function mutations in these proteins result in ataxia and epilepsy. Here, the authors combined mol. dynamics simulations with fluorescence spectroscopy of prokaryotic homolog glutamate transporter GltPh and patch-clamp recordings of mammalian EAATs to det. how these transporters conduct anions. Whereas outward- and inward-facing GltPh conformations were nonconductive, lateral movement of the glutamate transport domain from intermediate transporter conformations resulted in the formation of an anion-selective conduction pathway. Fluorescence quenching of inserted Trp residues indicated the entry of anions into this pathway, and mutations of homologous pore-forming residues had analogous effects on GltPh simulations and EAAT2/EAAT4 measurements of single-channel currents and anion/cation selectivities. These findings provided a mechanistic framework of how neurotransmitter transporters can operate as anion-selective and ligand-gated ion channels.
- 67Bastug, T.; Heinzelmann, G.; Kuyucak, S.; Salim, M.; Vandenberg, R. J.; Ryan, R. M. Position of the Third Na+ Site in the Aspartate Transporter GltPh and the Human Glutamate Transporter, EAAT1. PLoS One 2012, 7, e33058 DOI: 10.1371/journal.pone.003305867https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38Xks1Grs7Y%253D&md5=5b2a68466812e3ca7c2730b78dc4e0c0Position of the third Na+ site in the Aspartate transporter GltPh and the human glutamate transporter, EAAT1Bastug, Turgut; Heinzelmann, Germano; Kuyucak, Serdar; Salim, Marietta; Vandenberg, Robert J.; Ryan, Renae M.PLoS One (2012), 7 (3), e33058CODEN: POLNCL; ISSN:1932-6203. (Public Library of Science)Glutamate transport via the human excitatory amino acid transporters is coupled to the co-transport of three Na+ ions, one H+ and the counter-transport of one K+ ion. Transport by an archaeal homolog of the human glutamate transporters, GltPh, whose three dimensional structure is known is also coupled to three Na+ ions but only two Na+ ion binding sites have been obsd. in the crystal structure of GltPh. In order to fully utilize the GltPh structure in functional studies of the human glutamate transporters, it is essential to understand the transport mechanism of GltPh and accurately det. the no. and location of Na+ ions coupled to transport. Several sites have been proposed for the binding of a third Na+ ion from electrostatic calcns. and mol. dynamics simulations. In this study, we have performed detailed free energy simulations for GltPh and reveal a new site for the third Na+ ion involving the side chains of Threonine 92, Serine 93, Asparagine 310, Aspartate 312, and the backbone of Tyrosine 89. We have also studied the transport properties of alanine mutants of the coordinating residues Threonine 92 and Serine 93 in GltPh, and the corresponding residues in a human glutamate transporter, EAAT1. The mutant transporters have reduced affinity for Na+ compared to their wild type counterparts. These results confirm that Threonine 92 and Serine 93 are involved in the coordination of the third Na+ ion in GltPh and EAAT1.
- 68Kortzak, D.; Alleva, C.; Weyand, I.; Ewers, D.; Zimmermann, M. I.; Franzen, A.; Machtens, J.-P.; Fahlke, C. Allosteric gate modulation confers K+ coupling in glutamate transporters. EMBO J. 2019, 38, e101468 DOI: 10.15252/embj.2019101468There is no corresponding record for this reference.
- 69Setiadi, J.; Kuyucak, S. Free-Energy Simulations Resolve the Low-Affinity Na+-High-Affinity Asp Binding Paradox in GltPh. Biophys. J. 2019, 117, 780– 789, DOI: 10.1016/j.bpj.2019.07.01669https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1MXhsFSrs7nO&md5=730d8e0bc640faf763c5ebc5e004bcffFree-Energy Simulations Resolve the Low-Affinity Na+-High-Affinity Asp Binding Paradox in GltPhSetiadi, Jeffry; Kuyucak, SerdarBiophysical Journal (2019), 117 (4), 780-789CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)Glutamate transporters clear up excess extracellular glutamate by cotransporting three Na+ and one H+ with the countertransport of one K+. The archaeal homologs are selective to aspartate and only cotransport three Na+. The crystal structures of GltPh from archaea have been used in computational studies to understand the transport mechanism. Although some progress has been made with regard to the ligand-binding sites, a consistent picture of transport still eludes us. A major concern is the discrepancy between the computed binding free energies, which predict high-affinity Na+-low-affinity aspartate binding, and the exptl. results in which the opposite is obsd. Here, we show that the binding of the first two Na+ ions involves an intermediate state near the Na1 site, where two Na+ ions coexist and couple to aspartate with similar strengths, boosting its affinity. Binding free energies for Na+ and aspartate obtained using this intermediate state are in good agreement with the exptl. values. Thus, the paradox in binding affinities arises from the assumption that the ligands bind to the sites obsd. in the crystal structure following the order dictated by their binding free energies with no intermediate states. In fact, the presence of an intermediate state eliminates such a correlation between the binding free energies and the binding order. The intermediate state also facilitates transition of the first Na+ ion to its final binding site via a knock-on mechanism, which induces substantial conformational changes in the protein consistent with exptl. observations.
- 70Machtens, J.-P.; Briones, R.; Alleva, C.; de Groot, B. L.; Fahlke, C. Gating Charge Calculations by Computational Electrophysiology Simulations. Biophys. J. 2017, 112, 1396– 1405, DOI: 10.1016/j.bpj.2017.02.01670https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC2sXkt1Sqs7c%253D&md5=416dc846d784f5ce514486d29d376372Gating Charge Calculations by Computational Electrophysiology SimulationsMachtens, Jan-Philipp; Briones, Rodolfo; Alleva, Claudia; de Groot, Bert L.; Fahlke, ChristophBiophysical Journal (2017), 112 (7), 1396-1405CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)Elec. cell signaling requires adjustment of ion channel, receptor, or transporter function in response to changes in membrane potential. For the majority of such membrane proteins, the mol. details of voltage sensing remain insufficiently understood. Here, we present a mol. dynamics simulation-based method to det. the underlying charge movement across the membrane-the gating charge-by measuring elec. capacitor properties of membrane-embedded proteins. We illustrate the approach by calcg. the charge transfer upon membrane insertion of the HIV gp41 fusion peptide, and validate the method on two prototypical voltage-dependent proteins, the Kv1.2 K+ channel and the voltage sensor of the Ciona intestinalis voltage-sensitive phosphatase, against exptl. data. We then use the gating charge anal. to study how the T1 domain modifies voltage sensing in Kv1.2 channels and to investigate the voltage dependence of the initial binding of two Na+ ions in Na+-coupled glutamate transporters. Our simulation approach quantifies various mechanisms of voltage sensing, enables direct comparison with expts., and supports mechanistic interpretation of voltage sensitivity by fractional amino acid contributions.
- 71Groeneveld, M.; Slotboom, D.-J. Na+:Aspartate Coupling Stoichiometry in the Glutamate Transporter Homologue GltPh. Biochemistry 2010, 49, 3511– 3513, DOI: 10.1021/bi100430s71https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3cXksVantLY%253D&md5=45ca89490d62bfb11bfb2fde29c0db46Na+:Aspartate Coupling Stoichiometry in the Glutamate Transporter Homologue GltPhGroeneveld, Maarten; Slotboom, Dirk-JanBiochemistry (2010), 49 (17), 3511-3513CODEN: BICHAW; ISSN:0006-2960. (American Chemical Society)The Na+ aspartate symporter GltPh from Pyrococcus horikoshii is the only member of the glutamate transporter family for which crystal structures have been detd. The cation:aspartate coupling stoichiometry is unknown, thus hampering the elucidation of the ion coupling mechanism. Here we measure transport of 22Na+ and [14C]aspartate in proteoliposomes contg. purified GltPh and demonstrate that three Na+ ions are symported with aspartate.
- 72Samal, S. K.; Dash, M.; Van Vlierberghe, S.; Kaplan, D. L.; Chiellini, E.; van Blitterswijk, C.; Moroni, L.; Dubruel, P. Cationic polymers and their therapeutic potential. Chem. Soc. Rev. 2012, 41, 7147– 7194, DOI: 10.1039/c2cs35094g72https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38XhsVOnurfE&md5=418983b4d598ad31e8010d09cc1b60c5Cationic polymers and their therapeutic potentialSamal, Sangram Keshari; Dash, Mamoni; Van Vlierberghe, Sandra; Kaplan, David L.; Chiellini, Emo; van Blitterswijk, Clemens; Moroni, Lorenzo; Dubruel, PeterChemical Society Reviews (2012), 41 (21), 7147-7194CODEN: CSRVBR; ISSN:0306-0012. (Royal Society of Chemistry)A review. The last decade has witnessed enormous research focused on cationic polymers. Cationic polymers are the subject of intense research as non-viral gene delivery systems, due to their flexible properties, facile synthesis, robustness and proven gene delivery efficiency. Here, the authors review the most recent scientific advances in cationic polymers and their derivs. not only for gene delivery purposes but also for various alternative therapeutic applications. An overview of the synthesis and prepn. of cationic polymers is provided along with their inherent bioactive and intrinsic therapeutic potential. In addn., cationic polymer based biomedical materials are covered. Major progress in the fields of drug and gene delivery as well as tissue engineering applications is summarized in the present review.
- 73Kostritskii, A. Y.; Kondinskaia, D. A.; Nesterenko, A. M.; Gurtovenko, A. A. Adsorption of Synthetic Cationic Polymers on Model Phospholipid Membranes: Insight from Atomic-Scale Molecular Dynamics Simulations. Langmuir 2016, 32, 10402– 10414, DOI: 10.1021/acs.langmuir.6b0259373https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC28XhsFens7bP&md5=1d3d471cd08d1aaa4c10582314cbe459Adsorption of Synthetic Cationic Polymers on Model Phospholipid Membranes: Insight from Atomic-Scale Molecular Dynamics SimulationsKostritskii, Andrei Yu.; Kondinskaia, Diana A.; Nesterenko, Alexey M.; Gurtovenko, Andrey A.Langmuir (2016), 32 (40), 10402-10414CODEN: LANGD5; ISSN:0743-7463. (American Chemical Society)Although synthetic cationic polymers represent a promising class of effective antibacterial agents, the mol. mechanisms behind their antimicrobial activity remain poorly understood. To this end, we employ at.-scale mol. dynamics simulations to explore adsorption of several linear cationic polymers of different chem. structure and protonation (polyallylamine (PAA), polyethylenimine (PEI), polyvinylamine (PVA) and poly-L-lysine (PLL)) on model bacterial membranes (4:1 mixt. of zwitterionic phosphatidylethanolamine (PE) and anionic phosphatidylglycerol (PG) lipids). Overall, our findings show that binding of polycations to the anionic membrane surface effectively neutralizes its charge, leading to the reorientation of water mols. close to the lipid/water interface and to the partial release of counterions to the water phase. In certain cases one has even an overcharging of the membrane, which was shown to be a cooperative effect of polymer charges and lipid counterions. Protonated amine groups of polycations are found to interact preferably with head groups of anionic lipids, giving rise to formation of hydrogen bonds and to a noticeable lateral immobilization of the lipids. While all the above findings are mostly defined by the overall charge of a polymer, we found that the polymer architecture also matters. In particular, PVA and PEI are able to accumulate anionic PG lipids on the membrane surface, leading to lipid segregation. In turn, PLL whose charge twice exceeds charges of PVA/PEI does not induce such lipid segregation due to its considerably less compact architecture and relatively long side chains. We also show that partitioning of a polycation into the lipid/water interface is an interplay between its protonation level (the overall charge) and hydrophobicity of the backbone. Therefore, a possible strategy in creating highly efficient antimicrobial polymeric agents could be in tuning these polycation's properties through proper combination of protonated and hydrophobic blocks.
- 74Kondinskaia, D. A.; Gurtovenko, A. A. Supramolecular complexes of DNA with cationic polymers: The effect of polymer concentration. Polymer 2018, 142, 277– 284, DOI: 10.1016/j.polymer.2018.03.04874https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC1cXmslKlt7o%253D&md5=b803292ea83e9ebd5c210446b17035b6Supramolecular complexes of DNA with cationic polymers: The effect of polymer concentrationKondinskaia, Diana A.; Gurtovenko, Andrey A.Polymer (2018), 142 (), 277-284CODEN: POLMAG; ISSN:0032-3861. (Elsevier Ltd.)Supramol. complexes of DNA with cationic polymers are of tremendous importance as these polymers can successfully be used in gene therapy as delivery vectors of genetic material. In this work we employ atomistic mol. dynamics simulations to get an insight into the impact of polymer concn. on the structural and electrostatic properties of the complexes of DNA and cationic polymers. Four linear cationic polymers of different chem. structure and protonation state were studied: polyethylenimine (PEI), poly-L-lysine (PLL), polyvinylamine (PVA), and polyallylamine (PAA). For all considered polymers our computational findings clearly demonstrate that increasing the polymer concn. leads to the overcharging of the DNA mol. This is relevant as the overall pos. charge of the DNA-polycation complex is considered to be a prerequisite for efficient binding of the polyplexes to neg. charged cell membranes. However, the concn. effects themselves are found to be sensitive to the chem. structure of a polymer. In particular, flexible PEI chains, being characterized by the high affinity to DNA's phosphate groups, are able to bind to DNA in an independent manner. As a result, DNA and PEIs form compact complexes with the largest cumulative pos. charge among all four types of cationic polymers under study. In contrast, PVA chains show the affinity to the major groove of DNA due to the hydrophobic nature of their backbones. This leads to stronger interactions with the DNA mols. and to the competition between PVA chains for DNA binding sites. This competition is responsible for the obsd. satn. in the PVA-induced neutralizing of the DNA charges when the PVA concn. increases. As for PLL and PAA, both these polycations are found to be less effective in neutralizing the charge of the polyanionic DNA mols. as compared to PEI and PVA. PLL has relatively long side chains, so that some of them cannot access DNA phosphates and reside in aq. soln. In turn, PAA has a low protonation level, leading to a weak electrostatic binding to DNA. Overall, our findings can relate the properties of supramol. DNA-polycation complexes at elevated polymer concns. with the chem. structure of a polycation and therefore can be useful for the development of novel, highly efficient polycation-based delivery vectors of DNA.
- 75Rohs, R.; West, S. M.; Sosinsky, A.; Liu, P.; Mann, R. S.; Honig, B. The role of DNA shape in protein-DNA recognition. Nature 2009, 461, 1248– 1253, DOI: 10.1038/nature0847375https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BD1MXhtlOjtrvL&md5=2c314c36cc47209d7869e3fbda04f5e5The role of DNA shape in protein-DNA recognitionRohs, Remo; West, Sean M.; Sosinsky, Alona; Liu, Peng; Mann, Richard S.; Honig, BarryNature (London, United Kingdom) (2009), 461 (7268), 1248-1253CODEN: NATUAS; ISSN:0028-0836. (Nature Publishing Group)The recognition of specific DNA sequences by proteins is thought to depend on 2 types of mechanism: one that involves the formation of H-bonds with specific bases, primarily in the major groove, and one involving sequence-dependent deformations of the DNA helix. Here, by comprehensively analyzing the 3-dimensional structures of protein-DNA complexes, the authors show that the binding of Arg residues to narrow minor grooves is a widely used mode for protein-DNA recognition. This readout mechanism exploits the phenomenon that narrow minor grooves strongly enhance the neg. electrostatic potential of the DNA. The nucleosome core particle offers a prominent example of this effect. Minor-groove narrowing is often assocd. with the presence of A-tracts, AT-rich sequences that exclude the flexible TpA step. These findings indicate that the ability to detect local variations in DNA shape and electrostatic potential is a general mechanism that enables proteins to use information in the minor groove, which otherwise offers few opportunities for the formation of base-specific H-bonds, to achieve DNA-binding specificity.
- 76Andrews, A. J.; Luger, K. Nucleosome structure(s) and stability: variations on a theme. Annu. Rev. Biophys. 2011, 40, 99– 117, DOI: 10.1146/annurev-biophys-042910-15532976https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC3MXnvFaitbo%253D&md5=916cbdc89d17871e816f8e1b24926b90Nucleosome structure(s) and stability: variations on a themeAndrews, Andrew J.; Luger, KarolinAnnual Review of Biophysics (2011), 40 (), 99-117CODEN: ARBNCV; ISSN:1936-122X. (Annual Reviews Inc.)A review. Chromatin is a highly regulated, modular nucleoprotein complex that is central to many processes in eukaryotes. The organization of DNA into nucleosomes and higher-order structures has profound implications for DNA accessibility. Alternative structural states of the nucleosome, and the thermodn. parameters governing its assembly and disassembly, need to be considered in order to understand how access to nucleosomal DNA is regulated. Here, the authors provide a brief historical account of how the overriding perception regarding aspects of nucleosome structure has changed over the past 30 yr. The authors discuss recent tech. advances regarding nucleosome structure and its phys. characterization and review the evidence for alternative nucleosome conformations and their implications for nucleosome and chromatin dynamics.
- 77Hub, J. S.; de Groot, B. L. Detection of functional modes in protein dynamics. PLoS Comput. Biol. 2009, 5, e1000480, DOI: 10.1371/journal.pcbi.100048077https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BD1Mros1agug%253D%253D&md5=65120d75ba73b0d5fba27e45233f5ab2Detection of functional modes in protein dynamicsHub Jochen S; de Groot Bert LPLoS computational biology (2009), 5 (8), e1000480 ISSN:.Proteins frequently accomplish their biological function by collective atomic motions. Yet the identification of collective motions related to a specific protein function from, e.g., a molecular dynamics trajectory is often non-trivial. Here, we propose a novel technique termed "functional mode analysis" that aims to detect the collective motion that is directly related to a particular protein function. Based on an ensemble of structures, together with an arbitrary "functional quantity" that quantifies the functional state of the protein, the technique detects the collective motion that is maximally correlated to the functional quantity. The functional quantity could, e.g., correspond to a geometric, electrostatic, or chemical observable, or any other variable that is relevant to the function of the protein. In addition, the motion that displays the largest likelihood to induce a substantial change in the functional quantity is estimated from the given protein ensemble. Two different correlation measures are applied: first, the Pearson correlation coefficient that measures linear correlation only; and second, the mutual information that can assess any kind of interdependence. Detecting the maximally correlated motion allows one to derive a model for the functional state in terms of a single collective coordinate. The new approach is illustrated using a number of biomolecules, including a polyalanine-helix, T4 lysozyme, Trp-cage, and leucine-binding protein.
- 78Krivobokova, T.; Briones, R.; Hub, J. S.; Munk, A.; de Groot, B. L. Partial least-squares functional mode analysis: application to the membrane proteins AQP1, Aqy1, and CLC-ec1. Biophys. J. 2012, 103, 786– 796, DOI: 10.1016/j.bpj.2012.07.02278https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A528%3ADC%252BC38Xht1ektLbN&md5=b40f9fbb32b434d8a8b71f0956e4dd9bPartial Least-Squares Functional Mode Analysis: Application to the Membrane Proteins AQP1, Aqy1, and CLC-ec1Krivobokova, Tatyana; Briones, Rodolfo; Hub, Jochen S.; Munk, Axel; de Groot, Bert L.Biophysical Journal (2012), 103 (4), 786-796CODEN: BIOJAU; ISSN:0006-3495. (Cell Press)We introduce an approach based on the recently introduced functional mode anal. to identify collective modes of internal dynamics that maximally correlate to an external order parameter of functional interest. Input structural data can be either exptl. detd. structure ensembles or simulated ensembles, such as mol. dynamics trajectories. Partial least-squares regression is shown to yield a robust soln. to the multidimensional optimization problem, with a minimal and controllable risk of overfitting, as shown by extensive cross-validation. Several examples illustrate that the partial least-squares-based functional mode anal. successfully reveals the collective dynamics underlying the fluctuations in selected functional order parameters. Applications to T4 lysozyme, the Trp-cage, the aquaporin channels Aqy1 and hAQP1, and the CLC-ec1 chloride antiporter are presented in which the active site geometry, the hydrophobic solvent-accessible surface, channel gating dynamics, water permeability (pf), and a dihedral angle are defined as functional order parameters. The Aqy1 case reveals a gating mechanism that connects the inner channel gating residues with the protein surface, thereby providing an explanation of how the membrane may affect the channel. hAQP1 shows how the pf correlates with structural changes around the arom./arginine region of the pore. The CLC-ec1 application shows how local motions of the gating Glu148 couple to a collective motion that affects ion affinity in the pore.
- 79Bratanov, D.; Kovalev, K.; Machtens, J.-P.; Astashkin, R.; Chizhov, I.; Soloviov, D.; Volkov, D.; Polovinkin, V.; Zabelskii, D.; Mager, T.; Gushchin, I.; Rokitskaya, T.; Antonenko, Y.; Alekseev, A.; Shevchenko, V.; Yutin, N.; Rosselli, R.; Baeken, C.; Borshchevskiy, V.; Bourenkov, G.; Popov, A.; Balandin, T.; Büldt, G.; Manstein, D. J.; Rodriguez-Valera, F.; Fahlke, C.; Bamberg, E.; Koonin, E.; Gordeliy, V. Unique structure and function of viral rhodopsins. Nat. Commun. 2019, 10, 4939 DOI: 10.1038/s41467-019-12718-079https://chemport.cas.org/services/resolver?origin=ACS&resolution=options&coi=1%3ACAS%3A280%3ADC%252BB3MjjtFChsg%253D%253D&md5=f68e6b3fcfb1157bdb7cd373eaf7d8aeUnique structure and function of viral rhodopsinsBratanov Dmitry; Kovalev Kirill; Volkov Dmytro; Polovinkin Vitaly; Zabelskii Dmitrii; Alekseev Alexey; Shevchenko Vitaly; Baeken Christian; Borshchevskiy Valentin; Balandin Taras; Gordeliy Valentin; Bratanov Dmitry; Kovalev Kirill; Volkov Dmytro; Zabelskii Dmitrii; Alekseev Alexey; Shevchenko Vitaly; Baeken Christian; Borshchevskiy Valentin; Balandin Taras; Gordeliy Valentin; Bratanov Dmitry; Kovalev Kirill; Astashkin Roman; Polovinkin Vitaly; Gordeliy Valentin; Kovalev Kirill; Astashkin Roman; Soloviov Dmytro; Volkov Dmytro; Zabelskii Dmitrii; Gushchin Ivan; Alekseev Alexey; Shevchenko Vitaly; Borshchevskiy Valentin; Buldt Georg; Fahlke Christoph; Bamberg Ernst; Gordeliy Valentin; Kovalev Kirill; Alekseev Alexey; Shevchenko Vitaly; Kovalev Kirill; Machtens Jan-Philipp; Fahlke Christoph; Machtens Jan-Philipp; Chizhov Igor; Manstein Dietmar J; Soloviov Dmytro; Soloviov Dmytro; Soloviov Dmytro; Mager Thomas; Bamberg Ernst; Rokitskaya Tatyana; Antonenko Yuri; Yutin Natalya; Koonin Eugene; Rosselli Riccardo; Rodriguez-Valera Francisco; Bourenkov Gleb; Popov AlexanderNature communications (2019), 10 (1), 4939 ISSN:.Recently, two groups of rhodopsin genes were identified in large double-stranded DNA viruses. The structure and function of viral rhodopsins are unknown. We present functional characterization and high-resolution structure of an Organic Lake Phycodnavirus rhodopsin II (OLPVRII) of group 2. It forms a pentamer, with a symmetrical, bottle-like central channel with the narrow vestibule in the cytoplasmic part covered by a ring of 5 arginines, whereas 5 phenylalanines form a hydrophobic barrier in its exit. The proton donor E42 is placed in the helix B. The structure is unique among the known rhodopsins. Structural and functional data and molecular dynamics suggest that OLPVRII might be a light-gated pentameric ion channel analogous to pentameric ligand-gated ion channels, however, future patch clamp experiments should prove this directly. The data shed light on a fundamentally distinct branch of rhodopsins and may contribute to the understanding of virus-host interactions in ecologically important marine protists.
Supporting Information
Supporting Information
The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jctc.0c01246.
Dependence of the short-range part of the potential on β; convergence of P2X3 profiles calculated from the raw SPME potential; time convergence of the g_elpot calculations; raw data on the electrostatic potential and Na1 occupancy in the GltPh simulations; chemical structure of PVA and PAA; groove-specific division of space in DNA systems; and simulation times for DNA and GltPh systems (PDF)
Terms & Conditions
Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.