Analysis of ligand binding to proteins using molecular dynamics simulations

Abstract

This work aims to explore theoretically the molecular mechanisms of ligand binding to proteins through the use of molecular dynamics simulations. The binding of sodium dodecyl sulfate (SDS) to cobra cardio toxin A3 (CTX A3) and thiourea (TOU) to lysozyme have been chosen as the two model systems. Data acquisitions were made by Gromacs software. To begin with, the collisions of ligand molecules with every residue of CTX A3 and lysozyme were evaluated. With this information in hand, the average numbers of collisions with each residue was defined and then assessed. Next, a measure of the affinity of a residue, P_i, referred to as conformational factor, toward a ligand molecule was established. Based on the results provided, all site-making residues for CTX A3 and lysozyme were identified. The results are in good agreement with the experimental data. Finally, based on this method, all site-making residues of bovine carbonic anhydrase (BCA) toward the SDS ligand were predicted.