Protein grabs a ligand by extending anchor residues: Molecular simulation for Ca²⁺ binding to calmodulin loop

¹ Kobe University

^* To whom correspondence should be addressed. E-mail: stakada{at}kobe-u.ac.jp.

Submitted on November 18, 2005
Revised on December 20, 2005
Accepted on 10 January 2006

	Abstract

The structural difference in proteins between unbound and bound forms directly suggests the importance of the conformational plasticity of proteins. However, pathways that connect two-end structures and how they are coupled to the binding reaction are not well understood at atomic resolution. Here, we analyzed the free energy landscape, explicitly taking into account coupling between binding and conformational change by performing atomistic molecular dynamics simulations for Ca²⁺ binding to a calmodulin loop. Using the AMBER force field with explicit water solvent, we conducted umbrella sampling for the free energy surface and steered molecular dynamics for the pathway search. We found that, at an early stage of binding, some key residue sidechains extend their "arms" to catch Ca²⁺ and, after catching, they carry the Ca²⁺ to the center of the binding pocket. This grabbing motion resulted in smooth and stepwise exchange in coordination partners of Ca²⁺ from water oxygen to atoms in the calmodulin loop. The key residue that first caught the ion was one of the two acidic residues, which are highly conserved. In the pathway simulations, different pathways were observed between binding and dissociation reactions: The former was more diverse than the latter.

Key Words: conformational change, free energy landscape, molecular dynamics, steered MD