Application of QSAR Method in the Design of Enhanced Antimalarial Derivatives of Azetidine-2-carbonitriles, their Molecular Docking, Drug-likeness, and SwissADME Properties

The resistance of the P. falciparum strain to some of the antimalarial drugs has been a dominant dilemma facing the treatment of this fetid disease. This necessitates the detection and development of new antimalarial agents targeting the P. falciparum. Azetidine-2-carbonitriles reported for its antimalarial activities, could provide an alternative to the customized antimalarial drugs. Leading to the use of quantitative structure-activity relationship (QSAR) studies, which relates the structures of Azetidine-2-carbonitriles with their activities to generate predictive models. The structures were optimized using density functional theory (DFT) DFT/B3LYP/6-31G* basis set to generate their molecular descriptors, where five predictive models were constructed using the generated descriptors. The models were constructed using the genetic function algorithm component of a material studio, where the model with good statistical parameters, high coefficient of determination (R(2)) = 0.9465, cross-validated R(2) (Q(2)cv) = 0.8981, Q(2 )((L4O))cv = 0.9272, and highest external validated R(2) (R(2)(pred)) = 0.6915 was selected as the best model. These statistical results show the robustness, excellent power of prediction, and validity of the selected model. The descriptor, SpMax2_Bhp (the maximum absolute eigenvalue of Barysz matrix for n = 2 was weighted by polarizability), was revealed to be the most influential in the model due to its highest mean effect. The descriptor played a role in the design of sixteen (16) theoretical derivatives of Azetidine-2-carbonitriles using compound 25 as the design template by increasing polarizability of the compounds through substitution of the various group with electron deactivating groups (F, I, Cl, SO(3)H, CN, NO(2), etc.) at different position of the template. The designed compounds were docked with Plasmodium falciparum dihydroorotate dehydrogenase (Pf-DHODH), giving compound D9 the highest binding energy. The designed compounds were further screened for their drug-likeness, where they all pass Lipinski’s RO5. All the compounds show good skin permeability coefficient and have low Gastrointestinal absorption while few compounds D1, D2, D3, D14, and D15 inhibiting the CYP1A2.