A novel short neurotoxin, cobrotoxin c, from monocellate cobra (Naja kaouthia) venom: isolation and purification, primary and secondary structure determination, and tertiary structure modeling

A novel short neurotoxin, cobrotoxin c (CBT C) was isolated from the venom of monocellate cobra (Naja kaouthia) using a combination of ion-exchange chromatography and FPLC. Its primary structure was determined by Edman degradation. CBT C is composed of 61 amino acid residues. It differs from cobrotoxin b (CBT B) by only two amino acid substitutions, Thr/Ala11 and Arg/Thr56, which are not located on the functionally important regions by sequence similarity. However, the LD50 is 0.08 mg/g to mice, i.e. approximately five-fold higher than for CBT B. Strikingly, a structure-function relationship analysis suggests the existence of a functionally important domain on the outside of Loop III of CBT C. The functionally important basic residues on the outside of Loop III might have a pairwise interaction with alpha subunit, instead of gamma or delta subunits of the nicotinic acetylcholine receptor (nAChR). (C) 2002 Elsevier Science Inc. All rights reserved.

Structural analysis of human CCR2b and primate CCR2b by molecular modeling and molecular dynamics simulation

CCR2b, a chemokine receptor for MCP-1, -2, -3, -4, plays an important role in a variety of diseases involving infection, inflammation, and/or injury, as well as being a coreceptor for HIV-1 infection. Two models of human CCR2b (hCCR2b) were generated by h

Molecular modeling on human CCR5 receptors and complex with CD4 antigens and HIV-1 envelope glycoprotein gp120

AIM: To investigate the interaction between human CCR5 receptors (CCR5) and HIV-1 envelope glycoprotein gp120 (HIV-1 gp120) and HIV-1 receptor CD4 antigens (CD4). METHODS: The structurally con served regions (SCR) of human CCR5 was built by the SYBYL/Biopolymer module using the corresponding transmembrane (TM) domain of bacteriorhodopsin (bR) as the template. The coordinates for amino-ter minal residue sequence, and carboxyl-terminal residue sequence, extracellular and cytoplasmic loops were generated using LOOP SEARCH algorithm. Subsequently the structural model was merged into the complex with HIV-1 gp120 and CD4. RESULTS: Human CCR5 interacted with both an HIV-1 gp120 and CD4. The N-terminal residues (especially Met1 and Gln4) of human CCR5, contacted with CD4 residues, mainly 7Nith one span (56 - 59) of CD4 in electrostatic interaction and hydrogen-bonds. The binding sites of human CCR5 were buried in a hydrophobic center surrounded by a highly basic periphery. On the other hand, direct interatomic contacts were made between ? CCR5 residues and 6 gp120 amino-acid residues, which included van der Waals contacts, hydrophobic interaction, and hydrogen bonds. CONCLUSION: The interaction model should be helpful for rational design of novel anti-HIV drugs.

Characterize dynamic conformational space of human CCR5 extracellular domain by molecular modeling and molecular dynamics simulation

The chemokine receptor CCR5 is the receptor for several chemokines and major coreceptor for R5 human immunodeficiency virus type-1 strains entry into cell. Three-dimensional models of CCR5 were built by using homology modeling approach and 1 ns molecular dynamics (MD) simulation, because studies of site-directed mutagenesis and chimeric receptors have indicated that the N-terminus (Nt) and extracellular loops (ECLs) of CCR5 are important for ligands binding and viral fusion and entry, special attention was focused on disulfide bond function, conformational flexibility, hydrogen bonding, electrostatic interactions, and solvent-accessible surface area of Nt and ECLs of this protein part. We found that the extracellular segments of CCR5 formed a well-packet globular domain with complex interactions occurred between them in a majority of time of MID simulation, but Nt region could protrude from this domain sometimes. The disulfide bond Cys20-Cys269 is essential in controlling specific orientation of Nt region and maintaining conformational integrity of extracellular domain. RMS comparison analysis between conformers revealed the ECL1 of CCR5 stays relative rigid, whereas the ECL2 and Nt are rather flexible. Solvent-accessible surface area calculations indicated that the charged residues within Nt and ECL2 are often exposed to solvent. Integrating these results with available experimental data, a two-step gp120-CCR5 binding mechanism was proposed. The dynamic interaction of CCR5 extracellular domain with gp120 was emphasized. (C) 2004 Elsevier B.V. All rights reserved.

Impaired formation of beta-adrenergic receptor-nucleotide regulatory protein complexes in pseudohypoparathyroidism.

Decreased activity of the guanine nucleotide regulatory protein (N) of the adenylate cyclase system is present in cell membranes of some patients with pseudohypoparathyrodism (PHP-Ia) whereas others have normal activity of N (PHP-Ib). Low N activity in PHP-Ia results in a decrease in hormone (H)-stimulatable adenylate cyclase in various tissues, which might be due to decreased ability to form an agonist-specific high affinity complex composed of H, receptor (R), and N. To test this hypothesis, we compared beta-adrenergic agonist-specific binding properties in erythrocyte membranes from five patients with PHP-Ia (N = 45% of control), five patients with PHP-Ib (N = 97%), and five control subjects. Competition curves that were generated by increasing concentrations of the beta-agonist isoproterenol competing with [125I]pindolol were shallow (slope factors less than 1) and were computer fit to a two-state model with corresponding high and low affinity for the agonist. The agonist competition curves from the PHP-Ia patients were shifted significantly (P less than 0.02) to the right as a result of a significant (P less than 0.01) decrease in the percent of beta-adrenergic receptors in the high affinity state from 64 +/- 22% in PHP-Ib and 56 +/- 5% in controls to 10 +/- 8% in PHP-Ia. The agonist competition curves were computer fit to a "ternary complex" model for the two-step reaction: H + R + N in equilibrium HR + N in equilibrium HRN. The modeling was consistent with a 60% decrease in the functional concentration of N, and was in good agreement with the biochemically determined decrease in erythrocyte N protein activity. These in vitro findings in erythrocytes taken together with the recent observations that in vivo isoproterenol-stimulated adenylate cyclase activity is decreased in patients with PHP (Carlson, H. E., and A. S. Brickman, 1983, J. Clin. Endocrinol. Metab. 56:1323-1326) are consistent with the notion that N is a bifunctional protein interacting with both R and the adenylate cyclase. It may be that in patients with PHP-Ia a single molecular and genetic defect accounts for both decreased HRN formation and decreased adenylate cyclase activity, whereas in PHP-Ib the biochemical lesion(s) appear not to affect HRN complex formation.

Two Arginine-Glutamate Ionic Locks Near the Extracellular Surface of FFAR1 Gate Receptor Activation

Activation of a number of class A G protein-coupled receptors (GPCRs) is thought to involve two molecular switches, a rotamer toggle switch within the transmembrane domain and an ionic lock at the cytoplasmic surface of the receptor; however, the mechanism by which agonist binding changes these molecular interactions is not understood. Importantly, 80% of GPCRs including free fatty acid receptor 1 (FFAR1) lack the complement of amino acid residues implicated in either or both of these two switches; the mechanism of activation of these GPCRs is therefore less clear. By homology modeling, we identified two Glu residues (Glu-145 and Glu-172) in the second extracellular loop of FFAR1 that form putative interactions individually with two transmembrane Arg residues (Arg-183(5.39) and Arg-258(7.35)) to create two ionic locks. Molecular dynamics simulations showed that binding of agonists to FFAR1 leads to breakage of these Glu-Arg interactions. In mutagenesis experiments, breakage of these two putative interactions by substituting Ala for Glu-145 and Glu-172 caused constitutive receptor activation. Our results therefore reveal a molecular switch for receptor activation present on the extracellular surface of FFAR1 that is broken by agonist binding. Similar ionic locks between the transmembrane domains and the extracellular loops may constitute a mechanism common to other class A GPCRs also.

Evidence That Interspecies Polymorphism in the Human and Rat Cholecystokinin Receptor-2 Affects Structure of the Binding Site for the Endogenous Agonist Cholecystokinin

The cholecystokinin (CCK) receptor-2 exerts very important central and peripheral functions by binding the neuropeptides cholecystokinin or gastrin. Because this receptor is a potential therapeutic target, great interest has been devoted to the identification of efficient antagonists. However, interspecies genetic polymorphism that does not alter cholecystokinin-induced signaling was shown to markedly affect activity of synthetic ligands. In this context, precise structural study of the agonist binding site on the human cholecystokinin receptor-2 is a prerequisite to elucidating the molecular basis for its activation and to optimizing properties of synthetic ligands. In this study, using site-directed mutagenesis and molecular modeling, we delineated the binding site for CCK on the human cholecystokinin receptor-2 by mutating amino acids corresponding to that of the rat homolog. By doing so, we demonstrated that, although resembling that of rat homolog, the human cholecystokinin receptor-2 binding site also displays important distinct structural features that were demonstrated by susceptibility to several point mutations (F120A, Y189A, H207A). Furthermore, docking of CCK in the human and rat cholecystokinin receptor-2, followed by dynamic simulations, allowed us to propose a plausible structural explanation of the experimentally observed difference between rat and human cholecystokinin-2 receptors.

Chemically engineering ligand selectivity at the free fatty acid receptor 2 based on pharmacological variation between species orthologs

When it is difficult to develop selective ligands within a family of related G-protein-coupled receptors (GPCRs), chemically engineered receptors activated solely by synthetic ligands (RASSLs) are useful alternatives for probing receptor function. In the present work, we explored whether a RASSL of the free fatty acid receptor 2 (FFA2) could be developed on the basis of pharmacological variation between species orthologs. For this, bovine FFA2 was characterized, revealing distinct ligand selectivity compared with human FFA2. Homology modeling and mutational analysis demonstrated a single mutation in human FFA2 of C4.57G resulted in a human FFA2 receptor with ligand selectivity similar to the bovine receptor. This was exploited to generate human FFA2-RASSL by the addition of a second mutation at a known orthosteric ligand interaction site, H6.55Q. The resulting FFA2-RASSL displayed a >100-fold loss of activity to endogenous ligands, while responding to the distinct ligand sorbic acid with pEC(50) values for inhibition of cAMP, 5.83 ± 0.11; Ca(2+) mobilization, 4.63 ± 0.05; ERK phosphorylation, 5.61 ± 0.06; and dynamic mass redistribution, 5.35 ± 0.06. This FFA2-RASSL will be useful in future studies on this receptor and demonstrates that exploitation of pharmacological variation between species orthologs is a powerful method to generate novel chemically engineered GPCRs.-Hudson, B. D., Christiansen, E., Tikhonova, I. G., Grundmann, M., Kostenis, E., Adams, D. R., Ulven, T., Milligan, G. Chemically engineering ligand selectivity at the free fatty acid receptor 2 based on pharmacological variation between species orthologs.

Selective Orthosteric Free Fatty Acid Receptor 2 (FFA2) Agonists IDENTIFICATION OF THE STRUCTURAL AND CHEMICAL REQUIREMENTS FOR SELECTIVE ACTIVATION OF FFA2 VERSUS FFA3

Free fatty acid receptor 2 (FFA2; GPR43) is a G protein-coupled seven-transmembrane receptor for short-chain fatty acids (SCFAs) that is implicated in inflammatory and metabolic disorders. The SCFA propionate has close to optimal ligand efficiency for FFA2 and can hence be considered as highly potent given its size. Propionate, however, does not discriminate between FFA2 and the closely related receptor FFA3 (GPR41). To identify FFA2-selective ligands and understand the molecular basis for FFA2 selectivity, a targeted library of small carboxylic acids was examined using holistic, label-free dynamic mass redistribution technology for primary screening and the receptor-proximal G protein [S-35] guanosine 5'-(3-O-thio) triphosphate activation, inositol phosphate, and cAMP accumulation assays for hit confirmation. Structure-activity relationship analysis allowed formulation of a general rule to predict selectivity for small carboxylic acids at the orthosteric binding site where ligands with substituted sp(3)-hybridized alpha-carbons preferentially activate FFA3, whereas ligands with sp(2)- or sp-hybridized alpha-carbons prefer FFA2. The orthosteric binding mode was verified by site-directed mutagenesis: replacement of orthosteric site arginine residues by alanine in FFA2 prevented ligand binding, and molecular modeling predicted the detailed mode of binding. Based on this, selective mutation of three residues to their non-conserved counterparts in FFA3 was sufficient to transfer FFA3 selectivity to FFA2. Thus, selective activation of FFA2 via the orthosteric site is achievable with rather small ligands, a finding with significant implications for the rational design of therapeutic compounds selectively targeting the SCFA receptors.

CoMFA and homology-based models of the glycine binding site of N-methyl-D-aspartate receptor

Homology modeling was used to build 3D models of the N-methyl-D-aspartate (NMDA) receptor glycine binding site on the basis of an X-ray structure of the water-soluble AMPA-sensitive receptor. The docking of agonists and antagonists to these models was used to reveal binding modes of ligands and to explain known structure-activity relationships. Two types of quantitative models, 3D-QSAR/CoMFA and a regression model based on docking energies, were built for antagonists (derivatives of 4-hydroxy-2-quinolone, quinoxaline-2,3-dione, and related compounds). The CoMFA steric and electrostatic maps were superimposed on the homology-based model, and a close correspondence was marked. The derived computational models have permitted the evaluation of the structural features crucial for high glycine binding site affinity and are important for the design of new ligands.

Defining the molecular basis for the first potent and selective orthosteric agonists of the FFA2 free fatty acid receptor

FFA2 is a G protein-coupled receptor that responds to short chain fatty acids (SCFAs) and has generated interest as a therapeutic target for metabolic and inflammatory conditions. However, definition of its functions has been slowed by a dearth of selective ligands that can distinguish it from the closely related FFA3. At present, the only selective ligands described for FFA2 suffer from either poor potency, altered signaling due to allosteric modes of action, or a lack of function at non-human orthologs of the receptor. To address the need for novel selective ligands, we synthesized two compounds potentially having FFA2 activity and examined the molecular basis of their function. These compounds were confirmed to be potent and selective FFA2 agonists that interact with the orthosteric binding site. A combination of ligand structure-activity relationship, pharmacological analysis, homology modeling, species ortholog comparisons and mutagenesis studies were then employed to define the molecular basis of selectivity and function of these ligands. From this, we identified key residues within both extracellular loop 2 (ECL2) and the transmembrane domain (TM) regions of FFA2 critical for ligand function. One of these ligands was active with reasonable potency at rodent orthologs of FFA2 and demonstrated the role of FFA2 in the regulation of lipolysis in murine 3T3-L1 adipocytes. Together, these findings describe the first potent and selective FFA2 orthosteric agonists and demonstrate key aspects of ligand interaction within the orthosteric binding site of FFA2 that will be invaluable in future ligand development at this receptor.

Addressing Selective Polypharmacology of Antipsychotic Drugs Targeting the Bioaminergic Receptors through the Receptor Dynamic Conformational Ensembles

Selective polypharmacology, where a drug acts on multiple rather than single molecular targets involved in a disease, emerges to develop a structure-based system biology approach to design drugs selectively targeting a disease-active protein network. We focus on the bioaminergic receptors that belong to the group of integral membrane signalling proteins coupled to the G protein and represent targets for therapeutic agents against schizophrenia and depression. Among them, it has been shown that the serotonin (5-HT2A and 5-HT6), dopamine (D2 and D3) receptors induce a cognition-enhancing effect (group 1), while the histamine (H1) and serotonin (5-HT2C) receptors lead to metabolic side effects and the 5-HT2B serotonin receptor causes pulmonary hypertension (group 2). Thus, the problem arises to develop an approach that allows identifying drugs targeting only the disease-active receptors, i.e. group 1. The recent release of several crystal structures of the bioaminergic receptors, involving the D3 and H1 receptors provides the possibility to model the structures of all receptors and initiate a study of the structural and dynamic context of selective polypharmacology. In this work, we use molecular dynamics simulations to generate a conformational space of the receptors and subsequently characterize its binding properties applying molecular probe mapping. All-against-all comparison of the generated probe maps of the selected diverse conformations of all receptors with the Tanimoto similarity coefficient (Tc) enable to separate the receptors of group 1 from group 2. The pharmacophore built based on the Tc-selected receptor conformations, using the multiple probe maps discovers structural features that can be used to design molecules selective towards the receptors of group 1. The importance of several predicted residues to ligand selectivity is supported by the available mutagenesis and ligand structure-activity relationships studies. In addition, the Tc-selected conformations of the receptors for group 1 show good performance in isolation of known ligands from a random decoy. Our computational structure-based protocol to tackle selective polypharmacology of antipsychotic drugs could be applied for other diseases involving multiple drug targets, such as oncologic and infectious disorders.

Advances in GPCR modeling evaluated by the GPCR Dock 2013 assessment:Meeting new challenges

Despite tremendous successes of GPCR crystallography, the receptors with available structures represent only a small fraction of human GPCRs. An important role of the modeling community is to maximize structural insights for the remaining receptors and complexes. The community-wide GPCR Dock assessment was established to stimulate and monitor the progress in molecular modeling and ligand docking for GPCRs. The four targets in the present third assessment round presented new and diverse challenges for modelers, including prediction of allosteric ligand interaction and activation states in 5-hydroxytryptamine receptors 1B and 2B, and modeling by extremely distant homology for smoothened receptor. Forty-four modeling groups participated in the assessment. State-of-the-art modeling approaches achieved close-to-experimental accuracy for small rigid orthosteric ligands and models built by close homology, and they correctly predicted protein fold for distant homology targets. Predictions of long loops and GPCR activation states remain unsolved problems.

Status of GPCR modeling and docking as reflected by community-wide GPCR Dock 2010 assessment

The community-wide GPCR Dock assessment is conducted to evaluate the status of molecular modeling and ligand docking for human G protein-coupled receptors. The present round of the assessment was based on the recent structures of dopamine D3 and CXCR4 chemokine receptors bound to small molecule antagonists and CXCR4 with a synthetic cyclopeptide. Thirty-five groups submitted their receptor-ligand complex structure predictions prior to the release of the crystallographic coordinates. With closely related homology modeling templates, as for dopamine D3 receptor, and with incorporation of biochemical and QSAR data, modern computational techniques predicted complex details with accuracy approaching experimental. In contrast, CXCR4 complexes that had less-characterized interactions and only distant homology to the known GPCR structures still remained very challenging. The assessment results provide guidance for modeling and crystallographic communities in method development and target selection for further expansion of the structural coverage of the GPCR universe.

Quantitative single-molecule localization microscopy combined with rule-based modeling reveals ligand-induced TNF-R1 reorganization toward higher-order oligomers

We report on the assembly of tumor necrosis factor receptor 1 (TNF-R1) prior to ligand activation and its ligand-induced reorganization at the cell membrane. We apply single-molecule localization microscopy to obtain quantitative information on receptor cluster sizes and copy numbers. Our data suggest a dimeric pre-assembly of TNF-R1, as well as receptor reorganization toward higher oligomeric states with stable populations comprising three to six TNF-R1. Our experimental results directly serve as input parameters for computational modeling of the ligand-receptor interaction. Simulations corroborate the experimental finding of higher-order oligomeric states. This work is a first demonstration how quantitative, super-resolution and advanced microscopy can be used for systems biology approaches at the single-molecule and single-cell level.