Luca
Digiacomo‡
a,
Francesca
Giulimondi‡
a,
Anna Laura
Capriotti
b,
Susy
Piovesana
b,
Carmela Maria
Montone
b,
Riccardo Zenezini
Chiozzi
c,
Aldo
Laganà
b,
Morteza
Mahmoudi
d,
Daniela
Pozzi
a and
Giulio
Caracciolo
*a
aDepartment of Molecular Medicine, Sapienza University of Rome, Viale Regina Elena 291, 00161, Rome, Italy. E-mail: giulio.caracciolo@uniroma1.it
bDepartment of Chemistry, Sapienza University of Rome, P.le A. Moro 5, 00185, Rome, Italy
cBiomolecular Mass Spectrometry and Proteomics, Bijvoet Center for Biomolecular Research, Utrecht Institute for Pharmaceutical Sciences Utrecht University, Heidelberglaan 8, 3584 CS Utrecht, The Netherlands
dDepartment of Radiology, Precision Health Program, Michigan State University, MI, USA
First published on 17th May 2021
In the past few years, characterization of the protein corona (PC) that forms around liposomal systems has gained increasing interest for the development of novel therapeutic and diagnostic technologies. At the crossroads of fast-moving research fields, the interdisciplinarity of protein corona investigations poses challenges for experimental design and reporting. Isolation of liposome–protein complexes from biological fluids has been identified as a fundamental step of the entire workflow of PC characterization but exact specifications for conditions to optimize pelleting remain elusive. In the present work, key factors affecting precipitation of liposome–protein complexes by centrifugation, including time of centrifugation, total sample volume, lipid:protein ratio and contamination from biological NPs were comprehensively evaluated. Here we show that the total amount of isolated liposome–protein complexes and the extent of contamination from biological NPs may vary with influence factors. Our results provide protein corona researchers with precise indications to separate liposome–protein complexes from protein-rich fluids and include proper controls, thus they are anticipated to catalyze improved consistency of data mining and computational modelling of protein corona composition.
Over the last decade we have started to understand that several overlooked factors associated with the liposome–protein corona have long limited translation of liposome technology from benchtop to bedside.21,23,24 Although many relevant aspects of the liposome–protein corona have been clarified,25–28 there are still concerns to be addressed. Among them, separation of liposome–protein complexes from biological fluids has been identified as a fundamental step of the entire workflow of PC characterization. Commonly used isolation methods including size-exclusion chromatography (SEC), membrane ultrafiltration and centrifugation have been recently reviewed.29–31 While the first two methods have been shown to be prone to protein contamination,32 proper evaluation of factors affecting precipitation of liposome–protein complexes by centrifugation remains elusive. Experimental approaches based on centrifugation to study the PC of inorganic NPs were the first to be established.6 These methods were subsequently extended to liposomes and other lipid-based systems. However, even the field of inorganic NPs has yet to adopt a best practice approach mainly due to the enormous variability in nanomaterial characteristics.33 In particular, variations in material density make it impossible to provide an exact specification for centrifugation factors (e.g. time and speed) to optimize the pelleting of the nanomaterials with their associated proteins. Likewise, a systematic study of the factors influencing protein corona isolation from liposomes has not been reported so far thus impairing reproducibility and consistency among laboratories. To fulfil this gap, in the present work, via rational experimental design, key factors affecting isolation of liposome–protein complexes by centrifugation were comprehensively evaluated. These included time of incubation, total sample volume and lipid:protein ratio. The lipid amount in the supernatants was quantified by fluorescence measurements, while the protein composition of the coronas was determined by one-dimensional (1D) gel electrophoresis and liquid-chromatography tandem mass spectrometry (LC-MS/MS) experiments. In an attempt to generalize results, we extended our investigation to a wide range of lipid systems used in drug and gene delivery, i.e. un-PEGylated and PEGylated liposomes and cationic liposome/DNA complexes (lipoplexes). Our results indicate that the total amount of isolated proteins and the contamination from biological NPs may change with the influence factors, but the corona composition is only dependent on the liposomes' synthetic identity and the protein source. Reviewing the literature and the concerns associated with other isolation methods, it is clear that centrifugation should be the technique of choice for studies on the liposome–protein corona.
Protein corona isolation was performed by centrifuging all the samples for 15 minutes or 60 minutes at 4 °C, 21400g. This centrifugation speed corresponds to the highest speed of many benchtop centrifuge models and has been used in numerous previous investigations.36 After centrifugation, supernatants were collected and subjected to fluorescence measurements to quantify the residual amount of non-pelleted lipids.
As measured fluorescence depends on the sample volume, three distinct calibration curves for each total sample volume (i.e. 50 μL, 100 μL and 200 μL) were obtained (Fig. 1a). Recovery rate of lipids in the supernatant was represented by percentage ratio of the amount of lipid after precipitation to its initial amount (before precipitation) (Fig. 1b).
The lowest recovery rate of lipid in the supernatant (i.e. the highest efficacy of pelleting) was achieved when liposomes were exposed to low protein concentration (5% HP). This peculiar behavior is in agreement with the presence of large liposome–protein aggregates in solution at low HP concentration. In a previous publication,26 Giulimondi et al. have demonstrated that around 5% HP cationic liposome–protein complexes reach their isoelectric point. When neutrally charged complexes collide, short-range van der Waals forces dominate over electrostatic repulsions and induce formation of large-size aggregates. Such big clusters are more easily pelleted by centrifugation. Conversely, in solutions at 50% HP and 80% HP, the lipid amount in the supernatant ranged from 12% to 52% and followed general trends. As shown in Fig. 1b, increasing centrifugation time and decreasing sample volume generally led to smaller lipid recovery in the supernatant, i.e. to larger pelletizing of liposome–protein complexes. These findings clearly indicate that the amount of pelleted liposome–protein complexes is dependent on the technical parameters of the centrifugal procedure. Then, the question immediately arose whether the experimental protocol of centrifugation had an impact also on the corona composition. As relevance in characterization of protein corona has increased, the field has introduced a set of experimental methods to characterize it. This assortment of specializations poses concerns to investigation design and reporting.37 This is a point of great general interest as systematic changes in corona composition produced by the experimental protocol would make it impossible to compare results from different studies.
To address this issue, we studied the protein patterns of each sample by 1D SDS–PAGE. As a representative example, Fig. 2 summarizes the electrophoretic outcome for liposome–protein complexes obtained exposing liposomes to intermediate protein concentration (i.e. HP 50%). As shown in Fig. 2a (uncropped gel image is provided Fig. S1 in the ESI†), the lane intensity increased with total sample volume and centrifugation time and exhibited a very good correlation with the pelleted lipid amount evaluated by fluorescence measurements (Fig. 2b). This clearly indicates that the ratio of lipid to protein pelleted by centrifugation stays constant under different centrifugation conditions. The total lane intensities are also reported in the top panels of Fig. 2c, along with the corresponding intensity curves. An increasing trend of the intensity signal as function of the total sample volume can be recognized. However, the shape of the intensity curves seems to be the same for all the samples. This is confirmed by normalizing the curves, to obtain the molecular weight distributions of corona proteins (Fig. 2d). Interestingly, the observed increasing trend for the original intensity curves was found for the integral areas within different ranges of molecular weights between 20 kDa to 300 kDa (Fig. 2e). Conversely, no significant difference was found for the integral areas of the normalized patterns (Fig. 2f). Taken together, these results show that protein corona composition is not dependent on the centrifugation factors. This was confirmed for both low (i.e. 5% HP) and high (i.e. 80% HP) protein concentration (Fig. S1 in the ESI†). To our knowledge, this is the first demonstration that pelletizing liposome–protein complexes under different centrifuge conditions leads to the same protein corona composition.
Numerous previous investigations have clarified that protein binding to liposomes is controlled by several concomitant factors such as lipid composition,38 surface properties,39 incubation time,40 protein source41,42 and concentration.43 Among them, electrostatic attraction between proteins and lipids seems to be the dominant interaction. Thus, the effect of surface charge of liposomes on the efficacy of pelletizing was investigated next. Membrane charge density of lipid vesicles can be easily tuned by adjusting molar fraction of cationic lipids in the lipid bilayer.44 Fine-tuning molar fraction of DOTAP, we prepared (i) low charged liposomes (i.e. DOTAP–DOPE 1:3 (mol mol−1)) and (ii) highly charged liposomes (i.e. DOTAP–DOPE 3:1 (mol mol−1)).
As shown in Fig. S2,† the recovery rate of lipid in the supernatant was roughly the same for highly- and low charged liposomes. This result indicates that the efficacy of pelleting liposome–protein complexes is not affected by membrane charge density of liposomes. 1D SDS–PAGE experiments allowed us investigating the composition of the liposome–protein corona. Gel images corresponding to liposome–protein complexes for DOTAP–DOPE 1:1 (mol mol−1), DOTAP–DOPE 1:3 (mol mol−1) and DOTAP–DOPE 3:1 (mol mol−1) are shown in Fig. S1, S3 and S4† respectively. For each formulation, high correlation coefficients between the recovery rate of lipids in the pellet fluorescence and total lane intensity were measured (Fig. 3a–c). Average protein patterns for the three liposomal formulations upon incubation with human plasma at 5%, 50% and 80% HP are shown in Fig. 3d, e and f, respectively. We observe that protein patterns depended both on the employed formulation and protein concentration. This finding is in full agreement with previous results showing that the nanoparticle–protein corona is simultaneously shaped by (i) NPs' physiochemical properties and (ii) the protein source. However, for each combination of liposomal formulation and HP concentration, corona composition was not affected by the centrifugation protocol (Fig. S3 and S4†).
Fig. 3 1D SDS–PAGE analysis of the PC isolated from DOTAP–DOPE 1:3 (mol mol−1) (blue), DOTAP–DOPE 1:1 (mol mol−1) (green), DOTAP–DOPE 3:1 (mol mol−1) (yellow) under different centrifugation conditions (gel images are provided in Fig. S3, S1 and S4 in the ESI,† respectively). (a–c) Correlation analysis between total lane intensity from 1D SDS–PAGE and recovered lipid amount in pellets, measured by fluorescence experiments. (d–f) Average normalized 1D SDS–PAGE profiles grouped by HP concentration. |
As the protein corona affects manifold biological processes involved in the targeted delivery of nanomedicines, we included in our analysis two more relevant nanodelivery systems: PEGylated liposomes and PEGylated liposome/DNA complexes (lipoplexes). Grafting PEG to liposome surface represents one of the most used and effective surface functionalization to reduce protein binding, avoid opsonization and prolong bloodstream circulation of liposomes in vivo. On the other hand, PEGylated lipoplexes have been the object of numerous investigations for gene delivery purposes over the last two decades including cancer gene therapy.45
Fig. 4a depicts the recovery rates of lipid in pellets after incubation with HP followed by sample centrifugation, as a polar plot for unPEGylated DOTAP–DOPE liposomes, PEGylated DOTAP–DOPE liposomes and DOTAP–DOPE/DNA lipoplexes. Each radius of the plot represents one experimental condition, whereas the distance from the centre quantifies the relative lipid content in pellets to the original amount (i.e. pre-centrifuge), i.e. the bigger the distance from the centre, the higher the amount of pelleted liposomes with their associated proteins. The lipid content in pellets for the PEGylated liposome (violet) was much less abundant than that corresponding to the un-PEGylated formulation (green) as PEG reduced the adsorption of proteins. On the other hand, the pelleting of PEGylated lipoplexes (red) was intermediate with respect to the other two formulations. This suggests that plasma proteins bind specifically to plasmid DNA. This observation is in agreement with previous investigations that compared the protein profiles adsorbed to cationic liposomes and lipoplexes with the same lipid composition.46 Proteomics data showed that some complement proteins such as C5, C7 and C8 were exclusively found in the protein corona of lipoplexes. According to literature, this was attributed to the existence of a domain with DNA-binding activity in the complement proteins.47 1D SDS–PAGE experiments (gel images are provided in Fig. S1, S5 and S6†) confirmed that the protein corona is simultaneously shaped by physical–chemical properties of lipid vesicles and HP concentration (Fig. 4b). For small plasma concentration (5% HP), protein patterns are almost superimposed over the entire MW domain, but a peculiar peak within 20–30 kDa is exhibited by lipoplexes (red curve). A further peak at 30 kDa characterized its profile at 50% HP, whereas the main band for unPEGylated (green curve) and PEGylated (violet curve) liposomes was located at about 60 kDa and was much more prominent for the unPEGylated formulation. At 80% HP, that band represented the dominant peak for unPEGylated liposomes, whose pattern clearly differed from the other two, which in turn were found to be quite similar to each other. Second and foremost, electrophoretic results demonstrated that centrifugation only affected the total protein amount bound to vesicles but had no appreciable effect on the corona composition (Fig. S5 and S6 in the ESI†). As a next step, we performed nano liquid chromatography tandem mass spectrometry (nano LC-MS/MS) experiments to provide quantitative information on the protein corona composition (protein lists are reported in Table S1 in the ESI†). For each of the three formulations, we compared protein patterns isolated under the two experimental conditions leading to (i) the highest (i.e. total volume 50 μL and 60 min centrifugation time) and (ii) the lowest pelleting capacity (i.e. 200 μL total sample volume and 15 min centrifugation time). Results are represented as “Volcano plots” (Fig. 5), which aim to quantify the differences in relative protein abundances (RPA) of each single protein after corona isolation under the two explored experimental conditions. Proteins with similar RPA (i.e. p ≥ 0.05 or |log2(fold change)| ≤ 1) are indicated as grey dots and their cumulative RPA overcame 90% of the corona composition for all the studied formulations. Interestingly, this “overlap” percentage reached about 95% for PEGylated liposomes (Fig. 5b) and PEGylated lipoplexes (Fig. 5c). Corona proteins exhibiting relevant differences (i.e. p < 0.05 and |log2(fold change)| > 1) are indicated as colored dots, their cumulative RPA and total number are reported within the corresponding outer regions. Their contributions to the corona composition were found particularly small for PEGylated liposomes (Fig. 5b) and PEGylated lipoplexes (Fig. 5c), as they were always less than 4%. This value reached 8.4% for the unPEGylated formulation (Fig. 5a) and represented the highest portion of dissimilarity between the two experimental conditions of corona isolation. Of note, this portion is distributed over a remarkable number of protein (n = 59), thus the average RPA of the single proteins populating that region was quite small and read about 0.14%. In conclusion, nano LC-MS/MS experiments totally confirmed the findings of 1D SDS–PAGE demonstrating that the vast majority of the corona composition is not affected by the centrifugation protocol.
As recently pointed out experimental methods commonly used to isolate liposomes from biological fluids may be susceptible to protein contamination arising from unbound proteins and biological NPs.48 This could be caused by overlooked technical reasons and could help to rationalize data of different laboratories. It has been also postulated that different types of biological NPs present in human blood such as extracellular vesicles (EV) and lipoproteins need be properly considered in research field of bio–nano interactions. Fig. 6a shows a representative SDS–PAGE outcome of proteins isolated from liposome-free human plasma (indicated as “HP”) and liposomal samples (total volume 50 μL) isolated from DOTAP–DOPE 3:1 (mol mol−1), DOTAP–DOPE 1:1 (mol mol−1), DOTAP–DOPE 1:3 (mol mol−1), PEGylated DOTAP–DOPE 1:1 (mol mol−1) and PEGylated lipoplex, after 60 min centrifugation. Other two combinations of sample volumes and centrifugation times were studied, and results are reported in the ESI (Fig. S7 and S8†).
These three combinations of total sample volume and centrifugations time were chosen as they represent the conditions leading to highest (total volume 50 μL and 60 min centrifugation time), intermediate (total volume 50 μL and 15 min centrifugation time) and lowest (total volume 200 μL and 15 min centrifugation time) pelleting capacity. As a representative example of the comparison between the protein patterns of the investigated samples and the HP reference, we provide the corresponding absolute total lane intensity (Fig. 6b) and intensity profiles (Fig. 6c–g). The absolute profile of HP is superimposed to those of liposomal formulations as a black shaded area, in each panel. The gel image and the corresponding densitometric analysis indicate that the protein content of the “HP blank” sample (black shaded area) is negligible with respect to that truly associated to cationic unPEGylated liposomes and to PEGylated lipoplexes (Fig. 6c–e and g). These results are a direct evidence that the centrifugation-based isolation method resulted in minor, if any, contamination by unbound proteins and protein aggregates. However, the “HP blank” signal represented a small but considerable fraction of the protein patterns for non-loaded PEGylated liposomes (Fig. 6f). This is likely due to the lower recovery capacity by centrifugation for PEGylated liposomes than that achieved for the other systems (lipid recovery rates for PEGylated liposomes after centrifugation under all the studied conditions are reported in Fig. S9†). As a consequence, the HP baseline signal contributes more to the PC pattern. Adjusting the centrifugation protocol to enhance the recovery rates for PEGylated liposomes would be necessary to reduce the effects of plasma contamination and assess more accurately the PC composition. These trends are confirmed also for the other explored centrifugation conditions: total volume 50 μL, centrifugation time 15 min, and total volume 200 μL, centrifugation time 15 min (Fig. S7 and S8,† respectively).
The image processing procedure aims to evaluate the one-dimensional intensity distribution function of each sample, to obtain the corresponding one-dimensional molecular weight (MW) distribution. It is developed as follows. (i) On each of the gel images, background was preliminary removed to prevent unwanted baselines affect the resulting curves. (ii) The sampled intensity of the original image was converted into a new two-variable function where the vertical displacement was related to protein MWs (ladder lanes) and the intensity value was proportional to the detected protein amount. (iii) The projection of the two-variable function over a plane orthogonal to the image represented the absolute intensity profile of the generic lane. (iv) That intensity profile was finally normalized to 1 to obtain the corresponding protein pattern. Step (iii) and step (iv) were repeated for each lane of the studied gel. Further details can be found elsewhere.52
The acquired raw MS/MS data files from Xcalibur software (version 2.2 SP1.48, Thermo Fisher Scientific) were searched against Swiss-Prot human database by MaxQuant search engine, (v1.6.3.4)54 with the automatic setting for tryptic label free analysis. Fully tryptic digestion mode with maximum of two missed cleavages, multiplicity 1, minimum length of six amino acids and peptide tolerance of 20 ppm in first search and 4.5 ppm in the main search were set. For the modifications we set carbamidomethylation on cysteine as fixed, oxidation of methionine, acetylation of protein N-termini, and deamidation of glutamine and asparagine as variable. Recalibration and match between runs (with automatic settings) was enabled; protein identification were accepted with at least one unique + razor peptides for protein group and false discovery rate was set at 0.01 for peptide-spectrum matches, peptides and proteins identifications.
For data analysis, we used the raw intensity from MaxQuant output “proteinGroups.txt”. First, we deleted entries belonging to reverse, only identified by site and Contaminants (only the hit not deriving from the Human fasta file). Then the intensity of each protein was divided by its mass in Dalton and subsequently normalized to the sum of all the intensities of the same file; the normalized numbers were then used for the next analysis and graphs.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/d1na00211b |
‡ Equal contribution. |
This journal is © The Royal Society of Chemistry 2021 |