1. Introduction
The 491-residue-long human version of MDM2 (mouse double minute 2, EC. 2.3.2.27) protein [
1,
2] is a negative regulator of p53, which is a tumor suppressor involved in apoptosis, DNA repair, senescence, and activation of genes implicated in cell cycle arrest [
3,
4]. MDM2 can interact with p53 through different ways: (i) the N-terminal region of MDM2 (N-MDM2), comprising the first 125 residues, interacts with the p53 N-terminal transactivation (TA) domain (residues 15–29), altering p53-transcriptional activity; (ii) the full-length MDM2 exports p53 from the nucleus with the help of an export signal in its sequence [
2,
5,
6]; and (iii) the C-terminal region of MDM2 has a ubiquitin ligase function towards p53, activating p53 degradation via ubiquitin–proteasome pathways [
3,
7,
8]. Thus, the MDM2/p53 interaction hampers p53 transcriptional regulatory activity and accelerates degradation of p53 within the proteasome [
7]. On the other hand, there are various mechanisms that allow p53 to escape the interaction with MDM2 [
9]. For instance, the protein ARF, whose expression is induced by the oncogenic proteins cMyc and Ras, is capable of binding to MDM2, abolishing its ability to bind p53 and hampering its p53-regulatory properties [
10].
From a structural point of view, MDM2 is an intrinsically disordered protein (IDP) which possesses long unfolded regions connecting three relatively small and well-folded domains (
Figure 1). N-MDM2 includes the largest of them, which is the p53 binding domain (residues 25–109 [
11]), and is composed by four α-helices and a smaller fraction of secondary structures organized in β-strands [
6,
12] (inset in
Figure 1). The other two domains are a zinc finger (residues 300–330) and a C-terminal RING domain (starting from residue 438 onwards).
Studies aimed to understand the molecular mechanism behind the interplay between p53 and N-MDM2 have shown that the phosphorylation of residues in p53 is important in modulating such interaction and that the binding of different p53-derived ligands to N-MDM2 induces long-range conformational changes in the latter [
13,
14,
15]. Furthermore, interactions with p53 occur through the movement of a so-called “lid”, which shifts in the presence of p53. Therefore, it seems that N-MDM2, in spite of being largely endowed with a structure, has an inherent and high conformational plasticity which could allow this protein to bind several different partners. In fact, N-MDM2 not only binds the TA of p53 but, among others, also the S100 proteins [
16], the SAM domain of p73 [
17], and the citrullinating enzyme PADI4 [
18].
In this work, we studied the stability and conformational features of isolated N-MDM2 under different acidic pH conditions using several biophysical techniques, namely fluorescence, circular dichroism (CD), size exclusion chromatography (SEC), and differential scanning calorimetry (DSC). Furthermore, we modeled the simulated unfolding of N-MDM2 by using an analysis of the protein constraint network based on graph theory [
19]. All the findings obtained with the different experimental techniques indicate that N-MDM2 has a low conformational stability and that it acquires a native-like conformation in a narrow pH range (between pH ~7.0 and 10.0). Despite its small size, thermal or chemical unfolding of N-MDM2 was found to be a complex process, showing the presence of several intermediates with different amounts of secondary and tertiary structures. The presence of these intermediates was also observed by the in silico experiments, which captured the occurrence of multiple sub-structures organized in a hierarchical fashion during the simulated denaturation process.
3. Discussion
It has been shown, based on NMR measurements, that N-MDM2 is a functional domain that is very sensitive to conformational rearrangements, and even at 15 °C, an evident plasticity of this protein has been observed [
12] in some regions, which undergoes conformational re-adjustments to allow for the allocation of p53 [
11]. These observations seem to suggest a high inner mobility of the domain and a low stability. We have tried to address these questions in our study by investigating the structural features and the conformational stability of the domain from different in vitro and in silico approaches.
The first conclusion of our work is that N-MDM2 in solution had a rather low stability (with an apparent thermal denaturation midpoint of ~55 °C) in a narrow pH interval (between pH 7.0 and 10.0). A similar low stability has been reported in other prevalently α-helical proteins [
31,
32,
33,
34], where it has been associated with a large flexibility. It is interesting to compare our results with those obtained with model proteins such as chicken and equine lysozyme [
35,
36]. In fact, when the specific enthalpy of denaturation in the two transitions is combined and normalized with respect to the molecular mass MW (i.e., the quantity (ΔH
m1 + ΔH
m2)/MW is considered), the value obtained for N-MDM2 (9.4 cal g
–1) may seem comparable to that of chicken lysozyme (7.6 cal g
–1) [
35], which is considered a standard for the packing density of links in the native state. Nevertheless, the stability of N-MDM2 is essentially dictated by the first transition because it represents the loss of structural integrity and, very likely, will affect some of its interactions with other biomolecules. Such transition is characterized by an unfolding temperature of 52 °C and an unfolding enthalpy of 55 kcal mol
–1, and these parameters are much lower than those for chicken lysozyme, which shows a single transition characterized by an unfolding temperature of 75 °C and an unfolding enthalpy of about 100 kcal mol
–1. Thus, the unfolding features of N-MDM2 and chicken lysozyme are quite different: the unfolding of N-MDM2 occurs along two transitions at moderate-to-low temperatures, whereas in the case of chicken lysozyme, the unfolding takes place as a single transition at high temperatures. Furthermore, the stabilization energy is much lower for N-MDM2: at 52 °C, N-MDM2 has already lost a significant part of its structure whereas, at the same temperature, the structure of chicken lysozyme is little affected. In contrast to chicken lysozyme, when compared to N-MDM2, equine lysozyme shows a similar stability profile [
36], including two transitions with low unfolding temperatures and unfolding enthalpies. We suggest that the low stability of N-MDM2 would ease the interactions with other macromolecules and it is an essential feature of this domain to carry out its function(s).
The second conclusion of our work is that the loss of the secondary structure and the burial of solvent-exposed hydrophobic patches (as monitored by far-UV CD and ANS fluorescence), as the pH was varied, occurred in an identical manner (
Figure 2). On the other hand, the acquisition of tertiary structures (as monitored by intrinsic fluorescence,
Figure 2A) and compactness of the polypeptide chain (as shown by SEC,
Figure 3A) occurred at a different step, although in the latter technique we cannot rule out the possibility of protein–column interactions of the partially folded structures present at acidic pH values. At acidic pH, because of the similarity of the p
Ka values with those found in model peptides [
37,
38,
39], the loss of native-like structures could be associated with the presence of aspartic or glutamic acid residues. The fact that the populated species at low pH showed a large ANS fluorescence intensity suggests that they are probably molten globules [
20], as we could further confirm by the absence of thermal unfolding sigmoidal curves followed by fluorescence or CD at low pH values (
Figures S2 and S4). The narrow pH range (from 7.0 to 10.0) where N-MDM2 acquired a fully native-like conformation ensures that the domain is only functional under physiological conditions.
Furthermore, we also explored the unfolding of N-MDM2. As a preliminary remark, we cannot exclude at this stage that folding of the domain, and even its stability, might be partly modified due to its isolation from the whole MDM2. In particular, we cannot rule out that interactions with other domains of the intact protein could modulate its conformational and folding features. The chemical denaturation curves with different defined steps, depending on the probe used (intrinsic fluorescence, ANS fluorescence, SEC, or far-UV CD) and the broadening of some of the transitions (e.g., that from fluorescence with large
m-values,
Figure 5) that demonstrate a low cooperativity, all contribute to indicating that the unfolding of N-MDM2 was not a two-state process [
40]. These results were further supported by our SEC findings (
Figure 3B).
The simulation results confirmed the structural plasticity of N-MDM2, which was especially visible in the complex nature of its unfolding process. During this reaction, the increase in flexibility of the constraint network in the protein was visible in all the global indexes analyzed (
Figure 6). Moreover, the configuration entropy Σ clearly suggests a hierarchy in the melting of the different protein sub-units, which could be expected for larger proteins with a more complex architecture [
21] but is more unusual for a relatively small protein such as N-MDM2. The phase transition observed for Σ corresponds to the first and main event in the denaturation process of the protein, in which the melting of the structure leads to distinct and independent sub-units. In contrast, the rigidity order parameter Π appeared to be more sensitive to later events along the unfolding pathway of N-MDM2. Local indexes in the simulated unfolding process, obtained as a function of the residue number (
Figure 7) and mapped onto the native structure of N-MDM2 (
Figure 8), provided specific details on the different rigidity of the various regions of this domain. The variety observed in the inner flexibility of N-MDM2 determines the presence of multiple intermediates along the unfolding process, which could also be predicted by other less-sensitive computational tools (e.g., the folding kinetic prediction tool K-Fold [
41]) and, more importantly, whose presence was directly captured in our wet-laboratory experiments. The collective behavior of the sub-structures of the N-MDM2 domain, visible in the denaturing process, could be expected to be present, at least in part, also at physiological temperatures and to contribute to dictating its overall functional properties.
4. Materials and Methods
4.1. Materials
Imidazole, Trizma base, DNase, SIGMAFAST protease tablets, NaCl, Ni2+-resin, and Amicon centrifugal devices, with a molecular weight cut-off of 3 kDa, were from Sigma (Madrid, Spain). The β-mercaptoethanol was from BioRad (Madrid, Spain). Ampicillin and isopropyl-β-D-1-thiogalactopyranoside were obtained from Apollo Scientific (Stockport, UK). Triton X-100, Tris(2-carboxyethyl)phosphine (TCEP), dialysis tubing with a molecular weight cut-off of 3500 Da, and the SDS protein marker (PAGEmark Tricolor) were from VWR (Barcelona, Spain). The rest of the used materials were of analytical grade. Water was deionized and purified on a Millipore system.
4.2. Protein Expression and Purification
Monomeric His-tagged N-MDM2 (residues 6–125 of the intact protein) was purified as previously described [
17]. The first 5 residues of N-MDM2 were missing due to optimization of the cloning sequence. Protein concentration was determined by UV absorbance, employing an extinction coefficient at 280 nm estimated from the number of tyrosines (the construct of His-tagged N-MDM2 has eight tyrosines) in the sequence [
42].
4.3. Fluorescence
Fluorescence spectra were collected on a Cary Varian spectrofluorometer (Agilent, Santa Clara, CA, USA) interfaced with a Peltier unit. Following the standard protocols used in our laboratories, the samples were prepared the day before and left overnight at 5 °C; before experiments, samples were left for 1 h at 25 °C. Slit widths for the excitation and emission lights were 5 nm in all cases. A 1 cm-pathlength quartz cell (Hellma, Kruibeke, Belgium) was used. Concentration of N-MDM2 in the samples was 20 µM.
4.3.1. Intrinsic Fluorescence
Protein samples were excited at 280 nm in the pH range 2.0‒12.0. The other experimental parameters have been described elsewhere [
18]. Appropriate blank corrections were made in all spectra. Fluorescence experiments (pH or chemical denaturations) were repeated in triplicate with newly prepared samples. Variations of results among the experiments were lower than 10%.
The pH of each sample was measured after completion of pH denaturations with an ultra-thin Aldrich electrode in a Radiometer (Copenhagen) pH meter. The salts and acids used in the buffers have been described elsewhere [
34]. The wavelength-averaged emission intensity (also called the spectrum mass center), <λ>, was calculated as described [
34]. Briefly, the wavelength-averaged emission intensity, <λ>, is
, where
Ii is the intensity at wavelength λ
i. From its definition, this parameter is an integral of the value of the fluorescence spectrum, and thus, it allows one to obtain overall information on the intensities acquired in the spectrum (instead of using a single wavelength to monitor the fluorescence intensity). It is important to indicate that, conversely to the measurements of a single intensity value at a particular wavelength, <λ> is an intensive variable [
43]. We have reported the value of <λ> in units of µm
−1.
Urea denaturations at pH 7.0 (50 mM, sodium phosphate buffer), either followed by fluorescence or CD, were carried out by dilution of the proper amount of an 8 M urea stock solution. Fluorescence and CD experiments were also acquired in GdmCl (from a stock solution at concentration of 7 M). The denaturant concentration of the stock solutions was determined as described by using the refraction index of the solution [
34]. All the urea or GdmCl denaturations were shown to be irreversible by following <λ>, and therefore we cannot estimate the stability in terms of the Gibbs energy, ΔG, in the denaturation of N-MDM2. This parameter, as other thermodynamic features, should have been obtained from the curves shown in
Figure 5, where the variation in a spectroscopic property (raw ellipticity at a particular wavelength in CD technique; or fluorescence intensity at a particular wavelength or <λ> in fluorescence technique) is represented vs. the concentration of denaturant used.
4.3.2. Thermal Denaturations
Experiments were performed at constant heating rates of 60 °C/h with an average acquisition time of 1 s and collecting points every 0.2 °C. The “average time” is the “sampling time” of the instrument at each temperature. Ideally, in a thermal scan experiment, this time should be much lower than the scan rate of the experiment to ensure that the temperature is constant during the acquisition of fluorescence emissions at a particular temperature (in our conditions, the temperature change during the sampling time was 0.017 °C, and it could be considered negligible). Thermal scans were collected at 308 nm by excitation at 280 nm from 25 to 85 °C. The rest of the experimental set was the same as described above. Thermal denaturations were not reversible at any pH value.
4.3.3. ANS Binding
The excitation wavelength was 380 nm, and emission was measured from 400 to 600 nm. Slit widths were 5 nm for excitation and emission lights. ANS stock solutions were prepared in water and diluted to yield a final concentration of 100 µM. Blank solutions were subtracted from the corresponding spectra. Protein concentration, as in the intrinsic fluorescence experiments, was 20 µM for pH and GdmCl denaturations.
4.3.4. Thermal Unfolding by Using an External Probe (SYPRO Orange)
The thermal stability of N-MDM2 was assessed by differential scanning fluorimetry in a Mx3005P real-time qPCR (Agilent, Madrid, Spain). The fluorescence emission intensity of the extrinsic fluorophore SYPRO Orange (Thermo Fisher Scientific, Barcelona, Spain) was monitored as a function of the temperature at a scanning rate of 1 °C/min, using excitation and emission filters centered at 496 and 610 nm, respectively (the two available filters that were closest to the theoretical values estimated for the fluorophore excitation and emission, 491 and 586 nm). Assays were carried out at protein concentrations of 14 and 25 μM. Upon protein temperature unfolding, SYPRO Orange binds to solvent-exposed hydrophobic regions in the protein, undergoing a fluorescence quantum yield enhancement with an increase in the fluorescence intensity. The unfolding traces were analyzed according to a two-state (single-transition) unfolding model, from which the unfolding temperature and unfolding enthalpy could be estimated. Very briefly, the fluorescence intensity measured at any temperature, <
F>, is given by
where
PN and
PU are the populations of the native (N) and unfolded (U) states, and
FN and
FU are the intrinsic characteristic fluorescence signals for the native and unfolded states that are linear functions of the temperature:
The following standard relationships are considered to calculate <
F>:
where
Tm is the unfolding temperature, Δ
H and Δ
CP are the differences in enthalpy and heat capacity between the unfolded and the native states,
K is the conformational equilibrium constant, and Δ
G is the Gibbs energy associated with the protein stability. Fitting of Equations (1) to (3) was carried out by using user-defined routines in Origin 7.0 in order to estimate
Tm, Δ
H(
Tm), and Δ
CP.
4.4. Circular Dichroism (CD)
Circular dichroism spectra were collected on a Jasco J810 (Tokyo, Japan) spectropolarimeter fitted with a thermostated cell holder and interfaced with a Peltier unit. The instrument was periodically calibrated with (+)-10-camphorsulphonic acid.
4.4.1. Far-UV Spectra
Isothermal wavelength spectrum of N-MDM2 at different pH values or urea or GdmCl concentrations were acquired at a scan speed of 50 nm/min with a response time of 4 s and averaged over six scans at 25 °C in a 0.1 cm-pathlength cell. Step resolution was 0.2 nm and the bandwidth was 1 nm. Protein and buffer concentrations were the same used in intrinsic fluorescence experiments. Spectra were corrected by subtracting the baseline in all cases. The chemical and pH denaturations were repeated at least three times with new samples. The samples were prepared the day before and left overnight at 5 °C to allow for equilibration.
4.4.2. Thermal Denaturation Experiments
The experiments were performed at constant heating rates of 60 °C/h and a response time of 8 s. Thermal scans were collected in the far-UV region by following the changes in ellipticity at 222 nm from 25 to 80 °C in 0.1 cm pathlength cells, with a total protein concentration of 20 µM. Data were collected every 0.2 °C. Solution conditions were the same as those reported in the steady-state experiments. No difference was observed between the scans aimed to test a possible drift of the signal of the spectropolarimeter. Thermal-denaturations were not reversible at any pH for N-MDM2, as shown by (i) the comparison of spectra before and after the heating; and (ii) the changes in the voltage of the instrument [
44].
4.5. Size Exclusion Chromatography (SEC)
This technique was used to determine the compactness of N-MDM2 [
45,
46], at different pH values and GdmCl concentrations, by using a ~20 µM protein concentration. Samples were loaded in 50 mM of the corresponding buffer with 150 mM NaCl (to avoid interactions with the column) and 2 mM EDTA (to avoid polypeptide degradation) in a calibrated analytical Superdex 75 10/30 HR FPLC column (GE Healthcare, Barcelona, Spain) connected to an AKTA-FPLC (GE Healthcare) at 20 °C. For GdmCl denaturation experiments, the buffer was 50 mM (Tris, pH 8.0) with 150 mM NaCl and the corresponding concentration of GdmCl. The exclusion molecular weight of the column was 70 kDa for a globular protein, according to the manufacturer and our own measurements with blue dextran. The elution volumes,
Ve, of N-MDM2 under different pH values were obtained from chromatogram analyses with the UNICORN software version 5.01 (GE Healthcare) from three independent measurements with freshly prepared samples. Samples of N-MDM2 previously dissolved in the corresponding running buffer were eluted at 1 mL/min and continuously monitored with an on-line detector at a wavelength of 280 nm. The samples were loaded from ice, maintaining the minimum delay possible before the injection. The standards used for column calibration, at 20 µM of protein concentration, were ribonuclease A, chymotrypsinogen A, ovalbumin, and bovine serum albumin. The
Ve values of the standard proteins were measured at pH 7.0 (phosphate buffer, 50 mM) with 150 mM NaCl.
4.6. Differential Scanning Calorimetry (DSC)
The thermal stability of N-MDM2 was also assessed by differential scanning calorimetry in an Auto-PEAQ-DSC instrument (MicroCal, Malvern-Panalytical). The partial heat capacity of the protein solution was measured as a function of the temperature at a scanning rate of 1 °C/min under a pressure of 3.4 bars (50 psi) to avoid bubble formation and evaporation at high temperatures. Assays at protein concentrations of 14 and 25 μM were carried out, which correspond to 0.21 and 0.37 mg/mL, respectively. Several buffer-buffer scans were performed to ensure proper instrument and sample equilibration. After instrumental buffer baseline subtraction and concentration normalization, a baseline calculated from the progress of the unfolding process was subtracted to obtain the excess molar heat capacity of the protein. Because the unfolding traces showed two clear transitions, they were analyzed according to a two-transition unfolding model (i.e., considering two energetic domains), from which the unfolding temperature and enthalpy for each transition could be estimated. Very briefly, the partition function for a protein exhibiting two thermal unfolding transitions is given by
where
Ki(
T) is the conformational equilibrium constant for each of the two domains, and the average enthalpy, <Δ
H>, can be calculated:
where
Pi is the population of the corresponding conformational state (
P1: intermediate 1;
P2: intermediate 2;
P12: unfolded), and Δ
Hi is the difference in enthalpy between each conformational state and the native state (i.e., native state taken as a reference). The following standard relationships are considered to calculate <Δ
H>:
and the average excess molar heat capacity, <Δ
Cp>(
T), is calculated as the temperature derivative of <Δ
H>:
which is the final relationship employed for analyzing the calorimetric thermograms. Fitting was carried out by using user-defined routines in Origin 7.0 to estimate
Tm,i, Δ
Hi(
Tm,i), and Δ
CP,i. The populations of the four conformational states (native, intermediate 1, intermediate 2, and unfolded protein) were calculated by using the estimated parameters for each transition, as indicated above.
It could be thought that the model used is a three-state model, instead of a two-transition unfolding model. However, there is a marked conceptual difference between these two models: the three-state model considers the native state, the unfolded state, and an intermediate partially unfolded state, whereas the two-transition model considers two intermediate partially unfolded states. In mathematical terms, the three-state model is represented by a partition function containing three terms: 1 + K1 + K1K2, whereas the two-transition model is represented by a partition function containing four terms: 1 + K1 + K2 + K1K2. The two models may coincide when one of the intermediate states has a very low population, which occurs when the two unfolding temperatures (Tm1 and Tm2) are very far apart (i.e., K2 can be neglected compared to K1). However, both from a theoretical and operational point of view, we preferred to apply the two-transition model since it does not require any constraint on the intrinsic stability of the domains forming the protein.
4.7. Analysis of the pH, Thermal, and Chemical Denaturation Curves
The pH denaturations were analyzed, when possible, assuming that both protein species, protonated and deprotonated, contributed to the measured spectral properties:
where
X is the spectral property being observed (either ellipticity at a particular wavelength or fluorescence intensity at a particular wavelength, or else <λ>);
Xa is the intrinsic value of that property for the acidic species;
Xb is the intrinsic value of that property observed at high pH values; p
Ka is the apparent pH midpoint of the titrating group; and
n is the Hill coefficient (which was close to 1 in all the curves reported in this work).
The thermal and chemical denaturation data for N-MDM2 were fitted to
where
R is the gas constant and
T is the temperature (in K). The
XN and
XD correspond to the intrinsic values of a physical property
X of the native and denatured protein, respectively, being monitored. Both parameters showed a linear relationship with the temperature or the denaturant concentration.
Although chemical denaturations of N-MDM2 were not reversible, in order to have a qualitative way to compare the different curves obtained by the several techniques used, chemical denaturation curves for N-MDM2 were analyzed according to the linear extrapolation model [
34], in which the Gibbs energy is given by Δ
G =
m([
D]
1/2 − [
D]), where [
D] is the denaturant concentration, [
D]
1/2 is that at the midpoint of the transition, and
m is the slope of the denaturant dependence of the Gibbs energy.
The thermal denaturations, either followed by fluorescence (intrinsic or ANS) or CD, were irreversible as well. Nevertheless, we obtained an apparent thermal denaturation midpoint,
Tm, to allow for an estimate of the stability of N-MDM2 at the different pH values, where it was possible, and for a comparison with the DSC data from the change in Gibbs energy, Δ
G, given by
where ∆
Hm is the van’t Hoff unfolding enthalpy and ∆
Cp is the heat capacity change in the folding reaction. The shape of Equation (10) does not impose restrictions on the value of the Δ
Cp used in the fitting.
Fitting by non-linear least-squares analysis to Equations (8) and (9), when it was possible, was carried out by using the KaleidaGraph (Abelbeck software, Reading, PA, USA) on a PC computer.
4.8. Simulated Thermal Unfolding
The thermal unfolding of N-MDM2 was simulated by using the CNA algorithm for rigidity analysis of a protein structure [
19], implemented in the homonymous web server [
21]. Default simulation parameters were used, except when explicitly stated otherwise. With the aim of reducing the bias in the choice of the starting protein structure, we used residues 6–125 of the conformation modeled by AlphaFold (entry Q00987) [
22], which summarizes multiple structures of N-MDM2 deposited in the Protein Data Bank (PDB) [
47].
The protein in the CNA representation consists of a body-and-bar network, with such bodies and bars modeling the macromolecule atoms and their constraints, respectively. Each atom has a maximum of six degrees of freedom, which decrease depending on the interactions it is involved in. Non-covalent interactions such as hydrogen bonds and hydrophobic tethers correspond to five and two bars, respectively. By progressively removing these non-bonded constraints, the denaturation process of the protein is mimicked. In particular, hydrogen bonds are discarded according to their estimated energy [
19] from the weakest to the strongest. The removal proceeded at increasing energy E, with steps of 0.1 kcal mol
−1, so that only hydrogen bonds with a strength E
HB more favorable than such threshold value (i.e., |E
HB| ≥ E) are included in the network of each given protein state. Therefore, E represents the reaction coordinate of the simulated denaturation process.
We note that, in our definition, E has an opposite sign compared to the quantity given as output by the CNA algorithm, which is the cutoff energy E
cut (corresponding to removing hydrogen bonds with E
HB ≤ E
cut). Our convention has some advantages: (i) E is always positive; (ii) it can be more intuitively associated with a corresponding temperature T as the latter has also a positive value; and (iii) the interpretation of the graphs is clearer as the unfolding reaction proceeds from the left to the right along the
x-axis. The conversion from E to T is obtained by using the linear relationship by Radestock and Gohlke [
26], again considering the opposite sign between E and E
cut. Phase transition(s) at distinct temperature(s) can be identified while the network of constraints relaxes, going from more rigid to predominantly flexible states. This is visible in a variety of physical quantities that can be estimated, which reflect the degrees of freedom of the protein, its inner rigidity, the configuration entropy of the network, and other indexes that can be monitored both on a global and per-residue level (for details on how they are defined, we refer to [
19,
21,
23] and references therein).
The analysis of the constraint network performed by using CNA, including the temperature at which phase transitions take place, can be very sensitive to the input structure of the protein. This can be improved by analyzing a structural ensemble provided from another simulation technique such as classical molecular dynamics (MD) [
48] or, as an alternative, by producing an ensemble of networks starting from the single structure provided [
28]. The CNA web server offers an option to generate an ensemble of network topologies from the input structure by considering so-called “fuzzy” contributions that mimic non-covalent constraints randomly breaking or forming already in the native protein state [
28]. We exploited this option and performed our analysis on 500 network topologies, which is ten times larger than the number suggested as default.
5. Conclusions
MDM2 is a protein that can sample a large variety of conformations. Among the only three well-folded domains it contains, the N-terminal region N-MDM2 is the largest and plays a key role in the interaction of the whole protein with several molecular partners. In spite of being well-folded, N-MDM2 has its own distinctive features in terms of structural flexibility and plasticity. In this work, we have studied how the flexibility of N-MDM2 reflects its structural stability and, conversely, how the unfolding process encodes significant properties on the behaviour of this domain at physiological temperatures. The results obtained provide an overall picture on the sensitivity of N-MDM2 to changes in the surrounding environment, including temperature, pH, and the presence of chemical denaturants. In particular, N-MDM2 shows a very low stability and a hierarchical structure in the denaturation process that suggests a delicate organization of its sub-structures with the presence of various intermediates.
The overall findings suggest that a large conformational flexibility is an intrinsic feature of the N-terminal region of MDM2. Therefore, it is immediate to speculate that these properties are crucial for the function of N-MDM2 and especially for the interaction with p53 and the other proteins belonging to the same regulatory network. In addition to providing insights on an such isolated domain, our studies may set the basis to exploit this knowledge in some applications—as an example, to find new inhibitors of the interaction between p53 and N-MDM2 by using thermal stability measurements. More generally, the properties of N-MDM2 can also be fundamental for the interaction with additional molecular partners unrelated to p53, many of which have been uncovered in recent years.