Computational Predictions of G Protein-Coupled Receptor Structures and Binding Sites

G protein-coupled receptors (GPCRs) are the largest family of proteins within the human genome. They consist of seven transmembrane (TM) helices, with a N-terminal region of varying length and structure on the extracellular side, and a C-terminus on the intracellular side. GPCRs are involved in transmitting extracellular signals to cells, and as such are crucial drug targets. Designing pharmaceuticals to target GPCRs is greatly aided by full-atom structural information of the proteins. In particular, the TM region of GPCRs is where small molecule ligands (much more bioavailable than peptide ligands) typically bind to the receptors. In recent years nearly thirty distinct GPCR TM regions have been crystallized. However, there are more than 1,000 GPCRs, leaving the vast majority of GPCRs with limited structural information. Additionally, GPCRs are known to exist in a myriad of conformational states in the body, rendering the static x-ray crystal structures an incomplete reflection of GPCR structures. In order to obtain an ensemble of GPCR structures, we have developed the GEnSeMBLE procedure to rapidly sample a large number of variations of GPCR helix rotations and tilts. The lowest energy GEnSeMBLE structures are then docked to small molecule ligands and optimized. The GPCR family consists of five subfamilies with little to no sequence homology between them: class A, B1, B2, C, and Frizzled/Taste2. Almost all of the GPCR crystal structures have been of class A GPCRs, and much is known about their conserved interactions and binding sites. In this work we particularly focus on class B1 GPCRs, and aim to understand that family’s interactions and binding sites both to small molecules and their native peptide ligands. Specifically, we predict the full atom structure and peptide binding site of the glucagon-like peptide receptor and the TM region and small molecule binding sites for eight other class B1 GPCRs: CALRL, CRFR1, GIPR, GLR, PACR, PTH1R, VIPR1, and VIPR2. Our class B1 work reveals multiple conserved interactions across the B1 subfamily as well as a consistent small molecule binding site centrally located in the TM bundle. Both the interactions and the binding sites are distinct from those seen in the more well-characterized class A GPCRs, and as such our work provides a strong starting point for drug design targeting class B1 proteins. We also predict the full structure of CXCR4 bound to a small molecule, a class A GPCR that was not closely related to any of the class A GPCRs at the time of the work.

Binding Studies of Neuronal Nicotinic Acetylcholine Receptors Expressing Unnatural Amino Acids

Nicotinic acetylcholine receptors are pentameric ligand-gated ion channels mediating fast synaptic transmission throughout the peripheral and central nervous systems. They have been implicated in various processes related to cognitive functions, learning and memory, arousal, reward, motor control and analgesia. Therefore, these receptors present alluring potential therapeutic targets for the treatment of pain, epilepsy, Alzheimer’s disease, Parkinson’s disease, Tourette’s syndrome, schizophrenia, anxiety, depression and nicotine addiction. The work detailed in this thesis focuses on binding studies of neuronal nicotinic receptors and aims to further our knowledge of subtype specific functional and structural information.

Chapter 1 is an introductory chapter describing the structure and function of nicotinic acetylcholine receptors as well as the methodologies used for the dissertation work described herein. There are several different subtypes of nicotinic acetylcholine receptors known to date and the subtle variations in their structure and function present a challenging area of study. The work presented in this thesis deals specifically with the α4β2 subtype of nicotinic acetylcholine receptor. This subtype assembles into 2 closely related stoichiometries, termed throughout this thesis as A3B2 and A2B3 after their respective subunit composition. Chapter 2 describes binding studies of select nicotinic agonists on A3B2 and A2B3 receptors determined by whole-cell recording. Three key binding interactions, a cation-π and two hydrogen bonds, were probed for four nicotinic agonists, acetylcholine, nicotine, smoking cessation drug varenicline (Chantix®) and the related natural product cytisine.

Results from the binding studies presented in Chapter 2 show that the major difference in binding of these four agonists to A3B2 and A2B3 receptors lies in one of the two hydrogen bond interactions where the agonist acts as the hydrogen bond acceptor and the backbone NH of a conserved leucine residue in the receptor acts as the hydrogen bond donor. Chapter 3 focuses on studying the effect of modulating the hydrogen bond acceptor ability of nicotine and epibatidine on A3B2 receptor function determined by whole-cell recording. Finally, Chapter 4 describes single-channel recording studies of varenicline binding to A2B3 and A3B2 receptors.

Investigation of competitive antagonist binding to the nicotinic acetylcholine receptor using voltage-jump and light-flash techniques

1. The effect of 2,2’-bis-[α-(trimethylammonium)methyl]azobenzene (2BQ), a photoisomerizable competitive antagonist, was studied at the nicotinic acetycholine receptor of Electrophorus electroplaques using voltage-jump and light-flash techniques.

2. 2BQ, at concentrations below 3 μΜ, reduced the amplitude of voltage-jump relaxations but had little effect on the voltage-jump relaxation time constants under all experimental conditions. At higher concentrations and voltages more negative than -150 mV, 2BQ caused significant open channel blockade.

3. Dose-ratio studies showed that the cis and trans isomers of 2BQ have equilibrium binding constants (K _ᵢ) of .33 and 1.0 μΜ, respectively. The binding constants determined for both isomers are independent of temperature, voltage, agonist concentration, and the nature of the agonist.

4. In a solution of predominantly cis-2BQ, visible-light flashes led to a net cis→trans isomerization and caused an increase in the agonist-induced current. This increase had at least two exponential components; the larger amplitude component had the same time constant as a subsequent voltage-jump relaxation; the smaller amplitude component was investigated using ultraviolet light flashes.

5. In a solution of predominantly trans-2BQ, UV-light flashes led to a net trans→cis isomerization and caused a net decrease in the agonist-induced current. This effect had at least two exponential components. The smaller and faster component was an increase in agonist-induced current and had a similar time constant to the voltage-jump relaxation. The larger component was a slow decrease in the agonist-induced current with rate constant approximately an order of magnitude less than that of the voltage-jump relaxation. This slow component provided a measure of the rate constant for dissociation of cis-2BQ (k_ = 60/s at 20°C). Simple modelling of the slope of the dose-rate curves yields an association rate constant of 1.6 x 10⁸/M/s. This agrees with the association rate constant of 1.8 x 10⁸/M/s estimated from the binding constant (K_i). The Q₁₀ of the dissociation rate constant of cis-2BQ was 3.3 between 6° and 20°C. The rate constants for association and dissociation of cis-28Q at receptors are independent of voltage, agonist concentration, and the nature of the agonist.

6. We have measured the molecular rate constants of a competitive antagonist which has roughly the same K _ᵢ as d-tubocurarine but interacts more slowly with the receptor. This leads to the conclusion that curare itself has an association rate constant of 4 x 10⁹/M/s or roughly as fast as possible for an encounter-limited reaction.

Incremental Control Synthesis for Robotics in the Presence of Temporal Logic Specifications

This thesis presents methods for incrementally constructing controllers in the presence of uncertainty and nonlinear dynamics. The basic setting is motion planning subject to temporal logic specifications. Broadly, two categories of problems are treated. The first is reactive formal synthesis when so-called discrete abstractions are available. The fragment of linear-time temporal logic (LTL) known as GR(1) is used to express assumptions about an adversarial environment and requirements of the controller. Two problems of changes to a specification are posed that concern the two major aspects of GR(1): safety and liveness. Algorithms providing incremental updates to strategies are presented as solutions. In support of these, an annotation of strategies is developed that facilitates repeated modifications. A variety of properties are proven about it, including necessity of existence and sufficiency for a strategy to be winning. The second category of problems considered is non-reactive (open-loop) synthesis in the absence of a discrete abstraction. Instead, the presented stochastic optimization methods directly construct a control input sequence that achieves low cost and satisfies a LTL formula. Several relaxations are considered as heuristics to address the rarity of sampling trajectories that satisfy an LTL formula and demonstrated to improve convergence rates for Dubins car and single-integrators subject to a recurrence task.

The binding and catalytic properties of lysozyme

The binding and catalytic properties of hen's egg white lysozyme have been studied by a variety of techniques. These studies show that the enzyme has three contiguous binding subsites, A, B, and C. The application of nuclear magnetic resonance (NMR) spectroscopy to probe the binding environment of several saccharides to lysozyme has demonstrated that the reducing end sugar rings of chitotriose, chitobiose and the β-form of N-acetylglucosamine all bind in subsite C. The central sugar ring of chitotriose and the sugar ring at the nonreducing end of chitobiose were found to bind in subsite B, while the nonreducing end sugar residue of chitotriose occupied subsite A. The dynamics of the binding process has also been investigated by NMR. The formation rate constant of chitobiose--and chitotriose-enzyme complexes were found to be about 4 X 10^-6 M^-1 sec^-1 with small activation energies.

The stereochemical path of the lysozyme catalyzed hydrolysis of glycosidic bonds has been shown to proceed with at least 99.7% retention of configuration at C-1 of the sugar. The lysozyme catalyzed hydrolysis of glucosidic bonds has been shown to be largely carbonium ion in character by virtue of the α-deuterium kinetic isotope effect (k_H/k_D = 1.11) observed for the reaction. It is probable that the mechanism of action of the enzyme involves a carbonium ion intermediate which is stereospecifically quenched by solvent. However, acetamido group participation cannot be ruled out for natural substrates.