Frank E. Blaney

Definition of the G protein-coupled receptor transmembrane bundle binding pocket and calculation of receptor similarities for drug design.

Recent advances in structural biology for G-protein-coupled receptors (GPCRs) have provided new opportunities to improve the definition of the transmembrane binding pocket. Here a reference set of 44 residue positions accessible for ligand binding was defined through detailed analysis of all currently available crystal structures. This was used to characterize pharmacological relationships of Family A/Rhodopsin family GPCRs, minimizing evolutionary influence from parts of the receptor that do not generally affect ligand binding. The resultant dendogram tended to group receptors according to endogenous ligand types, although it revealed subdivision of certain classes, notably peptide and lipid receptors. The transmembrane binding site reference set, particularly when coupled with a means of identifying the subset of ligand binding residues, provides a general paradigm for understanding the pharmacology/selectivity profile of ligands at Family A GPCRs. This has wide applicability to GPCR drug design problems across many disease areas.

1,2,4-Triazolyl Azabicyclo[3.1.0]hexanes: A New Series of Potent and Selective Dopamine D3 Receptor Antagonists

The discovery of new highly potent and selective dopamine (DA) D(3) receptor antagonists has recently allowed the characterization of the DA D(3) receptor in a range of preclinical animal models of drug addiction. A novel series of 1,2,4-triazol-3-yl-azabicyclo[3.1.0]hexanes, members of which showed a high affinity and selectivity for the DA D(3) receptor and excellent pharmacokinetic profiles, is reported here. Members of a group of derivatives from this series showed good oral bioavailability and brain penetration and very high in vitro affinity and selectivity for the DA D(3) receptor, as well as high in vitro potency for antagonism at this receptor. Several members of this series also significantly attenuate the expression of conditioned place preference (CPP) to nicotine and cocaine.

Discovery of novel 1-azoniabicyclo[2.2.2]octane muscarinic acetylcholine receptor antagonists.

A novel 4-hydroxyl(diphenyl)methyl substituted quinuclidine series was discovered as a very promising class of muscarinic antagonists. The structure-activity relationships of the connectivity of the diphenyl moiety to the quinuclidine core and around the ring nitrogen side chain are described. Computational docking studies using an homology model of the M(3) receptor readily explained the observed structure-activity relationship of the various compounds. Compound 14o was identified as a very potent, slowly reversible M(3) antagonist with a very long in vivo duration of bronchoprotection.

Synthesis and pharmacological characterization of novel druglike corticotropin-releasing factor 1 antagonists.

To identify new CRF(1) receptor antagonists, an attempt to modify the bis-heterocycle moiety present in the top region of the dihydropyrrole[2,3]pyridine template was made following new pharmacophoric hypothesis on the CRF(1) receptor antagonists binding pocket. In particular, the 2-thiazole ring, present in the previous series of compounds, was replaced by more hydrophilic non aromatic heterocycles able to make appropriate H-bond interactions with amino acid residues Thr192 and Tyr195. This exploration, followed by an accurate analysis of the substitution of the pendant aryl ring, enabled to identify in vitro potent compounds showing excellent pharmacokinetics and outstanding in vivo activity in animal models of anxiety, both in rodents and primates.

Modeling GPCR active state conformations: The β2‐adrenergic receptor

The recent publication of several G protein‐coupled receptor (GPCR) structures has increased the information available for homology modeling inactive class A GPCRs. Moreover, the opsin crystal structure shows some active features. We have therefore combined information from these two sources to generate an extensively validated model of the active conformation of the β2‐adrenergic receptor. Experimental information on fully active GPCRs from zinc binding studies, site‐directed spin labeling, and other spectroscopic techniques has been used in molecular dynamics simulations. The observed conformational changes reside mainly in transmembrane helix 6 (TM6), with additional small but significant changes in TM5 and TM7. The active model has been validated by manual docking and is in agreement with a large amount of experimental work, including site‐directed mutagenesis information. Virtual screening experiments show that the models are selective for β‐adrenergic agonists over other GPCR ligands, for (R)‐ over (S)‐β‐hydroxy agonists and for β2‐selective agonists over β1‐selective agonists. The virtual screens reproduce interactions similar to those generated by manual docking. The C‐terminal peptide from a model of the stimulatory G protein, readily docks into the active model in a similar manner to which the C‐terminal peptide from transducin, docks into opsin, as shown in a recent opsin crystal structure. This GPCR‐G protein model has been used to explain site‐directed mutagenesis data on activation. The agreement with experiment suggests a robust model of an active state of the β2‐adrenergic receptor has been produced. The methodology used here should be transferable to modeling the active state of other GPCRs. Proteins 2011.

Structural determinants of drugs acting on the Nav1.8 channel.

The aim of this work is to study the role of pore residues on drug binding in the NaV1.8 channel. Alanine mutations were made in the S6 segments, chosen on the basis of their roles in other NaV subtypes; whole cell patch clamp recordings were made from mammalian ND7/23 cells. Mutations of some residues caused shifts in voltage dependence of activation and inactivation, and gave faster time course of inactivation, indicating that the residues mutated play important roles in both activation and inactivation in the NaV1.8 channel. The resting and inactivated state affinities of tetracaine for the channel were reduced by mutations I381A, F1710A, and Y1717A (for the latter only inactivated state affinity was measured), and by mutation F1710A for the NaV1.8-selective compound A-803467, showing the involvement of these residues for each compound, respectively. For both compounds, mutation L1410A caused the unexpected appearance of a complete resting block even at extremely low concentrations. Resting block of native channels by compound A-803467 could be partially removed (“disinhibition”) by repetitive stimulation or by a test pulse after recovery from inactivation; the magnitude of the latter effect was increased for all the mutants studied. Tetracaine did not show this effect for native channels, but disinhibition was seen particularly for mutants L1410A and F1710A. The data suggest differing, but partially overlapping, areas of binding of A-803467 and tetracaine. Docking of the ligands into a three-dimensional model of the NaV1.8 channel gave interesting insight as to how the ligands may interact with pore residues.

Functional analysis of CYP2d6.31 variant : Homology modeling suggests possible disruption of redox partner interaction by Arg440his substitution

Cytochrome P450 2D6 (CYP2D6) is an important human drug‐metabolizing enzyme that exhibits a marked genetic polymorphism. Numerous CYP2D6 alleles have been characterized at a functional level, although the consequences for expression and/or catalytic activity of a substantial number of rare variants remain to be investigated. One such allele, CYP2D6*31, is characterized by mutations encoding three amino acid substitutions: Arg296Cys, Arg440His and Ser486Thr. The identification of this allele in an individual with an apparent in vivo poor metabolizer phenotype prompted us to analyze the functional consequence of these substitutions on enzyme activity using yeast as a heterologous expression system. We demonstrated that the Arg440His substitution, alone or in combination with Arg296Cys and/or Ser486Thr, altered the respective kinetic parameters [Km (μM) and kcat (min−1)] of debrisoquine 4‐hydroxylation (wild‐type, 25; 0.92; variants, 43–68; 0.05–0.11) and dextromethorphan O‐demethylation (wild‐type, 1; 4.72; variants, 12–23; 0.64‐1.43), such that their specificity constants (kcat/Km) were decreased by more than 95% compared to those observed with the wild‐type enzyme. The rates of oxidation of rac‐metoprolol at single substrate concentrations of 40 and 400 μM were also markedly decreased by approximately 90% with each CYP2D6 variant containing the Arg440His substitution. These in vitro data confirm that the CYP2D6*31 allele encodes an enzyme with a severely impaired but residual catalytic activity and, furthermore, that the Arg440His exchange alone is the inactivating mutation. A homology model of CYP2D6 based on the crystal structure of rabbit CYP2C5 locates Arg440 on the proximal surface of the protein. Docking the structure of the FMN domain of human cytochrome P450 reductase to the CYP2D6 model suggests that Arg440 is a key member of a cluster of basic amino acid residues important for reductase binding. Proteins 2005.

Modeling active GPCR conformations.

The most significant advance in modeling GPCR active states has been the β(2)-adrenergic receptor-Gs complex as this essentially transforms active-state modeling into homology modeling. Various different molecular dynamics-based approaches for modeling active states are presented, and a number of key applications discussed. These simulations have given insights into the activation pathway, conformational changes, dimerization, hydration, the ionic lock, ligand binding, protonation, and sodium binding. Crystallography and simulations have shown that the presence of agonist alone is unlikely to be sufficient to form the active state and that restraints applied to the G protein-binding region are required. The role of various microswitches in activation is discussed, including the controversial rotamer toggle switch. The importance of explicitly simulating experimental molecular probes to understand activation is highlighted, along with the need to ensure that such molecules are well parameterized. Approaches to loop modeling are discussed. We argue that the role of successful virtual screening against active models should not be overestimated as the main conformational changes on activation occur in the intracellular region.

Novel 5-HT1A/1B/1D receptors antagonists with potent 5-HT reuptake inhibitory activity

Novel 2-methyl-5-quinolinyl-1-piperazinylalkyl-3,4-dihydro-2H-1,4-benzoxazin-3-ones showing high affinities for the 5-HT(1A/1B/1D) receptors coupled with potent 5-HT reuptake inhibitory activity have been discovered. This is the first report describing docking of the lead compound 6-{2-[4-(2-methyl-5-quinolinyl)-1-piperazinyl]ethyl}-2H-1,4-benzoxazin-3(4H)-one 1, into a model of the 5-HT transporter and the 5-HT(1A) receptor model.

G-protein-coupled receptor dynamics: dimerization and activation models compared with experiment

Our previously derived models of the active state of the β2-adrenergic receptor are compared with recently published X-ray crystallographic structures of activated GPCRs (G-protein-coupled receptors). These molecular dynamics-based models using experimental data derived from biophysical experiments on activation were used to restrain the receptor to an active state that gave high enrichment for agonists in virtual screening. The β2-adrenergic receptor active model and X-ray structures are in good agreement over both the transmembrane region and the orthosteric binding site, although in some regions the active model is more similar to the active rhodopsin X-ray structures. The general features of the microswitches were well reproduced, but with minor differences, partly because of the unexpected X-ray results for the rotamer toggle switch. In addition, most of the interacting residues between the receptor and the G-protein were identified. This analysis of the modelling has also given important additional insight into GPCR dimerization: re-analysis of results on photoaffinity analogues of rhodopsin provided additional evidence that TM4 (transmembrane helix 4) resides at the dimer interface and that ligands such as bivalent ligands may pass between the mobile helices. A comparison, and discussion, is also carried out between the use of implicit and explicit solvent for active-state modelling.