Next Article in Journal
Research on Molecular Structure and Electronic Properties of Ln3+ (Ce3+, Tb3+, Pr3+)/Li+ and Eu2+ Co-Doped Sr2Si5N8 via DFT Calculation
Next Article in Special Issue
Pure and Hybrid SCAN, rSCAN, and r2SCAN: Which One Is Preferred in KS- and HF-DFT Calculations, and How Does D4 Dispersion Correction Affect This Ranking?
Previous Article in Journal
Modifications of Plasma Membrane Organization in Cancer Cells for Targeted Therapy
Previous Article in Special Issue
Experimental and Density Functional Theory Studies on 1,1,1,4,4,4-Hexafluoro-2-Butene Pyrolysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Computational Study on the Mechanism of the Photouncaging Reaction of Vemurafenib: Toward an Enhanced Photoprotection Approach for Photosensitive Drugs

by
Abdulilah Dawoud Bani-Yaseen
Department of Chemistry & Earth Sciences, Faculty of Arts & Science, Qatar University, Doha 2713, Qatar
Molecules 2021, 26(7), 1846; https://doi.org/10.3390/molecules26071846
Submission received: 15 August 2020 / Revised: 21 September 2020 / Accepted: 1 October 2020 / Published: 25 March 2021

Abstract

:
The photochemical behavior of the photosensitive first-line anticancer drug vemurafenib (VFB) is of great interest due to the impact of such behavior on its pharmacological activity. In this work, we computationally elucidated the mechanism of the photoinduced release of VFB from the 4,5-dimethoxy-2-nitrobenzene (DMNB) photoprotecting group by employing various density functional theory (DFT)/time-dependent DFT (TD-DFT) approaches. The computational investigations included a comparative assessment of the influence of the position of the photoprotecting group as a substituent on the thermodynamics and kinetics of the photouncaging reactions of two VFB-DMNB prodrugs, namely pyrrole (NP) and sulfonamide (NS). With the aid of the DFT calculations concerning the activation energy barrier (∆G), the obtained results suggest that the step of the photoinduced intramolecular proton transfer of the DMNB moiety is not detrimental concerning the overall reaction profile of the photouncaging reaction of both prodrugs. However, the obtained results suggested that the position of the substitution position of the DMNB photoprotecting group within the prodrug structure has a substantial impact on the photouncaging reaction. In particular, the DMNB-Ns-VFB prodrug exhibited a notable increase in ∆G for the key step of ring opining within the DMNB moiety indicative of potentially hindered kinetics of the photouncaging process compared with DMNB-Np-VFB. Such an increase in ∆G may be attributed to the electronic influence of the NP fragment of the prodrug. The results reported herein elaborate on the mechanism of the photoinduced release of an important anticancer drug from photoprotecting groups with the aim of enhancing our understanding of the photochemical behavior of such photosensitive pharmaceutical materials at the molecular level.

1. Introduction

Examining the physicochemical properties of pharmaceutical materials at the molecular level can enhance our understanding of various associated biological processes that are important in drug discovery. Examining such molecular properties of pharmaceutical materials has recently attracted notable attention [1,2,3,4,5,6,7,8,9,10,11,12,13]. The phototoxicity of numerous pharmaceutical materials is a major therapeutic concern that can conspicuously constrain their clinical applications [14,15,16,17,18]. Several approaches have been recently developed toward reducing the impact of such undesired side effects, including photocaging [19,20,21,22,23,24]. Recently, the importance of the incorporation of photocaging groups into photoactive pharmaceutical materials has been increasingly recognized as an effectual approach toward reducing the phototoxicity of a wide spectrum of photoactive drugs [20,21,22,23,24].
In such an approach, a photoremovable photocaging moiety is selectively attached to specific positions of a molecule of interest via covalent bonding to afford a photocaged product of rendered photochemical activities. Subsequently, the functionality of the parent molecules can be released upon the elimination of the photocaging group via a photochemical process [25,26,27]. Such a process of uncaging is substantially important in the landscape of pharmaceutical and medicinal chemistry, where it requires irradiation with light of a specific wavelength (λirr) that is appropriate for uncaging the drug under study.
As such, several studies have been recently reported concerning the photocaging of materials of biological and pharmaceutical importance, such as the antimelanoma agent vemurafenib [20], anticancer drug doxazolidine [21], anticancer kinase inhibitor imatinib [22], small-molecule tubulin inhibitor [23], and photocaged materials for homeostasis [24]. Various types of photolabile-protecting groups (PPGs) have been developed and utilized in these works to afford the photocaged prodrugs of interest. Such photocaging groups have been mainly utilized to selectively photoprotect specific functional groups that are feasible for photoreactions; these include amines, hydroxyl, and carboxyl functional groups.
Interest in the applications of the antimelanoma agent vemurafenib (VFB) has recently grown. The chemical structure of VFB is shown in Figure 1. VFB has attracted such notable attention due to its therapeutic efficacy as an anticancer drug. However, although VFB is considered as a first-line anticancer drug, several studies have reported the severe phototoxicity associated with the administration of the drug [28,29,30,31,32,33]. On the other hand, investigating the physicochemical properties of materials of biological and medicinal applications at the molecular level is essential for a better understanding of their activities and functionalities [21,22,23,24,25,26,27,28,29,30,31]. As such, investigating the photochemical properties of photosensitive biological materials, such as the first-line anticancer drug VFB, at the molecular level is necessary to enhance the efficacy of the photocaging process and, consequently, reduce the potential phototoxicity.
In examining the structure, VFB bears two secondary amine functional groups that are feasible for photochemical instability. These two amines are within two moieties, namely the azaindole and sulfonamide moieties. Employing the photocaging approach for these two photoactive centers toward reducing the phototoxicity of VFB was recently reported by Peifer and co-workers using derivatives of 2-nitrotoluene, such as 4,5-dimethoxy-2-nitrobenzyl (DMNB) [20]. The reported results suggested a promising approach for photocaging and photouncaging of VFB using the DMNB moiety. However, a requested feature that is expected from an efficient photocaging process of biological materials is a reversible uncaging of the photoactive drug under ambient conditions and appropriate microenvironments. Importantly, such requirements necessitate investigating various mechanistic aspects of the uncaging process under biologically mimicked microenvironments. In the literature, there are several studies on the mechanism of photouncaging of 2-nitrotoluene derivatives [20,25,34,35,36]. In view of the suggested mechanisms in these references, a plausible mechanism for the photocleavage reaction of VFB from DMNB can be suggested, as illustrated in Scheme 1.
In light of the above-mentioned considerations, we provide, in this work, computational insights concerning selected mechanistic aspects of the photouncaging reaction of the DMNB-photoprotected version of the first-line anticancer drug VFB. By employing systematic density functional theory (DFT) calculations in implicit aqueous solutions, we address important scientific concerns at the molecular level related to the mechanism of the photouncaging reaction of VFB photoprotected by DMNB, which is induced by irradiation. The results presented in this work can help enhance our understanding of this important chemical process toward the development of efficient chemical approaches that can lead to reducing the phototoxicity of photosensitive drugs.

2. Results and Discussion

As can be noted from Figure 1, the structure of VFB bears two amine groups that are feasible for photoprotection with DMNB. These two amine groups are ascribed as pyrrole (NP) and sulfonamide (NS). These are ascribed with the correspondence of the moieties that bear these two amine groups, namely the sulfoxide and pyrrole moieties, respectively. Figure 2 displays the chemical structures of the DMNB-photocaged VFB prodrugs, namely 1 and 2. The ascriptions of 1 and 2 correspond to the photoprotected VFB at the NP and NS positions, respectively.
According to the mechanism illustrated in Scheme 1, the photouncaging reaction of the DMNB-protected substance is triggered by irradiation that induces a proton transfer in the excited state to afford I1. Such a step is plausibly considered as a fast step that has a minimal influence on the overall kinetics of the photouncaging reaction [34,35,36]. As such, this particular step was not considered in the DFT calculations concerning the transformation of 1 and 2 into their corresponding I1. Hence, we focused on the other three steps in this study. The other three steps of the photouncaging reaction can be briefly described. As such, following the photoinduced intramolecular proton transfer in 1 and 2, I1 undergoes an intramolecular cyclization reaction between C′′1 and O′6 to afford I2. Then, I2 undergoes a concurrent a ring-opening and proton transfer to yield I3 followed by the breakage of the N-C′1 bond to release the parent drug VFB. We also examined the transition state (TS) for each of these three steps of the photouncaging of 1 and 2 as well.
Geometry optimization. Figure 3 shows the DFT-optimized geometry of 1 and 2 and the corresponding intermediates that are generated via the photouncaging reaction. The accuracy of the optimized geometry was verified upon comparison with the previously reported X-ray structures obtained from the Cambridge Crystallographic Data Centre (CCDC). The key structural parameters of the examined species, namely the bond lengths of NP-C′1 and NS-C′1 are displayed in Figure 3. Comparing 1 vs. 2, the calculations revealed a shorter bond length for 1 with a ∆l of 0.018 Å. Additionally, comparing the corresponding intermediates of 1 vs. their counterparts of 2, this variation in the l(N-C′)key was retained for I1, whereas no significant differences were noted for I2 and I3. On the other hand, upon examining the optimized geometry of 1 and 1-I1, there was substantial change in the dihedral angle that comprises the NP-C′1 and the DMNB moiety. Hence, the formation of 1-I1 requires a ~100° rotation in this dihedral angle of the DMNB moiety around the NP-C′1 bond axis. This kind of rotation requirement may increase the energy barrier for the 11-I1 transformation, which in turn might be crucial for affording the photouncaging process of 1. The relative rotation requirement around the key-bond of the photouncaging reaction for I2 and I3 of both prodrugs is not significantly appreciable.
According to the mechanism illustrated in Scheme 1, we performed a geometry optimization of TSs for the three critical steps following the intramolecular proton transfer induced by irradiation employing the DFT method with the 6-31G basis set (DFT/6-31G) in a vacuum. The optimized geometry of TS1-3 of the photouncaging of 1 and 2 are displayed in Figure 4.
Figure 4 illustrates the distance (dA…B) in Å between the two key atoms involved in the key TS of the photouncaging reactions of 1 and 2. As stated above, TS1 refers to the transformation I1TS1I2, which corresponds to the formation of the C′1-O′6 bond ascribed as an intramolecular cyclization reaction. Comparing the 1-TS1 vs. 2-TS2, the DFT calculation revealed a dC′1…O′6 of 2.245 and 2.300 Å, respectively, revealing a difference of 0.055 Å. Such a small increase in the distance observed for 1-TS1 compared to 2-TS1 is indicative of the relatively comparable energy barrier that needs to be overcome for both intermediates. For the I2TS2I3 transformation, we found a dN′4…O′6 of 2.134 and 1.741 Å for 1-TS2 and 2-TS2, respectively, indicative of a relatively more feasible ring-opining for 1-TS2 compared to 2-TS2. As such, one may anticipate that a higher energy barrier is associated with the formation of 2-TS2 compared to 1-TS2.
Studying the photochemical behavior of a substance necessitates investigating the potential electronic transition associated with being irradiated with light of specific intensity and λ. On the other hand, for safe and efficient pharmacological applications, reducing the damage to the tissues that is caused by the irradiation of the uncaging process can be achieved by increasing the effective λ of the irradiation process. One important aspect of the photouncaging mechanism of the 1 and 2 prodrugs that we consider herein is to assess their photochemical behavior upon being irradiated with the same light source of specific λ.
Thus, investigating their spectral behavior is crucial toward providing better comparative rationalization of such potential discrepancies in their photochemical behaviors. As such, we used the time-dependent DFT (TD-DFT) approach to compute the UV–Vis absorption spectra of the parent compounds, namely VFB and DMNB, as well as prodrugs 1 and 2 in implicit aqueous solutions (Figure 5). First, the accuracy of the calculated spectra of VBF, 1, and 2 was verified by comparing the simulated spectra vs. the experimental counterparts [20,37]. As can be noted from Figure 5, VFB exhibits two major absorption bands positioned at λ of 253 and 315 nm. For 1 and 2, a third major absorption band appeared in the visible range at a λ of ~ 390 nm. However, comparing the absorption spectrum of VFB vs. DMNB, 1 and 2, one can infer that the absorption band of λ390 corresponds to the DMNB moiety.
The TD-DFT results revealed that all three major transition bands illustrated in Figure 5 correspond to π→π* electronic transitions. The key molecular orbitals (MOs) involved in such transitions are displayed in Figure 6. For VFB, both the highest molecular occupied orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) are delocalized over the π-conjugated part of the molecule comprising the phenyl, pyridine, and pyrrole moieties. The calculations revealed an energy gap (∆EHOMO-LUMO) of 4.46 and 3.68 eV for VFB and DMNB, respectively. One can anticipate that the lone-pair (lp) of the NP can participate in the charge delocalization that accounts for the shape and nature of these frontier MOs.
Comparing the MOs of VFB vs. 1 and 2, we infer that the HOMO and LUMO of VFB become a HOMO-1 and LUMO+1 in 1 and 2, respectively, whereas the HOMO and LUMO of 1 and 2 correspond to their counterparts of DMNB. However, the incorporation of DMNB into VFB via a photocaging reaction induced an increase in the ∆E of the MOs of 0.065 eV. This increase in the ∆E corresponds to the stabilization and destabilization of 0.025 and 0.040 in the energy of HOMO and LUMO, respectively, indicative of a blue-shift in the electronic transition and consequently reduced photoreactivity in the visible range of irradiation. Comparing the HOMO of 1 vs. the counterpart of 2, the HOMO comprises both HOMOs of VFB and DMNB, whereas for 2, it comprises only the HOMO of DMNB. In light of these changes in the nature of the MOs after photocaging, we suggest that 1 has to overcome higher energy barriers concerning the photouncaging reaction, which in turn is well in line with the results obtained from the DFT calculations obtained from the optimized geometry.
We further extended our investigations concerning the role of the position of photoprotection on the photouncaging reaction by conducting a natural bond orbital (NBO) analysis. In this approach, the contributions of NP and NS in the mechanistic pathway of 1 and 2 were evaluated. We focused herein on the energy of interaction between the MOs of the key atoms involved in the photouncaging reaction. Such types of interactions are presented in terms of an electronic transition from a donating MO to an accepting MO; i.e., DMO→AMO. The magnitude of interaction was calculated with respect to the second-order perturbation energy (E(2)) as revealed by the DFT-based NBO analysis.
The obtained results from the NBO analyses of all species involved in the photouncaging reaction of 1 and 2 are compiled in Table 1. The superpositions of the selected NBO and the corresponding energy of interactions are displayed in Figure 7. It is worth noting that NP and NS in VFB exhibited notable NBO discrepancies. NP exhibited a substantial two electron transition of the types lp(NP)→π*(C1-C2) and lp(NP)→π*(C3-C4) with E(2) of 42.0 and 35.9 kcal/mol, respectively; whereas NS exhibited two transitions of the types lp(NS)→σ*(S-O) and lp(NS)→σ*(C5-C6) with E(2) of 8.2 and 7.8 kcal/mol, respectively. These results revealed strong interactions between the lp(NP) and adjacent MOs indicative of enhanced charge delocalization, which in turn indicates a reduced contribution in the photouncaging reaction.
As such, we suggest that the lp(NS) can afford more contribution toward promoting the photouncaging reaction for 2 compared with 1. The formation of 1-I1 can trigger a substantial deformation in the delocalization of the NP moiety. In particular, upon the formation of 1-I1, a loss in the lone-pair character of NP is observed. Under these circumstances, the NP atom exhibits other types of stabilizations, namely π(NP-C1)→π*(C′1-C′2) with E(2) of 13.2 kcal/mol indicative of charge delocalization that is well in line with the frontier HOMO of prodrug 1. These results suggest that the change in the types of NBO electronic transition may cause a diminution in the relative reactivity of 1 compared with 2 toward the formation of I1.
The lp(NP) retained its characteristic upon the formation of I2 and I3. On the other hand, the main transitions of lp(NS) retained their characters, including the two main transitions. Concerning the transition between the lp(NP) and lp(NS) to the antibonding NBO of the key atoms of the DMNB moiety, namely σ* and π*, we suggest that such types of transition can promote bond formation breakage toward affording the species of interest. For example, 1-I2 and 2-I2 exhibited the transitions of lp(NP)→σ*(C′1-O′6) and lp(NS)→σ*(C′1-O′6), respectively, with stabilizing energies of 16.1 and 17.6 kcal/mol, respectively. These results suggest that such transitions can stimulate the breakage of the C′1-O′6 bond and consequently the ring-opening and proton transfer toward affording I3. Indeed, alike transitions can be noted in the TS2 of both intermediates indicative of the important contribution of lp(NP) and lp(NS) in the photouncaging reaction of 1 and 2.
We further extended our NBO analyses via examining the variation in the NBO charges of the key atoms involved in the photouncaging reaction of prodrugs 1 and 2; see Table 2. As anticipated, NP and NS exhibited negative charges for all species that can be attributed to the higher electronegativity of the nitrogen atom. However, comparing the NP vs. the NS for VFB, the calculations revealed NBO charges of −0.534 and −0.869 e, respectively. This relatively lower value for NP compared to NS can be attributed to the incorporation of the lp(NP) in the conjugation system as demonstrated by the frontier MO and NBO analyses.
However, both exhibited relatively similar behaviors concerning the decrease in their values upon the formation of the prodrugs 1 and 2, where a decrease in the NBO charge by 0.175 and 0.202 e was calculated for 1 and 2, respectively. On the other hand, NP and NS exhibited relatively negligible changes in their NBO charges in the corresponding intermediates and TSs. Importantly, the key atom C′1 exhibited a substantial variation in its NBO charge as the photouncaging reaction progressed. For example, C′1 exhibited approximately the same values of NBO charges for 1 and 2, whereas values of 0.019 and −0.005 were calculated for 1-I1 and 2-I1, respectively.
Similarly, for 1-TS1 and 2-TS1, values of 0.123 and 0.057 were calculated for the NBO charges, respectively, indicative of an ∆NBO charge of 0.066 e. Nevertheless, other species of 1 compared to 2 exhibited an insignificant type of ∆NBO. Indeed, these results indicate the significant influence of the first step, namely irradiation-induced prodrug→I1, in distinguishing the photochemical behavior of the prodrugs 1 and 2. Such observations were further illustrated as shown in the overall energy profile of the photouncaging reactions of 1 and 2.
Reaction Profile. As illustrated in Figure 8, we constructed the overall energy profile for the photouncaging reaction of 1 and 2 with the aid of DFT/B3LYP /6-31G calculations and took into consideration the corresponding key steps. The overall energy was constructed in terms of the relative change Gibbs free energy (∆G°) utilizing the results of the geometry optimization of intermediates and TSs. The displayed results were calculated relative to the energy of the reactants, namely prodrugs 1 and 2. For example, the value of 33.7 kcal/mol indicates the difference in Gibbs free energy between 1-I1 and 1; i.e., G°(1-I1) − G°(1).
As can be inferred from Figure 8, the DFT calculations revealed a ∆G° of 33.7 and 5.2 kcal/mol for the formation of 1-I1 and 2-I1, which is indicative of thermodynamically more favorable formation for 2-I1 compared with 1-I1 upon irradiation. This notable difference is indicative of a relatively more favorable photo-conversion of 2 compared to 1, which is in good agreement with the experimental results reported for these two prodrugs [20]. One can infer from Figure 8 that the I1TS1I2 transformations of both prodrugs are associated with comparable energy barriers of activation (∆G), where the ∆G values of 12.1 and 15.4 kcal/mol were calculated for the transformations of 1 and 2, respectively, which is in line with the results obtained from the geometry optimization.
On the other hand, one significant difference can be noted for the I2TS2I3. ∆G values of 1.9 and 40.3 kcal/mol were calculated for the transformations of 1 and 2, respectively. We suggest that such a significant increase in the ∆G for 2 compared to 1 might be attributed to both steric and electronic effects, as revealed by the geometry optimization calculations and NBO analyses. For the termination transformation, I3TS3Pr, ∆G values of 42.0 and 41.3 kcal/mol can be inferred from Figure 8 for 1 and 2, respectively, indicative of similar kinetic favoritism.
These results suggest that the intramolecular proton transfer step induced by irradiation is the key thermodynamic process that is responsible for discrepancies in the photochemical behavior observed for prodrugs 1 and 2. On the other hand, the second step is a key source of kinetic discrepancies between the two prodrugs, whereas the third step for both prodrugs exhibited a comparable increase in ∆G compared to the other steps. These results also suggest that both thermodynamic and kinetic factors can notably affect the photouncaging reactions of prodrugs 1 and 2.

3. Computational Methods

All calculations were performed with Gaussian 09 version D.01 [38]. Calculation of the Geometry optimization was conducted for all species comprised in the mechanism shown in Scheme 1 employing the DFT/B3LYP/6-31+G(d), unless otherwise noted. The optimized geometry was verified as minima utilizing frequency calculations. The optimized geometries for selected species were utilized as inputs to perform further computational tasks. The UV–Vis absorption spectra (company, city, abbreviated state of the parent compounds and prodrugs were simulated utilizing time-dependent DFT (TD-DFT) in implicit aqueous solutions. The implicit aqueous solution was incorporated into the DFT calculations employing the integral equation formalism polarizable continuum model (IEFPCM) [39]. Natural bond orbital (NBO) analyses were performed as implemented in Gaussian 09. The DFT method employed herein has been previously demonstrated to be efficient in reproducing the experimental results of the molecular properties of various substances in implicit solutions [8,40].

4. Conclusions

The photocaging of photosensitive pharmaceutical materials is a promising approach in pharmaceutical and biological applications toward the reduction in the potential associated phototoxicity of a drug molecule of interest. However, it is crucial to elucidate the mechanism of the photoinduced release of the parent molecule. The present study provides computational insights on the mechanism of the photouncaging reaction of the important first-line anticancer drug VFB from the photoprotecting group DMNB. Within a four-step mechanistic pathway, the obtained results concerning the thermodynamic factor and activation energy barrier revealed that the first step of the mechanism was not detrimental—namely, the photoinduced intramolecular proton transfer within the DMNB moiety.
We demonstrated that the position of the moiety DMNB as a substituent can have a substantial impact on the kinetics of the photouncaging reaction of VFB-DMNB prodrugs. Hence, these results suggest that the step of the concurrent ring-opening and intramolecular proton transfer might be the rate-determining step for prodrug 2 with an equivalent role for the termination step, whereas the termination step is the rate-determining step for prodrug 1. We anticipate that the results reported herein may provide insights on the development of photocaging/photouncaging approaches for reducing the phototoxicity of photosensitive pharmaceutical materials.

Funding

This work was supported by Qatar University [QUCD-CAS-2020-1].

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

Major aspects of the calculations were performed using the supercomputing facility at Texas A&M University at Qatar.

Conflicts of Interest

The author declares no conflict of interest.

Sample Availability

Samples of the compounds are not available from the author.

References

  1. Caron, G.; Ermondi, G. Updating molecular properties during early drug discovery. Drug Discov. Today 2017, 22, 835–840. [Google Scholar] [CrossRef] [PubMed]
  2. Caron, G.; Kihlberg, J.; Ermondi, G. Intramolecular hydrogen bonding: An opportunity for improved design in medicinal chemistry. Med. Res. Rev. 2019, 39, 1707–1729. [Google Scholar] [CrossRef] [PubMed]
  3. Dawoud Bani-Yaseen, A. Spectrofluorimetric study on the interaction between antimicrobial drug sulfamethazine and bovine serum albumin. J. Lumin. 2011, 131, 1042–1047. [Google Scholar] [CrossRef]
  4. Varma, M.V.; Lai, Y.; El-Kattan, A.F. Molecular properties associated with transporter-mediated drug disposition. Adv. Drug Deliv. Rev. 2017, 116, 92–99. [Google Scholar] [CrossRef] [PubMed]
  5. Bani-Yaseen, A.D. Synchronous spectrofluorimetric study of the supramolecular host-guest interaction of β-cyclodextrin with propranolol: A comparative study. Spectrochim. Acta. A Mol. Biomol. Spectrosc. 2015, 148, 93–98. [Google Scholar] [CrossRef]
  6. Tarcsay, Á.; Keserű, G.M. Contributions of Molecular Properties to Drug Promiscuity. J. Med. Chem. 2013, 56, 1789–1795. [Google Scholar] [CrossRef]
  7. Bani-Yaseen, A.D. Computational molecular perspectives on the interaction of propranolol with β-cyclodextrin in solution: Towards the drug-receptor mechanism of interaction. J. Mol. Liq. 2017, 227, 280–290. [Google Scholar] [CrossRef]
  8. Bani-Yaseen, A.D.; Al-Balawi, M. The solvatochromic, spectral, and geometrical properties of nifenazone: A DFT/TD-DFT and experimental study. Phys. Chem. Chem. Phys. 2014, 16, 15519–15526. [Google Scholar] [CrossRef]
  9. González, M.; Ellahioui, Y.; Álvarez, R.; Gallego-Yerga, L.; Caballero, E.; Vicente-Blázquez, A.; Ramudo, L.; Marín Folgado, M.; Sanz, C.; Medarde, M.; et al. The Masked Polar Group Incorporation (MPGI) Strategy in Drug Design: Effects of Nitrogen Substitutions on Combretastatin and Isocombretastatin Tubulin Inhibitors. Molecules 2019, 24, 4319. [Google Scholar] [CrossRef] [Green Version]
  10. Bani-Yaseen, A.D.; Hammad, F.; Ghanem, B.S.; Mohammad, E.G. On the Photophysicochemical Properties of Selected Fluoroquinolones: Solvatochromic and Fluorescence Spectroscopy Study. J. Fluoresc. 2013, 23, 93–101. [Google Scholar] [CrossRef]
  11. Das, S.S.; Alkahtani, S.; Bharadwaj, P.; Ansari, M.T.; ALKahtani, M.D.F.; Pang, Z.; Hasnain, M.S.; Nayak, A.K.; Aminabhavi, T.M. Molecular insights and novel approaches for targeting tumor metastasis. Int. J. Pharm. 2020, 585, 119556. [Google Scholar] [CrossRef] [PubMed]
  12. Steyn, S.J.; Varma, M.V.S. Cytochrome-P450-Mediated Drug – Drug Interactions of Substrate Drugs: Assessing Clinical Risk Based on Molecular Properties and an Extended Clearance Classification System. Mol. Pharm. 2020, 17, 3024–3032. [Google Scholar] [CrossRef] [PubMed]
  13. Bani-Yaseen, A.D. Solvatochromic and Fluorescence Behavior of Sulfisoxazole. J. Fluoresc. 2011, 21, 1061–1067. [Google Scholar] [CrossRef] [PubMed]
  14. Su, K.A.; Habel, L.A.; Achacoso, N.S.; Friedman, G.D.; Asgari, M.M. Photosensitizing antihypertensive drug use and risk of cutaneous squamous cell carcinoma. Br. J. Dermatol. 2018, 179, 1088–1094. [Google Scholar] [CrossRef]
  15. Monteiro, A.F.; Rato, M.; Martins, C. Drug-induced photosensitivity: Photoallergic and phototoxic reactions. Clin. Dermatol. 2016, 34, 571–581. [Google Scholar] [CrossRef]
  16. Peukert, S.; Nunez, J.; He, F.; Dai, M.; Yusuff, N.; DiPesa, A.; Miller-Moslin, K.; Karki, R.; Lagu, B.; Harwell, C.; et al. A method for estimating the risk of drug-induced phototoxicity and its application to smoothened inhibitors. Medchemcomm 2011, 2, 973. [Google Scholar] [CrossRef]
  17. Kim, W.B.; Shelley, A.J.; Novice, K.; Joo, J.; Lim, H.W.; Glassman, S.J. Drug-induced phototoxicity: A systematic review. J. Am. Acad. Dermatol. 2018, 79, 1069–1075. [Google Scholar] [CrossRef]
  18. Khandpur, S.; Porter, R.M.; Boulton, S.J.; Anstey, A. Drug-induced photosensitivity: New insights into pathomechanisms and clinical variation through basic and applied science. Br. J. Dermatol. 2017, 176, 902–909. [Google Scholar] [CrossRef]
  19. Wu, G.; Zhao, T.; Kang, D.; Zhang, J.; Song, Y.; Namasivayam, V.; Kongsted, J.; Pannecouque, C.; De Clercq, E.; Poongavanam, V.; et al. Overview of Recent Strategic Advances in Medicinal Chemistry. J. Med. Chem. 2019, 62, 9375–9414. [Google Scholar] [CrossRef]
  20. Horbert, R.; Pinchuk, B.; Davies, P.; Alessi, D.; Peifer, C. Photoactivatable Prodrugs of Antimelanoma Agent Vemurafenib. ACS Chem. Biol. 2015, 10, 2099–2107. [Google Scholar] [CrossRef]
  21. Tamura, R.; Balabanova, A.; Frakes, S.A.; Bargmann, A.; Grimm, J.; Koch, T.H.; Yin, H. Photoactivatable Prodrug of Doxazolidine Targeting Exosomes. J. Med. Chem. 2019, 62, 1959–1970. [Google Scholar] [CrossRef] [PubMed]
  22. Zindler, M.; Pinchuk, B.; Renn, C.; Horbert, R.; Döbber, A.; Peifer, C. Design, Synthesis, and Characterization of a Photoactivatable Caged Prodrug of Imatinib. ChemMedChem 2015, 10, 1335–1338. [Google Scholar] [CrossRef] [PubMed]
  23. Döbber, A.; Phoa, A.F.; Abbassi, R.H.; Stringer, B.W.; Day, B.W.; Johns, T.G.; Abadleh, M.; Peifer, C.; Munoz, L. Development and Biological Evaluation of a Photoactivatable Small Molecule Microtubule-Targeting Agent. ACS Med. Chem. Lett. 2017, 8, 395–400. [Google Scholar] [CrossRef] [Green Version]
  24. Basa, P.N.; Antala, S.; Dempski, R.E.; Burdette, S.C. A Zinc(II) Photocage Based on a Decarboxylation Metal Ion Release Mechanism for Investigating Homeostasis and Biological Signaling. Angew. Chem. 2015, 127, 13219–13223. [Google Scholar] [CrossRef] [Green Version]
  25. Klán, P.; Šolomek, T.; Bochet, C.G.; Blanc, A.; Givens, R.; Rubina, M.; Popik, V.; Kostikov, A.; Wirz, J. Photoremovable Protecting Groups in Chemistry and Biology: Reaction Mechanisms and Efficacy. Chem. Rev. 2013, 113, 119–191. [Google Scholar] [CrossRef] [PubMed]
  26. Pelliccioli, A.P.; Wirz, J. Photoremovable protecting groups: Reaction mechanisms and applications. Photochem. Photobiol. Sci. 2002, 1, 441–458. [Google Scholar] [CrossRef]
  27. Šolomek, T.; Wirz, J.; Klán, P. Searching for Improved Photoreleasing Abilities of Organic Molecules. Acc. Chem. Res. 2015, 48, 3064–3072. [Google Scholar] [CrossRef]
  28. Brugière, C.; Stefan, A.; Morice, C.; Cornet, E.; Moreau, A.; Allouche, S.; Verneuil, L. Vemurafenib skin phototoxicity is indirectly linked to ultraviolet A minimal erythema dose decrease. Br. J. Dermatol. 2014, 171, 1529–1532. [Google Scholar] [CrossRef]
  29. Woods, J.A.; Ferguson, J.S.; Kalra, S.; Degabriele, A.; Gardner, J.; Logan, P.; Ferguson, J. The phototoxicity of vemurafenib: An investigation of clinical monochromator phototesting and in vitro phototoxicity testing. J. Photochem. Photobiol. B Biol. 2015, 151, 233–238. [Google Scholar] [CrossRef]
  30. Sinha, R.; Larkin, J.; Gore, M.; Fearfield, L. Cutaneous toxicities associated with vemurafenib therapy in 107 patients with BRAF V600E mutation-positive metastatic melanoma, including recognition and management of rare presentations. Br. J. Dermatol. 2015, 173, 1024–1031. [Google Scholar] [CrossRef]
  31. Boudon, S.M.; Plappert-Helbig, U.; Odermatt, A.; Bauer, D. Characterization of Vemurafenib Phototoxicity in a Mouse Model. Toxicol. Sci. 2014, 137, 259–267. [Google Scholar] [CrossRef] [PubMed]
  32. Gelot, P.; Dutartre, H.; Khammari, A.; Boisrobert, A.; Schmitt, C.; Deybach, J.-C.; Nguyen, J.-M.; Seité, S.; Dréno, B. Vemurafenib: An unusual UVA-induced photosensitivity. Exp. Dermatol. 2013, 22, 297–298. [Google Scholar] [CrossRef] [PubMed]
  33. Boussemart, L.; Boivin, C.; Claveau, J.; Tao, Y.G.; Tomasic, G.; Routier, E.; Mateus, C.; Deutsch, E.; Robert, C. Vemurafenib and Radiosensitization. JAMA Dermatol. 2013, 149, 855. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Il’ichev, Y.V.; Schwörer, M.A.; Wirz, J. Photochemical Reaction Mechanisms of 2-Nitrobenzyl Compounds: Methyl Ethers and Caged ATP. J. Am. Chem. Soc. 2004, 126, 4581–4595. [Google Scholar] [CrossRef]
  35. Rhee, H.; Lee, J.-S.; Lee, J.; Joo, C.; Han, H.; Cho, M. Photolytic Control and Infrared Probing of Amide I Mode in the Dipeptide Backbone-Caged with the 4,5-Dimethoxy-2-nitrobenzyl Group. J. Phys. Chem. B 2008, 112, 2128–2135. [Google Scholar] [CrossRef]
  36. Sebej, P.; Solomek, T.; Hroudna, L.; Brancova, P.; Klan, P. Photochemistry of 2-Nitrobenzylidene Acetals. J. Org. Chem. 2009, 74, 8647–8658. [Google Scholar] [CrossRef]
  37. Morlière, P.; Boscá, F.; Silva, A.M.S.; Teixeira, A.; Galmiche, A.; Mazière, J.-C.; Nourry, V.; Ferreira, J.; Santus, R.; Filipe, P. A molecular insight into the phototoxic reactions observed with vemurafenib, a first-line drug against metastatic melanoma. Photochem. Photobiol. Sci. 2015, 14, 2119–2127. [Google Scholar] [CrossRef]
  38. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision D.01; Gaussian, Inc.: Wallingford, CT, USA, 2009. [Google Scholar]
  39. Tomasi, J.; Mennucci, B.; Cancès, E. The IEF version of the PCM solvation method: An overview of a new method addressed to study molecular solutes at the QM ab initio level. J. Mol. Struct. THEOCHEM 1999, 464, 211–226. [Google Scholar] [CrossRef]
  40. Shkoor, M.; Mehanna, H.; Shabana, A.; Farhat, T.; Bani-Yaseen, A.D. Experimental and DFT/TD-DFT computational investigations of the solvent effect on the spectral properties of nitro substituted pyridino [3,4-c]coumarins. J. Mol. Liq. 2020, 313, 113509. [Google Scholar] [CrossRef]
Figure 1. The chemical structures of vemurafenib (VFB) and 4,5-dimethoxy-2-nitrobenzene (DMNB).
Figure 1. The chemical structures of vemurafenib (VFB) and 4,5-dimethoxy-2-nitrobenzene (DMNB).
Molecules 26 01846 g001
Scheme 1. Mechanism of the photouncaging reaction of VBF protected by DMNB.
Scheme 1. Mechanism of the photouncaging reaction of VBF protected by DMNB.
Molecules 26 01846 sch001
Figure 2. Chemical structures of the DMNB-photocaged VFB pyrrole (NP) and sulfonamide (NS) amine groups considered in this work.
Figure 2. Chemical structures of the DMNB-photocaged VFB pyrrole (NP) and sulfonamide (NS) amine groups considered in this work.
Molecules 26 01846 g002
Figure 3. DFT/B3LYP/6-31+G(d)-optimized geometry of the DMNB-photocaged VFB and their corresponding photouncaging intermediates; uninvolved hydrogen atoms were omitted for clarity; highlighted regions indicate the center of the reaction; numbers indicate the bond length in Å.
Figure 3. DFT/B3LYP/6-31+G(d)-optimized geometry of the DMNB-photocaged VFB and their corresponding photouncaging intermediates; uninvolved hydrogen atoms were omitted for clarity; highlighted regions indicate the center of the reaction; numbers indicate the bond length in Å.
Molecules 26 01846 g003
Figure 4. DFT/B3LYP/6-31+G(d))-optimized geometry of the transition state (TS) of the photouncaging of 1 and 2; uninvolved hydrogen atoms are omitted for clarity; the given numbers correspond to the bond length of selected key bonds as indicated by arrows.
Figure 4. DFT/B3LYP/6-31+G(d))-optimized geometry of the transition state (TS) of the photouncaging of 1 and 2; uninvolved hydrogen atoms are omitted for clarity; the given numbers correspond to the bond length of selected key bonds as indicated by arrows.
Molecules 26 01846 g004
Figure 5. Time-dependent DFT (TD-DFT) (B3LYP/6-31 + G(d))-simulated UV–Vis absorption spectra of VFB in implicit aqueous solution; inset: simulated spectrum of DMNB under the same conditions; vertical lines indicate the position of the key electronic transition that is accountable for the corresponding absorption band.
Figure 5. Time-dependent DFT (TD-DFT) (B3LYP/6-31 + G(d))-simulated UV–Vis absorption spectra of VFB in implicit aqueous solution; inset: simulated spectrum of DMNB under the same conditions; vertical lines indicate the position of the key electronic transition that is accountable for the corresponding absorption band.
Molecules 26 01846 g005
Figure 6. The frontier molecular orbitals (Mos) of VFB, DMNB, 1, and 2; numbers in brackets refer to energy in hartree; uninvolved hydrogen atoms are omitted for clarity.
Figure 6. The frontier molecular orbitals (Mos) of VFB, DMNB, 1, and 2; numbers in brackets refer to energy in hartree; uninvolved hydrogen atoms are omitted for clarity.
Molecules 26 01846 g006
Figure 7. Superposition of NBOs of selected species; uninvolved hydrogen atoms are omitted for clarity.
Figure 7. Superposition of NBOs of selected species; uninvolved hydrogen atoms are omitted for clarity.
Molecules 26 01846 g007
Figure 8. Energy profile of the photouncaging reactions of 1 and 2. The Gibbs energy of the examined species was obtained from geometry optimization at the DFT/B3LYP /6-31G level of theory.
Figure 8. Energy profile of the photouncaging reactions of 1 and 2. The Gibbs energy of the examined species was obtained from geometry optimization at the DFT/B3LYP /6-31G level of theory.
Molecules 26 01846 g008
Table 1. Selected natural bond orbital (NBO) electronic transition of key elements (a) involved in the photouncaging reaction and their corresponding stabilization energies (E(2)).
Table 1. Selected natural bond orbital (NBO) electronic transition of key elements (a) involved in the photouncaging reaction and their corresponding stabilization energies (E(2)).
SpeciesTransition (Donor→Acceptor) E(2) kcal/molSpeciesTransition (Donor→Acceptor) E(2) kcal/mol
VFBlp(NP)→π*(C1-C2)42.02lp(NS)→σ*(S-O)10.2
lp(NP)→π*(C3-C4)35.9 lp(NS)→σ*(C5-C6)7.4
lp(NS)→σ*(S-O)8.2 lp(NS)→σ*(C5-C8)7.8
lp(NS)→σ*(C5-C6)7.8 lp(NS)→σ*(C′1-C′2)8.0
1lp(NP)→π*(C1-C2)41.32-I1lp(NS)→σ*(S-O)8.2
lp(NP)→π*(C3-C4)37.0 lp(NS)→σ*(C5-C6)6.3
lp(NP)→π*(C′1-C′2)6.7 lp(NS)→σ*(C5-C8)7.1
1-I1π(NP-C1)→π*(C′1-C′2)13.2 lp(NS)→σ*(C′1-C′2)33.8
π(C′3-N′4)→π*(C′1-C′2)5.82-I2lp(NS)→σ*(C′1-O′6)17.6
1-I2lp(NP)→σ*(C′1-O′6)16.1 lp(O′6)→σ*(NS-C′1)8.4
1-I3lp(NP)→σ*(C′1-O′6)5.62-I3lp(O′6)→σ*(NS-C′1)13.9
lp(O′6)→σ*(NP-C′1)8.52-TS1lp(NS)→σ*(C′1-C′2)26.1
lp(O′6)→s*(H)24.7 lp(O′6)→π*(C′1-C′2)9.4
1-TS1lp(NP)→π*(C′1-C′2)46.92-TS2lp(NS)→σ*(C′1-O′6)22.2
lp(O′6)→π*(C′1-C′2)6.72-TS3lp(NP)→σ*(O′6-H)7.1
1-TS2lp(NP)→σ*(C′1-O′6)16.4
1-TS3lp(NP)→σ*(O′6-H)6.6
(a) See Figure 1 for atom numbering.
Table 2. The NBO charges (a) of key atoms (b) of the intermediates and TSs of photouncaging reactions of prodrugs 1 and 2.
Table 2. The NBO charges (a) of key atoms (b) of the intermediates and TSs of photouncaging reactions of prodrugs 1 and 2.
Species12
NpNsC′1NpNsC′1
VFB −0.534−0.869
prodrug−0.359−0.849−0.282−0.534−0.667−0.281
I1−0.349−0.8360.019−0.533−0.598−0.005
TS1−0.336−0.8320.123−0.533−0.6000.057
I2−0.382−0.8360.231−0.534−0.6630.217
TS2−0.338−0.8420.247−0.533−0.6580.232
I3−0.384−0.8670.256−0.533−0.6870.242
TS3−0.591−0.8390.381−0.546−0.8320.361
(a) In units of electron (e); e = 1.6022 × 10−19 C. (b) See Figure 1 for atom numbering.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bani-Yaseen, A.D. Computational Study on the Mechanism of the Photouncaging Reaction of Vemurafenib: Toward an Enhanced Photoprotection Approach for Photosensitive Drugs. Molecules 2021, 26, 1846. https://doi.org/10.3390/molecules26071846

AMA Style

Bani-Yaseen AD. Computational Study on the Mechanism of the Photouncaging Reaction of Vemurafenib: Toward an Enhanced Photoprotection Approach for Photosensitive Drugs. Molecules. 2021; 26(7):1846. https://doi.org/10.3390/molecules26071846

Chicago/Turabian Style

Bani-Yaseen, Abdulilah Dawoud. 2021. "Computational Study on the Mechanism of the Photouncaging Reaction of Vemurafenib: Toward an Enhanced Photoprotection Approach for Photosensitive Drugs" Molecules 26, no. 7: 1846. https://doi.org/10.3390/molecules26071846

APA Style

Bani-Yaseen, A. D. (2021). Computational Study on the Mechanism of the Photouncaging Reaction of Vemurafenib: Toward an Enhanced Photoprotection Approach for Photosensitive Drugs. Molecules, 26(7), 1846. https://doi.org/10.3390/molecules26071846

Article Metrics

Back to TopTop