Allosteric enhancement of the BCR-Abl1 kinase inhibition activity of nilotinib by cobinding of asciminib.

Inhibitors that bind competitively to the ATP binding pocket in the kinase domain of the oncogenic fusion protein BCR-Abl1 are used successfully in targeted therapy of chronic myeloid leukemia (CML). Such inhibitors provided the first proof of concept that kinase inhibition can succeed in a clinical setting. However, emergence of drug resistance and dose-dependent toxicities limit the effectiveness of these drugs. Therefore, treatment with a combination of drugs without overlapping resistance mechanisms appears to be an appropriate strategy. In the present work, we explore the effectiveness of combination therapies of the recently developed allosteric inhibitor asciminib with the ATP-competitive inhibitors nilotinib and dasatinib in inhibiting the BCR-Abl1 kinase activity in CML cell lines. Through these experiments, we demonstrate that asciminib significantly enhances the inhibition activity of nilotinib, but not of dasatinib. Exploring molecular mechanisms for such allosteric enhancement via systematic computational investigation incorporating molecular dynamics, metadynamics simulations, and density functional theory calculations, we found two distinct contributions. First, binding of asciminib triggers conformational changes in the inactive state of the protein, thereby making the activation process less favorable by ∼4 kcal/mol. Second, the binding of asciminib decreases the binding free energies of nilotinib by ∼3 and ∼7 kcal/mol for the wildtype and T315I-mutated protein, respectively, suggesting the possibility of reducing nilotinib dosage and lowering risk of developing resistance in the treatment of CML.

Modulating the Photocyclization Reactivity of Diarylethenes through Changes in the Excited-State Aromaticity of the π-Linker.

Quantum chemical calculations are performed to explore if the reactivity of diarylethene switches toward photocyclization can be controlled by the excited-state aromaticity of their bridging π-linker. Using an archetypal diarylethene with a non-aromatic π-linker as a reference, completely different outcomes are found when the π-linker is allowed to become either aromatic (no reaction) or antiaromatic (fast reaction) upon photoexcitation. The results demonstrate a possibility to use the excited-state aromaticity concept for actual modulation of photochemical reactivity.

Photoswitching of Local (Anti)Aromaticity in Biphenylene-Based Diarylethene Molecular Switches.

Photoinduced tuning of (anti)aromaticity and associated molecular properties is currently in the focus of attention for both tailoring photochemical reactivity and designing new materials. Here, we report on the synthesis and spectroscopic characterization of diarylethene-based molecular switches embedded in a biphenylene structure composed of rings with different levels of local (anti)aromaticity. We show that it is possible to modulate and control the (anti)aromatic character of each ring through reversible photoswitching of the aryl units of the system between open and closed forms. Remarkably, it is shown that the irreversible formation of an annulated bis(dihydro-thiopyran) side-product that hampers the photoswitching can be efficiently suppressed when the aryl core formed by thienyl groups in one switch is replaced by thiazolyl groups in another.

Computational Comparison of Chemical and Isotopic Approaches to Control the Photoisomerization Dynamics of Light-Driven Molecular Motors.

Synthetic molecular motors driven by E/Z photoisomerization reactions are able to produce unidirectional rotary motion because of a structural asymmetry that makes one direction of rotation more probable than the other. In most such motors, this asymmetry is realized through the incorporation of a chemically asymmetric carbon atom. Here, we present molecular dynamics simulations based on multiconfigurational quantum chemistry to investigate whether the merits of this approach can be equaled by an alternative approach that instead exploits isotopic chirality. By first considering an N-methylpyrrolidine-cyclopentadiene motor design, it is shown that isotopically chiral variants of this design undergo faster photoisomerizations than a chemically chiral counterpart, while maintaining rotary photoisomerization quantum yields of similarly high magnitude. However, by subsequently considering a pyrrolinium-cyclopentene design, it is also found that the introduction of isotopic chirality does not provide any control of the directionality of the photoinduced rotations within this framework. Taken together, the results highlight both the potential usefulness of isotopic rather than chemical chirality for the design of light-driven molecular motors, and the need for further studies to establish the exact structural circumstances under which this asymmetry is best exploited.

Photoinduced Changes in Aromaticity Facilitate Electrocyclization of Dithienylbenzene Switches.

The concepts of excited-state aromaticity and antiaromaticity have in recent years with increasing frequency been invoked to rationalize the photochemistry of cyclic conjugated organic compounds, with the long-term goal of using these concepts to improve the reactivities of such compounds toward different photochemical transformations. In this regard, it is of particular interest to assess how the presence of a benzene motif affects photochemical reactivity, as benzene is well-known to completely change its aromatic character in its lowest excited states. Here, we investigate how a benzene motif influences the photoinduced electrocyclization of dithienylethenes, a major class of molecular switches. Specifically, we report on the synthesis of a dithienylbenzene switch where the typical nonaromatic, ethene-like motif bridging the two thienyl units is replaced by a benzene motif, and show that this compound undergoes electrocyclization upon irradiation with UV-light. Furthermore, through a detailed quantum chemical analysis, we demonstrate that the electrocyclization is driven jointly and synergistically by the loss of aromaticity in this motif from the formation of a reactive, antiaromatic excited state during the initial photoexcitation, and by the subsequent relief of this antiaromaticity as the reaction progresses from the Franck-Condon region. Overall, we conclude that photoinduced changes in aromaticity facilitate the electrocyclization of dithienylbenzene switches.

Light-driven rotary molecular motors without point chirality: a minimal design.

A fundamental requirement for achieving photoinduced unidirectional rotary motion about an olefinic bond in a molecular motor is that the potential energy surface of the excited state is asymmetric with respect to clockwise and counterclockwise rotations. In most available light-driven rotary molecular motors, such asymmetry is guaranteed by the presence of a stereocenter. Here, we present non-adiabatic molecular dynamics simulations based on multiconfigurational quantum chemistry to demonstrate that this chiral feature is not essential for inducing unidirectional rotary motion in molecules that incorporate a cyclohexenylidene moiety into a protonated Schiff-base framework. Rather, the simulations show that it is possible to exploit the intrinsic asymmetry of the puckered cyclohexenylidene to control the direction of photoinduced rotation.

Computational Insight to Improve the Thermal Isomerisation Performance of Overcrowded Alkene-Based Molecular Motors through Structural Redesign.

Synthetic overcrowded alkene-based molecular motors achieve 360° unidirectional rotary motion of one motor half (rotator) relative to the other (stator) through sequential photochemical and thermal isomerisation steps. In order to facilitate and expand the use of these motors for various applications, it is important to investigate ways to increase the rates and efficiencies of the reactions governing the rotary motion. Here, we use computational methods to explore whether the thermal isomerisation performance of some of the fastest available motors of this type can be further improved by reducing the sizes of the motor halves. Presenting three new redesigned motors that combine an indanylidene rotator with a cyclohexadiene, pyran or thiopyran stator, we first use multiconfigurational quantum chemical methods to verify that the photoisomerisations of these motors sustain unidirectional rotary motion. Then, by performing density functional calculations, we identify both stepwise and concerted mechanisms for the thermal isomerisations of the motors and show that the rate-determining free-energy barriers of these processes are up to 25 kJ mol-1 smaller than those of the original motors. Furthermore, the thermal isomerisations of the redesigned motors proceed in fewer steps. Altogether, the results suggest that the redesigned motors are useful templates for improving the thermal isomerisation performance of existing overcrowded alkene-based motors.

Computational design of faster rotating second-generation light-driven molecular motors by control of steric effects.

We report a systematic computational investigation of the possibility to accelerate the rate-limiting thermal isomerizations of the rotary cycles of synthetic light-driven overcrowded alkene-based molecular motors through modulation of steric interactions. Choosing as a reference system a second-generation motor known to accomplish rotary motion in the MHz regime and using density functional theory methods, we propose a three-step mechanism for the thermal isomerizations of this motor and show that variation of the steric bulkiness of the substituent at the stereocenter can reduce the (already small) free-energy barrier of the rate-determining step by a further 15-17 kJ mol(-1). This finding holds promise for future motors of this kind to reach beyond the MHz regime. Furthermore, we demonstrate and explain why one particular step is kinetically favored by decreasing and another step is kinetically favored by increasing the steric bulkiness of this substituent, and identify a possible back reaction capable of impeding the rotary rate.

How method-dependent are calculated differences between vertical, adiabatic, and 0-0 excitation energies?

Through a large number of benchmark studies, the performance of different quantum chemical methods in calculating vertical excitation energies is today quite well established. Furthermore, these efforts have in recent years been complemented by a few benchmarks focusing instead on adiabatic excitation energies. However, it is much less well established how calculated differences between vertical, adiabatic and 0-0 excitation energies vary between methods, which may be due to the cost of evaluating zero-point vibrational energy corrections for excited states. To fill this gap, we have calculated vertical, adiabatic, and 0-0 excitation energies for a benchmark set of molecules covering both organic and inorganic systems. Considering in total 96 excited states and using both TD-DFT with a variety of exchange-correlation functionals and the ab initio CIS and CC2 methods, it is found that while the vertical excitation energies obtained with the various methods show an average (over the 96 states) standard deviation of 0.39 eV, the corresponding standard deviations for the differences between vertical, adiabatic, and 0-0 excitation energies are much smaller: 0.10 (difference between adiabatic and vertical) and 0.02 eV (difference between 0-0 and adiabatic). These results provide a quantitative measure showing that the calculation of such quantities in photochemical modeling is well amenable to low-level methods. In addition, we also report on how these energy differences vary between chemical systems and assess the performance of TD-DFT, CIS, and CC2 in reproducing experimental 0-0 excitation energies.

Molecular dynamics and absolute binding free energy studies of piperine derivatives as potential inhibitors of SARS-CoV-2 main protease.

The work presents a library of piperine derivatives as potential inhibitors of the main protease protein (Mpro) functionality using Docking Studies, Molecular Dynamics (MD) Simulations and Absolute Binding Free-Energy calculations. 342 ligands were selected for this work and docked with Mpro protein. Among all the ligands studied, PIPC270, PIPC299, PIPC252, PIPC63, PIPC311 were the top five docked conformations having significant hydrogen bonding and hydrophobic interactions inside the active pocket of Mpro. These top five ligands were subjected to MD simulations for 100 ns using GROMACS. Root Mean Square Deviation (RMSD), Root Mean Square Fluctuation (RMSF), Radius of Gyration (Rg), Solvent Accessible Surface Area (SASA) and hydrogen bond analysis revealed that the ligands bounded to protein remain stable without significant deviations during the course of MD simulations. Absolute binding free energy (ΔGb) was calculated for theses complexes and found that the ligand PIPC299 shows the prevalent binding affinity with binding free-energy of about -113.05 Kcal/mol. Thus, these molecules can be further tested by in vitro and in vivo studies on Mpro. This study lays a path to explore the new functionality of piperine derivatives as novel drug like molecules.Communicated by Ramaswamy H. Sarma.

Activation of Abl1 Kinase Explored Using Well-Tempered Metadynamics Simulations on an Essential Dynamics Sampled Path.

Well-tempered metadynamics (wT-metaD) simulations using path collective variables (CVs) have been successfully applied in recent years to explore conformational transitions in protein kinases and other biomolecular systems. While this methodology has the advantage of describing the transitions with a limited number of predefined path CVs, it requires as an input a reference path connecting the initial and target states of the system. It is desirable to automate the path generation using approaches that do not rely on the choice of geometric CVs to describe the transition of interest. To this end, we developed an approach that couples essential dynamics sampling with wT-metaD simulations. We used this newly developed procedure to explore the activation mechanism of Abl1 kinase and compute the associated free energy barriers. Through these simulations, we identified a three-step mechanism for the activation that involved two metastable intermediates that possessed a partially open activation loop and differed primarily in the "in" or "out" conformation of the aspartate residue of the DFG motif. One of these states is similar to a conformation that was detected in previous spectroscopic studies of Abl1 kinase, albeit its mechanistic role in the activation was hitherto not well understood. The present study establishes its intermediary role in the activation and predicts a rate-determining free energy barrier of 13.8 kcal/mol that is in good agreement with previous experimental and computational estimates. Overall, our study demonstrates the usability of essential dynamics sampling as a path CV in wT-metaD to conveniently study conformational transitions and accurately calculate the associated barriers.

Unidirectional Rotary Motion in Isotopically Chiral Molecular Motors: A Computational Analysis.

Molecular dynamics simulations are performed to explore if isotopic chirality can induce unidirectional rotary motion in molecular motors operated through double-bond photoisomerizations. Using a high-quantum yield motor featuring a chemically asymmetric carbon atom as reference, it is found that isotopically chiral counterparts of this motor sustain such motion almost equally well. Overall, the study reveals a previously unexplored role for isotopic chirality in the design of rotary molecular motors.

Molecular Photoswitching Aided by Excited-State Aromaticity.

Central to the development of optoelectronic devices is the availability of efficient synthetic molecular photoswitches, the design of which is an arena where the evolving concept of excited-state aromaticity (ESA) is yet to make a big impact. The aim of this minireview is to illustrate the potential of this concept to become a key tool for the future design of photoswitches. The paper starts with a discussion of challenges facing the use of photoswitches for applications and continues with an account of how the ESA concept has progressed since its inception. Then, following some brief remarks on computational modeling of photoswitches and ESA, the paper describes two different approaches to improve the quantum yields and response times of switches driven by E/Z photoisomerization or photoinduced H-atom/proton transfer reactions through simple ESA considerations. It is our hope that these approaches, verified by quantum chemical calculations and molecular dynamics simulations, will help stimulate the application of the ESA concept as a general tool for designing more efficient photoswitches and other functional molecules used in optoelectronic devices.

Excited-State Aromaticity Improves Molecular Motors: A Computational Analysis.

A new approach to the design of more efficient light-driven rotary molecular motors is presented and evaluated computationally based on molecular dynamics simulations. The approach involves enabling part of the motor to become aromatic in the photoactive excited state, and is found to sharply increase the rotary quantum yields of the photoisomerizations that underlie the motor function. Excited-state aromaticity thus holds promise as a guiding principle toward better-performing molecular motors.

On the possibility to accelerate the thermal isomerizations of overcrowded alkene-based rotary molecular motors with electron-donating or electron-withdrawing substituents.

We employ computational methods to investigate the possibility of using electron-donating or electron-withdrawing substituents to reduce the free-energy barriers of the thermal isomerizations that limit the rotational frequencies achievable by synthetic overcrowded alkene-based molecular motors. Choosing as reference systems one of the fastest motors known to date and two variants thereof, we consider six new motors obtained by introducing electron-donating methoxy and dimethylamino or electron-withdrawing nitro and cyano substituents in conjugation with the central olefinic bond connecting the two (stator and rotator) motor halves. Performing density functional theory calculations, we then show that electron-donating (but not electron-withdrawing) groups at the stator are able to reduce the already small barriers of the reference motors by up to 18 kJ mol(-1). This result outlines a possible strategy for improving the rotational frequencies of motors of this kind. Furthermore, exploring the origin of the catalytic effect, it is found that electron-donating groups exert a favorable steric influence on the thermal isomerizations, which is not manifested by electron-withdrawing groups. This finding suggests a new mechanism for controlling the critical steric interactions of these motors. Graphical Abstract The introduction of electron-donating groups in one of the fastest rotary molecular motors known to date is found to reduce the free-energy barriers of the thermal steps that limit the rotational frequencies by up to 18 kJ mol(-1).