Homology modeling and dynamics of the extracellular domain of rat and human neuronal nicotinic acetylcholine receptor subtypes α4β2 and α7

In recent years, it has become clear that the neuronal nicotinic acetylcholine receptor (nAChR) is a valid target in the treatment of a variety of diseases, including Alzheimer’s disease, anxiety, and nicotine addiction. As with most membrane proteins, information on the three-dimensional (3D) structure of nAChR is limited to data from electron microscopy, at a resolution that makes the application of structure-based design approaches to develop specific ligands difficult. Based on a high-resolution crystal structure of AChBP, homology models of the extracellular domain of the neuronal rat and human nAChR subtypes α4β2 and α7 (the subtypes most abundant in brain) were built, and their stability assessed with molecular dynamics (MD). All models built showed conformational stability over time, confirming the quality of the starting 3D model. Lipophilicity and electrostatic potential studies performed on the rat and human α4β2 and α7 nicotinic models were compared to AChBP, revealing the importance of the hydrophobic aromatic pocket and the critical role of the α-subunit Trp—the homolog of AChBP-Trp 143—for ligand binding. The models presented provide a valuable framework for the structure-based design of specific α4β2 nAChR subtype ligands aimed at improving therapeutic and diagnostic applications. Figure Electrostatic surface potential of the binding site cavity of the neuronal nicotinic acetylcholine receptor (nAChR). Nicotinic models performed with the MOLCAD program: a rat α7, b rat α4β2, c human α7, d human α4β2. All residues labeled are part of the α7 (a,c) or α4 (b,d) subunit with the exception of Phe 117, which belongs to subunit β2 (d). Violet Very negative, blue negative, yellow neutral, red very positive     

The structures of the human neuronal nicotinic acetylcholine receptor β2- and α3-subunit genes (CHRNB2 and CHRNA3)

The α4-subunit gene (CHRNA4) of the neuronal nicotinic acetylcholine receptor (nAChR) subunit family has recently been identified in two families as the gene responsible for autosomal dominant nocturnal frontal lobe epilepsy (ADNFLE), a rare monogenic idiopathic epilepsy. As a result of this finding, other subunits of the neuronal nAChR gene family are being considered as candidate genes for ADNFLE in families not linked to CHRNA4 and for other idiopathic epilepsies. α4-subunitsoften assemble together with β2-subunits (gene symbol CHRNB2) to build heteromeric nAChRs. The gene encoding another abundant AChR subunit, the α3-subunit gene (CHRNA3), is present with those encoding two other subunits, CHRNB4 and CHRNA5, in a gene cluster whose functional role is still unclear. Here we provide the information on the genomic structures of both the CHRNB2 and the CHRNA3 genes that is necessary for comprehensive mutational analyses, and we refine the genomic assignment of CHRNB2 on chromosome 1. Received: 5 August 1998 / Accepted: 13 October 1998     

Functional and structural interaction of (−)-lobeline with human α4β2 and α4β4 nicotinic acetylcholine receptor subtypes

To determine the pharmacologic activity of (−)-lobeline between human (h)α4β2 and hα4β4 nicotinic acetylcholine receptors (AChRs), functional and structural experiments were performed. The Ca²⁺ influx results established that (−)-lobeline neither actives nor enhances the function of the studied AChR subtypes, but competitively inhibits hα4β4 AChRs with potency ∼10-fold higher than that for hα4β2 AChRs. This difference is due to a higher binding affinity for the [³H]cytisine sites at hα4β4 compared to hα4β2 AChRs, which, in turn, can be explained by our molecular dynamics (MD) results: (1) higher stability of (−)-lobeline and its hydrogen bonds within the α4β4 pocket compared to the α4β2 pocket, (2) (−)-lobeline promotes Loop C to cap the binding site at the α4β4 pocket, but forces Loop C to get apart from the α4β2 pocket, precluding the gating process elicited by agonists, and (3) the orientation of (−)-lobeline within the α4β4, but not the α4β2, subpocket, promoted by the t− (or t+) rotameric state of α4-Tyr98, remains unchanged during the whole MD simulation. This study gives a detailed view of the molecular and dynamics events evoked by (−)-lobeline supporting the differential binding affinity and subsequent inhibitory potency between hα4β2 and hα4β4 AChRs, and supports the possibility that the latter subtype is also involved in its activity.     

Homology modeling and flex-ligand docking studies on the guinea pig β2 adrenoceptor: structural and experimental similarities/ differences with the human β2

The trachea of a guinea pig is widely used in drug development assays focused on the treatment of pulmonary diseases. Some of these drugs relax the airways by binding to the guinea pig β₂-adrenoceptor (Gβ₂AR). In this work, the amino acid sequence of the Gβ₂AR was searched to carry out homology modeling, using the Swiss-Model server, with the human β₂AR as the parent template. The Gβ₂AR 3-D structure was structurally and energetically optimized in vacuo using NAMD 2.6 program. The refined 3-D model obtained was used for further study. Molecular docking simulations were performed by testing a set of well-known β₂AR ligands using the AutoDock 3.0.5 program. The results show that the homology model of Gβ₂AR has a 3-D structure very similar to the crystal structure of recently studied human β₂AR. This was also corroborated by identity (94.23%), Ramachandran map, and docking results. The theoretical simulation showed that the ligands bind at sites that are similar to those reported for the human β₂AR. The R-enantiomer ligands showed correlation with in vitro data. We have obtained a Gβ₂AR 3-D model which can be used to carry out computational screening as a complementary tool during the drug design and experimental tests under guinea pig models.     

Analysis of L-DOPA and droxidopa binding to human β2-adrenergic receptor

Over the last two decades, an increasing number of studies has been devoted to a deeper understanding of the molecular process involved in the binding of various agonists and antagonists to active and inactive conformations of β₂-adrenergic receptor (β₂AR). The 3.2 Å x-ray crystal structure of human β₂AR active state in combination with the endogenous low affinity agonist adrenaline offers an ideal starting structure for studying the binding of various catecholamines to adrenergic receptors. We show that molecular docking of levodopa (L-DOPA) and droxidopa into rigid and flexible β₂AR models leads for both ligands to binding anchor sites comparable to those experimentally reported for adrenaline, namely D113/N312 and S203/S204/S207 side chains. Both ligands have a hydrogen bond network that is extremely similar to those of noradrenaline and dopamine. Interestingly, redocking neutral and protonated versions of adrenaline to rigid and flexible β₂AR models results in binding poses that are more energetically stable and distinct from the x-ray crystal structure. Similarly, lowest energy conformations of noradrenaline and dopamine generated by docking into flexible β₂AR models had binding free energies lower than those of best poses in rigid receptor models. Furthermore, our findings show that L-DOPA and droxidopa molecules have binding affinities comparable to those predicted for adrenaline, noradrenaline, and dopamine, which are consistent with previous experimental and computational findings and supported by the molecular dynamics simulations of β₂AR-ligand complexes performed here.     

Distinct and shared roles of β-arrestin-1 and β-arrestin-2 on the regulation of C3a receptor signaling in human mast cells

Background

The complement component C3a induces degranulation in human mast cells via the activation of cell surface G protein coupled receptors (GPCR; C3aR). For most GPCRs, agonist-induced receptor phosphorylation leads to the recruitment of β-arrestin-1/β-arrestin-2; resulting in receptor desensitization and internalization. Activation of GPCRs also leads to ERK1/2 phosphorylation via two temporally distinct pathways; an early response that reflects G protein activation and a delayed response that is G protein independent but requires β-arrestins. The role of β-arrestins on C3aR activation/regulation in human mast cells, however, remains unknown.

Methodology/Principal Findings

We utilized lentivirus short hairpin (sh)RNA to stably knockdown the expression of β-arrestin-1 and β-arrrestin-2 in human mast cell lines, HMC-1 and LAD2 that endogenously expresses C3aR. Silencing β-arrestin-2 attenuated C3aR desensitization, blocked agonist-induced receptor internalization and rendered the cells responsive to C3a for enhanced NF-κB activity as well as chemokine generation. By contrast, silencing β-arrestin-1 had no effect on these responses but resulted in a significant decrease in C3a-induced mast cell degranulation. In shRNA control cells, C3a caused a transient ERK1/2 phosphorylation, which peaked at 5 min but disappeared by 10 min. Knockdown of β-arrestin-1, β-arrestin-2 or both enhanced the early response to C3a and rendered the cells responsive for ERK1/2 phosphorylation at later time points (10–30 min). Treatment of cells with pertussis toxin almost completely blocked both early and delayed C3a-induced ERK1/2 phosphorylation in β-arrestin1/2 knockdown cells.

Conclusion/Significance

This study demonstrates distinct roles for β-arrestins-1 and β-arrestins-2 on C3aR desensitization, internalization, degranulation, NF-κB activation and chemokine generation in human mast cells. It also shows that both β-arrestin-1 and β-arrestin-2 play a novel and shared role in inhibiting G protein-dependent ERK1/2 phosphorylation. These findings reveal a new level of complexity for C3aR regulation by β-arrestins in human mast cells.     

Voltage-dependent and -independent block of α-amino-3-hydroxy-5-methylisoxazole-4-propionate receptor channels

Polyamine-containing toxins and synthetic dicationic derivatives of adamantane and phenylcyclohexyl selectively antagonize Ca(2+)-permeable α-amino-3-hydroxy-5-methylisoxazole-4-propionate (AMPA) receptor channels. These compounds demonstrate voltage-dependent open-channel block and are trapped by closed channels. In this study, we describe an alternative mechanism of non-competitive AMPA receptor inhibition caused by 9-aminoacridine and some of its derivatives. These compounds exhibit similar potency against Ca(2+)-permeable and Ca(2+)-impermeable AMPA receptors. The inhibition is largely voltage-independent, binding and unbinding do not require presence of agonist. We conclude that 9-aminoacridine binds to a shallow site in the AMPA receptor, which is located above the activation gate. A comparison of three-dimensional structures of the antagonists suggests that the 'V-like' shape of the hydrophobic headgroup favors voltage-dependent binding to the deep site in the channel pore, whereas the compounds possessing flat aromatic headgroups preferably bind to the shallow site. The characterization of the novel mechanism of AMPA receptor channel antagonism opens a way to develop a new family of pharmacological agents, which can be of scientific and practical importance.     

Unique assignment of inter-subunit association in GABAA α1β3γ2 receptors determined by molecular modeling

Recent publications defined requirements for inter-subunit contacts in a benzodiazepine-sensitive GABA_A receptor (GABA_ARα1β3γ2). There is strong evidence that the heteropentameric receptor contains two α1, two β3, and one γ2 subunit. However, the available data do not distinguish two possibilities: When viewed clockwise from an extracellular viewpoint the subunits could be arranged in either γ2β3α1β3α1 or γ2α1β3α1β3 configurations. Here we use molecular modeling to thread the relevant GABA_AR subunit sequences onto a template of homopentameric subunits in the crystal structure of the acetylcholine binding protein (AChBP). The GABA_A sequences are known to have 15-18% identity with the acetylcholine binding protein and nearly all residues that are conserved within the nAChR family are present in AChBP. The correctly aligned GABA_A sequences were threaded onto the AChBP template in the γ2β3α1β3α1 or γ2α1β3α1β3  arrangements. Only the γ2α1β3α1β3 arrangement satisfied three known criteria: (1) α1 His¹⁰² binds at the γ2 subunit interface in proximity to γ2 residues Thr¹⁴², Phe⁷⁷, and Met¹³⁰; (2) α1 residues 80-100 bind near γ2 residues 91-104; and (3) α1 residues 58-67 bind near the β3 subunit interface. In addition to predicting the most likely inter-subunit arrangement, the model predicts which residues form the GABA and benzodiazepine binding sites.     

Distribution of GDNF family receptor α3 and RET in rat and human non-neural tissues

The neurotrophic growth factor artemin binds selectively to GDNF family receptor α3 (GFRα3), forming a molecular complex with the co-receptor RET which mediates downstream signaling. This signaling pathway has been demonstrated to play an important role in the survival and maintenance of nociceptive sensory neurons and in the development of sympathetic neurons. However, the presence and potential role of this artemin-responsive pathway in non-neural tissues has not been fully explored to-date. To study the distribution of GFRα3 and RET in adult rat and human non-neural tissues, we carried out a comprehensive immunohistochemical study. We stained major organs from the digestive, urinary, reproductive, immune, respiratory and endocrine systems, and from other systems (cardiovascular, skeletal muscle), as well as regions of the nervous system for comparison. In both rat and human, the majority of non-neural cells did not exhibit detectable GFRα3-like immunoreactivity. In the rat, GFRα3- and RET-like staining were found in the same non-neural cell type only in kidney. In the human digestive and reproductive systems, a subset of epithelial cells exhibited GFRα3- and RET-like staining, suggesting co-localization. In other tissues, sub-populations of cells expressed either GFRα3- or RET-like immunoreactivity. The functional consequences of GFRα3 expression in non-neural cells remain to be determined.     

Defined concatenated α6α1β3γ2 GABAA receptor constructs reveal dual action of pyrazoloquinolinone allosteric modulators

Pyrazoloquinolinones (PQs) have been extensively studied as modulators of GABA_A receptors with different subunit composition, exerting modulatory effects by binding at α+/β- interfaces of GABA_A receptors. PQs with a substituent in position R7 have been reported to preferentially modulate α6- subunit containing GABA_A receptors which are mostly expressed in the cerebellum but were also found in the olfactory bulb, in the cochlear nucleus, in the hippocampus and in the trigeminal sensory pathway. They are considered potentially interesting in the context of sensori-motor gating deficits, depressive-like behavior, migraine and orofacial pain. Here we explored the option to modify the lead ligands’ R7 position. In the compound series we observed two different patterns of allosteric modulation in recombinantly expressed α6β3γ2 receptors, namely monophasic and biphasic positive modulation. In the latter case the additional phase occurred in the nanomolar range, while all compounds displayed robust modulation in the micromolar range. Nanomolar, near silent binding has been reported to occur at benzodiazepine binding sites, but was not investigated at the diazepam insensitive α6+/γ2- interface. To clarify the mechanism underlying the biphasic effect we tested one of the compounds in concatenated receptors. In these constructs the subunits are covalently linked, allowing to form either the α6+/γ2- interface, or the α6+/β3- interface, to study the resulting modulation. With this approach we were able to ascribe the nanomolar modulation to the α6+/γ2- interface. While not all compounds display the nanomolar phase, the strong modulation at the α6+/β3 interface proved to be tolerant for all tested R7 groups. This provides the future option to introduce e.g. isotope labelled or fluorescent moieties or substituents that enhance solubility and bioavailability.     

Production and purification of recombinant human α2C2 adrenergic receptor using Saccharomyces cerevisiae

The objective is to generate milligram quantities of recombinant human ₂C2 adrenergic receptor for X-ray crystallographic studies. It has been cloned in Saccharomyces cerevisiae, and the production level is at best about 13 pmol/mg of membrane protein, as estimated by radio-ligand binding assay. The receptor is solubilized with sucrose monolaurate followed by immunoaffinity purification and reconstitution into phospholipid vesicles. The efficiency of solubilization and immuno-purification are 60% and 91%, respectively.     

The TF1-ATPase and ATPase activities of assembled α3β3γ, α3β3γδ, and α3β3γε complexes are stimulated by low and inhibited by high concentrations of rhodamine 6G whereas the dye only inhibits the α3β3, and α3β3δ complexes

The ATPase activity of the F₁-ATPase from the thermophilic bacterium PS3 is stimulated at concentrations of rhodamine 6G up to about 10 µM where 70% stimulation is observed at 36°C. Half maximal stimulation is observed at about 3 µM dye. At rhodamine 6G concentrations greater than 10 µM, ATPase activity declines with 50% inhibition observed at about 75 µM dye. The ATPase activities of the ₃₃ and ₃₃ complexes assembled from isolated subunits of TF₁ expressed inE. coli deleted of theunc operon respond to increasing concentrations of rhodamine 6G nearly identically to the response of TF₁. In contrast, the ATPase activities of the ₃₃ and ₃₃ complexes are only inhibited by rhodamine 6G with 50% inhibition observed, respectively, at 35 and 75 µM dye at 36°C. The ATPase activity of TF₁ is stimulated up to 4-fold by the neutral detergent, LDAO. In the presence of stimulating concentrations of LDAO, the ATPase activity of TF₁ is no longer stimulated by rhodamine 6G, but rather, it is inhibited with 50% inhibition observed at about 30 µM dye at 30°C. One interpretation of these results is that binding of rhodamine 6G to a high-affinity site on TF₁ stimulates ATPase activity and unmasks a low-affinity, inhibitory site for the dye which is also exposed by LDAO.     

Discovery of highly potent and selective biphenylacylsulfonamide-based β3-adrenergic receptor agonists and molecular modeling based on the solved X-ray structure of the β2-adrenergic receptor: Part 6

As an extension of research, we have investigated modification of left-hand side (LHS) of biphenyl analogues containing an acylsulfonamide moiety in the development of potent and selective human β₃-adrenergic receptor (AR) agonists. Result of structure–activity relationships (SAR) and cassette-dosing evaluation in dogs showed that the hydroxynorephedrine analogue 16 had an excellent balance of in vitro and in vivo potency with pharmacokinetic profiles. In addition, to facilitate structure-based drug design (SBDD), we also have performed a docking study of biphenyl analogues based on the X-ray structure of the β₂-adrenergic receptor.     

Identification of a negative allosteric site on human α4β2 and α3β4 neuronal nicotinic acetylcholine receptors

Acetylcholine-based neurotransmission is regulated by cationic, ligand-gated ion channels called nicotinic acetylcholine receptors (nAChRs). These receptors have been linked to numerous neurological diseases and disorders such as Alzheimer's disease, Parkinson's disease, and nicotine addiction. Recently, a class of compounds has been discovered that antagonize nAChR function in an allosteric fashion. Models of human α4β2 and α3β4 nicotinic acetylcholine receptor (nAChR) extracellular domains have been developed to computationally explore the binding of these compounds, including the dynamics and free energy changes associated with ligand binding. Through a blind docking study to multiple receptor conformations, the models were used to determine a putative binding mode for the negative allosteric modulators. This mode, in close proximity to the agonist binding site, is presented in addition to a hypothetical mode of antagonism that involves obstruction of C loop closure. Molecular dynamics simulations and MM-PBSA free energy of binding calculations were used as computational validation of the predicted binding mode, while functional assays on wild-type and mutated receptors provided experimental support. Based on the proposed binding mode, two residues on the β2 subunit were independently mutated to the corresponding residues found on the β4 subunit. The T58K mutation resulted in an eight-fold decrease in the potency of KAB-18, a compound that exhibits preferential antagonism for human α4β2 over α3β4 nAChRs, while the F118L mutation resulted in a loss of inhibitory activity for KAB-18 at concentrations up to 100 μM. These results demonstrate the selectivity of KAB-18 for human α4β2 nAChRs and validate the methods used for identifying the nAChR modulator binding site. Exploitation of this site may lead to the development of more potent and subtype-selective nAChR antagonists which may be used in the treatment of a number of neurological diseases and disorders.     

Genomic organization and mapping of the human and mouse neuronal β2-nicotinic acetylcholine receptor genes

As a first step in determining whether there are polymorphisms in the nicotinic acetylcholine receptor (nAChR) genes that are associated with nicotine addiction, we isolated genomic clones of the β2-nAChR genes from human and mouse BAC libraries. Although cDNA sequences were available for the human gene, only the promoter sequence had been reported for the mouse gene. We determined the genomic structures by sequencing 12 kb of the human gene and over 7 kb of the mouse gene. While the sizes of exons in the mouse and human genes are the same, the introns differ in size. Both promoters have a high GC content (60–80%) proximal to the AUG and share a neural-restrictive silencer element (NRSE), but overall sequence identity is only 72%. Using a 6-Mb YAC contig of Chr 1, we mapped the human β2-nAChR gene, CHRNB2, to 1q21.3 with the order of markers cen, FLG, IVL, LOR, CHRNB2, tel. The mouse gene, Acrb2, had previously been mapped to Chr 3 in a region orthologous to human Chr 1. We refined mapping of the mouse gene and other markers on a radiation hybrid panel of Chr 3 and found the order cen, Acrb2, Lor, Iv1, Flg, tel. Our results indicate that this cluster of markers on human Chr 1 is inverted with respect to its orientation on the chromosome compared with markers in the orthologous region of mouse Chr 3. Received: 26 January 1999 / Accepted: 10 May 1999     

Physiological and Pathological Roles of α3β1 Integrin

α3β1 integrin has been considered to be a mysterious adhesion molecule due to the pleiotropy in its ligand-binding specificity. However, recent studies have identified laminin isoforms as high-affinity ligands for this integrin, and demonstrated that α3β1 integrin plays a number of essential roles in development and differentiation, mainly by mediating the establishment and maintenance of epithelial tissues. Furthermore, α3β1 integrin is also implicated in many other biological phenomena, including cell growth and apoptosis, angiogenesis and neural functions. This integrin receptor forms complexes with various other membrane proteins, such as the transmembrane-4 superfamily proteins (tetraspanins), cytoskeletal proteins and signaling molecules. Recently, lines of evidence have been reported showing that complex formation regulates integrin functions in cell adhesion and migration, signal transduction across cell membranes, and cytoskeletal organization. In addition to these roles in physiological processes, α3β1 integrin performs crucial functions in various pathological processes, especially in wound healing, tumor invasion and metastasis, and infection by pathogenic microorganisms.This revised version was published online in June 2005 with a corrected cover date.     

Expression of integrins α3β1 and α5β1 and GlcNAc β1,6 glycan branching influences metastatic melanoma cell migration on fibronectin

Acquisition of metastatic potential is accompanied by changes in cell surface N-glycosylation. One of the best-studied changes is increased expression of N-acetylglucosaminyltransferase V enzyme (GnT-V) and its products, β1,6-branched N-linked oligosaccharides, observed in the tumorigenesis of many cancers. In this study we demonstrate that during the transition from the vertical growth phase (VGP) (WM793 cell line) to the metastatic stage (WM1205Lu line), β1,6 glycosylation of melanoma cell surface proteins increases as a consequence of elevated expression of the GnT-V-encoding Mgat-5 gene. Treatment with swainsonine led to reduced cell motility on fibronectin in both cell lines; the effect was stronger in metastatic cells, probably due to the higher content of GlcNAc β1,6-branched glycans on the main fibronectin receptors – integrins α5β1 and α3β1. Our results show that GlcNAc β1,6 N-glycosylation of cell surface receptors, which increases with the aggressiveness of melanoma cells, is an important factor influencing melanoma cell migration.