Deducing the Conformational Properties of a Tyrosine Kinase Inhibitor in Solution by Optical Spectroscopy and Computational Chemistry

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 anilinequinazoline 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).
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 offtarget 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 SCHEME 1 | Chemical structures of Dacomitinib and Gefitinib and their relationship.
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 wildtype 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 IC 50 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 ultravioletvisible (UV-Vis) absorption of Dacomitinib in dimethylsulfoxide (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.

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 (C 24 H 25 ClFN 5 O 2 ) 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).

Dacomitinib Conformers in Crystal Phase and in Isolation
The chemical structure and nomenclature of Dacomitinib (C 24 H 25 ClFN 5 O 2 , 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-ylbuten-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-4fluoroaniline ring (C) with respect to the quinazoline ring (Rings A and B). The O-bridge provides small tunes to the methoxy group (-OCH 3 ) 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-1yl-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.
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 . 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.
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 . 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) 17), the former controls the relevant orientation of the methoxy group to the quinazoline ring and the latter controls the 3-chloro-4fluoroaniline 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 . 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. 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.

Spectral Signatures of Dacomitinib Conformers in the UV-Vis Absorption Spectra
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 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. 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 C 1 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). 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) Hlinker toward the piperidine side chain in the Dacomitinib conformers. Folding the 3-chloro-4-fluoroaniline ring (Ring C) 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 . 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) H 3 ) 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 lonepair electron containing O (4) and the electron less populated in N (6) , which are represented as intense red protrusion and deep blue notch, respectively.
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 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.
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 = ε lm i -ε i gm , in which ε lm i is the orbital energies of a conformer and ε i gm 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).  Frontiers in Chemistry | www.frontiersin.org FIGURE 7 | Excess orbital energy spectrum (EOES) of Dacomitinib in DMSO at the B3LYP/6-311G** level of theory: ε i = ε lm i -ε i gm , where ε lm i = ith orbital energy of a local minimum conformer (i.e., Daco-P2, Daco-T1, or Daco-T2) and ε i gm = 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 (C 24 H 25 ClFN 5 O 2 ) has a total of 123 doubly occupied orbitals, which consists of 33 1s core orbitals. While most of the core 1s EOES ( ε i = ε lm i -ε i gm ) 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  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 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.
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.