Fluorescence Imaging of Neural Activity, Neurochemical Dynamics, and Drug-Specific Receptor Conformation with Genetically Encoded Sensors.

Recent advances in fluorescence imaging permit large-scale recording of neural activity and dynamics of neurochemical release with unprecedented resolution in behaving animals. Calcium imaging with highly optimized genetically encoded indicators provides a mesoscopic view of neural activity from genetically defined populations at cellular and subcellular resolutions. Rigorously improved voltage sensors and microscopy allow for robust spike imaging of populational neurons in various brain regions. In addition, recent protein engineering efforts in the past few years have led to the development of sensors for neurotransmitters and neuromodulators. Here, we discuss the development and applications of these genetically encoded fluorescent indicators in reporting neural activity in response to various behaviors in different biological systems as well as in drug discovery. We also report a simple model to guide sensor selection and optimization.

A high-affinity cocaine binding site associated with the brain acid soluble protein 1.

Cocaine exerts its stimulant effect by inhibiting dopamine (DA) reuptake, leading to increased dopamine signaling. This action is thought to reflect the binding of cocaine to the dopamine transporter (DAT) to inhibit its function. However, cocaine is a relatively weak inhibitor of DAT, and many DAT inhibitors do not share cocaine's behavioral actions. Further, recent reports show more potent actions of the drug, implying the existence of a high-affinity receptor for cocaine. We now report high-affinity binding of cocaine associated with the brain acid soluble protein 1 (BASP1) with a dissociation constant (Kd) of 7 nM. Knocking down BASP1 in the striatum inhibits [3H]cocaine binding to striatal synaptosomes. Depleting BASP1 in the nucleus accumbens but not the dorsal striatum diminishes locomotor stimulation in mice. Our findings imply that BASP1 is a pharmacologically relevant receptor for cocaine.

Functional characterization of inactivating ABCC8 variants causing congenital hyperinsulinism.

Congenital hyperinsulinism (CHI; OMIM: 256450) is characterized by persistent insulin secretion despite severe hypoglycemia. The most common causes are variants in the ATP-binding cassette subfamily C member 8(ABCC8) and potassium inwardly-rectifying channel subfamily J member 11(KCNJ11) genes. These encode ATP-sensitive potassium (KATP) channel subunit sulfonylurea receptor 1 (SUR1) and inwardly rectifying potassium channel (Kir6.2) proteins. A 7-day-old male infant presented with frequent hypoglycemic episodes and was clinically diagnosed with CHI, underwent trio-whole-exome sequencing, revealing compound heterozygous ABCC8 variants (c.307C>T, p.His103Tyr; and c.3313_3315del, p.Ile1105del) were identified. In human embryonic kidney 293 (HEK293) and rat insulinoma cells (INS-1) transfected with wild-type and variant plasmids, KATP channels formed by p.His103Tyr were delivered to the plasma membrane, whereas p.Ile1105del or double variants (p.His103Tyr coupled with p.Ile1105del) failed to be transported to the plasma membrane. Compared to wild-type channels, the channels formed by the variants (p.His103Tyr; p.Ile1105del) had elevated basal [Ca2+]i, but did not respond to stimulation by glucose. Our results provide evidence that the two ABCC8 variants may be related to CHI owing to defective trafficking and dysfunction of KATP channels.

Cryo-Electron Microscopy: Moving Beyond X-Ray Crystal Structures for Drug Receptors and Drug Development.

Electron cryo-microscopy (cryo-EM) has revolutionized structure determination of membrane proteins and holds great potential for structure-based drug discovery. Here we discuss the potential of cryo-EM in the rational design of therapeutics for membrane proteins compared to X-ray crystallography. We also detail recent progress in the field of drug receptors, focusing on cryo-EM of two protein families with established therapeutic value, the γ-aminobutyric acid A receptors (GABAARs) and G protein-coupled receptors (GPCRs). GABAARs are pentameric ion channels, and cryo-EM structures of physiological heteromeric receptors in a lipid environment have uncovered the molecular basis of receptor modulation by drugs such as diazepam. The structures of ten GPCR-G protein complexes from three different classes of GPCRs have now been determined by cryo-EM. These structures give detailed insights into molecular interactions with drugs, GPCR-G protein selectivity, how accessory membrane proteins alter receptor-ligand pharmacology, and the mechanism by which HIV uses GPCRs to enter host cells.

Structural Insights into ATP-Sensitive Potassium Channel Mechanics: A Role of Intrinsically Disordered Regions.

Commonly used techniques, such as CryoEM or X-ray, are not able to capture the structural reorganizations of disordered regions of proteins (IDR); therefore, it is difficult to assess their functions in proteins based exclusively on experiments. To fill this gap, we used computational molecular dynamics (MD) simulation methods to capture IDR dynamics and trace biological function-related interactions in the Kir6.2/SUR1 potassium channel. This ATP-sensitive octameric complex, one of the critical elements in the insulin secretion process in human pancreatic ß-cells, has four to five large, disordered fragments. Using unique MD simulations of the full Kir6.2/SUR1 channel complex, we present an in-depth analysis of the dynamics of the disordered regions and discuss the possible functions they could have in this system. Our MD results confirmed the crucial role of the N-terminus of the Kir6.2 fragment and the L0-loop of the SUR1 protein in the transfer of mechanical signals between domains that trigger insulin release. Moreover, we show that the presence of IDRs affects natural ligand binding. Our research takes us one step further toward understanding the action of this vital complex.

The Emerging Structural Pharmacology of ATP-Sensitive Potassium Channels.

ATP-sensitive potassium channels (KATP) are energy sensors that participate in a range of physiologic processes. These channels are also clinically validated drug targets. For decades, KATP inhibitors have been prescribed for diabetes and KATP activators have been used for the treatment of hypoglycemia, hypertension, and hair loss. In this Emerging Concepts article, we highlight our current knowledge about the drug binding modes observed using cryogenic electron microscopy techniques. The inhibitors and activators bind to two distinct sites in the transmembrane domain of the sulfonylurea receptor (SUR) subunit. We also discuss the possible mechanism of how these drugs allosterically modulate the dimerization of SUR nucleotide-binding domains (NBDs) and thus KATP channel activity. SIGNIFICANCE STATEMENT: ATP-sensitive potassium channels (KATP) are fundamental to energy homeostasis, and they participate in many vital physiological processes. KATP channels are important drug targets. Both KATP inhibitors (insulin secretagogues) and KATP activators are broadly used clinically for the treatment of related diseases. Recent cryogenic electron microscopy studies allow us to understand the emerging concept of KATP structural pharmacology.

WNK4 kinase: from structure to physiology.

With no lysine kinase-4 (WNK4) belongs to a serine-threonine kinase family characterized by the atypical positioning of its catalytic lysine. Despite the fact that WNK4 has been found in many tissues, the majority of its study has revolved around its function in the kidney, specifically as a positive regulator of the thiazide-sensitive NaCl cotransporter (NCC) in the distal convoluted tubule of the nephron. This is explained by the description of gain-of-function mutations in the gene encoding WNK4 that causes familial hyperkalemic hypertension. This disease is mainly driven by increased downstream activation of the Ste20/SPS1-related proline-alanine-rich kinase/oxidative stress responsive kinase-1-NCC pathway, which increases salt reabsorption in the distal convoluted tubule and indirectly impairs renal K+ secretion. Here, we review the large volume of information that has accumulated about different aspects of WNK4 function. We first review the knowledge on WNK4 structure and enumerate the functional domains and motifs that have been characterized. Then, we discuss WNK4 physiological functions based on the information obtained from in vitro studies and from a diverse set of genetically modified mouse models with altered WNK4 function. We then review in vitro and in vivo evidence on the different levels of regulation of WNK4. Finally, we go through the evidence that has suggested how different physiological conditions act through WNK4 to modulate NCC activity.

SAMPL7 blind predictions using nonequilibrium alchemical approaches.

In the context of the SAMPL7 challenge, we computed, employing a non-equilibrium (NE) alchemical technique, the standard binding free energy of two series of host-guest systems, involving as a host the Isaac's TrimerTrip, a Cucurbituril-like open cavitand, and the Gilson's Cyclodextrin derivatives. The adopted NE alchemy combines enhanced sampling molecular dynamics simulations with driven fast out-of-equilibrium alchemical trajectories to recover the free energy via the Jarzynski and Crooks NE theorems. The GAFF2 non-polarizable force field was used for the parametrization. Performances were acceptable and similar in accuracy to those we submitted for Gibb's Deep Cavity Cavitands in the previous SAMPL6 host-guest challenge, confirming the reliability of the computational approach and exposing, in some cases, some important deficiencies of the GAFF2 non-polarizable force field.

Estimating drug potency in the competitive target mediated drug disposition (TMDD) system when the endogenous ligand is included.

Predictions for target engagement are often used to guide drug development. In particular, when selecting the recommended phase 2 dose of a drug that is very safe, and where good biomarkers for response may not exist (e.g. in immuno-oncology), a receptor occupancy prediction could even be the main determinant in justifying the approved dose, as was the case for atezolizumab. The underlying assumption in these models is that when the drug binds its target, it disrupts the interaction between the target and its endogenous ligand, thereby disrupting downstream signaling. However, the interaction between the target and its endogenous binding partner is almost never included in the model. In this work, we take a deeper look at the in vivo system where a drug binds to its target and disrupts the target's interaction with an endogenous ligand. We derive two simple steady state inhibition metrics (SSIMs) for the system, which provides intuition for when the competition between drug and endogenous ligand should be taken into account for guiding drug development.

The Effect of DREADD Activation of Leptin Receptor Positive Neurons in the Nucleus of the Solitary Tract on Sleep Disordered Breathing.

Obstructive sleep apnea (OSA) is recurrent obstruction of the upper airway due to the loss of upper airway muscle tone during sleep. OSA is highly prevalent, especially in obesity. There is no pharmacotherapy for OSA. Previous studies have demonstrated the role of leptin, an adipose-tissue-produced hormone, as a potent respiratory stimulant. Leptin signaling via a long functional isoform of leptin receptor, LEPRb, in the nucleus of the solitary tract (NTS), has been implicated in control of breathing. We hypothesized that leptin acts on LEPRb positive neurons in the NTS to increase ventilation and maintain upper airway patency during sleep in obese mice. We expressed designer receptors exclusively activated by designer drugs (DREADD) selectively in the LEPRb positive neurons of the NTS of Leprb-Cre-GFP mice with diet-induced obesity (DIO) and examined the effect of DREADD ligand, J60, on tongue muscle activity and breathing during sleep. J60 was a potent activator of LEPRb positive NTS neurons, but did not stimulate breathing or upper airway muscles during NREM and REM sleep. We conclude that, in DIO mice, the stimulating effects of leptin on breathing during sleep are independent of LEPRb signaling in the NTS.

Receptor-Based Pharmacophore Modeling in the Search for Natural Products for COVID-19 Mpro.

Considering the urgency of the COVID-19 pandemic, we developed a receptor-based pharmacophore model for identifying FDA-approved drugs and hits from natural products. The COVID-19 main protease (Mpro) was selected for the development of the pharmacophore model. The model consisted of a hydrogen bond acceptor, donor, and hydrophobic features. These features demonstrated good corroboration with a previously reported model that was used to validate the present model, showing an RMSD value of 0.32. The virtual screening was carried out using the ZINC database. A set of 208,000 hits was extracted and filtered using the ligand pharmacophore mapping, applying the lead-like properties. Lipinski's filter and the fit value filter were used to minimize hits to the top 2000. Simultaneous docking was carried out for 200 hits for natural drugs belonging to the FDA-approved drug database. The top 28 hits from these experiments, with promising predicted pharmacodynamic and pharmacokinetic properties, are reported here. To optimize these hits as Mpro inhibitors and potential treatment options for COVID-19, bench work investigations are needed.

PGF2α stimulates the 10-pS Cl- channel and thiazide-sensitive Na+-Cl- cotransporter in the distal convoluted tubule.

We used patch-clamp and Western blot analysis to test whether PGF2α stimulates the basolateral 10-pS Cl- channel and thiazide-sensitive Na+-Cl- cotransporter (NCC) in the distal convoluted tubule (DCT) via a prostaglandin F receptor (FP-R). Single channel and whole cell recordings demonstrated that PGF2α stimulated the 10-pS Cl- channel in the DCT. The stimulatory effect of PGF2α on the Cl- channel was mimicked by a FP-R agonist, latanoprost, but was abrogated by blocking FP-R with AL8810. Also, the effect of PGF2α on the Cl- channel in the DCT was recapitulated by stimulating PKC but was blocked by inhibiting PKC. Furthermore, inhibition of p38 MAPK but not ERK blocked the effect of PGF2α on the 10-pS Cl- channel. Inhibition of NADPH oxidase also abrogated the stimulatory effect of PGF2α on the 10-pS Cl- channel, while the addition of 10 µM H2O2 mimicked the stimulatory effect of PGF2α on the 10-pS Cl- channel. Moreover, superoxide-related species may mediate the stimulatory effect of PGF2α on the 10-pS Cl- channel because the stimulatory effect of PGF2α and H2O2 was not additive. Western blot analysis showed that infusion of PGF2α in vivo not only increased the expression of FP-R but also increased the expression of total NCC and phosphorylated NCC. We conclude that PGF2α stimulates the basolateral 10-pS Cl- channel in the DCT by activating FP-R through PKC/p38 MAPK and NADPH oxidase-dependent pathways. The stimulatory effects of PGF2α on the Cl- channel and NCC may contribute to PGF2α-induced increases in NaCl reabsorption in the DCT.

Determination of Real Time in Vivo Drug Receptor Occupancy for a Covalent Binding Drug as a Clinical Pharmacodynamic Biomarker by Immunocapture-LC-MS/MS.

We report a novel immunocapture (IC)-LC-MS/MS methodology to directly measure real time in vivo receptor occupancy (RO) for a covalent binding drug in blood lysate. A small molecule quencher was added immediately after sample collection to convert the free receptor to a quencher-bound receptor (QB-R) which was measured with the drug-bound receptor (DB-R) simultaneously by LC-MS/MS after immunocapture enrichment, followed by trypsin digestion. Addition of the quencher is necessary to prevent the free receptor from ex vivo binding with the drug. The real time RO was calculated based on the concentrations of DB-R and the free receptor (which is now QB-R) that were obtained from each sample. This strategy has been successfully applied to the measurement of the RO for Bruton's tyrosine kinase (BTK) in the blood lysate of monkeys after dosing with branebrutinib (BMS-986195), a covalent BTK inhibitor being evaluated to treat rheumatoid arthritis. A custom-made quencher, which is more reactive to BTK than branebrutinib, was added in excess amount to bind with all available free BTK to form quencher-bound BTK (QB-BTK) during blood sample collection. To measure a wide range of % BTK RO, including those of <5% or >95%, the required LLOQ at 0.125 nM for QB-BTK and 0.250 nM for drug-bound BTK (DB-BTK) in blood lysate were successfully achieved by using this IC-LC-MS/MS strategy. This proof-of-concept assay demonstrated its suitability with high throughput for real time in vivo BTK RO measurement as a pharmacodynamic (PD) biomarker for clinical drug development.

Binding energies of the drugs capreomycin and streptomycin in complex with tuberculosis bacterial ribosome subunits.

Despite advances, tuberculosis remains a significant infectious disease, whose mortality presents alarming numbers. Although it can be cured, the number of cases of antimicrobial resistant strains is increasing, requiring the use of less efficient second-line drugs. Capreomycin and streptomycin are part of this group, being antibiotics whose mechanism of action is the inhibition of protein synthesis when interacting with the tuberculosis bacterial ribosome. Their binding mechanisms are distinct: capreomycin is able to bind to both ribosomal (30S and 50S) subunits, whereas streptomycin binds only to the smaller one (30S). In this context, the biochemical characterization of these binding sites for a proper understanding of their complex interactions is of crucial importance to increase their efficacy. Through crystallographic data and computer simulations, in this work we calculated the interaction binding energies of capreomycin and streptomycin in complex with the tuberculosis bacterial ribosome subunits, by using density functional theory (DFT) within the molecular fractionation with conjugated caps (MFCC) approach. For capreomycin in the 30S (50S) subunit, we investigated the binding energies of 44 (30) residues presented within a pocket radius of 14 Å (30 Å). Regarding streptomycin, 60 nucleotide (25 amino acid) residues distributed up to 12.5 Å (15 Å) away from the drug in the 30S subunit (S12 protein) were taken into account. We also identify the contributions of hydrogen bonds and hydrophobic interactions in the drug-receptor complex, and the regions of the drugs that most contributed to the anchorages of them in their binding sites, as well as identify residues that are most associated with mutations.

Quantitative PET Imaging in Drug Development: Estimation of Target Occupancy.

Positron emission tomography, an imaging tool using radiolabeled tracers in humans and preclinical species, has been widely used in recent years in drug development, particularly in the central nervous system. One important goal of PET in drug development is assessing the occupancy of various molecular targets (e.g., receptors, transporters, enzymes) by exogenous drugs. The current linear mathematical approaches used to determine occupancy using PET imaging experiments are presented. These algorithms use results from multiple regions with different target content in two scans, a baseline (pre-drug) scan and a post-drug scan. New mathematical estimation approaches to determine target occupancy, using maximum likelihood, are presented. A major challenge in these methods is the proper definition of the covariance matrix of the regional binding measures, accounting for different variance of the individual regional measures and their nonzero covariance, factors that have been ignored by conventional methods. The novel methods are compared to standard methods using simulation and real human occupancy data. The simulation data showed the expected reduction in variance and bias using the proper maximum likelihood methods, when the assumptions of the estimation method matched those in simulation. Between-method differences for data from human occupancy studies were less obvious, in part due to small dataset sizes. These maximum likelihood methods form the basis for development of improved PET covariance models, in order to minimize bias and variance in PET occupancy studies.

WNK signalling pathways in blood pressure regulation.

Hypertension (high blood pressure) is a major public health problem affecting more than a billion people worldwide with complications, including stroke, heart failure and kidney failure. The regulation of blood pressure is multifactorial reflecting genetic susceptibility, in utero environment and external factors such as obesity and salt intake. In keeping with Arthur Guyton's hypothesis, the kidney plays a key role in blood pressure control and data from clinical studies; physiology and genetics have shown that hypertension is driven a failure of the kidney to excrete excess salt at normal levels of blood pressure. There is a number of rare Mendelian blood pressure syndromes, which have shed light on the molecular mechanisms involved in dysregulated ion transport in the distal kidney. One in particular is Familial hyperkalemic hypertension (FHHt), an autosomal dominant monogenic form of hypertension characterised by high blood pressure, hyperkalemia, hyperchloremic metabolic acidosis, and hypercalciuria. The clinical signs of FHHt are treated by low doses of thiazide diuretic, and it mirrors Gitelman syndrome which features the inverse phenotype of hypotension, hypokalemic metabolic alkalosis, and hypocalciuria. Gitelman syndrome is caused by loss of function mutations in the thiazide-sensitive Na/Cl cotransporter (NCC); however, FHHt patients do not have mutations in the SCL12A3 locus encoding NCC. Instead, mutations have been identified in genes that have revealed a key signalling pathway that regulates NCC and several other key transporters and ion channels in the kidney that are critical for BP regulation. This is the WNK kinase signalling pathway that is the subject of this review.

Asymptotic analysis of a TMDD model: when a reaction contributes to the destruction of its product.

The multi-scale dynamics of a two-compartment with first order absorption Target-Mediated Drug Disposition (TMDD) pharmacokinetics model is analysed, using the Computational Singular Perturbation (CSP) algorithm. It is shown that the process evolves along two Slow Invariant Manifolds (SIMs), on which the most intense components of the model are equilibrated, so that the less intensive are the driving ones. The CSP tools allow for the identification of the components of the TMDD model that (i) constrain the evolution of the process on the SIMs, (ii) drive the system along the SIMs and (iii) generate the fast time scales. Among others, such diagnostics identify (i) the factors that determine the start and the duration of the period in which the ligand-receptor complex acts and (ii) the processes that determine its degradation rate. The counterintuitive influence of the process that transfers the ligand from the tissue to the main compartment, as it is manifested during the final stage of the process, is studied in detail.

Disrupted Prediction Error Links Excessive Amygdala Activation to Excessive Fear.

Basolateral amygdala (BLA) is critical for fear learning, and its heightened activation is widely thought to underpin a variety of anxiety disorders. Here we used chemogenetic techniques in rats to study the consequences of heightened BLA activation for fear learning and memory, and to specifically identify a mechanism linking increased activity of BLA glutamatergic neurons to aberrant fear. We expressed the excitatory hM3Dq DREADD in rat BLA glutamatergic neurons and showed that CNO acted selectively to increase their activity, depolarizing these neurons and increasing their firing rates. This chemogenetic excitation of BLA glutamatergic neurons had no effect on the acquisition of simple fear learning, regardless of whether this learning led to a weak or strong fear memory. However, in an associative blocking task, chemogenetic excitation of BLA glutamatergic neurons yielded significant learning to a blocked conditioned stimulus, which otherwise should not have been learned about. Moreover, in an overexpectation task, chemogenetic manipulation of BLA glutamatergic neurons prevented use of negative prediction error to reduce fear learning, leading to significant impairments in fear inhibition. These effects were not attributable to the chemogenetic manipulation enhancing arousal, increasing asymptotic levels of fear learning or fear memory consolidation. Instead, chemogenetic excitation of BLA glutamatergic neurons disrupted use of prediction error to regulate fear learning. SIGNIFICANCE STATEMENT: Several neuropsychiatric disorders are characterized by heightened activation of the amygdala. This heightened activation has been hypothesized to underlie increased emotional reactivity, fear over generalization, and deficits in fear inhibition. Yet the mechanisms linking heightened amygdala activation to heightened emotional learning are elusive. Here we combined chemogenetic excitation of rat basolateral amygdala glutamatergic neurons with a variety of behavioral approaches to show that, although simple fear learning is unaffected, the use of prediction error to regulate this learning is profoundly disrupted, leading to formation of inappropriate fear associations and impaired fear inhibition.

Muscle KATP channels: recent insights to energy sensing and myoprotection.

ATP-sensitive potassium (K(ATP)) channels are present in the surface and internal membranes of cardiac, skeletal, and smooth muscle cells and provide a unique feedback between muscle cell metabolism and electrical activity. In so doing, they can play an important role in the control of contractility, particularly when cellular energetics are compromised, protecting the tissue against calcium overload and fiber damage, but the cost of this protection may be enhanced arrhythmic activity. Generated as complexes of Kir6.1 or Kir6.2 pore-forming subunits with regulatory sulfonylurea receptor subunits, SUR1 or SUR2, the differential assembly of K(ATP) channels in different tissues gives rise to tissue-specific physiological and pharmacological regulation, and hence to the tissue-specific pharmacological control of contractility. The last 10 years have provided insights into the regulation and role of muscle K(ATP) channels, in large part driven by studies of mice in which the protein determinants of channel activity have been deleted or modified. As yet, few human diseases have been correlated with altered muscle K(ATP) activity, but genetically modified animals give important insights to likely pathological roles of aberrant channel activity in different muscle types.

Comparative transcriptomic analysis identifies evolutionarily conserved gene products in the vertebrate renal distal convoluted tubule.

Understanding the molecular basis of the complex regulatory networks controlling renal ion transports is of major physiological and clinical importance. In this study, we aimed to identify evolutionarily conserved critical players in the function of the renal distal convoluted tubule (DCT) by a comparative transcriptomic approach. We generated a transgenic zebrafish line with expression of the red fluorescent mCherry protein under the control of the zebrafish DCT-specific promoter of the thiazide-sensitive NaCl cotransporter (NCC). The mCherry expression was then used to isolate from the zebrafish mesonephric kidneys the distal late (DL) segments, the equivalent of the mammalian DCT, for subsequent RNA-seq analysis. We next compared this zebrafish DL transcriptome to the previously established mouse DCT transcriptome and identified a subset of gene products significantly enriched in both the teleost DL and the mammalian DCT, including SLCs and nuclear transcription factors. Surprisingly, several of the previously described regulators of NCC (e.g., SPAK, KLHL3, ppp1r1a) in the mouse were not found enriched in the zebrafish DL. Nevertheless, the zebrafish DL expressed enriched levels of related homologues. Functional knockdown of one of these genes, ppp1r1b, reduced the phosphorylation of NCC in the zebrafish pronephros, similar to what was seen previously in knockout mice for its homologue, Ppp1r1a. The present work is the first report on global gene expression profiling in a specific nephron portion of the zebrafish kidney, an increasingly used model system for kidney research. Our study suggests that comparative analysis of gene expression between phylogenetically distant species may be an effective approach to identify novel regulators of renal function.