In order to contribute to the structural basis for rational design of calmodulin (CaM) inhibitors, we analyzed the interaction of CaM with 14 classic antagonists and two compounds that do not affect CaM, using docking and molecular dynamics (MD) simulations, and the data were compared to available experimental data. The Ca(2+)-CaM-Ligands complexes were simulated 20 ns, with CaM starting in the "open" and "closed" conformations. The analysis of the MD simulations provided insight into the conformational changes undergone by CaM during its interaction with these ligands. These simulations were used to predict the binding free energies (ΔG) from contributions ΔH and ΔS, giving useful information about CaM ligand binding thermodynamics. The ΔG predicted for the CaM's inhibitors correlated well with available experimental data as the r(2) obtained was 0.76 and 0.82 for the group of xanthones. Additionally, valuable information is presented here: I) CaM has two preferred ligand binding sites in the open conformation known as site 1 and 4, II) CaM can bind ligands of diverse structural nature, III) the flexibility of CaM is reduced by the union of its ligands, leading to a reduction in the Ca(2+)-CaM entropy, IV) enthalpy dominates the molecular recognition process in the system Ca(2+)-CaM-Ligand, and V) the ligands making more extensive contact with the protein have higher affinity for Ca(2+)-CaM. Despite their limitations, docking and MD simulations in combination with experimental data continue to be excellent tools for research in pharmacology, toward a rational design of new drugs.
Keywords: MMPBSA; calmodulin; conformational flexibility; docking; molecular dynamics; molecular recognition; protein-ligand interactions; protein-peptide interactions.