Open source simulation tool for electrophoretic stacking, focusing, and separation

A neural network model for rapid prediction of analyte focusing in isotachophoresis

We present the development and demonstration of a neural network (NN) model for fast and accurate prediction of whether or not a chosen analyte is focused by an isotachophoresis (ITP) buffer system. The NN model is useful in the rapid evaluation of possible ITP chemistries applicable to analytes of interest. We trained and tested the NN model for univalent species based on extensive data sets of over 10,000 anionic and 10,000 cationic ITP simulations. The NN model uses as inputs the mobilities and the acid dissociation constants of leading electrolyte ion, trailing electrolyte ion, counterion, and a single analyte as well as the leading-to-counterion concentration ratio of the leading zone. The output then indicates whether the chosen electrolyte system yields stable ITP focusing of the analyte. The prediction accuracy of the NN model is over 97.7%. We demonstrate the applicability of the NN by validating its predictions with reported experimental data for anionic and cationic ITP. We have packaged the NN model in a free, web-based application named IONN (isotachophoresis on neural network), which can be used to rapidly screen ITP electrolyte systems.

Fertile Goeppertella from the Jurassic of Patagonia: mosaic evolution in the Dipteridaceae-Matoniaceae lineage

Goeppertella has been postulated as a monophyletic group, whose precise position within the Gleichenoid families Dipteriaceae and Matoniaceae, remains poorly understood. Previously described Goeppertella specimens are based on frond fragments and its fertile morphology is represented by a few, poorly preserved specimens. We describe a new species based on the largest collection of fertile specimens known to date and discuss the evolutionary history of the genus based on the additional reproductive characters provided by the fossils described. Plant impressions were collected in Early Jurassic sediments of Patagonia, Argentina. The specimens were described, and silicone rubber casts were developed to examine in detail vegetative and reproductive features. The new species was compared with other Goeppertella species. Finally, a backbone analysis was performed in the context of a previously published combined matrix of Dipteridaceae, using the maximum parsimony criterion. The new species is described based on a combination of features that have not been previously reported. The vegetative morphology shows affinities with most fossil and extant Dipteriaceae, contrasting with the reproductive morphology which is more comparable with the scarce number of fossil dipteridaceous forms and it is more spread in the sister family, Matoniaceae. The backbone analysis indicates that the position of the new species vary among different positions among Dipteridaceae and Matoniaceae. Additional analyses, discriminating the signal of reproductive and vegetative character, are provided to discuss the base of this uncertainty. We consider Goeppertella as a member of the family Dipteridaceae since we interpret most shared features with Matoniaceae as plesiomorphic conditions for the family. In contrast, most shared features with Dipteridaceae represent apomorphies for the group. Thus, Goeppertella would represent an early diverging genus in Dipteridaceae, considering the venation characters as the most important in order to define the family.

Picomolar detection limits for glyphosate by two-dimensional column-coupled isotachophoresis/capillary zone electrophoresis-mass spectrometry

Capillary electrophoresis-mass spectrometry often lacks sufficient limits of detection for trace substances in the environment due to its low loadability. To overcome this problem, we conducted a feasibility study for column-coupling isotachophoresis to capillary electrophoresis-mass spectrometry. The first dimension isotachophoresis preconcentrated the analytes. The column-coupling of both dimensions was achieved by a hybrid capillary microfluidic chip setup. Reliable analyte transfer by voltage switching was enabled by an in-chip capacitively coupled contactless conductivity detector placed around the channel of the common section between two T-shaped crossings in the chip connecting both dimensions. This eliminated the need to calculate the moment of analyte transfer. A commercial capillary electrophoresis-mass spectrometry instrument with easily installable adaptations operated the setup. Prior to coupling isotachophoresis with capillary zone electrophoresis-mass spectrometry, both dimensions were optimized individually by simulations and verified experimentally. Both dimensions were able to stack/separate all degradation products of glyphosate, the most important herbicide applied worldwide. The first dimension isotachophoresis also removed phosphate, which is a critical matrix component in many environmental samples. Enrichment and separation of glyphosate and its main degradation product aminomethylphosphonic acid by the 2D setup provided an excellent limit of detection of 150 pM (25 ng/L) for glyphosate. This article is protected by copyright. All rights reserved.

Isotachophoresis: Theory and Microfluidic Applications

Isotachophoresis (ITP) is a versatile electrophoretic technique that can be used for sample preconcentration, separation, purification, and mixing, and to control and accelerate chemical reactions. Although the basic technique is nearly a century old and widely used, there is a persistent need for an easily approachable, succinct, and rigorous review of ITP theory and analysis. This is important because the interest and adoption of the technique has grown over the last two decades, especially with its implementation in microfluidics and integration with on-chip chemical and biochemical assays. We here provide a review of ITP theory starting from physicochemical first-principles, including conservation of species, conservation of current, approximation of charge neutrality, pH equilibrium of weak electrolytes, and so-called regulating functions that govern transport dynamics, with a strong emphasis on steady and unsteady transport. We combine these generally applicable (to all types of ITP) theoretical discussions with applications of ITP in the field of microfluidic systems, particularly on-chip biochemical analyses. Our discussion includes principles that govern the ITP focusing of weak and strong electrolytes; ITP dynamics in peak and plateau modes; a review of simulation tools, experimental tools, and detection methods; applications of ITP for on-chip separations and trace analyte manipulation; and design considerations and challenges for microfluidic ITP systems. We conclude with remarks on possible future research directions. The intent of this review is to help make ITP analysis and design principles more accessible to the scientific and engineering communities and to provide a rigorous basis for the increased adoption of ITP in microfluidics.

Electrophoresis simulations using Chebyshev pseudo-spectral method on a moving mesh

We present the implementation and demonstration of the Chebyshev pseudo-spectral method coupled with an adaptive mesh method for performing fast and highly accurate electrophoresis simulations. The Chebyshev pseudo-spectral method offers higher numerical accuracy than all other finite difference methods and is applicable for simulating all electrophoresis techniques in channels with open or closed boundaries. To improve the computational efficiency, we use a novel moving mesh scheme that clusters the grid points in the regions with poor numerical resolution. We demonstrate the application of the Chebyshev pseudo-spectral method on a moving mesh for simulating nonlinear electrophoretic processes through examples of isotachophoresis (ITP), isoelectric focusing (IEF), and electromigration-dispersion in capillary zone electrophoresis (CZE) at current densities as high as 1000 A/m 2 . We also show the efficacy of our moving mesh method over existing methods that cluster the grid points in the regions with large concentration gradients. We have integrated the adaptive Chebyshev pseudo-spectral method in the open-source SPYCE simulator and verified its implementation with other electrophoresis simulators.

Web-Based Open-Source Tool for Isotachophoresis

We present the development of a client-side web-based simulator for complex electrophoresis phenomena, including isotachophoresis. The simulation tool is called Client-based Application for Fast Electrophoresis Simulation (CAFES). CAFES uses the broad cross-browser compatibility of JavaScript to provide a rapid and easy-to-use tool for coupled unsteady electromigration, diffusion, and equilibrium electrolyte reactions among multiple weak electrolytes. The code uses a stationary grid (for simplicity) and an adaptive time step to provide reliable estimates of ion concentration dynamics (including pH profile evolution), requiring no prior installation nor compilation. CAFES also offers a large database of commonly used species and their relevant physicochemical properties. We present a validation of predictions from CAFES by comparing them to experimental data of peak- and plateau-mode isotachophoresis experiments. The code yields accurate estimates of interface velocity, plateau length and relative intensity, and pH variations while significantly reducing the computation time compared to existing codes. The tool is open-source and available for free at https://microfluidics.stanford.edu/cafes.

Electrochemical removal of amphoteric ions

Several harmful or valuable ionic species present in seawater, brackish water, and wastewater are amphoteric, weak acids or weak bases, and, thus, their properties depend on local water pH. Effective removal of these species can be challenging for conventional membrane technologies, necessitating chemical dosing of the feedwater to adjust pH. A prominent example is boron, which is considered toxic in high concentrations and often requires additional membrane passes to remove during seawater desalination. Capacitive deionization (CDI) is an emerging membraneless technique for water treatment and desalination, based on electrosorption of salt ions into charging microporous electrodes. CDI cells show strong internally generated pH variations during operation, and, thus, CDI can potentially remove pH-dependent species without chemical dosing. However, development of this technique is inhibited by the complexities inherent to the coupling of pH dynamics and ion properties in a charging CDI cell. Here, we present a theoretical framework predicting the electrosorption of pH-dependent species in flow-through electrode CDI cells. We demonstrate that such a model enables insight into factors affecting species electrosorption and conclude that important design rules for such systems are highly counterintuitive. For example, we show both theoretically and experimentally that for boron removal, the anode should be placed upstream and the cathode downstream, an electrode order that runs counter to the accepted wisdom in the CDI field. Overall, we show that to achieve target separations relying on coupled, complex phenomena, such as in the removal of amphoteric species, a theoretical CDI model is essential.

Dynamic computer simulations of electrophoresis: 2010-2020

The transport of components in liquid media under the influence of an applied electric field can be described with the continuity equation. It represents a nonlinear conservation law that is based upon the balance laws of continuous transport processes and can be solved in time and space numerically. This procedure is referred to as dynamic computer simulation. Since its inception four decades ago, the state of dynamic computer simulation software and its use has progressed significantly. Dynamic models are the most versatile tools to explore the fundamentals of electrokinetic separations and provide insights into the behavior of buffer systems and sample components of all electrophoretic separation methods, including moving boundary electrophoresis, CZE, CGE, ITP, IEF, EKC, ACE, and CEC. This article is a continuation of previous reviews (Electrophoresis 2009, 30, S16–S26 and Electrophoresis 2010, 31, 726–754) and summarizes the progress and achievements made during the 2010 to 2020 time period in which some of the existing dynamic simulators were extended and new simulation packages were developed. This review presents the basics and extensions of the three most used one‐dimensional simulators, provides a survey of new one‐dimensional simulators, outlines an overview of multi‐dimensional models, and mentions models that were briefly reported in the literature. A comprehensive discussion of simulation applications and achievements of the 2010 to 2020 time period is also included.

Numerical simulations of paper-based electrophoretic separations with open-source tools

A new tool for the solution of electromigrative separations in paper-based microfluidics devices is presented. The implementation is based on a recently published complete mathematical model for describing these types of separations, and was developed on top of the open-source toolbox electroMicroTransport, based on OpenFOAM^® , inheriting all its features as native 3D problem handling, support for parallel computation, and a GNU GPL license. The presented tool includes full support for paper-based electromigrative separations (including EOF and the novel mechanical and electrical dispersion effects), compatibility with a well-recognized electrolyte database, and a novel algorithm for computing and controlling the electric current in arbitrary geometries. Additionally, the installation on any operating system is available due to its novel installation option in the form of a Docker image. A validation example with data from literature is included, and two extra application examples are provided, including a 2D free-flow IEF problem, which demonstrates the capabilities of the toolbox for dealing with computational and physicochemical modeling challenges simultaneously. This tool will enable efficient and reliable numerical prototypes of paper-based electrophoretic devices to accompany the contemporary fast growth in paper-based microfluidics.

HIV detection from human serum with paper-based isotachophoretic RNA extraction and reverse transcription recombinase polymerase amplification

The number of people living with HIV continues to increase with the current total near 38 million, of which about 26 million are receiving antiretroviral therapy (ART). These treatment regimens are highly effective when properly managed, requiring routine viral load monitoring to assess successful viral suppression. Efforts to expand access by decentralizing HIV nucleic acid testing in low- and middle-income countries (LMICs) has been hampered by the cost and complexity of current tests. Sample preparation of blood samples has traditionally relied on cumbersome RNA extraction methods, and it continues to be a key bottleneck for developing low-cost POC nucleic acid tests. We present a microfluidic paper-based analytical device (μPAD) for extracting RNA and detecting HIV in serum, leveraging low-cost materials, simple buffers, and an electric field. We detect HIV virions and MS2 bacteriophage internal control in human serum using a novel lysis and RNase inactivation method, paper-based isotachophoresis (ITP) for RNA extraction, and duplexed reverse transcription recombinase polymerase amplification (RT-RPA) for nucleic acid amplification. We design a specialized ITP system to extract and concentrate RNA, while excluding harsh reagents used for lysis and RNase inactivation. We found the ITP μPAD can extract and purify 5000 HIV RNA copies per mL of serum. We then demonstrate detection of HIV virions and MS2 bacteriophage in human serum within 45-minutes.

Simul 6: A fast dynamic simulator of electromigration

Simul 6 is a 1D dynamic simulator of electromigration based on the mathematical model of electromigration in free solutions. The model consists of continuity equations for the movement of electrolytes in a separation channel, acid-base equilibria of weak electrolytes, and the electroneutrality condition. It accounts for any number of multivalent electrolytes or ampholytes and provides a complete picture about dynamics of electromigration and diffusion in the separation channel. The equations are solved numerically using software means which allow for parallelization and multithreaded computation. Simul 6 has a user-friendly graphical interface. It is typically used for inspection of system peaks (zones) in electrophoresis, stacking and preconcentrating analytes, optimization of separation conditions, method development in either capillary zone electrophoresis, isotachophoresis, and isoelectric focusing. Simul 6 is the successor of Simul 5, and has been launched as a free software available for download at https://simul6.app/.

High-resolution numerical simulations of electrophoresis using the Fourier pseudo-spectral method

We present the formulation, implementation, and performance evaluation of the Fourier pseudo-spectral method for performing fast and accurate simulations of electrophoresis. We demonstrate the applicability of this method for simulating a wide variety of electrophoretic processes such as capillary zone electrophoresis, transient-isotachophoresis, field amplified sample stacking, and oscillating electrolytes. Through these simulations, we show that the Fourier pseudo-spectral method yields accurate and stable solutions on coarser computational grids compared with other nondissipative spatial discretization schemes. Moreover, due to the use of coarser grids, the Fourier pseudo-spectral method requires lower computational time to achieve the same degree of accuracy. We have demonstrated the application of the Fourier pseudo-spectral method for simulating realistic electrophoresis problems with current densities as high as 5000 A/m² with over tenfold speed-up compared to the commonly used second-order central difference scheme, to achieve a given degree of accuracy. The Fourier pseudo-spectral method is also suitable for simulating electrophoretic processes involving a large number of concentration gradients, which render the adaptive grid-refinement techniques ineffective. We have integrated the numerical scheme in a new electrophoresis simulator named SPYCE, which we offer to the community as open-source code.

Nucleic acid sample preparation from whole blood in a paper microfluidic device using isotachophoresis

Nucleic acid amplification tests (NAATs) are a crucial diagnostic and monitoring tool for infectious diseases. A key procedural step for NAATs is sample preparation: separating and purifying target nucleic acids from crude biological samples prior to nucleic acid amplification and detection. Traditionally, sample preparation has been performed with liquid- or solid-phase extraction, both of which require multiple trained user steps and significant laboratory equipment. The challenges associated with sample preparation have limited the dissemination of NAAT point-of-care diagnostics in low resource environments, including low- and middle-income countries. We report on a paper-based device for purification of nucleic acids from whole blood using isotachophoresis (ITP) for point-of-care NAATs. We show successful extraction and purification of target nucleic acids from large volumes (33 µL) of whole human blood samples with no moving parts and few user steps. Our device utilizes paper-based buffer reservoirs to fully contain the liquid ITP buffers and does not require complex filling procedures, instead relying on the natural wicking of integrated paper membranes. We perform on-device blood fractionation via filtration to remove leukocytes and erythrocytes from our sample, followed by integrated on-paper proteolytic digestion of endogenous plasma proteins to allow for successful isotachophoretic extraction. Paper-based isotachophoresis purifies and concentrates target nucleic acids that are added directly to recombinase polymerase amplification (RPA) reactions. We show consistent amplification of input copy concentrations of as low as 3 × 10³ copies nucleic acid per mL input blood with extraction and purification taking only 30 min. By employing a paper architecture, we are able to incorporate these processes in a single, robust, low-cost design, enabling the direct processing of large volumes of blood, with the only intermediate user steps being the removal and addition of tape. Our device represents a step towards a simple, fully integrated sample preparation system for nucleic acid amplification tests at the point-of-care.

Instabilities of the pH gradient in carrier ampholyte-based isoelectric focusing: Elucidation of the contributing electrokinetic processes by computer simulation

Electrokinetic processes that lead to pH gradient instabilities in carrier ampholyte-based IEF are reviewed. In addition to electroosmosis, there are four of electrophoretic nature, namely (i) the stabilizing phase with the plateau phenomenon, (ii) the gradual isotachophoretic loss of carrier ampholytes at the two column ends in presence of electrode solutions, (iii) the inequality of the mobilities of positively and negatively charged species of ampholytes, and (iv) the continuous penetration of carbonate from the catholyte into the focusing column. The impact of these factors to cathodic and anodic drifts was analyzed by simulation of carrier ampholyte-based focusing in closed and open columns. Focusing under realistic conditions within a 5 cm long capillary in which three amphoteric low molecular mass dyes were focused in a pH 3-10 gradient formed by 140 carrier ampholytes was investigated. In open columns, electroosmosis displaces the entire gradient toward the cathode or anode whereas the electrophoretic processes act bidirectionally with a transition around pH 4 (drifts for pI > 4 and pI < 4 typically toward the cathode and anode, respectively). The data illustrate that focused zones of carrier ampholytes have an electrophoretic flux and that dynamic simulation can be effectively used to assess the magnitude of each of the electrokinetic destabilizing factors and the resulting drift for a combination of these effects. Predicted drifts of focused marker dyes are compared to those observed experimentally in a setup with coated capillary and whole column optical imaging.

Electrorotation of leaky-dielectric and conducting microspheres in asymmetric electrolytes and angular velocity reversal

We consider the central problem of polarizable and leaky-dielectric uncharged spherical particle freely suspended in an unbounded nonsymmetric binary electrolyte, which is forced by an ambient time-harmonic uniform electric field. Under the assumption of a "weak field," we employ the linearized standard electrokinetic model of binary electrolytes to account for such anion/cation asymmetry. A simplified generalized asymmetric dipole-term approximation, valid for a dielectric/conducting microsphere, is analytically derived for an arbitrary Debye scale and for any mismatch between ion diffusivities and valances. A two-peak unified dispersion spectrum covering all range of practical frequencies (KHz to MHz), is found for the case of a rotating electric field (ROT). The angular velocity of a free polarized particle is composed of dielectrophoretic contribution, resulting from the electrical torque (dipole term) as well as from the induced electroosmotic (ICEO) flow field. The two effects usually act in opposite directions. Under ROT excitation, we obtain a cofield rotation at high frequencies (MHz) and a counter-field behavior at low frequencies (KHz). The low-frequency dispersion is generally governed by electric double-layer charging and the high frequency by a Maxwell-Wagner relaxation process. ICEO generally dominates the low-frequency cofield response; however, it can be shown that depending on the electrolyte asymmetry, yet another dielectrophoretic related switching (reversal) point might exist. Furthermore, for large frequencies and depending on the complex permittivity ratio between the particle and electrolyte, we find a second switching point. Explicit expressions for the above two frequency reversal values are obtained in terms of the problem physical parameters and are compared against experimental results. Finally, we provide an analytical solution for the ROT ICEO velocity field of a microsphere as a function of electrolyte asymmetry and Debye length and compare it with numerical simulations.

Computer simulation of the isotachophoretic migration and separation of norpseudoephedrine stereoisomers with a free or immobilized neutral chiral selector

A detailed computer simulation study of the isotachophoretic migration and separation of norpseudoephedrine stereoisomers for cases with the neutral selector added to the leader, immobilized to the capillary wall or support, or partially present in the separation column is presented. The electrophoretic transport of the analytes from the sampling compartment into the separation medium with the selector, the formation of a transient mixed zone, the separation dynamics of the stereoisomers with a free or immobilized selector, the dependence of the leader pH, the ionic mobility of norpseudoephedrine, the complexation constant and selector immobilization on steady-state plateau zone properties, and zone changes occurring during the transition from the chiral environment into a selector free leader are thereby visualized in a hitherto unexplored way. For the case with the selector dissolved in the leading electrolyte, simulation data are compared to those observed in experimental setups with coated fused-silica capillaries that feature minimized electroosmosis and zone detection with conductivity and absorbance detectors.

Contemporary chiral simulators for capillary zone electrophoresis

For separation of enantiomers in presence of a chiral selector, data obtained with the 1D dynamic simulators SIMUL5complex and GENTRANS are compared to data predicted by PeakMaster 6, a recently released generalized model of the linear theory of electromigration. Four electrophoretic systems with stereoisomers of weak bases were investigated. They deal with the estimation of input data for complexation together with the elucidation of the origin of observed system peaks, the interference of analyte and system peak migration, the change of enantiomer migration order as function of the selector concentration and the inversion of analyte migration direction in presence of a multiply negatively charged selector. For all systems, data predicted with PeakMaster 6 are in agreement with those of the dynamic simulators and simulation data compare well with experimental data that were monitored with setups featuring conductivity and/or UV absorbance detection along the capillary. SIMUL5complex and GENTRANS provide the full dynamics of any buffer and sample arrangement and require very long execution time intervals. PeakMaster 6 is restricted to conventional CZE, is based on an approximate solution of the transport equations, provides data for realistic experimental conditions within seconds and represents a practical tool for an experimentalist.

Open source capillary electrophoresis

Open source paradigm is becoming widely accepted in scientific communities and open source hardware is finding its steady place in chemistry research. In this review article, we provide the reader with the most up-to-date information on open source hardware and software resources enabling the construction and utilization of an "open source capillary electrophoresis instrument". While CE is still underused as a separation technique, it offers unique flexibility, low-cost, and high efficiency and is particularly suitable for open source instrumental development. We overview the major parts of CE instruments, such as high voltage power supplies, detectors, data acquisition systems, and CE software resources with emphasis on availability of the open source information on the web and in the scientific literature. This review is the first of its kind, revealing accessible blueprints of most parts from which a fully functional open source CE system can be built. By collecting the extensive information on open source capillary electrophoresis in this review article, the authors aim at facilitating the dissemination of knowledge on CE within and outside the scientific community, fosters innovation and inspire other researchers to improve the shared CE blueprints.

Amplification-free detection of DNA in a paper-based microfluidic device using electroosmotically balanced isotachophoresis

We present a novel microfluidic paper-based analytical device (μPAD) which utilizes the native high electroosmotic flow (EOF) in nitrocellulose to achieve stationary isotachophoresis (ITP) focusing. This approach decouples sample accumulation from the length of the channel, resulting in significant focusing over short channel lengths. We provide a brief theory for EOF-balanced ITP focusing under continuous injection from a depleting reservoir and present the design of a short (7 mm) paper-based microfluidic channel, which allows a 200 μL sample to be processed in approximately 6 min, resulting in a 20 000-fold increase in concentration - a full order of magnitude improvement compared to previous paper-based ITP devices. We show the stability of the assay over longer (40 min) durations of time, and using Morpholino probes, we present the applicability of the device for amplification-free detection of nucleic acids, with a limit-of-detection (LoD) of 5 pM in 10 min. Finally, we utilize the small footprint of the channel and show a multiplexed platform in which 12 assays operate in parallel in a 24-well plate format.

Isotachophoresis applied to biomolecular reactions

This review discusses research developments and applications of isotachophoresis (ITP) to the initiation, control, and acceleration of chemical reactions, emphasizing reactions involving biomolecular reactants such as nucleic acids, proteins, and live cells. ITP is a versatile technique which requires no specific geometric design or material, and is compatible with a wide range of microfluidic and automated platforms. Though ITP has traditionally been used as a purification and separation technique, recent years have seen its emergence as a method to automate and speed up chemical reactions. ITP has been used to demonstrate up to 14 000-fold acceleration of nucleic acid assays, and has been used to enhance lateral flow and other immunoassays, and even whole bacterial cell detection assays. We here classify these studies into two categories: homogeneous (all reactants in solution) and heterogeneous (at least one reactant immobilized on a solid surface) assay configurations. For each category, we review and describe physical modeling and scaling of ITP-aided reaction assays, and elucidate key principles in ITP assay design. We summarize experimental advances, and identify common threads and approaches which researchers have used to optimize assay performance. Lastly, we propose unaddressed challenges and opportunities that could further improve these applications of ITP.

High-resolution dynamic computer simulation of electrophoresis using a multiphysics software platform

The modeling and simulation software COMSOL Multiphysics^® was recently extended with an electrophoretic transport interface. Its performance was investigated by comparison to results obtained using the 1D dynamic electrophoresis simulators GENTRANS and SIMUL5. Simulations of zone electrophoresis, isotachophoresis, isoelectric focusing and of an oscillating electrolyte system were performed. Smooth profiles were essentially identical indicating that the COMSOL electrophoretic transport interface is able to reproduce results of the 1D simulators. Differences in the way the respective numerical schemes handle steep concentration gradients and associated instabilities were observed. The COMSOL electrophoretic transport interface is expected to be useful as a general model for simulations in 1D, 2D or 3D geometries, as well as for simulations combining electrophoresis with other physical phenomena.

Electrohydrodynamic instability of ion-concentration shock wave in electrophoresis

Capillary electrophoresis techniques often involve ion-concentration shock waves in an electrolyte solution, propagating under the effect of an external electric field. These shock waves are characterized by self-sharpening gradients in ion concentrations and electrical conductivity that are collinear with the electric field. The coupling of electric field and fluid motion at the shock interface sometimes leads to an undesirable electrohydrodynamic (EHD) instability. Using linear stability analysis, we describe the motion of small-amplitude disturbances of an electrophoretic shock wave. Our analysis shows that the EHD instability results due to the competition between destabilizing electroviscous flow and stabilizing electromigration of the shock wave. The ratio of timescales corresponding to electroviscous flow and electromigration yields a threshold criterion for the onset of instability. We present a validation of this threshold criterion with published experimental data and also describe the physical mechanism underlying the EHD instability of the electrophoretic shock wave.

Modeling of formation and prevention of a pure water zone in capillary isoelectric focusing with narrow pH range carrier ampholytes

This paper comprises a continuation of computer simulation studies dealing with carrier ampholyte based CIEF in presence of narrow pH gradients. With this technique, amphoteric sample components with pI values outside the pH gradient are migrating isotachophoretically toward the cathode or anode whereas components with pI values within the gradient become focused. In order to understand the processes occurring in presence of the electric field, the behavior of both carrier ampholytes and amphoteric sample components is investigated by computer modeling. Characteristics of two pH unit gradients with end components having pI values at or around 7.00 and conditions that lead to the formation of a water zone at neutrality were investigated. Data obtained reveal that a zone of water is formed in focusing with carrier ampholytes when the applied pH range does not cover the neutral region, ends at pH 7.00 or begins at pH 7.00. The presence of additional amphoteric components that cover the neutrality region prevent water zone formation under current flow. This situation is met in experiments with narrow pH gradients that end or begin around neutrality. Simulation data reveal that no water zone evolves when atmospheric carbon dioxide dissolved in the catholyte causes the migration of carbonic acid (in the form of carbonate and/or hydrogen carbonate ions) from the catholyte through the focusing structure. An electrolyte change in the electrode solution does not have an impact on the focusing part but does change the isotachophoretic pattern migrating behind the leading ion.

Simulation and experimental study of enzyme and reactant mixing in capillary electrophoresis based on-line methods

The establishment of an efficient reaction mixture represents a crucial part of capillary electrophoresis based on-line enzymatic assays. For ketamine N-demethylation to norketamine mediated by the cytochrome P450 3A4 enzyme, mixing of enzyme and reactants in the incubation buffer at physiological pH was studied by computer simulation. A dynamic electrophoretic simulator that encompasses Taylor-Aris diffusivity which accounts for dispersion due to the parabolic flow profile associated with pressure driven flow was utilized. The simulator in the diffusion mode was used to predict transverse diffusional reactant mixing occurring during hydrodynamic plug injection of configurations featuring four and seven plugs. The same simulator in the electrophoretic mode was applied to study electrophoretic reactant mixing caused by voltage application in absence of buffer flow. Resulting conclusions were experimentally verified with enantioselective analysis of norketamine in a background electrolyte at low pH. Furthermore, simulations visualize buffer changes that occur upon power application between incubation buffer and background electrolyte and have an influence on the reaction mixture.

Validation of CE modeling with a contactless conductivity array detector

Dynamic computer simulation data are compared for the first time with CE data obtained with a laboratory made system comprising an array of 8 contactless conductivity detectors (C(4) Ds). The experimental setup featured a 50 μm id linear polyacrylamide (LPA) coated fused-silica capillary of 70 cm length and a purpose built sequential injection analysis manifold for fluid handling of continuous or discontinuous buffer configurations and sample injection. The LPA coated capillary exhibits a low EOF and the manifold allows the placement of the first detector at about 2.7 cm from the sample inlet. Agreement of simulated electropherograms with experimental data was obtained for the migration and separation of cationic and anionic analyte and system zones in CZE configurations in which EOF and other column properties are constant. For configurations with discontinuous buffer systems, including ITP, experimental data obtained with the array detector revealed that the EOF is not constant. Comparison of simulation and experimental data of ITP systems provided the insight that the EOF can be estimated with an ionic strength dependent model similar to that previously used to describe EOF in fused-silica capillaries dynamically double coated with Polybrene and poly(vinylsulfonate). For the LPA coated capillaries, the electroosmotic mobility was determined to be 17-fold smaller compared to the case with the charged double coating. Simulation and array detection provide means for quickly investigating electrophoretic transport and separation properties. Without realistic input parameters, modeling alone is not providing data that match CE results.

Theory of multi-species electrophoresis in the presence of surface conduction

Electrophoresis techniques are characterized by concentration disturbances (or waves) propagating under the effect of an electric field. These techniques are usually performed in microchannels where surface conduction through the electric double layer (EDL) at channel walls is negligible compared with bulk conduction. However, when electrophoresis techniques are integrated in nanochannels, shallow microchannels or charged porous media, surface conduction can alter bulk electrophoretic transport. The existing mathematical models for electrophoretic transport in multi-species electrolytes do not account for the competing effects of surface and bulk conduction. We present a mathematical model of multi-species electrophoretic transport incorporating the effects of surface conduction on bulk ion-transport and provide a methodology to derive analytical solutions using the method of characteristics. Based on the analytical solutions, we elucidate the propagation of nonlinear concentration waves, such as shock and rarefaction waves, and provide the necessary and sufficient conditions for their existence. Our results show that the presence of surface conduction alters the propagation speed of nonlinear concentration waves and the composition of various zones. Importantly, we highlight the role of surface conduction in formation of additional shock and rarefaction waves which are otherwise not present in conventional electrophoresis.

Open source and DIY hardware for DNA nanotechnology labs

A set of instruments and specialized equipment is necessary to equip a laboratory to work with DNA. Reducing the barrier to entry for DNA manipulation should enable and encourage new labs to enter the field. We present three examples of open source/DIY technology with significantly reduced costs relative to commercial equipment. This includes a gel scanner, a horizontal PAGE gel mold, and a homogenizer for generating DNA-coated particles. The overall cost savings obtained by using open source/DIY equipment was between 50 and 90%.

Impact of Taylor-Aris diffusivity on analyte and system zone dispersion in CZE assessed by computer simulation and experimental validation

Application of pressure-driven laminar flow has an impact on zone and boundary dispersion in open tubular CE. The GENTRANS dynamic simulator for electrophoresis was extended with Taylor-Aris diffusivity which accounts for dispersion due to the parabolic flow profile associated with pressure-driven flow. Effective diffusivity of analyte and system zones as functions of the capillary diameter and the amount of flow in comparison to molecular diffusion alone were studied for configurations with concomitant action of imposed hydrodynamic flow and electroosmosis. For selected examples under realistic experimental conditions, simulation data are compared with those monitored experimentally using modular CE setups featuring both capacitively coupled contactless conductivity and UV absorbance detection along a 50 μm id fused-silica capillary of 90 cm total length. The data presented indicate that inclusion of flow profile based Taylor-Aris diffusivity provides realistic simulation data for analyte and system peaks, particularly those monitored in CE with conductivity detection.

NAIL: Nucleic Acid detection using Isotachophoresis and Loop-mediated isothermal amplification

Nucleic acid amplification tests are the gold standard for many infectious disease diagnoses due to high sensitivity and specificity, rapid operation, and low limits of detection. Despite the advantages of nucleic acid amplification tests, they currently offer limited point-of-care (POC) utility due to the need for complex instruments and laborious sample preparation. We report the development of the Nucleic Acid Isotachophoresis LAMP (NAIL) diagnostic device. NAIL uses isotachophoresis (ITP) and loop-mediated isothermal amplification (LAMP) to extract and amplify nucleic acids from complex matrices in less than one hour inside of an integrated chip. ITP is an electrokinetic separation technique that uses an electric field and two buffers to extract and purify nucleic acids in a single step. LAMP amplifies nucleic acids at constant temperature and produces large amounts of DNA that can be easily detected. A mobile phone images the amplification results to eliminate the need for laser fluorescent detection. The device requires minimal user intervention because capillary valves and heated air chambers act as passive valves and pumps for automated fluid actuation. In this paper, we describe NAIL device design and operation, and demonstrate the extraction and detection of pathogenic E. coli O157:H7 cells from whole milk samples. We use the Clinical and Laboratory Standards Institute (CLSI) limit of detection (LoD) definitions that take into account the variance from both positive and negative samples to determine the diagnostic LoD. According to the CLSI definition, the NAIL device has a limit of detection (LoD) of 1000 CFU mL(-1) for E. coli cells artificially inoculated into whole milk, which is two orders of magnitude improvement to standard tube-LAMP reactions with diluted milk samples and comparable to lab-based methods. The NAIL device potentially offers significant reductions in the complexity and cost of traditional nucleic acid diagnostics for POC applications.

Taylor dispersion in equilibrium gradient focusing at steady state

An analytic expression is presented for the effective dispersion coefficient in the case where a solute is focused in a parabolic flow against a linear gradient in a restoring force. This expression was derived by employing a minor variation on the method of moments used by Aris in his development of the dispersion coefficients for a time-dependent, isocratic system. In the present case, dispersion is controlled by two dimensionless groups, a Peclet number which is proportional to the parabolic component of the flow, and a gradient number which is proportional to the slope of the restoring force. These results confirm that the Aris-Taylor expression for the dispersion coefficient should not be applied in cases where a solute is focused to a stationary steady state.

1000-fold sample focusing on paper-based microfluidic devices

We present an experimental and analytical study of a novel paper-based analytical device (μPAD) for isotachophoretic sample focusing. Guided by a simple heat transfer model, we further developed wax printing fabrication to enable the creation of shallow channels, which are critical in providing sufficient dissipation of Joule heat, and thus enable the use of high electric fields and short analysis time. This results in a device that is self-contained on a simple piece of filter paper and does not require any specialized enclosures or cooling devices to combat evaporation at high temperatures. Furthermore, we provide an analytical model for isotachophoretic sample accumulation in porous media, introduce a simple figure of merit for evaluating and comparing the efficiency of such devices, and present experimental validation in both paper and glass channels. Using this device we demonstrate the processing of 30 μL of sample achieving 1000-fold increase in peak concentration in 6 min. We believe that this method and device can serve as a guide to the design of low-cost, rapid and highly sensitive paper-based diagnostic platforms.

Theory of ion transport with fast acid-base equilibrations in bioelectrochemical systems

Bioelectrochemical systems recover valuable components and energy in the form of hydrogen or electricity from aqueous organic streams. We derive a one-dimensional steady-state model for ion transport in a bioelectrochemical system, with the ions subject to diffusional and electrical forces. Since most of the ionic species can undergo acid-base reactions, ion transport is combined in our model with infinitely fast ion acid-base equilibrations. The model describes the current-induced ammonia evaporation and recovery at the cathode side of a bioelectrochemical system that runs on an organic stream containing ammonium ions. We identify that the rate of ammonia evaporation depends not only on the current but also on the flow rate of gas in the cathode chamber, the diffusion of ammonia from the cathode back into the anode chamber, through the ion exchange membrane placed in between, and the membrane charge density.

Isotachophoresis with ionic spacer and two-stage separation for high sensitivity DNA hybridization assay

We present an on-chip electrophoretic assay for rapid and high sensitivity nucleic acid (NA) detection. The assay uses isotachophoresis (ITP) to enhance NA hybridization and an ionic spacer molecule to subsequently separate reaction products. In the first stage, the probe and target focus and mix rapidly in free solution under ITP. The reaction mixture then enters a region containing a sieving matrix, which allows the spacer ion to overtake and separate the slower probe-target complex from free, unhybridized probes. This results in the formation of two focused ITP peaks corresponding to probe and probe-target complex signals. For a 149 nt DNA target, we achieve a 220 fM limit of detection (LOD) within 10 min, with a 3.5 decade dynamic range. This LOD constitutes a 12× improvement over previous ITP-based hybridization assays. The technique offers an alternative to traditional DNA hybridization assays, and can be multiplexed and extended to detect other biomolecules.