Analysis of α4 β1integrin specific antagonists binding modes: structural insights by molecular docking, molecular dynamics and linear interaction energy method for free energy calculations

The VLA-4 antigen (α4β1 integrin) is involved in the pathophysiology of a variety of diseases including asthma, multiple sclerosis, rheumatoid arthritis and diabetes. The ligand selectivity toward this integrin remains a difficult problem, mainly due to the fact that 3D structures of most integrins are still unknown. We initially built a 3D computational model of the α4β1 ligand binding site, taking the crystal structure of the integrin αVβ3 as template. Then, we performed a computational study on a set of seven α4β1 antagonists, evaluating the binding modes of 4-[N'-(2-methylphenyl)ureido]phenylacetyl and derivatives by molecular docking. Molecular dynamics simulations were used to improve the receptor-ligand energy landscape exploration by the docking algorithm. The compounds were systematically arranged in two main binding modes, and in all cases, pointed out that these antagonists preferably bind to the α4β1 integrin active site in an extended conformation that resembles the one in solution. LIE (linear interaction energy) calculations also confirmed this statement given that the most prevailing binding mode is also the energetically most favored one. This study benefits the comprehension of the mechanism of this family of antagonists and may provide useful information for rational drug design.