Advances in Ligand Binding Predictions using Molecular Dynamics Simulations

Abstract: Biochemical processes all involve associations and dissociations of chemical entities. Understanding these is of substantial importance for many modern pharmaceutical applications. In this thesis, longstanding problems with regard to ligand binding are treated with computational methods, applied to proteins of key pharmaceutical importance. Homology modeling, docking, molecular dynamics simulations and free-energy calculations are used here for quantitative characterization of ligand binding to proteins. By combining computational tools, valuable contributions have been made for pharmaceutically relevant areas: a neglected tropical disease, an ion channel anti-drug-target, and GPCR drug-targets.We report three compounds inhibiting cruzain, the main cysteine protease of the protozoa causing Chagas’ disease. The compounds were found through an extensive virtual screening study and validated with experimental enzymatic assays. The compounds inhibit the enzyme in the ?M-range and are therefore valuable in further lead optimization studies.A high-resolution crystal structure of the BRICHOS domain is reported, together with molecular dynamics simulations and hydrogen-deuterium exchange mass spectrometry studies. This work revealed a plausible mechanism for how the chaperone activity of the domain may operate.Rationalization of structure-activity relationships for a set of analogous blockers of the hERG potassium channel is given. A homology model of the ion channel was used for docking compounds and molecular dynamics simulations together with the linear interaction energy method employed for calculating the binding free-energies.The three-dimensional coordinates of two GPCRs, 5HT1B and 5HT2B, were derived from homology modeling and evaluated in the GPCR Dock 2013 assessment. Our models were in good correlation with the experimental structures and all of them placed among the top quarter of all models assessed. Finally, a computational method, based on molecular dynamics free-energy calculations, for performing alanine scanning was validated with the A2A adenosine receptor bound to either agonist or antagonist. The calculated binding free-energies were found to be in good agreement with experimental data and the method was subsequently extended to non-alanine mutations. With extensive experimental mutation data, this scheme is a valuable tool for quantitative understanding of ligand binding and can ultimately be used for structure-based drug design.