Computational modeling and target synthesis of monomethoxy-substituted o-diphenylisoxazoles with unexpectedly high antimitotic microtubule destabilizing activity.

The ability of monomethoxy-substituted o-diphenylisoxazoles 2a-d to interact with the colchicine site of tubulin was predicted using computational modeling, docking studies, and calculation of binding affinity. The respective molecules were synthesized in high yields by three steps reaction using easily available benzaldehydes, acetophenones, and arylnitromethanes as starting material. The calculated antitubulin effect was confirmed in vivo in a sea urchin embryo model. Compounds 2a and 2c showed high antimitotic microtubule destabilizing activity compared to that of CA4. Isoxazole 2a also exhibited significant cytotoxicity against human cancer cells in NCI60 screen. For the first time, isoxazole-linked CA4 derivatives 2a and 2c with only one methoxy substituent were identified as potent antimitotic microtubule destabilizing agents. These molecules could be considered as promising structures for further optimization.

The role of human in the loop: lessons from D3R challenge 4.

The rapid development of new machine learning techniques led to significant progress in the area of computer-aided drug design. However, despite the enormous predictive power of new methods, they lack explainability and are often used as black boxes. The most important decisions in drug discovery are still made by human experts who rely on intuitions and simplified representation of the field. We used D3R Grand Challenge 4 to model contributions of human experts during the prediction of the structure of protein-ligand complexes, and prediction of binding affinities for series of ligands in the context of absence or abundance of experimental data. We demonstrated that human decisions have a series of biases: a tendency to focus on easily identifiable protein-ligand interactions such as hydrogen bonds, and neglect for a more distributed and complex electrostatic interactions and solvation effects. While these biases still allow human experts to compete with blind algorithms in some areas, the underutilization of the information leads to significantly worse performance in data-rich tasks such as binding affinity prediction.

Quantifying Possible Routes for SpnF-Catalyzed Formal Diels-Alder Cycloaddition.

The Diels-Alder reaction is a cornerstone of modern organic synthesis. Despite this, it remains essentially inaccessible to biosynthetic approaches. Only a few natural enzymes catalyze even a formal [4 + 2] cycloaddition, and it remains uncertain if any of them proceed via the Diels-Alder mechanism. In this study, we focus on the [4 + 2] cycloaddition step in the biosynthesis of spinosyn A, a reaction catalyzed by SpnF enzyme, one of the most promising "true Diels-Alderase" candidates. The four currently proposed mechanisms (including the Diels-Alder one) for this reaction in water (as a first-order approximation of the enzymatic reaction) are evaluated by an exhaustive quantum mechanical search for possible transition states (728 were found in total). We find that the line between the recently proposed bis-pericyclic [J. Am. Chem. Soc. 2016, 138 (11), 3631] and Diels-Alder routes is blurred, and favorable transition states of both types may coexist. Application of the Curtin-Hammett principle, however, reveals that the bis-pericyclic mechanism accounts for ∼83% of the reaction flow in water, while the classical Diels-Alder mechanism contributes only ∼17%. The current findings provide a route for modeling this reaction inside the SpnF active site and inferring the catalytic architecture of possible Diels-Alderases.

An explicit account of solvation is essential for modeling Suzuki-Miyaura coupling in protic solvents.

We compared explicit and implicit solvation approaches in modeling the free energy profile of the final step of Suzuki-Miyaura coupling. Both approaches produced similar ΔG(≠) in all the studied solvents (benzene, toluene, DMF, ethanol, and water). Solvation free energies of individual reaction components reasonably correlated for explicit and implicit models in aprotic solvents (RMSE = 30-50 kJ mol(-1), R(2) > 0.71). However for ethanol and water the correlation was poor. We attributed this difference to the formation of the PdH-O hydrogen bond with Pd(PPh3)2 which was surprisingly observed in explicit modeling. Further QM calculations of the Pd(PPh3)2-H2O system confirmed the direction (PdH) and stability of this bonding. Therefore we stress the need for considering explicit solvation for modeling Pd-catalyzed reactions in protic solvents.