Modeling disordered protein interactions from biophysical principles

<div><p>Disordered protein-protein interactions (PPIs), those involving a folded protein and an intrinsically disordered protein (IDP), are prevalent in the cell, including important signaling and regulatory pathways. IDPs do not adopt a single dominant structure in isolation but often become ordered upon binding. To aid understanding of the molecular mechanisms of disordered PPIs, it is crucial to obtain the tertiary structure of the PPIs. However, experimental methods have difficulty in solving disordered PPIs and existing protein-protein and protein-peptide docking methods are not able to model them. Here we present a novel computational method, IDP-LZerD, which models the conformation of a disordered PPI by considering the biophysical binding mechanism of an IDP to a structured protein, whereby a local segment of the IDP initiates the interaction and subsequently the remaining IDP regions explore and coalesce around the initial binding site. On a dataset of 22 disordered PPIs with IDPs up to 69 amino acids, successful predictions were made for 21 bound and 18 unbound receptors. The successful modeling provides additional support for biophysical principles. Moreover, the new technique significantly expands the capability of protein structure modeling and provides crucial insights into the molecular mechanisms of disordered PPIs.</p></div

L-RMSD vs Model Score and IDP RMSD.

<p>Inc: incorrect; Acc: acceptable; Med: medium. PDB ID: 2bzw.</p

IDP-LZerD consists of four steps.

<p>1. fragment structure prediction, 2. fragment docking, 3. path assembly, and 4. refinement. Steps 1 and 2 correspond to “dock” and Steps 3 and 4 correspond to “coalesce.”</p

Biological case studies.

<p><b>A</b>: <i>β</i>-catenin in complex with TCF7L2. Green: native TCF7L2; orange: rank 1 model of TCF7L2; <i>f</i>_<i>nat</i> 0.38, I-RMSD 2.85 Å, L-RMSD 7.94 Å. PDB ID: 1jpw. <b>B-E</b>: Human and mouse Cbp/p300 TAZ1 domain in complex with CITED2 and Hif-1<i>α</i>. Green/cyan: native CITED2/Hif-1<i>α</i>; orange/yellow: model CITED2/Hif-1<i>α</i>. Ball and stick: LPXL motif. <b>B-C</b>: Human TAZ1 and CITED2. <b>B</b>: bound (1p4qB); rank 5 model; <i>f</i>_<i>nat</i> 0.27, I-RMSD: 4.2 Å, L-RMSD: 7.6 Å. <b>C</b>: unbound (1l3eB); rank 9 model; <i>f</i>_<i>nat</i> 0.17, I-RMSD: 7.1 Å, L-RMSD: 9.6 Å. <b>D-E</b>: Mouse TAZ1 and Hif-1<i>α</i>. <b>D</b>: bound (1l8cA); rank 16 model; <i>f</i>_<i>nat</i> 0.05, I-RMSD 11.7 Å, L-RMSD 20.1 Å. <b>E</b>: unbound (1u2nA); rank 9 model; <i>f</i>_<i>nat</i> 0.20, I-RMSD 6.4 Å, L-RMSD 10.4 Å. <b>F</b>: Unbound complex between BoNT/A and sn2. Green: native sn2; orange: rank 1 model of sn2; <i>f</i>_<i>nat</i> 0.00, I-RMSD 15.7 Å, L-RMSD 38.2 Å. Receptor PDB ID 1xtfA (unbound).</p

Correlation between the docking score (DI) and the RMSD of the fragments.

<p>Data for sequence windows 1 and 2 of 4i9o. Green: top 30 docked fragments by DI score.</p

Complex between Bcl2-L-1 and BAD.

<p><b>(A)</b>: A model of the bound structure (2bzw) before (purple) and after (orange) refinement vs. native (green). <b>(B)</b>: Unbound (1pq0); blue-to-red (N-terminus on the left): native BAD; rainbow: top 7 models of BAD.</p

Fragment modeling and docking accuracy for 9-residue IDR complexes from ELM.

<p>Fragment modeling and docking accuracy for 9-residue IDR complexes from ELM.</p