Concentration-Dependent Binding of Small Ligands to Multiple Saturable Sites in Membrane Proteins

The binding and mechanism of a positive allosteric modulator of Kv3 channels

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.

The Modulation by Anesthetics and Analgesics of Respiratory Rhythm in the Nervous System

Rhythmic eupneic breathing in mammals depends on the coordinated activities of the neural system that sends cranial and spinal motor outputs to respiratory muscles. These outputs modulate lung ventilation and adjust respiratory airflow, which depends on the upper airway patency and ventilatory musculature. Anesthetics are widely used in clinical practice worldwide. In addition to clinically necessary pharmacological effects, respiratory depression is a critical side effect induced by most general anesthetics. Therefore, understanding how general anesthetics modulate the respiratory system is important for the development of safer general anesthetics. Currently used volatile anesthetics and most intravenous anesthetics induce inhibitory effects on respiratory outputs. Various general anesthetics produce differential effects on respiratory characteristics, including the respiratory rate, tidal volume, airway resistance, and ventilatory response. At the cellular and molecular levels, the mechanisms underlying anesthetic-induced breathing depression mainly include modulation of synaptic transmission of ligand-gated ionotropic receptors (e.g., γ-aminobutyric acid, N-methyl-D-aspartate, and nicotinic acetylcholine receptors) and ion channels (e.g., voltage-gated sodium, calcium, and potassium channels, two-pore domain potassium channels, and sodium leak channels), which affect neuronal firing in brainstem respiratory and peripheral chemoreceptor areas. The present review comprehensively summarizes the modulation of the respiratory system by clinically used general anesthetics, including the effects at the molecular, cellular, anatomic, and behavioral levels. Specifically, analgesics, such as opioids, which cause respiratory depression and the "opioid crisis", are discussed. Finally, underlying strategies of respiratory stimulation that target general anesthetics and/or analgesics are summarized.

Alchemical Free Energy Calculations on Membrane-Associated Proteins

Membrane proteins have diverse functions within cells and are well-established drug targets. The advances in membrane protein structural biology have revealed drug and lipid binding sites on membrane proteins, while computational methods such as molecular simulations can resolve the thermodynamic basis of these interactions. Particularly, alchemical free energy calculations have shown promise in the calculation of reliable and reproducible binding free energies of protein-ligand and protein-lipid complexes in membrane-associated systems. In this review, we present an overview of representative alchemical free energy studies on G-protein-coupled receptors, ion channels, transporters as well as protein-lipid interactions, with emphasis on best practices and critical aspects of running these simulations. Additionally, we analyze challenges and successes when running alchemical free energy calculations on membrane-associated proteins. Finally, we highlight the value of alchemical free energy calculations calculations in drug discovery and their applicability in the pharmaceutical industry.

Allosteric Modulation of Membrane Proteins by Small Low-Affinity Ligands

Membrane proteins may respond to a variety of ligands under an applied external stimulus. These ligands include small low-affinity molecules that account for functional effects in the mM range. Understanding the modulation of protein function by low-affinity ligands requires characterizing their atomic-level interactions under dilution, challenging the current resolution of theoretical and experimental routines. Part of the problem derives from the fact that small low-affinity ligands may interact with multiple sites of a membrane protein in a highly degenerate manner to a degree that it is better conceived as a partition phenomenon, hard to track at the molecular interface of the protein. Looking for new developments in the field, we rely on the classic two-state Boltzmann model to devise a novel theoretical description of the allosteric modulation mechanism of membrane proteins in the presence of small low-affinity ligands and external stimuli. Free energy stability of the partition process and its energetic influence on the protein coupling with the external stimulus are quantified. The outcome is a simple formulation that allows the description of the equilibrium shifts of the protein in terms of the grand-canonical partition function of the ligand at dilute concentrations. The model's predictions of the spatial distribution and response probability shift across a variety of ligand concentrations, and thermodynamic conjugates can be directly compared to macroscopic measurements, making it especially useful to interpret experimental data at the atomic level. Illustration and discussion of the theory is shown in the context of general anesthetics and voltage-gated channels for which structural data are available.

Metabolites of Cannabigerol Generated by Human Cytochrome P450s Are Bioactive

The phytocannabinoid cannabigerol (CBG) is the central biosynthetic precursor to many cannabinoids, including Δ⁹-tetrahydrocannabinol (THC) and cannabidiol (CBD). Though the use of CBG has recently witnessed a widespread surge because of its beneficial health effects and lack of psychoactivity, its metabolism by human cytochrome P450s is largely unknown. Herein, we describe comprehensive in vitro and in vivo cytochrome P450 (CYP)-mediated metabolic studies of CBG, ranging from liquid chromatography tandem mass spectrometry-based primary metabolic site determination, synthetic validation, and kinetic behavior using targeted mass spectrometry. These investigations revealed that cyclo-CBG, a recently isolated phytocannabinoid, is the major metabolite that is rapidly formed by selected human cytochrome P450s (CYP2J2, CYP3A4, CYP2D6, CYP2C8, and CYP2C9). Additionally, in vivo studies with mice administered with CBG supported these studies, where cyclo-CBG is the major metabolite as well. Spectroscopic binding studies along with docking and modeling of the CBG molecule near the heme in the active site of P450s confirmed these observations, pointing at the preferred site selectivity of CBG metabolism at the prenyl chain over other positions. Importantly, we found out that CBG and its oxidized CBG metabolites reduced inflammation in BV2 microglial cells stimulated with LPS. Overall, combining enzymological studies, mass spectrometry, and chemical synthesis, we showcase that CBG is rapidly metabolized by human P450s to form oxidized metabolites that are bioactive.

Concentration-Dependent Thermodynamic Analysis of the Partition Process of Small Ligands into Proteins

In the category of functional low-affinity interactions, small ligands may interact with multiple protein sites in a highly degenerate manner. Better conceived as a partition phenomenon at the molecular interface of proteins, such low-affinity interactions appear to be hidden to our current experimental resolution making their structural and functional characterization difficult in the low concentration regime of physiological processes. Characterization of the partition phenomenon under higher chemical forces could be a relevant strategy to tackle the problem provided the results can be scaled back to the low concentration range. Far from being trivial, such scaling demands a concentration-dependent understanding of self-interactions of the ligands, structural perturbations of the protein, among other molecular effects. Accordingly, we elaborate a novel and detailed concentration-dependent thermodynamic analysis of the partition process of small ligands aiming at characterizing the stability and structure of the dilute phenomenon from high concentrations. In analogy to an “aggregate” binding constant of a small molecule over multiple sites of a protein receptor, the model defines the stability of the process as a macroscopic equilibrium constant for the partition number of ligands that can be used to analyze biochemical and functional data of two-component systems driven by low-affinity interactions. Acquisition of such modeling-based structural information is expected to be highly welcome by revealing more traceable protein-binding spots for non-specific ligands.

A Potential Mechanism of Sodium Channel Mediating the General Anesthesia Induced by Propofol

General anesthesia has revolutionized healthcare over the past 200 years and continues to show advancements. However, many phenomena induced by general anesthetics including paradoxical excitation are still poorly understood. Voltage-gated sodium channels (Na_V) were believed to be one of the proteins targeted during general anesthesia. Based on electrophysiological measurements before and after propofol treatments of different concentrations, we mathematically modified the Hodgkin–Huxley sodium channel formulations and constructed a thalamocortical model to investigate the potential roles of Na_V. The ion channels of individual neurons were modeled using the Hodgkin–Huxley type equations. The enhancement of propofol-induced GABAa current was simulated by increasing the maximal conductance and the time-constant of decay. Electroencephalogram (EEG) was evaluated as the post-synaptic potential from pyramidal (PY) cells. We found that a left shift in activation of Na_V was induced primarily by a low concentration of propofol (0.3–10 μM), while a left shift in inactivation of Na_V was induced by an increasing concentration (0.3–30 μM). Mathematical simulation indicated that a left shift of Na_V activation produced a Hopf bifurcation, leading to cell oscillations. Left shift of Na_V activation around a value of 5.5 mV in the thalamocortical models suppressed normal bursting of thalamocortical (TC) cells by triggering its chaotic oscillations. This led to irregular spiking of PY cells and an increased frequency in EEG readings. This observation suggests a mechanism leading to paradoxical excitation during general anesthesia. While a left shift in inactivation led to light hyperpolarization in individual cells, it inhibited the activity of the thalamocortical model after a certain depth of anesthesia. This finding implies that high doses of propofol inhibit the network partly by accelerating Na_V toward inactivation. Additionally, this result explains why the application of sodium channel blockers decreases the requirement for general anesthetics. Our study provides an insight into the roles that Na_V plays in the mechanism of general anesthesia. Since the activation and inactivation of Na_V are structurally independent, it should be possible to avoid side effects by state-dependent binding to the Na_V to achieve precision medicine in the future.

Open-channel blocking action of volatile anaesthetics desflurane and sevoflurane on human voltage-gated Kv 1.5 channel

BACKGROUND AND PURPOSE

Volatile anaesthetics have been shown to differentially modulate mammalian Shaker-related voltage-gated potassium (K_v 1.x) channels. This study was designed to investigate molecular and cellular mechanisms underlying the modulatory effects of desflurane or sevoflurane on human K_v 1.5 (hK_v 1.5) channels.

EXPERIMENTAL APPROACH

Thirteen single-point mutations were constructed within pore domain of hK_v 1.5 channel using site-directed mutagenesis. The effects of desflurane or sevoflurane on heterologously expressed wild-type and mutant hK_v 1.5 channels were examined by whole-cell patch-clamp technique. A computer simulation was conducted to predict the docking pose of desflurane or sevoflurane within hK_v 1.5 channel.

KEY RESULTS

Both desflurane and sevoflurane increased hK_v 1.5 current at mild depolarizations but decreased it at strong depolarizations, indicating that these anaesthetics produce both stimulatory and inhibitory actions on hK_v 1.5 channels. The inhibitory effect of desflurane or sevoflurane on hK_v 1.5 channels arose primarily from its open-channel blocking action. The inhibitory action of desflurane or sevoflurane on hK_v 1.5 channels was significantly attenuated in T480A, V505A, and I508A mutant channels, compared with wild-type channel. Computational docking simulation predicted that desflurane or sevoflurane resides within the inner cavity of channel pore and has contact with Thr479, Thr480, Val505, and Ile508.

CONCLUSION AND IMPLICATIONS

Desflurane and sevoflurane exert an open-channel blocking action on hK_v 1.5 channels by functionally interacting with specific amino acids located within the channel pore. This study thus identifies a novel molecular basis mediating inhibitory modulation of hK_v 1.5 channels by desflurane and sevoflurane.

Donepezil Inhibits Acetylcholinesterase via Multiple Binding Modes at Room Temperature

Donepezil is a second generation acetylcholinesterase (AChE) inhibitor for treatment of Alzheimer's disease (AD). AChE is important for neurotransmission at neuromuscular junctions and cholinergic brain synapses by hydrolyzing acetylcholine into acetate and choline. In vitro data support that donepezil is a reversible, mixed competitive and noncompetitive inhibitor of AChE. The experimental fact then suggests a more complex binding mechanism beyond the molecular view in X-ray models resolved at cryogenic temperatures that show a unique binding mode of donepezil in the active site of the enzyme. Aiming at clarifying the mechanism behind that mixed competitive and noncompetitive nature of the inhibitor, we have applied molecular dynamics (MD) simulations and docking and free-energy calculations to investigate microscopic details and energetics of donepezil association for conditions of substrate-free and -bound states of the enzyme. Liquid-phase MD simulation at room temperature shows AChE transits between "open" and "closed" conformations to control accessibility to the active site and ligand binding. As shown by docking and free-energy calculations, association of donepezil involves its reversible axial displacement and reorientation in the active site of the enzyme, assisted by water molecules. Donepezil binds equally well the main-door anionic binding site PAS, the acyl pocket, and the catalytic site CAS by respectively adopting outward-inward-inward orientations regardless of substrate occupancy-the overall stability of that reaction process depends however on co-occupancy of the enzyme being preferential for its substrate-free state. All together, our findings support a physiologically relevant mechanism of AChE inhibition by donepezil involving multistable interactions modes at the molecular origin of the inhibitor's activity.

Atomistic Model for Simulations of the Sedative Hypnotic Drug 2,2,2-Trichloroethanol

2,2,2-Trichloroethanol (TCE) is the active form of the sedative hypnotic drug chloral hydrate, one of the oldest sleep medications in the market. Understanding of TCE's action mechanisms to its many targets, particularly within the ion channel family, could benefit from the state-of-the-art computational molecular studies. In this direction, we employed de novo modeling aided by the force field toolkit to develop CHARMM36-compatible TCE parameters. The classical potential energy function was calibrated targeting molecular conformations, local interactions with water molecules, and liquid bulk properties. Reference data comes from both tabulated thermodynamic properties and ab initio calculations at the MP2 level. TCE solvation free energy calculations in water and oil reproduce a lipophilic, yet nonhydrophobic, behavior. Indeed, the potential mean force profile for TCE partition through the phospholipid bilayer reveals the sedative's preference for the interfacial region. The calculated partition coefficient also matches experimental measures. Further validation of the proposed parameters is supported by the model's ability to recapitulate quenching experiments demonstrating TCE binding to bovine serum albumin.

Binding of the general anesthetic sevoflurane to ion channels

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.

General anesthetics are central to modern medicine, yet their microscopic mechanism of action is still unknown. Here, we demonstrate that a clinically used anesthetic, sevoflurane, binds the mammalian voltage-gated potassium channel Kv1.2 effecting a shift in its open probability, even at low concentrations. The results, supported by recent experimental measurements, are promising as they demonstrate that the molecular process of direct binding of anesthetic to ion channels play a relevant role in anesthesia.

Location and Character of Volatile General Anesthetics Binding Sites in the Transmembrane Domain of TRPV1

It has been proposed that general anesthesia results from direct multisite interactions with multiple and diverse ion channels in the brain. An understanding of the mechanisms by which general anesthetics modulate ion channels is essential to clarify their underlying behavior and their role in reversible immobilization and amnesia. Despite the fact that volatile general anesthetics are drugs that primarily induce insensitivity to pain, they have been reported to sensitize and active the vanilloid-1 receptor, TRPV1, which is known to mediate the response of the nervous system to certain harmful stimuli and which plays a crucial role in the pain pathway. Currently, the mechanism of action of anesthetics is unknown and the precise molecular sites of interaction have not been identified. Here, using ∼2.5 μs of classical molecular dynamics simulations and metadynamics, we explore these enigmas. Binding sites are identified and the strength of the association is further characterized using alchemical free-energy calculations. Anesthetic binding/unbinding proceeds primarily through a membrane-embedded pathway, and subsequently, a complex scenario is established involving multiple binding sites featuring single or multiple occupancy states of two small volatile drugs. One of the five anesthetic binding sites reported was previously identified experimentally, and another one, importantly, is identical to that of capsaicin, one of the chemical stimuli that activate TRPV1. However, in contrast to capsaicin, isoflurane and chloroform binding free-energies render modest to no association compared to capsaicin, suggesting a different activation mechanism. Uncovering chloroform and isoflurane modulatory sites will further our understanding of the TRPV1 molecular machinery and open the possibility of developing site-specific drugs.

Electrophysiological Analysis of Voltage-Gated Ion Channel Modulation by General Anesthetics

Voltage-gated ion channels (VGICs) of excitable tissues are emerging as targets likely involved in both the therapeutic and toxic effects of inhaled and intravenous general anesthetics. Whereas sevoflurane and propofol inhibit voltage-gated Na⁺ channels (Navs), sevoflurane potentiates certain voltage-gated K⁺ channels (Kvs). The combination of these effects would dampen neural excitability and, therefore, might contribute to the clinical endpoints of general anesthesia. As the body of work regarding the interaction of general anesthetics with VGICs continues to grow, a multidisciplinary approach involving functional, biochemical, structural, and computational techniques, many of which are detailed in other chapters, has increasingly become necessary to solve the molecular mechanism of general anesthetic action on VGICs. Here, we focus on electrophysiological and modeling approaches and methodologies to describe how our work has elucidated the biophysical basis of the inhibition Navs by propofol and the potentiation of Kvs by sevoflurane.