Abstract

Epidermal growth factor receptor (EGFR), a transmembrane glycoprotein, is overexpressed in many cancers such as head-neck, breast, prostate, and skin cancers for this reason it is a good target in cancer therapy and diagnosis. In nanobody-based cancer diagnosis and treatment, nanobodies with high affinity toward receptor (e.g. EGFR) results in effective treatment or diagnosis of cancer. In this regard, the main aim of this study is to develop a method based on molecular dynamic (MD) simulations for designing of 7D12 based nanobody with high affinity compared with wild-type nanobody. By surveying electrostatic and desolvation interactions between different residues of 7D12 and EGFR, the critical residues of 7D12 that play the main role in the binding of 7D12 to EGFR were elucidated and based on these residues, five logical variants were designed. Following the 50ns MD simulations, pull and umbrella sampling simulation were performed for 7D12 and all its variants in complex with EGFR. Binding free energy of 7D12 (and all its variants) with EGFR was obtained by weighted histogram analysis method. According to binding free energy results, GLY101 to GLU mutation showed the highest binding affinity but this variant is unstable after 50ns MD simulations. ALA100 to GLU mutation shows suitable binding enhancement with acceptable structural stability. Suitable position and orientation of GLU in residue 100 of 7D12 against related amino acids of EGFR formed some extra hydrogen and electrostatic interactions which resulted in binding enhancement.