Emma H.
Yee
a,
Seunghyeon
Kim
a and
Hadley D.
Sikes
*abc
aDepartment of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
bProgram in Polymers and Soft Matter, Massachusetts Institute of Technology, Cambridge, MA 02139, USA
cAntimicrobial Resistance Interdisciplinary Research Group, Singapore-MIT Alliance for Research and Technology, 1 CREATE Way, Singapore 138602. E-mail: sikes@mit.edu
First published on 26th April 2021
Oxygen-tolerant radical polymerization has demonstrated applications in biosensors as a signal amplification method for molecular recognition events. In particular, eosin-mediated photo-redox polymerization, a visible light-initiated radical copolymerization method using N-vinyl pyrrolidone and PEDGA monomers, can be performed in aqueous microliter-scale droplets under atmospheric conditions, and has been used for rapid (≤90 s) signal amplification in several diagnostic assays. In recent years, significant progress has been made in understanding the reaction mechanism, and here we assess the accuracy of the proposed mechanism via experimental validation in an assay format. A 2D reaction–diffusion model was developed and compared to experimental behavior of eosin photopolymerization for paper-based signal amplification. For 3 and 4 mm test zones, the model predicted, within an order of magnitude, experimentally observed effects of oxygen exposure and eosin photoinitiator concentration on polymer formation in a droplet. Both model and experimental results demonstrated that high oxygen exposure and low eosin concentration restrict polymer formation to the center of circular wells. Decreasing the surface-area-to-volume ratio of the reaction droplet and increasing eosin concentration allow polymerization throughout the zone, initially forming in radially intermediate zones due to oxygen's role not just as a reaction inhibitor but also a promoter via photoinitiator regeneration. Reaction volumes as low as 20 μL on 3, 4, and 5 mm diameter reaction zones enabled sensitive signal amplification, although higher oxygen exposure (3–10 μL droplets) showed greatly reduced sensitivity. Observing oxygen tolerance limits and experimentally validating the reaction mechanism can help better understand the eosin photopolymerization system and its applications in diagnostic assay signal amplification.
Many oxygen tolerant methods of performing ATRP,8 RAFT,9–11 redox-mediated polymerization,3 and photoinitiated polymerization15–19 have been demonstrated, enhancing their ability for use in many laboratory-based and point-of-care assays. In particular, visible light-initiated free radical polymerization with eosin Y (eosin dianion) and a tertiary amine has been used as a rapid, low-cost signal amplification method to detect biomarkers of several infectious diseases.16,20–23 In eosin photopolymerization, the absorption of a photon of wavelength in the vicinity of 517 nm excites eosin into a triplet intermediate state where it can undergo electron transfer with a triethanolamine (TEOA) coinitiator.24,25 The electron donation to eosin's triplet excited state via a loosely bound charge-transfer complex yields an eosin radical trianion and a TEOA radical cation26 (Fig. 1A). TEOA radical cations either decay to TEOA radicals, which propagate the radical polymerization of monomers, or react with oxygen to form stable iminium cations and superoxide.27 Eosin radical trianions can also react with oxygen to regenerate eosin dianion (Fig. 1A). Like most radical polymerization reactions, eosin photopolymerization is inhibited by oxygen, which hinders propagation by forming peroxy radicals which are less reactive than carbon-centered free radicals. For polymerization to occur at ambient conditions, a small amount of free eosin Y is added to the monomer solution15 such that oxygen is sufficiently scavenged by eosin radical trianion and TEOA radical cation species in solution.27,28 However, oxygen also plays a uniquely promotional role by regenerating eosin via oxidation of eosin radical trianions, allowing the photointiation process to continue even with low concentrations of eosin photoinitiator.
TEOA radicals catalyze the copolymerization of polyethylene glycol diacrylate (PEGDA) and N-vinylpyrrolidone (VP). VP is commonly employed in acrylate photopolymerization reactions to reduce oxygen inhibition and enhance the rate of radical polymerization and final conversion.25,27,29,30
When used for signal amplification, eosin photoinitiators are coupled to nucleic acids16,22 or proteins16,20,21,23 that specifically recognize captured biomarkers. Local photoinitiator concentration will increase when biodetection events occur, and polymer forms when the local concentration exceeds the minimum initiator threshold where propagation reactions become competitive with inhibition reactions. Polymer formation can enable a colorimetric readout through the use of the pH indicating dye, phenolphthalein, in the monomer solution. Eosin photopolymerization has been previously used as an amplification method for colorimetric detection in diagnostic tests for tuberculosis,20,22 malaria,21 and periodontal disease.23
A thorough understanding of the reaction kinetics and reaction mechanism of eosin photopolymerization promotes its effective use for signal amplification. An ability to predict how polymerization is affected by changes in geometry or changes in key species like inhibitory oxygen and eosin photointiator is crucial in understanding how the signal amplification method can best be applied.
Models have previously been developed to study eosin-mediated photopolymerization at ambient conditions,25,31,32 and some even hypothesized a regeneration mechanism for eosin.31–33 Many insights into the eosin regeneration mechanism have since been made,26,28,34 but have yet to be combined into a kinetic model and validated with experimental results in actual diagnostic test conditions.
Here, we developed a 2D axisymmetric reaction–diffusion model to simulate diagnostic test conditions where a droplet of monomer solution is applied to a paper biodetection zone and irradiated with green light for photoinitiated signal amplification. The reaction–diffusion system models polymerization throughout the droplet, initiated by the free eosin Y in the monomer solution rather than by bound initiators at the sensing interface. It did not incorporate surface-immobilized eosin that occurs in molecular recognition events. The model incorporated key reactions in the proposed mechanism and solved for species concentration throughout the system over time. The model's predictions of when and where in the test zone polymer forms were compared with experimental results. The effects on polymer formation of changing eosin Y photoinitiator concentration and oxygen exposure in the monomer solution were investigated. We then investigated if the model's predicted trends for solution-phase polymerization behavior could also be used to predict and understand polymer formation in biodetection assays where eosin-conjugated streptavidin was captured on paper-based tests by an engineered binding protein.
Test zones were imaged with an iPhone 7 camera. Because ambient lighting and operator positioning affect image color,35 those variables were controlled in the setup. The smartphone was clamped a set distance directly above the test well to be imaged. Lighting conditions were controlled by imaging in a windowless room with constant overhead lighting in a ceiling-less box of white paper to minimize glare.
Ji = −Di∇ci. |
Species balances were solved for each chemical species using reaction terms from the rate equations in the kinetic reaction mechanism. All reactions included in the model, kinetic rate constants, diffusion coefficients, and other parameters were listed in ESI.† To simulate oxygen exposure to droplets at ambient conditions, a flux boundary condition was applied at all outer surfaces of the model using an approximation of Fick's Law at each surface point in the mesh. Here, cO2,0 is the initial dissolved oxygen concentration in the monomer solution.34 The approximate oxygen boundary layer size, δ, was set to 0.1 mm.36
Initial concentration of eosin and all other monomer components in the model was set identical to experimental conditions with which it was compared (see ESI†).
The biodetection assay was created by first immobilizing cellulose binding domain (CBD)-rcSso7d streptavidin binder constructs.38 We previously engineered a charge-neutralized variant of the Sulfolobus solfataricus DNA binding protein, rcSso7d, to specifically bind streptavidin. Fusing it to CBD enables rapid, high efficiency immobilization of the binder on cellulose paper.38,39 Eosin-5-isothiocyanate was conjugated to streptavidin (ESI†) and applied to the assay for specific molecular capture that could be detected via eosin photopolymerization.
In examining the model, we looked not only at concentration of monomer solution species and double bond conversion over time, but also at the Damköhler number (DaO2), a ratio of oxygen reaction rate to diffusion rate throughout the droplet. This helps understand what is happening to oxygen throughout the droplet over time.
Like experimental results (Fig. 2C), the model predicts that under the set conditions, the highest rates of double bond formation and initial polymer formation occur in intermediate areas not at the droplet center but also not at the edge (Fig. 2D). Why conversion is highest in the area can be explained by the roles of oxygen in the system. While oxygen plays a major inhibitory role in the eosin photopolymerization reaction, it also reacts with eosin radicals to regenerate eosin. Polymerization can only occur when oxygen concentration is essentially zero, but that polymerization also relies on the eosin regenerated by oxygen to form initiating TEOA radicals. For this reason we see that the area when oxygen first drops to zero (Fig. 2E) is also the area where oxygen is reacting the fastest (high DaO2) (Fig. 2F), regenerating eosin that then reacts with TEOA to form high concentrations of TEOA and eosin Y radicals (Fig. 2G and H).
Both the model and experiments found that increasing eosin concentration decreases t0.2 and allows polymer to form further from center (Fig. 3). In further validation of our proposed reaction mechanism, the model predicted t0.2 quite accurately, and captured its nonlinear dependence on eosin concentration.
![]() | ||
Fig. 3 In both model and experimental results, increasing eosin Y concentration increased the area in which polymer will form and reduced the irradiation time necessary to form polymer. |
Eosin and the TEOA radicals generated by photoinitiation both react with oxygen to deplete it and allow polymerization to occur. Therefore, higher eosin concentrations increase the rate of oxygen depletion, which reduces t0.2 and allow oxygen-depleting reactions to occur closer to the oxygen-exposed edge of the droplet. However, the area of highest polymerization rate also moves closer to the droplet edge, following the oxygen diffusion front, as it depends on oxygen to regenerate eosin initiators.
It is important to note that most of the calculated t0.2 or “gelation time” is taken up by oxygen inhibition, when radical polymerization cannot occur due to the presence of oxygen in the reaction droplet.18,27 For example, in modeling the case of 0.4 μM eosin Y (Fig. 3), it took 90 seconds for eosin radical trianion and TEOA radical cation reactions to deplete oxygen completely in at least one point of the volume, but it is 115 seconds before double bond conversion reaches 0.2. Because t0.2 is so heavily dependent on oxygen inhibition, the scaling of t0.2 with catalyst (eosin) loading is very high and does not appear to follow the usual propagation–initiation rate scaling of free radical polymerization reactions.40
As expected, due to oxygen inhibition, larger test zones with thinner droplets took longer to polymerize and polymerized in an area closer to the well center, as the region where oxygen could be sufficiently depleted by initiators decreased. This was observed both experimentally and from the model (Fig. 4 and Table 1). It was observed experimentally that significantly less polymer formed in 5 mm-diameter wells (Fig. 4B).
![]() | ||
Fig. 4 Model predictions and experimental effects of increasing test zone diameter and droplet surface-area-to-volume ratio. |
Well diameter (mm) | 3 | 4 | 5 |
Experimental t0.2 (s) | 85 | 95 | 105 |
Model predicted t0.2 (s) | 115 | 157 | Never |
While the model did predict no polymerization (even after 300 s irradiation) for 5 mm-diameter wells, it is likely that radical formation is barely outcompeting oxygen inhibition in the high oxygen exposure conditions of the 5 mm-diameter wells, and the model margin of error prevented accurate prediction in an environment so close to total inhibition. The model could be made more accurate by approximating the droplet as an ellipsoid and measuring more kinetic parameters for the system rather than using approximations from literature. We used sensitivity analysis to quantify the impact of possible inaccuracies in several rate constants and the oxygen boundary layer (Fig. S1†). For example, we calculated that t0.2 would decrease by 5% if kisc in our system is 10% larger than the literature value. However, in spite of these approximations, the model accurately predicted trends in polymer spatial and temporal formation with well diameter and oxygen exposure of reaction droplets.
To simulate actual biodetection applications of eosin photopolymerization and determine if the trends observed in our model would apply to this scenario, biodetection assays were created by immobilizing streptavidin-binding proteins on 3, 4, and 5 mm test zones and capturing eosin-conjugated streptavidin. Volumes of capture binder and streptavidin eosin were adjusted depending on the test zone size such that all test sizes and all replicates had approximately the same surface concentration of eosin (Fig. S2†). Eosin photopolymerization was then used for signal amplification on test zones of various sizes, all with 20 μL monomer solution droplets. Zones were irradiated for a time 5 seconds below t0.2 such that negative samples (no immobilized streptavidin binder) returned no polymer.
In this biodetection format, the same trends previously predicted by the model and observed experimentally were also seen. As the test zone diameter increases, and the droplet volume becomes more exposed to oxygen, polymer formation becomes restricted to areas closer to the droplet center (Fig. 5).
In some cases, it could be advantageous to use shallow droplets that are not sensitive to photoactivation (irradiation) time, and can be irradiated for long periods of time without risk of false positives, or formation of colored polymer due to the free eosin Y in the monomer solution, not due to biorecognition events. However, shallow droplets would be useful only for assays where high sensitivity is not required. Here, it was observed that high oxygen exposure in very shallow droplets greatly reduced sensitivity to molecular recognition events (Fig. 6). Eosin photopolymerization was again performed on biodetection assay test zones where eosin-streptavidin had been specifically captured at a uniform surface concentration for both 3 mm and 5 mm diameter zones. However, this time, 3 μL and 10 μL monomer solution was applied to the zones, respectively. Both 3 mm and 5 mm-diameter zones demonstrated greatly reduced polymer response when a smaller droplet of monomer solution was applied.
![]() | ||
Fig. 6 Signal amplification sensitivity of eosin photopolymerization in 3 and 5 mm biodetection assays was greatly decreased when very shallow reaction droplets were used. |
When eosin photopolymerization is used to amplify signal from eosin-coupled biodetection events, it is a threshold process: polymer forms when local concentration of eosin photoinitiator exceeds the minimum initiator threshold (due to macrophotoinitiator-linked molecular recognition events) where propagation reactions become competitive with inhibition reactions. For maximum sensitivity, the reaction with added solution-based initiator and no surface-bound eosin should be very close to that threshold, such that even a small number of surface molecular recognition events that bind eosin to the zone surface will bump the initiator concentration over the threshold. It is evident that under conditions of high oxygen exposure, the system is not at all close to that threshold, and lots of initiator must first be immobilized before polymerization can occur—greatly reducing sensitivity.
These results can guide application of eosin photopolymerization, and other radical polymerization reactions, in paper-based diagnostic assays. A variety of platforms exist, such as lateral and vertical flow assays, where solutions flow via capillary action through paper channels where reactions occur. Sliding chip assays, where paper reaction zones are manually slid in and out of fluidic paths to enable incubation with reagent droplets, have also been proposed. Another format that allows for reactions to occur in droplets is 96-well printed sheets,41 similar to traditional 96-well plates commonly used for ELISA and other high-throughput screening assays.
The results presented in this paper indicate that eosin photopolymerization sensitivity is promoted in monomer droplets with lower surface-area-to-volume ratios, which minimize oxygen exposure. This indicates that eosin photopolymerization signal amplification is most applicable in assay platforms where droplets of at least 10 s of μL are feasible, such as high throughput screening assays in 96-well printed formats. For more portable, point-of-care applications, devices like sliding chip assays,42 where solutions can be localized in small areas, could be compatible with this amplification method. This is likely also true for other oxygen-inhibited radical polymerization reactions that can be performed at ambient conditions, such as enzyme- and redox pair-mediated free radical polymerization,1,3 and some modified forms of ATRP8 and RAFT.9–11
Eosin photopolymerization and other radical polymerization methods would likely be incompatible with flow-through formats (lateral and vertical flow assays) where monomer solutions would be highly exposed to oxygen.43 It is possible that radical polymerization could work in flow-through formats for assays targeting molecules in such high abundance that low sensitivity is acceptable. Additionally, in some cases, greatly increasing initiator concentration or light intensity could overcome inhibitory reactions.
Trends observed in model and experimental results demonstrate that careful consideration of assay platform compatibility with radical polymerization reactions—even relatively oxygen-tolerant ones—is necessary when determining if or how radical polymerization can be applied for signal amplification. While this paper specifically investigated eosin-mediated photo-redox polymerization reactions, its conclusions regarding the effects of oxygen exposure and initiator concentration on spatial polymer formation and sensitivity could be generalized to many other oxygen-inhibited free radical polymerization reactions, such at various ATRP, RAFT, and redox-initiated radical polymerization methods, which are only oxygen tolerant to a certain point. While concentration of initiator or other methods of suppressing oxygen inhibition can be varied for specific reaction conditions, in some applications radical polymerization may not allow achievement of desired sensitivity and/or reaction speed. This work illustrates the value of appropriately matching signal amplification chemistry to the desired assay format and performance needs.
Footnote |
† Electronic supplementary information (ESI) available: Includes model equations, variables, and initial conditions, materials, eosin conjugation and imaging methods, and supplementary model vs. experimental results. See DOI: 10.1039/d1py00413a |
This journal is © The Royal Society of Chemistry 2021 |