A Comparative Study of Molecular Structure, Pka, Lipophilicity, Solubility, Absorption and Polar Surface Area of Some Antiplatelet Drugs
Total Page:16
File Type:pdf, Size:1020Kb
International Journal of Molecular Sciences Article A Comparative Study of Molecular Structure, pKa, Lipophilicity, Solubility, Absorption and Polar Surface Area of Some Antiplatelet Drugs Milan Remko 1,*, Anna Remková 2 and Ria Broer 3 1 Department of Pharmaceutical Chemistry, Faculty of Pharmacy, Comenius University in Bratislava, Odbojarov 10, SK-832 32 Bratislava, Slovakia 2 Department of Internal Medicine, Faculty of Medicine, Slovak Medical University, Limbová 12, SK–833 03 Bratislava, Slovakia; [email protected] 3 Department of Theoretical Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands; [email protected] * Correspondence: [email protected]; Tel.: +421-2-5011-7291 Academic Editor: Michael Henein Received: 18 February 2016; Accepted: 11 March 2016; Published: 19 March 2016 Abstract: Theoretical chemistry methods have been used to study the molecular properties of antiplatelet agents (ticlopidine, clopidogrel, prasugrel, elinogrel, ticagrelor and cangrelor) and several thiol-containing active metabolites. The geometries and energies of most stable conformers of these drugs have been computed at the Becke3LYP/6-311++G(d,p) level of density functional theory. Computed dissociation constants show that the active metabolites of prodrugs (ticlopidine, clopidogrel and prasugrel) and drugs elinogrel and cangrelor are completely ionized at pH 7.4. Both ticagrelor and its active metabolite are present at pH = 7.4 in neutral undissociated form. The thienopyridine prodrugs ticlopidine, clopidogrel and prasugrel are lipophilic and insoluble in water. Their lipophilicity is very high (about 2.5–3.5 logP values). The polar surface area, with regard to the structurally-heterogeneous character of these antiplatelet drugs, is from very large interval of values of 3–255 Å2. Thienopyridine prodrugs, like ticlopidine, clopidogrel and prasugrel, with the lowest polar surface area (PSA) values, exhibit the largest absorption. A high value of polar surface area (PSA) of cangrelor (255 Å2) results in substantial worsening of the absorption in comparison with thienopyridine drugs. Keywords: antiplatelet agents; molecular structure; solvent effect; pKa; lipophilicity; solubility; absorption; polar surface area 1. Introduction Arterial thrombosis is a common cause of ischemic stroke, myocardial infarction and limb gangrene. This disease is a leading cause of morbidity and mortality in developed countries [1]. Atherothrombosis is essentially a platelet-driven process [2–4]. Pharmaceutical research has focused on the development of antiplatelet drugs targeting key pathways of platelet activation [2–6]. Three classes of antiplatelet agents are currently approved for clinical use: (i) cyclooxygenase-1 (COX-1) inhibitors (aspirin); (ii) adenosine diphosphate (ADP) P2Y12 receptor antagonists (clopidogrel); and (iii) glycoprotein IIb/IIIa inhibitors (abciximab) [2,3]. Besides the cornerstone of antiplatelet therapy, aspirin, also thienopyridines (ticlopidine, clopidogrel, prasugrel), especially in dual therapy with aspirin, have become standard antiplatelet therapies [7]. Owing to known clinical limitations of traditional drugs, like clopidogrel, a number of novel drugs are currently at various stages of clinical development. Three agents, prasugrel, cangrelor and ticagrelor, were already approved for antiplatelet therapy. Int. J. Mol. Sci. 2016, 17, 388; doi:10.3390/ijms17030388 www.mdpi.com/journal/ijms Int. J. Mol. Sci. 2016, 17, 388 2 of 18 Int. J. Mol. Sci. 2016, 17, 388 2 of 17 Ample experimental biological evidence is available for the relationship between the chemical and biologicalAmple experimental properties ofbiological new antiplatelets evidence is avai targetinglable for the relationship P2Y12 platelet between receptor the chemical and their biologicaland biological activity. properties Nevertheless, of new experimental antiplatelets studiestargeting engaged the P2Y1 in2 theplatelet systematic receptor comparisonand their of biological activity. Nevertheless, experimental studies engaged in the systematic comparison of the the physicochemical and pharmacokinetic parameters of these antiplatelet agents are not very well physicochemical and pharmacokinetic parameters of these antiplatelet agents are not very well characterized.characterized. The The absence absence of experimentalof experimental structural structural data data of of new new syntheticsynthetic antiplatelets targeting targeting the P2Y12the platelet P2Y12 platelet receptor receptor presents presents a challenge a challenge to the to applicationthe application of molecularof molecular modeling modeling methods methods with the aimwith to the enhance aim to our enhance understanding our understanding of the subtle of pharmacologicalthe subtle pharmacological effects of these effects antithrombotics. of these antithrombotics.Here, we use density functional methods [8–10] to investigate the structure of six P2Y12 platelet receptor inhibitorsHere, we use (approved density functional drugs ticlopidine, methods [8–10] clopidogrel, to investigate prasugrel, the structure cangrelor of six and P2Y12 ticagrelor) platelet and experimentalreceptor inhibitors drug elinogrel, (approved respectively drugs ticlopidine, (Figure clopidogrel,1). The polarizable prasugrel, cangrelor continuum and methodticagrelor) was and used to accountexperimental for the drug effect elinogrel, of solvent respectively on the (Figure equilibrium 1). The conformation polarizable continuum of these method drugs. was Our used interest to is also focusedaccount for on the the effect evaluation of solvent of on the the important equilibrium physicochemical conformation of these properties drugs. ofOur drugs interest studied is also (pKa, focused on the evaluation of the important physicochemical properties of drugs studied (pKa, lipophilicity, solubility, absorption and polar surface area). The results of the molecular modeling lipophilicity, solubility, absorption and polar surface area). The results of the molecular modeling investigationsinvestigations of these of these antiplatelets antiplatelets were were compared compared to to the the available available experimental experimental structural structural data dataand and discusseddiscussed in relation in relation to theto the present present theories theories ofof actionaction of of these these agents. agents. Cl 1 3 5 2 4 N S Cl Ticlopidine O 9 7 1 3 5 8 2 4 6 N OH SH Ticlopidine metabolite DBMET01202 CH3 7 O O Cl 6 1 3 5 2 4 N S Clopidogrel CH3 7 O O 11 Cl 6 O 1 3 5 9 10 2 4 8 N OH SH Clopidogrel metabolite H4 Figure 1. Cont. Int. J. Mol. Sci. 2016, 17, 388 3 of 18 Int. J. Mol. Sci. 2016, 17, 388 3 of 17 1 F 2 3 5 8 4 N 6 9 10 O 11 CH O 7 3 S 12 Prasugrel O 1 F 2 5 3 b 9 4 8 10 OH 6 N O 7 a SH O 11 Prasugrel metabolite R-138727 CH3 H HN N O 1 2 4 N 3 5 F O O O 6 12 O 7 9 S 11 S N 8 N 10 Cl Elinogrel Hb Ha 5 4 3NH 2 1 N F N N 6 OH F N7 N 8 H C S 9 3 OH O Ticagrelor OH 5 4 3 NH 2 N 1 F N N6 OH F N 7 N 8 H C S 9 3 OH Ticagrelor metabolite AR-C124910XX OH Figure 1. Cont. Int. J. Mol. Sci. 2016, 17, 388 4 of 18 Int. J. Mol. Sci. 2016, 17, 388 4 of 17 H3C S 4 3 NH 2 N 1 5 N 6 N 7 8 S N O 9 10 O OH HO O O P OH Cangrelor 11 12 OH O P HO CF3 OH 13 14 15 P16 Cl Cl O 17 FigureFigure 1. Structure 1. Structure and and atom atomlabeling labeling in th thee antiplatelet antiplatelet drugs drugs investigated. investigated. 2. Results2. Results and and Discussion Discussion 2.1. Molecular Structures 2.1. Molecular Structures The prediction of the equilibrium structure of drugs is one of the main goals of modern quantum Thechemistry prediction [11]. Quantum of theequilibrium chemical calculations structure are of now drugs an iseffective one of tool the in main drug goals design of and modern quantumdevelopment. chemistry They [11 are]. Quantum used for an chemical accurate calculationsdetermination are of molecular now an effective structures tool and in properties drug design for and development.use in a wide They variety are usedof computer-aided for an accurate drug determination design (CADD) of studies molecular [12,13]. structures and properties for use in a wideThe relative variety orientations of computer-aided of important drug functional design (CADD)groups of studies the antiplatelet [12,13]. agents studied here Theare defined relative by orientations dihedral angles of important (α, β, γ, functionalδ, ε, ζ, η, ν, groups θ, μ and of ξ the; Table antiplatelet 1). The overall agents studiedview of the here are definedoptimized by dihedral geometry angles of the (α ,drugsβ, γ, isδ, presented", ζ, η, ν, inθ, thµeand Supplementaryξ; Table1). TheMaterial. overall The view relevant of the values optimized of geometrythe torsion of the angles drugs are is presentedreported in inTable the 1. Supplementary In all cases studied, Material. except The for cangrelor relevant valuestetrasodium, of the the torsion angleswater are reportedhad only a in small Table effect1. In on all the cases equilibrium studied, ge exceptometry for of cangrelor the antiplatel tetrasodium,et agents investigated the water had (Table 1). Table 2 shows the solvation energies obtained from combining calculations performed in a only a small effect on the equilibrium geometry of the antiplatelet agents investigated (Table1). Table2 vacuum with calculations based on the solvation method used. B3LYP/6-311++G (d,p) calculations in showswater the yielded solvation conformers energies not obtained significantly from different combining