Daniele Granata

Fluorine-19 NMR and computational quantification of isoflurane binding to the voltage-gated sodium channel NaChBac

Significance How general anesthetics modulate the function of voltage-gated sodium (NaV) channels remains a mystery. Here, strategic placements of 19F probes, guided by molecular dynamics simulations, allowed for high-resolution NMR quantitation of the volatile anesthetic isoflurane binding to the bacterial Nav channel NaChBac. The data provided experimental evidence showing that channel blockade at the base of the ion selectivity filter and the restricted pivot motion at the S4–S5 linker and the P2–S6 helix hinge underlie the action of isoflurane on NaChBac. The results contribute to a better understanding of the molecular mechanisms of general anesthesia. Voltage-gated sodium channels (NaV) play an important role in general anesthesia. Electrophysiology measurements suggest that volatile anesthetics such as isoflurane inhibit NaV by stabilizing the inactivated state or altering the inactivation kinetics. Recent computational studies suggested the existence of multiple isoflurane binding sites in NaV, but experimental binding data are lacking. Here we use site-directed placement of 19F probes in NMR experiments to quantify isoflurane binding to the bacterial voltage-gated sodium channel NaChBac. 19F probes were introduced individually to S129 and L150 near the S4–S5 linker, L179 and S208 at the extracellular surface, T189 in the ion selectivity filter, and all phenylalanine residues. Quantitative analyses of 19F NMR saturation transfer difference (STD) spectroscopy showed a strong interaction of isoflurane with S129, T189, and S208; relatively weakly with L150; and almost undetectable with L179 and phenylalanine residues. An orientation preference was observed for isoflurane bound to T189 and S208, but not to S129 and L150. We conclude that isoflurane inhibits NaChBac by two distinct mechanisms: (i) as a channel blocker at the base of the selectivity filter, and (ii) as a modulator to restrict the pivot motion at the S4–S5 linker and at a critical hinge that controls the gating and inactivation motion of S6.

Voltage-Gated Sodium Channels: Evolutionary History and Distinctive Sequence Features.

Voltage-gated sodium channels (Nav) are responsible for the rising phase of the action potential. Their role in electrical signal transmission is so relevant that their emergence is believed to be one of the crucial factors enabling development of nervous system. The presence of voltage-gated sodium-selective channels in bacteria (BacNav) has raised questions concerning the evolutionary history of the ones in animals. Here we review some of the milestones in the field of Nav phylogenetic analysis and discuss some of the most important sequence features that distinguish these channels from voltage-gated potassium channels and transient receptor potential channels.

Voltage-Gated Sodium Channels

Voltage-gated sodium channels (Nav) are responsible for the rising phase of the action potential. Their role in electrical signal transmission is so relevant that their emergence is believed to be one of the crucial factors enabling development of nervous system. The presence of voltage-gated sodium-selective channels in bacteria (BacNav) has raised questions concerning the evolutionary history of the ones in animals. Here we review some of the milestones in the field of Nav phylogenetic analysis and discuss some of the most important sequence features that distinguish these channels from voltage-gated potassium channels and transient receptor potential channels.

Native fold and docking pose discrimination by the same residue‐based scoring function

Structure prediction and quality assessment are crucial steps in modeling native protein conformations. Statistical potentials are widely used in related algorithms, with different parametrizations typically developed for different contexts such as folding protein monomers or docking protein complexes. Here, we describe BACH‐SixthSense, a single residue‐based statistical potential that can be successfully employed in both contexts. BACH‐SixthSense shares the same approach as BACH, a knowledge‐based potential originally developed to score monomeric protein structures. A term that penalizes steric clashes as well as the distinction between polar and apolar sidechain‐sidechain contacts are crucial novel features of BACH‐SixthSense. The performance of BACH‐SixthSense in discriminating correctly the native structure among a competing set of decoys is significantly higher than other state‐of‐the‐art scoring functions, that were specifically trained for a single context, for both monomeric proteins (QMEAN, Rosetta, RF_CB_SRS_OD, benchmarked on CASP targets) and protein dimers (IRAD, Rosetta, PIE*PISA, HADDOCK, FireDock, benchmarked on 14 CAPRI targets). The performance of BACH‐SixthSense in recognizing near‐native docking poses within CAPRI decoy sets is good as well. Proteins 2015; 83:621–630.

Patterns of coevolving amino acids unveil structural and dynamical domains

Significance Patterns of pairwise correlations in sequence alignments can be used to reconstruct the network of residue-residue contacts and thus the three-dimensional structure of proteins. Less explored, and yet extremely intriguing, is the functional relevance of such coevolving networks: Do they encode for the collective motions occurring in proteins at thermal equilibrium? Here, by combining coevolutionary coupling analysis with a state-of-the-art dimensionality reduction approach, we show that the network of pairwise evolutionary couplings can be analyzed to reveal communities of amino acids, which we term “evolutionary domains,” that are in striking agreement with the quasi-rigid protein domains obtained from elastic network models and molecular dynamics simulations. Patterns of interacting amino acids are so preserved within protein families that the sole analysis of evolutionary comutations can identify pairs of contacting residues. It is also known that evolution conserves functional dynamics, i.e., the concerted motion or displacement of large protein regions or domains. Is it, therefore, possible to use a pure sequence-based analysis to identify these dynamical domains? To address this question, we introduce here a general coevolutionary coupling analysis strategy and apply it to a curated sequence database of hundreds of protein families. For most families, the sequence-based method partitions amino acids into a few clusters. When viewed in the context of the native structure, these clusters have the signature characteristics of viable protein domains: They are spatially separated but individually compact. They have a direct functional bearing too, as shown for various reference cases. We conclude that even large-scale structural and functionally related properties can be recovered from inference methods applied to evolutionary-related sequences. The method introduced here is available as a software package and web server (spectrus.sissa.it/spectrus-evo_webserver).

Propofol inhibits the voltage-gated sodium channel NaChBac at multiple sites

Voltage-gated sodium (NaV) channels are important targets of general anesthetics, including the intravenous anesthetic propofol. Electrophysiology studies on the prokaryotic NaV channel NaChBac have demonstrated that propofol promotes channel activation and accelerates activation-coupled inactivation, but the molecular mechanisms of these effects are unclear. Here, guided by computational docking and molecular dynamics simulations, we predict several propofol-binding sites in NaChBac. We then strategically place small fluorinated probes at these putative binding sites and experimentally quantify the interaction strengths with a fluorinated propofol analogue, 4-fluoropropofol. In vitro and in vivo measurements show that 4-fluoropropofol and propofol have similar effects on NaChBac function and nearly identical anesthetizing effects on tadpole mobility. Using quantitative analysis by 19F-NMR saturation transfer difference spectroscopy, we reveal strong intermolecular cross-relaxation rate constants between 4-fluoropropofol and four different regions of NaChBac, including the activation gate and selectivity filter in the pore, the voltage sensing domain, and the S4–S5 linker. Unlike volatile anesthetics, 4-fluoropropofol does not bind to the extracellular interface of the pore domain. Collectively, our results show that propofol inhibits NaChBac at multiple sites, likely with distinct modes of action. This study provides a molecular basis for understanding the net inhibitory action of propofol on NaV channels.

Ion Channel Sensing: Are Fluctuations the Crux of the Matter?

The nonselective cation channel TRPV1 is responsible for transducing noxious stimuli into action potentials propagating through peripheral nerves. It is activated by temperatures greater than 43 °C, while remaining completely nonconductive at temperatures lower than this threshold. The origin of this sharp response, which makes TRPV1 a biological temperature sensor, is not understood. Here we used molecular dynamics simulations and free energy calculations to characterize the molecular determinants of the transition between nonconductive and conductive states. We found that hydration of the pore and thus ion permeation depends critically on the polar character of its molecular surface: in this narrow hydrophobic enclosure, the motion of a polar side-chain is sufficient to stabilize either the dry or wet state. The conformation of this side-chain is in turn coupled to the hydration state of four peripheral cavities, which undergo a dewetting transition at the activation temperature.

TRPV1 activation relies on hydration/dehydration of nonpolar cavities

TRPV1 promotes cationic currents across cellular membranes in response to multiple stimuli such as increased temperature or pressure, binding of chemicals, low pH and voltage. The molecular underpinnings of TRPV1 gating, in particular the mechanism of temperature sensitivity, are still largely unknown. Here, we used molecular simulations and electrophysiology to shed light on the closed to open transition. Specifically, we found that gating of TRPV1 relies on the motion of an evolutionarily conserved amino acid (N676) in the middle of the S6 helix. On rotation, the side chain of this asparagine faces either the central pore or the S4-S5 linker. Only in the former case is the central pore hydrated and thus conductive. Interestingly, when N676 rotates toward the linker, we observe hydration of four so far unreported small nonpolar cavities. Based on these findings, we propose a model for TRPV1 gating involving the dynamic hydration of these four cavities. Free energy calculations indicate that this gating mechanisms is markedly temperature dependent favoring the open state at high temperature. On the basis of this model, which is able to rationalize a wealth of seemingly conflicting and/or unrelated experimental observations, we predicted the behavior of a single residue mutant, F580Y, the consequences of which are confirmed experimentally and give support to the model.

Cross-kingdom auxiliary subunit modulation of a voltage-gated sodium channel

Voltage-gated, sodium ion–selective channels (NaV) generate electrical signals contributing to the upstroke of the action potential in animals. NaVs are also found in bacteria and are members of a larger family of tetrameric voltage-gated channels that includes CaVs, KVs, and NaVs. Prokaryotic NaVs likely emerged from a homotetrameric Ca2+-selective voltage-gated progenerator, and later developed Na+ selectivity independently. The NaV signaling complex in eukaryotes contains auxiliary proteins, termed beta (β) subunits, which are potent modulators of the expression profiles and voltage-gated properties of the NaV pore, but it is unknown whether they can functionally interact with prokaryotic NaV channels. Herein, we report that the eukaryotic NaVβ1-subunit isoform interacts with and enhances the surface expression as well as the voltage-dependent gating properties of the bacterial NaV, NaChBac in Xenopus oocytes. A phylogenetic analysis of the β-subunit gene family proteins confirms that these proteins appeared roughly 420 million years ago and that they have no clear homologues in bacterial phyla. However, a comparison between eukaryotic and bacterial NaV structures highlighted the presence of a conserved fold, which could support interactions with the β-subunit. Our electrophysiological, biochemical, structural, and bioinformatics results suggests that the prerequisites for β-subunit regulation are an evolutionarily stable and intrinsic property of some voltage-gated channels.

A consistent picture of TRPV1 activation emerges from molecular simulations and experiments

Although experimental structures of both the TRPV1 closed and open states are available, the conformational changes occurring in the pore domain and resulting in ionic conduction remain elusive. Here, we use bioinformatics analyses, molecular dynamics simulations and site-directed mutagenesis to shed light on this issue and suggest a possible molecular mechanism for TRPV1 activation. In light of our hypothesis, we re-examine the results of the previously published water accessibility and mutagenesis experiments, and analyze the newly available structures of TRP and other evolutionary related ion channels. Overall, we show that several independent lines of evidence corroborate our hypothesis, which highlights the rotation of a conserved asparagine toward the pore and the dehydration of hydrophobic cavities as key components of TRPV1 activation. Importantly, this molecular mechanism provides also a rationale for the coupling between the TRPV1 upper and lower gates.