Ab initio calculations are performed to elucidate the mechanism of the photoisomerization of azobenzene. We obtain the excitation energies of the $S_1{\mathrm{(n\to \pi }}^{*}\mathrm{),}$ $S_2{\mathrm{(\pi \to \pi }}^{*}\mathrm{),}$ and $S_3{\mathrm{(n}}^2{\mathrm{\to \pi }}^{\text{*2}})$ states by complete active space self-consistent field (CASSCF) and multireference single double configuration interaction (MRSDCI) calculations. Two-dimensional potential surfaces of the ground- and excited states are obtained at the CASSCF level in order to investigate the isomerization pathways. A conical intersection between the ground state and the $S_1$ state is found near the midpoint of the rotation pathway, and causes a radiationless transition. On the other hand, the $S_2$ state has local minima at the cis and trans structures, so that the isomerization proceeds at the $S_2$ surface following the deexcitation.

Alectinib has shown a greater efficacy to ALK-rearranged non-small-cell lung cancers in first-line setting; however, most patients relapse due to acquired resistance, such as secondary mutations in ALK including I1171N and G1202R. Although ceritinib or lorlatinib was shown to be effective to these resistant mutants, further resistance often emerges due to ALK-compound mutations in relapse patients following the use of ceritinib or lorlatinib. However, the drug for overcoming resistance has not been established yet. We established lorlatinib-resistant cells harboring ALK-I1171N or -G1202R compound mutations by performing ENU mutagenesis screening or using an in vivo mouse model. We performed drug screening to overcome the lorlatinib-resistant ALK-compound mutations. To evaluate these resistances in silico, we developed a modified computational molecular dynamic simulation (MP-CAFEE). We identified 14 lorlatinib-resistant ALK-compound mutants, including several mutants that were recently discovered in lorlatinib-resistant patients. Some of these compound mutants were found to be sensitive to early generation ALK-TKIs and several BCR-ABL inhibitors. Using our original computational simulation, we succeeded in demonstrating a clear linear correlation between binding free energy and in vitro experimental IC50 value of several ALK-TKIs to single- or compound-mutated EML4-ALK expressing Ba/F3 cells and in recapitulating the tendency of the binding affinity reduction by double mutations found in this study. Computational simulation revealed that ALK-L1256F single mutant conferred resistance to lorlatinib but increased the sensitivity to alectinib. We discovered lorlatinib-resistant multiple ALK-compound mutations and an L1256F single mutation as well as the potential therapeutic strategies for these ALK mutations. Our original computational simulation to calculate the binding affinity may be applicable for predicting resistant mutations and for overcoming drug resistance in silico. FUND: This work was mainly supported by MEXT/JSPS KAKENHI Grants and AMED Grants.