Christian
Bortolini
ab,
Tadas
Kartanas
a,
Davor
Copic
c,
Itzel
Condado Morales
a,
Yuewen
Zhang
a,
Pavan K.
Challa
a,
Quentin
Peter
a,
Tamás
Jávorfi
d,
Rohanah
Hussain
d,
Mingdong
Dong
b,
Giuliano
Siligardi
*d,
Tuomas P. J.
Knowles
*ae and
Jérôme
Charmet
*f
aChemistry Department, University of Cambridge, Lensfield Road, Cambridge, CB3 0FF, UK
bInterdisciplinary Nanoscience Center, Aarhus University, Aarhus, 8000, Denmark
cInstitute for Manufacturing, Engineering Department, University of Cambridge, Charles Babbage Road, Cambridge, CB3 0FS, UK
dDiamond Light Source, Harwell Science and Innovation Campus, Didcot, OX11 0DE, UK. E-mail: giuliano.siligardi@diamond.ac.uk
eCavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge, CB3 0HE, UK. E-mail: tpjk2@cam.ac.uk
fInstitute of Digital Healthcare, WMG, University of Warwick, Coventry, CV4 7AL, UK. E-mail: j.charmet@warwick.ac.uk
First published on 22nd November 2018
Circular dichroism spectroscopy has become a powerful tool to characterise proteins and other biomolecules. For heterogeneous samples such as those present for interacting proteins, typically only average spectroscopic features can be resolved. Here we overcome this limitation by using free-flow microfluidic size separation in-line with synchrotron radiation circular dichroism to resolve the secondary structure of each component of a model protein mixture containing monomers and fibrils. To enable this objective, we have integrated far-UV compatible measurement chambers into PDMS-based microfluidic devices. Two architectures are proposed so as to accommodate for a wide range of concentrations. The approach, which can be used in combination with other bulk measurement techniques, paves the way to the study of complex mixtures such as the ones associated with protein misfolding and aggregation diseases including Alzheimer's and Parkinson's diseases.
Bulk measurement techniques such as circular dichroism (CD),2,3 infrared (IR) spectroscopy,4,5 nuclear magnetic resonance (NMR)6,7 or thioflavin T (ThT) fluorescence8–11 have been used to study protein folding and aggregation. However, such techniques only report on the ensemble average and therefore do not allow for a precise resolution of interactions and changes at the molecular level.
Single molecule techniques, such as atomic force microscopy (AFM),12,13 electron microscopy (EM),14,15 infrared nanospectroscopy16 and high resolution imaging17,18 have received a considerable interest in recent years to study protein folding and aggregates at the molecular level. For example AFM and transmission electron microscopy (TEM) are commonly used to study peptide self-assembly through the mapping of sample morphology at different aggregation stages.19–21 To date, however, the level of structural information available from such approaches has been more limited than from bulk spectroscopy. Besides, such measurement often require operation in non-native environment.
Even though conventional separation techniques, such as liquid chromatography, have been combined with bulk measurement techniques to resolve complex mixtures22–26 their operation usually results in a loss of the temporal information and native environment. For example, it has been acknowledged that the interaction with the stationary phase and the dilution of the samples in size exclusion chromatography may modify the state of the sample.27,28
In recent years, microfluidic approaches have opened up new opportunities to study complex biological processes.29–32 Here, we take advantage of the laminar flow properties inherent to operation in microfluidic devices to separate a complex mixture into well-resolved size-dependent fractions, using an H-filter configuration.33 In brief, we flow the solution of interest containing the protein mixture, alongside a buffer solution. The different components of the mixture diffuse in a size-dependent manner, into the buffer solution and the free-flowing solution is then separated into well resolved fractions (see Fig. 1). The concentration of the isolated fraction is obtained based on the concentration profile (Fig. 1c) calculated by numerically solving the problem of particles diffusing in a fully developed Poiseuille flow in a rectangular microfluidics channel (see Calculation of the diffusion profile section). The fractions are then studied with a label-free bulk measurement technique, in this case synchrotron radiation circular dichroism (SRCD), a powerful technique to study the secondary structure of chiral molecules (see Synchrotron radiation circular dichroism section). In particular, we use the highly collimated microbeam generated at Diamond B23 beamline for SRCD,34 which enables on-chip measurement. We show that this combination gives information that could not be obtained by studying the complex mixture without separation. In particular, we demonstrate that by excluding larger species from one fraction (fraction 2 in Fig. 1), it is possible to resolve precisely its structure and reconstruct the structure of the other species (fraction 2) by subtracting the spectrum of the resolved fraction from the that of the overall mixture. Even though the concept proposed herein is used with SRCD, it is can also be adapted to other bulk measurement techniques, such as UV/vis, IR absorption and fluorescence microscopy.
A challenge encountered in interfacing microfluidics with CD is the incompatibility of conventional polydimethylsiloxane (PDMS) based microfluidic devices with Dd-UV measurement (see Fig. S1†). Even though the combination of SRCD and microfluidics has been reported previously to study protein refolding kinetics of cytochrome C from 4 M to 0.8 M GuHCL,35 the mixing devices used in these studies were made of fused silica with the beamlight focused on a masked slit of 60 μm × 15 mm. The fabrication of such devices requires access to specialised microfabrication equipment and expertise, which is usually not readily available in conventional biophysical laboratories. Since soft-lithography is one of the most widely used technique to fabricate microfluidic devices,36,37 and in an effort to make our finding available to a broad scientific community, we propose here microfluidic devices fabricated using conventional PDMS-based soft-lithography, compatible with far-UV measurement, including SRCD. In this manuscript, we describe proof-of-principle microfluidic devices based on two architectures enabling the integration of measurement chambers of different height, thus allowing the measurement of a wide range of concentration (see Fig. 2).
The first architecture consists of a two-layer device that comprise, at appropriate positions along the microfluidic circuit, measurement chambers confined vertically between two quartz slide windows (Fig. 2a). The fabrication process consists in aligning the 2 complementary PDMS devices, with channels facing out, before plasma bonding. The connections between the 2 layers and to the inlet and outlet ports are made using a biopsy punch. Finally, quartz windows are plasma bonded on each side of the device to seal the channels. The second set of devices is fabricated using a single layer architecture and a one-mould process (Fig. 2b). The fabrication steps consist in pouring uncured PDMS in the master mold and pressing a 5 × 5 mm quartz window (cut from a quartz slide using a diamond scriber) onto the measurement chamber area. Once the device is cured, the structured PDMS is peeled gently, making sure the quartz window stays in place, and the ports are punched. Finally, a quartz slide is plasma bonded onto the device to seal the channels (see Fig. S4†). The second architecture was also tested successfully using a 3D printed master (see Fig. S5.a†).
(1) |
(2) |
This equation can be solved analytically. The following formulation of the solution is symmetrical and converges quickly:39
(3) |
The general convection–diffusion equation for the local concentration c(x, y, z) and diffusion coefficient D is given by:
(4) |
It is similarly simplified by translational invariance of the flow and time invariance:
(5) |
Assuming that diffusion in the x direction is negligible (D∂x2c ≪ vx∂xc):
(6) |
This equation is numerically integrated with a trapezoid method and Neumann boundary conditions. The space step dx is chosen with the Courant–Friedrichs–Lewy condition:
(7) |
This choice of dx means the step matrix S is independent on D and Q, the flow rate, given by . Defining the dimensionless step size dx′ = dx × D/Q, the number of steps to reach a position L is found to be:
(8) |
The evolution of an initial concentration distribution c0 can therefore be quickly calculated by repeated matrix squaring. Indeed using Si = S2i, only at most log2(Nsteps) matrix multiplications are necessary. For example 50 is 110010 in binary, so c50 = S6·(S5·(S1·c0)). Fig. S6† shows concentration profiles obtained using the above for insulin monomers.
In the first set of separation experiments, solutions of insulin monomers (0.2 mg ml−1) or fibres (0.4 mg ml−1) were injected alongside a buffer solution (water) and SRCD spectra were measured in the two measurement chambers, each collecting 50% of the solution, after a diffusion length of 90 mm (Fig. 3a). The high diffusion chamber collects a fraction of the small, high diffusion coefficient molecules, while the low diffusion measurement chambers retains most of the larger molecules, with a lower diffusion coefficient (Fig. 3a and b). The flow rates were 30 and 300 μl h−1 for the insulin and the buffer solutions respectively. The devices were made using the first architecture and the heights of the measurement chambers, corresponding to the path length, ranged typically between 3 and 5 mm. The spectra obtained for the monomers, with a typical high α-helical content (see Fig. 3c), and fibres, with a stronger β-sheet signal (see Fig. S7†), are in good agreement with the results obtained using benchtop CD instrument in static mode and confirm that neither the devices, nor the flow rate, induce any artefact (see Fig. S8† for the spectra and secondary structure of undiluted solutions using a benchtop CD instrument).
In order to verify the efficiency of the diffusion-based separation, we have measured the concentration of each fraction and compared it with the expected theoretical value. The concentration ratios were measured by dividing the amplitude of the CD signal at 208 nm (monomer) and 222 nm (fibre) in each chamber by the corresponding amplitude of the total concentration. These values were compared to the concentration expected due to diffusion, as calculated using the area under the theoretical concentration profile (Fig. 3b). The measured and expected (calculated) fraction in each chamber show excellent agreement as detailed in Table 1.
Measured | Calculated | |
---|---|---|
Amyloid fibrils (222 nm) | 2.3% | 1.1% |
Monomers (208 nm) | 13.8% | 11.7% |
We also measured the concentration profile of a solution of insulin monomers (2 mg ml−1) across a 2000 μm wide microfluidic channel, using the second architecture with a 50 μm high channel, and flow rates of 30 and 300 μl h−1 for the monomers and buffer solution respectively. The SRCD signal amplitude at 208 nm was measured continuously for 1 minute in 5 positions across the channel width. The normalised average amplitude points (and standard deviation) show a good match with the expected concentration profile as shown in Fig. 4.
Finally, we performed a time-dependent dilution experiment, using BSA (0.2 mg ml−1, 60 μl h−1) in a Y-junction channel fabricated with a 3D printed master mould based on the second architecture. The measurement, limited by the time needed to acquire a full spectrum, is performed at the end of the channel, across its entire width, where the original sample gets diluted due to the progressive introduction of the aqueous buffer solution from 0 μl h−1 to 60 μl h−1. The spectra shown in Fig. 5 do not exhibit any distortions, as confirmed after normalisation (Fig. S5.b†). For faster experiments, continuous measurement at a single wavelength can be performed.
The validation results presented above, obtained using simple protein solutions, constitute a proof that the different microfluidics device architectures proposed are compatible with SRCD.
Using the H-filter microfluidic separation device architecture presented herein, we isolate a fraction of monomers from the mixture. Using an 18 cm long diffusion channel, one can direct 22% of monomers and no fibres (0%) in the high-diffusion chamber by collecting one third of the flow (see Fig. 6a, right hand-side spectrum). The amplitude of the spectrum collected is then adjusted by a multiplication factor (in this case, 100/22 = 4.55) to account for the 100% of monomer present in the mixture and subtracted from the spectrum obtained for the mixture. The resulting reconstructed spectrum, shown in Fig. 6b, compares well with the spectrum of a fibril solution (of the same concentration) measured in the same chamber. This observation is confirmed by the analysis of the secondary structure, calculated using BeStSel, a method for the secondary structure determination and fold recognition from protein circular dichroism spectra.42Fig. 6c shows that the α-helical content for the reconstructed and measured spectra are identical (6 and 7% for the reconstructed and measured spectra respectively). The β-sheet and unordered contents, however, show small discrepancies. The β-sheet content is of 36% for the reconstructed and 42% for the measured spectra while the unordered content is of 42 and 38% for the reconstructed and measured spectra respectively. The differences observed can be attributed to small distortions of the spectrum that can arise due to variation of the microfluidics chip position with respect to the photon beam between the measurements, the measurement in dynamic mode and resulting flow fluctuations, as well as the low concentration collected, imposed by the necessity to exclude fibres from the high-diffusion chamber. It is noted also that even though monomers are stable at room temperature and neutral pH, they can start to aggregate with the fibrils, which act as catalysts.43 However, in this case, such conversion is very slow and, therefore, the content of the solutions studied is expected to be stable for the duration of the experiments presented herein (see Fig. S9† for details).
The approach described above can also be adapted to resolve more complex mixtures. Microfluidic diffusion-based separation enables the exclusion of species above a given hydrodynamic radius, in a single separation step. Even though it would theoretically be possible to isolate the smallest molecule in the complex mixture (e.g. a monomer), it has been shown that it is not possible to precisely resolve biomolecules unless their hydrodynamic radii differ by at least a factor three.44 Therefore, in order to increase the separation resolution, one should combine the approach with other high resolution microfluidic separation techniques, such as free-flow electrophoresis that enables the selection of biomolecules based on their electrophoretic mobility.45–47 However, ultimately, the detection of a single constituent from a mixture will be limited by the sensitivity of the measurement technique. Nevertheless, the separation of a heterogeneous mixture into less complex, yet well-resolved, fractions is expected to provide further insights into complex biological phenomena.
In order to make the results presented in this study available to a broad scientific community, we have developed new fabrication methods to integrate far-UV compatible measurement chambers (confined between quartz windows) into PDMS based microfluidic devices fabricated using conventional soft-lithography approaches. Two device architectures, which enable the measurement of a wide range of concentrations, are presented, characterised and validated using the highly collimated and high photon-flux microbeam light available at Diamond Light Source B23 beamline for SRCD.
The device architecture presented herein can also be used in combination with other sensing modalities requiring far-UV transparency and more generally, the principle can also be adapted to other bulk measurement techniques. In summary, the possibility to precisely separate a heterogeneous mixture into well-resolved, simpler fractions opens up a range of opportunities for the study of complex biological phenomena. In particular, such developments open up interesting perspective to study protein misfolding and aggregation diseases such as Alzheimer's and Parkinson's diseases and may provide time-resolved information about the protein aggregation pathway.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c8lc00757h |
This journal is © The Royal Society of Chemistry 2019 |