Efficiency considerations in the construction of interpolated potential energy surfaces for the calculation of quantum observables by diffusion Monte Carlo

A modified Shepard interpolation scheme is used to construct global potential energy surfaces (PES) in order to calculate quantum observables--vibrationally averaged internal coordinates, fully anharmonic zero-point energies and nuclear radial distribution functions--for a prototypical loosely bound molecular system, the water dimer. The efficiency of PES construction is examined with respect to (a) the method used to sample configurational space, (b) the method used to choose which points to add to the PES data set, and (c) the use of either a one- or two-part weight function. The most efficient method for constructing the PES is found to require a quantum sampling regime, a combination of both h-weight and rms methods for choosing data points and use of the two-part weight function in the interpolation. Using this regime, the quantum diffusion Monte Carlo zero-point energy converges to the exact result within addition of 50 data points. The vibrationally averaged O-O distance and O-O radial distribution function, however, converge more slowly and require addition of over 500 data points. The methods presented here are expected to be applicable to both other loosely bound complexes as well as tightly bound molecular species. When combined with high quality ab initio calculations, these methods should be able to accurately characterize the PES of such species.D.L.C. would like to acknowledge the financial support of an Australian Postgraduate Research Award. This work has also been supported by Large Grant No. A00104447 from the Australian Research Council and by grants of computer time from the Australian Partnership in Advanced Computing (APAC) National Merit Allocation Scheme

Investigating the role of loop c hydrophilic residue 'T244' in the binding site of ρ1 GABAC receptors via site mutation and partial agonism

The loop C hydrophilic residue, threonine 244 lines the orthosteric binding site of ρ1 GABAC receptors was studied by point mutation into serine, alanine and cysteine, and tested with GABA, some representative partial agonists and antagonists. Thr244 has a hydroxyl group essential for GABA activity that is constrained by the threonine methyl group, orienting it toward the binding site. Significant decreases in activation effects of the studied ligands at ρ1 T244S mutant receptors, suggests a critical role for this residue. Results of aliphatic and heteroaromatic partial agonists demonstrate different pharmacological effects at ρ1 T244S mutant receptors when co-applied with GABA EC50 responses. ρ1 T244A and ρ1 T244C mutant receptors have minimal sensitivity to GABA at high mM concentrations, whereas, the ρ1 WT partial agonists, β-alanine and MTSEA demonstrate more efficacy and potency, respectively, than GABA at these mutant receptors. This study explores the role of Thr244 in the binding of agonists as an initial step during channel gating by moving loop C towards the ligand

Zolpidem is a potent stoichiometry-selective modulator of α1β3 GABAA receptors : evidence of a novel benzodiazepine site in the α1-α1 interface

Zolpidem is not a typical GABAA receptor hypnotic. Unlike benzodiazepines, zolpidem modulates tonic GABA currents in the rat dorsal motor nucleus of the vagus, exhibits residual effects in mice lacking the benzodiazepine binding site, and improves speech, cognitive and motor function in human patients with severe brain injury. The receptor by which zolpidem mediates these effects is not known. In this study we evaluated binary α1β3 GABAA receptors in either the 3α1:2β3 or 2α1:3β3 subunit stoichiometry, which differ by the existence of either an α1-α1 interface, or a β3-β3 interface, respectively. Both receptor stoichiometries are readily expressed in Xenopus oocytes, distinguished from each other by using GABA, zolpidem, diazepam and Zn2+. At the 3α1:2β3 receptor, clinically relevant concentrations of zolpidem enhanced GABA in a flumazenil-sensitive manner. The efficacy of diazepam was significantly lower compared to zolpidem. No modulation by either zolpidem or diazepam was detected at the 2α1:3β3 receptor, indicating that the binding site for zolpidem is at the α1-α1 interface, a site mimicking the classical α1-γ2 benzodiazepine site. Activating α1β3 (3α1:2β3) receptors may, in part, mediate the physiological effects of zolpidem observed under distinct physiological and clinical conditions, constituting a potentially attractive drug target

Cannabigerolic acid, a major biosynthetic precursor molecule in cannabis, exhibits divergent effects on seizures in mouse models of epilepsy

Background and Purpose: Cannabis has been used to treat epilepsy for millennia, with such use validated by regulatory approval of cannabidiol (CBD) for Dravet syndrome. Unregulated artisanal cannabis-based products used to treat children with intractable epilepsies often contain relatively low doses of CBD but are enriched in other phytocannabinoids. This raises the possibility that other cannabis constituents might have anticonvulsant properties. Experimental Approach: We used the Scn1a+/− mouse model of Dravet syndrome to investigate the cannabis plant for phytocannabinoids with anticonvulsant effects against hyperthermia-induced seizures. The most promising, cannabigerolic acid (CBGA), was further examined against spontaneous seizures and survival in Scn1a+/− mice and in electroshock seizure models. Pharmacological effects of CBGA were surveyed across multiple drug targets. Key Results: The initial screen identified three phytocannabinoids with novel anticonvulsant properties: CBGA, cannabidivarinic acid (CBDVA) and cannabigerovarinic acid (CBGVA). CBGA was most potent and potentiated the anticonvulsant effects of clobazam against hyperthermia-induced and spontaneous seizures, and was anticonvulsant in the MES threshold test. However, CBGA was proconvulsant in the 6-Hz threshold test and a high dose increased spontaneous seizure frequency in Scn1a+/− mice. CBGA was found to interact with numerous epilepsy-relevant targets including GPR55, TRPV1 channels and GABAA receptors. Conclusion and Implications: These results suggest that CBGA, CBDVA and CBGVA may contribute to the effects of cannabis-based products in childhood epilepsy. Although these phytocannabinoids have anticonvulsant potential and could be lead compounds for drug development programmes, several liabilities would need to be overcome before CBD is superseded by another in this class

Ligand-gated ion channels in genetic disorders and the question of efficacy

Whole-genome sequencing has unearthed a substantial number of individual variants in ion channels associated with genetic disorders. Ligand-gated ion channels, including glycine, γ-aminobutyric acid type A and nicotinic acetylcholine receptors, have long been known to harbour genetic variants associated with hyperekplexia and different forms of epilepsy. In some of these cases, missense variants enhance or impair the intrinsic ability of the receptor to convert ligand binding to channel opening, or the efficacy of receptor activation. We review the current understanding of how ligand-gated ion channels are activated and the properties that define the efficacy of an agonist, and how these properties can be altered by disease-causing variants. Additionally, we consider the mechanisms defining drug modulation of receptors and consider how this may differ in genetic variants. This fundamental knowledge is likely to be essential in understanding how effective treatments will be for patients with genetic variants in ligand-gated ion channels

High and low GABA sensitivity α4β2δ GABAA receptors are expressed in Xenopus laevis oocytes with divergent stoichiometries

GABAA receptors that contain the a4 and d subunits are thought to be located extrasynaptically, mediating tonic currents elicited by low concentrations of GABA. These a4bd receptors are modulated by neurosteroids and certain anesthetics, identifying them as important drug targets in research. However, pharmacological studies on these receptors have often yielded variable results, possibly due to the expression of receptors in different stoichiometries or arrangements. In this study, we injected different ratios of a4, b2 and d cRNA into Xenopus oocytes and measured the sensitivity to GABA and DS2 activation of the resulting receptor populations. By creating a matrix of RNA injection ratios from stock RNA concentrations, we were able to compare the changes in pharmacology between injection ratios where the ratio of only one subunit was altered. We identified two distinct populations of receptors, the first with an EC50 value of approximately 100 nM to GABA, a low Hill slope of approximately 0.3 and substantial direct activation by DS2. The second population had an EC50 value of approximately 1 lM to GABA, a steeper Hill slope of 1 and little direct activation, but substantial potentiation, by DS2. The second population was formed with high a4 ratios and low b2 ratios, but altering the ratio of d subunit injected had little effect. We propose that receptors with high sensitivity to GABA and direct activation by DS2 are the result of a greater number of b2 subunits being incorporated into the receptor

Coadministered cannabidiol and clobazam : preclinical evidence for both pharmacodynamic and pharmacokinetic interactions

Objective: Cannabidiol (CBD) has been approved by the US Food and Drug Administration (FDA) to treat intractable childhood epilepsies, such as Dravet syn-drome and Lennox‐Gastaut syndrome. However, the intrinsic anticonvulsant activity of CBD has been questioned due to a pharmacokinetic interaction between CBD and a first‐line medication, clobazam. This recognized interaction has led to speculation that the anticonvulsant efficacy of CBD may simply reflect CBD augmenting clobazam exposure. The present study aimed to address the nature of the interaction between CBD and clobazam. Methods: We examined whether CBD inhibits human CYP3A4 and CYP2C19 mediated metabolism of clobazam and N‐desmethylclobazam (N‐CLB), respectively, and performed studies assessing the effects of CBD on brain and plasma pharma-cokinetics of clobazam in mice. We then used the Scn1a+/− mouse model of Dravet syndrome to examine how CBD and clobazam interact. We compared anticonvulsant effects of CBD‐clobazam combination therapy to monotherapy against thermally‐in-duced seizures, spontaneous seizures and mortality in Scn1a+/− mice. In addition, we used Xenopus oocytes expressing γ‐aminobutyric acid (GABA)A receptors to investigate the activity of GABAA receptors when treated with CBD and clobazam together. Results: CBD potently inhibited CYP3A4 mediated metabolism of clobazam and CYP2C19 mediated metabolism of N‐CLB. Combination CBD‐clobazam treatment resulted in greater anticonvulsant efficacy in Scn1a+/− mice, but only when an anticonvulsant dose of CBD was used. It is important to note that a sub‐anticonvulsant dose of CBD did not promote greater anticonvulsant effects despite increasing plasma clobazam concentrations. In addition, we delineated a novel pharmacodynamic mechanism where CBD and clobazam together enhanced inhibitory GABAA receptor activation

A pharmacological characterization of GABA, THIP and DS2 at binary α4β3 and β38 receptors : GABA activates β38 receptors via the β3(+)8(-) interface

There is growing evidence that GABA (γ-aminobutyric acid) can activate GABAA receptors (GABAARs) in the absence of an α subunit. In this study, we compared the pharmacology of homomeric and binary α4, β3 or δ subunits with ternary α4β3δ to identify subunit interfaces that contribute to the pharmacology of GABA, THIP, and DS2, and the antagonists, Zn2þ, gabazine and bicuculline. β3δ receptors form functional GABA-gated channels when expressed in Xenopus oocytes with a pharmacology that differs to homomeric β3, binary α4β3 and ternary α4β3δ receptors. GABA had similar potency at α4β3 and β3δ receptors (25 mM and 26 mM, respectively) but differed at α4β3δ receptors where GABA exhibited a biphasic concentration-response (EC50 (1)¼12.6 nM; EC50 (2)¼6.3 μM). THIP activated β3δ receptors (EC50¼456 μM) but was a more potent activator of α4β3 (EC50¼27 μM) and α4β3δ receptors (EC50 (1)¼ 27.5 nM; EC50 (2)¼29.5 μΜ), indicating that the α4 subunit significantly contribute to its potency. The δpreferring modulator, DS2 had marginal or no effect at β3δ and α4β3 receptors, indicating a role for both the α4 and δ subunits for its potency. Gabazine inhibited GABA-elicited currents at β3δ receptors whereas bicuculline activated these receptors. Mutational analysis verified that GABA binds to the β3(þ) δ() interface formed by the β3 and δ subunits. In conclusion, evaluating agents against binary GABAARs such as β3δ and α4β3 receptors enables identification of interfaces that may contribute to the pharmacology of the more complex ternary α4β3δ receptors

AE Succinimide, an Analogue of Methyllycaconitine, When Bound Generates a Nonconducting Conformation of the α4β2 Nicotinic Acetylcholine Receptor

Nicotinic acetylcholine (nACh) receptors are pentameric ligand-gated ion channels that mediate fast synaptic transmission. The α4β2 nACh receptor is highly expressed in the brain and exists in two functional stoichiometries: the (α4)2(β2)3 and (α4)3(β2)2 that differ by an ACh-binding site at the α4−α4 interface of (α4)3(β2)2 receptors. Methyllycaconitine (MLA) is an nACh receptor antagonist, and while potent at both α7 and α4β2 nACh receptors, it has a higher selectivity for the α7 nACh receptor. The anthranilate-succinimide ester side-chain is important for its activity and selectivity. Here we identify a simplified MLA analogue that contains only the A and E ring skeleton of MLA, AE succinimide, that binds close to the channel lumen to display insurmountable inhibition at α4β2 nACh receptors. Although inhibition by AE succinimide was found to be voltage dependent indicating a possible pore channel blocker, substituted-cysteine accessibility experiments indicated it did not bind between 2′−16′ region of the channel pore. Instead, we found that upon binding and in the presence of ACh, there is a conformational change to the channel membrane that was identified when the compound was assessed against (α4 V13′C)β2 nACh receptors. It was found that in the 3:2 stoichiometry the two adjacent α4 subunits containing 13′ cysteine mutations formed a disulfide bond and occluded ion conductance. This was reversed by treatment with the reducing agent, dithiothreitol. Thus, AE succinimide has a different mechanism of inhibition to both MLA and other AE analogues, such as AE bicyclic alcohol, in that upon binding to an as yet unidentified site, AE succinimide in the presence of ACh induces a conformational change to the channel that generates a ligand-bound closed stateThis research was supported by a Project (APP1069417) from the Australian National Health and Medical Research Council (M.D.M. and M.C.) and by a Discovery Project (DP0986469) from the Australian Research Council (M.D.M. and M.C.). T.Q., G.Q., J.I.H. were supported by Australian Postgraduate Award. D.I. was supported by International Postgraduate Research Scholarship, and T.Q., G.X.Q., J.I.H., and D.I. were also supported by the John A. Lamberton scholarshi