Abstract
Integrin α(IIb)β(3) mediates platelet aggregation by binding the Arginyl-Glycyl-Aspartic acid (RGD) sequence of fibrinogen. RGD binding occurs at a site topographically proximal to the α(IIb) and β(3) subunits, promoting the conformational activation of the receptor from bent to extended states. While several experimental approaches have characterized RGD binding to α(IIb)β(3) integrin, applying computational methods has been significantly more challenging due to limited sampling and the need for a priori information regarding the interactions between the RGD peptide and integrin. In this study, we employed all-atom simulations using funnel metadynamics (FM) to evaluate the interactions of an RGD peptide with the α(IIb) and β(3) subunits of integrin. FM incorporates an external history-dependent potential on selected degrees of freedom while applying a funnel-shaped restraint potential to limit RGD exploration of the unbound state. Furthermore, it does not require a priori information about the interactions, enhancing the sampling at a low computational cost. Our FM simulations reveal significant molecular changes in the β(3) subunit of integrin upon RGD binding and provide a free-energy landscape with a low-energy binding mode surrounded by higher-energy prebinding states. The strong agreement between previous experimental and computational data and our results highlights the reliability of FM as a method for studying dynamic interactions of complex systems such as integrin.