Simul 6: A fast dynamic simulator of electromigration

Methamphetamine detection using portable capillary electrophoresis coupled with a swab-based extraction device

Methamphetamine (MA) is one of the most virulent illicit drugs that can be synthesized from household materials leading to its prevalent trafficking and local manufacturing in clandestine drug laboratories (clan labs). The significant problems of tracing MA in clan labs and monitoring drug abusers lie in the lag time between sample collection and analysis and the number of tests done. Capillary electrophoresis (CE) is a rapid separation technique amenable to miniaturization and field testing. Herein, we developed a simple transient isotachophoretic (tITP)-CE method to detect MA and its precursor pseudoephedrine (PSE) in clan labs and non-invasive biological fluids. The method was implemented on the ETD-100, a commercial fully automated portable CE instrument with an integrated swab-based extraction system. Within 2 min of insertion of the swab, MA and PSE were automatically extracted with a leading electrolyte (LE) and then separated on covalently modified capillaries. The ETD-100 showed a limit of detection (LOD) and quantification (LOQ) of MA 0.02 and 0.05 μg/swab and 0.02 and 0.06 μg/swab of PSE, with an enhancement factor of 118 and 328, respectively, when compared to a normal non-tITP injection. The intra and inter-day relative standard deviation in terms of migration time were in the range of 0.75-1.93 % for both MA and PSE and were 2.0-2.4 % for both MA and PSE peak height. The method was demonstrated with the detection of spiked MA and PSE on different household materials as well as in non-invasive biological fluids with a recovery above 60 %.

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.

Mobilization in two-step capillary isoelectric focusing: Concepts assessed by computer simulation

The mobilization step in a two-step capillary isoelectric focusing protocol is discussed by means of dynamic computer simulation data for systems without and with spacer compounds that establish their zones at the beginning and end of the focusing column. After focusing in an electroosmosis-free environment (first step), mobilization (second step) can be induced electrophoretically, by the application of a hydrodynamic flow, or by a combination of both means. Dynamic simulations provide insight into the complexity of the various modes of electrophoretic mobilization and dispersion associated with hydrodynamic mobilization. The data are discussed together with the relevant literature.

Relevance of the calculation of the diffusion coefficient in a capillary electrophoresis experiment

This paper analyzes the role of the diffusion coefficient in the movement of analytes that can reversibly react with a selector given a product in the presence of drift. The problem mimics the movement of enantiomers in a capillary electrophoresis experiment. As is well known, the signal in the capillary must be sharp enough to make a good determination of the effective mobility of the analytes being analyzed. The essence of the technique is based on fast interconversion rates. Therefore, the effective diffusion coefficient must be negligible during the experiment. In the present work, an exact expression for both the apparent mobility and the diffusion coefficient is obtained. This is done by writing the rate equations governing the process and solving them using the generating function technique. The effective mobility coincides with the Wren and Rowe equation, whereas the diffusion coefficient allows us to determine the values of the parameters to be taken into account so that this quantity is minimal or close to zero. On the other hand, the numerical solution of the kinetic equations and Monte Carlo simulations allow us to follow the signal in the capillary and to determine its space-time evolution.

Streamlined Data Processing for Determination of Equilibrium Dissociation Constants with Accurate Constant via Transient Incomplete Separation (ACTIS)

The determination of accurate equilibrium dissociation constants, K_d, of protein-small molecule complexes is important but challenging as all established methods have inherent sources of inaccuracy. Accurate Constant via Transient Incomplete Separation (ACTIS) is a new method for K_d determination using transient incomplete separation of the complex from the unbound small molecule in a pressure-driven flow inside a capillary. ACTIS is accurate, and its accuracy is invariant to variations in geometries of both the fluidic system and the flow. Furthermore, ACTIS is implemented using a simple fluidic system supporting its accuracy and providing a simple-to-follow/copy template for instrumentation. Despite the simple and robust instrumentation/acquisition, the current data processing workflow is cumbersome, time consuming, and prone to hard-to-trace human errors therefore hindering ACTIS' ability to become a practical reference method for K_d determination. This technical note describes a streamlined workflow for processing ACTIS data; the workflow is implemented as a set of open-source software tools called prACTISed (https://github.com/prACTISedProgram/prACTISed). These tools allow all steps of data processing to be performed in a fast and straightforward fashion. These practical software tools complement the simple instrumentation serving both developers and users of ACTIS.

Evolution of the theoretical description of the isoelectric focusing experiment: II. An open system isoelectric focusing experiment is a transient, bidirectional isotachophoretic experiment

The carrier ampholytes-based (CA-based) isoelectric focusing (IEF) experiment evolved from Svensson's closed system IEF (constant spatial current density, absence of convective mixing, counter-balancing electrophoretic and diffusive fluxes yielding a steady state pH gradient) to the contemporary open system IEF (absence of convective mixing, large cross-sectional area electrode vessels, lack of counter-balancing electrophoretic- and diffusive fluxes leading to transient pH gradients). Open system IEF currently is described by a two-stage model: In the first stage, a rapid IEF process forms the pH gradient which, in the second stage, is slowly degraded by isotachophoretic processes that move the most acidic and most basic CAs into the electrode vessels. An analysis of the effective mobilities and the effective mobility to conductivity ratios of the anolyte, catholyte, and the CAs indicates that in open system IEF experiments a single process, transient bidirectional isotachophoresis (tbdITP) operates from the moment current is turned on until it is turned off. In tbdITP, the anolyte and catholyte provide the leading ions and the pI 7 CA or the reactive boundary of the counter-migrating H₃ O⁺ and OH^- ions serves as the shared terminator. The outcome of the tbdITP process is determined by the ionic mobilities, pK_a values, and loaded amounts of all ionic and ionizable components: It is constrained by both the transmitted amount of charge and the migration space available for the leading ions. tbdITP and the resulting pH gradient can never reach steady state with respect to the spatial coordinate of the separation channel.

Evolution of the theoretical description of the isoelectric focusing experiment: III. Carrier ampholyte behavior in transient, bidirectional isotachophoresis

In modern isoelectric focusing (IEF) systems, where (i) convective mixing is prevented by gels or small cross-sectional area separation channels, (ii) current densities vary spatially due to the presence of electrode vessels with much larger cross-sectional areas than those of the gels or separation channels, and (iii) electrophoretic and diffusive fluxes do not balance each other, stationary, steady-state pH gradients cannot form (open-system IEF). Open-system IEF is currently described as a two-stage process: A rapid IEF process forms the pH gradient from the carrier ampholytes (CAs) in the first stage, then isotachophoretic processes degrade the pH gradient in the second stage as the extreme pI CAs are moved into the electrode vessels where they become diluted. Based on the ratios of the local effective mobilities and the local conductivities ( μ L eff ( x ) $\mu _{\rm{L}}^{{\rm{eff}}}( x )$ / κ ( x ) $\kappa ( x )$ values) of the anolyte, catholyte, and the CAs, we pointed out in the preceding paper (Vigh G, Gas B, Electrophoresis 2023, 44, x x x) that in open-system IEF, a single process, transient, bidirectional isotachophoresis (tbdITP) operates from the moment current is turned on. In this paper, we demonstrate some of the operational features of the tbdITP model using the new ITP/IEF version of Simul 6.

On-capillary fluorescent labeling of saccharides for capillary electrophoresis

The feasibility of on-capillary derivatization of saccharides by aromatic amine-based fluorescent labeling agents was tested. To avoid the problematic evolution of gaseous hydrogen cyanide, the Schiff base reduction by sodium cyanoborohydride, as the second step of the standard reductive amination protocol, was omitted. Glucose was used as a model analyte and 7-amino-1,3-naphthalenedisulfonic acid as the labeling agent. Our experiments showed that the direct reaction of the saccharide with the labeling agent in 2.5-M acetic acid yields a labeled product that is sufficiently stable to be separated from the labeling agent in 20-mM phosphate buffer, pH 3.5, and detected using UV detection. The glucose and label zones were introduced separately into the capillary and mixed using a negative voltage. Mixing voltage, its duration, the concentration of acetic acid in the reaction zone, and the waiting time between mixing and separation were optimized. To show the applicability of the procedure to a broader range of analytes, a mixture of different types of saccharides, that is, xylose (pentose), fucose (hexose), glucose (hexose), N-acetylglucosamine (N-acetylaminosaccharide), and lactose (disaccharide), was subjected to derivatization and analysis under the optimal conditions. The linearity and repeatability of the process were evaluated as critical parameters for its analytical applications. Six-point calibration dependences in the 1-50 mM range showed excellent determination coefficients of 0.9992 or higher for all five saccharides tested. The repeatability of the labeled saccharide peak areas was between 2.2% and 4.3%.

Recent progress in analytical capillary isotachophoresis (2018 - March 2022)

This review brings a survey of papers on analytical capillary and microchip isotachophoresis published since 2018 until the first quarter of 2022. Theoretical papers extending fundamental knowledge include those on computer simulations that remain an important research tool useful in the design of electrolyte systems. Many papers are focused on instrumental aspects where new media including microfluidic devices and their hyphenation to various detection techniques bring remarkable results. Papers reporting analytical applications demonstrate the potential of contemporary analytical isotachophoresis. Although it is not being used on a mass scale, its special features are attracting continued interest resulting in applications of isotachophoresis both as a stand-alone analytical method and as a part of multidimensional separation techniques.

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.

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.