# Intrinsic dynamics study identifies two amino acids of TIMP-1 critical for its LRP-1-mediated endocytosis in neurons

Laurie Verzeaux, Nicolas Belloy, Jessica Thevenard-Devy, Jérome Devy, Géraldine Ferracci, Laurent Martiny, Stéphane Dedieu, Manuel Dauchez, Hervé Emonard, Nicolas Etique, et al.

<figure><img src="image_1.png" /></figure>

## To cite this version:

Laurie Verzeaux, Nicolas Belloy, Jessica Thevenard-Devy, Jérome Devy, Géraldine Ferracci, et al. Intrinsic dynamics study identifies two amino acids of TIMP-1 critical for its LRP-1-mediated endocytosis in neurons. Scientific Reports, 2017, 7 (1), pp.5375. 10.1038/s41598-017-05039-z. hal-02184512

HAL Id: hal-02184512

https://hal.univ-reims.fr/hal-02184512

Submitted on 13 Dec 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

OPEN

# Intrinsic dynamics study identifies two amino acids of TIMP-1 critical for its LRP-1-mediated endocytosis in neurons

Laurie Verzeaux<fnref n="1" />, Nicolas Belloy<fnref n="1" />,², Jessica Thevenard-Devy<fnref n="1" />, Jérôme Devy<fnref n="1" />, Géraldine Ferracci³, Laurent Martiny<fnref n="1" />, Stéphane Dedieu<fnref n="1" />, Manuel Dauchez<fnref n="1" />,², Hervé Emonard<fnref n="1" />, Nicolas Etique<fnref n="1" /> & Emmanuelle Devarenne-Charpentier<fnref n="1" />

The tissue inhibitor of metalloproteinases-1 (TIMP-1) exerts inhibitory activity against matrix metalloproteinases and cytokine-like effects. We previously showed that TIMP-1 reduces neurite outgrowth in mouse cortical neurons and that this cytokine-like effect depends on TIMP-1 endocytosis mediated by the low-density lipoprotein receptor-related protein-1 (LRP-1). To gain insight into the interaction between TIMP-1 and LRP-1, we considered conformational changes that occur when a ligand binds to its receptor. TIMP-1 conformational changes have been studied using biomolecular simulations, and our results provide evidence for a hinge region that is critical for the protein movement between the N- and C-terminal TIMP-1 domains. *In silico* mutants have been proposed on residues F12 and K47, which are located in the hinge region. Biological analyses of these mutants show that F12A or K47A mutation does not alter MMP inhibitory activity but impairs the effect of TIMP-1 on neurite outgrowth. Interestingly, these mutants bind to LRP-1 but are not endocytosed. We conclude that the intrinsic dynamics of TIMP-1 are not involved in its binding to LRP-1 but rather in the initiation of endocytosis and associated biological effects.

The tissue inhibitor of metalloproteinase-1 (TIMP-1) is a natural inhibitor of matrix metalloproteinases (MMP) and several a disintegrin and metalloproteinases (ADAM)<fnref n="1" />. Due to its inhibitory functions, TIMP-1 is largely involved in the control of extracellular matrix remodelling in both physiological and pathological conditions². Moreover, TIMP-1 has been widely depicted as a cytokine-like effector that triggers various cellular responses independently of its MMP inhibitory activity³⁻⁵. For instance, we recently showed that TIMP-1 decreased neurite outgrowth in cortical neurons and that this effect was mainly dependent on its endocytosis mediated by the low-density lipoprotein receptor-related protein-1 (LRP-1)⁶. LRP-1 exhibits important endocytic and signalling functions that regulate the behaviour of many cell types⁷. In neurons, LRP-1 is abundantly expressed and mediates the endocytosis of various extracellular ligands including TIMP-1⁸. High levels of TIMP-1 are secreted by hyperactive astrocytes, and TIMP-1 expression is highly correlated with various neurological diseases⁹,¹⁰. Characterising TIMP-1/LRP-1 interaction could thus be of physiological relevance in the treatment of certain neurodegenerative disorders.

Molecular docking is the traditional method for predicting how a ligand binds a receptor¹¹, ¹². Nevertheless, this method requires the determination of both partner structures, which is usually done by X-ray diffraction or nuclear magnetic resonance (NMR) spectroscopy. The well-characterised structure of TIMP-1 consists of six disulfide-bonded loops forming two structurally distinct domains, and the three-dimensional TIMP-1 structure has a wedge-shaped appearance¹³, ¹⁴. The N-terminal domain (N-TIMP-1), which is composed of 126 amino

Received: 4 October 2016

Accepted: 24 May 2017

Published online: 14 July 2017

<footnote n="1">¹CNRS UMR 7369: Matrice Extracellulaire et Dynamique Cellulaire (MEDyC), UFR Sciences Exactes et Naturelles, Université de Reims Champagne-Ardenne (URCA), Laboratoire SiRMa - Campus Moulin de la Housse, BP 1039, 51687, Reims cedex, France. ²Plate-forme de Modélisation Moléculaire Multi-échelle (P3M), Université de Reims Champagne-Ardenne, Reims, France. ³Aix-Marseille Université, CNRS, Centre de Recherche en Neurobiologie et Neurophysiologie de Marseille (CRN2M), UMR 7286, Plate-forme de Recherche en Neurosciences (PFRN), Marseille, France. Nicolas Etique and Emmanuelle Devarenne-Charpentier contributed equally to this work. Correspondence and requests for materials should be addressed to N.E. (email: nicolas.etique@univ-reims.fr)</footnote>
acids, carries the inhibitory activity against MMPs by forming a non-covalent 1:1 stoichiometric complex with the proteinase¹⁵. The C-terminal domain (C-TIMP-1), which is composed of 58 amino acids, is structurally less characterised, but it has been shown to interact with the proMMP-9 hemopexin domain¹⁴.

LRP-1 is a large receptor composed of a long extracellular α-chain non-covalently associated to a short transmembrane β-chain⁷. The α-chain contains four ligand-binding domains composed of cysteine-rich complement-type repeats, and domains II and IV are the major binding regions, interacting with more than forty ligands⁷, ¹⁶. Unfortunately, the high molecular mass of LRP-1 and the presence of a predicted unordered region are major obstacles to elucidate the entire LRP-1 structure. Consequently, the use of molecular docking tools is not appropriate for studying the TIMP-1/LRP-1 interaction.

Protein conformational changes defined by protein flexibility and dynamics play a crucial role in ligand/receptor interaction¹⁷⁻¹⁹. We have thus hypothesised that the alteration of these properties could modify the TIMP-1/LRP-1 interaction and associated cellular effects. Protein dynamics can be evaluated in silico by normal mode analysis (NMA)²⁰, ²¹ and/or principal component analysis (PCA) of molecular dynamics (MD) simulations²². For instance, combining these approaches helped us to obtain reliable results consistent with the experimental data in the case of CD47/TSP-1 interaction by identifying large amplitude motions of the TSP-1 C-terminal domain²³, ²⁴. We have thus combined NMA and MD simulations with biological experiments to characterise the TIMP-1/LRP-1 interaction. NMA performed on the structure energy of TIMP-1 showed movement between the N- and C-terminal domains of TIMP-1 and indicated regions with high deformation energy and low carbon alpha atomic fluctuation. NMA and MD clearly pointed out that these regions are located in a hinge region that could be essential for protein movement. Interestingly, a single mutation of residue F12 or K47 (numbering of residues in the mature secreted protein) located in this region inhibits TIMP-1 cytokine-like activity in neurons but surprisingly does not alter TIMP-1 binding to LRP-1 domains II and IV. The data obtained using in silico simulations and biological experiments highlight the relevance of protein dynamics in the TIMP-1/LRP-1 interaction and associated biological effects.

## Results and Discussion

**Determination of TIMP-1 intrinsic dynamics using molecular modelling.** TIMP-1 has been described as a protein whose function could be controlled by its intrinsic dynamics²⁵. Protein dynamics, which represent intrinsic subregional motions, could also be a factor in ligand binding to its receptor¹⁷, ¹⁸. We have thus hypothesised that alteration of these protein dynamics could modify the TIMP-1/LRP-1 interaction. We first studied TIMP-1 intrinsic molecular motion using NMA. The energy of the TIMP-1 structure (PDB 1UEA) after the addition of a hydrogen atom was minimised by successively combining the Steepest Descent and Adopted Basis Newton-Raphson (ABNR) methods. The NMA of this structure generated a set of conformations. We have excluded the first six modes characterised by a null frequency corresponding to rigid-body rotation and translation movements²⁶. Of the next following modes, modes 7 to 9 were characterised by the lowest frequencies corresponding to large movements, which are the most functionally relevant for proteins²⁶ (Fig. 1a). Moreover, large region displacements are associated with a high degree of collectivity¹⁹. Taking into account the intermediate structures of TIMP-1 along the three modes, the flexibility profile along mode 7 exhibited the highest degree of collectivity (numerous arrows with similar lengths). This mode, which is presented in the supplementary Movie 1, showed large interdomain displacements comparable to a “wirecutter” movement between the N- and C-terminal domains.

Along mode 7, we determined the associated deformation energy and alpha carbon (Cα) atomic fluctuation to identify the TIMP-1 regions that govern this movement. Deformation energy analysis was used to measure the amount of local flexibility in the protein structure (i.e. atomic motion relative to neighbouring atoms)²⁷, ²⁸. High deformation indicates flexible regions, and our data showed that residues with the highest deformation energy (red bar) were present throughout the primary structure (Fig. 1b). In the TIMP-1 three-dimensional (3D) structure, these residues are located between the N- and C-terminal domains (Fig. 1c). Cα atomic fluctuations provide amplitudes of the absolute atomic motion²⁸, and our data (Fig. 1d) showed that residues with high fluctuation (red bar) were dispersed in the TIMP-1 primary structure but located outside the 3D structure (Fig. 1d and e). Interestingly, regions with the highest deformation energy seemed to be localised in the regions exhibiting the lowest atomic fluctuation (Fig. 1c and e). To confirm this observation, the deformation energy and atomic fluctuation data have been normalised and then superposed (Fig. 2a). The graph shows nine regions (green arrow) in which the residues exhibited both high deformation energy and low Cα atomic fluctuation. Interestingly, these regions were grouped in an area of the TIMP-1 3D structure that could be linked to a hinge essential for the interdomain displacements (Fig. 2b). Among the residues present in the nine regions mentioned above, we sought to identify those with mutations that could disturb the movement. Here, we excluded the amino acids present in secondary structures in order to prevent protein misfolding (Fig. 2c). Residues K47 and F12 of the mature secreted protein seemed to be of particular interest because the K47 side chain forms a hydrogen bond with the carbonyl oxygen of the F12 backbone (distance < 3.2 Å). Moreover, this bond was maintained in all intermediate TIMP-1 structures along mode 7 (Fig. 2d), suggesting that these two amino acids could contribute to chain flexibility and N- and C-terminal interdomain displacement. To support this finding, a 100 ns MD simulation was carried out and subjected to PCA to identify collective motions after fitting the trajectory on the first reference frame (Fig. 3). Each PCA showed that F12 and K47 are located outside regions displaying large atomic displacements. Moreover, F12 and K47 along the third PC made the lowest contribution to the total atomic displacement. This additional analysis confirmed that F12 and K47 could be of interest for TIMP-1 movement. We hypothesised that disrupting the interaction between F12 and K47 could affect the interdomain movements of TIMP-1 and its subsequent biological activity.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Normal mode analysis of TIMP-1 structure. (a) Intrinsic dynamic of TIMP-1 along modes 7, 8 and 9. Open structures of TIMP-1 are in the upper part of the panel and closed structures are in the lower part. Arrows (blue and red) indicate the direction and the deviation (length of arrows) of each residue returning to the 3D reference structure of TIMP-1. (b) Deformation energy of each TIMP-1 residue (numbering of the mature secreted protein) along mode 7. The lowest deformation energies are coloured blue and the highest are red. (c) Visualisation of the TIMP-1 3D-structure of each amino acid deformation energy. Residues with the lowest deformation energy are thin and coloured blue and those with the highest atomic fluctuation are thick and coloured red. (d) Cα atomic fluctuation of each TIMP-1 residue (numbering of the mature secreted protein) along mode 7. The lowest atomic fluctuations are coloured blue and the highest are red. (e) Visualisation of the TIMP-1 3D-structure of each amino acid atomic fluctuation. Residues with the lowest atomic fluctuations are thin and coloured blue and those with the highest atomic fluctuations are thick and coloured red.</figcaption></figure>

**Generation of TIMP-1 mutants.** To disrupt the interaction between F12 and K47, these residues were mutated into alanine. The choice of alanine was based on its low steric hindrance and/or its apolar character. This mutation should prevent hydrogen bond formation between K47 and F12 and theoretically disturb the movement described by mode 7. The two mutated TIMP-1 sequences have been inserted into the p3X-FLAG-CMV-14
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Identification of residues involved in TIMP-1 interdomain movements. (a) Superposition of the deformation energy (dotted line) and atomic fluctuation (coloured bars) data from TIMP-1 residues (numbering of the mature secreted protein) along mode 7. The green arrows point out residues or sets of residues with both high deformation energy and low $C\alpha$ atomic fluctuation. (b) Localisation in the TIMP-1 3D structure of the residues identified in (a) with both high deformation energy and low $C\alpha$ atomic fluctuation. These residues are in an area surrounded by a dotted line. (c) Localisation in the TIMP-1 secondary structure of the residues identified in (a) with both high deformation energy and low $C\alpha$ atomic fluctuation (green arrows). $\alpha$ helices are defined by the purple area and $\beta$ sheets by the yellow area. (d) Left: localisation in the TIMP-1 3D structure of F12 and K47 residues. Right: hydrogen bond formed between the carbonyl oxygen of the F12 backbone and the K47 side chain.</figcaption></figure>

vector to produce proteins with a C-terminal 3X-FLAG tag, facilitating their purification (Fig. 4a). Each mutant (T1-F12A and T1-K47A) and a wild-type TIMP-1 (T1-WT) were produced in Chinese hamster ovary (CHO) cells to be correctly folded and glycosylated. Western blot analysis using antibodies directed against TIMP-1 or 3X-FLAG showed that T1-WT and TIMP-1 mutants exhibited migration profiles that correspond to the expected 37 kDa molecular mass of the 3X-FLAG glycosylated proteins (Fig. 4b, left panel). Non-transfected CHO cells did not express endogenous TIMP-1 (data not shown) while CHO transfected cells expressed and secreted recombinant T1-WT, T1-K47A or T1-F12A with homogenous expression levels (Fig. 4b, right panel). To ensure proper folding of the recombinant proteins, we tested their ability to inhibit MMP-1, -2, -3 and -9. Indeed, it has been demonstrated that correct folding of TIMP-1 is necessary for its MMP inhibitory activity and that a single point mutation in the TIMP-1 sequence may alter this activity¹⁴, ²⁵, ²⁹. Our data showed that the Ki values exhibited by mutants were similar to those of T1-WT (Table 1).
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Principal component analysis derived from molecular dynamics simulation of TIMP-1. For each PC, the contribution of each residue to the overall atomic displacement was determined. The green arrows point out residues F12 and K47.</figcaption></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /><figcaption>**Figure 4.** Production of TIMP-1 mutants. (a) Schematic representation of T1-WT, T1-F12A and T1-K47A proteins. (b) A representative western blot with the indicated antibodies of cell extracts and conditioned medium (cond. med) of stably transfected CHO cells with p3X-FLAG-CMV-14 expressing T1-WT, T1-F12A or T1-K47A. The relative protein level based on three independent experiments is indicated under each band. The gels were run under the same experimental conditions and were shown as cropped gels/blots (Cropped gels/blots are shown in Supplementary Figure S1).</figcaption></figure>

<table><caption>**Table 1.** Inhibition constants (Ki) of TIMP-1 mutants against MMP-1, MMP-2, MMP-3 and MMP-9. Mean values ± standard error are presented. Data are based on three independent experiments.</caption><thead><tr><th rowspan="2"></th><th colspan="4">Ki (nM)</th></tr><tr><th>MMP-1</th><th>MMP-2</th><th>MMP-3</th><th>MMP-9</th></tr></thead><tbody><tr><td>T1-WT</td><td>4.03 ± 1.57</td><td>0.09 ± 0.04</td><td>1.18 ± 0.47</td><td>4.25 ± 0.28</td></tr><tr><td>T1-F12A</td><td>3.82 ± 1.33</td><td>0.06 ± 0.02</td><td>0.58 ± 0.60</td><td>4.76 ± 0.57</td></tr><tr><td>T1-K47A</td><td>3.38 ± 2.92</td><td>0.03 ± 0.01</td><td>0.55 ± 0.19</td><td>4.40 ± 0.78</td></tr></tbody></table>
These data indicate that individual single point mutations (K47A and F12A) do not modify the inhibitory activity of TIMP-1 towards these MMPs and confirms that the folding of the two TIMP-1 mutants was correct.

**Effect of the TIMP-1 mutants on neurite outgrowth.** We then evaluated the effect of F12A and K47A mutations on the cytokine-like TIMP-1 activity in mouse cortical neurons. Indeed, we and others have reported that TIMP-1 directly bound to LRP-1 at the plasma membrane of mouse cortical neurons, leading to cell morphology changes<sup>6,30</sup>. We first evaluated the effects of TIMP-1 mutants on neuron morphology in similar experimental conditions, as previously described<sup>6</sup>. Cortical neurons treated or not treated with T1-WT, T1-F12A or T1-K47A were labelled with anti-βIII tubulin antibodies to visualise the microtubule cytoskeleton and to measure neurite outgrowth (Fig. 5a). As expected, T1-WT significantly reduced the neurite length after 30 min treatment (Fig. 5a and b). In contrast, T1-F12A and T1-K47A exerted no effect on cortical neuron morphology (Fig. 5a and b). To explain this result, it was hypothesised that the TIMP-1 mutant binding to LRP-1 was impaired. To test this, we evaluated the ability of T1-F12A and T1-K47A to interact with LRP-1 domains II and IV by surface plasmon resonance (SPR). Surprisingly, this analysis demonstrated a direct interaction between T1-F12A or T1-K47A and LRP-1 domains II (DII) and IV (DIV), with a nanomolar range of affinity similar to that of T1-WT (Table 2). T1-F12A and T1-K47A were associated with DII about 10 times faster than with DIV, while they dissociated with DII slower compared to DIV (Table 2). The measured association ($k_{on}$) and dissociation ($k_{off}$) rates were similar to those previously obtained for wild-type TIMP-1<sup>6</sup>. Our SPR data demonstrated that the lack of neurite outgrowth inhibition could not be explained by a difference of affinity towards LRP-1 between the mutants and T1-WT.

To support these unexpected results, we performed a competition assay by co-treating cortical neurons with different T1-WT/T1-F12A or T1-WT/T1-K47A ratios (Fig. 5c). The equimolar ratio of T1-WT/T1-F12A or T1-WT/T1-K47A reduced the effect of T1-WT on neurite length by about twofold. Moreover, a five-fold excess of mutants compared to TIMP-1 completely abrogated the TIMP-1-WT effect. These data indicated that wild-type TIMP-1 and its mutants interacted with identical, or at least spatially similar, amino acids of the LRP-1 α-chain. Intriguingly, although both TIMP-1 mutants bound to LRP-1, F12A and K47A mutations prevented the biological effects associated with the TIMP-1/LRP-1 interaction.

**Study of LRP-1-mediated endocytosis of the TIMP-1 mutants.** LRP-1-mediated endocytosis and signalling are associated events<sup>31</sup>. Since TIMP-1 mutants bind to LRP-1 without triggering a biological effect in neurons, we investigated whether these mutants were internalised by LRP-1 utilising two experimental methods. First, we used eFluor®-conjugated T1-WT (fluoT1-WT), T1-F12A (fluoT1-F12A) and T1-K47A (fluoT1-K47A) according to a previously reported experimental protocol (6) to discriminate cell surface bound from internalised proteins (Fig. 6). At 4 °C, when endocytosis was blocked, we observed that the fluoT1-WT bound to the cell surface of the neurons. This binding was greatly inhibited by RAP, a potent antagonist of ligand binding to LRP receptors<sup>32,33</sup> or R2629, a previously validated LRP-1 blocking antibody<sup>34</sup>. The fluorescence at the cell surface was quite similar at about 1500 A.U. for the T1-WT and mutants and confirmed the SPR data (Table 2). RAP similarly impeded the binding of the two eFluor®-conjugated TIMP-1 mutants to the cell surface of neurons. The very low fluorescence intensities measured at 4 °C in the intracellular compartment are considered to be nonspecific. At 37 °C, when endocytosis was initiated, about 30% of the fluoT1-WT bound to the surface of neurons was internalised. Interestingly, the fluorescence intensity at the cell surface remained equal at 4 °C and 37 °C for the two mutants. This result suggests that these mutants were not endocytosed by LRP-1.

To support these data, we then examined the T1-F12A and T1-K47A internalisation by confocal microscopy imaging (Fig. 7). For this purpose, cortical neurons were incubated with T1-WT, T1-F12A or T1-K47A for 1 h at 4 °C, and then transferred to 37 °C for 10 min to allow endocytosis. Since LRP-1 ligands are internalised in the endosomal compartments<sup>7,35</sup>, cortical neurons were stained with antibodies directed to TIMP-1 and EEA1, an early endosome membrane marker<sup>36</sup> (Fig. 7). Confocal imaging analysis revealed that EEA1 labelling was statistically constant in all treated cells (Fig. 8a). T1-WT accumulated inside the cells (Fig. 8b) and co-located with EEA1 (Fig. 8c), demonstrating that wild-type TIMP-1 was actively endocytosed in cortical neurons. RAP inhibited this effect, confirming that TIMP-1 endocytosis was mediated by LRP-1. Interestingly, TIMP-1 labelling was significantly lessened in mutant-treated cells as compared to T1-WT-treated cells (Fig. 8b). Moreover, TIMP-1/EEA1 colocalisation was also decreased in T1-F12A- or T1-K47A-treated cells (Fig. 8c).

Since LRP-1 has been described as a fast endocytic receptor<sup>37</sup>, a rapid cycle of TIMP-1 uptake, trafficking to the cell surface and exocytosis could explain our data. Since our SPR data show that the association and dissociation rates were similar between the T1-WT and mutants (Table 2), this hypothesis is barely conceivable. Another explanation could be a delay of TIMP-1 mutant internalisation. To test this, the endocytosis assay was extended up to one hour, and in this condition, the TIMP-1 mutants were not endocytosed (data not shown). Together, these data highlighted that F12A or K47A mutation prevents TIMP-1 endocytosis in primary cortical neurons.

Several studies have revealed that phosphorylation of the tyrosine residue in the NPXY motifs present in the β-chain intracellular tail represents a molecular mechanism for switching LRP-1 function from that of an endocytic receptor to that of a signalling receptor by modulating the class and type of adaptor proteins that will associate with these motifs<sup>38,39</sup>. The mechanism leading to the phosphorylation of YXXL and NPXY motifs following ligand binding is still poorly understood. Therefore, it is difficult to explain why TIMP-1 mutants that bind to LRP-1 are not internalised and do not trigger cellular effects in neurons.

Our biological results are in accordance with those of several studies showing that proteins may undergo conformational changes following their binding<sup>40-42</sup>. They confirm our hypothesis that disrupting the interaction between F12 and K47 affects TIMP-1 dynamics and its subsequent biological activity. It would be interesting to use other approaches to confirm that the intrinsic dynamics of F12A and K47A mutants are disturbed. For example, NMR spectroscopy is a powerful technique for measuring protein dynamics. Nevertheless, this technique is restricted by some practical limitations (time cost, quantity of material, protein size). Recent advancements
a

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

b

<figure><img src="image_12.png" /></figure>

c

<figure><img src="image_13.png" /></figure>

**Figure 5.** Effect of TIMP-1 mutants on neurite outgrowth. (a) Immunofluorescence of cortical neurons from mouse embryos cultured for 24 h on poly-L-lysine-coated coverslips and then treated for 30 min with PBS (CTRL) or 5 nM T1-WT, T1-F12A and T1K47A. Cells were labelled with anti-βIII-tubulin monoclonal antibody and observed by fluorescent microscopy. The images are representative of more than three distinct experiments. Scale bar: 10 µm. (b) Quantification of neurite mean length per cell performed using the ImageJ plugin NeuronJ. For each group, at least 20 randomly selected fields were analysed in three independent experiments (n = 20). The data are representative of at least three independent experiments with at least five mice embryos. Error bars indicate mean ± SD, *** indicates P-value < 0.001 and NS indicates a value that is not statistically significant. (c) Neurite mean length of cortical neurons from mouse embryos cultured for 24 h on poly-L-lysine-coated coverslips and then treated for 30 min with PBS (CTRL) or T1-WT or co-treated with different T1-WT/T1-F12A or T1-WT/T1-K47A ratios (0:1, 1:1 and 1:5). Cells were labelled with anti-βIII-tubulin monoclonal antibody, observed by fluorescent microscopy and analysed with the ImageJ plugin NeuronJ Quantification. For each group, at least 10 randomly selected fields were analysed (n = 10). The data are representative of at least three independent experiments with at least five mice embryos. Error bars indicate mean ± SEM, *** indicates P-value < 0.001 and NS indicates a value that is not statistically significant.
<table><caption>**Table 2.** Surface plasmon resonance data for direct binding of TIMP-1 mutants to immobilised LRP-1 domains II (DII) and IV (DIV). Data are based on three measurements using five different concentrations for each measurement. The equilibrium constants of dissociation ($K_D$) were calculated from the association ($k_{on}$) and dissociation ($k_{off}$) rate constants. Mean values ± standard error are presented.</caption><thead><tr><th></th><th>k<sub>on</sub> (0.106 M<sup>-1</sup>s<sup>-1</sup>)</th><th>k<sub>off</sub> (0.10<sup>-2</sup> s<sup>-1</sup>)</th><th>K<sub>D</sub> (nM)</th></tr></thead><tbody><tr><td colspan="4">DII</td></tr><tr><td>T1-WT</td><td>1.28 ± 0.35</td><td>1.36 ± 0.20</td><td>10.8 ± 2.1</td></tr><tr><td>T1-F12A</td><td>1.07 ± 0.25</td><td>1.92 ± 0.18</td><td>18.5 ± 4.6</td></tr><tr><td>T1-K47A</td><td>0.70 ± 0.44</td><td>8.30 ± 0.78</td><td>10.4 ± 4.5</td></tr><tr><td colspan="4">DIV</td></tr><tr><td>T1-WT</td><td>0.148 ± 0.001</td><td>0.56 ± 0.020</td><td>37.4 ± 1.7</td></tr><tr><td>T1-F12A</td><td>0.143 ± 0.006</td><td>0.55 ± 0.039</td><td>38.7 ± 3.8</td></tr><tr><td>T1-K47A</td><td>0.427 ± 0.300</td><td>1.14 ± 0.830</td><td>26.5 ± 0.9</td></tr></tbody></table>

<figure><img src="image_14.png" /><figcaption>**Figure 6.** Biochemical analysis of LRP-1-mediated endocytosis of wild-type TIMP-1 and mutants. Cortical neurons from mouse embryos cultured for 48 h on poly-L-lysine-coated coverslips were incubated at 4°C for 1 hour with 5 nM fluoT1-WT, fluoT1-F12A or fluoT1K47A in the presence or absence of 500 nM RAP or 30 µg/mL of blocking LRP-1 polyclonal antibodies (R2629). After careful washing, part of the cells was used to quantify the surface bound and intracellular signal at 4°C. The other part was then incubated at 37°C for an additional 10 min to quantify the surface bound and intracellular signal. Fluorescence intensity was quantified by spectrophotometry and expressed as arbitrary units (A.U.). Values below 300 A.U. are considered to be nonspecific. Error bars indicate mean ± SD, *** indicates significantly different with a P-value < 0.001 and NS indicates a value that is not statistically significant.</figcaption></figure>

in computational methods allow biomolecular simulation to be considered as an alternative to experimental structural procedures. Consequently, *in silico* approaches remain a reliable way to describe and predict protein movements in association with biological experiments.
<figure><img src="image_15.png" /><figcaption>**Figure 7.** Confocal analysis of LRP-1-mediated endocytosis. Cortical neurons from mouse embryos cultured for 24 h on poly-L-lysine-coated coverslips were incubated at 4 °C for 30 min with 5 nM T1-WT, T1-F12A or T1K47A in the presence or absence of 500 nM RAP and then incubated at 37 °C for 10 min. Cells were stained with Alexa Fluor 488 for TIMP-1 and Alexa Fluor 568 for EEA1 and analysed by confocal microscopy. Images were treated with AMIRA software. Maximal Intensity Projection (MIP) in greyscale of TIMP-1 labelling (first panel) and EEA1 labelling (second panel), TIMP-1/EEA1 merge (third panel, TIMP-1 in red and EEA1 in green), differential interference contrast (DIC) and isosurface representation of TIMP-1/EEA1 colocalisation (fourth panel) are shown. Images are representative of more than three distinct experiments. Scale bar: 10 µm.</figcaption></figure>

The characterisation of the TIMP-1/LRP-1 interaction could be physiologically relevant in some neurodegenerative disorders. Molecular docking is the classical method used to predict how a ligand binds a receptor. However, this method necessitates the structural determination of both partners, and the high molecular weight of LRP-1 excludes the use of this method. In this study, we combined NMA and MD simulations with biological experiments to gain insight into the TIMP-1/LRP-1 interaction. We observed large displacement between the N- and C-terminal TIMP-1 domains and identified a hinge region that could be essential to this movement. Residues F12 and K47 seem to play a role in the intrinsic dynamics of TIMP-1. Interestingly, these residues are highly conserved in numerous mammal species (mouse, rabbit, horse...). This suggests the importance of these residues in TIMP-1 functions. *In silico* mutants have been proposed on these residues located in the hinge region, and their biological activity has been tested in mouse cortical neurons. We show that F12A or K47A mutations do not prevent TIMP-1 mutants from binding to LRP-1. Nevertheless, these mutations block their endocytosis and ability to reduce neurite length. Taken together, these results suggest that TIMP-1 binding to LRP-1 is followed by a TIMP-1 conformational change, which is required for its endocytosis and cellular effects.

In the absence of structural data regarding one of the two partners, considering the intrinsic dynamics of the proteins could provide additional information for studying ligand/receptor interactions and their associated cellular effects.

## Methods

**Molecular Modelling.** NMA of the TIMP-1 protein was performed using the PDB 1UEA crystallographic structure. After the addition of hydrogen atoms, the structure was successively minimised using a combination of the Steepest Descent (1000 steps) and ABNR (50000 steps) methods. NMA was carried out using the Bio3D package within the Rstudio environment²⁷. We used the calphax force field to account for stronger beta bridges and helix 1-4 interactions, specifying the six disulfide bridges (1-70, 3-99, 132-137, 145-166, 127-174, 13-124) when building the force constants. Graphs were generated with Rstudio.

A 100-ns-long molecular dynamics simulation in explicit solvent (TIP3P) was performed on the TIMP-1 minimised structure used for NMA with the NAMD software, along with the CHARMM27 force field, in the
<figure><img src="image_16.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>**Figure 8.** Quantification of TIMP-1 and EEA1 labelling. (a) Quantification of average intensity of TIMP-1 labelling per µm² for each cell (n > 10) using ImageJ software. The data are representative of at least three independent experiments with at least five mice embryos. Error bars indicate mean ± SD and **** indicates significantly different from T1-WT without RAP and with a P-value < 0.0001. (b) Quantification of average intensity of EEA1 labelling per µm² for each cell (n > 10) using ImageJ software. Error bars indicate mean ± SD and NS indicates a value that is not statistically significant. (c) Quantification of TIMP-1/EEA1 colocalisation (n > 10) calculated using AMIRA software. Error bars indicate mean ± SD and *** indicates significantly different from T1-WT without RAP with a P-value < 0.001.</figcaption></figure>

NPT ensemble at 293 K - 1 atm, with a cut-off of 11 Å for long-range interactions. Trajectory analysis and graphs were produced using the Bio3D package within the Rstudio environment²⁸. Pictures were generated with PyMol software.
**Plasmid construction, production and purification of recombinant proteins.** The pCMV6-XL5 vector containing TIMP-1 cDNA was obtained from Origene (Origene, Rockville, USA). Mutated TIMP-1 forms were generated at the F12 and K47 positions using the QuickChange II Site-directed Mutagenesis kit (Stratagene, Agilent Technologies, Les Ulis, France) with the following primers: F12A (forward): 5'-CCACCACAGACGGCCGCCTGCAATTCCGACCTCG-3'; (reverse): 5'-CGAGGT CGGAATTGCAGGCGGCCGTCTGTGGGTGG-3', K47A (forward): 5'-GATGACCAAGATGTATGCAGGGT TCCAAGCCTTAGGG-3'; (reverse): 5'-CCCTAAGGCTTGGAACCCTGCATACATCTTGTCATC-3'. TIMP-1 and mutated cDNA were cloned into the p3XFLAG-CMV14 vector (Sigma-Aldrich, Saint Quentin Fallavier, France) to produce TIMP-1 with fusion FLAG-tag in C-terminal using the following primers: 5'-TACAGAATTCCACCATGGCCCC-3' (forward) and 5'-GAGAGCTAGCGGCTATCTGGGACC-3' (reverse). All constructs were confirmed by DNA sequencing analysis (Beckman Coulter Genomics, France). CHO cells were stably transfected with these constructs using Lipofectamine 2000 (Invitrogen, Cergy Pontoise, France) to produce FLAG-tagged TIMP-1 (T1-WT), FLAG-tagged TIMP-1 F12A (T1-F12A) and FLAG-tagged TIMP-1 K47A (T1-K47A). Cells were grown in DMEM-Ham's F-12 medium with 10% (v/v) foetal bovine serum and 500 ng/ml geneticin at 37°C in a humidified chamber containing 5% (v/v) CO₂. The conditioned media were collected and incubated for 1 h with 10 mM EDTA in acidic conditions (pH 3.5) to remove the MMPs complexed to wild-type or mutated TIMP-1. After equilibration to pH 7.4 with NaOH, the conditioned media were incubated at 4°C for 3 h with anti-FLAG M2 affinity gel (Sigma-Aldrich). The resin was then collected and washed 5 times with Tris-buffered saline (TBS, 50 mM Tris-HCl pH 7.4 containing 150 mM NaCl). The recombinant proteins were eluted under acidic conditions (0.1 M glycine-HCl pH 3.5), equilibrated at pH 7.4 by 1 M Tris pH 8.0 and dialysed against PBS. TIMP-1 wild-type and mutant concentrations were determined with the Human TIMP-1 Quantikine ELISA kit (R&D Systems, Lille, France). Recombinant protein purity and quantification were controlled by SDS-PAGE followed by silver nitrate staining and immunoblotting. Recombinant RAP was prepared as previously described⁶.

**Western blotting.** Conditioned media and whole-cell protein extracts were prepared and analysed by western blotting using anti-TIMP-1 (1/500, Merck Biosciences), anti-FLAG (1/1000, mouse, clone M2, Sigma-Aldrich) or anti-actin (1/1000, Santa-Cruz Biotechnology), as previously described⁶.

**MMP inhibition assays.** The inhibitory activities of T1-WT, T1-F12A and T1-K47A against MMP-1, -2, -3 and -9 were assayed as follows. All MMPs and substrates were obtained from Merck Biosciences.

**MMP-1:** proMMP-1 was activated by incubation in activation buffer (25 mM Hepes, 5 mM CaCl₂, 20% (v/v) glycerol, 0.01% (v/v) Brij 35 and 4.5 mM APMA) for 2 h at 37°C. More than 80% of the proMMP-1 was activated in these conditions (data not shown). The activated MMP-1 (2 nM) was then incubated for 1 h at 27°C in Tris-test buffer (50 mM Tris-HCl pH 7.5 containing 150 mM NaCl and 5 mM CaCl₂) containing 0.92 to 6.44 nM TIMP-1 or its mutated forms. The assay was initiated by adding MMP-1/-9 fluorogenic substrate DNP-Pro-Cha-Gly-Cys(Me)-His-Ala-Lys(N-Me-Abz)-NH₂ ($\lambda_{exc}$ = 365 nm/$\lambda_{em}$ = 450 nm) at a final concentration of 5 µM.

**MMP-2:** 0.76 nM MMP-2 was incubated in Tris-test buffer containing 0.17 to 13.5 nM TIMP-1 or its mutated forms for 2 h at 27°C before adding MMP-2/-7 fluorogenic substrate MCA-Pro-Leu-Gly-Leu-Dpa-Ala-Arg-NH TFA ($\lambda_{exc}$ = 325 nm/$\lambda_{em}$ = 393 nm) at a final concentration of 2.5 µM.

**MMP-3:** 2.3 nM MMP-3 catalytic domain was incubated with 0.1 to 5.4 nM TIMP-1 or its mutated forms in Tris-test buffer for 1 h at 27°C. MMP-3 fluorogenic substrate NBD-Arg-Pro-Lys-Pro-Leu-Ala-Nva-Trp-Lys(DMC)-NH₂ ($\lambda_{exc}$ = 350 nm/$\lambda_{em}$ = 465 nm) was added at a final concentration of 3 µM.

**MMP-9:** 0.3 nM MMP-9 was incubated in Tris-test buffer containing 0.7 to 11 nM TIMP-1 or its mutated forms for 2 h at 27°C. The assay was initiated by adding MMP-1/-9 fluorogenic substrate (DNP-Pro-Cha-Gly-Cys(Me)-His-Ala-Lys(N-Me-Abz)-NH₂ ($\lambda_{exc}$ = 365 nm/$\lambda_{em}$ = 450 nm) at a final concentration of 5 µM.

The rate of each substrate cleavage was measured in triplicate for each concentration examined, using an Infinite F200 PRO spectrofluorimeter (Tecan, Lyon, France) with one measure per minute for 20 min. Non-linear regression analysis with Graphpad software (La Jolla, USA) allowed us to calculate the Kᵢ using the Morrison equation⁴³.

**Surface plasmon resonance analysis.** Surface plasmon resonance experiments were performed using Biacore T200 (GE Healthcare Life Sciences, Velizy-Villacoublay, France) as previously described⁶. Injections were performed with HBS-N buffer (50 mM Hepes, 150 mM NaCl, Tween 0.005% (v/v), pH 7.4). Recombinant LRP-1 mini-domains DII and DIV (R&D) with an Fc-tag were immobilised on a CMD500m sensor chip (Xantech) functionalised with anti-Fc antibody. Different concentrations of TIMP-1 and of its mutated forms (5–160 nM) were injected in multi-cycle kinetics to allow determination of the association ($k_{on}$) and dissociation ($k_{off}$) constants. EGF was used as the negative control.

**Neuron morphological analysis.** Primary cultures of cortical neurons were prepared from CD1 mice embryos as previously described⁶ and cultured as previously described³⁰. These experiments involving mice were carried out in accordance with the French and European guidelines for the care and use of laboratory animals and were approved by the ethic committee of the University of Reims Champagne-Ardenne (Protocol number: 1A08248094720). According to the specific biological activity of TIMP-1 in neurons⁶,³⁰, we chose to use a concentration of 5 nM of TIMP-1, which allowed TIMP-1 quantification in the endocytosis assay. Five nM of T1-WT, T1-F12A or T1-K47A or mixtures of T1-WT and mutants with different ratios were added to the cultures in serum-free media for 30 min. The neurons were rinsed with PBS and fixed in 4% (v/v) paraformaldehyde for 10 min. Immunofluorescence studies with anti-βIII-tubulin antibodies (1/100, clone SDL.3D10, Sigma-Aldrich) were performed as previously described³⁰. The cells were observed using an Olympus BH2-RFCA fluorescent
microscope with a 40x objective lens and photographed with a U-MPTV Camera using the DP Controller software (all from Olympus). In each experiment, fluorescence microphotographs were taken from at least 5 randomly selected fields per slide, with 4 slides per experimental condition (total 20 fields). The neurite length corresponded to the average length of the neurites of each cell and was obtained by dividing the total length of the neurites in a field by the number of neurite segments, as previously described⁶. The neurite length of all βIII-tubulin positive cells was measured using the ImageJ plugin NeuronJ. A neurite segment was defined as the distance between branching points or the distance between the branching point and the top of a neurite.

**Fluorescent TIMP-1 preparation and endocytosis experiment.** Ten µg of T1-WT, T1-F12A or T1-K47A were biotinylated with the ProtOn Biotin Labelling Kit according to the manufacturer's instructions. The biotinylated products were incubated with 0.5 µg eFluor® 450 Streptavidin (eBioscience) for 30 min at RT. The neuron-surface binding and endocytosis of 5 nM fluoT1-WT, fluoT1-F12A or fluoT1-K47A were studied as previously described⁶. Briefly, cortical neurons plated in 24-well plates were incubated in assay medium (culture medium containing 0.1% (v/v) bovine serum albumin, BSA) at 4 °C for 1 h. For the binding assays, cells were incubated with 5 nM fluo-T1 (-WT, -F12A and -K47A) at 4 °C for 2 h in the absence or presence of 500 nM RAP or 30 µg/mL blocking LRP-1 polyclonal antibodies (R2629)³⁴. The cells were then carefully rinsed five times with cold PBS and surface-digested with 0.1% (w/v) pronase in culture medium at 4 °C to degrade the surface-bound ligands and to detach cells. After cell collection by centrifugation, fluorescence in the supernatant corresponding to the surface-bound ligand and fluorescence in the pelleted cells corresponding to the internalised ligand were measured by spectrofluorimetry (λexc: 405 nm/λem: 445 nm) (Infinite F200 PRO, Tecan, Lyon, France). Appropriate controls using only eFluor were performed. For the endocytosis assays, after fluo-T1 binding (as described above), the cells were carefully rinsed with PBS and further cultured in assay medium pre-warmed at 37 °C for 15 min. To distinguish surface binding from intracellular accumulation, the cells were washed twice with cold PBS and surface-digested with pronase as described above.

**Confocal microscopy analysis of LRP-1 mediated endocytosis.** Neurons were first equilibrated at 4 °C for 30 min with cold serum-free defined medium supplemented with 1 mg/ml BSA to saturate the cell surface. Five nM of T1-WT, T1-F12A or T1-K47A were then added for 60 min at 4 °C in the presence or absence of 500 nM RAP. Endocytosis was allowed by incubation at 37 °C for 10 min or 1 hour and stopped by replacing the cells at 4 °C. The neurons were extensively rinsed with cold PBS and fixed in 4% (v/v) paraformaldehyde for 10 min. For immunofluorescence, cells were incubated with anti-TIMP-1 (1/300, Merck Biosciences) and anti-EEA1 (1/300, Abcam) at 4 °C overnight, followed by secondary antibodies conjugated to Alexa Fluor 488 (1/1000, green, Molecular Probes) or to Alexa Fluor 568 (1/1000, red, Molecular Probes) at room temperature for 60 min. After rinsing with PBS, the slides were mounted in fluorescence mounting medium. Acquisitions were made with a confocal microscope LSM710 with a 63x objective and the Zeiss operating system associated with the ZEN software program (all from Carl Zeiss MicroImaging GmbH, Germany). Acquisitions were performed by exciting Alexa Fluor 488 and Alexa Fluor 568 with an argon laser, a helium-neon laser and a chameleon infrared laser tuned at 730 nm. Emitted fluorescence was detected through the appropriate wavelength window. In each experiment, fluorescence microphotographs were taken from at least 10 randomly selected fields per slide, with 3 slides per experimental condition (total 30 fields). In each field, 20 images were captured with a 0.3 µm z-step. Maximal intensity projection representations were realised using the AMIRA software program (v5.4.2, Visualization Sciences Group, Burlington, MA, USA). For the colocalisation studies, a multichannel field module was used, followed by a correlation plot treatment (subrange values of 15 to 255; gamma, 0.5). Then, an isosurface representation of the correlation plot was realised using the same threshold for each treatment. The TIMP-1 average intensity of fluorescence per µm² for each cell was performed using the particle analysis function of the ImageJ software. Quantification of the TIMP-1/EEA1 colocalisation was performed using the AMIRA software program.

**Statistical analysis.** Statistical analysis was performed using a Student's t-test or one-way ANOVA with a Newman Keuls post-test (GraphPad Prism version 6.01 for Windows, GraphPad Software, San Diego California USA, www.graphpad.com). The values represent the mean ± SD of at least three independent experiments. P-values of <0.05, <0.01, <0.001 and <0.0001 are identified with *, **, *** and ****, respectively. NS represents a value that is not statistically significant.

**References**

1. Murphy, G. Tissue inhibitors of metalloproteinases. *Genome Biology* **12**, 233 (2011).

2. Arpino, V., Brock, M. & Gill, S. E. The role of TIMPs in regulation of extracellular matrix proteolysis. *Matrix Biology: Journal of the International Society for Matrix Biology* **44–46**, 247–254 (2015).

3. Ries, C. Cytokine functions of TIMP-1. *Cellular and Molecular Life Science: CMLS* **71**, 659–672 (2014).

4. Stetler-Stevenson, W. G. Tissue inhibitors of metalloproteinases in cell signaling: metalloproteinase-independent biological activities. *Science Signaling* **1**, re6 (2008).

5. Bridoux, L. et al. A crucial role for Lyn in TIMP-1 erythroid cell survival signalling pathway. *FEBS Letters* **587**, 1524–1528 (2013).

6. Thevenard, J. et al. Low-density lipoprotein receptor-related protein-1 mediates endocytic clearance of tissue inhibitor of metalloproteinases-1 and promotes its cytokine-like activities. *PloS One* **9**, e103839 (2014).

7. Lillis, A. P., Van Duyn, L. B., Murphy-Ullrich, J. E. & Strickland, D. K. LDL receptor-related protein 1: unique tissue-specific functions revealed by selective gene knockout studies. *Physiological Reviews* **88**, 887–918 (2008).

8. Bu, G., Maksymovitch, E. A., Geuze, H. & Schwartz, A. L. Subcellular localization and endocytic function of low density lipoprotein receptor-related protein in human glioblastoma cells. *The Journal of Biological Chemistry* **269**, 298742–9882 (1994).

9. Johnson, K. M. & Crocker, S. J. TIMP-1 couples RhoK activation to IL-1beta-induced astrocyte responses. *Neuroscience Letters* **609**, 165–170 (2015).
10. Mroczko, B., Groblewska, M. & Barcikowska, M. The role of matrix metalloproteinases and tissue inhibitors of metalloproteinases in the pathophysiology of neurodegeneration: a literature study. *Journal of Alzheimer's disease: JAD* **37**, 273–283 (2013).

11. Guedes, I. A., de Magalhães, C. S. & Dardenne, L. E. Receptor–ligand molecular docking. *Biophysical Reviews* **6**, 75–87 (2014).

12. Du, X. *et al.* Insights into protein-ligand interactions: mechanisms, models, and methods. *International Journal of Molecular Sciences* **17** (2016).

13. Gomis-Ruth, F. X. *et al.* Mechanism of inhibition of the human matrix metalloproteinase stromelysin-1 by TIMP-1. *Nature* **389**, 77–81 (1997).

14. Brew, K. & Nagase, H. The tissue inhibitors of metalloproteinases (TIMPs): an ancient family with structural and functional diversity. *Biochimica et Biophysica Acta* **1803**, 55–71 (2010).

15. Murphy, G. *et al.* The N-terminal domain of tissue inhibitor of metalloproteinases retains metalloproteinase inhibitory activity. *Biochemistry* **30**, 8097–8102 (1991).

16. Etique, N., Verzeaux, L., Dedieu, S. & Emonard, H. LRP-1: a checkpoint for the extracellular matrix proteolysis. *BioMed Research International* **2013**, 152163 (2013).

17. Jiang, X., Dias, J. A. & He, X. Structural biology of glycoprotein hormones and their receptors: insights to signaling. *Molecular and Cellular Endocrinology* **382**, 424–451 (2014).

18. Mittag, T., Kay, L. E. & Forman-Kay, J. D. Protein dynamics and conformational disorder in molecular recognition. *Journal of Molecular Recognition: JMR* **23**, 105–116 (2010).

19. Skjaerven, L., Reuter, N. & Martinez, A. Dynamics, flexibility and ligand-induced conformational changes in biological macromolecules: a computational approach. *Future Medicinal Chemistry* **3**, 2079–2100 (2011).

20. Tama, F. Normal mode analysis with simplified models to investigate the global dynamics of biological systems. *Protein and Peptide Letters* **10**, 119–132 (2003).

21. Ma, J. Usefulness and limitations of normal mode analysis in modeling dynamics of biomolecular complexes. *Structure* **13**, 373–380 (2005).

22. Skjaerven, L., Martinez, A. & Reuter, N. Principal component and normal mode analysis of proteins: a quantitative comparison using the GroEL subunit. *Proteins* **79**, 232–243 (2011).

23. Floquet, N., Dedieu, S., Martiny, L., Dauchez, M. & Perahia, D. Human thrombospondin's (TSP-1) C-terminal domain opens to interact with the CD-47 receptor: a molecular modeling study. *Archives of Biochemistry and Biophysics* **478**, 103–109 (2008).

24. Jeanne, A., Schneider, C., Martiny, L. & Dedieu, S. Original insights on thrombospondin-1-related antireceptor strategies in cancer. *Frontiers in Pharmacology* **6**, 252 (2015).

25. Grossman, M. *et al.* The intrinsic protein flexibility of endogenous protease inhibitor TIMP-1 controls its binding interface and affects its function. *Biochemistry* **49**, 6184–6192 (2010).

26. Hollup, S. M., Salensminde, G. & Reuter, N. WEBnm@: a web application for normal mode analyses of proteins. *BMC Bioinformatics* **6**, 52 (2005).

27. Hinsen, K. Analysis of domain motions by approximate normal mode calculations. *Proteins* **33**, 417–429 (1998).

28. Grant, B. J., Rodrigues, A. P., ElSawy, K. M., McCammon, J. A. & Caves, L. S. Bio3d: an R package for the comparative analysis of protein structures. *Bioinformatics* **22**, 2695–2696 (2006).

29. Hamze, A. B., Wei, S., Bahudhanapati, H., Kota, S., Acharya, K.R., Brew K. Constraining specificity in the N-domain of tissue inhibitor of metalloproteinases-1; gelatinase-selective inhibitors. *Protein science: a publication of the Protein Society* **16**, 1905–1913 (2007).

30. Ould-yahoui, A. *et al.* A new role for TIMP-1 in modulating neurite outgrowth and morphology of cortical neurons. *PloS One* **4**, e8289 (2009).

31. Van Gool, B., Dedieu, S., Emonard, H. & Roebroek, A. J. The matricellular receptor LRP1 forms an interface for signaling and endocytosis in modulation of the extracellular tumor environment. *Frontiers in Pharmacology* **6**, 271 (2015).

32. Obermoeller-McCormick, L. M. *et al.* Dissection of receptor folding and ligand-binding property with functional minireceptors of LDL receptor-related protein. *Journal of Cell Science* **114**, 899–908 (2001).

33. Herz, J., Goldstein, J. L., Strickland, D. K., Ho, Y. K. & Brown, M. S. 39-kDa protein modulates binding of ligands to low density lipoprotein receptor-related protein/alpha 2-macroglobulin receptor. *The Journal of Biological Chemistry* **266**, 21232–21238 (1991).

34. Mikhailenko, I. *et al.* Recognition of alpha 2-macroglobulin by the low density lipoprotein receptor-related protein requires the cooperation of two ligand binding cluster regions. *The Journal of Biological Chemistry* **276**, 39484–39491 (2001).

35. Laatsch, A. *et al.* Low density lipoprotein receptor-related protein 1 dependent endosomal trapping and recycling of apolipoprotein E. *PloS One* **7**, e29385 (2012).

36. Mu, F. T. *et al.* EEA1, an early endosome-associated protein. EEA1 is a conserved alpha-helical peripheral membrane protein flanked by cysteine "fingers" and contains a calmodulin-binding IQ motif. *The Journal of Biological Chemistry* **270**, 13503–13511 (1995).

37. Li, Y., Marzolo, M. P., van Kerkhof, P., Strous, G. J. & Bu, G. The YXXL motif, but not the two NPXY motifs, serves as the dominant endocytosis signal for low density lipoprotein receptor-related protein. *The Journal of Biological Chemistry* **275**, 17187–17194 (2000).

38. Ranganathan, S. *et al.* Serine and threonine phosphorylation of the low density lipoprotein receptor-related protein by protein kinase C alpha regulates endocytosis and association with adaptor molecules. *The Journal of Biological Chemistry* **279**, 40536–40544 (2004).

39. Boucher, P. & Gotthardt, M. LRP and PDGF signaling: a pathway to atherosclerosis. *Trends in Cardiovascular Medicine* **14**, 55–60 (2004).

40. Keskin, O. Binding induced conformational changes of proteins correlate with their intrinsic fluctuations: a case study of antibodies. *BMC Structural Biology* **7**, 31 (2007).

41. Plattner, N. & Noe, F. Protein conformational plasticity and complex ligand-binding kinetics explored by atomistic simulations and Markov models. *Nature Communications* **6**, 7653 (2015).

42. Goh, C. S., Milburn, D. & Gerstein, M. Conformational changes associated with protein-protein interactions. *Current Opinion in Structural Biology* **14**, 104–109 (2004).

43. Morrison, J. F. Kinetics of the reversible inhibition of enzyme-catalysed reactions by tight-binding inhibitors. *Biochimica et Biophysica Acta* **185**, 269–286 (1969).

## Acknowledgements

The authors thank Dr. Nicolas Floquet for preliminary studies in this project, Olivier Boquet for his technical assistance in generating mutants, Jean-Luc Breda for animal care, Christine Terryn for technical support with confocal imaging and Sebastien Almagro for the ImageJ software plugins. They are grateful to Dr. D.K. Strickland (Department of Surgery, University of Maryland School of Medicine, Baltimore, MD, USA) for providing R2629 antibodies. They thank the HPC Regional Center ROMEO and the PICT (Cell and Tissue Imaging Platform) from the University of Reims Champagne-Ardenne (France) for providing time and support. This work was supported by the Université de Reims Champagne-Ardenne, the Region Champagne-Ardenne and the Centre National de la Recherche Scientifique (CNRS), France. Laurie Verzeaux was the recipient of a grant from the Region Champagne-Ardenne.
## Author Contributions

N.E. and E.D.-C. contributed to the conception and design of the study; N.B. performed the molecular modelling experiments. L.V., J.T.-D., J.D., G.F. and N.E. performed biological experiments. L.V., M.D., H.E., S.D., E.D.-C. and N.E. analysed the data and performed statistical analyses. H.E., E.D.-C. and N.E. drafted the manuscript and L.M. and S.D. revised the manuscript. All authors approved the final version of the manuscript.

## Additional Information

**Supplementary information** accompanies this paper at doi:10.1038/s41598-017-05039-z

**Competing Interests:** The authors declare that they have no competing interests.

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_19.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2017