scholarly journals Potency-Enhancing Mutations of Gating Modifier Toxins for the Voltage-Gated Sodium Channel NaV1.7 Can Be Predicted Using Accurate Free-Energy Calculations

Toxins ◽  
2021 ◽  
Vol 13 (3) ◽  
pp. 193
Author(s):  
Dana Katz ◽  
Dan Sindhikara ◽  
Michael DiMattia ◽  
Abba E. Leffler

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 ± 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 ± 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.

2020 ◽  
Vol 295 (15) ◽  
pp. 5067-5080 ◽  
Author(s):  
Akello J. Agwa ◽  
Poanna Tran ◽  
Alexander Mueller ◽  
Hue N. T. Tran ◽  
Jennifer R. Deuis ◽  
...  

Huwentoxin-IV (HwTx-IV) is a gating modifier peptide toxin from spiders that has weak affinity for the lipid bilayer. As some gating modifier toxins have affinity for model lipid bilayers, a tripartite relationship among gating modifier toxins, voltage-gated ion channels, and the lipid membrane surrounding the channels has been proposed. We previously designed an HwTx-IV analogue (gHwTx-IV) with reduced negative charge and increased hydrophobic surface profile, which displays increased lipid bilayer affinity and in vitro activity at the voltage-gated sodium channel subtype 1.7 (NaV1.7), a channel targeted in pain management. Here, we show that replacements of the positively-charged residues that contribute to the activity of the peptide can improve gHwTx-IV's potency and selectivity for NaV1.7. Using HwTx-IV, gHwTx-IV, [R26A]gHwTx-IV, [K27A]gHwTx-IV, and [R29A]gHwTx-IV variants, we examined their potency and selectivity at human NaV1.7 and their affinity for the lipid bilayer. [R26A]gHwTx-IV consistently displayed the most improved potency and selectivity for NaV1.7, examined alongside off-target NaVs, compared with HwTx-IV and gHwTx-IV. The lipid affinity of each of the three novel analogues was weaker than that of gHwTx-IV, but stronger than that of HwTx-IV, suggesting a possible relationship between in vitro potency at NaV1.7 and affinity for lipid bilayers. In a murine NaV1.7 engagement model, [R26A]gHwTx-IV exhibited an efficacy comparable with that of native HwTx-IV. In summary, this study reports the development of an HwTx-IV analogue with improved in vitro selectivity for the pain target NaV1.7 and with an in vivo efficacy similar to that of native HwTx-IV.


Science ◽  
2019 ◽  
Vol 363 (6433) ◽  
pp. 1303-1308 ◽  
Author(s):  
Huaizong Shen ◽  
Dongliang Liu ◽  
Kun Wu ◽  
Jianlin Lei ◽  
Nieng Yan

Voltage-gated sodium channel Nav1.7 represents a promising target for pain relief. Here we report the cryo–electron microscopy structures of the human Nav1.7-β1-β2 complex bound to two combinations of pore blockers and gating modifier toxins (GMTs), tetrodotoxin with protoxin-II and saxitoxin with huwentoxin-IV, both determined at overall resolutions of 3.2 angstroms. The two structures are nearly identical except for minor shifts of voltage-sensing domain II (VSDII), whose S3-S4 linker accommodates the two GMTs in a similar manner. One additional protoxin-II sits on top of the S3-S4 linker in VSDIV. The structures may represent an inactivated state with all four VSDs “up” and the intracellular gate closed. The structures illuminate the path toward mechanistic understanding of the function and disease of Nav1.7 and establish the foundation for structure-aided development of analgesics.


2016 ◽  
Vol 47 (S 01) ◽  
Author(s):  
W. Fazeli ◽  
B. Schattling ◽  
B. Engeland ◽  
M. Friese ◽  
D. Isbrand

Author(s):  
Christina Schindler ◽  
Hannah Baumann ◽  
Andreas Blum ◽  
Dietrich Böse ◽  
Hans-Peter Buchstaller ◽  
...  

Here we present an evaluation of the binding affinity prediction accuracy of the free energy calculation method FEP+ on internal active drug discovery projects and on a large new public benchmark set.<br>


2020 ◽  
Author(s):  
E. Prabhu Raman ◽  
Thomas J. Paul ◽  
Ryan L. Hayes ◽  
Charles L. Brooks III

<p>Accurate predictions of changes to protein-ligand binding affinity in response to chemical modifications are of utility in small molecule lead optimization. Relative free energy perturbation (FEP) approaches are one of the most widely utilized for this goal, but involve significant computational cost, thus limiting their application to small sets of compounds. Lambda dynamics, also rigorously based on the principles of statistical mechanics, provides a more efficient alternative. In this paper, we describe the development of a workflow to setup, execute, and analyze Multi-Site Lambda Dynamics (MSLD) calculations run on GPUs with CHARMm implemented in BIOVIA Discovery Studio and Pipeline Pilot. The workflow establishes a framework for setting up simulation systems for exploratory screening of modifications to a lead compound, enabling the calculation of relative binding affinities of combinatorial libraries. To validate the workflow, a diverse dataset of congeneric ligands for seven proteins with experimental binding affinity data is examined. A protocol to automatically tailor fit biasing potentials iteratively to flatten the free energy landscape of any MSLD system is developed that enhances sampling and allows for efficient estimation of free energy differences. The protocol is first validated on a large number of ligand subsets that model diverse substituents, which shows accurate and reliable performance. The scalability of the workflow is also tested to screen more than a hundred ligands modeled in a single system, which also resulted in accurate predictions. With a cumulative sampling time of 150ns or less, the method results in average unsigned errors of under 1 kcal/mol in most cases for both small and large combinatorial libraries. For the multi-site systems examined, the method is estimated to be more than an order of magnitude more efficient than contemporary FEP applications. The results thus demonstrate the utility of the presented MSLD workflow to efficiently screen combinatorial libraries and explore chemical space around a lead compound, and thus are of utility in lead optimization.</p>


Sign in / Sign up

Export Citation Format

Share Document