The Scientific Method as a Scaffold to Enhance Communication Skills in Chemistry.

Scientific success in the field of chemistry depends upon the mastery of a wide range of soft skills, most notably scientific writing and speaking. However, training for scientific communication is typically limited at the undergraduate level, where students struggle to express themselves in a clear and logical manner. The underlying issue is deeper than basic technical skills; rather, it is a problem of students' unawareness of a fundamental and strategic framework for writing and speaking with a purpose. The methodology has been implemented for individual mentorship and in our regional summer research program to deliver a blueprint of thought and reasoning that endows students with the confidence and skills to become more effective communicators. Our didactic process intertwines undergraduate research with the scientific method and is partitioned into six steps, referred to as "phases", to allow for focused and deep thinking on the essential components of the scientific method. The phases are designed to challenge the student in their zone of proximal development so they learn to extract and ultimately comprehend the elements of the scientific method through focused written and oral assignments. Students then compile their newly acquired knowledge to create a compelling and logical story, using their persuasive written and oral presentations to complete a research proposal, final report, and formal 20 min presentation. We find that such an approach delivers the necessary guidance to promote the logical framework that improves writing and speaking skills. Over the past decade, we have witnessed both qualitative and quantitative gains in the students' confidence in their abilities and skills (developed by this process), preparing them for future careers as young scientists.

Parameterization of the GPR119 Receptor Agonist AR231453.

The GPR119 receptor is a class A G protein-coupled receptor expressed mainly in pancreatic beta cells. Since GPR119 receptor activation ameliorates Type 2 Diabetes through an increase in glucose-dependent insulin release, the development of new GPR119 receptor agonists would be worthwhile. A better understanding of the way agonists interact with the receptor would help to design better ligands for the receptor. It also would help to better understand the agonist mechanism of action. An understanding of how agonists interact with the receptor can be acquired using molecular dynamics simulations, which cannot be performed without having force field parameters for the ligand molecule. This study presents the development of CHARMM force field parameters for AR231453, the prototypical first potent and orally available GPR119 agonist, using the Force Field Tool Kit. The parameters are validated through Normal Mode Analysis calculations and molecular dynamics simulations in combination with infrared spectroscopy. © 2017 Wiley Periodicals, Inc.

Design, synthesis, and analysis of antagonists of GPR55: Piperidine-substituted 1,3,4-oxadiazol-2-ones.

A series of 1,3,4-oxadiazol-2-ones was synthesized and tested for activity as antagonists at GPR55 in cellular beta-arrestin redistribution assays. The synthesis was designed to be modular in nature so that a sufficient number of analogues could be rapidly accessed to explore initial structure-activity relationships. The design of analogues was guided by the docking of potential compounds into a model of the inactive form of GPR55. The results of the assays were used to learn more about the binding pocket of GPR55. With this oxadiazolone scaffold, it was determined that modification of the aryl group adjacent to the oxadiazolone ring was often detrimental and that the distal cyclopropane was beneficial for activity. These results will guide further exploration of this receptor.

Crucial positively charged residues for ligand activation of the GPR35 receptor.

GPR35 is a G protein-coupled receptor expressed in the immune, gastrointestinal, and nervous systems in gastric carcinomas and is implicated in heart failure and pain perception. We investigated residues in GPR35 responsible for ligand activation and the receptor structure in the active state. GPR35 contains numerous positively charged amino acids that face into the binding pocket that cluster in two distinct receptor regions, TMH3-4-5-6 and TMH1-2-7. Computer modeling implicated TMH3-4-5-6 for activation by the GPR35 agonists zaprinast and pamoic acid. Mutation results for the TMH1-2-7 region of GPR35 showed no change in ligand efficacies at the K1.32A, R2.65A, R7.33A, and K7.40A mutants. However, mutation of arginine residues in the TMH3-4-5-6 region (R4.60, R6.58, R3.36, R(164), and R(167) in the EC2 loop) had effects on signaling for one or both agonists tested. R4.60A resulted in a total ablation of agonist-induced activation in both the ß-arrestin trafficking and ERK1/2 activation assays. R6.58A increased the potency of zaprinast 30-fold in the pERK assay. The R(167)A mutant decreased the potency of pamoic acid in the ß-arrestin trafficking assay. The R(164)A and R(164)L mutants decreased potencies of both agonists. Similar trends for R6.58A and R(167)A were observed in calcium responses. Computer modeling showed that the R6.58A mutant has additional interactions with zaprinast. R3.36A did not express on the cell surface but was trapped in the cytoplasm. The lack of surface expression of R3.36A was rescued by a GPR35 antagonist, CID2745687. These results clearly show that R4.60, R(164), R(167), and R6.58 play crucial roles in the agonist initiated activation of GPR35.

Identification of the GPR55 antagonist binding site using a novel set of high-potency GPR55 selective ligands.

GPR55 is a class A G protein-coupled receptor (GPCR) that has been implicated in inflammatory pain, neuropathic pain, metabolic disorder, bone development, and cancer. Initially deorphanized as a cannabinoid receptor, GPR55 has been shown to be activated by non-cannabinoid ligands such as l-α-lysophosphatidylinositol (LPI). While there is a growing body of evidence of physiological and pathophysiological roles for GPR55, the paucity of specific antagonists has limited its study. In collaboration with the Molecular Libraries Probe Production Centers Network initiative, we identified a series of GPR55 antagonists using a ß-arrestin, high-throughput, high-content screen of ~300000 compounds. This screen yielded novel, GPR55 antagonist chemotypes with IC50 values in the range of 0.16-2.72 µM [Heynen-Genel, S., et al. (2010) Screening for Selective Ligands for GPR55: Antagonists (ML191, ML192, ML193) (Bookshelf ID NBK66153; PMID entry 22091481)]. Importantly, many of the GPR55 antagonists were completely selective, with no agonism or antagonism against GPR35, CB1, or CB2 up to 20 µM. Using a model of the GPR55 inactive state, we studied the binding of an antagonist series that emerged from this screen. These studies suggest that GPR55 antagonists possess a head region that occupies a horizontal binding pocket extending into the extracellular loop region, a central ligand portion that fits vertically in the receptor binding pocket and terminates with a pendant aromatic or heterocyclic ring that juts out. Both the region that extends extracellularly and the pendant ring are features associated with antagonism. Taken together, our results provide a set of design rules for the development of second-generation GPR55 selective antagonists.

The importance of hydrogen bonding and aromatic stacking to the affinity and efficacy of cannabinoid receptor CB2 antagonist, 5-(4-chloro-3-methylphenyl)-1-[(4-methylphenyl)methyl]-N-[(1S,2S,4R)-1,3,3-trimethylbicyclo[2.2.1]hept-2-yl]-1H-pyrazole-3-carboxamide (SR144528).

Despite the therapeutic promise of the subnanomolar affinity cannabinoid CB2 antagonist, 5-(4-chloro-3-methylphenyl)-1-[(4-methylphenyl)methyl]-N-[(1S,2S,4R)-1,3,3-trimethylbicyclo[2.2.1]hept-2-yl]-1H-pyrazole-3-carboxamide (SR144528, 1), little is known about its binding site interactions and no primary interaction site for 1 at CB2 has been identified. We report here the results of Glide docking studies in our cannabinoid CB2 inactive state model that were then tested via compound synthesis, binding, and functional assays. Our results show that the amide functional group of 1 is critical to its CB2 affinity and efficacy and that aromatic stacking interactions in the TMH5/6 aromatic cluster of CB2 are also important. Molecular modifications that increased the positive electrostatic potential in the region between the fenchyl and aromatic rings led to more efficacious compounds. This result is consistent with the EC-3 loop negatively charged amino acid, D275 (identified via Glide docking studies) acting as the primary interaction site for 1 and its analogues.

Identification of the GPR55 agonist binding site using a novel set of high-potency GPR55 selective ligands.

Marijuana is the most widely abused illegal drug, and its spectrum of effects suggests that several receptors are responsible for the activity. Two cannabinoid receptor subtypes, CB1 and CB2, have been identified, but the complex pharmacological properties of exogenous cannabinoids and endocannabinoids are not fully explained by their signaling. The orphan receptor GPR55 binds a subset of CB1 and CB2 ligands and has been proposed as a cannabinoid receptor. This designation, however, is controversial as a result of recent studies in which lysophosphatidylinositol (LPI) was identified as a GPR55 agonist. Defining a biological role for GPR55 requires GPR55 selective ligands that have been unavailable. From a ß-arrestin, high-throughput, high-content screen of 300000 compounds run in collaboration with the Molecular Libraries Probe Production Centers Network initiative (PubChem AID1965), we identified potent GPR55 selective agonists. By modeling of the GPR55 activated state, we compared the GPR55 binding conformations of three of the novel agonists obtained from the screen, CID1792197, CID1172084, and CID2440433 (PubChem Compound IDs), with that of LPI. Our modeling indicates the molecular shapes and electrostatic potential distributions of these agonists mimic those of LPI; the GPR55 binding site accommodates ligands that have inverted-L or T shapes with long, thin profiles that can fit vertically deep in the receptor binding pocket while their broad head regions occupy a horizontal binding pocket near the GPR55 extracellular loops. Our results will allow the optimization and design of second-generation GPR55 ligands and provide a means for distinguishing GPR55 selective ligands from those interacting with cannabinoid receptors.

Parameterization of Org27569: an allosteric modulator of the cannabinoid CB1 G protein-coupled receptor.

The cannabinoid CB1 receptor is a class A G protein-coupled receptor (GPCR) that is the most widely expressed GPCR in the brain. Many GPCRs contain allosteric binding sites for endogenous and/or synthetic ligands, which are topographically distinct from the agonist-binding site that is known as the orthosteric site. While both endogenous and synthetic ligands that act at the CB1 orthosteric site have been known for some time, compounds that act at a CB1 allosteric site have only recently been discovered. The most studied of these is 5-chloro-3-ethyl-1H-indole-2-carboxylic acid [2-(4-piperidin-1-ylphenyl)ethyl]amide (Org27569). Because allosteric ligands are thought to act through conformational changes in the receptor that are transmitted from the allosteric to the orthosteric site, computational studies of the structural and dynamic interactions of Org27569 with the CB1 receptor are crucial to achieve a molecular level understanding of the basis of action of this important new class of compounds. To date, such computational studies have not been possible due to the lack of a complete set of molecular mechanics force field parameters for Org27569. Here, we present the development of missing CHARMM force field parameters for Org27569 using previously published methods and the validation and application of these new parameters using normal mode analysis and molecular dynamics simulations combined with experimental infrared measurements.

Lipid bilayer molecular dynamics study of lipid-derived agonists of the putative cannabinoid receptor, GPR55.

Both L-α-lysophosphatidylinositol (LPI) and 2-arachidonoyl-sn-glycero-3-phosphoinositol (2-AGPI) have been reported to activate the putative cannabinoid receptor, GPR55. Recent microsecond time-scale molecular dynamics (MD) simulations and isothiocyanate covalent labeling studies have suggested that a transmembrane helix 6/7 (TMH6/7) lipid pathway for ligand entry may be necessary for interaction with cannabinoid receptors. Because LPI and 2-AGPI are lipid-derived ligands, conformations that each assumes in the lipid bilayer are therefore likely important for their interaction with GPR55. We report here the results of 70 ns NAMD molecular dynamics (MD) simulations of LPI and of 2-AGPI in a fully hydrated bilayer of 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC). These simulations are compared with a 70 ns simulation of the cannabinoid CB1 receptor endogenous ligand, N-arachidonoylethanolamine (anandamide, AEA) in a POPC bilayer. These simulations revealed that (1) LPI and 2-AGPI sit much higher in the bilayer than AEA, with inositol headgroups that can at times be solvated completely by water; (2) the behavior of the acyl chains of AEA and 2-AGPI are similar in their flexibilities in the bilayer, while the acyl chain of LPI has reduced flexibility; and (3) both 2-AGPI and LPI can adopt a tilted headgroup orientation by hydrogen bonding to the phospholipid phosphate/glycerol groups or via intramolecular hydrogen bonding. This tilted head group conformation (which represents over 40% of the conformer population of LPI (42.2 ± 3.3%) and 2-AGPI (43.7 ± 1.4%)) may provide a low enough profile in the lipid bilayer for LPI and 2-AGPI to enter GPR55 via the putative TMH6/7 entry port.

GEC1-kappa opioid receptor binding involves hydrophobic interactions: GEC1 has chaperone-like effect.

We demonstrated previously that the protein GEC1 (glandular epithelial cell 1) bound to the human kappa opioid receptor (hKOPR) and promoted cell surface expression of the receptor by facilitating its trafficking along the secretory pathway. Here we showed that three hKOPR residues (Phe345, Pro346, and Met350) and seven GEC1 residues (Tyr49, Val51, Leu55, Thr56, Val57, Phe60, and Ile64) are indispensable for the interaction. Modeling studies revealed that the interaction was mediated via direct contacts between the kinked hydrophobic fragment in hKOPR C-tail and the curved hydrophobic surface in GEC1 around the S2 beta-strand. Intramolecular Leu44-Tyr109 interaction in GEC1 was important, likely by maintaining its structural integrity. Microtubule binding mediated by the GEC1 N-terminal domain was essential for the GEC1 effect. Expression of GEC1 also increased cell surface levels of the GluR1 subunit and the prostaglandin EP3.f receptor, which have FPXXM and FPXM sequences, respectively. With its widespread distribution in the nervous system and its predominantly hydrophobic interactions, GEC1 may have chaperone-like effects for many cell surface proteins along the biosynthesis pathway.

Analogs of JHU75528, a PET ligand for imaging of cerebral cannabinoid receptors (CB1): development of ligands with optimized lipophilicity and binding affinity.

Cyano analogs of Rimonabant with high binding affinity for the cerebral cannabinoid receptor (CB1) and with optimized lipophilicity have been synthesized as potential positron emission tomography (PET) ligands. The best ligands of the series are optimal targets for the future radiolabeling with PET isotopes and in vivo evaluation as radioligands with enhanced properties for PET imaging of CB1 receptors in human subjects. Extracellular electrophysiological recordings in rodent brain slices demonstrated that JHU75528, 4, the lead compound of the new series, has functional CB antagonist properties that are consistent with its structural relationship to Rimonabant. Molecular modeling analysis revealed an important role of the binding of the cyano group with the CB1 binding pocket.

Bisphosphonate inhibition of phosphoglycerate kinase: quantitative structure-activity relationship and pharmacophore modeling investigation.

We report the results of a three-dimensional quantitative structure-activity relationship (3D-QSAR) and pharmacophore modeling investigation of the interaction of the enzyme 3-phosphoglycerate kinase (PGK) with aryl and alkyl bisphosphonates. For the human enzyme, the IC50 values are predicted within a factor of 2 over the 240x experimental range in activity, while for the yeast enzyme, binding of the more flexible alkyl bisphosphonates is predicted within a factor of approximately 4 (over a 2500x range in activity). Pharmacophore models indicate the importance of two negative ionizable features, one hydrophobic feature, and one halogen feature, and docking studies indicate that bisphosphonates bind in a manner similar to the 3-phosphoglycerate molecule identified crystallographically. The results give a good account of the activities of a diverse range of bisphosphonate inhibitors and are of interest in the context of developing inhibitors of glycolysis in organisms that are totally reliant on glycolysis for ATP production, such as trypanosomatid parasites.

Bisphosphonate inhibition of the exopolyphosphatase activity of the Trypanosoma brucei soluble vacuolar pyrophosphatase.

Trypanosoma brucei, the causative agent of African trypanosomiasis, contains a soluble, vacuolar pyrophosphatase, TbVSP1, not present in humans, which is essential for the growth of bloodstream forms in their mammalian host. Here, we report the inhibition of a recombinant TbVSP1 expressed in Escherichia coli by a panel of 81 bisphosphonates. The IC50 values were found to vary from approximately 2 to 850 microM. We then used 3D QSAR (comparative molecular field and comparative molecular similarity index; CoMFA and CoMSIA) methods to analyze the enzyme inhibition results. The R2 values for the experimental versus the QSAR-predicted activities were 0.78 or 0.61 for CoMFA and 0.79 or 0.68 for CoMSIA, for two different alignments. The root-mean-square (rms) pIC50 error for the best CoMFA model was 0.41 for five test sets of five activity predictions, which translates to a factor of approximately 2.6 error in IC50 prediction. For CoMSIA, the rms pIC50 error and error factors were 0.35 and 2.2, respectively. In general, the most active compounds contained both a single aromatic ring and a hydrogen bond donor feature. Thirteen of the more potent compounds were then tested in vivo in a mouse model of T. brucei infection. The most active compound in vivo provided a 40% protection from death with no apparent side effects, suggesting that further development of such compounds may be of interest.

Pyridinium-1-yl bisphosphonates are potent inhibitors of farnesyl diphosphate synthase and bone resorption.

We report the design, synthesis and testing of a series of novel bisphosphonates, pyridinium-1-yl-hydroxy-bisphosphonates, based on the results of comparative molecular similarity indices analysis and pharmacophore modeling studies of farnesyl diphosphate synthase (FPPS) inhibition, human Vgamma2Vdelta2 T cell activation and bone resorption inhibition. The most potent molecules have high activity against an expressed FPPS from Leishmania major, in Dictyostelium discoideum growth inhibition, in gammadelta T cell activation and in an in vitro bone resorption assay. As such, they represent useful new leads for the discovery of new bone resorption, antiinfective and anticancer drugs.

A quantitative structure-activity relationship and pharmacophore modeling investigation of aryl-X and heterocyclic bisphosphonates as bone resorption agents.

We have used quantitative structure-activity relationship (QSAR) techniques, together with pharmacophore modeling, to investigate the relationships between the structures of a wide variety of geminal bisphosphonates and their activity in inhibiting osteoclastic bone resorption. For aryl-X (X = alkyl, oxyalkyl, and sulfanylalkyl) derivatives of pamidronate and one alendronate, a molecular field analysis (MFA) yielded an R(2) value of 0.900 and an F-test of 54 for a training set of 29 compounds. Using reduced training sets, the activities of 20 such compounds were predicted with an average error of 2.1 over a 4000x range in activity. Such good results were only obtained when using the X-ray crystallographic structure of farnesyl pyrophosphate (FPP) bound to the target enzyme, farnesyl pyrophosphate synthase (FPP synthase), to guide the initial molecular alignment. For a series of heterocyclic bisphosphonates, use of the MFA method yielded an R(2) of 0.873 and an F-test of 36 for a training set of 26 compounds. Using a reduced training set, the activities of 20 compounds were predicted with an average error of 2.5 over a 2000x range in activity. With the heterocyclic compounds, test calculations indicated the importance of correct choice of protonation of the heterocyclic rings. For example, thiazoles, pyrazoles, and triazoles have low ( approximately 2-3) pK(a) values and the derived bisphosphonates are inactive in bone resorption since they cannot readily be side chain protonated and are thus poor carbocation reactive intermediate analogues. On the other hand, aminothiazoles, imidazoles, pyridyl, and aminopyridyl species typically have pK(a) values in the range approximately 5-9 and, in the absence of unfavorable steric interactions, the corresponding bisphosphonates are generally good inhibitors. However, aminoimidazole bisphosphonates are generally less active, since their pK(a)s ( approximately 11) are so high, due to guanidinium-like resonance, that they cannot readily be deprotonated, which we propose results in poor cellular uptake. The results of pharmacophore modeling using the Catalyst program revealed the importance of two negative ionizable and one positive charge feature for both aryl-X and heterocyclic pharmacophores, together with the presence of a distal hydrophobic feature in the aryl bisphosphonate and a more proximal aromatic feature in the heterocyclic bisphosphonate pharmacophores. When taken together, these results show that it is now possible to predict the activity, within a factor of about 2.3, of a wide range of aryl-X and heterocyclic bisphosphonates. The results emphasize the importance of utilizing crystallographic structural information to guide the initial alignment of extended bisphosphonates, and in the case of heterocyclic bisphosphonates, the importance of side chain protonation state. These simple ideas may facilitate the design of other, novel bisphosphonates, of use in bone resorption therapy, and as antiparasitic and immunotherapeutic agents.