{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,3]],"date-time":"2026-10-03T22:08:13Z","timestamp":1791065293457,"version":"4.1.0"},"reference-count":36,"publisher":"MDPI AG","issue":"3","license":[{"start":{"date-parts":[[2021,3,7]],"date-time":"2021-03-07T00:00:00Z","timestamp":1615075200000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"content-domain":{"domain":[],"crossmark-restriction":false},"short-container-title":["Toxins"],"abstract":"<jats:p>Gating modifier toxins (GMTs) isolated from venomous organisms such as Protoxin-II (ProTx-II) and Huwentoxin-IV (HwTx-IV) that inhibit the voltage-gated sodium channel NaV1.7 by binding to its voltage-sensing domain II (VSDII) have been extensively investigated as non-opioid analgesics. However, reliably predicting how a mutation to a GMT will affect its potency for NaV1.7 has been challenging. Here, we hypothesize that structure-based computational methods can be used to predict such changes. We employ free-energy perturbation (FEP), a physics-based simulation method for predicting the relative binding free energy (RBFE) between molecules, and the cryo electron microscopy (cryo-EM) structures of ProTx-II and HwTx-IV bound to VSDII of NaV1.7 to re-predict the relative potencies of forty-seven point mutants of these GMTs for NaV1.7. First, FEP predicted these relative potencies with an overall root mean square error (RMSE) of 1.0 \u00b1 0.1 kcal\/mol and an R2 value of 0.66, equivalent to experimental uncertainty and an improvement over the widely used molecular-mechanics\/generalized born-surface area (MM-GB\/SA) RBFE method that had an RMSE of 3.9 \u00b1 0.8 kcal\/mol. Second, inclusion of an explicit membrane model was needed for the GMTs to maintain stable binding poses during the FEP simulations. Third, MM-GB\/SA and FEP were used to identify fifteen non-standard tryptophan mutants at ProTx-II[W24] predicted in silico to have a at least a 1 kcal\/mol gain in potency. These predicted potency gains are likely due to the displacement of high-energy waters as identified by the WaterMap algorithm for calculating the positions and thermodynamic properties of water molecules in protein binding sites. Our results expand the domain of applicability of FEP and set the stage for its prospective use in biologics drug discovery programs involving GMTs and NaV1.7.<\/jats:p>","DOI":"10.3390\/toxins13030193","type":"journal-article","created":{"date-parts":[[2021,3,7]],"date-time":"2021-03-07T21:52:15Z","timestamp":1615153935000},"page":"193","update-policy":"https:\/\/doi.org\/10.3390\/mdpi_crossmark_policy","source":"Crossref","is-referenced-by-count":19,"title":["Potency-Enhancing Mutations of Gating Modifier Toxins for the Voltage-Gated Sodium Channel NaV1.7 Can Be Predicted Using Accurate Free-Energy Calculations"],"prefix":"10.3390","volume":"13","author":[{"given":"Dana","family":"Katz","sequence":"first","affiliation":[{"name":"Schr\u04e7dinger, Inc., 120 West 45th St., New York, NY 10036, USA"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Dan","family":"Sindhikara","sequence":"additional","affiliation":[{"name":"Schr\u04e7dinger, Inc., 120 West 45th St., New York, NY 10036, USA"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Michael","family":"DiMattia","sequence":"additional","affiliation":[{"name":"Schr\u04e7dinger, Inc., 120 West 45th St., New York, NY 10036, USA"}],"role":[{"vocabulary":"crossref","role":"author"}]},{"given":"Abba E.","family":"Leffler","sequence":"additional","affiliation":[{"name":"Schr\u04e7dinger, Inc., 120 West 45th St., New York, NY 10036, USA"}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"1968","published-online":{"date-parts":[[2021,3,7]]},"reference":[{"key":"ref_1","doi-asserted-by":"crossref","first-page":"E3051","DOI":"10.1073\/pnas.1700761114","article-title":"Structures of closed and open states of a voltage-gated sodium channel","volume":"114","author":"Lenaeus","year":"2017","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"ref_2","doi-asserted-by":"crossref","first-page":"7093","DOI":"10.1021\/jm501981g","article-title":"Voltage-Gated Sodium Channels: Structure, Function, Pharmacology, and Clinical Indications","volume":"58","author":"Kraus","year":"2015","journal-title":"J. Med. Chem."},{"key":"ref_3","doi-asserted-by":"crossref","unstructured":"Yamakawa, K. (2016). Mutations of Voltage-Gated Sodium Channel Genes SCN1A and SCN2A in Epilepsy, Intellectual Disability, and Autism. Neuronal and Synaptic Dysfunction in Autism Spectrum Disorder and Intellectual Disability, Academic Press.","DOI":"10.1016\/B978-0-12-800109-7.00015-7"},{"key":"ref_4","doi-asserted-by":"crossref","first-page":"49","DOI":"10.1038\/nrn3404","article-title":"The NaV1.7 sodium channel: From molecule to man","volume":"14","author":"Yang","year":"2013","journal-title":"Nat. Rev. Neurosci."},{"key":"ref_5","doi-asserted-by":"crossref","first-page":"2299","DOI":"10.1021\/jm501765v","article-title":"Engineering potent and selective analogues of GpTx-1, a tarantula venom peptide antagonist of the Na(V)1.7 sodium channel","volume":"58","author":"Murray","year":"2015","journal-title":"J. Med. Chem."},{"key":"ref_6","doi-asserted-by":"crossref","first-page":"70","DOI":"10.1016\/j.cbpa.2017.03.004","article-title":"Progress and challenges in the optimization of toxin peptides for development as pain therapeutics","volume":"38","author":"Netirojjanakul","year":"2017","journal-title":"Curr. Opin. Chem. Biol."},{"key":"ref_7","doi-asserted-by":"crossref","first-page":"175","DOI":"10.1016\/j.tips.2009.12.007","article-title":"Targeting voltage sensors in sodium channels with spider toxins","volume":"31","author":"Bosmans","year":"2010","journal-title":"Trends Pharm. Sci."},{"key":"ref_8","doi-asserted-by":"crossref","first-page":"2704","DOI":"10.1021\/acs.jmedchem.5b01947","article-title":"Single Residue Substitutions That Confer Voltage-Gated Sodium Ion Channel Subtype Selectivity in the NaV1.7 Inhibitory Peptide GpTx-1","volume":"59","author":"Murray","year":"2016","journal-title":"J. Med. Chem."},{"key":"ref_9","doi-asserted-by":"crossref","first-page":"637","DOI":"10.1085\/jgp.116.5.637","article-title":"A hot spot for the interaction of gating modifier toxins with voltage-dependent ion channels","volume":"116","author":"Winterfield","year":"2000","journal-title":"J. Gen. Physiol."},{"key":"ref_10","doi-asserted-by":"crossref","first-page":"22707","DOI":"10.1074\/jbc.M113.461392","article-title":"Analysis of the structural and molecular basis of voltage-sensitive sodium channel inhibition by the spider toxin huwentoxin-IV (\u03bc-TRTX-Hh2a)","volume":"288","author":"Minassian","year":"2013","journal-title":"J. Biol. Chem."},{"key":"ref_11","doi-asserted-by":"crossref","first-page":"702","DOI":"10.1016\/j.cell.2018.12.018","article-title":"Structural Basis of Nav1.7 Inhibition by a Gating-Modifier Spider Toxin","volume":"176","author":"Xu","year":"2019","journal-title":"Cell"},{"key":"ref_12","doi-asserted-by":"crossref","first-page":"14734","DOI":"10.1021\/bi026546a","article-title":"Two tarantula peptides inhibit activation of multiple sodium channels","volume":"41","author":"Middleton","year":"2002","journal-title":"Biochemistry"},{"key":"ref_13","doi-asserted-by":"crossref","first-page":"27300","DOI":"10.1074\/jbc.M708447200","article-title":"Tarantula huwentoxin-IV inhibits neuronal sodium channels by binding to receptor site 4 and trapping the domain ii voltage sensor in the closed configuration","volume":"283","author":"Xiao","year":"2008","journal-title":"J. Biol. Chem."},{"key":"ref_14","doi-asserted-by":"crossref","first-page":"40","DOI":"10.1016\/j.peptides.2013.03.011","article-title":"Potency optimization of Huwentoxin-IV on hNav1.7: A neurotoxin TTX-S sodium-channel antagonist from the venom of the Chinese bird-eating spider Selenocosmia huwena","volume":"44","author":"Revell","year":"2013","journal-title":"Peptides"},{"key":"ref_15","doi-asserted-by":"crossref","first-page":"14187","DOI":"10.1073\/pnas.1922903117","article-title":"Employing NaChBac for cryo-EM analysis of toxin action on voltage-gated Na(+) channels in nanodisc","volume":"117","author":"Gao","year":"2020","journal-title":"Proc. Natl. Acad. Sci. USA"},{"key":"ref_16","doi-asserted-by":"crossref","first-page":"9041","DOI":"10.1074\/jbc.RA118.002553","article-title":"Gating modifier toxins isolated from spider venom: Modulation of voltage-gated sodium channels and the role of lipid membranes","volume":"293","author":"Agwa","year":"2018","journal-title":"J. Biol. Chem."},{"key":"ref_17","doi-asserted-by":"crossref","first-page":"6015","DOI":"10.1021\/bi0477034","article-title":"Solution Structure and Lipid Membrane Partitioning of VSTx1, an Inhibitor of the KvAP Potassium Channel","volume":"44","author":"Jung","year":"2005","journal-title":"Biochemistry"},{"key":"ref_18","doi-asserted-by":"crossref","first-page":"34333","DOI":"10.1038\/srep34333","article-title":"Molecular basis of the interaction between gating modifier spider toxins and the voltage sensor of voltage-gated ion channels","volume":"6","author":"Lau","year":"2016","journal-title":"Sci. Rep."},{"key":"ref_19","doi-asserted-by":"crossref","first-page":"17049","DOI":"10.1074\/jbc.M116.729095","article-title":"Interaction of Tarantula Venom Peptide ProTx-II with Lipid Membranes Is a Prerequisite for Its Inhibition of Human Voltage-gated Sodium Channel NaV1.7","volume":"291","author":"Henriques","year":"2016","journal-title":"J. Biol. Chem."},{"key":"ref_20","doi-asserted-by":"crossref","unstructured":"Shah, B., Sindhikara, D., Borrelli, K., and Leffler, A.E. (2020). Water Thermodynamics of Peptide Toxin Binding Sites on Ion Channels. Toxins, 12.","DOI":"10.3390\/toxins12100652"},{"key":"ref_21","doi-asserted-by":"crossref","first-page":"930","DOI":"10.1016\/j.jmb.2016.11.021","article-title":"Free Energy Perturbation Calculation of Relative Binding Free Energy between Broadly Neutralizing Antibodies and the gp120 Glycoprotein of HIV-1","volume":"429","author":"Clark","year":"2017","journal-title":"J. Mol. Biol."},{"key":"ref_22","doi-asserted-by":"crossref","first-page":"201","DOI":"10.1007\/978-1-4939-9608-7_9","article-title":"Protein-Ligand Binding Free Energy Calculations with FEP","volume":"2022","author":"Wang","year":"2019","journal-title":"Methods Mol. Biol."},{"key":"ref_23","doi-asserted-by":"crossref","first-page":"6061","DOI":"10.1021\/acs.jctc.0c00660","article-title":"Enhancing Water Sampling in Free Energy Calculations with Grand Canonical Monte Carlo","volume":"16","author":"Ross","year":"2020","journal-title":"J. Chem. Theory Comput."},{"key":"ref_24","doi-asserted-by":"crossref","unstructured":"Beard, H., Cholleti, A., Pearlman, D., Sherman, W., and Loving, K.A. (2013). Applying Physics-Based Scoring to Calculate Free Energies of Binding for Single Amino Acid Mutations in Protein-Protein Complexes. PLoS ONE, 8.","DOI":"10.1371\/journal.pone.0082849"},{"key":"ref_25","doi-asserted-by":"crossref","first-page":"1481","DOI":"10.1016\/j.jmb.2019.02.003","article-title":"Relative Binding Affinity Prediction of Charge-Changing Sequence Mutations with FEP in Protein-Protein Interfaces","volume":"431","author":"Clark","year":"2019","journal-title":"J. Mol. Biol."},{"key":"ref_26","doi-asserted-by":"crossref","unstructured":"Deplazes, E. (2017). Molecular Simulations of Disulfide-Rich Venom Peptides with Ion Channels and Membranes. Molecules, 22.","DOI":"10.3390\/molecules22030362"},{"key":"ref_27","doi-asserted-by":"crossref","first-page":"2695","DOI":"10.1021\/ja512751q","article-title":"Accurate and Reliable Prediction of Relative Ligand Binding Potency in Prospective Drug Discovery by Way of a Modern Free-Energy Calculation Protocol and Force Field","volume":"137","author":"Wang","year":"2015","journal-title":"J. Am. Chem. Soc."},{"key":"ref_28","doi-asserted-by":"crossref","first-page":"1550023","DOI":"10.1142\/S0219633615500236","article-title":"Calculation of free energy changes due to mutations from alchemical free energy simulations","volume":"14","author":"Rashid","year":"2015","journal-title":"J. Theor. Comput. Chem."},{"key":"ref_29","doi-asserted-by":"crossref","first-page":"2804","DOI":"10.1021\/acschemneuro.0c00293","article-title":"Rational Design of \u03b1-Conotoxin RegIIA Analogues Selectively Inhibiting the Human \u03b13\u03b22 Nicotinic Acetylcholine Receptor through Computational Scanning","volume":"11","author":"Xu","year":"2020","journal-title":"ACS Chem. Neurosci."},{"key":"ref_30","doi-asserted-by":"crossref","first-page":"420","DOI":"10.1016\/j.drudis.2009.01.012","article-title":"Healthy skepticism: Assessing realistic model performance","volume":"14","author":"Brown","year":"2009","journal-title":"Drug Discov. Today"},{"key":"ref_31","doi-asserted-by":"crossref","first-page":"38","DOI":"10.1016\/j.molcel.2020.10.039","article-title":"Structural Basis for High-Affinity Trapping of the NaV1.7 Channel in Its Resting State by Tarantula Toxin","volume":"81","author":"Wisedchaisri","year":"2021","journal-title":"Mol. Cell"},{"key":"ref_32","doi-asserted-by":"crossref","unstructured":"Negron, C., Pearlman, D.A., and Del Angel, G. (2019). Predicting mutations deleterious to function in beta-lactamase TEM1 using MM-GBSA. PLoS ONE, 14.","DOI":"10.1371\/journal.pone.0214015"},{"key":"ref_33","doi-asserted-by":"crossref","unstructured":"Gao, S., Na, R., Yang, L., Yu, H., Zhao, X., and Huang, X. (2020). Investigation of binding modes of spider toxin-human voltage-gated sodium channel subtybe 1.7. J. Biomol. Struct. Dyn., 1\u20139.","DOI":"10.1080\/07391102.2020.1783363"},{"key":"ref_34","doi-asserted-by":"crossref","first-page":"1315","DOI":"10.1016\/S0021-9258(17)49888-7","article-title":"Comprehensive engineering of the tarantula venom peptide huwentoxin-IV to inhibit the human voltage-gated sodium channel hNa(v)1.7","volume":"295","author":"Neff","year":"2020","journal-title":"J. Biol. Chem."},{"key":"ref_35","doi-asserted-by":"crossref","first-page":"1863","DOI":"10.1021\/acs.jctc.8b01026","article-title":"OPLS3e: Extending Force Field Coverage for Drug-Like Small Molecules","volume":"15","author":"Roos","year":"2019","journal-title":"J. Chem. Theory Comput."},{"key":"ref_36","doi-asserted-by":"crossref","first-page":"18378","DOI":"10.33011\/livecoms.2.1.18378","article-title":"Best Practices for Alchemical Free Energy Calculations","volume":"2","author":"Mey","year":"2020","journal-title":"Living J. Comput. Mol. Sci."}],"container-title":["Toxins"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/www.mdpi.com\/2072-6651\/13\/3\/193\/pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2025,10,11]],"date-time":"2025-10-11T05:34:28Z","timestamp":1760160868000},"score":1,"resource":{"primary":{"URL":"https:\/\/www.mdpi.com\/2072-6651\/13\/3\/193"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2021,3,7]]},"references-count":36,"journal-issue":{"issue":"3","published-online":{"date-parts":[[2021,3]]}},"alternative-id":["toxins13030193"],"URL":"https:\/\/doi.org\/10.3390\/toxins13030193","relation":{},"ISSN":["2072-6651"],"issn-type":[{"value":"2072-6651","type":"electronic"}],"subject":[],"published":{"date-parts":[[2021,3,7]]}}}