Bioequivalence Study of Velpatasvir/Sofosbuvir Oral Coated Tablets in Healthy Volunteers Under Fasting Conditions.

This study was conducted as a single-site, open-label, randomized, replicated crossover trial with 4 treatment periods. The aim was to evaluate the bioequivalence of a generic test drug containing velpatasvir and sofosbuvir compared to an established brand-name medication in healthy White subjects under fasting conditions. Blood samples were collected at specified intervals up to 72 hours after dosing to measure the concentrations of velpatasvir and sofosbuvir using a certified high-performance liquid chromatography with tandem mass spectrometry method. The bioequivalence of the 2 formulations was confirmed when statistical analysis showed that confidence intervals for the log-transformed peak concentration and area under the concentration-time curve from time 0 to the last quantifiable sample were within an acceptable range from 80% to 125%. Criteria for bioequivalence were met for both area under the concentration-time curve from time 0 until the last quantifiable sample and peak concentration parameters. No adverse effects were reported during this trial in both groups.

Clinical Pharmacology of GP40321 (Insulin Glulisine Biosimilar): Pharmacokinetic and Pharmacodynamic Comparability in a Hyperinsulinemic-Euglycemic Clamp Procedure.

The aim of the study was to compare the pharmacokinetics (PK) and pharmacodynamics (PD) of T-glu (GP40321, test drug), and reference insulin glulisine in a hyperinsulinemic-euglycemic clamp procedure. During this study, 34 healthy male volunteers underwent the hyperinsulinemic-euglycemic clamp procedure following subcutaneous 0.3 U/kg injection of T-glu or reference insulin glulisine in a randomized, double-blind, crossover study. Plasma glucose levels were monitored every 5 minutes for 8 hours. Glucose infusion rate adjustment was based on the blood glucose measurements. Evaluation of PD was performed using the glucose infusion rate values, while PK was calculated using insulin concentrations measured via enzyme-linked immunosorbent assay. The study results showed that the 90% CI for the geometric mean ratios of primary PK and PD of T-glu and reference insulin glulisine were within 80%-125% comparability limits, and that the safety profiles were comparable. PK, PD, and safety similarity of T-glu and reference insulin glulisine was demonstrated.

Endocannabinoids enhance hKV7.1/KCNE1 channel function and shorten the cardiac action potential and QT interval.

BACKGROUND: Genotype-positive patients who suffer from the cardiac channelopathy Long QT Syndrome (LQTS) may display a spectrum of clinical phenotypes, with often unknown causes. Therefore, there is a need to identify factors influencing disease severity to move towards an individualized clinical management of LQTS. One possible factor influencing the disease phenotype is the endocannabinoid system, which has emerged as a modulator of cardiovascular function. In this study, we aim to elucidate whether endocannabinoids target the cardiac voltage-gated potassium channel KV7.1/KCNE1, which is the most frequently mutated ion channel in LQTS. METHODS: We used two-electrode voltage clamp, molecular dynamics simulations and the E4031 drug-induced LQT2 model of ex-vivo guinea pig hearts. FINDINGS: We found a set of endocannabinoids that facilitate channel activation, seen as a shifted voltage-dependence of channel opening and increased overall current amplitude and conductance. We propose that negatively charged endocannabinoids interact with known lipid binding sites at positively charged amino acids on the channel, providing structural insights into why only specific endocannabinoids modulate KV7.1/KCNE1. Using the endocannabinoid ARA-S as a prototype, we show that the effect is not dependent on the KCNE1 subunit or the phosphorylation state of the channel. In guinea pig hearts, ARA-S was found to reverse the E4031-prolonged action potential duration and QT interval. INTERPRETATION: We consider the endocannabinoids as an interesting class of hKV7.1/KCNE1 channel modulators with putative protective effects in LQTS contexts. FUNDING: ERC (No. 850622), Canadian Institutes of Health Research, Canada Research Chairs and Compute Canada, Swedish National Infrastructure for Computing.

CryoEM structures of anion exchanger 1 capture multiple states of inward- and outward-facing conformations.

Anion exchanger 1 (AE1, band 3) is a major membrane protein of red blood cells and plays a key role in acid-base homeostasis, urine acidification, red blood cell shape regulation, and removal of carbon dioxide during respiration. Though structures of the transmembrane domain (TMD) of three SLC4 transporters, including AE1, have been resolved previously in their outward-facing (OF) state, no mammalian SLC4 structure has been reported in the inward-facing (IF) conformation. Here we present the cryoEM structures of full-length bovine AE1 with its TMD captured in both IF and OF conformations. Remarkably, both IF-IF homodimers and IF-OF heterodimers were detected. The IF structures feature downward movement in the core domain with significant unexpected elongation of TM11. Molecular modeling and structure guided mutagenesis confirmed the functional significance of residues involved in TM11 elongation. Our data provide direct evidence for an elevator-like mechanism of ion transport by an SLC4 family member.

The Single Residue K12 Governs the Exceptional Voltage Sensitivity of Mitochondrial Voltage-Dependent Anion Channel Gating.

The voltage-dependent anion channel (VDAC) is a ß-barrel channel of the mitochondrial outer membrane (MOM) that passively transports ions, metabolites, polypeptides, and single-stranded DNA. VDAC responds to a transmembrane potential by "gating," i.e. transitioning to one of a variety of low-conducting states of unknown structure. The gated state results in nearly complete suppression of multivalent mitochondrial metabolite (such as ATP and ADP) transport, while enhancing calcium transport. Voltage gating is a universal property of ß-barrel channels, but VDAC gating is anomalously sensitive to transmembrane potential. Here, we show that a single residue in the pore interior, K12, is responsible for most of VDAC's voltage sensitivity. Using the analysis of over 40 µs of atomistic molecular dynamics (MD) simulations, we explore correlations between motions of charged residues inside the VDAC pore and geometric deformations of the ß-barrel. Residue K12 is bistable; its motions between two widely separated positions along the pore axis enhance the fluctuations of the ß-barrel and augment the likelihood of gating. Single channel electrophysiology of various K12 mutants reveals a dramatic reduction of the voltage-induced gating transitions. The crystal structure of the K12E mutant at a resolution of 2.6 Å indicates a similar architecture of the K12E mutant to the wild type; however, 60 µs of atomistic MD simulations using the K12E mutant show restricted motion of residue 12, due to enhanced connectivity with neighboring residues, and diminished amplitude of barrel motions. We conclude that ß-barrel fluctuations, governed particularly by residue K12, drive VDAC gating transitions.

A molecular switch controls the impact of cholesterol on a Kir channel.

SignificanceCholesterol is one of the main components found in plasma membranes and is involved in lipid-dependent signaling enabled by integral membrane proteins such as inwardly rectifying potassium (Kir) channels. Similar to other ion channels, most of the Kir channels are down-regulated by cholesterol. One of the very few notable exceptions is Kir3.4, which is up-regulated by this important lipid. Here, we discovered and characterized a molecular switch that controls the impact (up-regulation vs. down-regulation) of cholesterol on Kir3.4. Our results provide a detailed molecular mechanism of tunable cholesterol regulation of a potassium channel.

The voltage-sensing domain of a hERG1 mutant is a cation-selective channel.

A cationic leak current known as an "omega current" may arise from mutations of the first charged residue in the S4 of the voltage sensor domains of sodium and potassium voltage-gated channels. The voltage-sensing domains (VSDs) in these mutated channels act as pores allowing nonspecific passage of cations, such as Li+, K+, Cs+, and guanidinium. Interestingly, no omega currents have been previously detected in the nonswapped voltage-gated potassium channels such as the human-ether-a-go-go-related (hERG1), hyperpolarization-activated cyclic nucleotide-gated, and ether-a-go-go channels. In this work, we discovered a novel omega current by mutating the first charged residue of the S4 of the hERG1, K525 to serine. To characterize this omega current, we used various probes, including the hERG1 pore domain blocker, dofetilide, to show that the omega current does not require cation flux via the canonical pore domain. In addition, the omega flux does not cross the conventional selectivity filter. We also show that the mutated channel (K525S hERG1) conducts guanidinium. These data are indicative of the formation of an omega current channel within the VSD. Using molecular dynamics simulations with replica-exchange umbrella sampling simulations of the wild-type hERG1 and the K525S hERG1, we explored the molecular underpinnings governing the cation flow in the VSD of the mutant. We also show that the wild-type hERG1 may form water crevices supported by the biophysical surface accessibility data. Overall, our multidisciplinary study demonstrates that the VSD of hERG1 may act as a cation-selective channel wherein a mutation of the first charged residue in the S4 generates an omega current. Our simulation uncovers the atomistic underpinning of this mechanism.

A second S4 movement opens hyperpolarization-activated HCN channels.

Rhythmic activity in pacemaker cells, as in the sino-atrial node in the heart, depends on the activation of hyperpolarization-activated cyclic nucleotide-gated (HCN) channels. As in depolarization-activated K+ channels, the fourth transmembrane segment S4 functions as the voltage sensor in hyperpolarization-activated HCN channels. But how the inward movement of S4 in HCN channels at hyperpolarized voltages couples to channel opening is not understood. Using voltage clamp fluorometry, we found here that S4 in HCN channels moves in two steps in response to hyperpolarizations and that the second S4 step correlates with gate opening. We found a mutation in sea urchin HCN channels that separate the two S4 steps in voltage dependence. The E356A mutation in S4 shifts the main S4 movement to positive voltages, but channel opening remains at negative voltages. In addition, E356A reveals a second S4 movement at negative voltages that correlates with gate opening. Cysteine accessibility and molecular models suggest that the second S4 movement opens up an intracellular crevice between S4 and S5 that would allow radial movement of the intracellular ends of S5 and S6 to open HCN channels.

Cryo-EM structure of the sodium-driven chloride/bicarbonate exchanger NDCBE.

SLC4 transporters play significant roles in pH regulation and cellular sodium transport. The previously solved structures of the outward facing (OF) conformation for AE1 (SLC4A1) and NBCe1 (SLC4A4) transporters revealed an identical overall fold despite their different transport modes (chloride/bicarbonate exchange versus sodium-carbonate cotransport). However, the exact mechanism determining the different transport modes in the SLC4 family remains unknown. In this work, we report the cryo-EM 3.4 Å structure of the OF conformation of NDCBE (SLC4A8), which shares transport properties with both AE1 and NBCe1 by mediating the electroneutral exchange of sodium-carbonate with chloride. This structure features a fully resolved extracellular loop 3 and well-defined densities corresponding to sodium and carbonate ions in the tentative substrate binding pocket. Further, we combine computational modeling with functional studies to unravel the molecular determinants involved in NDCBE and SLC4 transport.

Toward Reducing hERG Affinities for DAT Inhibitors with a Combined Machine Learning and Molecular Modeling Approach.

Psychostimulant drugs, such as cocaine, inhibit dopamine reuptake via blockading the dopamine transporter (DAT), which is the primary mechanism underpinning their abuse. Atypical DAT inhibitors are dissimilar to cocaine and can block cocaine- or methamphetamine-induced behaviors, supporting their development as part of a treatment regimen for psychostimulant use disorders. When developing these atypical DAT inhibitors as medications, it is necessary to avoid off-target binding that can produce unwanted side effects or toxicities. In particular, the blockade of a potassium channel, human ether-a-go-go (hERG), can lead to potentially lethal ventricular tachycardia. In this study, we established a counter screening platform for DAT and against hERG binding by combining machine learning-based quantitative structure-activity relationship (QSAR) modeling, experimental validation, and molecular modeling and simulations. Our results show that the available data are adequate to establish robust QSAR models, as validated by chemical synthesis and pharmacological evaluation of a validation set of DAT inhibitors. Furthermore, the QSAR models based on subsets of the data according to experimental approaches used have predictive power as well, which opens the door to target specific functional states of a protein. Complementarily, our molecular modeling and simulations identified the structural elements responsible for a pair of DAT inhibitors having opposite binding affinity trends at DAT and hERG, which can be leveraged for rational optimization of lead atypical DAT inhibitors with desired pharmacological properties.

The origin and radiation of the phosphoprotein phosphatase (PPP) enzymes of Eukaryotes.

Phosphoprotein phosphatase (PPP) enzymes are ubiquitous proteins involved in cellular signaling pathways and other functions. Here we have traced the origin of the PPP sequences of Eukaryotes and their radiation. Using a bacterial PPP Hidden Markov Model (HMM) we uncovered "BacterialPPP-Like" sequences in Archaea. A HMM derived from eukaryotic PPP enzymes revealed additional, unique sequences in Archaea and Bacteria that were more like the eukaryotic PPP enzymes then the bacterial PPPs. These sequences formed the basis of phylogenetic tree inference and sequence structural analysis allowing the history of these sequence types to be elucidated. Our phylogenetic tree data strongly suggest that eukaryotic PPPs ultimately arose from ancestors in the Asgard archaea. We have clarified the radiation of PPPs within Eukaryotes, substantially expanding the range of known organisms with PPP subtypes (Bsu1, PP7, PPEF/RdgC) previously thought to have a more restricted distribution. Surprisingly, sequences from the Methanosarcinaceae (Euryarchaeota) form a strongly supported sister group to eukaryotic PPPs in our phylogenetic analysis. This strongly suggests an intimate association between an Asgard ancestor and that of the Methanosarcinaceae. This is highly reminiscent of the syntrophic association recently demonstrated between the cultured Lokiarchaeal species Prometheoarchaeum and a methanogenic bacterial species.

A general mechanism of KCNE1 modulation of KCNQ1 channels involving non-canonical VSD-PD coupling.

Voltage-gated KCNQ1 channels contain four separate voltage-sensing domains (VSDs) and a pore domain (PD). KCNQ1 expressed alone opens when the VSDs are in an intermediate state. In cardiomyocytes, KCNQ1 co-expressed with KCNE1 opens mainly when the VSDs are in a fully activated state. KCNE1 also drastically slows the opening of KCNQ1 channels and shifts the voltage dependence of opening by >40 mV. We here show that mutations of conserved residues at the VSD-PD interface alter the VSD-PD coupling so that the mutant KCNQ1/KCNE1 channels open in the intermediate VSD state. Using recent structures of KCNQ1 and KCNE beta subunits in different states, we present a mechanism by which KCNE1 rotates the VSD relative to the PD and affects the VSD-PD coupling of KCNQ1 channels in a non-canonical way, forcing KCNQ1/KCNE1 channels to open in the fully-activated VSD state. This would explain many of the KCNE1-induced effects on KCNQ1 channels.

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.

Identification of PUFA interaction sites on the cardiac potassium channel KCNQ1.

Polyunsaturated fatty acids (PUFAs), but not saturated fatty acids, modulate ion channels such as the cardiac KCNQ1 channel, although the mechanism is not completely understood. Using both simulations and experiments, we find that PUFAs interact directly with the KCNQ1 channel via two different binding sites: one at the voltage sensor and one at the pore. These two amphiphilic binding pockets stabilize the negatively charged PUFA head group by electrostatic interactions with R218, R221, and K316, while the hydrophobic PUFA tail is selectively stabilized by cassettes of hydrophobic residues. The rigid saturated tail of stearic acid prevents close contacts with KCNQ1. By contrast, the mobile tail of PUFA linoleic acid can be accommodated in the crevice of the hydrophobic cassette, a defining feature of PUFA selectivity in KCNQ1. In addition, we identify Y268 as a critical PUFA anchor point underlying fatty acid selectivity. Combined, this study provides molecular models of direct interactions between PUFAs and KCNQ1 and identifies selectivity mechanisms. Long term, this understanding may open new avenues for drug development based on PUFA mechanisms.

Identification of multiple substrate binding sites in SLC4 transporters in the outward-facing conformation: Insights into the transport mechanism.

Solute carrier family 4 (SLC4) transporters mediate the transmembrane transport of HCO3-, CO32-, and Cl- necessary for pH regulation, transepithelial H+/base transport, and ion homeostasis. Substrate transport with varying stoichiometry and specificity is achieved through an exchange mechanism and/or through coupling of the uptake of anionic substrates to typically co-transported Na+. Recently solved outward-facing structures of two SLC4 members (human anion exchanger 1 [hAE1] and human electrogenic sodium bicarbonate cotransporter 1 [hNBCe1]) with different transport modes (Cl-/HCO3- exchange versus Na+-CO32- symport) revealed highly conserved three-dimensional organization of their transmembrane domains. However, the exact location of the ion binding sites and their protein-ion coordination motifs are still unclear. In the present work, we combined site identification by ligand competitive saturation mapping and extensive molecular dynamics sampling with functional mutagenesis studies which led to the identification of two substrate binding sites (entry and central) in the outward-facing states of hAE1 and hNBCe1. Mutation of residues in the identified binding sites led to impaired transport in both proteins. We also showed that R730 in hAE1 is crucial for anion binding in both entry and central sites, whereas in hNBCe1, a Na+ acts as an anchor for CO32- binding to the central site. Additionally, protonation of the central acidic residues (E681 in hAE1 and D754 in hNBCe1) alters the ion dynamics in the permeation cavity and may contribute to the transport mode differences in SLC4 proteins. These results provide a basis for understanding the functional differences between hAE1 and hNBCe1 and may facilitate potential drug development for diseases such as proximal and distal renal tubular acidosis.

Lipid regulation of hERG1 channel function.

The lipid regulation of mammalian ion channel function has emerged as a fundamental mechanism in the control of electrical signalling and transport specificity in various cell types. In this work, we combine molecular dynamics simulations, mutagenesis, and electrophysiology to provide mechanistic insights into how lipophilic molecules (ceramide-sphingolipid probe) alter gating kinetics and K+ currents of hERG1. We show that the sphingolipid probe induced a significant left shift of activation voltage, faster deactivation rates, and current blockade comparable to traditional hERG1 blockers. Microseconds-long MD simulations followed by experimental mutagenesis elucidated ceramide specific binding locations at the interface between the pore and voltage sensing domains. This region constitutes a unique crevice present in mammalian channels with a non-swapped topology. The combined experimental and simulation data provide evidence for ceramide-induced allosteric modulation of the channel by a conformational selection mechanism.

Polarization Effects in Water-Mediated Selective Cation Transport across a Narrow Transmembrane Channel.

Despite the progress in modeling complex molecular systems of ever-increasing complexity, a quantitatively accurate computational treatment of ion permeation through narrow membrane channels remains challenging. An important factor to reach this goal is induced electronic polarization, which is likely to impact the permeation rate of small ions through narrow molecular pores. In this work, we extended the recently developed polarizable force field based on the classical Drude oscillators to assess the role of induced polarization effects on the energetics of sodium and potassium ion transport across the gramicidin A (gA) ion channel. The inclusion of induced polarization lowers barriers present in 1D potential of mean force (PMF) for cation permeation by ∼50% compared to those obtained with the additive force field. Conductance properties calculated with 1D PMFs from Drude simulations are in better agreement with experimental results. Polarization of single-file water molecules and protein atoms forming the narrow pore has a direct impact on the free-energy barriers and cation-specific solid-state NMR chemical shifts. Sensitivity analysis indicates that small changes to water-channel interactions can alter the free energy barrier for ion permeation. These results, illustrating polarization effects present in the complex electrostatic environment of the gA channel, have broad implications for revising proposed mechanisms of ion permeation and selectivity in a variety of ion channels.

Refinement of a cryo-EM structure of hERG: Bridging structure and function.

The human-ether-a-go-go-related gene (hERG) encodes the voltage-gated potassium channel (KCNH2 or Kv11.1, commonly known as hERG). This channel plays a pivotal role in the stability of phase 3 repolarization of the cardiac action potential. Although a high-resolution cryo-EM structure is available for its depolarized (open) state, the structure surprisingly did not feature many functionally important interactions established by previous biochemical and electrophysiology experiments. Using molecular dynamics flexible fitting (MDFF), we refined the structure and recovered the missing functionally relevant salt bridges in hERG in its depolarized state. We also performed electrophysiology experiments to confirm the functional relevance of a novel salt bridge predicted by our refinement protocol. Our work shows how refinement of a high-resolution cryo-EM structure helps to bridge the existing gap between the structure and function in the voltage-sensing domain (VSD) of hERG.