- 1Department of Physics and Astronomy, Optical Sciences Centre, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, VIC, Australia
- 2Department of Chemistry and Biotechnology, Centre for Translatonal Atomaterials, Faculty of Science, Engineering and Technology, Swinburne University of Technology, Melbourne, VIC, Australia
Dacomitinib (PF-00299804) was recently approved by the Food and Drug Administration (FDA) as a tyrosine kinase inhibitor (TKI). Unfortunately, side effects and disease resistance eventually result from its use. Off-target effects in some kinase inhibitors have arisen from drug conformational plasticity; however, the conformational states of Dacomitinib in solution are presently unknown. To fill this gap, we have used computational chemistry to explore optimized molecular geometry, properties, and ultraviolet-visible (UV-Vis) absorption spectra of Dacomitinib in dimethyl sulfoxide (DMSO) solution. Potential energy scans led to the discovery of two planar and two twisted conformers of Dacomitinib. The simulated UV-Vis spectral signatures of the planar conformers reproduced the two experimental spectral bands at 275 and 343 nm in solution. It was further discovered that Dacomitinib forms conformers through its three flexible linkers of two C–NH–C bridges, which control the orientations of the 3-chloro-4-fluoroaniline ring (Ring C) and the quinazoline ring (Rings A and B) and the 4-piperidin-1-yl-buten-2-nal side chain, and one C–O–C local bridge which controls the methoxy group locally. When in isolation, these flexible linkers form close hexagon and pentagon loops through strong intramolecular hydrogen bonding so that the “planar” conformers Daco-P1 and Daco-P2 are more stable in isolation. Such flexibility of the ligand and its ability to dock and bind with protein also depend on their interaction with the environment, in addition to their energy and spectra in isolation. However, an accurate quantum mechanical study on drug/ligand conformers in isolation provides necessary reference information for the ability to form a complex with proteins.
Introduction
Lung cancer is the most lethal malignancy among all kinds of cancer worldwide, and non-small-cell lung carcinoma (NSCLC) is counted as the most prevalent (Abdelhameed et al., 2019; Lau et al., 2019). Overexpression of epidermal growth factor receptors (EGFR) and abnormal signaling through the human epidermal growth factor receptor family (HERs) are considered as contributors to NSCLC (Huang et al., 2009; Costa and Kobayashi, 2015). Quinazoline derivatives designed to block HER activation are popular antitumor agents (Selvam and Kumar, 2011; Shagufta, 2017; Qiu et al., 2019; Solyanik, 2019). Dacomitinib (Reed and Smaill, 2016) (trade name: Vizimpro), an anilinoquinazoline derivative, recently received FDA approval (Shirley, 2018; Roskoski, 2019) as a tyrosine kinase inhibitor (TKI) for the treatment of locally advanced or metastatic NSCLC, as it irreversibly inhibits three out of four HER families: HER-1, HER-2, and HER-4 tyrosine kinases (Gonzales et al., 2008). This highly selective, small-molecule TKI of the epidermal growth factor receptor (EGFR) (Bello et al., 2013) has the aniline–quinazoline ring structures with an acrylamide branch known as the Michael acceptor. This Michael acceptor, a chemically reactive electrophilic warhead, targets a cysteine nucleophile in the adenine pocket to form an irreversible covalent adduct (Gajiwala et al., 2013).
Despite high controversy (Caruso, 2018; Passaro and de Marinis, 2018) and earlier recommendations for use in second-line treatment (Gonzales et al., 2008; Zugazagoitia et al., 2017), dacomitinib was recently suggested as a first-line clinical therapy worldwide (Wu et al., 2017; Decoster et al., 2018; Shirley, 2018; Abdelhameed et al., 2019; Chustecka, 2019; Roeper and Griesinger, 2019). In addition, this more potent TKI could possibly find other applications such as suppressing brain tumor (Chen et al., 2019), chemoresistant ovarian cancer (Momeny et al., 2017), and breast cancer (Kalous et al., 2012) in future.
Compared to the first-generation anticancer agents, Dacomitinib has prolonged progression-free survival and overcame the resistance from the EGFR-muted NSCLC (Gonzales et al., 2008; Lau et al., 2019) but at the expense of increased cytotoxicity (Zugazagoitia et al., 2017). Still, NSCLC eventually acquired resistance to Dacomitinib by either EGFR T790M or C797S mutations (Kobayashi et al., 2018). Commonly found adverse reactions of this EGFR antagonist are diarrhea, rash, nail changes (paronychia and onycholysis), fatigue, stomatitis, loss of appetite, conjunctivitis, weight loss, alopecia, cough, itching, nausea, etc. (Shirley, 2018; Chustecka, 2019). These side effects may possibly result from its two reactive metabolites (Attwa et al., 2018). Alternatively, unexpected off-target effects can also arise from drug conformational plasticity (Hantschel, 2015), i.e., distinct drug conformers, which bind on-target and off-target proteins. It is therefore important to determine the relevant conformational states in solution.
As Dacomitinib is a very newly FDA-approved drug (Shirley, 2018; Roskoski, 2019), to the best of our knowledge, there are no reported detailed electronic structural studies of the drug and its conformation in solution, although the crystal structure of the drug is available and drug conformation in the drug–protein complex (Gajiwala et al., 2013) was also reported. Without detailed molecular structural information and knowledge about how the drug responds to the environment, it is difficult to understand the mechanism of the ligand–protein interaction. The mechanism of the ligand–protein interaction is important for the development of new drugs. For example, Dacomitinib (PF00299804) is considered in the same class of a NSCLC drug Gefitinib (ZD1839) as shown in Scheme 1, which differs only in the side chain. However, Dacomitinib (PF00299804) and Gefitinib (DZ1839) exhibit significant drug potent activities. The former effectively inhibits the in vitro kinase activity of wild-type EGFR (IC50 = 6 nM) (Engelman et al., 2007), whereas the latter is an inhibitor that specifically binds and inhibits the EGFR tyrosine kinase, with the IC50 value of 2–37 nM in NR6wtEGFR cells (Pedersen et al., 2005).
In order to understand how the chemical structural changes in Gefitinib and Dacomitinib affect the drug potency and ability to bind with proteins, conformation(s) in solution is required. The present study provides the first accurate quantum mechanical electronic structures and conformations of Dacomitinib in solution. In this study, the ultraviolet-visible (UV-Vis) absorption of Dacomitinib in dimethyl-sulfoxide (DMSO) solution and in silico, for the first time, is investigated applying the time-dependent density functional theory (TD-DFT). The experimental measurements together with computational calculations elucidate the interrelationship of the molecular structure and spectral properties of Dacomitinib.
Materials and Methods
Experimental Details
Dacomitinib and DMSO were purchased from Focus Bioscience Pty Ltd. (MedChemExpress) and Sigma-Aldrich Pty Ltd., respectively. All UV-Vis absorption experiments were carried out on the same day, using a pair of matched quartz cuvettes (purchased from Stama Pty Ltd.) with a 1-cm path length. The electronic absorption spectra (absorbance of λmax < 0.1) of the Dacomitinib solution in DMSO (5–25 μL) in the range of 250–400 nm was recorded on a Perkin-Elmer LAMBDA 1050 UV/Vis/NIR spectrophotometer at room temperature (293 K).
Computational Details
The IUPAC name of Dacomitinib (C24H25ClFN5O2) is (2E)-N-{4-[(3-chloro-4-fluorophenyl)amino]-7-methoxyquinazolin-6-yl}-4-(piperidin-1-yl)but-2-en amide. The initial 3D geometry of the crystal structure of Dacomitinib was downloaded from the PubChem website (https://pubchem.ncbi.nlm.nih.gov/compound/Dacomitinib#section=3D-Conformer) with a ligand ID code of 1C9 in Figure 1. Optimization of this moiety in DMSO was performed using the density-functional theory (DFT)-based Becke three-parameter Lee–Yang–Parr hybrid functional (B3LYP) (Becke, 1993) in combination with the 6-311+G(d,p) basis set and the conductor-like polarizable continuum model (CPCM) (Cossi et al., 2003). Potential energy scan (PES) was performed by relaxed scan through the C(17)-N(5) bond using B3LYP/6-311G. The resulting local minima geometries were further optimized at B3LYP/6-311+G(d,p) in the DMSO solvent. Absorption UV-Vis spectra in DMSO solution were calculated using the time-dependent (TD)-DFT method for the lowest 40 excited states of singlet–singlet transitions using the same model. The UV-Vis spectrum of the most stable structure of the drug Daco-P1 in DMSO solution was calculated using both B3LYP/6-311+G(d,p) and CAM-B3LYP/6-311+G(d,p) methods, respectively. The method has been established in our previous studies of this class of TKIs such as AG-1478 (Khattab et al., 2016) and SKF86002 (Van Dongen et al., 2017). The excess orbital energy spectrum (EOES) proposed by Islam and Wang (2015) was employed to examine the core electron energies between the Dacomitinib conformers. The EOES recognizes which orbital sites change in energy in response to structural changes through conformation (Wang and Chatterjee, 2017); EOES was performed by analyzing the orbital energies of the higher-energy conformers with respect to those of the global minimum structure (Khattab et al., 2016). All calculations were carried out using the Gaussian 16 Revision A.03 computational chemistry package (Frisch, 2016).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/538256/fchem-08-00596-HTML/image_m/fchem-08-00596-g001.jpg)
Figure 1. Nomenclature and chemical structure of Dacomitinib ((2E)-N-(4-((3-chloro-4-fluorophenyl)amino)-7-methoxyquinazolin-6-yl)-4-piperidin-1-ylbut-2- enamide). The Arrows indicate the most likely rotations. The crystal structure of (1C9) (Lau et al., 2019) and the co-crystallized structures of Dacomitinib from the Protein Data Bank (PDB), 4I23 and 4I24 (Abdelhameed et al., 2019). The dihedral angles of C(12)-C(11)-O(3)-C(24) and C(19)-C(18)-N(5)-C(17) define the flexibility of the ligand structures, which vary according to the conditions of crystallization of the drug (ligand).
Results and Discussion
Dacomitinib Conformers in Crystal Phase and in Isolation
The chemical structure and nomenclature of Dacomitinib (C24H25ClFN5O2, CAS 1110813-31-4) are given in Figure 1. A single-bond NH linker, C(17)-N(5)(H)-C(18), connects the 3-chloro-4-fluoroaniline ring (Ring C) and the quinazoline ring (Rings A and B). The initial crystal structure of Dacomitinib was downloaded from the PubChem website (https://pubchem.ncbi.nlm.nih.gov/compound/Dacomitinib#section=3D-Conformer) with a ligand ID code of 1C9. The Dacomitinib co-crystal structures by Gajiwala et al. (2013) with IDs of 4I23 and 4I24 were downloaded from the Protein Data Bank (PDB) for more structural information. The three crystal structures of Dacomitinib are given in Figure 1. The co-crystal ligand structures (4I23 and 4I24) are not the complete chemical structure of Dacomitinib but are the core structure of 3-chloro-4-fluoroaniline ring and the quinazoline ring present in 4I24 (without piperidin-1-yl) and in 4I23 (without 4-piperidin-1-yl-buten-2-nal). As can be seen in this figure, the conformations of the Dacomitinib ligand depend on both the crystallization conditions and the complexes that the ligand forms with the proteins (polymers). The most apparent conformational changes are in three major bridges of C(10)-N(6)H-C(25), C(11)-O(3)-C(24) and C(17)-N(5)-C(18), which control the flexibility of the ligand (drug molecule) through intramolecular hydrogen bonding (HB). The two N-bridges control the relevant orientations of the 3-4-piperidin-1-yl-buten-2-nal tail and the chloro-4-fluoroaniline ring (C) with respect to the quinazoline ring (Rings A and B). The O-bridge provides small tunes to the methoxy group (-OCH3) and the quinazoline ring (Rings A and B) locally.
The core structure of Dacomitinib is the two-ring backbone system [3-chloro-4-fluoroaniline ring (Ring C)] and the quinazoline ring (Rings A and B) as shown in Scheme 1. Variations of the side R-chain and results in Gefitinib and Dacomitinib are shown in Scheme 1. The crystal structures (1C9, 4I24, and 4I23) reported in Figure 1 have different side R-chains, and only the 1C9 structure contains the full Dacomitinib ligand structure (piperidin-1-yl is missing in 4I24 and 4-piperidin-1-yl-buten-2-nal is missing in 4I23). Table 1 collects the selected structural properties of the ligand in Figure 1 from different sources, that is, the same ligand when being confined in crystal or co-crystal complexes is different from that in isolation (solution may be considered in this case). Moreover, under different crystallization conditions, the ligand results in different conformations with respect to the flexible regions.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/538256/fchem-08-00596-HTML/image_m/fchem-08-00596-t001.jpg)
Table 1. Selected geometric properties of crystal structure (literature) and calculated Dacomitinib conformers obtained.
Full geometry optimization was performed to locate a planar minimum energy structure of Dacomitinib (Daco-P1). Here, the “planar” structure Daco-P1 only means that the dihedral angle of the plane formed by Ring C (3-chloro-4-fluoroaniline ring) and Rings A and B (the quinazoline ring) is very small when the drug molecule is in isolation. Based on the obtained structure (Daco-P1), three more local minima (conformers) were found by a potential energy scan (PES) calculation through rotation of the C(17)-N(5) bond (i.e., dihedral angle ∠N(7)-C(17)-N(5)-H), which evidently alter neither any bond lengths nor most of the bond angles outside of the selected parameters. The perimeters of the three aromatic rings in four conformers, for example, are altered by <0.01 Å. Figure 2 presents the PES and the obtained local minima and their structures, in which the differences of the bond angles related to the NH linker are negligible too, except that the NH linker local to dihedral angles altered considerably [all include N(5)]. In addition, our previous study on a TKI drug AG-1478 (Khattab et al., 2016) in the same class discovered that rotation of the C(17)-N(5) bond produced possibly potent conformers of the drug. Therefore, the present study concentrates on the rotation of the C(17)-N(5) bond for conformer searches. As an illustration, the dihedral angle of ∠N(7)-C(17)-N(5)-H showed substantial conformational alterations from the planar structure of Daco-P1 to the twisted structure of Daco-T1 by 160°. The obtained four stable local minima (conformers) are in two pairs of “planar” Daco-P1 and Daco-P2 and twisted Daco-T1 and Daco-T2 conformers. A small energy barrier of ~1.5 kcal·mol−1 between Daco-T1 and Daco-T2 is as shown in Figure 2. However, for Dacomitinib conformers (Daco-P1 and Daco-P2) to become more stable, the twisted conformers need to cross a higher-energy barrier of ~5 kcal·mol−1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/538256/fchem-08-00596-HTML/image_m/fchem-08-00596-g002.jpg)
Figure 2. Potential energy surface (PES) scan of Dacomitinib conformers in DMSO through local rotation of the C(17)-N(5) bond (dihedral angle ∠N(7)-C(17)-N(5)-H) using the B3LYP/6-311G level of theory. Pointed arrows at a particular point on PES refer to the corresponding 3D structures of Dacomitinib. In the 3D structures of Dacomitinib conformers, atoms of carbon, hydrogen, nitrogen, oxygen, chlorine, and fluorine are denoted by gray-, white-, blue-, red-, green-, and cyan-colored atoms, respectively.
Figure 3 displays the obtained local minimum structures of Dacomitinib obtained from calculations in isolation. Daco-P1 is evidently the most energetically stable, i.e., global, minimum structure. The other local minimum structures such as Daco-P2, Daco-T1, and Daco-T2 are 0.51, 4.44, and 4.60 kcal mol−1, respectively, above the Dac-P1 conformer. Table 1 compares some selective structural parameters of the crystal structures (1C9, 4I23, and 4I24) in Figure 1 and the calculated structures of Dacomitinib in isolation in Figure 3. Table S1 in Supplementary Materials contains the Cartesian coordinates of the four conformers. Figure S1 provides the 3D crystal structure (1C9) and four theoretically obtained Dacomitinib conformers.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/538256/fchem-08-00596-HTML/image_m/fchem-08-00596-g003.jpg)
Figure 3. The optimized three-dimensional (3D) structures of the Dacomitinib conformers obtained from the search around the crystal structures are given along with their relative energies. In the 3D structures of Dacomitinib conformers, atoms of carbon, hydrogen, nitrogen, oxygen, chlorine, and fluorine are denoted by gray-, white-, blue-, red-, green-, and cyan-colored atoms, respectively. The energy difference between the most table minimum energy structure (Daco-P1) and the crystal structure (1C9) is 91.15 kcal·mol−1.
As can be seen in Table 1, almost all geometric parameters except for the one highlighted (i.e., the flexible regions and their HB) are very close in values. In particular, the Daco-P1 and Daco-P2 structures exhibit large similarities to the 1C9 crystal structure and the root mean square deviation (RMSD) is smaller than 0.62, which is nearly half the RMSD of Daco-T1 (1.144) and Daco-T2 (1.631). In addition, the dipole moment of the crystal structure is calculated as 11.967 D, which is in close agreement to 12.721 D for Daco-P1 (12.721 D) and Daco-P2 (11.132 D). When in isolation, the “planar” Daco-P1 and Daco-P2 are more stable than the twist conformer pairs. The planar drug (ligand) conformers in dilute solution engage with stronger intramolecular HBs (see Table 1), which contribute to a reduction in the energy. That is, the atoms in these regions are able to form pentagon or hexagon loops through intramolecular hydrogen bonding (HB). For example, the H at the left-hand side forms intramolecular HB of O or N at the right-hand side to close the loop of Hex-1 of H-C(15)-C(10)-N(6)-C(25)-O(4), Hex-2 of H-C(19)-C(18)-N(5)-C(17)-N(7), and Pent H-N(6)-C(10)-C(11)-O(3). The planar structures Daco-P1 and Daco-P2 enhances the intramolecular HB and therefore stabilizes the structures. The three HB of C(15)H···O(4), N(6)H···O(3), and C(19)H···N(7) of Daco-P1 and Daco-P1 are stronger than Daco-T1 and Daco-T2 as their H···O(N) distances are shorter.
The crystal structures are confined structures, which most likely influence the bond angles. The most significant differences of the crystal structures of Dacomitinib are two of the three flexible regions, Hex-2 of H-C(19)-C(18)-N(5)-C(17)-N(7) and Pent H-N(6)-C(10)-C(11)-O(3). In the pair of related dihedral angles C(12)-C(11)-O(3)-C(24) and C(19)-C(18)-N(5)-C(17), the former controls the relevant orientation of the methoxy group to the quinazoline ring and the latter controls the 3-chloro-4-fluoroaniline ring (Ring C) and the quinazoline ring (Rings A and B). As indicated in Table 1, this pair of dihedral angles deviates from planarity (0° or 180°) in the crystal structures (1C9, 4I24, and 4I23) due to confinement in a crystal or a complex. When in isolation, this pair of dihedral angles adopt planar configuration (0° or 180°) (except for the twist conformers) in order to form strong intramolecular HB to reduce the energy as much as possible (Steiner, 2002). As a result, the free Daco-PI conformer in the DMSO solution is more stable than the crystal structure in the same solution by ~9 kcal·mol−1.
Spectral Signatures of Dacomitinib Conformers in the UV-Vis Absorption Spectra
Figure 4 compares the calculated UV-Vis spectra of the Daco conformers with the measurement in DMSO solution. Our previous studies on this class of drugs, such as AG-1478 (Khattab et al., 2016), revealed that the TD-DFT method with the B3LYP/6-311+G(d,p) method produced UV-Vis spectra of the drug in excellent agreement with the measurements, which is in agreement with the present study. The calculated UV-Vis spectra of Daco-P1 using B3LYP/6-311+G(d,p) reproduce the maximum UV-Vis transition of the drug at 345.18 nm, which is in excellent agreement with the measurement of 343 nm in the same solvent of DMSO. While the CAM-B3LYP/6-311+G(d,p) method produces the same transition at 302.74 nm which is approximately a −40-nm shift from the measurement. As a result, B3LYP/6-311+G(d,p) is employed to produce the UV-Vis spectra for other Daco drugs. The resultant UV-Vis spectra are compared in Figure S2 of the Supplementary Materials.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/538256/fchem-08-00596-HTML/image_m/fchem-08-00596-g004.jpg)
Figure 4. Comparison of the measured UV-Vis absorption spectra of Dacomitinib in the DMSO spectrum (solid black line) with that calculated (colored) in the same solvent. All spectra including both the calculated and the measured are normalized to their own individual maximum absorbances. Please note that the wavelength cutoff due to solvent is 250 nm, and the maximum absorbances of some Dacomitinib conformers lie in the region smaller than 250 nm.
In Figure 4, the experimental UV-Vis spectrum in the region of 250–450 nm of the Dacomitinib drug in DMSO is compared with the calculated UV-Vis absorption spectra of the four conformers of Dacomitinib as well as the crystal structure (1C9 without further optimization) under the same conditions. The experimental UV-Vis spectrum exhibits two major bands, Band A at 261 nm and Band B at 343 nm; both bands (black) have shoulders at 275 nm of Band A and shoulders at 320 and 358 nm at both sides of Band B. UV-Vis spectra, measured at different Dacomitinib concentrations, are given in Figure S3 of the Supplementary Materials.
Inspection of the simulated UV-Vis spectra (see Figure 4) revealed that the planar conformers (i.e., Daco-P1 and Daco-P2) are more likely the dominant Dacomitinib conformations in solution, over the twisted conformers of Daco-T1 and Daco-T2. The simulated UV-Vis spectra of the planar pair of Daco-P1 and Daco-P2 agree better with the measured spectrum, exhibiting two bands at about 280 and 345 nm; the UV-Vis spectrum of the crystal structure (1C9) is the next with two major bands at 262.22 and 390.1 nm. The major HOMO–LUMO transition at 390.1 nm shifts 47 nm from the measured 343 nm. The twisted dacomitinib conformers display one major band at ~283.27 (278.75) nm for Daco-T1 (Daco-T2) and a shoulder band at ~344.03 nm (336.75 nm) for Daco-T1 (Daco-T2). Although the shoulder band above 300 nm stems from the HOMO–LUMO transition, the strength of the band is not significant.
The ratio of the peak intensities at 275 and 343 nm of the experimental UV-Vis spectrum (peak ratio of 1.44) has a close match with the calculated spectra of the planar conformers (peak ratio of 1.29 and 1.33 for Daco-P1 and Daco-P2, respectively), rather than the crystal structure of ~3.0, and the twisted conformers with a peak ratio of 2.99 and 2.64 for Daco-T1 and Daco-T2, respectively. This suggests that in solution the planar conformers are predominant, which is supported by the calculated Boltzmann distribution at the experimental temperature of 293.15 K. That is, the drug in isolation is predominately populated by Daco-P1 (51.08%) and Daco-P2 (48.91%) but the twisted conformers Daco-T2 and Daco-T1 can be neglected. When in isolation including in solution, the Dacomitinib ligand does not take the same conformation as it is in the crystal phase.
Figure 5 shows the frontier orbital diagrams ranging from HOMO-3 to LUMO+3 of the low-energy-lying conformers, Daco-P1, Daco-P2, Daco-T1, and Daco-T2, calculated in DMSO solution, along with their corresponding orbital density distributions of HOMO and LUMO. Considering subtle differences in energy gap, Daco-P1 is the lowest-energy configuration and hence most stable structure. The conformers have a small energy difference in HOMO → LUMO energy gap (4.12 eV, 4.12 eV, 4.19 eV, and 4.25 eV for Daco-P1, Daco-P2, Daco-T1, and Daco-T2, respectively). All conformers exhibit a C1 point group symmetry including the “planar” Daco-P1 and Daco-P2. As a result, the HOMO–LUMO energy gap or the H → L transitions will be dominant (>96% see Table S2). The HOMO → LUMO gap contributes the major band of UV-Vis (λ > 300 nm in Figure 3), and the planar conformers have slightly smaller HOMO → LUMO gaps and therefore larger λ in the UV-Vis spectra compared with the twisted ones (see Figure 3). The probability of the transitions of electrons increases with the decrease in the energy gap (Δε) between the occupied orbitals (HOMO-n) and the virtual orbitals (LUMO+m).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/538256/fchem-08-00596-HTML/image_m/fchem-08-00596-g005.jpg)
Figure 5. Energy gaps (in eV) among different energy levels (from HOMO-2 to LUMO+2) and HOMO–LUMO orbital charge densities for Dacomitinib conformers in DMSO at the B3LYP/6-311G(d,p) level of theory. The positive and negative electron densities are shown in green and red, respectively.
Figure 5 also provides the orbitals of HOMO and LUMO for all conformers obtained in isolation. Two apparent features of the frontier orbitals are noticeable: firstly, the orbital electron density of HOMOs and LUMOs concentrates on the aromatic quinazoline and aniline rings through the –N(5)H– linker rather than on the piperidine side chain. The planar structure of Daco-P1 and Daco-P2 helps delocalize the HOMO electrons in the conjugate system, whereas the HOMOs of Daco-T1 and Daco-T2 push the HOMO electrons to diffuse into the piperidine side chain through the other –N(6)H– bridge. Secondly, the LUMOs of the 3-chloro-4-fluoroaniline ring (Ring C) of Dacomitinib are less populated and the electrons diffuse through the –N(6)H– linker toward the piperidine side chain in the Dacomitinib conformers. Folding the 3-chloro-4-fluoroaniline ring (Ring C) over toward the quinazoline ring (Rings A and B) in Daco-T1 and Daco-T2 helps the LUMO electron density to concentrate in the –N(6)H– bridge region.
Electronic Structures of the Dacomitinib Conformers
Molecular electrostatic potentials (MEP) of a drug determine approximately the shape of the drug. Shape complementarity is important in drug–target interactions, and so a less stable conformer of a drug can sometimes be more potent if it has the right shape (Chen and Wang, 2009). Hence, it is crucial to realize the differences in electronic structures and shapes among the conformers, which may explain their characteristics like drug potency. In Figure 6, the molecular electrostatic potential (MEP) portrays the distribution of electron density of four conformers along with the scale of the color spectrum having minimum (negative extreme) and maximum (positive extreme) electrostatic potential energy ends marked by red and blue colors, respectively. Overall, in the four conformers, multiple electron donating spots (red) can be identified, which are dominated by the large electronegative atoms such as F, Cl, O, and N as shown in Figure 6. For example, the electron density is withdrawn from the methoxy (O(3)-C(24)H3) group by the quinazoline ring indicated by a light blue color. The electron density of piperidine containing the crotonamide Michael acceptor on C(10) in general is neither high nor low except the amide group (HN(6)-C(25)O(4)). In this amide group, high electron density is observed over a lone-pair electron containing O(4) and the electron less populated in N(6), which are represented as intense red protrusion and deep blue notch, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/538256/fchem-08-00596-HTML/image_m/fchem-08-00596-g006.jpg)
Figure 6. Comparative molecular electrostatic potential (MEP) maps for four conformers of Dacomitinib: (A) Daco-P1, (B) Daco-P2, (C) Daco-T1, and (D) Daco-T2 in DMSO solution.
The MEPs of planar Daco-P1 and Daco-P2 depict that they both have similar electron density distributions other than the orientation of the yellowish protrusion from the chlorine atom. In contrast to the planar conformers, the twisted ones (Daco-T1 and Daco-T2) have a significant difference in electron density indicated by comparative intense red and blue colors because of their packed structure. A drop in electron density over the N(6) atom and a significant increase in electron localization over the regions of O(4), N(7), N(8), and F(1) atoms are noticeable from Daco-P1 to Daco-T2 [refer to Figure 6 from left (a) to right (b)], respectively. Hence, the planar structures have more electron delocalization and thus are more stable compared to the twisted structures. The shape is the most crucial alteration among the four conformers of Dacomitinib—when the 3-chloro-4-fluoroaniline ring rotates over the single C(17)-N(5) bond and folds over the quinazoline ring, which makes the nitrogen atoms (i.e., N(7) and N(8)) of the core ring in twisted structures more reactive to form a bond in space, as the MEP map indicates the electrostatic potential of the molecular surface, and it can help in the prediction of intermolecular electrostatic bond formation such as hydrogen bond (an electrostatic force) (Andree and Aakeröy, 2018). The MEP of Daco-P1 and Daco-P2 clearly indicates a possibility of intermolecular hydrogen bond formation between the N(7) atom (high negative potential reflected as reddish MEP surface) and closely aligned hydrogen atom at C(19) (relatively positive potential reflected as a greenish MEP surface).
Another electronic property which helps to identify specific atoms of the drug conformers is the identification of responsible atomic sites of the drug conformers. The experimental UV-Vis spectrum implies that more than one such Dacomitinib conformer may exist in the DMSO solution. Recent X-ray photoemission spectroscopy (XPS) studies have revealed that conformers experience core electron energy changes (Islam et al., 2018). The excess orbital energy spectrum (EOES) (Islam and Wang, 2015) was therefore introduced to reveal the orbital response to the conformational changes from the reference (which is often the global minimum structure) using a simple method. It is quite helpful in the study of AG-1478 drug conformers (Khattab et al., 2016). The EOES in Figure 6 are three spectra of excess orbital energies Δεi = -εigm, in which is the orbital energies of a conformer and εigm is the orbital energies of the reference which is usually the most stable conformer (Daco-P1). In Figure 7, only atom-specific core 1s energies are considered. The gray solid circles are the excess core 1s energies of Daco-P2 with respect to Daco-P1. The red squares represent the excess core 1s energies of Daco-T1 with respect to Daco-P1, and the blue triangle is the excess core 1s energies of Daco-T2 with respect to Daco-P1. If the excess core 1s orbital energies Δεi are outside of the red horizontal dashed lines at a cutoff, for example, ±1.75 kcal mol−1 in Figure 6, the conformational impact is significant on that atomic site. As a result, using the core 1s EOES, one can identify the atomic sites which respond significantly to the conformational changes from the most stable conformer (Daco-P1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/538256/fchem-08-00596-HTML/image_m/fchem-08-00596-g007.jpg)
Figure 7. Excess orbital energy spectrum (EOES) of Dacomitinib in DMSO at the B3LYP/6-311G** level of theory: Δεi = -εigm, where = ith orbital energy of a local minimum conformer (i.e., Daco-P2, Daco-T1, or Daco-T2) and εigm = ith orbital energy of a global minimum conformer, Daco-P1. The X axis is the orbital order from the innermost orbital as Cl(1) to the HOMO as MO(123). The valence orbitals from MO(34) to MO(122) are omitted for simplicity.
In its ground electronic state, Dacomitinib (C24H25ClFN5O2) has a total of 123 doubly occupied orbitals, which consists of 33 1s core orbitals. While most of the core 1s EOES (Δεi = -εigm) of Dacomitinib conformers in Figure 7 distribute within the red horizontal dashed lines at ±1.75 kcal mol−1 (i.e., a cutoff), a few points, however, locate outside of the two parallel red dash lines. This suggests that most of the atoms in conformer Daco-P2 do not change their energies apparently (those within the parallel red dash lines), and a few atoms, which are local to the particular bond rotation and breaking or formation of hydrogen bonding in conformer Daco-P2, respond to the conformational changes significantly. As shown by the gray solid circles in Figure 6, the most significant core 1s electron-energy changes between Daco-P2 with respect to Daco-P1 are C(19)1s and C(23)1s, as this pair of carbons, C(19) and C(23), switch their roles in Daco-P1 and in Daco-P2. In Daco-P1, atom C(19) is involved in hydrogen bonding C(19) ···HN(7), and atom C(23) does not. However, in Daco-P2, the hydrogen bonding becomes C(23) ···HN(7), and atom C(19) is not involved in hydrogen bonding.
The EOES shows that in the planar conformer pair of Daco-P1 and Daco-P2, only a small group of atoms associate with apparent core 1s energy changes. In the EOES of the Daco-T1 and Daco-P1 pair (red squares in Figure 7), more atoms are disturbed with significant changes. The most significant atomic 1s energy changes (outside of the parallel red dash line cutoff) are Cl(1)1s, F(2)1s, N(5)1s, C(10)1s, C(15)1s, C(17)1s, C(18)1s, C(20)1s, and C(21)1s. However, the core 1s electron energy changes of Daco-T2 with respect to Daco-P1 (blue triangles) are Cl(1)1s, F(2)1s, N(5)1s, C(10)1s, C(14)1s, C(15)1s, C(17)1s, C(18)1s, C(20)1s, and C(21)1s. All of the atomic sites affected significantly due to Daco-T1 and Daco-T2 with respect to Daco-P1 are local to either the halogen atoms, C(20)-Cl and C(21)-F, or the local atomic site strained with the C(17)-N(5) bond rotation, that is, the chain formed by C(15)-C(14)-C(17)-N(5)-C(18).
In solution, evidently, multiple conformers of Dacomitinib present and planar conformers are predominant. However, none of our conformers in solution are identical with the previously reported crystal structure (Gajiwala et al., 2013) of Dacomitinib in the complex with its cognate kinase domain. Hence, it is possible that the flexible Dacomitinib changes its conformation when binding with the proteins or confining in crystals, in addition to different phases of the compound.
Conclusion
In this work, we have used the TD-DFT method to explore the molecular geometry, structural features, molecular properties, and UV-Vis absorption spectra of TKI Dacomitinib in the DMSO solution. Potential energy scans locate two stable planar conformers (i.e., Daco-P1 and Daco-P2) and two less stable twisted conformers (i.e., Daco-T1 and Daco-T2), which all possess lower total energy in isolation than the crystal Dacomitinib structure by over 8 kcal·mol−1. The Daco-P1 and Daco-P2 conformers may engage with intramolecular hydrogen bonding, which are C(19)-H…N(7), and C(23)-H…N(7), respectively, forming a hexagon and pentagon structures to lock the flexibility, hereby contributing to stabilization of the planar conformation. The simulated UV-Vis absorption spectrum of Dacomitinib generated from the contribution of the planar conformers closely reproduces the two measured spectral bands at 275 and 343 nm. In particular, the core 1s energy EOES reveals the atomic site-specific changes of the conformers with respect to Daco-P1. The identified local atomic sites between other conformers and Daco-P1 and Daco-P2 may contribute to the drug potency study. This is because the drug potency significantly depends on the shape/conformation of a drug molecule, in which the flexibility of the ligand and its ability to dock and bind with protein are also important, in addition to their energy. However, an accurate quantum mechanical study on drug/ligand conformers in isolation provides necessary reference information for the ability to form complexes with proteins (Yuan et al., 2017). The crystal structure clearly shows large similarities with Daco-P, the only major difference being the twisting around the N5-C17 bond, which is line with our previous study on a related TKI (Khattab et al., 2016). One potential application of this investigation is to extend the study to identification of the conformer(s) when binding with proteins using molecular dynamic studies, which may contribute to reveal the side effects of the Dacomitinib class of drugs.
Data Availability Statement
All datasets generated for this study are included in the article/Supplementary Material.
Author Contributions
MK conducted the experiments, carried out the computational calculations, and writing the manuscript. FB assisted in computational calculations. AC and FW designed and supervised the research. The corresponding authors edited and revised the manuscript.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
MK acknowledges the Swinburne University Postgraduate Research Award (SUPRA). FB acknowledges the Australian Government Research Training Program Scheme (RTPS) Award. All calculations were performed with thanks to the National Computer Infrastructure (NCI) in addition to the Swinburne University Supercomputer Facilities.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fchem.2020.00596/full#supplementary-material
References
Abdelhameed, A. S., Kadi, A. A., Attwa, M. W., and AlRabiah, H. (2019). Validated LC-MS/MS assay for quantification of the newly approved tyrosine kinase inhibitor, dacomitinib, and application to investigating its metabolic stability. PLoS ONE 14:e0214598. doi: 10.1371/journal.pone.0214598
Andree, S. N. L., and Aakeröy, C. B. (2018). Molecular electrostatic potentials as a quantitative measure of hydrogen bonding preferences in solution. Supramol. Chem. 30, 455–463. doi: 10.1080/10610278.2017.1418876
Attwa, M. W., Kadi, A. A., and Abdelhameed, A. S. (2018). Characterization of reactive intermediates formation in dacomitinib metabolism and bioactivation pathways elucidation by LC-MS/MS: in vitro phase I metabolic investigation. RSC Adv. 8, 38733–38744. doi: 10.1039/C8RA06709K
Becke, A. D. (1993). Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98, 5648–5652. doi: 10.1063/1.464913
Bello, C. L., Smith, E., Ruiz-Garcia, A., Ni, G., Alvey, C., and Loi, C.-M. (2013). A phase I, open-label, mass balance study of [(14)C] dacomitinib (PF-00299804) in healthy male volunteers. Cancer Chemother. Pharmacol. 72, 379–385. doi: 10.1007/s00280-013-2207-9
Caruso, C. (2018). Dacomitinib approved, but might not be used. Cancer Discov. 8:1500. doi: 10.1158/2159-8290.CD-NB2018-138
Chen, F. F., and Wang, F. (2009). Electronic structure of the azide group in 3c-azido-3c-deoxythymidine (AZT) compared to small azide compounds. Molecules 14, 2656–2668. doi: 10.3390/molecules14072656
Chen, W. S., Hong, L., Wang, F., and Li, J. J. (2019). Investigation of dacomitinib on reducing cell necrosis and enhancing cell apoptosis in C6 glioma rat model by MRI. Biosci. Rep. 39:BSR20190006. doi: 10.1042/BSR20190006
Chustecka, Z. (2019). Dacomitinib for EGFR+ lung cancer gets EU OK. Medscape. Available online at: https://www.medscape.com/viewarticle/908538
Cossi, M., Rega, N., Scalmani, G., and Barone, V. (2003). Energies, structures, and electronic properties of molecules in solution with the C-PCM solvation model. J. Comput. Chem. 24, 669–681. doi: 10.1002/jcc.10189
Costa, D. B., and Kobayashi, S. S. (2015). Whacking a mole-cule: clinical activity and mechanisms of resistance to third generation EGFR inhibitors in EGFR mutated lung cancers with EGFR-T790M. Trans. Lung Cancer Res. 4, 809–815. doi: 10.3978/j.issn.2218-6751.2015.05.05
Decoster, L., Giron, P., Mignon, S., and De Grève, J. (2018). The evolving first-line treatment of advanced non-small cell lung cancer harbouring epidermal growth factor receptor mutations. Transl. Lung Cancer Res. 7, S134–S137. doi: 10.21037/tlcr.2018.03.08
Engelman, J. A., Zejnullahu, K., Gale, C. M., Lifshits, E., Gonzales, A. J., Shimamura, T., et al. (2007). PF00299804, an irreversible pan-ERBB inhibitor, is effective in lung cancer models with EGFR and ERBB2 mutations that are resistant to gefitinib. Cancer Res. 67, 11924–11932. doi: 10.1158/0008-5472.CAN-07-1885
Frisch, M. J., Trucks, G. W., Schlegel, H. B., Scuseria, G. E., Robb, M. A., Cheeseman, J. R., et al. (2016). Gaussian 16, Revision C.01. Wallingford, CT: Gaussian, Inc.
Gajiwala, K. S., Feng, J., Ferre, R., Ryan, K., Brodsky, O., Weinrich, S., et al. (2013). Insights into the aberrant activity of mutant EGFR kinase domain and drug recognition. Structure 21, 209–219. doi: 10.1016/j.str.2012.11.014
Gonzales, A. J., Hook, K. E., Althaus, I. W., Ellis, P. A., Trachet, E., Delaney, A. M., et al. (2008). Antitumor activity and pharmacokinetic properties of PF-00299804, a second-generation irreversible pan-erbB receptor tyrosine kinase inhibitor. Mol. Cancer. Ther. 7, 1880–1889. doi: 10.1158/1535-7163.MCT-07-2232
Hantschel, O. (2015). Unexpected off-targets and paradoxical pathway activation by kinase inhibitors. ACS Chem. Biol. 10, 234–245. doi: 10.1021/cb500886n
Huang, Z., Brdlik, C., Jin, P., and Shepard, H. M. (2009). A pan-HER approach for cancer therapy: background, current status and future development. Expert Opin. Biol. Ther. 9, 97–110. doi: 10.1517/14712590802630427
Islam, S., Ganesan, A., Auchettl, R., Plekan, O., Acres, R. G., Wang, F., et al. (2018). Electronic structure and intramolecular interactions in three methoxyphenol isomers, J. Chem. Phys. 149:134312. doi: 10.1063/1.5048691
Islam, S., and Wang, F. (2015). The d-electrons of Fe in ferrocene: the excess orbital energy spectrum (EOES). RSC Adv. 5, 11933–11943. doi: 10.1039/C4RA14506B
Kalous, O., Conklin, D., Desai, A. J., O'Brien, N. A., Ginther, C., Anderson, L., et al. (2012). Dacomitinib (PF-00299804), an irreversible Pan-HER inhibitor, inhibits proliferation of HER2-amplified breast cancer cell lines resistant to trastuzumab and lapatinib. Mol. Cancer Ther. 11, 1978–1987. doi: 10.1158/1535-7163.MCT-11-0730
Khattab, M., Chatterjee, C., Clayton, A. H. A, and Wang, F. (2016). Tyrosine kinase inhibitor (AG-1478): two conformers of a tyrosine kinase inhibitor (AG-1478) disclosed using simulated UV-vis absorption spectroscopy. New J. Chem. 40, 8296–8304. doi: 10.1039/C6NJ01909A
Kobayashi, Y., Fujino, T., Nishino, M., Koga, T., Chiba, M., Sesumi, Y., et al. (2018). EGFR T790M and C797S mutations as mechanisms of acquired resistance to dacomitinib. J. Thorac. Oncol. 13, 727–731. doi: 10.1016/j.jtho.2018.01.009
Lau, S. C. M., Batra, U., Mok, T. S. K., and Loong, H. H. (2019). Dacomitinib in the management of advanced non-small-cell lung cancer. Drugs 79, 823–831. doi: 10.1007/s40265-019-01115-y
Momeny, M., Zarrinrad, G., Moghaddaskho, F., Poursheikhani, A., Sankanian, G., Zaghal, A., et al. (2017). Dacomitinib, a pan-inhibitor of ErbB receptors, suppresses growth and invasive capacity of chemoresistant ovarian carcinoma cells. Sci. Rep. 7:4204. doi: 10.1038/s41598-017-04147-0
Passaro, A., and de Marinis, F. (2018). Dacomitinib in EGFR-positive non-small cell lung cancer: an attractive but broken option. Transl. Lung Cancer Res. 7, S100–S102. doi: 10.21037/tlcr.2018.02.09
Pedersen, M. W., Pedersen, N., Ottesen, L. H., and Poulsen, H. S. (2005). Differential response to gefitinib of cells expressing normal EGFR and the mutant EGFRvIII. Br J Cancer. 93, 915–923. doi: 10.1038/sj.bjc.6602793
Qiu, X., Lin, Q., Ning, Z., Qian, X., Li, P., Ye, L., et al. (2019). Quantitative bioanalytical assay for the human epidermal growth factor receptor (HER) inhibitor dacomitinib in rat plasma by UPLC-MS/MS. J. Pharm. Biomed. Anal. 166, 66–70. doi: 10.1016/j.jpba.2018.12.041
Reed, J. E., and Smaill, J. B. (2016). “The discovery of dacomitinib, a potent irreversible EGFR inhibitor,” in Comprehensive Accounts of Pharmaceutical Research and Development: From Discovery to Late-Stage Process Development, eds A. F. Abdel-Magid, J. Pesti, and R. Vaidyanathan (Washington, DC: ACS Publication), 207–233. doi: 10.1021/bk-2016-1239.ch008
Roeper, J., and Griesinger, F. (2019). Epidermal growth factor receptor tyrosine kinase inhibitors in advanced nonsmall cell lung cancer: what is the preferred first-line therapy? Curr. Opin. Oncol. 31, 1–7. doi: 10.1097/CCO.0000000000000495
Roskoski, R. (2019). Properties of FDA-approved small molecule protein kinase inhibitors. Pharmacol Res. 144, 19–50. doi: 10.1016/j.phrs.2019.03.006
Shagufta, and Ahmad, I. (2017). An insight into the therapeutic potential of quinazoline derivatives as anticancer agents. Med. Chem. Comm. 8, 871–885. doi: 10.1039/C7MD00097A
Shirley, M. (2018). Dacomitinib: first global approval. Drugs 78, 1947–1953. doi: 10.1007/s40265-018-1028-x
Solyanik, G. I. (2019). Quinazoline compounds for antitumor treatment. Exp. Oncol. 41, 3–6. doi: 10.32471/exp-oncology.2312-8852.vol-41-no-1.12414
Steiner, T. (2002). The hydrogen bond in the solid state. Angew. Chem. Int. Ed. 41, 48–76. doi: 10.1002/1521-3773(20020104)41:1<48::AID-ANIE48>3.0.CO;2-U
Van Dongen, M., Khattab, M., Clayton, A. H. A., and Wang, F. (2017). Exploring the optical reporting characteristics of drugs—UV-Vis spectra and conformation of tyrosine kinase inhibitor SKF86002. N. J. Chem. 41, 14567–14573. doi: 10.1039/C7NJ03361C
Wang, F., and Chatterjee, S. (2017). Dominant carbons in trans- and cis-resveratrol isomerization. J. Phys. Chem. B 121, 4745–4755. doi: 10.1021/acs.jpcb.7b02115
Wu, Y. L., Cheng, Y., Zhou, X., Lee, K. H., Nakagawa, K., Niho, S., et al. (2017). Dacomitinib versus gefitinib as first-line treatment for patients with EGFR-mutation-positive non-small-cell lung cancer (ARCHER 1050): a randomised, open-label, phase 3 trial. Lancet Oncol. 18, 1454–1466. doi: 10.1016/S1470-2045(17)30608-3
Yuan, J., Yuan, Y., Tian, X., Liu, Y., and Sun, J. (2017). Insights into the photobehavior of fluorescent oxazinone, quinazoline, and difluoroboron derivatives: molecular design based on the structure–property relationships. J. Phys. Chem. C 121, 8091–8108. doi: 10.1021/acs.jpcc.7b01360
Keywords: tyrosine kinase inhibitor (EGFR-TKI), drug structures, density functional theory (DFT) calculations, spectroscopy, off-target
Citation: Kabir ML, Backler F, Clayton AHA and Wang F (2020) Deducing the Conformational Properties of a Tyrosine Kinase Inhibitor in Solution by Optical Spectroscopy and Computational Chemistry. Front. Chem. 8:596. doi: 10.3389/fchem.2020.00596
Received: 27 February 2020; Accepted: 08 June 2020;
Published: 28 July 2020.
Edited by:
Jamie Platts, Cardiff University, United KingdomReviewed by:
Mathieu Linares, Linköping University, SwedenArtur Nenov, University of Bologna, Italy
Copyright © 2020 Kabir, Backler, Clayton and Wang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Andrew H. A. Clayton, aclayton@swin.edu.au; Feng Wang, fwang@swin.edu.au