Explore the vast range of titles available.

The problem of finding the correct set of partners for a given pair of interacting protein families based on multi-sequence alignments (MSAs) has received great attention over the years. Recently, the native contacts of two interacting proteins were shown to store the strongest mutual information (MI) signal to discriminate MSA concatenations with the largest fraction of correct pairings. Although that signal might be of practical relevance in the search for an effective heuristic to solve the problem, the number of MSA concatenations with near-native MI is large, imposing severe limitations. Here, a Genetic Algorithm that explores possible MSA concatenations according to a MI maximization criteria is shown to find degenerate solutions with two error sources, arising from mismatches among (i) similar and (ii) non-similar sequences. If mistakes made among similar sequences are disregarded, type-(i) solutions are found to resolve correct pairings at best true positive (TP) rates of 70%—far above the very same estimates in type-(ii) solutions. A machine learning classification algorithm helps to show further that differences between optimized solutions based on TP rates are not artificial and may have biological meaning associated with the three-dimensional distribution of the MI signal. Type-(i) solutions may therefore correspond to reliable results for predictive purposes, found here to be more likely obtained via MI maximization across protein systems having a minimum critical number of amino acid contacts on their interaction surfaces (N > 200).

The response of a membrane-bound Kv1.2 ion channel to an applied transmembrane potential has been studied using molecular dynamics simulations. Channel deactivation is shown to involve three intermediate states of the voltage sensor domain (VSD), and concomitant movement of helix S4 charges 10-15 Å along the bilayer normal; the latter being enabled by zipper-like sequential pairing of S4 basic residues with neighboring VSD acidic residues and membrane-lipid head groups. During the observed sequential transitions S4 basic residues pass through the recently discovered charge transfer center with its conserved phenylalanine residue, F²³³. Analysis indicates that the local electric field within the VSD is focused near the F²³³ residue and that it remains essentially unaltered during the entire process. Overall, the present computations provide an atomistic description of VSD response to hyperpolarization, add support to the sliding helix model, and capture essential features inferred from a variety of recent experiments.

Small-molecule modulators of diverse voltage-gated K + (Kv) channels may help treat a wide range of neurological disorders. However, developing effective modulators requires understanding of their mechanism of action. We apply an orthogonal approach to elucidate the mechanism of action of an imidazolidinedione derivative (AUT5), a highly selective positive allosteric modulator of Kv3.1 and Kv3.2 channels. AUT5 modulation involves positive cooperativity and preferential stabilization of the open state. The cryo-EM structure of the Kv3.1/AUT5 complex at a resolution of 2.5 Å reveals four equivalent AUT5 binding sites at the extracellular inter-subunit interface between the voltage-sensing and pore domains of the channel’s tetrameric assembly. Furthermore, we show that the unique extracellular turret regions of Kv3.1 and Kv3.2 essentially govern the selective positive modulation by AUT5. High-resolution apo and bound structures of Kv3.1 demonstrate how AUT5 binding promotes turret rearrangements and interactions with the voltage-sensing domain to favor the open conformation. To promote the development of effective small molecule modulators that may help treat diverse neuropsychiatric disorders, this study elucidates the mechanism of a specific positive modulator of neuronal potassium channels at near-atomic resolution.

The direct-site hypothesis assumes general anesthetics bind ion channels to impact protein equilibrium and function, inducing anesthesia. Despite advancements in the field, a first principle all-atom demonstration of this structure-function premise is still missing. We focus on the clinically used sevoflurane interaction to anesthetic-sensitive Kv1.2 mammalian channel to resolve if sevoflurane binds protein's well-characterized open and closed structures in a conformation-dependent manner to shift channel equilibrium. We employ an innovative approach relying on extensive docking calculations and free-energy perturbation of all potential binding sites revealed by the latter, and find sevoflurane binds open and closed structures at multiple sites under complex saturation and concentration effects. Results point to a non-trivial interplay of site and conformation-dependent modes of action involving distinct binding sites that increase channel open-probability at diluted ligand concentrations. Given the challenge in exploring more complex processes potentially impacting channel-anesthetic interaction, the result is revealing as it demonstrates the process of multiple anesthetic binding events alone may account for open-probability shifts recorded in measurements.

Membrane proteins are primary targets for most therapeutic indications in cancer and neurological diseases, binding over 50% of all known small molecule drugs. Understanding how such ligands impact membrane proteins requires knowledge on the molecular structure of ligand binding, a reasoning that has driven relentless efforts in drug discovery and translational research. Binding of small ligands appears however highly complex involving interaction to multiple transmembrane protein sites featuring single or multiple occupancy states. Within this scenario, looking for new developments in the field, we investigate the concentration-dependent binding of ligands to multiple saturable sites in membrane proteins. The study relying on docking and free-energy perturbation provides us with an extensive description of the probability density of protein-ligand states that allows for computation of thermodynamic properties of interest. It also provides one- and three-dimensional spatial descriptions for the ligand density across the protein-membrane system which can be of interest for structural purposes. Illustration and discussion of the results are shown for binding of the general anesthetic sevoflurane against Kv1.2, a mammalian ion channel for which experimental data are available.

The recent crystal structure of Orai, the pore unit of a calcium release-activated calcium (CRAC) channel, is used as the starting point for molecular dynamics and free-energy calculations designed to probe this channel’s conduction properties. In free molecular dynamics simulations, cations localize preferentially at the extracellular channel entrance near the ring of Glu residues identified in the crystal structure, whereas anions localize in the basic intracellular half of the pore. To begin to understand ion permeation, the potential of mean force (PMF) was calculated for displacing a single Na ⁺ ion along the pore of the CRAC channel. The computed PMF indicates that the central hydrophobic region provides the major hindrance for ion diffusion along the permeation pathway, thereby illustrating the nonconducting nature of the crystal structure conformation. Strikingly, further PMF calculations demonstrate that the mutation V174A decreases the free energy barrier for conduction, rendering the channel effectively open. This seemingly dramatic effect of mutating a nonpolar residue for a smaller nonpolar residue in the pore hydrophobic region suggests an important role for the latter in conduction. Indeed, our computations show that even without significant channel-gating motions, a subtle change in the number of pore waters is sufficient to reshape the local electrostatic field and modulate the energetics of conduction, a result that rationalizes recent experimental findings. The present work suggests the activation mechanism for the wild-type CRAC channel is likely regulated by the number of pore waters and hence pore hydration governs the conductance.

The X-ray structure of the bacterial voltage-gated sodium channel NavAb has been reported in a conformation with a closed conduction pore. Comparison between this structure and the activated-open and resting-closed structures of the voltage-gated Kv1.2 potassium channel suggests that the voltage-sensor domains (VSDs) of the reported structure are not fully activated. Using the aforementioned structures of Kv1.2 as templates, molecular dynamics simulations are used to identify analogous functional conformations of NavAb. Specifically, starting from the NavAb crystal structure, conformations of the membrane-bound channel are sampled along likely pathways for activation of the VSD and opening of the pore domain. Gating charge computations suggest that a structural rearrangement comparable to that occurring between activated-open and resting-closed states is required to explain experimental values of the gating charge, thereby confirming that the reported VSD structure is likely an intermediate along the channel activation pathway. Our observation that the X-ray structure exhibits a low pore domain-opening propensity further supports this notion. The present molecular dynamics study also identifies conformations of NavAb that are seemingly related to the resting-closed and activated-open states. Our findings are consistent with recent structural and functional studies of the orthologous channels NavRh, NaChBac, and NavMs and offer possible structures for the functionally relevant conformations of NavAb.

Significance A central paradigm in the subject of cellular excitability is whether or not a dynamic membrane voltage field contributes to the operation of voltage sensor proteins. This issue has challenged an entire generation of scientists because static and dynamic contributions of the field cannot be discriminated from electrophysiology measurements. The dilemma imposes that structure-based calculations are required to solve the issue. By following this approach, we report that reshaping or dynamics of the field does not impact significantly the domain energetics or operation. As such, the functional diversity of voltage sensors observed in nature originates from variations in the detailed atom–atom interactions of the domain (chemical free energy) rather than in its coupling with the membrane voltage. A dynamic transmembrane voltage field has been suggested as an intrinsic element in voltage sensor (VS) domains. Here, the dynamic field contribution to the VS energetics was analyzed via electrostatic calculations applied to a number of atomistic structures made available recently. We find that the field is largely static along with the molecular motions of the domain, and more importantly, it is minimally modified across VS variants. This finding implies that sensor domains transfer approximately the same amount of gating charges when moving the electrically charged S4 helix between fixed microscopic configurations. Remarkably, the result means that the observed operational diversity of the domain, including the extension, rate, and voltage dependence of the S4 motion, as dictated by the free energy landscape theory, must be rationalized in terms of dominant variations of its chemical free energy.

Many voltage-gated potassium (Kv) channels display a time-dependent phenomenon called C-type inactivation, whereby prolonged activation by voltage leads to the inhibition of ionic conduction, a process that involves a conformational change at the selectivity filter toward a non-conductive state. Recently, a high-resolution structure of a strongly inactivated triple-mutant channel kv1.2-kv2.1-3m revealed a novel conformation of the selectivity filter that is dilated at its outer end, distinct from the well-characterized conductive state. While the experimental structure was interpreted as the elusive non-conductive state, our molecular dynamics simulations and electrophysiological measurements show that the dilated filter of kv1.2-kv2.1-3m is conductive and, as such, cannot completely account for the inactivation of the channel observed in the structural experiments. The simulation shows that an additional conformational change, implicating isoleucine residues at position 398 along the pore lining segment S6, is required to effectively block ion conduction. The I398 residues from the four subunits act as a state-dependent hydrophobic gate located immediately beneath the selectivity filter. These observations are corroborated by electrophysiological experiments showing that ion permeation can be resumed in the kv1.2-kv2.1-3m channel when I398 is mutated to an asparagine—a mutation that does not abolish C-type inactivation since digitoxin (AgTxII) fails to block the ionic permeation of kv1.2-kv2.1-3m_I398N. As a critical piece of the C-type inactivation machinery, this structural feature is the potential target of a broad class of quaternary ammonium (QA) blockers and negatively charged activators thus opening new research directions toward the development of drugs that specifically modulate gating states of Kv channels.

Studies on the interaction sites of peptide toxins and ion channels typically involve site-directed mutations in toxins. However, natural mutant toxins exist among them, offering insights into how the evolutionary process has conserved crucial sequences for activities and molecular target selection. In this study, we present a comparative investigation using electrophysiological approaches and computational analysis between two alpha toxins from evolutionarily close scorpion species of the genus Tityus, namely, Tst3 and Ts3 from T. stigmurus and T. serrulatus, respectively. These toxins exhibit three natural substitutions near the C-terminal region, which is directly involved in the interaction between alpha toxins and Nav channels. Additionally, we characterized the activity of the Tst3 toxin on Nav1.1-Nav1.7 channels. The three natural changes between the toxins did not alter sensitivity to Nav1.4, maintaining similar intensities regarding their ability to alter opening probabilities, delay fast inactivation, and induce persistent currents. Computational analysis demonstrated a preference for the down conformation of VSD4 and a shift in the conformational equilibrium towards this state. This illustrates that the sequence of these toxins retained the necessary information, even with alterations in the interaction site region. Through electrophysiological and computational analyses, screening of the Tst3 toxin on sodium isoform revealed its classification as a classic α-NaTx with a broad spectrum of activity. It effectively delays fast inactivation across all tested isoforms. Structural analysis of molecular energetics at the interface of the VSD4-Tst3 complex further confirmed this effect.