RESUMO
Unsaturated ketone derivatives are known as inhibitors of monoamine oxidase B (MAO-B), a potential drug target of Parkinson's disease. Here, docking-based alignment, 3 D-QSAR (three-dimensional quantitative structure-activity relationship) studies, ADMET (absorption, distribution, metabolism, excretion, and toxicity) prediction, molecular dynamics (MD) simulation, and MM_GBSA binding free energy were performed on a novel series of MAO-B inhibitors. The objective is to predict new MAO-B inhibitors with high potency activity. The 3 D-QSAR models were created using comparative molecular field analysis (CoMFA) and comparative molecular similarity index analysis (CoMSIA). Molecular docking findings indicated that compounds with strong inhibitory efficacy also had a high binding affinity. 3 D-QSAR studies showed the importance of steric, electrostatic, and H-bond acceptor fields on the inhibitory activity of MAO-B. Based on the appropriate 3 D-QSAR model, a new series of MAO-B inhibitors were predicted and their pharmacokinetic characteristics were evaluated using in silico ADMET prediction. All screened compounds show good oral bioavailability without any side effects. Moreover, the dynamic behavior and stability of the most active compounds were evaluated using MD simulations. The results showed that unsaturated ketone derivatives are stable and compact during the 100 ns of MD simulation. Finally, the binding free energy of complexes was determined using the MM_GBSA method; the findings indicated that the T1 compound is more stable (ΔGbinding = -409.506 KJ/mol) than the data set's highest active compound (ΔGbinding = -31.883 KJ/mol).Communicated by Ramaswamy H. Sarma.
Assuntos
Simulação de Dinâmica Molecular , Relação Quantitativa Estrutura-Atividade , Simulação de Acoplamento Molecular , Monoaminoxidase/química , Disponibilidade BiológicaRESUMO
Quantitative Structure Activity Relationship (QSAR) analysis techniques are tools largely utilized in many research fields, including drug discovery processes. In this work electronic descriptors are calculated with the Gaussian 03W software using the DFT method with the BecKe 3-parameters exchange functional and Lee-Yang-Parr correlation functional, with Kohn and Sham orbitals (KS) developed on a Gaussian Basis of type 6-31G (d), in combination with five Lipinski parameters that have been calculated with ChemOffice software, in order to develop a statistically verified 2D-QSAR model able to predict the biological activity of new molecules belonging to the same range of coumarins rather than chemical synthesis and biological evaluations that require more time and resources. Two QSAR models against both MCF-7 and HepG-2 cell lines are obtained using the multiple linear regression method. The predictive power of these models has been confirmed by internal and external validation. The Leverage method was used to determine the domain of applicability of the 2D-QSAR models developed. The results indicate that the best QSAR model is the one that links the 2D descriptors with the CDK inhibitory activity of the cell line (HepG-2) R2 = 0.748, R2cv = 0.618, MSE = 0.03 for the learning series and R2 = 0.73, MSE = 0.18 for the test series. This model implies that coumarin inhibitory activity is strongly related to dipole moment and the number of hydrogen bond donors. The results obtained suggest the importance of studying structure-activity relationships as a principal axis in drug design. The docking procedure using AutoDOCK Tools was also used to understand the mechanisms of molecular interactions and consequently, to develop new inhibitors.