Samantha
Hume
a,
Gordon
Hithell
a,
Gregory M.
Greetham
b,
Paul M.
Donaldson
b,
Michael
Towrie
b,
Anthony W.
Parker
b,
Matthew J.
Baker
c and
Neil T.
Hunt
*d
aDepartment of Physics, University of Strathclyde, SUPA, 107 Rottenrow East, Glasgow, G4 0NG, UK
bSTFC Central Laser Facility, Research Complex at Harwell, Rutherford Appleton Laboratory, Harwell Campus, Didcot, OX11 0QX, UK
cWestCHEM, Department of Pure and Applied Chemistry, University of Strathclyde, Technology and Innovation Centre, 99 George Street, Glasgow, G1 1RD, UK
dDepartment of Chemistry, York Biomedical Research Institute, University of York, Heslington, York, YO10 5DD, UK. E-mail: neil.hunt@york.ac.uk
First published on 14th May 2019
The amide I infrared band of proteins is highly sensitive to secondary structure, but studies under physiological conditions are prevented by strong, overlapping water absorptions, motivating the widespread use of deuterated solutions. H/D exchange raises fundamental questions regarding the impact of increased mass on protein dynamics, while deuteration is impractical for biomedical or commercial applications of protein IR spectroscopy. We show that 2D-IR spectroscopy can avoid this problem because the 2D-IR amide I signature of proteins dominates that of water even at sub-millimolar protein concentrations. Using equine blood serum as a test system, we investigate the significant implications of being able to measure the spectroscopy and dynamics of proteins in water, demonstrating relevance in areas ranging from fundamental science to the clinic. Measurements of vibrational relaxation dynamics of serum proteins reveals that deuteration slows down the rate of amide I vibrational relaxation by >10%, indicating a dynamic impact of isotopic exchange in some proteins. The unique link between protein secondary structure and 2D-IR amide I lineshape allows differentiation of signals due to albumin and globulin protein fractions in serum leading to measurements of the biomedically-important albumin to globulin ratio (AGR) with an accuracy of ±4% across a clinically-relevant range. Furthermore, we demonstrate that 2D-IR spectroscopy enables differentiation of the structurally similar globulin proteins IgG, IgA and IgM, opening up a straightforward spectroscopic approach to measuring levels of serum proteins that are currently only accessible via biomedical laboratory testing.
In the biomedical arena, spectroscopic interrogation of biofluids holds attractions as a label-free, minimally-invasive screening technology.8 Blood serum is easily obtained, with minimal patient discomfort, and contains a range of potentially diagnostic chemical markers by virtue of contact with most of the major organs.9,10 Current technologies use antibody assays to enhance the signal associated with a target biomolecule, relying critically on the availability of specific antibodies for proteins of interest and requiring significant sample preparation. Moreover, the heterogeneous nature of disease means that single-metabolite detection may be inferior to a broad biomolecular fingerprint of metabolic function as an early warning of deteriorating patient health9,10 or, for example, to indicate the presence of cancers.10,11 The protein content of blood serum represents an ideal substrate for holistic analysis. Human serum contains ∼70 mg mL−1 of proteins composed of albumin (∼35–50 mg mL−1) and the globulins (∼25–35 mg mL−1). Diagnostically, measurement of the albumin to globulin ratio (AGR) is valuable. Changes in the AGR are linked to an inflammatory response and can correlate with the ability of the patient to survive cancer therapy, complementing progress towards personalized treatment and precision medicine.12–16 Currently, the globulin level is derived indirectly, from laboratory-based measurements of total protein and albumin content rather than by direct measurement. Moreover, the globulins encompass a huge number of proteins. The γ-globulins constitute the bulk of the serum globulin fraction. Of these, immunoglobulin-G (IgG) is the most abundant, accounting for ∼80% of the γ-globulins, while IgA (∼13%) and IgM (∼6%) are the next most abundant.16 As well as bulk changes in globulin concentration, changes in serum-levels of each of these individual globulin components are associated with health-related issues. IgG levels are found to increase in cases of liver disease or chronic infection. IgA is linked to cirrhosis while changes in IgM levels can warn of the presence of trypanosomiasis or antibody deficiency syndrome.16
There is thus considerable benefit in a straightforward spectroscopic measurement that can not only deliver the AGR directly but also differentiate between the major globulin components. Current IR spectroscopic studies of blood serum employ dried samples to avoid the problem of water absorption, which can introduce artefacts from the drying process. Attenuated total reflection methods enable studies in aqueous liquids, but such linear spectroscopic methods cannot separate albumin and globulin signals or identify contributions from individual globulin proteins.8
Using blood serum as our exemplar system, we show that H/D exchange is not necessary for label-free protein IR studies in H2O-based media. Typical serum protein concentrations in humans are in the sub mM range (35–50 mg mL−1 albumin corresponds to 0.5–0.7 mM; 25–30 mg mL−1 γ-globulins ∼0.15–0.25 mM), which correspond closely to those used for 2D-IR spectroscopic studies of proteins in D2O. Thus, our approach extends beyond blood serum, allowing detailed amide I studies of proteins in physiological solvents for the first time. Our measurements show that deuteration slows down the vibrational relaxation dynamics of serum proteins, while the enhanced spectral resolution of 2D-IR relative to IR absorption enables accurate differentiation of protein signals in the complex aqueous serum environment, even when their secondary structure composition is similar. These results indicate that the ability to measure the 2D-IR spectrum of proteins in water has relevance in areas ranging from fundamental molecular science to healthcare. Applications of 2D-IR spectroscopy in the biomedical arena are particularly timely in light of recent advances enabling 2D-IR spectral-acquisition in a few seconds17–20 alongside demonstrations of high throughput screening methodologies.21
By contrast to the IR absorption spectra (Fig. 1 and 2(a)), the 2D-IR spectrum of pure serum shows considerable structure (Fig. 2(b)). The negative feature (red) located on the 2D-IR diagonal near 1650 cm−1 is assigned to the v = 0–1 transitions of modes observed in the IR absorption spectrum and contains two distinct contributions with pump frequencies of 1639 and 1656 cm−1 (arrows). Positive (blue) peaks due to the accompanying v = 1–2 transitions are shifted to lower probe frequencies by vibrational anharmonicity.
Comparison of 2D-IR spectra of serum and pure water (Fig. 2(d)) shows that the 2D-IR signal of water is significantly weaker than that of the serum under the same sample conditions. Though weak, the measured 2D-IR response of water was found to be in good agreement with previous observations.24–26 The corresponding IR pump–probe spectrum is also shown for comparison (Fig. 2(c)).27
At 1650 cm−1, the molar extinction coefficient of the serum proteins is at least two orders of magnitude larger than that of water.28 As 2D-IR signals are dependent upon the 4th power of the vibrational transition dipole moment, this leads to enhancement of the strong amide I mode of the biological macromolecules relative to the more plentiful, but weakly-absorbing, water molecules.1,5 As a result, the 2D-IR response of proteins is the dominant feature in the serum spectrum despite the large absorbance of the δH–O–H mode of water.
Comparing the 2D-IR spectrum of pure serum to the spectra of serum albumin and γ-globulins obtained individually in water (Fig. 3) allows assignment of the peaks at 1656 and 1639 cm−1 in the serum spectrum to the albumin and the globulin components respectively. The difference in frequency of the two protein signals arises from the fact that serum albumin has a largely α-helical secondary structure while the globulins have a higher proportion of β-sheet, which shifts the center of mass of the amide I band to lower frequency.3 The higher-order dependence of the 2D-IR signal upon the transition dipole moment means that the lineshapes appearing on the diagonal of a 2D-IR spectrum are narrower than those found in the IR absorption spectrum. This leads to the appearance of two well-resolved peaks along the diagonal of the 2D-IR plot, where only one broad signal was observed in the IR absorption spectrum (Fig. 1(c) and 2(a)).1,5
Fig. 3 IR absorption spectra of (a) serum (b) serum albumin (c) γ-globulins. 2D-IR spectra of (d) serum (e) serum albumin (f) γ-globulins. Dashed grey horizontal lines show peak positions of albumin and γ-globulins. The color scale is as shown in Fig. 2. |
In the case of serum albumin, the γ-globulins and IgG in D2O, the dynamics of the peak of the amide I v = 0–1 transition were well-represented by single exponential decays with lifetimes of 0.89, 0.90 and 0.93 ps respectively (Fig. 4(d)–(f), red). In the case of samples made in H2O, the amide I v = 0–1 peak (Fig. 4(a)) is overlapped by the bleach of the v = 0–1 transition of the H2O δH–O–H mode (Fig. 4(b)). The latter was significantly smaller in amplitude than the protein response (<20% of the total amplitude at Tw = 0 ps) and was well-represented by a bi-exponential function featuring a ∼220 ± 40 fs decay and a ∼1.2 ± 0.2 ps rise-time due to the effects of residual sample heating, which persisted to Tw values longer than 5 ps (Fig. 4(b) and (c)). This behaviour of the water band is in agreement with previous work.27 In the data shown in Fig. 4(d)–(f), the water response has been subtracted by scaling the signals of water and the protein/serum to the Tw signal at 5 ps where no protein contribution is observed. Using this method, single-exponential vibrational relaxation times of 0.78, and 0.74 and 0.78 ps were observed for the amide I v = 0–1 transition for albumin, the γ-globulins and IgG in water. These values were found to be robust using other data analysis approaches, including fitting to tri-exponential functions to account for the bi-exponential relaxation of the water signal and the protein relaxation behaviour.
It is important to note that the vibrational relaxation time of the amide I band of a large protein is, by definition, a weighted average over a large number of coupled amide I oscillators. In the case of the γ-globulins, this is a mixture of proteins. It has been shown previously using fibrillar aggregates of short chain peptides that the lifetime of the amide I band is sensitive to secondary structure and to the level of solvation of a given residue.29 However, the consistent observation of a reduction in vibrational lifetime by around 10% upon moving from D2O to H2O indicates that the average lifetime of the amide I mode of these proteins is being reduced. It is perhaps to be expected that this may be arising from a combination of responses from solvent-exposed and buried residues, which are perturbed differently by H/D exchange, but this is a topic for further study. These findings are however consistent with previous studies of solvent isotope-dependent vibrational dynamics, suggesting that isotopic exchange of the solvent may be responsible for altering the observed protein dynamics.30 In addition to pump–probe data, the spectral diffusion of the 2D-IR amide I lineshape of the proteins in water was compared to that in D2O. Indications of altered linewidths and spectral diffusion processes in H2O were present in the data (ESI Fig. S1†) but as these observations relate to the entire amide I lineshape, which is underpinned by considerable structure arising from different secondary structural contributions a numerical analysis of this as a whole is not physically-relevant.
In relation to analytical studies of complex protein mixtures such as serum, the relative dynamics of the signals due to water and proteins can be exploited to optimise the contrast between the 2D-IR protein response and that of water. The faster relaxation time of the δH–O–H mode than the protein amide I response of serum (∼0.83 ± 0.1 ps at 1650 cm−1) means that, at Tw = 250 fs, the water signal is at a minimum prior to the onset of the small rising signal due to water heating (Fig. 4(c)). The result is that the spectrum shows only the protein signature at this waiting time and on this basis, the following spectral analysis of protein samples was carried out using a Tw of 250 fs, though it is stressed that the water signal at other values of Tw is not sufficiently large as to prevent measurement of protein relaxation dynamics.
(i) Using the relative amplitudes of the peaks assigned to albumin and globulins on the 2D-IR spectrum diagonal (Fig. 5(a) and (b)). In this approach, the ratio of the absolute values of the amplitudes of two distinct peaks at 1656 cm−1 and 1639 cm−1 on the diagonal of each 2D spectrum, assigned to the albumin and globulin fractions respectively, was used to determine the AGR. Scaling of the globulin amplitude by a factor of 1.8 was performed to account for the measured differences in signal amplitude between albumin and the γ-globulins per unit concentration (see Fig. S3†).
(ii) Using the amplitudes of the v = 0–1 peaks due to albumin and globulins taken from pump-frequency slices through the 2D spectra (Fig. 5(c) and (d)). This method utilised slices through the 2D-IR spectrum at pump frequencies of 1656 cm−1 (albumin) and 1639 cm−1 (globulin). The ratio of the absolute values of the amplitude of the globulin pump slice and that of the albumin slice was used to determine the AGR following application of the scaling factor (1.8) to the globulin signal.
(iii) Analysis using a linear combination of the 2D-IR spectra of albumin and γ-globulins (Fig. 5(e) and (f)). All 2D-IR spectra were normalised to the albumin peak at 1656 cm−1. A linear combination analysis fitted the serum 2D-IR spectrum to the linear sum of the independent 2D-IR spectra of albumin and the globulins. The coefficients of the relative contributions of the two protein spectra were then used to evaluate the AGR, following scaling of the globulin fraction by 1.8. The methods are described in detail in the ESI† alongside all 2D-IR spectra (Fig. S2†).
All three methods for determining the AGR spectroscopically produced a linear relationship when the measured value (points and dashed lines, Fig. 5(a), (c) and (e)) was plotted against the actual AGR. The latter was obtained via sending a sample of the as-received equine blood serum for standard laboratory testing at the Glasgow School of Veterinary Medicine and adding the quantity of the known γ-globulin spike to the globulin component. Ideal agreement between the known and 2D-IR-measured AGR values is represented by the solid black line in Fig. 5(a), (c) and (e). Of the three 2D-IR methods used, the pump slice approach (Fig. 5(c) and (d)) was most accurate at the higher values of the AGR, which correspond most closely to the expected human clinical range of 1–2 (horse serum AGR values are slightly lower than human levels). At lower AGR values, the agreement obtained with the pump-slice method was less effective, possibly owing to the very large γ-globulin spike distorting the albumin response. The results obtained from the 2D-IR diagonals were good across the full range of the samples studied (Fig. 5(a) and (b)), with most 2D-IR-derived values being within the measurement error of the actual AGR value, though a constant offset from the actual AGR value was noted. This is attributed to differences in the anharmonicities of the proteins. Finally, the linear combination yielded excellent agreement over the mid-range of the spiked samples (AGR = 0.5–0.7), but was less effective at the extremities. Taking an average of the three analysis approaches (Fig. 5(g)) produced agreement with actual AGR values across the full range of samples, within the experimental uncertainty and is the best approach. Leave one out-type tests of the analysis protocol also showed accuracy to within the expected error of the measurement (Fig. S4†). Overall, the 2D-IR measurements tested here show accuracy over a clinically-relevant range.
Based on sample-to-sample variation, the accuracy of the 2D-IR-derived AGR measurement was ±0.03 (∼4%). Direct comparisons with the current wet assay technique are not possible because these tests derive the AGR value from the difference in total protein and albumin concentrations and so do not directly measure globulin content, however typical quoted accuracies are ∼1%. Although the spectroscopic approach is less accurate than current technologies, this is a first demonstration and there is considerable scope for improvement of the accuracy through engineering approaches to sample path length repeatability and improved data collection protocols.
2D-IR spectra of IgG, IgA and IgM are shown in Fig. 6(a)–(c). As expected, the spectral features are all very similar due to their comparable secondary protein structures and compositions. In order to highlight the subtle discrepancies between the signals, difference 2D-IR spectra (Fig. 6(d)–(f)) were constructed by subtracting the average of triplicate measurements of IgG, IgA and IgM at a known concentration (Fig. 6(a)–(c)) from the spectrum of the equivalent concentration of γ-globulins to account for different maximum solubility levels of the immunoglobulin proteins studied. The results of this process shows that there is little difference between the spectrum of IgG and that of the γ-globulins (Fig. 6(d)). Although somewhat trivial, this expected result acts as an effective control for the process. By contrast, the difference spectra obtained for IgA and IgM do show regions of spectral differences with the γ-globulins. IgA in particular (Fig. 6(e)) clearly shows regions of decreased negative (red) spectral density in the diagonal region near ∼1640 cm−1 and increases in the diagonal part of the amide I band (blue) near 1657 cm−1. IgM also shows spectral differences to the γ-globulins in the diagonal region of the spectrum (Fig. 6(f)), though the effect is less marked than for IgA and the amplitude of the difference signal is reduced.
In order to determine whether the signals observed in the difference spectra are sufficient to quantify changes in serum levels of IgA and IgM, measurements were carried out on a range of serum samples spiked with additional quantities of IgA and IgM. The positions of the peaks in the (γ-globulins–IgA) difference spectrum (Fig. 6(e), 1640 and 1657 cm−1) are not the same as those used to measure the γ-globulin fraction for the above AGR analysis (Fig. 3(f)). Furthermore, the opposing sign of the two components gives two points of reference that can be used to separate the contributions from IgA and γ-globulins.
The 2D-IR spectra of serum samples spiked with concentrations of IgA from 0–15 mg mL−1 are shown in Fig. S5.† To determine the ability of the 2D-IR spectrum to determine the IgA content, the ratio of the amplitudes on the spectrum diagonal at 1657:1640 cm−1 (the positive and negative peaks in the (γ-globulin–IgA) difference spectrum) was plotted as a function of IgA concentration (Fig. 6(g)). The measurement is difficult because the IgA and γ-globulin signals overlap strongly and there is no portion of the spectrum that is unique to either the γ-globulins or to IgA. However, calculation of the spectra expected from this experiment using reconstructions from the individually-measured albumin, γ-globulin, and IgA spectra (Fig. 6(h)) show that the 1657:1640 cm−1 amplitude ratio should decrease as the IgA level was increased if IgA is influencing the signal. Importantly, the gradient of the decrease in this ratio would be significantly shallower than that observed if the generic γ-globulin response increased (Fig. 6(g)). It can be seen from a comparison of the calculated (Fig. 6(h)) and measured (Fig. 6(g)) data that the 2D-IR response recovered matches well with that expected for an IgA-specific signal increase. The gradient of −0.0015 from the experimental data matches well with the value of −0.0012 derived from calculated data. Furthermore, the gradient of the measured amplitude ratio is much closer to that predicted for a change in IgA levels than for a change in γ-globulin fraction. The fact that the correlation with IgA levels persists down to ∼1 mg mL−1 compares well with expected serum levels of IgA, which are in the range of 13% of ∼30 mg mL−1 (4 mg mL−1).16
Repeating the exercise for IgM also showed a linear relationship between the measured 2D-IR signal at the peaks of the (IgM–γ-globulin) difference spectrum and the IgM concentration (Fig. S6†). In the case of IgM, the smaller magnitude of the spectral differences between IgM and γ-globulins led to a significantly more noisy correlation, while the calculated differences show that a less clear separation of the IgM and γ-globulin responses would be anticipated.
The reported proof of concept experiments illustrate the promising potential for development of 2D-IR as a tool for differentiating IgA and IgM protein contributions from albumin and γ-globulins based on their 2D-IR responses. Our approach utilizes two characteristic points of spectral difference that separate IgA/M from albumin and the γ-globulins. Whilst acknowledging that applying this approach to samples of unknown protein levels would be challenging, we believe that this first study gives a clear indicator that the potential exists to sub-divide the γ-globulin spectral component into its major constituents. This result encourages further developments in sample handling protocols and data analysis strategies, including absolute calibration methods needed to extend the capabilities of this approach.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9sc01590f |
This journal is © The Royal Society of Chemistry 2019 |