On the determinants of residence times and dissociation mechanisms of complexes of interleukin-13 with its low and high affinity receptors

Interleukin-13 (IL-13) is an immunomodulatory cell signaling cytokine that has been implicated in neurodegenerative disease and chronic inflammation. IL-13 binds to its low and high affinity receptors, IL-13 receptor 1 (IL-13R1) and IL-13 receptor 2 (IL-13R2), respectively, with residence times that vary accordingly. As the binding kinetics of the cytokine-receptor complexes influence cellular responses, we employed the molecular dynamics (MD) simulation-based{tau} -random acceleration molecular dynamics method ({tau}RAMD) to compute relative residence times for wild-type (WT) IL-13 and 19 IL-13 mutants to the two receptors. Comparison with experimental kinetic data shows that the{tau} RAMD computations capture the trends in residence times. Analysis of simulated dissociation trajectories of the cytokine-receptor complexes reveals two distinct dissociation pathways of IL-13 from each of the receptors. This study thus pinpoints key determinants of the interaction of IL-13 with its receptors which could be targeted for therapeutic design. Statement of SignificanceCytokines are regulatory proteins that bind to cell surface receptors and thereby send signals to the cellular interior. Interleukin-13 (IL-13) is a cytokine that has a low and a high affinity receptor. It has important physiological roles, and its deregulation is involved in diseases such as atopic dermatitis and asthma. We employed a molecular dynamics simulation-based method to compute the effects of changes in the sequence of IL-13 on the lifetimes of complexes of IL-13 and its receptors. Comparison with experiments supports the validity of the computational approach and analysis of the simulations reveals two distinct ways in which IL-13 dissociates from each receptor. These results thus provide a map for targeting IL-13 - receptor interactions for the design of therapeutics.