Structural modeling of ion channels using AlphaFold2, RoseTTAFold2, and ESMFold

Ion channels play key roles in human physiology and are important targets in drug discovery. The atomic-scale structures of ion channels provide invaluable insights into a fundamental understanding of the molecular mechanisms of channel gating and modulation. Recent breakthroughs in deep learning-based computational methods, such as AlphaFold, RoseTTAFold, and ESMFold have transformed research in protein structure prediction and design. We review the application of AlphaFold, RoseTTAFold, and ESMFold to structural modeling of ion channels using representative voltage-gated ion channels, including human voltage-gated sodium (NaV) channel - NaV1.8, human voltage-gated calcium (CaV) channel - CaV1.1, and human voltage-gated potassium (KV) channel - KV1.3. We compared AlphaFold, RoseTTAFold, and ESMFold structural models of NaV1.8, CaV1.1, and KV1.3 with corresponding cryo-EM structures to assess details of their similarities and differences. Our findings shed light on the strengths and limitations of the current state-of-the-art deep learning-based computational methods for modeling ion channel structures, offering valuable insights to guide their future applications for ion channel research.

H3K4me1 recruits DNA repair proteins in plants

DNA repair proteins can be recruited by their histone reader domains to specific epigenomic features, with consequences on intragenomic mutation rate variation. Here, we investigated H3K4me1-associated hypomutation in plants. We first examined two proteins which, in plants, contain Tudor histone reader domains: PRECOCIOUS DISSOCIATION OF SISTERS 5 (PDS5C), involved in homology-directed repair, and MUTS HOMOLOG 6 (MSH6), a mismatch repair protein. The MSH6 Tudor domain of Arabidopsis (Arabidopsis thaliana) binds to H3K4me1 as previously demonstrated for PDS5C, which localizes to H3K4me1-rich gene bodies and essential genes. Mutations revealed by ultradeep sequencing of wild-type and msh6 knockout lines in Arabidopsis show that functional MSH6 is critical for the reduced rate of single base substitution mutations in gene bodies and H3K4me1-rich regions. We explored the breadth of these mechanisms among plants by examining a large rice (Oryza sativa) mutation dataset. H3K4me1-associated hypomutation is conserved in rice as are the H3K4me1 binding residues of MSH6 and PDS5C Tudor domains. Recruitment of DNA repair proteins by H3K4me1 in plants reveals convergent, but distinct, epigenome-recruited DNA repair mechanisms from those well described in humans. The emergent model of H3K4me1-recruited repair in plants is consistent with evolutionary theory regarding mutation modifier systems and offers mechanistic insight into intragenomic mutation rate variation in plants.

Elucidating molecular mechanisms of protoxin-II state-specific binding to the human NaV1.7 channel

Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells, and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-II (PTx2) has high selectivity for hNaV1.7 and is a valuable scaffold for designing novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD and the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.

Toward Digital Twin Technology for Precision Pharmacology

The authors demonstrate the feasibility of technological innovation for personalized medicine in the context of drug-induced arrhythmia. The authors use atomistic-scale structural models to predict rates of drug interaction with ion channels and make predictions of their effects in digital twins of induced pluripotent stem cell-derived cardiac myocytes. The authors construct a simplified multilayer, 1-dimensional ring model with sufficient path length to enable the prediction of arrhythmogenic dispersion of repolarization. Finally, the authors validate the computational pipeline prediction of drug effects with data and quantify drug-induced propensity to repolarization abnormalities in cardiac tissue. The technology is high throughput, computationally efficient, and low cost toward personalized pharmacologic prediction.

Harnessing AlphaFold to reveal state secrets: Prediction of hERG closed and inactivated states

To design safe, selective, and effective new therapies, there must be a deep understanding of the structure and function of the drug target. One of the most difficult problems to solve has been resolution of discrete conformational states of transmembrane ion channel proteins. An example is KV11.1 (hERG), comprising the primary cardiac repolarizing current, IKr. hERG is a notorious drug anti-target against which all promising drugs are screened to determine potential for arrhythmia. Drug interactions with the hERG inactivated state are linked to elevated arrhythmia risk, and drugs may become trapped during channel closure. However, the structural details of multiple conformational states have remained elusive. Here, we guided AlphaFold2 to predict plausible hERG inactivated and closed conformations, obtaining results consistent with myriad available experimental data. Drug docking simulations demonstrated hERG state-specific drug interactions aligning well with experimental results, revealing that most drugs bind more effectively in the inactivated state and are trapped in the closed state. Molecular dynamics simulations demonstrated ion conduction that aligned with earlier studies. Finally, we identified key molecular determinants of state transitions by analyzing interaction networks across closed, open, and inactivated states in agreement with earlier mutagenesis studies. Here, we demonstrate a readily generalizable application of AlphaFold2 as a novel method to predict discrete protein conformations and novel linkages from structure to function.

Structural modeling of hERG channel-drug interactions using Rosetta

The human ether-a-go-go-related gene (hERG) not only encodes a potassium-selective voltage-gated ion channel essential for normal electrical activity in the heart but is also a major drug anti-target. Genetic hERG mutations and blockage of the channel pore by drugs can cause long QT syndrome, which predisposes individuals to potentially deadly arrhythmias. However, not all hERG-blocking drugs are proarrhythmic, and their differential affinities to discrete channel conformational states have been suggested to contribute to arrhythmogenicity. We used Rosetta electron density refinement and homology modeling to build structural models of open-state hERG channel wild-type and mutant variants (Y652A, F656A, and Y652A/F656 A) and a closed-state wild-type channel based on cryo-electron microscopy structures of hERG and EAG1 channels. These models were used as protein targets for molecular docking of charged and neutral forms of amiodarone, nifekalant, dofetilide, d/l-sotalol, flecainide, and moxifloxacin. We selected these drugs based on their different arrhythmogenic potentials and abilities to facilitate hERG current. Our docking studies and clustering provided atomistic structural insights into state-dependent drug-channel interactions that play a key role in differentiating safe and harmful hERG blockers and can explain hERG channel facilitation through drug interactions with its open-state hydrophobic pockets.

A Novel Peptide Prevents Enterotoxin- and Inflammation-Induced Intestinal Fluid Secretion by Stimulating Sodium-Hydrogen Exchanger 3 Activity

BACKGROUND & AIMS

Acute diarrheal diseases are the second most common cause of infant mortality in developing countries. This is contributed to by lack of effective drug therapy that shortens the duration or lessens the volume of diarrhea. The epithelial brush border sodium (Na+)/hydrogen (H+) exchanger 3 (NHE3) accounts for a major component of intestinal Na+ absorption and is inhibited in most diarrheas. Because increased intestinal Na+ absorption can rehydrate patients with diarrhea, NHE3 has been suggested as a potential druggable target for drug therapy for diarrhea.

METHODS

A peptide (sodium-hydrogen exchanger 3 stimulatory peptide [N3SP]) was synthesized to mimic the part of the NHE3 C-terminus that forms a multiprotein complex that inhibits NHE3 activity. The effect of N3SP on NHE3 activity was evaluated in NHE3-transfected fibroblasts null for other plasma membrane NHEs, a human colon cancer cell line that models intestinal absorptive enterocytes (Caco-2/BBe), human enteroids, and mouse intestine in vitro and in vivo. N3SP was delivered into cells via a hydrophobic fluorescent maleimide or nanoparticles.

RESULTS

N3SP uptake stimulated NHE3 activity at nmol/L concentrations under basal conditions and partially reversed the reduced NHE3 activity caused by elevated adenosine 3',5'-cyclic monophosphate, guanosine 3',5'-cyclic monophosphate, and Ca2+ in cell lines and in in vitro mouse intestine. N3SP also stimulated intestinal fluid absorption in the mouse small intestine in vivo and prevented cholera toxin-, Escherichia coli heat-stable enterotoxin-, and cluster of differentiation 3 inflammation-induced fluid secretion in a live mouse intestinal loop model.

CONCLUSIONS

These findings suggest pharmacologic stimulation of NHE3 activity as an efficacious approach for the treatment of moderate/severe diarrheal diseases.

Elucidating Molecular Mechanisms of Protoxin-2 State-specific Binding to the Human NaV1.7 Channel

Human voltage-gated sodium (hNaV) channels are responsible for initiating and propagating action potentials in excitable cells and mutations have been associated with numerous cardiac and neurological disorders. hNaV1.7 channels are expressed in peripheral neurons and are promising targets for pain therapy. The tarantula venom peptide protoxin-2 (PTx2) has high selectivity for hNaV1.7 and serves as a valuable scaffold to design novel therapeutics to treat pain. Here, we used computational modeling to study the molecular mechanisms of the state-dependent binding of PTx2 to hNaV1.7 voltage-sensing domains (VSDs). Using Rosetta structural modeling methods, we constructed atomistic models of the hNaV1.7 VSD II and IV in the activated and deactivated states with docked PTx2. We then performed microsecond-long all-atom molecular dynamics (MD) simulations of the systems in hydrated lipid bilayers. Our simulations revealed that PTx2 binds most favorably to the deactivated VSD II and activated VSD IV. These state-specific interactions are mediated primarily by PTx2's residues R22, K26, K27, K28, and W30 with VSD as well as the surrounding membrane lipids. Our work revealed important protein-protein and protein-lipid contacts that contribute to high-affinity state-dependent toxin interaction with the channel. The workflow presented will prove useful for designing novel peptides with improved selectivity and potency for more effective and safe treatment of pain.

Opening of capsaicin receptor TRPV1 is stabilized equally by its four subunits

Capsaicin receptor TRPV1 is a nociceptor for vanilloid molecules, such as capsaicin and resiniferatoxin (RTX). Even though cryo-EM structures of TRPV1 in complex with these molecules are available, how their binding energetically favors the open conformation is not known. Here, we report an approach to control the number of bound RTX molecules (0-4) in functional rat TRPV1. The approach allowed direct measurements of each of the intermediate open states under equilibrium conditions at both macroscopic and single-molecule levels. We found that RTX binding to each of the four subunits contributes virtually the same activation energy, which we estimated to be 1.70 to 1.86 kcal/mol and found to arise predominately from destabilizing the closed conformation. We further showed that sequential bindings of RTX increase open probability without altering single-channel conductance, confirming that there is likely a single open-pore conformation for TRPV1 activated by RTX.

Energetic landscape of polycystin channel gating

Members of the polycystin family (PKD2 and PKD2L1) of transient receptor potential (TRP) channels conduct Ca²⁺ and depolarizing monovalent cations. Variants in PKD2 cause autosomal dominant polycystic kidney disease (ADPKD) in humans, whereas loss of PKD2L1 expression causes seizure susceptibility in mice. Understanding structural and functional regulation of these channels will provide the basis for interpreting their molecular dysregulation in disease states. However, the complete structures of polycystins are unresolved, as are the conformational changes regulating their conductive states. To provide a holistic understanding of the polycystin gating cycle, we use computational prediction tools to model missing PKD2L1 structural motifs and evaluate more than 150 mutations in an unbiased mutagenic functional screen of the entire pore module. Our results provide an energetic landscape of the polycystin pore, which enumerates gating sensitive sites and interactions required for opening, inactivation, and subsequent desensitization. These findings identify the external pore helices and specific cross-domain interactions as critical structural regulators controlling the polycystin ion channel conductive and nonconductive states.

Molecular determinants of μ-conotoxin KIIIA interaction with the human voltage-gated sodium channel NaV1.7

The voltage-gated sodium (NaV) channel subtype NaV1.7 plays a critical role in pain signaling, making it an important drug target. Here we studied the molecular interactions between μ-Conotoxin KIIIA (KIIIA) and the human NaV1.7 channel (hNaV1.7). We developed a structural model of hNaV1.7 using Rosetta computational modeling and performed in silico docking of KIIIA using RosettaDock to predict residues forming specific pairwise contacts between KIIIA and hNaV1.7. We experimentally validated these contacts using mutant cycle analysis. Comparison between our KIIIA-hNaV1.7 model and the cryo-EM structure of KIIIA-hNaV1.2 revealed key similarities and differences between NaV channel subtypes with potential implications for the molecular mechanism of toxin block. The accuracy of our integrative approach, combining structural data with computational modeling, experimental validation, and molecular dynamics simulations, suggests that Rosetta structural predictions will be useful for rational design of novel biologics targeting specific NaV channels.

Elucidation of a dynamic interplay between a beta-2 adrenergic receptor, its agonist, and stimulatory G protein

G protein-coupled receptors (GPCRs) represent the largest group of membrane receptors for transmembrane signal transduction. Ligand-induced activation of GPCRs triggers G protein activation followed by various signaling cascades. Understanding the structural and energetic determinants of ligand binding to GPCRs and GPCRs to G proteins is crucial to the design of pharmacological treatments targeting specific conformations of these proteins to precisely control their signaling properties. In this study, we focused on interactions of a prototypical GPCR, beta-2 adrenergic receptor (β2AR), with its endogenous agonist, norepinephrine (NE), and the stimulatory G protein (Gs). Using molecular dynamics (MD) simulations, we demonstrated the stabilization of cationic NE, NE(+), binding to β2AR by Gs protein recruitment, in line with experimental observations. We also captured the partial dissociation of the ligand from β2AR and the conformational interconversions of Gs between closed and open conformations in the NE(+)-β2AR-Gs ternary complex while it is still bound to the receptor. The variation of NE(+) binding poses was found to alter Gs α subunit (Gsα) conformational transitions. Our simulations showed that the interdomain movement and the stacking of Gsα α1 and α5 helices are significant for increasing the distance between the Gsα and β2AR, which may indicate a partial dissociation of Gsα The distance increase commences when Gsα is predominantly in an open state and can be triggered by the intracellular loop 3 (ICL3) of β2AR interacting with Gsα, causing conformational changes of the α5 helix. Our results help explain molecular mechanisms of ligand and GPCR-mediated modulation of G protein activation.

Structural modeling of peptide toxin-ion channel interactions using RosettaDock

Voltage-gated ion channels play essential physiological roles in action potential generation and propagation. Peptidic toxins from animal venoms target ion channels and provide useful scaffolds for the rational design of novel channel modulators with enhanced potency and subtype selectivity. Despite recent progress in obtaining experimental structures of peptide toxin-ion channel complexes, structural determination of peptide toxins bound to ion channels in physiologically important states remains challenging. Here we describe an application of RosettaDock approach to the structural modeling of peptide toxins interactions with ion channels. We tested this approach on 10 structures of peptide toxin-ion channel complexes and demonstrated that it can sample near-native structures in all tested cases. Our approach will be useful for improving the understanding of the molecular mechanism of natural peptide toxin modulation of ion channel gating and for the structural modeling of novel peptide-based ion channel modulators.

TRPV1 Opening is Stabilized Equally by Its Four Subunits

Capsaicin receptor TRPV1 is a nociceptor for vanilloid molecules such as capsaicin and resiniferatoxin (RTX). Even though cryo-EM structures of TRPV1 in complex with these molecules are available, how their binding energetically favors the open conformation is not known. Here we report an approach to control the number of bound RTX molecules (0-to-4) in functional mouse TRPV1. The approach allowed direct measurements of each of the intermediate open states under equilibrium conditions at both macroscopic and single-molecule levels. We found that RTX binding to each of the four subunits contributes virtually the same activation energy, which we estimated to be 1.86 kcal/mol and found to arise predominately from destabilizing the closed conformation. We further showed that sequential bindings of RTX increase open probability without altering single-channel conductance, confirming that there is likely a single open-pore conformation for TRPV1 activated by RTX.

Adenylyl cyclase isoform 1 contributes to sinoatrial node automaticity via functional microdomains

Sinoatrial node (SAN) cells are the heart's primary pacemaker. Their activity is tightly regulated by β-adrenergic receptor (β-AR) signaling. Adenylyl cyclase (AC) is a key enzyme in the β-AR pathway that catalyzes the production of cAMP. There are current gaps in our knowledge regarding the dominant AC isoforms and the specific roles of Ca2+-activated ACs in the SAN. The current study tests the hypothesis that distinct AC isoforms are preferentially expressed in the SAN and compartmentalize within microdomains to orchestrate heart rate regulation during β-AR signaling. In contrast to atrial and ventricular myocytes, SAN cells express a diverse repertoire of ACs, with ACI as the predominant Ca2+-activated isoform. Although ACI-KO (ACI-/-) mice exhibit normal cardiac systolic or diastolic function, they experience SAN dysfunction. Similarly, SAN-specific CRISPR/Cas9-mediated gene silencing of ACI results in sinus node dysfunction. Mechanistically, hyperpolarization-activated cyclic nucleotide-gated 4 (HCN4) channels form functional microdomains almost exclusively with ACI, while ryanodine receptor and L-type Ca2+ channels likely compartmentalize with ACI and other AC isoforms. In contrast, there were no significant differences in T-type Ca2+ and Na+ currents at baseline or after β-AR stimulation between WT and ACI-/- SAN cells. Due to its central characteristic feature as a Ca2+-activated isoform, ACI plays a unique role in sustaining the rise of local cAMP and heart rates during β-AR stimulation. The findings provide insights into the critical roles of the Ca2+-activated isoform of AC in sustaining SAN automaticity that is distinct from contractile cardiomyocytes.

Structural modeling of the hERG potassium channel and associated drug interactions

The voltage-gated potassium channel, KV11.1, encoded by the human Ether-à-go-go-Related Gene (hERG), is expressed in cardiac myocytes, where it is crucial for the membrane repolarization of the action potential. Gating of the hERG channel is characterized by rapid, voltage-dependent, C-type inactivation, which blocks ion conduction and is suggested to involve constriction of the selectivity filter. Mutations S620T and S641A/T within the selectivity filter region of hERG have been shown to alter the voltage dependence of channel inactivation. Because hERG channel blockade is implicated in drug-induced arrhythmias associated with both the open and inactivated states, we used Rosetta to simulate the effects of hERG S620T and S641A/T mutations to elucidate conformational changes associated with hERG channel inactivation and differences in drug binding between the two states. Rosetta modeling of the S641A fast-inactivating mutation revealed a lateral shift of the F627 side chain in the selectivity filter into the central channel axis along the ion conduction pathway and the formation of four lateral fenestrations in the pore. Rosetta modeling of the non-inactivating mutations S620T and S641T suggested a potential molecular mechanism preventing F627 side chain from shifting into the ion conduction pathway during the proposed inactivation process. Furthermore, we used Rosetta docking to explore the binding mechanism of highly selective and potent hERG blockers - dofetilide, terfenadine, and E4031. Our structural modeling correlates well with much, but not all, existing experimental evidence involving interactions of hERG blockers with key residues in hERG pore and reveals potential molecular mechanisms of ligand interactions with hERG in an inactivated state.

Common Structural Pattern for Flecainide Binding in Atrial-Selective Kv1.5 and Nav1.5 Channels: A Computational Approach

Atrial fibrillation (AF) is the most common cardiac arrhythmia. Its treatment includes antiarrhythmic drugs (AADs) to modulate the function of cardiac ion channels. However, AADs have been limited by proarrhythmic effects, non-cardiovascular toxicities as well as often modest antiarrhythmic efficacy. Theoretical models showed that a combined blockade of Na_v1.5 (and its current, I_Na) and K_v1.5 (and its current, I_Kur) ion channels yield a synergistic anti-arrhythmic effect without alterations in ventricles. We focused on K_v1.5 and Na_v1.5 to search for structural similarities in their binding site (BS) for flecainide (a common blocker and widely prescribed AAD) as a first step for prospective rational multi-target directed ligand (MTDL) design strategies. We present a computational workflow for a flecainide BS comparison in a flecainide-K_v1.5 docking model and a solved structure of the flecainide-Na_v1.5 complex. The workflow includes docking, molecular dynamics, BS characterization and pattern matching. We identified a common structural pattern in flecainide BS for these channels. The latter belongs to the central cavity and consists of a hydrophobic patch and a polar region, involving residues from the S6 helix and P-loop. Since the rational MTDL design for AF is still incipient, our findings could advance multi-target atrial-selective strategies for AF treatment.

Structural and Functional Characterization of a Novel Scorpion Toxin that Inhibits NaV1.8 via Interactions With the DI Voltage Sensor and DII Pore Module

Voltage-gated sodium channel Na_V1.8 regulates transmission of pain signals to the brain. While Na_V1.8 has the potential to serve as a drug target, the molecular mechanisms that shape Na_V1.8 gating are not completely understood, particularly mechanisms that couple activation to inactivation. Interactions between toxin producing animals and their predators provide a novel approach for investigating Na_V structure-function relationships. Arizona bark scorpions produce Na⁺ channel toxins that initiate pain signaling. However, in predatory grasshopper mice, toxins inhibit Na_V1.8 currents and block pain signals. A screen of synthetic peptide toxins predicted from bark scorpion venom showed that peptide NaTx36 inhibited Na⁺ current recorded from a recombinant grasshopper mouse Na_V1.8 channel (OtNa_V1.8). Toxin NaTx36 hyperpolarized OtNa_V1.8 activation, steady-state fast inactivation, and slow inactivation. Mutagenesis revealed that the first gating charge in the domain I (DI) S4 voltage sensor and an acidic amino acid (E) in the DII SS2 – S6 pore loop are critical for the inhibitory effects of NaTx36. Computational modeling showed that a DI S1 – S2 asparagine (N) stabilizes the NaTx36 – OtNa_V1.8 complex while residues in the DI S3 – S4 linker and S4 voltage sensor form electrostatic interactions that allow a toxin glutamine (Q) to contact the first S4 gating charge. Surprisingly, the models predicted that NaTx36 contacts amino acids in the DII S5 – SS1 pore loop instead of the SS2 – S6 loop; the DII SS2 – S6 loop motif (QVSE) alters the conformation of the DII S5 – SS1 pore loop, enhancing allosteric interactions between toxin and the DII S5 – SS1 pore loop. Few toxins have been identified that modify Na_V1.8 gating. Moreover, few toxins have been described that modify sodium channel gating via the DI S4 voltage sensor. Thus, NaTx36 and OtNa_V1.8 provide tools for investigating the structure-activity relationship between channel activation and inactivation gating, and the connection to alternative pain phenotypes.

Structural insights into TRPV2 activation by small molecules

Transient receptor potential vanilloid 2 (TRPV2) is involved in many critical physiological and pathophysiological processes, making it a promising drug target. Here we present cryo-electron microscopy (cryo-EM) structures of rat TRPV2 in lipid nanodiscs activated by 2-aminoethoxydiphenyl borate (2-APB) and propose a TRPV2-specific 2-ABP binding site at the interface of S5 of one monomer and the S4-S5 linker of the adjacent monomer. In silico docking and electrophysiological studies confirm the key role of His521 and Arg539 in 2-APB activation of TRPV2. Additionally, electrophysiological experiments show that the combination of 2-APB and cannabidiol has a synergetic effect on TRPV2 activation, and cryo-EM structures demonstrate that both drugs were able to bind simultaneously. Together, our cryo-EM structures represent multiple functional states of the channel, providing a native picture of TRPV2 activation by small molecules and a structural framework for the development of TRPV2-specific activators.

Towards Structure-Guided Development of Pain Therapeutics Targeting Voltage-Gated Sodium Channels

Voltage-gated sodium (Na_V) channels are critical molecular determinants of action potential generation and propagation in excitable cells. Normal Na_V channel function disruption can affect physiological neuronal signaling and lead to increased sensitivity to pain, congenital indifference to pain, uncoordinated movement, seizures, or paralysis. Human genetic studies have identified human Na_V1.7 (hNa_V1.7), hNa_V1.8, and hNa_V1.9 channel subtypes as crucial players in pain signaling. The premise that subtype selective Na_V inhibitors can reduce pain has been reinforced through intensive target validation and therapeutic development efforts. However, an ideal therapeutic has yet to emerge. This review is focused on recent progress, current challenges, and future opportunities to develop Na_V channel targeting small molecules and peptides as non-addictive therapeutics to treat pain.

Computational design of peptides to target NaV1.7 channel with high potency and selectivity for the treatment of pain

The voltage-gated sodium Na_V1.7 channel plays a key role as a mediator of action potential propagation in C-fiber nociceptors and is an established molecular target for pain therapy. ProTx-II is a potent and moderately selective peptide toxin from tarantula venom that inhibits human Na_V1.7 activation. Here we used available structural and experimental data to guide Rosetta design of potent and selective ProTx-II-based peptide inhibitors of human Na_V1.7 channels. Functional testing of designed peptides using electrophysiology identified the PTx2-3127 and PTx2-3258 peptides with IC₅₀s of 7 nM and 4 nM for hNa_V1.7 and more than 1000-fold selectivity over human Na_V1.1, Na_V1.3, Na_V1.4, Na_V1.5, Na_V1.8, and Na_V1.9 channels. PTx2-3127 inhibits Na_V1.7 currents in mouse and human sensory neurons and shows efficacy in rat models of chronic and thermal pain when administered intrathecally. Rationally designed peptide inhibitors of human Na_V1.7 channels have transformative potential to define a new class of biologics to treat pain.

Single amino acid residue mediates reciprocal specificity in two mosquito odorant receptors

The southern house mosquito, Culex quinquefasciatus, utilizes two odorant receptors, CquiOR10 and CquiOR2, narrowly tuned to oviposition attractants and well conserved among mosquito species. They detect skatole and indole, respectively, with reciprocal specificity. We swapped the transmembrane (TM) domains of CquiOR10 and CquiOR2 and identified TM2 as a specificity determinant. With additional mutations, we showed that CquiOR10A73L behaved like CquiOR2. Conversely, CquiOR2L74A recapitulated CquiOR10 specificity. Next, we generated structural models of CquiOR10 and CquiOR10A73L using RoseTTAFold and AlphaFold and docked skatole and indole using RosettaLigand. These modeling studies suggested space-filling constraints around A73. Consistent with this hypothesis, CquiOR10 mutants with a bulkier residue (Ile, Val) were insensitive to skatole and indole, whereas CquiOR10A73G retained the specificity to skatole and showed a more robust response than the wildtype receptor CquiOR10. On the other hand, Leu to Gly mutation of the indole receptor CquiOR2 reverted the specificity to skatole. Lastly, CquiOR10A73L, CquiOR2, and CquiOR2L74I were insensitive to 3-ethylindole, whereas CquiOR2L74A and CquiOR2L74G gained activity. Additionally, CquiOR10A73G gave more robust responses to 3-ethylindole than CquiOR10. Thus, we suggest the specificity of these receptors is mediated by a single amino acid substitution, leading to finely tuned volumetric space to accommodate specific oviposition attractants.

Ensuring scientific reproducibility in bio-macromolecular modeling via extensive, automated benchmarks

Each year vast international resources are wasted on irreproducible research. The scientific community has been slow to adopt standard software engineering practices, despite the increases in high-dimensional data, complexities of workflows, and computational environments. Here we show how scientific software applications can be created in a reproducible manner when simple design goals for reproducibility are met. We describe the implementation of a test server framework and 40 scientific benchmarks, covering numerous applications in Rosetta bio-macromolecular modeling. High performance computing cluster integration allows these benchmarks to run continuously and automatically. Detailed protocol captures are useful for developers and users of Rosetta and other macromolecular modeling tools. The framework and design concepts presented here are valuable for developers and users of any type of scientific software and for the scientific community to create reproducible methods. Specific examples highlight the utility of this framework, and the comprehensive documentation illustrates the ease of adding new tests in a matter of hours.

Disruption of protein quality control of the human ether-à-go-go related gene K+ channel results in profound long QT syndrome

BACKGROUND

Long QT syndrome (LQTS) is a hereditary disease that predisposes patients to life-threatening cardiac arrhythmias and sudden cardiac death. Our previous study of the human ether-à-go-go related gene (hERG)-encoded K⁺ channel (K_v11.1) supports an association between hERG and RING finger protein 207 (RNF207) variants in aggravating the onset and severity of LQTS, specifically T613M hERG (hERG_T613M) and RNF207 frameshift (RNF207_G603fs) mutations. However, the underlying mechanistic underpinning remains unknown.

OBJECTIVE

The purpose of the present study was to test the role of RNF207 in the function of hERG-encoded K⁺ channel subunits.

METHODS

Whole-cell patch-clamp experiments were performed in human embryonic kidney (HEK 293) cells and human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) together with immunofluorescent confocal and high resolution microscopy, auto-ubiquitinylation assays, and co-immunoprecipitation experiments to test the functional interactions between hERG and RNF207.

RESULTS

Here, we demonstrated that RNF207 serves as an E3 ubiquitin ligase and targets misfolded hERG_T613M proteins for degradation. RNF207_G603fs exhibits decreased activity and hinders the normal degradation pathway; this increases the levels of hERG_T613M subunits and their dominant-negative effect on the wild-type subunits, ultimately resulting in decreased current density. Similar findings are shown for hERG_A614V, a known dominant-negative mutant subunit. Finally, the presence of RNF207_G603fs with hERG_T613M results in significantly prolonged action potential durations and reduced hERG current in human-induced pluripotent stem cell-derived cardiomyocytes.

CONCLUSION

Our study establishes RNF207 as an interacting protein serving as a ubiquitin ligase for hERG-encoded K⁺ channel subunits. Normal function of RNF207 is critical for the quality control of hERG subunits and consequently cardiac repolarization. Moreover, our study provides evidence for protein quality control as a new paradigm in life-threatening cardiac arrhythmias in patients with LQTS.

Different arrhythmia-associated calmodulin mutations have distinct effects on cardiac SK channel regulation

Calmodulin (CaM) plays a critical role in intracellular signaling and regulation of Ca²⁺-dependent proteins and ion channels. Mutations in CaM cause life-threatening cardiac arrhythmias. Among the known CaM targets, small-conductance Ca²⁺-activated K⁺ (SK) channels are unique, since they are gated solely by beat-to-beat changes in intracellular Ca²⁺. However, the molecular mechanisms of how CaM mutations may affect the function of SK channels remain incompletely understood. To address the structural and functional effects of these mutations, we introduced prototypical human CaM mutations in human induced pluripotent stem cell–derived cardiomyocyte-like cells (hiPSC-CMs). Using structural modeling and molecular dynamics simulation, we demonstrate that human calmodulinopathy-associated CaM mutations disrupt cardiac SK channel function via distinct mechanisms. CaM_D96V and CaM_D130G mutants reduce SK currents through a dominant-negative fashion. By contrast, specific mutations replacing phenylalanine with leucine result in conformational changes that affect helix packing in the C-lobe, which disengage the interactions between apo-CaM and the CaM-binding domain of SK channels. Distinct mutant CaMs may result in a significant reduction in the activation of the SK channels, leading to a decrease in the key Ca²⁺-dependent repolarization currents these channels mediate. The findings in this study may be generalizable to other interactions of mutant CaMs with Ca²⁺-dependent proteins within cardiac myocytes.

Molecular determinants of pro-arrhythmia proclivity of d- and l-sotalol via a multi-scale modeling pipeline

Drug isomers may differ in their proarrhythmia risk. An interesting example is the drug sotalol, an antiarrhythmic drug comprising d- and l- enantiomers that both block the hERG cardiac potassium channel and confer differing degrees of proarrhythmic risk. We developed a multi-scale in silico pipeline focusing on hERG channel – drug interactions and used it to probe and predict the mechanisms of pro-arrhythmia risks of the two enantiomers of sotalol. Molecular dynamics (MD) simulations predicted comparable hERG channel binding affinities for d- and l-sotalol, which were validated with electrophysiology experiments. MD derived thermodynamic and kinetic parameters were used to build multi-scale functional computational models of cardiac electrophysiology at the cell and tissue scales. Functional models were used to predict inactivated state binding affinities to recapitulate electrocardiogram (ECG) QT interval prolongation observed in clinical data. Our study demonstrates how modeling and simulation can be applied to predict drug effects from the atom to the rhythm for dl-sotalol and also increased proarrhythmia proclivity of d- vs. l-sotalol when accounting for stereospecific beta-adrenergic receptor blocking.

Toggle switch residues control allosteric transitions in bacterial adhesins by participating in a concerted repacking of the protein core

Critical molecular events that control conformational transitions in most allosteric proteins are ill-defined. The mannose-specific FimH protein of Escherichia coli is a prototypic bacterial adhesin that switches from an 'inactive' low-affinity state (LAS) to an 'active' high-affinity state (HAS) conformation allosterically upon mannose binding and mediates shear-dependent catch bond adhesion. Here we identify a novel type of antibody that acts as a kinetic trap and prevents the transition between conformations in both directions. Disruption of the allosteric transitions significantly slows FimH's ability to associate with mannose and blocks bacterial adhesion under dynamic conditions. FimH residues critical for antibody binding form a compact epitope that is located away from the mannose-binding pocket and is structurally conserved in both states. A larger antibody-FimH contact area is identified by NMR and contains residues Leu-34 and Val-35 that move between core-buried and surface-exposed orientations in opposing directions during the transition. Replacement of Leu-34 with a charged glutamic acid stabilizes FimH in the LAS conformation and replacement of Val-35 with glutamic acid traps FimH in the HAS conformation. The antibody is unable to trap the conformations if Leu-34 and Val-35 are replaced with a less bulky alanine. We propose that these residues act as molecular toggle switches and that the bound antibody imposes a steric block to their reorientation in either direction, thereby restricting concerted repacking of side chains that must occur to enable the conformational transition. Residues homologous to the FimH toggle switches are highly conserved across a diverse family of fimbrial adhesins. Replacement of predicted switch residues reveals that another E. coli adhesin, galactose-specific FmlH, is allosteric and can shift from an inactive to an active state. Our study shows that allosteric transitions in bacterial adhesins depend on toggle switch residues and that an antibody that blocks the switch effectively disables adhesive protein function.

Directed Evolution of a Selective and Sensitive Serotonin Sensor via Machine Learning

Serotonin plays a central role in cognition and is the target of most pharmaceuticals for psychiatric disorders. Existing drugs have limited efficacy; creation of improved versions will require better understanding of serotonergic circuitry, which has been hampered by our inability to monitor serotonin release and transport with high spatial and temporal resolution. We developed and applied a binding-pocket redesign strategy, guided by machine learning, to create a high-performance, soluble, fluorescent serotonin sensor (iSeroSnFR), enabling optical detection of millisecond-scale serotonin transients. We demonstrate that iSeroSnFR can be used to detect serotonin release in freely behaving mice during fear conditioning, social interaction, and sleep/wake transitions. We also developed a robust assay of serotonin transporter function and modulation by drugs. We expect that both machine-learning-guided binding-pocket redesign and iSeroSnFR will have broad utility for the development of other sensors and in vitro and in vivo serotonin detection, respectively.