Prerpints.org logo
Preprint
Article

A FRET Based Assay for the Identification of PCNA Inhibitors

Altmetrics

Downloads

137

Views

62

Comments

0

A peer-reviewed article of this preprint also exists.

supplementary.zip (1.48MB )

Submitted:

06 July 2023

Posted:

06 July 2023

You are already at the latest version

Alerts
Abstract
Proliferating cell nuclear antigen (PCNA) is the key regulator of human DNA metabolism. One important interaction partner is p15, involved in DNA replication and repair. Targeting the PCNA-p15 interaction is a promising therapeutic strategy against cancer. Here, a Förster resonance energy transfer (FRET)-based assay for the analysis of PCNA-p15 interaction was developed. Next to the application as screening tool for the identification and characterization of PCNA-p15 interaction inhibitors, the assay is also suitable for the investigation of mutation-induced changes in their affinity. This is particularly useful for analysing disease associated PCNA or p15 variants at the molecular level. Recently, the PCNA variant C148S has been associated with Ataxia-telangiectasia-like disorder type 2 (ATLD2). ATLD2 is a neurodegenerative disease based on defects in DNA repair due to an impaired PCNA. Incubation time dependent FRET measurements indicated no effect on PCNAC148S-p15 affinity, but on PCNA stability. The impaired stability and increased aggregation behaviour of PCNAC148S was confirmed by intrinsic tryptophan fluorescence, DSF and AF4 measurements. The analysis of the disease associated PCNA variant demonstrated the versatility of the interaction assay as developed.
Keywords: 
Subject: Medicine and Pharmacology  -   Medicine and Pharmacology

1. Introduction

Controlled DNA metabolism is crucial for cell survival. The key regulators in DNA metabolism are sliding clamps that encircle the DNA helix. Sliding clamps are ubiquitously present in all living organisms and in some viruses [1]. Despite their low sequence homology and varying oligomeric states, all sliding clamps show a high overall structure similarity with a pseudo six-fold symmetry [2,3,4,5]. They function as binding platform for proteins that are involved in DNA metabolism and tether them to the DNA at the replication fork. Additionally, sliding clamps increase the processivity of DNA polymerases [6,7].
Human DNA metabolism is controlled by the 86 kDa homotrimeric proliferating cell nuclear antigen (PCNA) [8,9]. Through the interaction with a variety of binding partners, PCNA is involved in DNA synthesis, replication, repair, chromatin assembly, chromatin maintenance and cell cycle control [9,10]. Many proteins bind PCNA via their PCNA-interacting protein (PIP)–box or PIP degron that targets PCNA for degradation. The binding occurs at the front face of PCNA between two domains of one subunit [11,12,13]. There are three binding sites for PIP-box interactions, one at each monomer. Protein binding is strongly regulated by post-translational modifications, binding affinities and local protein concentrations [9]. Other interaction motifs are the KA-box [14] or the AlkB homologue 2 PCNA-interacting motif (APIM), which overlap with the binding site of PIP-box proteins [15,16].
Due to the importance of PCNA for DNA metabolism, an altered function can cause drastic consequences for the cell. The expression of PCNA is elevated in many types of tumours and is associated with high proliferation rates [17,18]. Targeting PCNA by interfering with the PCNA interaction network, preventing trimerization or reduction of chromatin stability are promising therapeutic strategies. The clinical relevance of PCNA led to the identification of peptidic, small molecule, and aptamer inhibitors of PCNA interactions, which efficiently reduce cell growth and induce cell death [18,19]. A promising example of a PCNA interaction inhibitor is ATX-101, which is currently in phase II clinical studies. ATX-101 is a cell-penetrating peptide based on the APIM binding motif. It reduces primary metabolism, inhibits DNA repair, induces apoptosis and increases the effect of chemotherapeutic agents in many tumour types and cancer cell lines, while healthy cells were less affected [20,21,22]. Another example of a PCNA inhibitor targeting the PCNA binding groove is T2 amino alcohol (T2AA) [23]. T2AA disrupt the interaction of PCNA to PIP-box containing proteins and increases the effects of DNA damaging agents [23,24,25]. Molecules that bind to the hydrophobic binding pocket of PCNA affect a large network of protein-protein interactions in multiple signalling pathways. This can lead to undesired side effects. The inhibition of one specific therapeutically relevant interaction might circumvent this problem. The PCNA-p15 interaction is particularly well suited for this purpose [26]. p15 is an intrinsically disordered protein involved in DNA replication and repair by its interaction with PCNA [27]. Increased expression of p15 is associated with cancer and chemotherapeutic resistance [28,29,30,31]. p15 interacts with its canonical PIP-box sequence (62QKGIGEFF69) at the front site of PCNA. But what makes the p15 interaction special is that the p15 residues 59-62 form a β-turn, which direct the N-terminus of p15 towards the channel of the ring. This allows the interaction of the residues 51NPVCVRPT58 to the inner side of the PCNA ring [26]. The unique binding mode allows the specific targeting of the therapeutically relevant PCNA-p15 interaction.
Although an elevated PCNA expression correlates with highly proliferative cells and is associated with tumours, low PCNA amounts or an impaired functionality can also be fatal for the cell. PCNA knockout in mice is embryonically lethal [32], indicating the essential role of PCNA for the cell. It is hardly surprising that almost no disease associated PCNA mutations have been discovered so far. However, in 2014, the PCNAS228I variant was linked to ataxia-telangiectasia-like disorder type 2 (ATLD2) [33]. ATLD2 is a neurodegenerative disorder with underlying defects in DNA repair based on an impaired PCNA. The main clinical features are telangiectasia, sensorineural hearing loss, premature aging, short stature, photosensitivity and neurodegeneration [33]. Further symptoms include developmental delay, foot deformity, speech and swallowing deficits, muscle weakness and mild oculomotor apraxia [33,34]. Previous studies showed that in patient-derived cells an unchanged PCNAS228I protein level and an unaffected DNA replication is observed. However, the cells show an enhanced sensitivity against UV radiation indicating defects in DNA repair [33]. Structural studies with PCNAS228I observed a large conformational change at the interdomain connection loop (IDCL) and at the PIP-box binding groove [35], which likely impairs the interaction of several PIP-box containing proteins [33,35]. Recently, the PCNA variant NM_002592.2(PCNA): c.443G > C(p.C148S) was associated with ATDL2 [34]. The molecular consequences of this mutation were not understood during the experiments of this publication. In the meantime, however, a drastic decrease in the thermal stability of PCNAC148S has been demonstrated [36].

2. Results

2.1. Development of a FRET assay to analyse the PCNA-p15 interaction.

The PCNA-p15 interaction presents a promising target for cancer therapy. For screening PCNA-p15 inhibitors targeting the unique interaction site, a FRET based binding assay was established. For this purpose, p15 and PCNA were genetically conjugated with a FRET donor-acceptor pair. The p15-mNeonGreen and mScarlet-I-PCNA fusion proteins were expressed in E. coli BL21(DE3) cells and purified by affinity chromatography afterwards. For analysis of PCNA-p15 interaction by FRET, 1 µM p15-mNeonGreen was mixed with varying concentrations of mScarlet-I-PCNA (0 - 2.125 µM) and incubated at 37 °C for 1 h. Fluorescence was then excited at 488 nm and the emission spectra were measured at 500 – 700 nm (Figure 1B). The detected emission maxima correspond to the emission of mNeonGreen at 517 nm and mScarlet-I at 594 nm [37,38]. Due to the wide excitation range of mScarlet-I, the acceptor is also directly excited by donor excitation. In addition, the emission spectrum of the donor overlaps with that of the acceptor. Accordingly, three emissions were detected for mScarlet-I-PCNA-p15-mNeonGreen when excited at 488 nm: direct donor molecule emission, direct acceptor molecule emission, and sensitized emission. The emission at 594 nm is elevated for mScarlet-I-PCNA/mNeonGreen-p15 in comparison to the non-binding control (mNeonGreen/mScarlet-I-PCNA). Simultaneously, the donor emission is reduced with increasing acceptor concentrations (donor quenching). The occurrence of both observations proves the energy transfer from mNeonGreen to mScarlet-I upon PCNA-p15 interaction. Thus, the assay setup is suitable for quantitative analysis of PCNA-p15 interaction and inhibition.
For the analysis of the PCNA-p15 affinity, the samples were prepared as in the experiment before. Binding was evaluated by analysing the sensitized emission (SE) using the three-channel method, which considers the spectral crosstalk of the fluorophores [39,40]. By plotting the SE against the acceptor concentration, a dissociation constant of 430.8 ± 59.9 nM for the PCNA-p15 interaction was determined (Figure 2). This is in good agreement with the binding affinity already published for PCNA-p15 interaction (isothermal titration calorimetry, 1.1 µM) [26].
To investigate the ability to detect inhibitory molecules in this assay, the effect of nonfluorescent PCNA, the peptide p1551-70 and the small molecule inhibitor T2AA were tested. Unlabeled PCNA competes with mScarlet-I-PCNA for the binding to p15-mNeonGreen. Therefore, a reduced binding signal of p15-mNeonGreen/mScarlet-I-PCNA in presence of unlabeled PCNA is expected. p1551-70 is a peptide based on the sequence of p15, which binds competitively to p15 at PCNA. T2AA is an already published PCNA inhibitor interacting with the PIP-box binding groove [23] and should inhibit the PCNA-p15 interaction. To test the inhibitory effect of the molecules, 1 µM mScarlet-I-PCNA and 1 µM p15-mNeonGreen were incubated with varying concentrations of the inhibitors. By plotting the inhibitor concentration against the sensitized emission, IC50 values of 1.23 ± 0.06 (unlabeled PCNA), 9.81 ± 0.99 (p1551-70) and 13.81 ± 2.0 µM (T2AA) were determined. Ki values were calculated using the Cheng Prusoff equation and turned out to be 0.38 ± 0.02 (unlabeled PCNA), 2.99 ± 0.26 (p1551-70) and 4.21 ± 0.55 µM (T2AA), respectively. Assuming that the conjugation of mScarlet-I does not affect the binding affinity of PCNA-p15, the expected Ki value is equal to the determined KD value [41]. The observed KD value of 0.43 µM is quite similar to the Ki value actually obtained in this assay, demonstrating the accuracy of the method. T2AA inhibits the interaction of the PIP-box containing PL-peptide with PCNA at an IC50 of 1-1.5 µM [24,25] and a Ki value of 0.91-1.36 µM was calculated. Thus, T2AA-induced inhibition of the PCNA-p15 interaction is in the same concentration range as the inhibition of the PCNA-PL peptide interaction. Accordingly, this assay is suitable for the identification and quantification of inhibitory peptides and small molecules targeting the PCNA-p15 interaction.

2.2. Screening for PCNA inhibitors

To identify new classes of PCNA inhibitors, compounds of the in-house database were analysed regarding their effect on the PCNA-p15 interaction. Among them derivates of pyranobezochinones, aminomethylfuropyranones, totarols and benzofuran-4,7-diones were tested. However, no promising results were obtained.

2.3. Analysing the relevance of the mutation C148S for PCNA function

Next to the application as inhibitor screening tool, the FRET assay was also accurate in determining the binding affinity of PCNA-p15. Thus, mutational effects on this interaction can be investigated. Recently, the mutation C148S was associated with the disease ATLD2 [34,36]. However, nothing was known about the consequences of this mutation for PCNA function. C148 is highly conserved in different eukaryotic species (Table S1), indicating the importance of the cysteine at position 148 for PCNA. C148S is located at the inner site of the PCNA ring, far away from the typical binding site of most PCNA interaction partners (Figure 4). Its neighboring amino acid R149 is involved in sliding of PCNA along the DNA helix as well as for p15 binding [26,43]. We employed the established FRET assay to determine the influence of the C148S mutation on the interaction between PCNA and p15.

2.3.1. C148S reduces the apparent binding of PCNA to p15

The influence of C148S on PCNA-p15 binding was analysed by the FRET assay. p15-mNeonGreen was used as donor and varying concentrations of mScarlet-I-PCNAWT or mScarlet-I-PCNAC148S were used as acceptor molecule. The sensitized emission was plotted against the acceptor concentration and KD values of 0.43 ± 0.06 µM for p15/PCNAWT and 3.09 ± 1.11 µM for p15/PCNAC148S were ascertained. Thus, an approximately 7 times lower apparent binding of PCNAC148S in comparison to the wild-type protein was observed. However, saturation must be approached to accurately determine KD values from saturation binding experiments [44]. This was not the case for the p15/PCNAC148S interaction, resulting in an uncertainty of the obtained KD value. Nevertheless, it can be concluded that the incorporation of C148S strongly impairs the apparent binding of PCNA with p15.
Preprints 78749 g005
Figure 5. Analysis of the effect of C148S on the PCNA-p15 binding by FRET. 1 µM p15-mNeonGreen was incubated with 0 – 3.09 µM mScarlet-I-PCNA (black circles) or mScarlet-I-PCNAC148S (grey squares) for 1 h at 37 °C. The sensitised emission (SE) of three independent experiments was plotted against the acceptor concentrations, resulting in KD values of 0.43 ± 0.06 µM for PCNAWT and 3.09 ± 1.11 µM for PCNAC148S. The given error is the standard error of the fit. RFU = relative fluorescence units.
Figure 5. Analysis of the effect of C148S on the PCNA-p15 binding by FRET. 1 µM p15-mNeonGreen was incubated with 0 – 3.09 µM mScarlet-I-PCNA (black circles) or mScarlet-I-PCNAC148S (grey squares) for 1 h at 37 °C. The sensitised emission (SE) of three independent experiments was plotted against the acceptor concentrations, resulting in KD values of 0.43 ± 0.06 µM for PCNAWT and 3.09 ± 1.11 µM for PCNAC148S. The given error is the standard error of the fit. RFU = relative fluorescence units.
Preprints 78749 g005

2.3.2. Alchemical free-energy calculations predict a destabilizing effect of the C148S mutation on PCNA

We next asked if the reduced apparent binding signal was caused by a direct impairment of the p15/PCNA binding site or by an affected structure or stability of PCNA. Therefore, we conducted alchemical free energy calculations. Alchemical free-energy calculations are a useful in silico method for studying the effect of mutations, as they can provide remarkable prediction accuracy, if properly set up [45,46]. It is possible to estimate the differences in the free energy (∆G) between two defined sets of states, i.e. the unfolded and folded state of a protein or the bound and free state of a complex. The effect of mutations on this free energy difference can also be given as the relative change in ∆G, which is commonly referred to as the relative folding free energy (∆∆Gfolding) or the relative binding free energy (∆∆Gaffinity). In case of PCNA, this allows to study the effect of the C148S mutation on protein stability (∆∆Gfolding) and PCNA-p15 affinity (∆∆Gaffinity).
Alchemical free-energy methods work by gradually transforming one chemical state into another following an unphysical pathway. The alchemical pathway does not need not follow the physical folding or binding pathway and is therefore more accessible by molecular dynamics (MD) simulations. In our case we made use of a non-equilibrium approach based on the Jarzynski equality and the Crooks fluctuation theorem [47,48,49]. In this approach two standard MD simulations are run, one for each state (WT and C148S) of the system. Snapshots from these equilibrium simulations are extracted and used to perform alchemical transition by driving the Hamiltonian H of the system from one state to another using a parameter called λ. This parameter is gradually increased from 0 to 1 during the simulation.
To estimate both the ∆∆Gfolding of PCNA and the ∆∆Gaffinity between PCNA and p15 upon amino acid mutation, the wild-type amino acid was alchemically morphed from a cysteine at position 148, to a serine. To keep the systems to a reasonable size, we included only one PCNA monomer, instead of the homotrimer. This drastically reduced computational cost. The stability of monomeric PCNA during 100 ns of MD simulation was tested. No great conformational change, and a root mean square deviation (RMSD) below 3 Å for the protein was observed. We therefore concluded that it was viable to only use the monomer going onward. For the calculation of ∆∆Gfolding the transition for the unbound state, in this case approximated by a Gly-Cys-Gly tripeptide, and the folded state were performed [50]. For the calculation of ∆∆Gaffinity the transitions were performed for free PCNA and the PCNA-p15 complex. The resulting ∆G values for each transition were used to calculate the ∆∆G values according to the thermodynamic cycle shown in Figure S4. We also performed calculations for two further mutations to construct a closed thermodynamic cycle with the respective ∆∆G values at its edges to confirm convergence of the simulations (Figure 6). This cycle includes mutations for serine to alanine and alanine to cysteine, forming a cycle with three edges. Furthermore, the simulation length for the transitions were optimized by running these simulations with different transition times. This optimization was performed for the mutation involving the biggest perturbation, in this case A148C. ∆∆Gfolding was monitored as a function of transition time and no significant change in ∆∆Gfolding was observed when increasing the transition time beyond 100 ps (Figure S5). We therefore concluded that a transition time of 100 ps was sufficient for the mutations and the system under investigation. Following this optimized setup, cycle closure for the calculation of ∆∆Gfolding and ∆∆Gaffinity was achieved, with all the single legs of each cycle summing up to 0.1 kcal/mol and -0.19 kcal/mol. These values are within the expected range of error for alchemical free energy calculations [45,46], and within the uncertainty of our approach as seen by the standard deviation over three independent replicas. The estimated ∆∆Gfolding value of 2.39 ± 0.47 kcal/mol for the mutation C148S suggest a destabilizing effect of the mutation. An increase in ∆∆Gfolding is equivalent to a more positive folding free energy of PCNAC148S. This leads to a decreased thermostability and/or chemical stability [51]. It furthermore increases the probability of the protein to adapt intermediate folding states, which are normally poorly populated. These intermediate states often display large hydrophobic patches on their surface which can initiate aggregation [52]. The calculated ∆∆Gaffinity of 0.21 ± 0.22 kcal/mol for the C148S mutation, suggest no direct change in binding affinity of PCNAC148S to p15 compared to PCNAWT. One likely conclusion is, that the mutation leads to a less stable variant of PCNA, which might be more prone to aggregation or denaturation. When this protein is used as the titrator in a binding assay, less functional protein than expected might be present in solution. The measured affinity would therefore be apparently lower than it actually is. This might explain the apparently reduced binding of p15 to PCNA measured with the FRET assay, even if the binding free energy is not changed by the mutation. We therefore decided to experimentally test the thermal and chemical stability of PCNAC148S, as well as its aggregation behavior.

2.3.3. The amino acid substitution C148S reduces the chemical stability of PCNA

For analysing the influence of C148S on the chemical stability of PCNA, we measured the intrinsic tryptophan fluorescence at 320 nm (Figure 7). The PCNA trimer contains three tryptophans (W28), one per monomer, allowing the analysis of the chemical stability by the intrinsic tryptophan fluorescence. The tryptophans of PCNA are buried in a hydrophobic environment, which get more exposed to the aqueous solvent during unfolding. The aqueous environment results in a blue-shifted tryptophan emission, reducing the emission at 320 nm [53].
Our results show a three-state behaviour for the urea induced denaturation of the wild-type protein. Such three-state curve indicates the presence of an intermediate state where the tryptophans are more exposed to the solvent. In comparison to the wild-type PCNA, the denaturation of PCNAC148S arises in only two states. The overlay of both curves reveals a high similarity of the tryptophan fluorescence of PCNAWT and PCNAC148S between 2 – 8 M urea (Figure 7A). Since the intrinsic fluorescence of tryptophans is strongly dependent on their environment, the tryptophans (W28) probably undergo similar conformational states between 2 - 8 M urea in both proteins.
Without the addition of urea, the tryptophan fluorescence is higher for PCNAWT than for PCNAC148S. Accordingly, the tryptophans are probably more buried in the wild-type protein under this condition. The tryptophan fluorescence intensity of PCNAC148S is constant between 0-2 M urea and similar to the fluorescence intensity of the wild-type’s intermediate state. One possible explanation might be, that PCNAC148S predominantly populates this intermediate conformational state between 0-2 M urea, even in the absence of urea. Based on these results, a decrease in chemical stability due to the incorporation of C148S was concluded.
Furthermore, the midpoint denaturation (Cm) for the complete denaturation were determined, which were at 3.95 ± 0.18 M for PCNAWT and at 3.32 ± 0.15 M for PCNAC148S. Accordingly, the higher Cm value of PCNAC148S and the different denaturation states indicate a decrease in PCNA stability due to the incorporation of C148S.

2.3.4. C148S affects the thermal stability of PCNA

The intrinsic tryptophan fluorescence measurements revealed a negative effect of C148S on the stability of PCNA. To confirm this result, the thermal stability of PCNA was analysed by differential scanning fluorimetry (DSF). In the assay, an extrinsic fluorescent dye is used whose fluorescence is quenched in an aqueous environment. Upon protein unfolding, more hydrophobic sites are accessible for the interaction with the fluorescent dye, resulting in an elevated fluorescence intensity. This allows the analysis of thermal protein denaturation based on the fluorescence intensity of the dye.
The DSF experiment revealed aberrations in the denaturation profile of PCNAC148S and PCNAWT (Figure 8). PCNAC148S shows a higher initial fluorescence signal in comparison to PCNAWT, which implies more hydrophobic patches on the exposed surface of PCNAC148S prior to heating. Such exposed hydrophobic patches might be explained by denatured or unproperly folded protein samples [54]. This observation is consistent with the interpretation of the previous experiment, where we assumed a stronger exposure of W28 to the aqueous solvent for PCNAC148S than for PCNAWT.
The fluorescence intensity of PCNAC148S peaks at lower temperatures in the normalized plot, indicating the denaturation of PCNAC148S at lower temperatures (Figure 8B). By plotting the first derivative of the melting curves against the temperature, two melting points for PCNAWT (Tm1 = 49.1 ± 0.1 and Tm2 = 54.54 ± 0.66 °C) were determined. Accordingly, similar to the previously performed chemical denaturation, a three-state denaturation profile for PCNAWT was detected. Due to the poor resolution of the PCNAC148S‘s derivative, the determination of Tm-values was not possible here.
ΔG-values were estimated from the DSF melting curves, as explained in detail in the methods section. A ΔG-value for protein unfolding of 6.54 ± 0.48 kcal/mol for PCNAWT and 4.89 ± 0.26 kcal/mol for PCNAC148S was ascertained. Thus, the experimentally determined ΔΔGunfolding (1.65 ± 0.74 kcal/mol) is similar to the ΔΔGunfolding values predicted by the alchemical free energy calculations (2.39 ± 0.47 kcal/mol). The measured positive ΔΔGunfolding value as well as the normalized melting curves strongly supports our hypothesis that C148S impairs the thermal stability of PCNA.

2.3.5. PCNAC148S exhibits an enhanced aggregation behavior

Reduced protein stability leads to more frequent folding/unfolding intermediates with exposed hydrophobic sites on the protein surface. This can initiate the formation of protein aggregates. Reduced protein stability is therefore often also associated with an increased aggregation behaviour [52]. Thus, we asked whether the introduction of the mutation C148S also has an influence on the aggregation behaviour of PCNA. The use of asymmetrical flow field-flow fractionation (AF4) enables the separation of proteins based on their hydrodynamic diameter and is suitable to detect protein aggregates [55]. Aggregation was induced by incubating PCNAWT or PCNAC148S at 37 °C for 0, 1, or 16 h or at 65 °C for 10 min (Figure 9). The molecular weight of the samples was estimated using the co-running size standard.
Based on the analysis of the molecular weight standard, both PCNAWT and PCNAC148S were mainly present in their trimeric forms after 0 h at 37°C (elution time around 10 -15 min). These data are consistent with crosslinking experiments, where trimers for both PCNAWT and PCNAC148S were detected (Figure S3). When comparing PCNAWT and PCNAC148S after 0 h at 37°C, an additional early elution peak (elution time around 5 – 10 min) is more pronounced for PCNAC148S, which might correlate to monomeric PCNA.
High temperatures induce protein unfolding and often promotes aggregation [56]. Since a Tm-value of 54.54 °C was determined for PCNAWT unfolding by DSF, we induced the aggregation of PCNA at 65 °C for 10 min. This condition resulted in later molecule elution times (above 15 min) and broader signal peaks and was used as positive control for protein aggregation (Figure 9). In AF4, later elution times correlate with higher hydrodynamic diameters of the molecules. Increasing the thermal stress by elongation of the incubation time at 37 °C elevated the protein elution time (above 15 min) for both, PCNAWT and PCNAC148S. This increase of the elution time occurred more rapidly for PCNAC148S than for the PCNAWT, which is particularly evident after protein incubation at 37 °C for 1 h (Figure 9). These high hydrodynamic diameter molecules can be explained by the formation of protein aggregates.
Accordingly, PCNAC148S exhibits a lower thermal stability and an enhanced aggregation behaviour. Since both PCNAWT and PCNAC148S were mainly present in their trimeric form after 0 h at 37 °C, we suggest that the mutation does not affect the trimerization of the protein.

2.3.6. The impaired thermal stability of PCNAC148S is responsible for its reduced binding signal to p15 measured by FRET.

The prediction of ΔΔGaffinity suggests no impairment of PCNAC148S affinity towards p15. The method proved to be quite accurate in predicting the ΔΔGunfolding. However, a higher KD-value for p15 binding to PCNAC148S in comparison to PCNAWT was determined. We therefore decided to test whether this discrepancy might be a result of reduced protein stability. The PCNA binding assays include an incubation step at 37°C for 1 h. AF4 measurements revealed that such protein incubation elevated the hydrodynamic diameter of PCNAC148S. In contrast, the hydrodynamic diameter of PCNAWT remained almost unchanged. This is also true for mScarlet-I-PCNAWT and mScarlet-I-PCNAC148S (Figure S2), which were used for the FRET bindings assay. Accordingly, a lower amount of mScarlet-I-PCNAC148S might be available for the binding assays, which might affect the apparent binding.
To check this assumption, FRET measurements were performed, in which the affinities were determined after varying incubation times at 37 °C (Figure 10). The KD-values of mScarlet-I-PCNAWT and mScarlet-I-PCNAC148S are similar after 15 min at 37 °C. However, after 30 min incubation, the KD-values of mScarlet-I-PCNAC148S to p15-mNeonGreen are increased, while the affinity of mScarlet-I-PCNAWT to p15-mNeonGreen remains constant over 60 min. The consistent KD-values of mScarlet-I-PCNAWT indicate that the binding of p15 and PCNA is already in equilibrium after 15 min incubation. The dramatic increase in the KD-values of mScarlet-I-PCNAC148S from 30 min at 37 °C suggests an effect of the incubation on the calculated binding affinities. Based on the lower stability of PCNAC148S, we conclude that the continued incubation at 37 °C results in less functional mScarlet-I-PCNAC148S that is available for p15-mNeonGreen binding. These results strengthen our confidence that the measured decreased binding affinity of PCNAC148S after 60 min at 37 °C is caused by the impaired thermal stability of PCNAC148S.

3. Discussion

3.1. Development of a FRET based assay to analyse the PCNA-p15 interaction

In this study, a FRET-based assay for the identification of specific PCNA-p15 inhibitors was developed. This assay is a convenient method for analysing the effect of proteins, peptides and small molecules on the PCNA-p15 interaction. The determined binding affinity is in the same range as the already published KD value of PCNA and p15 measured by isothermal titration calorimetry (1.1 µM) [33], indicating the accuracy of the method. However, as the p15-mNeonGreen concentration is similar to the dissociation constant, a titration effect on the determined KD value cannot be excluded [57].
The assay was performed in 384 well plates and can easily be scaled up to a 1536 format [58]. The conjugation with fluorescent proteins facilitates protein purification, obviates the need for additional labeling steps, and can increases protein stability and thus the robustness of the assay. In addition, fluorescent fusion proteins are more environmentally friendly in comparison to other hazards like radioisotopes. The PCNA-p15 FRET assay is similar to fluorescence polarization methods, that are already reported for screening of PCNA-PIP-box inhibitors [24,25,59,60]. However, the FRET assay can be applied using plate readers with simpler optics that are widely distributed in many laboratories. The assay focusses on the in vitro characterization of PCNA-p15 interaction inhibitors. In addition to the identification of PIP-box binders, this also allows the identification of specific inhibitors of the PCNA-p15 interaction. These specific inhibitors may circumvent undesired side effects due to interference with various DNA metabolism processes. However, the FRET assay does not allow the direct differentiation between PCNA-PIP box inhibitors and specific PCNA-p15. This question can be investigated by testing the inhibitory activity of a promising compound for further PIP-box binding partners.

3.2. The amino acid substitution C148S reduces the amount of functional PCNA

In addition to the application as screening tool for PCNA inhibitors, we also successfully used the FRET assay for analysis of mutational effects on PCNA-p15 binding. In this study we investigated the functional consequences of the newly described PCNA mutation C148S, which is related to ATLD2. Although we initially assumed an effect of the C148S incorporation on PCNA-p15 binding, we now suggest that the mutation does not affect PCNA-p15 affinity. This conclusion is consistent with the recently published data of Magrino et al., in which no impairment of the binding affinity of PCNAC148S to PIP-box interaction partners was detected [36]. The incubation time dependent analysis of the KD-values and the free energy calculation for PCNA-p15 affinity indicate that the previously measured decreased binding of PCNAC148S is rather related to a reduced amount of natively folded PCNAC148S than to its lower affinity. When a protein with a reduced stability is used as the titrator in a binding assay and were incubated for 1 h at 37 °C, less functional protein is present in solution. The measured affinity would therefore be apparently lower than it actually is. The incubation time dependent FRET measurements support this hypothesis as the KD value of PCNAWT was constant over 1 h, while the KD value of PCNAC148S increased after continuous incubation at 37 °C. The decreased amount of functional PCNA can be explained by the reduced chemical and thermal stability of PCNAC148S, which was also confirmed by Magrino et al [36]. Destabilizing mutations are more likely to have intermediate folding states with external hydrophobic regions, which are normally marginally populated. An increased frequency of an intermediate state of PCNAC148S was indicated by the chemical denaturation experiment. There, W28 of PCNAC148S was predominantly present in the intermediate conformation even without the addition of urea. Further, the high initial fluorescence values of PCNAC148S in the DSF measurement also indicated an increased exposure of hydrophobic amino acids on the protein surface. Such external hydrophobic sites often initiate protein aggregation [52,61], explaining the dramatically increased aggregation behaviour of PCNAC148S measured by AF4. However, the affinity measurements, the detection of PCNAC148S-trimers by AF4 and the unchanged crystal structure of PCNAC148S prove that the intermediate PCNA state is at least partially functional.
Now, the question arises why a substitution of cysteine to serine impairs the stability of PCNA. PCNA does not exert any stabilizing disulfide bonds [62] that could be prevented by the amino acid exchange. However, it has been shown that a cysteine to serine substitution reduces the stability and increases the aggregation behaviour of proteins even if they are not involved in disulfide bonds [63,64]. C148 is buried in a hydrophobic environment and its sulphur atom forms a Π-bond with the backbone of F144 [36,62]. Although there is a high structure similarity between serine and cysteine, buried cysteines are more hydrophobic, which might explain the destabilizing effect of serine at this position. The free energy calculations indicate that the exchange to both serine and alanine destabilizes PCNA. Accordingly, the highly conserved cysteine at position 148 is important for PCNA stability.

3.3. Consequences of the reduced PCNAC148S stability for the cell

The functional consequences of PCNAC148S for the cell are not fully understood yet. Based on our data and the reduced amount of chromatin bound PCNAC148S in patient derived cells [36], we are suggesting that the mutation is hypomorphic in nature. PCNA downregulation causes drastic consequences for the cell. It abolishes DNA synthesis and cell proliferation and impairs the genomic integrity and DNA repair [65,66,67]. However, it is unlikely that PCNAC148S causes such strong effects, because it would probably not be compatible with life. To date, only one other human disease-related PCNA variant has been described, which also underlies ATLD2 [33]. The common component that unites both ATLD2 associated PCNA variants is the dramatic loss of the stability [36]. The two ATLD2 associated PCNA variants PCNAS228I and PCNAC148S differ with respect to their binding affinities to PIP-box binding partner. The impairment of the binding affinity of PCNAS228I is probably caused by the altered structure of PCNAS228I, while the structure of PCNA is not affected by C148S incorporation [35,36]. PCNAS228I has no effect on DNA replication, but on DNA repair [33]. PCNA and its binding partners are present at high levels during DNA replication [68,69]. Green et al. suggested that the high level of PCNA at the replication fork might compensate the reduced binding affinity of PCNAS228I during DNA replication [70]. It could be similar for PCNAC148S. The lower stability of PCNAC148S might be less important during DNA replication due to its high expression level. The effect of the C148S amino acid substitution might also be fatal for DNA repair. Considering that symptoms of ATLD2 patients are already linked with defects in DNA repair [33], it might be reasonable that the C148S also affects DNA repair. The impact of C148S on DNA repair was already suggested by the higher sensitivity of patient cells to IR or zeocine induced DNA damage [36].
Beyond that, the increased aggregation behaviour of PCNAC148S might contribute to the disease. In addition to the loss of protein function, aggregation can affect cellular protein homeostasis [71]. Many disease-associated protein variants show an increased aggregation, indicating the relevance of aggregation on various human diseases [72]. Future work on analysing the effects of PCNA aggregation on the cells, might give further insights for understanding of the disease.

3.4. PCNA stabilizing molecules as therapeutic option

In accordance with Magrino et al. [36],we suspect that the disease is essentially driven by the lower stability of PCNAC148S. A possible therapeutic approach might be the identification of a PCNA stabilizing molecule. Such pharmaceutical chaperones continue gaining importance for the treatment of destabilizing missense mutations [73,74]. Due to the quantities of PCNA interaction partners [9,10], a stabilizer might interfere with the interaction of other binding partners. However, most PCNA interaction partners bind at similar sites at the front site of PCNA [13,14,15]. Accordingly, the identification of an allosteric stabilizer might be a promising approach.

4. Materials and Methods

4.1. Plasmid construction and mutagenesis

A list of all used oligonucleotides is given in Table S2. The construct PET28b-PCNA, pET11d-p15 have been kindly provided by Alfredo de Biasio (King Abdulla University of Science and Technology, Thuwal, Saudi Arabia). For the generation of the p15-mNeonGreen plasmid, the sequence of mNeonGreen from Branchiostoma lanceolatum [38] was integrated C-terminal of the p15 sequence, connected by a flexible linker (GGGGS). The incorporation of the mNeonGreen sequence was done by In-Fusion Cloning (Clontech Laboratories, Takara, Saint- Germain-en-Laye, France). Analogously, the sequence of mScarlet-I [37] was fused to the N-terminus of PCNA for the construction of the mScarlet-I-PCNA expression plasmid. mScarlet-I and PCNA were connected by the hydrophilic flexible linker (GEGQGQGPGRGYAYRS) [75]. The C148S substitution was introduced in the PCNA sequence by site directed mutagenesis.

4.2. Protein purification

E. coli BL21(DE3) cells containing the corresponding plasmid were cultivated at 37 °C to an optical density (OD578nm) of 0.6. For induction of PCNA or mScarlet-I-PCNA expression, 1 mM IPTG were added and the cells were incubated for 16 hours at 23 °C. After cultivation, the cell pellet was resuspended in buffer A (50 mM Tris-HCl, pH 7.6, 150 mM NaCl, 10% glycerol, 1 mM benzamidine, 1 mM PMSF). Cells were disrupted by sonification and the solution was centrifuged at 50,000 x g for 1 h at 4 °C afterwards. The lysate was added to the nickel NTA column, previously equilibrated with buffer B (50 mM Tris-HCl, pH 7.6, 150 mM NaCl, 10% glycerol). The protein elution was induced by increasing concentrations of imidazole in buffer B (50-500 mM). The quality of the purification was analysed by SDS-PAGE and Coomassie brilliant-blue staining. The PCNA containing fractions of interest were dialyzed against phosphate-buffered saline (PBS), aliquoted, frozen with liquid nitrogen and stored at -80 °C.
The expression of p15-mNeonGreen was induced by the addition of 1 mM IPTG and the cells were incubated at 30 °C for additional 4 hours. The cell pellet was resuspended in buffer C (20 mM Tris-HCl, pH 8, 150 mM NaCl, 10% glycerol, 1 mM benzamidine, 1 mM PMSF) and the cells were disrupted by sonification on ice. The solution was centrifuged at 50,000 x g for 1 h at 4 °C and the supernatant was transferred to a column equilibrated with buffer D (20 mM Tris-HCl, pH 8, 150 mM NaCl, 10% glycerol). p15-mNeonGreen was eluted with increasing concentrations of imidazole (50-500 mM) in buffer D. The quality of the purification was analysed by SDS-PAGE and Coomassie brilliant-blue staining. The protein fractions of interest were dialyzed against PBS containing 1 mM DTT, aliquoted, frozen with liquid nitrogen and stored at -80 °C.
The protein concentration of p15-mNeonGreen was determined using a nanophotometer (Pearl, Implen GmbH, Munich, Germany) and the extinction coefficient 52830 M-1cm-1 at 280 nm. The concentration of PCNA was determined as promoter concentration using the extinction coefficient at 280 nm (15930 M-1cm-1 for His-PCNA and 53290 M-1cm-1 for His-mScarlet-I-PCNA).

4.3. FRET binding assay

For FRET-KD-assay 1 µM p15-mNeonGreen and variable concentrations of mScarlet-I-PCNA (0-3.1 µM) were mixed in PBS containing 1 mM DTT. For initial testing of inhibitors, 1 µM p15-mNeonGreen, 1 µM mScarlet-I-PCNA and 50 µM of the inhibitors were used. IC50 values were determined using 1 µM p15-mNeonGreen, 1 µM mScarlet-I-PCNA and varying inhibitor concentrations. The assay was performed in a final volume of 20 µL per well in a 384 well microtiter plate (3766, Corning, NY, USA). After incubation at 37 °C for 60 min at 300 rpm, the fluorescence of the samples was detected using the microplate reader infinite 200Pro (Tecan, Männedorf, Switzerland). The fluorescence was measured in three channels: 488/525 nm, 561/610 nm and 488/610 nm and the FRET effect was determined by evaluating the sensitized emission [39]. mNeonGreen as well as mScarlet-I were used as non-binding controls, respectively (Figure S1). The sensitized emission of mNeonGreen/mScarlet-I-PCNA was subtracted from the sensitive emission of the p15-mNeonGreen/mScarlet-I-PCNA interaction. For the determination of the binding affinities, the sensitized emission of three experiments was measured in technical triplicates and plotted against the acceptor concentration. The KD-values were calculated using following KD-Fit [39] (equation 1) and GraphPad Prism 5 (GraphPad Software, La Jolla, USA):
E m F R E T = E m F R E T m a x ( 1 − 2 K D X − a + K D + ( X − a − K D ) 2 + 4 K D X ) ,
where EmFRET is the FRET signal intensity, EmFRETmax is the maximum of EmFRET and a is the constant donor concentration.
For the determination of the inhibitory activities, the sensitized emission of three experiments was measured in technical triplicates, normalized and plotted against the inhibitor concentration. The IC50-values were calculated using the log (inhibitor) vs. response- variable slope (four parameters) of GraphPad Prism 5 (GraphPad Software, La Jolla, USA):
The small molecule positive control T2AA was purchased from Sigma Aldrich (St. Louis, USA). The peptide p1551-70 containing the amino acids NPVCVRPTPKWQKGIGEFFR was commercially acquired from GenicBio Limited (Shanghai, China).

4.6. Differential Scanning Fluorimetry (DSF)

The thermal stability of PCNAWT and PCNAC148S was determined by DSF [76]. 4 µM PCNA were mixed with 20 x SYPRO Orange Protein Gel Stain (Sigma Aldrich, St. Louis, USA) in PBS. The samples were heated from 30-95 °C by 0.5 °C/5 s in the Rotor-Gene Q 2 Plex HRM (Qiagen, Hilden, Germany). SYPRO orange was excited at 470 ± 10 nm and the emission was detected at 610 ± 5 nm. The measurement was conducted three times in quadruplicates. The melting points were obtained using the Rotor-Gene Q software by determining the extreme points of the first derivative of the fluorescence dF/dT(T) formed from the melting curves as function F(T). ∆Gunfolding was estimated from the melting curves with the method described by Wright et al. [77]. In short, the fraction of folded protein, Pf, at any given temperature is calculated with equation 2. Where Fmin is the minimal measured fluorescence, Fmax is the maximal fluorescence and F is the measured fluorescence at the given temperature.
P f = 1 − F − F m i n F m a x − F m i n
Next the fraction of unfolded protein, Pu is calculated with equation 3 at each temperature.
P u = 1 − P f
Ku, the equlibrium constant of unfolding is calculated with equation 4.
K u = P u P f
The natural logarithm of Ku was than plotted against 1/T, for all values corresponing to the protein folding between 10-50% with T, being the absolute temperatur. From the linear fit ∆Hunfolding and ∆Sunfolding can be derived. ∆Hunfolding, being the slope and ∆Sunfolding being the y intercept. ∆Gunfolding is then calculated using equation 5.
Δ G u n f o l d i n g = Δ H u n f o l d i n g − T Δ S u n f o l d i n g

4.7. Asymmetrical flow field-flow fractionation (AF4)

AF4 was used for analysis of the aggregation behaviour of PCNAWT and PCNAC148S. Thermal stress was induced by incubating the proteins at 37 °C for 1 or 16 h or at 65 °C for 15 min. Afterwards, the samples were cooled to 4 °C until they were injected into AF4. The AF2000 Flow FFF system (Postnova Analytics, Landsberg, Germany) coupled with an UV/Vis detector as well as a fluorescence detector (Shimadzu, Kyōto, Japan) was used. The UV/Vis detector was set to a wavelength of 280 nm. The fluorescence detector was set to an excitation wavelength of 569 nm and emission wavelength of 610 nm with a gain of 128. The Bio-Rad size standard (MW 1.35 - 670 kDa, Bio-Rad Laboratories, Hercules, USA) served for the determination of the molecular weight of the analysed proteins. The separation channel was equipped with a regenerated cellulose membrane (MWCO 10 kDa) and a trapezoidal spacer with tapered ends (height: 350 µm). All analyses were performed at 25 °C in phosphate buffer pH 6.6 (1.675 g/L NaH2PO4 x 1 H2O, 0.875 g/L Na2HPO4 x 2 H2O, 8.765 g/L NaCl) as mobile phase and samples were eluted at a detector flow rate of 0.5 mL/min. The samples (injection volume of 20 µL) were injected into the separation channel during focusing step (4 min) at an injection flow of 0.2 mL/min, a cross-flow of 2.5 mL/min and a focus-flow of 2.8 mL/min. After a transition time of 0.2 min, the elution step was started: the cross-flow was kept constant (2.5 mL/min) for 10 min and then decreased linearly over 10 min to 0 mL/min. To ensure complete elution of the samples, the elution step was continued for 40 min with a detector flow of 0.5 mL/min (cross flow 0 mL/min). This was followed by a rinsing step with a detector flow of 0.2 mL/min (4 min). The data were evaluated using GraphPad Prism 5 (GraphPad Software, La Jolla, USA) and a baseline correction was performed for a better data visualisation.

4.8. Intrinsic tryptophan fluorescence

The intrinsic tryptophan fluorescence measurements were carried out in PBS containing 20 µL of 0.2 mg/mL PCNAWT or PCNAC148S solution in a 384 well microplate (788876, Greiner Bio-One, Kremsmünster, Austria). Chemical denaturation was induced by increasing urea concentrations (0-8 M) and a subsequent incubation at RT for 5 h. The intrinsic tryptophan fluorescence was excited at 284 nm and the emission were detected at 320 nm with a bandwidth of 9 nm at RT by using the plate reader Synergy MX (BioTek Instruments GmbH, Bad Friedrichshall, Germany). The normalized fluorescence intensities of three measurements in duplicates were plotted against the urea concentration. The Cm-values were determined by using equation 6 and OriginPro 2023 (OriginLab Corporation, Northampton, USA), where A2 is the final value, A1 is the initial value, x0 is the center and p is the power.
y = A 2 + ( A 1 − A 2 ) / ( 1 + ( x / x 0 ) ^ p )
Due to the biphasic denaturation profile of PCNAWT, the graph was split into two graphs for the determination of the Cm-values. The given deviation of the Cm values is the standard error of the fit.

4.9. Alchemical free energy predictions

System preparation

The crystal structure of PCNA in complex with p15 (PDB ID: 6GWS) was retrieved from the PDB server. Only one PCNA monomer in complex with p15 (chain A and chain D respectively) were kept for all further calculations. The remaining structure was checked for missing atoms using the Structure Preparation panel from MOE (Molecular Operating Environment, Chemical Computing Group, Canada version: 2022.02). Protonation states were assigned with the Protonate 3D function using default settings at pH=7.0. Mutations were generated using the Protein Builder from MOE. Sidechains were repacked and minimized using the Minimize tool, while tethering the backbone. The structures were exported from MOE in PDB format and used for MD system preparation within GROMACS 2019.3. [78].

Free energy calculations and hybrid topologies

To estimate the relative free energy difference for the folding of PCNA or the binding of p15 to PCNA, we used a well-established thermodynamic cycle shown in Figure S4 [79]. The amino acids were morphed along the ΔG1 and ΔG2 arrows for the unfolded and folded, or the free PCNA and PCNA-p15 complex. The unfolded state for the estimation of the free folding energy difference, ΔΔGfolding, was modelled by a glycine-x-glycine tripeptide, as described earlier [45,80]. All simulations were conducted separately. Hybrid topologies were generated using the pmx package after running the equilibrium simulations [80]. This way the equilibrium trajectories can be used to extract coordinates for different transitions later on.

MD Simulations

All MD simulations were carried out using GROMACS 2019.3. Proteins were modelled with the Amber99SB*ILDN forcefield [81], while water molecules were represented by the TIP3P water model [81]. Proteins, or tripeptides were solvated in a cubic box with 150 mM of Na+ and Cl- ions. The total charge of the system was neutralized with either Na+ or Cl-. The system was energy minimized using the steepest decent algorithm, then equilibrate for 500 ps in the NVT ensemble and then for 500 ps in the NPT ensemble. Production simulations were run for 100 ns in the NPT ensemble. During simulation H-Bonds were constrained using LINCS algorithm [82]. Particle Mesh Ewald was used to treat Electrostatic interactions in the simulation [83]. The cutoff for electrostatic interactions was set to 1.2 nm. Simulation temperature was controlled by the velocity rescaling thermostat at 300 K with a time constant of 0.1 ps [84]. Pressure coupling was done, using the Parinello-Rahman barostat at 1 atm [85].
From these equilibrium trajectories, snapshots were extracted to construct hybrid topologies and perform the fast-alchemical switching reactions. The first 5 ns of every trajectory were discarded and the rest was used to generate 95 snapshots (one every 1000 ps). Hybrid structures and topologies were than generated using the pmx package for each frame [80]. After energy minimization and 50 ps of equilibration, alchemical transitions were performed in 100 ps. This resulted in a change in λ value of 5*10-5/step. A softcore potential was used for the electrostatic interaction and the Lennard-Jones interactions during the transition, with default parameters. The work values for each transition were extracted and used to estimate the free energy at 300 K using the scripts provided with the pmx package. The free energy estimation was based on the Crooks Fluctuation Theorem [47]. Bennett acceptance ratio was used as a maximum likelihood estimator [86]. The uncertainty given is the standard deviation of the free energy estimation over three independent replicas.

4.10. Glutaraldehyde crosslinking of PCNA

46 µM PCNAWT or PCNAC148S were incubated with 0.004% glutaraldehyde in PBS for 10 min at RT. The reaction was stopped by the addition of 0.1 M glycine. For analysis of the crosslinking, 5 µg protein were loaded onto the SDS-PAGE gel (10%), followed by a Coomassie brilliant- blue staining or a western blot analysis.

4.11. Western Blot analysis

The same volume of SDS sample buffer (100 mM Tris-HCl, 200 mM dithiothreitol, 4 % (w/v) SDS, 0.2 % (w/v) bromophenol blue, 20 % (v/v) glycerol, pH 6.8) was added to the protein samples and heated for 10 min at 95 °C. The proteins were separated by SDS PAGE and transferred to a polyvinylidene fluoride membrane using semidry electroblotting (Trans-Blot SD, Bio-Rad, Hercules, USA). Afterwards, the membrane was blocked with 5% BSA in PBS over night at 4 °C under shaking. The immunodetection was conducted using first a His6 antibody (1:1000 in PBS with 0.01% NaN3, 1 h at RT, mouse anti-6x-His, Thermo Fisher Scientific, Waltham, USA), followed by a horseradish peroxidase (HRP) coupled antibody (1:6000 in PBS with 0.1% Tween20, 1.5 h at RT, rabbit anti-mouse, Thermo Fisher Scientific, Waltham, USA). The washed membrane was treated with Immuno Cruz Western Blot Luminol Reagent (Santa Cruz Biotechnology, Dallas, USA) and the luminicence was detected by chemiluminescence reader (ChemoCam ECL imager, Intas, Göttingen, Germany).

Author Contributions

Conceptualization, J.J. and S.H.; methodology, S.H., S.S., A.A., K.L.; validation, S.H., S.S. and K.L.; formal analysis, S.H., S.S. and K.L.; investigation, S.H., S.S.; writing—original draft preparation, S.H., S.S.; writing—review and editing, J.J.; supervision, J.J.; project administration, J.J.; funding acquisition, J.J.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available in the manuscript and in the supplementary material as provided.

Acknowledgments

We thank Anke Ostmann for technical assistance. We acknowledge the support of HPC facility of the University of Münster for providing computational resources.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Indiani, C.; O’Donnell, M. The replication clamp-loading machine at work in the three domains of life. Nature Reviews Molecular Cell Biology 2006, 7, 751–761. [Google Scholar] [CrossRef]
  2. Kelman, Z.; O’Donnell, M. Structural and functional similarities of prokaryotic and eukaryotic DNA polymerase sliding clamps. Nucleic Acids Res 1995, 23, 3613–3620. [Google Scholar] [CrossRef]
  3. Krishna, T.S.; Kong, X.P.; Gary, S.; Burgers, P.M.; Kuriyan, J. Crystal structure of the eukaryotic DNA polymerase processivity factor PCNA. Cell 1994, 79, 1233–1243. [Google Scholar] [CrossRef]
  4. Matsumiya, S.; Ishino, Y.; Morikawa, K. Crystal structure of an archaeal DNA sliding clamp: proliferating cell nuclear antigen from Pyrococcus furiosus. Protein Sci 2001, 10, 17–23. [Google Scholar] [CrossRef]
  5. Williams, G.J.; Johnson, K.; Rudolf, J.; McMahon, S.A.; Carter, L.; Oke, M.; Liu, H.; Taylor, G.L.; White, M.F.; Naismith, J.H. Structure of the heterotrimeric PCNA from Sulfolobus solfataricus. Acta Crystallogr Sect F Struct Biol Cryst Commun 2006, 62, 944–948. [Google Scholar] [CrossRef]
  6. Mulye, M.; Singh, M.I.; Jain, V. From Processivity to Genome Maintenance: The Many Roles of Sliding Clamps. Genes (Basel) 2022, 13. [Google Scholar] [CrossRef]
  7. Stukenberg, P.T.; Studwell-Vaughan, P.S.; O’Donnell, M. Mechanism of the sliding beta-clamp of DNA polymerase III holoenzyme. The Journal of biological chemistry 1991, 266, 11328–11334. [Google Scholar] [CrossRef]
  8. Gulbis, J.M.; Kelman, Z.; Hurwitz, J.; O’Donnell, M.; Kuriyan, J. Structure of the C-terminal region of p21(WAF1/CIP1) complexed with human PCNA. Cell 1996, 87, 297–306. [Google Scholar] [CrossRef]
  9. Moldovan, G.L.; Pfander, B.; Jentsch, S. PCNA, the maestro of the replication fork. Cell 2007, 129, 665–679. [Google Scholar] [CrossRef]
  10. Choe, K.N.; Moldovan, G.-L. Forging ahead through darkness: PCNA, still the principal conductor at the replication fork. Mol Cell 2017, 65, 380–392. [Google Scholar] [CrossRef]
  11. Havens, C.G.; Walter, J.C. Docking of a specialized PIP Box onto chromatin-bound PCNA creates a degron for the ubiquitin ligase CRL4Cdt2. Mol Cell 2009, 35, 93–104. [Google Scholar] [CrossRef] [PubMed]
  12. Prestel, A.; Wichmann, N.; Martins, J.M.; Marabini, R.; Kassem, N.; Broendum, S.S.; Otterlei, M.; Nielsen, O.; Willemoës, M.; Ploug, M.; et al. The PCNA interaction motifs revisited: thinking outside the PIP-box. Cellular and Molecular Life Sciences 2019, 76, 4923–4943. [Google Scholar] [CrossRef]
  13. Warbrick, E. PCNA binding through a conserved motif. Bioessays 1998, 20, 195–199. [Google Scholar] [CrossRef]
  14. Xu, H.; Zhang, P.; Liu, L.; Lee, M.Y. A novel PCNA-binding motif identified by the panning of a random peptide display library. Biochemistry 2001, 40, 4512–4520. [Google Scholar] [CrossRef]
  15. Gilljam, K.M.; Feyzi, E.; Aas, P.A.; Sousa, M.M.; Müller, R.; Vågbø, C.B.; Catterall, T.C.; Liabakk, N.B.; Slupphaug, G.; Drabløs, F.; et al. Identification of a novel, widespread, and functionally important PCNA-binding motif. J Cell Biol 2009, 186, 645–654. [Google Scholar] [CrossRef]
  16. Hara, K.; Uchida, M.; Tagata, R.; Yokoyama, H.; Ishikawa, Y.; Hishiki, A.; Hashimoto, H. Structure of proliferating cell nuclear antigen (PCNA) bound to an APIM peptide reveals the universality of PCNA interaction. Acta Crystallogr F Struct Biol Commun 2018, 74, 214–221. [Google Scholar] [CrossRef]
  17. Stoimenov, I.; Helleday, T. PCNA on the crossroad of cancer. Biochemical Society Transactions 2009, 37, 605–613. [Google Scholar] [CrossRef]
  18. Cardano, M.; Tribioli, C.; Prosperi, E. Targeting Proliferating Cell Nuclear Antigen (PCNA) as an Effective Strategy to Inhibit Tumor Cell Proliferation. Curr Cancer Drug Targets 2020, 20, 240–252. [Google Scholar] [CrossRef] [PubMed]
  19. Wang, S.C. PCNA: a silent housekeeper or a potential therapeutic target? Trends Pharmacol Sci 2014, 35, 178–186. [Google Scholar] [CrossRef]
  20. Bose, S.; Ingham, M.; Singh-Kandah, S.V.; Magana, W.; Schwartz, G.K. A phase II study, with a safety lead-in, to evaluate ATX-101, a peptide drug targeting PCNA, in advanced dedifferentiated liposarcoma and leiomyosarcoma. Journal of Clinical Oncology 2022, 40, TPS11587–TPS11587. [Google Scholar] [CrossRef]
  21. Müller, R.; Misund, K.; Holien, T.; Bachke, S.; Gilljam, K.M.; Våtsveen, T.K.; Rø, T.B.; Bellacchio, E.; Sundan, A.; Otterlei, M. Targeting proliferating cell nuclear antigen and its protein interactions induces apoptosis in multiple myeloma cells. PloS one 2013, 8, e70430. [Google Scholar] [CrossRef] [PubMed]
  22. Lemech, C.R.; Kichenadasse, G.; Marschner, J.-P.; Alevizopoulos, K.; Otterlei, M.; Millward, M. ATX-101, a cell-penetrating protein targeting PCNA, can be safely administered as intravenous infusion in patients and shows clinical activity in a Phase 1 study. Oncogene 2023, 42, 541–544. [Google Scholar] [CrossRef] [PubMed]
  23. Inoue, A.; Kikuchi, S.; Hishiki, A.; Shao, Y.; Heath, R.; Evison, B.J.; Actis, M.; Canman, C.E.; Hashimoto, H.; Fujii, N. A small molecule inhibitor of monoubiquitinated Proliferating Cell Nuclear Antigen (PCNA) inhibits repair of interstrand DNA cross-link, enhances DNA double strand break, and sensitizes cancer cells to cisplatin. The Journal of biological chemistry 2014, 289, 7109–7120. [Google Scholar] [CrossRef]
  24. Punchihewa, C.; Inoue, A.; Hishiki, A.; Fujikawa, Y.; Connelly, M.; Evison, B.; Shao, Y.; Heath, R.; Kuraoka, I.; Rodrigues, P.; et al. Identification of small molecule proliferating cell nuclear antigen (PCNA) inhibitor that disrupts interactions with PIP-box proteins and inhibits DNA replication. The Journal of biological chemistry 2012, 287, 14289–14300. [Google Scholar] [CrossRef]
  25. Evison, B.J.; Actis, M.L.; Wu, S.Z.; Shao, Y.; Heath, R.J.; Yang, L.; Fujii, N. A site-selective, irreversible inhibitor of the DNA replication auxiliary factor proliferating cell nuclear antigen (PCNA). Bioorg Med Chem 2014, 22, 6333–6343. [Google Scholar] [CrossRef]
  26. De Biasio, A.; de Opakua, A.I.; Mortuza, G.B.; Molina, R.; Cordeiro, T.N.; Castillo, F.; Villate, M.; Merino, N.; Delgado, S.; Gil-Carton, D.; et al. Structure of p15(PAF)-PCNA complex and implications for clamp sliding during DNA replication and repair. Nat Commun 2015, 6, 6439. [Google Scholar] [CrossRef] [PubMed]
  27. González-Magaña, A.; Blanco, F.J. Human PCNA Structure, Function and Interactions. Biomolecules 2020, 10. [Google Scholar] [CrossRef] [PubMed]
  28. Xie, C.; Yao, M.; Dong, Q. Proliferating cell unclear antigen-associated factor (PAF15): a novel oncogene. Int J Biochem Cell Biol 2014, 50, 127–131. [Google Scholar] [CrossRef]
  29. Bao, S.; Wu, Q.; McLendon, R.E.; Hao, Y.; Shi, Q.; Hjelmeland, A.B.; Dewhirst, M.W.; Bigner, D.D.; Rich, J.N. Glioma stem cells promote radioresistance by preferential activation of the DNA damage response. Nature 2006, 444, 756–760. [Google Scholar] [CrossRef]
  30. Ong, D.S.T.; Hu, B.; Ho, Y.W.; Sauvé, C.-E.G.; Bristow, C.A.; Wang, Q.; Multani, A.S.; Chen, P.; Nezi, L.; Jiang, S.; et al. PAF promotes stemness and radioresistance of glioma stem cells. Proceedings of the National Academy of Sciences of the United States of America 2017, 114, E9086–E9095. [Google Scholar] [CrossRef]
  31. Yan, R.; Zhu, K.; Dang, C.; Lan, K.; Wang, H.; Yuan, D.; Chen, W.; Meltzer, S.J.; Li, K. Paf15 expression correlates with rectal cancer prognosis, cell proliferation and radiation response. Oncotarget 2016, 7, 38750–38761. [Google Scholar] [CrossRef] [PubMed]
  32. Roa, S.; Avdievich, E.; Peled, J.U.; Maccarthy, T.; Werling, U.; Kuang, F.L.; Kan, R.; Zhao, C.; Bergman, A.; Cohen, P.E.; et al. Ubiquitylated PCNA plays a role in somatic hypermutation and class-switch recombination and is required for meiotic progression. Proceedings of the National Academy of Sciences of the United States of America 2008, 105, 16248–16253. [Google Scholar] [CrossRef]
  33. Baple, E.L.; Chambers, H.; Cross, H.E.; Fawcett, H.; Nakazawa, Y.; Chioza, B.A.; Harlalka, G.V.; Mansour, S.; Sreekantan-Nair, A.; Patton, M.A.; et al. Hypomorphic PCNA mutation underlies a human DNA repair disorder. J Clin Invest 2014, 124, 3137–3146. [Google Scholar] [CrossRef] [PubMed]
  34. Raslan, I.R.; de Assis Pereira Matos, P.C.A.; Boaratti Ciarlariello, V.; Daghastanli, K.H.; Rosa, A.B.R.; Arita, J.H.; Aranda, C.S.; Barsottini, O.G.P.; Pedroso, J.L. Beyond Typical Ataxia Telangiectasia: How to Identify the Ataxia Telangiectasia-Like Disorders. Mov Disord Clin Pract 2021, 8, 118–125. [Google Scholar] [CrossRef]
  35. Duffy, C.M.; Hilbert, B.J.; Kelch, B.A. A Disease-Causing Variant in PCNA Disrupts a Promiscuous Protein Binding Site. J Mol Biol 2016, 428, 1023–1040. [Google Scholar] [CrossRef] [PubMed]
  36. Magrino, J.; Munford, V.; Martins, D.J.; Homma, T.K.; Page, B.; Gaubitz, C.; Freire, B.L.; Lerario, A.M.; Vilar, J.B.; Amorin, A.; et al. A thermosensitive PCNA allele underlies an ataxia-telangiectasia-like disorder. The Journal of biological chemistry 2023, 299, 104656. [Google Scholar] [CrossRef]
  37. Bindels, D.S.; Haarbosch, L.; van Weeren, L.; Postma, M.; Wiese, K.E.; Mastop, M.; Aumonier, S.; Gotthard, G.; Royant, A.; Hink, M.A.; et al. mScarlet: a bright monomeric red fluorescent protein for cellular imaging. Nature Methods 2017, 14, 53–56. [Google Scholar] [CrossRef]
  38. Shaner, N.C.; Lambert, G.G.; Chammas, A.; Ni, Y.; Cranfill, P.J.; Baird, M.A.; Sell, B.R.; Allen, J.R.; Day, R.N.; Israelsson, M.; et al. A bright monomeric green fluorescent protein derived from Branchiostoma lanceolatum. Nature Methods 2013, 10, 407–409. [Google Scholar] [CrossRef]
  39. Song, Y.; Madahar, V.; Liao, J. Development of FRET assay into quantitative and high-throughput screening technology platforms for protein-protein interactions. Ann Biomed Eng 2011, 39, 1224–1234. [Google Scholar] [CrossRef]
  40. Liao, J.-y.; Song, Y.; Liu, Y. A new trend to determine biochemical parameters by quantitative FRET assays. Acta Pharmacologica Sinica 2015, 36, 1408–1415. [Google Scholar] [CrossRef]
  41. Lazareno, S. Quantification of receptor interactions using binding methods. J Recept Signal Transduct Res 2001, 21, 139–165. [Google Scholar] [CrossRef]
  42. Cheng, Y.; Prusoff, W.H. Relationship between the inhibition constant (K1) and the concentration of inhibitor which causes 50 per cent inhibition (I50) of an enzymatic reaction. Biochemical pharmacology 1973, 22, 3099–3108. [Google Scholar] [CrossRef] [PubMed]
  43. De March, M.; Merino, N.; Barrera-Vilarmau, S.; Crehuet, R.; Onesti, S.; Blanco, F.J.; De Biasio, A. Structural basis of human PCNA sliding on DNA. Nature Communications 2017, 8, 13935. [Google Scholar] [CrossRef] [PubMed]
  44. Motulsky, H.J.; Neubig, R.R. Analyzing binding data. Curr Protoc Neurosci, 2010, Chapter 7, Unit 7.5,. [CrossRef]
  45. Gapsys, V.; Michielssens, S.; Seeliger, D.; de Groot, B.L. Accurate and Rigorous Prediction of the Changes in Protein Free Energies in a Large-Scale Mutation Scan. Angewandte Chemie (International ed. in English) 2016, 55, 7364–7368. [Google Scholar] [CrossRef] [PubMed]
  46. Patel, D.; Patel, J.S.; Ytreberg, F.M. Implementing and Assessing an Alchemical Method for Calculating Protein-Protein Binding Free Energy. Journal of chemical theory and computation 2021, 17, 2457–2464. [Google Scholar] [CrossRef]
  47. Crooks, G.E. Entropy production fluctuation theorem and the nonequilibrium work relation for free energy differences. Physical Review E 1999, 60, 2721–2726. [Google Scholar] [CrossRef]
  48. Jarzynski, C. Nonequilibrium Equality for Free Energy Differences. Physical Review Letters 1997, 78, 2690–2693. [Google Scholar] [CrossRef]
  49. Jarzynski, C. Equilibrium free-energy differences from nonequilibrium measurements: A master-equation approach. Physical Review E 1997, 56, 5018–5035. [Google Scholar] [CrossRef]
  50. Seeliger, D.; de Groot, B.L. Protein Thermostability Calculations Using Alchemical Free Energy Simulations. Biophysical journal 2010, 98, 2309–2316. [Google Scholar] [CrossRef]
  51. Huang, P.; Chu, S.K.S.; Frizzo, H.N.; Connolly, M.P.; Caster, R.W.; Siegel, J.B. Evaluating Protein Engineering Thermostability Prediction Tools Using an Independently Generated Dataset. ACS Omega 2020, 5, 6487–6493. [Google Scholar] [CrossRef]
  52. Wang, W. Protein aggregation and its inhibition in biopharmaceutics. International journal of pharmaceutics 2005, 289, 1–30. [Google Scholar] [CrossRef]
  53. Royer, C.A. Probing protein folding and conformational transitions with fluorescence. Chem Rev 2006, 106, 1769–1784. [Google Scholar] [CrossRef] [PubMed]
  54. Gao, K.; Oerlemans, R.; Groves, M.R. Theory and applications of differential scanning fluorimetry in early-stage drug discovery. Biophys Rev 2020, 12, 85–104. [Google Scholar] [CrossRef] [PubMed]
  55. Leeman, M.; Nilsson, L.; Storm, M.U. Practical applications of asymmetrical flow field-flow fractionation (AF4): a review. 2015.
  56. Mahler, H.C.; Friess, W.; Grauschopf, U.; Kiese, S. Protein aggregation: pathways, induction factors and analysis. J Pharm Sci 2009, 98, 2909–2934. [Google Scholar] [CrossRef] [PubMed]
  57. Jarmoskaite, I.; AlSadhan, I.; Vaidyanathan, P.P.; Herschlag, D. How to measure and evaluate binding affinities. Elife 2020, 9. [Google Scholar] [CrossRef]
  58. Sowa, S.T.; Vela-Rodríguez, C.; Galera-Prat, A.; Cázares-Olivera, M.; Prunskaite-Hyyryläinen, R.; Ignatev, A.; Lehtiö, L. A FRET-based high-throughput screening platform for the discovery of chemical probes targeting the scaffolding functions of human tankyrases. Scientific Reports 2020, 10, 12357. [Google Scholar] [CrossRef]
  59. Bartolowits, M.D.; Gast, J.M.; Hasler, A.J.; Cirrincione, A.M.; O’Connor, R.J.; Mahmoud, A.H.; Lill, M.A.; Davisson, V.J. Discovery of Inhibitors for Proliferating Cell Nuclear Antigen Using a Computational-Based Linked-Multiple-Fragment Screen. ACS Omega 2019, 4, 15181–15196. [Google Scholar] [CrossRef]
  60. Wegerer, A.; Sun, T.; Altenbuchner, J. Optimization of an E. coli L-rhamnose-inducible expression vector: test of various genetic module combinations. BMC Biotechnology 2008, 8, 2–2. [Google Scholar] [CrossRef]
  61. Jahn, T.R.; Radford, S.E. Folding versus aggregation: polypeptide conformations on competing pathways. Arch Biochem Biophys 2008, 469, 100–117. [Google Scholar] [CrossRef]
  62. Kontopidis, G.; Wu, S.-Y.; Zheleva, D.I.; Taylor, P.; McInnes, C.; Lane, D.P.; Fischer, P.M.; Walkinshaw, M.D. Structural and biochemical studies of human proliferating cell nuclear antigen complexes provide a rationale for cyclin association and inhibitor design. Proceedings of the National Academy of Sciences of the United States of America 2005, 102, 1871–1876. [Google Scholar] [CrossRef]
  63. González-Mondragón, E.; Zubillaga, R.A.; Saavedra, E.; Chánez-Cárdenas, M.E.; Pérez-Montfort, R.; Hernández-Arana, A. Conserved Cysteine 126 in Triosephosphate Isomerase Is Required Not for Enzymatic Activity but for Proper Folding and Stability. Biochemistry 2004, 43, 3255–3263. [Google Scholar] [CrossRef]
  64. Pavlin, M.; Qasem, Z.; Sameach, H.; Gevorkyan-Airapetov, L.; Ritacco, I.; Ruthstein, S.; Magistrato, A. Unraveling the Impact of Cysteine-to-Serine Mutations on the Structural and Functional Properties of Cu(I)-Binding Proteins. Int J Mol Sci 2019, 20. [Google Scholar] [CrossRef] [PubMed]
  65. Tsai, Y.-C.; Wang, Y.-H.; Liu, Y.-C. Overexpression of PCNA Attenuates Oxidative Stress-Caused Delay of Gap-Filling during Repair of UV-Induced DNA Damage. Journal of Nucleic Acids 2017, 2017, 8154646. [Google Scholar] [CrossRef] [PubMed]
  66. Jaskulski, D.; deRiel, J.K.; Mercer, W.E.; Calabretta, B.; Baserga, R. Inhibition of cellular proliferation by antisense oligodeoxynucleotides to PCNA cyclin. Science 1988, 240, 1544–1546. [Google Scholar] [CrossRef] [PubMed]
  67. Groehler, A.L.; Lannigan, D.A. A chromatin-bound kinase, ERK8, protects genomic integrity by inhibiting HDM2-mediated degradation of the DNA clamp PCNA. J Cell Biol 2010, 190, 575–586. [Google Scholar] [CrossRef]
  68. Sporbert, A.; Domaing, P.; Leonhardt, H.; Cardoso, M.C. PCNA acts as a stationary loading platform for transiently interacting Okazaki fragment maturation proteins. Nucleic Acids Res 2005, 33, 3521–3528. [Google Scholar] [CrossRef]
  69. Kurki, P.; Vanderlaan, M.; Dolbeare, F.; Gray, J.; Tan, E.M. Expression of proliferating cell nuclear antigen (PCNA)/cyclin during the cell cycle. Exp Cell Res 1986, 166, 209–219. [Google Scholar] [CrossRef]
  70. Green, C.M.; Baple, E.L.; Crosby, A.H. PCNA mutation affects DNA repair not replication. Cell Cycle 2014, 13, 3157–3158. [Google Scholar] [CrossRef]
  71. Invernizzi, G.; Papaleo, E.; Sabate, R.; Ventura, S. Protein aggregation: mechanisms and functional consequences. Int J Biochem Cell Biol 2012, 44, 1541–1554. [Google Scholar] [CrossRef]
  72. De Baets, G.; Van Doorn, L.; Rousseau, F.; Schymkowitz, J. Increased Aggregation Is More Frequently Associated to Human Disease-Associated Mutations Than to Neutral Polymorphisms. PLoS Comput Biol 2015, 11, e1004374. [Google Scholar] [CrossRef]
  73. Gil-Martínez, J.; Bernardo-Seisdedos, G.; Mato, J.M.; Millet, O. The use of pharmacological chaperones in rare diseases caused by reduced protein stability. Proteomics 2022, 22, e2200222. [Google Scholar] [CrossRef]
  74. Liguori, L.; Monticelli, M.; Allocca, M.; Hay Mele, B.; Lukas, J.; Cubellis, M.V.; Andreotti, G. Pharmacological Chaperones: A Therapeutic Approach for Diseases Caused by Destabilizing Missense Mutations. Int J Mol Sci 2020, 21. [Google Scholar] [CrossRef]
  75. Leonhardt, H.; Rahn, H.P.; Weinzierl, P.; Sporbert, A.; Cremer, T.; Zink, D.; Cardoso, M.C. Dynamics of DNA replication factories in living cells. J Cell Biol 2000, 149, 271–280. [Google Scholar] [CrossRef] [PubMed]
  76. Boivin, S.; Kozak, S.; Meijers, R. Optimization of protein purification and characterization using Thermofluor screens. Protein Expr Purif 2013, 91, 192–206. [Google Scholar] [CrossRef] [PubMed]
  77. Wright, T.A.; Stewart, J.M.; Page, R.C.; Konkolewicz, D. Extraction of Thermodynamic Parameters of Protein Unfolding Using Parallelized Differential Scanning Fluorimetry. The Journal of Physical Chemistry Letters 2017, 8, 553–558. [Google Scholar] [CrossRef] [PubMed]
  78. Abraham, M.J.; Murtola, T.; Schulz, R.; Páll, S.; Smith, J.C.; Hess, B.; Lindahl, E. GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX 2015, 1-2, 19–25. [Google Scholar] [CrossRef]
  79. Aldeghi, M.; de Groot, B.L.; Gapsys, V. Accurate Calculation of Free Energy Changes upon Amino Acid Mutation. Methods Mol Biol 2019, 1851, 19–47. [Google Scholar] [CrossRef] [PubMed]
  80. Gapsys, V.; Michielssens, S.; Seeliger, D.; de Groot, B.L. pmx: Automated protein structure and topology generation for alchemical perturbations. 2015; 36, 348–354. [Google Scholar] [CrossRef]
  81. Lindorff-Larsen, K.; Piana, S.; Palmo, K.; Maragakis, P.; Klepeis, J.L.; Dror, R.O.; Shaw, D.E. Improved side-chain torsion potentials for the Amber ff99SB protein force field. Proteins 2010, 78, 1950–1958. [Google Scholar] [CrossRef]
  82. Hess, B.; Bekker, H.; Berendsen, H.J.C.; Fraaije, J.G.E.M. LINCS: A linear constraint solver for molecular simulations. 1997; 18, 1463–1472. [Google Scholar] [CrossRef]
  83. Essmann, U.; Perera, L.; Berkowitz, M.L.; Darden, T.; Lee, H.; Pedersen, L.G. A smooth particle mesh Ewald method. 1995; 103, 8577–8593. [Google Scholar] [CrossRef]
  84. Bussi, G.; Donadio, D.; Parrinello, M. Canonical sampling through velocity rescaling. The Journal of chemical physics 2007, 126, 014101. [Google Scholar] [CrossRef]
  85. Parrinello, M.; Rahman, A. Polymorphic transitions in single crystals: A new molecular dynamics method. 1981; 52, 7182–7190. [Google Scholar] [CrossRef]
  86. Shirts, M.R.; Bair, E.; Hooker, G.; Pande, V.S. Equilibrium free energies from nonequilibrium measurements using maximum-likelihood methods. Phys Rev Lett 2003, 91, 140601. [Google Scholar] [CrossRef]
  87. HomoloGene [Internet]. Bethesda (MD): National Library of Medicine (US), National Center for Biotechnology Information; 2022 12 01. Available online: https://www.ncbi.nlm.nih.gov/homologene.
Preprints 78749 g001
Figure 1. Analysis of PCNA-p15 interaction by FRET. (A) Schematic representation of the FRET based assay for analysis of PCNA-p15 interaction. PCNA-p15 interaction brings the fluorescent proteins mNeonGreen (donor) and mScarlet-I (acceptor) in close proximity. This allows energy transfer from the donor molecule to the acceptor, resulting in a reduced donor emission (donor quenching) and an increased acceptor emission (sensitized emission). Binding is evaluated by analysing the sensitized emission using the three-channel method. (B) 1 µM p15-mNeonGreen (left) or mNeonGreen (right) were incubated with 0 - 2.125 µM mScarlet-I-PCNA for 1 h at 37 °C. The fluorescence was excited at 488 nm and the emission was detected at 517 – 700 nm. For p15-mNeonGreen/mScarlet-I, a quenched donor emission and an increased acceptor emission was observed.
Figure 1. Analysis of PCNA-p15 interaction by FRET. (A) Schematic representation of the FRET based assay for analysis of PCNA-p15 interaction. PCNA-p15 interaction brings the fluorescent proteins mNeonGreen (donor) and mScarlet-I (acceptor) in close proximity. This allows energy transfer from the donor molecule to the acceptor, resulting in a reduced donor emission (donor quenching) and an increased acceptor emission (sensitized emission). Binding is evaluated by analysing the sensitized emission using the three-channel method. (B) 1 µM p15-mNeonGreen (left) or mNeonGreen (right) were incubated with 0 - 2.125 µM mScarlet-I-PCNA for 1 h at 37 °C. The fluorescence was excited at 488 nm and the emission was detected at 517 – 700 nm. For p15-mNeonGreen/mScarlet-I, a quenched donor emission and an increased acceptor emission was observed.
Preprints 78749 g001
Preprints 78749 g002
Figure 2. Affinity analysis of PCNA-p15 interaction by FRET. 1 µM p15-mNeonGreen was incubated with 0 – 3.09 µM mScarlet-I-PCNA for 1 h at 37 °C. Binding is evaluated by analysing the sensitized emission using the three-channel method. The sensitised emission (SE) of three independent experiments in triplicates was plotted against the mScarlet-I-PCNA concentration, resulting in a KD value of 430.8 ± 59.9 nM for PCNA-p15 interaction. The given error is the standard error of the fit.
Figure 2. Affinity analysis of PCNA-p15 interaction by FRET. 1 µM p15-mNeonGreen was incubated with 0 – 3.09 µM mScarlet-I-PCNA for 1 h at 37 °C. Binding is evaluated by analysing the sensitized emission using the three-channel method. The sensitised emission (SE) of three independent experiments in triplicates was plotted against the mScarlet-I-PCNA concentration, resulting in a KD value of 430.8 ± 59.9 nM for PCNA-p15 interaction. The given error is the standard error of the fit.
Preprints 78749 g002
Preprints 78749 g003
Figure 3. Analysis of PCNA-p15 inhibition by unlabeled PCNA, peptide p1551-70 and small molecule inhibitor T2AA. 1 µM p15-mNeonGreen and 1 µM mScarlet-I-PCNA were incubated with varying inhibitor concentrations (0-100 µM for T2AA and p1551-70, 0-20 µM for unlabeled PCNA) for 1 h at 37 °C. Relative binding was plotted against PCNA, p1551-70, or T2AA concentrations, resulting in IC50 values of 1.24 ± 0.07, 9.81 ± 0.99, 13.81 ± 2.0 µM. The given error is the standard error of the fit. Ki values were determined using Cheng-Prusoff equation [42] and turned out to be at 0.37 ± 0.03 (unlabeled PCNA), 2.94 ± 0.31 (p1551-70) and 4.12 ± 0.61 µM (T2AA), respectively.
Figure 3. Analysis of PCNA-p15 inhibition by unlabeled PCNA, peptide p1551-70 and small molecule inhibitor T2AA. 1 µM p15-mNeonGreen and 1 µM mScarlet-I-PCNA were incubated with varying inhibitor concentrations (0-100 µM for T2AA and p1551-70, 0-20 µM for unlabeled PCNA) for 1 h at 37 °C. Relative binding was plotted against PCNA, p1551-70, or T2AA concentrations, resulting in IC50 values of 1.24 ± 0.07, 9.81 ± 0.99, 13.81 ± 2.0 µM. The given error is the standard error of the fit. Ki values were determined using Cheng-Prusoff equation [42] and turned out to be at 0.37 ± 0.03 (unlabeled PCNA), 2.94 ± 0.31 (p1551-70) and 4.12 ± 0.61 µM (T2AA), respectively.
Preprints 78749 g003
Preprints 78749 g004
Figure 4. Structure of human PCNA in complex with p1552-71 (PDB code 6EHT). Human PCNA is shown in gray with the position of the mutation at C148 highlighted in orange. p1552-71 is presented in blue and interacts with its C-terminal site with the hydrophobic pocket of PCNA. The N-terminal site binds to the inner site of the PCNA ring near the position of the mutation C148S. The position of the C148S mutation is emphasized in the window.
Figure 4. Structure of human PCNA in complex with p1552-71 (PDB code 6EHT). Human PCNA is shown in gray with the position of the mutation at C148 highlighted in orange. p1552-71 is presented in blue and interacts with its C-terminal site with the hydrophobic pocket of PCNA. The N-terminal site binds to the inner site of the PCNA ring near the position of the mutation C148S. The position of the C148S mutation is emphasized in the window.
Preprints 78749 g004
Preprints 78749 g006
Figure 6. Closed thermodynamic cycles of the ΔΔGfolding of PCNA (A) and ΔΔGaffinity of PCNA-p15 (B) upon amino acid mutation. ΔΔG values were obtained using the thermodynamic cycle shown in the methods section. The sum of all ΔΔG values within on cycle was calculated and resulted in values close to zero for both experiments, implicating good convergence of the simulations. The calculated ΔΔG for the mutation PCNAC148S suggest a destabilizing effect of the mutation on PCNA (ΔΔGC148S = 2.39 kcal/mol). The same mutation is predicted to have no effects on the binding free energy between PCNA and p15 (ΔΔGC148S = 0.21 kcal/mol). The uncertainty is the standard deviation of three independent replicas.
Figure 6. Closed thermodynamic cycles of the ΔΔGfolding of PCNA (A) and ΔΔGaffinity of PCNA-p15 (B) upon amino acid mutation. ΔΔG values were obtained using the thermodynamic cycle shown in the methods section. The sum of all ΔΔG values within on cycle was calculated and resulted in values close to zero for both experiments, implicating good convergence of the simulations. The calculated ΔΔG for the mutation PCNAC148S suggest a destabilizing effect of the mutation on PCNA (ΔΔGC148S = 2.39 kcal/mol). The same mutation is predicted to have no effects on the binding free energy between PCNA and p15 (ΔΔGC148S = 0.21 kcal/mol). The uncertainty is the standard deviation of three independent replicas.
Preprints 78749 g006
Preprints 78749 g007
Figure 7. Analysing the influence of C148S on the chemical stability of PCNA by intrinsic tryptophan fluorescence measurements. The urea induced denaturation of PCNAWT (black) and PCNAC148S (grey) were measured by analysing the intrinsic tryptophan fluorescence at 320 nm. (A) The fluorescence intensity of PCNA was plotted against varying urea concentrations (0-8 M) in an overlay of one representative experiment. (B) The normalized tryptophan fluorescence of PCNAC148S of three experiments was plotted against the urea concentration. PCNAC148S denatures at a Cm of 3.32 ± 0.15 M urea. (C) The normalized tryptophan fluorescence of PCNAWT of three experiments was plotted against the urea concentration. Due to the three-state denaturation of PCNAWT, the graph was split into two parts for the determination of the Cm values. A Cm1 at 1.56 ± 0.12 M and a Cm2 at 3.95 ± 0.18 M urea were ascertained. The standard deviation of Cm is the standard error of the fit.
Figure 7. Analysing the influence of C148S on the chemical stability of PCNA by intrinsic tryptophan fluorescence measurements. The urea induced denaturation of PCNAWT (black) and PCNAC148S (grey) were measured by analysing the intrinsic tryptophan fluorescence at 320 nm. (A) The fluorescence intensity of PCNA was plotted against varying urea concentrations (0-8 M) in an overlay of one representative experiment. (B) The normalized tryptophan fluorescence of PCNAC148S of three experiments was plotted against the urea concentration. PCNAC148S denatures at a Cm of 3.32 ± 0.15 M urea. (C) The normalized tryptophan fluorescence of PCNAWT of three experiments was plotted against the urea concentration. Due to the three-state denaturation of PCNAWT, the graph was split into two parts for the determination of the Cm values. A Cm1 at 1.56 ± 0.12 M and a Cm2 at 3.95 ± 0.18 M urea were ascertained. The standard deviation of Cm is the standard error of the fit.
Preprints 78749 g007
Preprints 78749 g008
Figure 8. Analysing the influence of C148S on the thermal stability of PCNA by DSF. (A) Representative melting curves of PCNAWT (black line), PCNAC148S (grey line), denatured PCNAWT (black dashed line) and denatured PCNAC148S (grey dashed line). (B) Normalized melting curves of PCNAWT and PCNAC148S indicate a reduced thermal stability of PCNAC148S. (C) First derivative of the melting curve of PCNAWT and PCNAC148S. For PCNAWT two melting points were determined, one at 49.5 ± 0.1 °C and one at 54.54 ± 0.66 °C. The calculated ΔG values are 5.65 ± 0.48 kcal/mol for PCNAWT and 4.89 ± 0.26 kcal/mol for PCNAC148S.
Figure 8. Analysing the influence of C148S on the thermal stability of PCNA by DSF. (A) Representative melting curves of PCNAWT (black line), PCNAC148S (grey line), denatured PCNAWT (black dashed line) and denatured PCNAC148S (grey dashed line). (B) Normalized melting curves of PCNAWT and PCNAC148S indicate a reduced thermal stability of PCNAC148S. (C) First derivative of the melting curve of PCNAWT and PCNAC148S. For PCNAWT two melting points were determined, one at 49.5 ± 0.1 °C and one at 54.54 ± 0.66 °C. The calculated ΔG values are 5.65 ± 0.48 kcal/mol for PCNAWT and 4.89 ± 0.26 kcal/mol for PCNAC148S.
Preprints 78749 g008
Preprints 78749 g009
Figure 9. Analysis of the aggregation behaviour of PCNA by AF4. Proteins were incubated at 37 °C for 0 h (black), 1 h (grey) or 16 h (blue) or at 65 °C for 10 min (orange) and were injected to AF4 afterwards. The induction of aggregates was evaluated by plotting the UV absorbance at 280 nm against the elution time. The increase of thermal stress correlates with elevated elution times, which can be explained by the formation of higher molecular weight complexes.
Figure 9. Analysis of the aggregation behaviour of PCNA by AF4. Proteins were incubated at 37 °C for 0 h (black), 1 h (grey) or 16 h (blue) or at 65 °C for 10 min (orange) and were injected to AF4 afterwards. The induction of aggregates was evaluated by plotting the UV absorbance at 280 nm against the elution time. The increase of thermal stress correlates with elevated elution times, which can be explained by the formation of higher molecular weight complexes.
Preprints 78749 g009
Preprints 78749 g010
Figure 10. Incubation time dependent analysis of PCNA-p15 affinity. The affinities of 1 µM p15 to PCNAWT or PCNAC148S (0 - 3.09 µM) were determined by FRET assay after 15, 30, 45, 60 min incubation at 37 °C. The KD-values of PCNAWT are constant at different incubation times. The KD-values of the PCNAC148S/p15 interaction increases during the incubation at 37 °C, indicating the stability caused decrease in the apparent binding.
Figure 10. Incubation time dependent analysis of PCNA-p15 affinity. The affinities of 1 µM p15 to PCNAWT or PCNAC148S (0 - 3.09 µM) were determined by FRET assay after 15, 30, 45, 60 min incubation at 37 °C. The KD-values of PCNAWT are constant at different incubation times. The KD-values of the PCNAC148S/p15 interaction increases during the incubation at 37 °C, indicating the stability caused decrease in the apparent binding.
Preprints 78749 g010
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated