Capillary Electrophoresis Sensitivity Enhancement Based on Adaptive

Apr 30, 2018 - The Frcutoff is required for calculation of optimal averaging window size for a .... Strong correlation (Supporting Information Figure ...
0 downloads 3 Views 3MB Size
Subscriber access provided by Kaohsiung Medical University

Article

Capillary Electrophoresis Sensitivity Enhancement Based on Adaptive Moving Average Method Tomas Drevinskas, Laimutis Telksnys, Audrius Maruska, Jelena Gorbatsova, and Mihkel Kaljurand Anal. Chem., Just Accepted Manuscript • DOI: 10.1021/acs.analchem.8b00664 • Publication Date (Web): 30 Apr 2018 Downloaded from http://pubs.acs.org on May 3, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Capillary Electrophoresis Sensitivity Enhancement Based on Adaptive Moving Average Method

Tomas Drevinskas1,2, Laimutis Telksnys2,3, Audrius Maruška1*, Jelena Gorbatsova4, Mihkel Kaljurand4,1

1

Instrumental Analysis Open Access Centre, Faculty of Natural Sciences, Vytautas Magnus University, Vileikos 8, LT44404 Kaunas, Lithuania 2

Department of Systems’ Analysis, Faculty of Informatics, Vytautas Magnus University, Vileikos 8, LT44404 Kaunas, Lithuania 3

Recognition Processes Department, Institute of Mathematics and Informatics, Goštauto 12, LT01108 Vilnius, Lithuania 4

Department of Chemistry, Faculty of Sciences, Tallinn University of Technology, Akadeemia tee 15, 12618 Tallinn, Estonia

Keywords: Capillary electrophoresis, Noise filtering, Contactless conductivity detection, Instrumentation, Software, Migration velocity

Correspondence to: Prof. Audrius Maruška [email protected]

Abstract: In the present work we demonstrate the novel approach to improve the sensitivity of the “out of lab” portable capillary electrophoretic measurements. Nowadays, many enhancement methods are: (i) underused (non-optimal), (ii) overused (distorts the data), or (iii) inapplicable in field-portable instrumentation due to lack of computational power. Described innovative migration velocity-adaptive moving average method uses optimal averaging window size and can be easily implemented with microcontroller. The contactless conductivity detection was used as a model for the development of a signal processing method and the demonstration of its impact on the sensitivity. The frequency characteristics of the recorded electropherograms and peaks were clarified. Higher electrophoretic mobility analytes exhibit higher frequency peaks, while lower electrophoretic mobility analytes exhibit lower frequency peaks. Based on obtained data, a migration velocity-adaptive moving average algorithm was created, adapted and programmed into capillary electrophoresis data processing software. Employing the developed algorithm, each data point is processed depending on a certain migration time of the analyte. Because of the implemented migration velocityadaptive moving average method the signal-to-noise ratio improved up to 11 times for sampling frequency of 4.6 Hz and up to 22 times for sampling frequency of 25 Hz. This paper could potentially be used as a methodological guideline for the development of new smoothing algorithms that require adaptive conditions in capillary electrophoresis, and other separation methods.

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 23

Introduction In the era of high speed and high-resolution data acquisition, there is a need for novel yet simple and efficient computational algorithms.1,2 The good example of such modernization is the last generation of analog-to-digital converters3,4. Affordable price, quality and size ratio allows to integrate such models (including field programmable gate array-based converters) as the main data acquisition modules5 not only into the expensive commercial laboratory equipment, but into the portable analytical systems as well. This is important for rapid microseparation methods, where little to none of the peak in the chromatogram/electropherogram should be lost in the filter. Whereas filtering of chromatographic, or electrophoretic data, in real-time requires significant computational power.6 Unfortunately, the microprocessors that are responsible for the computation in autonomous and portable instruments are not capable of handling high amount of data with computationally intensive algorithms. Therefore, current in situ data analysis systems are restricted only to low-speed data acquisition.7,8 The current trends in analytical method and instrumentation development includes miniaturization and performance.9 Miniaturization of the instrument could decrease the sensitivity.10 However, miniaturized techniques offer numerous advantages: miniaturized instruments (i) can be easily integrated into other hardware, (ii) can be coupled with other methods, so that in-line and on-line techniques can be realized, (iii) and can be used for in situ analysis in portable and autonomous instrumentation.7,11 The design of such instrumentation is not a trivial task due to the requirements for operation in a harsh and changing environment.12 The capillary electrophoresis (CE) separation method can be used for field-portable analytical investigations.13 Electrophoretic techniques offer several advantages over chromatographic techniques: they are (i) compact14, (ii) require fewer moving parts,15,16 (iii) are less susceptible to column/ capillary clogging17 and (iv) do not require complex sample preparation.18,19 Nevertheless, electrophoretic techniques are considered less sensitive (up to 2 orders of magnitude) compared to the chromatographic techniques, which pose problems detecting low amounts of chemicals.20,21 This is due to a minute dimension of a detection cell, which is mostly determined by an inner diameter of the separation capillary (25-75µm). There are three main strategies of increasing the sensitivity of the miniaturized and conventional separation instrumentation: (i) hardware improvement (such as Z or bubble cells for UV/VIS detection, higher excitation voltage used for contactless conductivity detection), (ii) software improvement and (iii) analysis condition adjustment.22,23 Signal conditioning using special software with previously developed algorithms is not a new approach; it has already been used in the field of analytical chemistry for several decades20,24,25. Software improvements mainly deal with signal compensation and signal filtering techniques26. There are several signal conditioning methods used with separation techniques: (i) low pass filtering (that includes hardware analog operational amplifier filtering, switched capacitor filtering, digital filtering, fast Fourier transform filtering), (ii) moving average and its variations (simple moving average filtering, weighted average filtering) and (iii) more complex mathematical algorithms such as Savitzky-Golay, or a finite impulse response method that requires more computer processing power for signal recording and processing than simple moving average method.24,26,27 In conventional separation systems, the detector signals are oversampled: the data is acquired at the higher rate than is required to record the peaks; therefore, noise is visible in the electropherograms. Random process theory implies that the information of peaks in the electropherograms is contained in the low-frequency region of the frequency spectrum, and the noise is encoded in the high-frequency region. This work clarifies the optimal smoothing

ACS Paragon Plus Environment

Page 3 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

parameters for each peak in the electropherogram, and the development of a smoothing algorithm that reduces the noise at each peak, depending on the optimal smoothing parameters. In this work, contactless conductivity detection (C4D), was selected to test the model due to its ease of use and modification possibilities.28,29 This technique is modifiable opposed to instrumentation supplied by official suppliers, which is closed-source and, unmodifiable since the algorithms are pre-programmed and restricted to a certain platform.30 The method described in this paper utilizes a posteriori information: (i) electrophoretic velocity and (i) frequency characteristics of the electropherograms. The method is simple, easily implemented, computationally efficient, and applicable to portable and stand-alone instrumentation in real-time and post-processing applications. The aim of this work was to design and test a migration velocity-adaptive moving average method that improves the sensitivity of acquired data in capillary electrophoresis.

Experimental Section Chemicals Acetic acid (99.8%) was purchased from Reachem (Slovakia). Sodium chloride (99.5%), potassium chloride (>99%) and 2-amino-2-methyl-1-propanol (AMP) (90%) were from Sigma-Aldrich (Germany). Tris(hydroxymethyl)aminomethane (TRIS) (99.8%) was from Merck (Germany). Bidistilled water was produced in the laboratory using Fistreem Cyclon bidistillator (UK).

Instrumentation The capillary electrophoresis system was donated by Agilent Technologies (Germany). Detection: previously optimized wireless, capacitance-to-digital technology based on differential measurement capability and, contactless conductivity detection system was used.31,32 Detection was performed using 32 kHz square wave signal with an amplitude of 3.3 V. Each electrode width was 24 mm, and the detection gap was 0.3 mm. The detector was completely shielded from the outer environment. The detector was connected to an Arduino Nano R3 microcontroller board, and data was sent using two-wire interface (TWI). The microcontroller board pre-processed the data so that values corresponding to picofarads (pF), or femtofarads (fF) were obtained and sent to the computer via transceiver module at the sampling frequencies of 4.6 Hz, 13 Hz, and 25 Hz. The function, named as “CAPCHOP”, which increases the sensitivity of the detector and in the specification of integrated circuit, was switched on during the operation.33 UV detection was used for comparative reasons.

Workflow of the research The research work was comprised of 4 stages (supporting information Figure S1): (i) capillary electrophoresis separations, (ii) signal analysis, (iii) design of migration velocity adaptive moving average method (algorithm) and (iv) ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 23

testing. First, the capillary electrophoresis procedure was performed and electropherogram recorded. The frequency characteristics of the peaks were determined, the cutoff frequency was found and the relationship between peak migration time and cutoff frequency was clarified by fitting the logarithmic equation. The migration velocity adaptive moving average algorithm was designed as follows: (a) the averaging window was calculated for each data-point of the electropherogram according to the logarithmic equation obtained, (b) the moving average was calculated for each point in the electropherogram and (c) the data processing software was upgraded with the new algorithm. Last, the developed algorithm was tested and processed, and unprocessed data was compared, showing significant improvement of signal-to-noise ratio (S/N ratio) and limits of quantification.

Sample preparation For the preparation of the stock solution, an adequate amount of chemical compound was dissolved in bidistilled water and frozen at -20 °C. For further use, solutions were defrosted and mixed together in adequate ratios. The obtained mixture was diluted, and 8 levels of different concentrations were obtained in the range of 0.23 – 417 μM, which were used for calibration and determination of detection and quantification limits.

Capillary electrophoresis procedure The procedure for the capillary electrophoretic analysis of cation was optimized and indirect mode was applied to detect negative peaks19. Analyses were performed in 0.5 M acetic acid background electrolyte (BGE), where 13 kV voltage potential was applied in 40 °C. The separation capillary was 59 cm total length (Ltot), 47 cm effective length (Leff), 50 μm inner diameter (I.D.) and 365 μm outer diameter (O.D.). Prior to each analysis, the capillary was flushed with background electrolyte for 3 min and electro-conditioned with 13 kV for 5 min. Samples were injected hydrodynamically at 50 mbar for 30 s. Calibration was performed in accordance with method validation guidelines (ICHQ2R1). 8 levels of different concentrations were used for calibration. At least 3 repetitions were performed for determination of concentration. Precision was expressed providing a percentage of relative standard deviation (RSD (%)). For each calibration curve, (y = a * b; y – peak area, a – slope of the calibration curve, b - concentration) the coefficient of determination (R2) was calculated. The quantification limit (LOQ) was calculated in accordance with the equation indicated in supporting information (equation S1).

Electrophoresis signal analysis Electrophoresis signal is observed (Figure 1A): x(j)=xr(j) (j=1,…,J;r=1,…,R); sampling period Δ = 0.216 s; r = 1,..,R – the number of the local peak of electrophoresis signal; xr(j)(j=rb, rb+1,…,re); rb is the starting point of xr(j). re – the end point of the xr(j).

ACS Paragon Plus Environment

Page 5 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Figure 1. Representation of signal analysis using STADIA3 software (A) imported electropherogram signal into STADIA3 software: x(j) = xr(j) (j=1,…,J; r=1,…,R). Sampling period Δ = 0.216 s34. r – the number of the local peak of electrophoresis signal. xr(j) (j=rb, rb+1,…,re); rb is the starting point of xr(j). re – the end point of the xr(j). (B) Calculated estimation of spectral density function Sr(f) of the rth peak and Frcutoff – cutoff frequency of the rth peak

The beginning (rb) and the end (re) of the peaks in the electropherogram were manually selected according to chromatographic peak determination procedure fixing starting and ending points. Since the peak selection stage was considered critical; therefore, each peak was inspected manually. Each peak in the electropherogram has its own frequency characteristics; therefore, the optimal averaging window size for each peak in the electropherogram is different. In order to clarify frequency characteristics of the peaks, spectral density functions were calculated for each selected peak using STADIA3 software following the equation indicated in supporting information (equation S2 and S3).34 Estimation of covariance function was calculated using equation indicated in supporting information (equation S4). It must be mentioned that such calculations can be performed using conventional software such as Matlab, or Rstudio if equations S2, S3 and S4 are followed.

Model of smoothing method Calculated spectral density functions help to determine the cutoff frequency of the peak, which is a threshold value (Hz) between information-carrying and noise frequencies in STADIA3 software. The cutoff frequency for each analyzed peak (Frcutoff) was manually determined at the peak and baseline breakpoint (Figure 1B). The Frcutoff is required for calculation of optimal averaging window size for a particular peak in the electropherogram following equation (1):  =

 (1) 2

Where Wr – the averaging window size of the rth peak, Fs – sampling frequency (Hz) and Frcutoff – the cutoff frequency (Hz) of the rth peak. If averaging window size is calculated for each peak, it leads to the development of adaptive ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 23

averaging algorithm, where the migration velocity-adaptive moving average algorithm is defined by following statements: (i) during electrophoresis process, the C4D detector provides an array of values, where analysis time – t (t = Δ * j, Δ – sampling period (0.216 s), j – point number) and amplitude x(j) is measured; (ii) depending on the data point number j, averaging window size Wr(j) is calculated at each data point following equation (2):  () = − 3.12 ln   + 16.8 (2)  Where Wr(j) – averaging window size, Fs – sampling frequency, ln – natural logarithm, Leff – effective length of the separation capillary (m), j – point number (j = t / Δ, t – analysis time and Δ – sampling period (0.216 s)); (iii) at each point – j estimation of smoothed signal zr(j) is calculated from original electropherogram (amplitude values) – x(j) following equation (3) (Supporting information Figure S2): () 1   − = 2 ()

#$(%)&(

! ')*

"()#$(%)&' ; , = 1, 2 … , (); () = − (3.12 ln 

  + 16.8) (3) 

Where zr(j-Wr(j)/2) (r = 1, 2,…, R) – calculated estimation of smoothed signal, Wr(j) – averaging window size, x(j) – original signal (amplitude values), i – step, Fs – sampling frequency (Hz), ln – natural logarithm, Leff – effective length of the separation capillary (m), j – point number (j = t / Δ, t – analysis time and Δ – sampling period (0.216 s)). The result of this procedure is an array of values providing migration velocity adaptively averaged electropherogram.

Results and discussion Frequency characteristics of the peaks The electro-migration velocity of the analytes depend on the charge, ionic volume and viscosity of the background electrolyte.35 Analyses of the mixture of the model analytes and real samples were performed. It was observed that during steady-state CE analysis, each later peak is wider than the former, providing a higher peak width value (supporting information Figure S3). These observations comply with the theoretical background stating that analyte migration velocity is a vector sum of the electrophoretic velocity of the substance and electroosmotic velocity of the BGE.35,36 The same concentration zones of various analytes do not migrate at equal speed, resulting in the difference of time required to pass the detection gap and result in uneven peak widths.36 A similar observation of analyte elution profiles is observed in isocratic chromatographic separations due to chromatographic zone broadening. For the electrophoretic and chromatographic signal smoothing finite response filtering method, which is classified as a weighted averaging method, can be used.27 Literature states that finite response filtering outperforms moving average filtering technique, or Savitzky-Golay method, because the recorded signals are transported to frequency domain, analyzed and transported back to time domain filtered.27,37 Finite response filtering is appropriate for chromatographic techniques, if gradient elution programs are used, due to the fact that (i) the peak shape and frequency characteristics that need signal analysis are hardly predictable, and (ii) chromatography systems are connected to computers that can analyze acquired data and readily process it. In electrophoretic techniques, the peak shape and frequency characteristics depend on the migration velocity and can be predicted (as stated below), which makes the signal filtering algorithm simpler and computationally efficient

ACS Paragon Plus Environment

Page 7 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

without the signal analysis process. As electrophoretic separation results show (Supporting information Figure S-3), the lowest peak widths were found for potassium ions (15.4 s and 20.1 s), which have highest electrophoretic mobility in BGE of 0.5 M acetic acid. Peak widths of other inorganic cations ranged between 18.2 s and 20.1 s, and organic cations provided a higher peak (AMP – 40.0 s and TRIS – 46.2 s). The highest peak width was observed for the unidentified compound in a real sample (ion-exchanger bleeding products) (69.7 s). Strong relationship between peak width and migration time was observed providing correlation of 0.86 (R2 0.74). Oversampled signals provide electropherograms with noise. In order to clarify if the signal is oversampled, the frequency characteristics must be determined.38 The main parameter is the cutoff frequency, which is the threshold in the frequency spectrum between noise and electropherogram peak information. The spectral density function is a measure that is used to determine the cutoff frequency of the peaks. Spectral density functions were calculated for each peak and the cutoff frequency was determined (Supporting information Figure S-4). All the calculated spectral density functions were found to have a single peak in the lowfrequency region (Supporting information Figure S-4 B). The peak in the low-frequency region of the spectral density function was considered as significant frequency carrying the peak information. The estimates in spectral density function of the high-frequency region that were at the level of the baseline were considered as a noise carrying frequency. It was observed that each later peak provides lower cutoff frequency. The same procedure was repeated for other electropherograms investigated in this work. The maximum frequency that was significant in recording the peak was determined, as depicted in Figure 1B (and supporting information Figure S-4 B). It was noticed that a single dataset from a single analysis, where concentrations of the analytes are similar, providing a comparable peak area, results in a strong relationship between frequency and analyte migration parameters (Supporting information Figure S-5 A). For a single dataset of the real sample (ion exchanger bleeding products), the highest relationship between migration velocity (m/s) and cutoff frequency (Hz) that is significant in recording the information of the peak, was found fitting linear curve where R2 was 0.99. According to Nyquist–Shannon sampling theorem, in order to not lose the information, the sampling frequency must be at least 2 times higher than the cutoff frequency of the signal.39,40 If the cutoff frequency of the peak is known, then the minimum sampling frequency is calculated by multiplying the cutoff frequency by 2. From the minimum sampling frequency, the period can be calculated when multiplied by the sampling frequency of the data acquisition system, provides a value that corresponds to maximum averaging window size (the equation (1) in Model of smoothing method). The relationship in the same dataset (ion exchanger bleeding product peaks) between peak migration velocity (Leff / mt (m/s); Leff – effective length of the capillary, mt – migration time of the peak) and averaging window size was clarified. Strong correlation (Supporting information Figure S-5) was observed providing R2 0.99 of the fitted logarithmic curve.

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 8 of 23

Figure 2. Determination of cutoff frequency (A) and averaging window size (B) for peaks

Scatterplots were created for multiple datasets of 27 data points, containing real and modeled samples, separated and unseparated, symmetric and asymmetric peaks of high, medium and low concentrations covering more cases and profiles of electrophoretic separations (Figure 2 A and B). In this case, a more scattered relationship was determined providing R2 0.88 of linear trend-line fitting for a scatterplot in Figure 2 A (cutoff frequency versus against migration velocity). The logarithmic fitting represented in Figure 2 B (averaging window size plotted against migration velocity (m/s)) provided a high negative correlation (-0.92 (R2 0.92)) of multiple electrophoretic profiles and analytical cases. Findings indicated that all investigated peaks provided similar tendency, which was the following: lower electrophoretic mobility analytes provided lower cutoff frequency.

Migration velocity-adaptive moving average algorithm Using the techniques represented in Figure 2 B, the empirical equation was obtained that was used for calculation of optimal averaging window size (supporting information, equation (S-5)), which was modified to equation (2) (Model of smoothing method section). The equation (2) instructs that a peak, depending on its migration time, can be averaged using moving average function when certain sampling frequency is used in data acquisition system and the effective length of the capillary is known. The obtained trivial equation serves as a constraint in a migration velocity– adaptive moving average function and can be explained by equation (3) (Electrophoresis signal analysis section). Considering the obtained migration velocity-adaptive moving average equation, the function was programmed in the software and further analyses were performed with an upgraded version of the software. The algorithm (Supporting information Figure S-6) performs specific sequence of the functions. At each point in the electropherogram the optimal averaging window is calculated and moving average function, according to the optimal averaging window, is calculated for each point of the electropherogram. At the beginning the averaging window size is 1, meaning that no averaging is performed at the start. Later, after a defined amount of time has passed, the averaging window size is incremented and kept unchanged until further incrementing. Average values are calculated from the original signal data-points, according to the averaging window size. This technique does not change the discreteness level, or information density of the recorded data, and no points in the electropherogram are lost. The result of this procedure

ACS Paragon Plus Environment

Page 9 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

is a smoothing of the noise in the electropherograms with no alterations to the shape of the peak. The code of the algorithm is represented in supporting information Annex S-5. The algorithm was tested with simulated repetitive (at different moments of simulated analysis time) step unit functions (Figure 3 A, B, and C) and simulated repetitive Gaussian peaks (Figure 3 D, E and F) following the methodology described in the literature41. Additionally, simulated white noise filtering was investigated (Figure 3 G), Gaussian noise (Figure 3 H) and filtering of simulated electrophoretic spikes (mainly caused by air bubbles) (Figure 3 I).

Figure 3. Smoothing of simulated signals. A, B and C – repetitive step unit functions (baseline at 0 and step at 1.0), window size 10 s, sampling frequency 4.6 Hz. Step position at (A) 30 s, (B) 220.5 s, (C) 602.5 s. D, E and F – Gaussian functions (baseline at 0) a = 1.0, c = 1. Step position at (D) 30 s, (E) 220.5 s, (F) 602.5 s. (G) White noise smoothing (sampling frequency 4.6 Hz), (H) Gaussian noise (sampling frequency 4.6 Hz) smoothing, (I) smoothing of simulated spikes in the electropherogram

Each later step or Gaussian peak is smoothed using bigger averaging window size.

Improvement of sensitivity

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 23

The algorithm that uses original electropherogram and processes it into a smoothed signal was employed in three different ways: (i) real-time processing in Arduino nano microcontroller, (ii) real-time processing in data acquisition software Universal Acquisition and (iii) post-processing in data visualization software Viewer. All alternatives provided identical analytical results stating that this method can be used: (i) in real-time data acquisition in the case of data storage restrictions and (ii) as a post-processing method for later analysis if there are no data storage restrictions. The first alternative is beneficial for the in situ autonomous analytical systems, such as on rovers.7,12,42 The second alternative is beneficial in clinical/ pre-clinical studies, process and quality control and other analytical applications, where data storage is not a concern. When used, the algorithm significantly improves the signal to noise (S/N) ratio of the electropherograms (Figure 4). The improvement of the S/N ratio of peaks, which have higher electrophoretic mobility, was found to be smaller than for those peaks, which migrate slower. This is due to the logarithmic equation used for the determination of averaging window size. Since the peaks of higher electrophoretic migration velocity are narrower and encoded by higher frequency information than low electrophoretic mobility analytes, the optimal averaging window size is lower. This means, that just after the start of the electrophoretic separation, the moving average function must not be used; after some time has passed, it should be used with utmost care so that peak shape would not be altered. After the defined amount of time has passed, the averaging window size is incremented. For K+ ions, which are fundamentally one of the fastest electro-migrating cations (migrate at 283 s) at the current conditions (Fs 4.6 Hz, 13 kV Ltot 59 cm, Leff 47 cm, temperature 40 °C, BGE – 0.5 M acetic acid solution) used, the averaging window size is 16 and TRIS peak, which migrates at mt of 948 s has the averaging window size of 31 points. For varying sampling frequency, capillary length, voltage, or temperature and electroosmotic instability the averaging window sizes are different and they only depend on the apparent analyte velocity. As observed in Figure 4 A, B and C, baseline drift is visible in original and processed electropherograms; therefore, it must be stressed that developed algorithm has no effect on baseline drift. For the compensation of the baseline drift due to temperature fluctuations, a different algorithm has been developed.43 If a uniform averaging window size is used for the whole electropherogram, non-optimal signal filtering is obtained if: (i) the averaging window size is too small, then the noise is still visible in the electropherograms (as demonstrated for moving average, FFT and Savitzky-Golay filters (Figure 4 G, H and I)) and (ii) if the averaging window size is too high, moving average algorithm distorts peaks and corrupts the information (as demonstrated for moving average, FFT and Savitzky-Golay filters (Figure 4 D, E and F)). It was observed that the proposed algorithm filters the noise of the electropherograms well if a higher Fs is used (Figure 4 A, B, and C). The sampling frequency should be selected with care since signal performance degrades at a higher Fs. Particularly for the experimental setup described in this paper, the effective dynamic range (effective resolution) of the detector operated at Fs of 4.6 Hz is 20.9 bits (the range is encoded by 220.9 (1956712) steps, detector operated at Fs of 13 Hz provides effective resolution of 20.7 bit and a detector operated at Fs of 25 Hz provides effective resolution of 20.1 bits.33 The tradeoff between the signal-to-noise ratio over the sampling frequency and averaging should be considered and carefully calculated.

ACS Paragon Plus Environment

Page 11 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Figure 4. Visual comparison of original and differently smoothed signals. Electropherograms: (A) collected at a sampling frequency of 4.6 Hz; (B) at 13 Hz; (C) at 25 Hz. Peaks: 1 – K+, 2 – Ca2+, 3 – Na+, 4 – AMP and 5 – TRIS. CE Conditions: 13 kV, 40 °C, Ltot 59 cm, Leff 47 cm, I.D. 50 μm, O.D. 365 μm. Detection: contactless conductivity at 3.3 V square wave of 32 kHz. K+ peak smoothing (D, E, and F) compared to different filters: proposed adaptive averaging, conventional averaging with a window size of 7 and 31, FFT filter with 7 and 31 points and Savitzky-Golay filter with 7 (4 + 3) and 31 (16 + 15) points

Different signal smoothing methods (moving average with the window size of 7 and 31, FFT filter with the window size of 7 and 31 and Savitzky-Golay smoothing method with the window size of 7 and 31) were compared by calculating main peak parameters, including a UV absorbance peak. The electropherograms were smoothed and the parameters (peak height, width and baseline noise) of each peak were measured. The increase of S/N ratio, height change, and width change was calculated, comparing them to original unsmoothed peaks; this is represented in Table 1 (full details represented in supporting information Table S-1). The methods with higher averaging window sizes (W = 31) provided higher S/N gain; however, those methods also distorted the peaks, which are of higher migration velocity. The methods with an averaging window of 7 provided insignificant peak distortions; however, the S/N gain was relatively low. The adaptive moving average method provided insignificant peak distortions and optimal S/N ratio.

Table 1. Comparison of main peak parameters after data processing using different methods

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Smoothing method

Page 12 of 23

Height Width S/N change change increase (%) (%)

Analyte: K+ Original Adaptive averaging Averaging, W=7 Averaging, W=31 FFT, W=7 FFT, W=31 Savitzky-Golay, (4+3) Savitzky-Golay, (16+15) Analyte: TRIS Original Adaptive averaging Averaging, W=7 Averaging, W=31 FFT, W=7 FFT, W=31 Savitzky-Golay, (4+3) Savitzky-Golay, (16+15)

100.0 98.0 97.8 72.0 94.3 54.4 100.0 95.3

100.0 103.0 102.9 128.3 113.4 165.7 107.4 108.9

1.0 4.6 3.0 8.3 3.9 6.9 1.8 5.3

100.0 101.8 112.8 99.0 104.0 90.1 97.3 95.4

100.0 102.9 97.4 110.4 98.8 117.9 100.4 100.8

1.0 11.1 3.4 11.4 4.3 11.4 1.8 5.3

Calibration parameters for each model analyte were calculated (Table 2). The comparison of calibration data using different filtering methods is represented in supporting information Table S-2. The precision of the method used did not exceed 3 % of relative standard deviation. An Insignificant difference in the calibration data between the processed and unprocessed electropherograms was observed. For calibration curves of (i) K+ cations, it was found that the difference was < 0.6 %, (ii) Na+ cations the difference was < 0.1 %, (iii) AMP < 0.2 % and (iv) TRIS < 0.2 %. The difference of R2 for (i) K+ was 0.02 %, (ii) Na+ was 0.17 %, (iii) AMP was 0.014 % and for (iv) TRIS the difference was 0.02 %. These differences can be considered as insignificant and related only to peak integration repeatability. It must be stressed that in the analytical process (sampling, sample preparation, pipetting, sample injection and measurement) errors that are introduced are usually up to two orders higher than calibration curve differences obtained in this work44.

Table 2. Comparison of calibration parameters for different analytes of processed and unprocessed data

Data Analyte K+ Na+ AMP

Unprocessed Equation

R2

y = 0.000185x y = 0.000398x y = 0.001176x

0.9961 0.9948 0.9970

Processed

LOQ (μM) 0.223 0.140 0.107

Equation y = 0.000184x y = 0.000398x y = 0.001174x

LOQ (μM) 0.9964 0.013 0.9965 0.008 0.9969 0.003 R2

ACS Paragon Plus Environment

S/N increase (times) Fs 4.6 Hz

Fs 13 Hz

Fs 25 Hz

4.6 7.6 7.2

8.3 8.5 13.6

9.8 11.8 16.4

Page 13 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

TRIS

Analytical Chemistry

y = 0.001327x

0.9988

0.094

y = 0.001329x

0.9986 0.002

11.1

14.4

22.3

The integrated area of the peaks in the electropherograms are two-dimensional values, while the S/N ratio is a onedimensional value, therefore the increase of LOQ is higher than the increase of the S/N ratio. Such tendency is observed and represented in Table 2. The literature states that the increase of S/N ratio is a square root of averaging window size26. Empirical determination of the increase of the S/N ratio was found slightly higher than the theoretical calculations. This observation can be explained by the fact that electropherograms contain less spectral components than white noise signal, which is typically used in theoretical calculations. For peaks: K+, Na+, AMP, and TRIS, the S/N ratio increased. 4.6, 7.6, 7.2 and 11.1 times at Fs of 4.6 Hz and LOQ decreased 16.8, 16.9, 35.6 and 47.0 times at Fs of 4.6 Hz. The significant decrease of LOQ can be explained by the fact that the peaks, which electro-migrate at slower velocity have more data-points and, are therefore averaged using a higher number of averaging points. Although concentrations at improved LOQ were not investigated, concentration close to the LOQ of TRIS (5 nM) was analyzed using capillary electrophoresis and is represented in supporting information Figure S-7. Although the adaptive smoothing approach has been used in the past for several applications, it has not been applied for capillary electrophoresis or similar separation method since the allowable size of the averaging window for different segments in electrophoretic data is unknown. This research work provides an answer to the main question raised: what is the allowable averaging window for a peak (actual and theoretical) in the electrophoretic data series? Knowing the allowable averaging window for different peaks in the electropherograms’ adaptive (or segmented) smoothing can be applied, resulting in significant sensitivity improvement. In many separations situations exist, when peaks are unseparated indicating overlapping. The averaging method is not directly targeted for resolution of overlapping peaks. However, if overlapped peaks are of a different cutoff frequency, they will be averaged using a different averaging window size, depending on their migration velocity. If overlapped peaks of a similar migration velocity fall into the same averaging window, they will be averaged using the same averaging window size. The developed algorithm is intended for use with portable and autonomous instrumentation, which normally have limited computational capabilities and cannot handle high speed and high-resolution data. In addition to this, the proposed algorithm is suitable not only for C4D detector data, but also for UV (supporting information table S-1 and Figure S-8) and other type detector data or separation methods. As represented in the figure, minor deviations from the original data was observed (peak no. 3 - phospholipase A2 glycosylated forms). This can be explained by the fact that the method, optimized for C4D detector and original UV detector, has a different detection window size. Moreover, the proposed methodology should not be restricted to adaptive FFT low pass, or Savitzky-Golay filtering methods. This paper can be used as a methodological guideline for the development of new smoothing algorithms that require adaptive conditions in capillary electrophoresis, or similar peak migration tendency providing separation methods.

Conclusions The migration velocity adaptive moving average algorithm, which significantly improves the sensitivity of detectors and the electropherograms, was developed. The algorithm improves the signal-to-noise ratio from 4.6 to 11.1 times, and the limit of quantification from 16 to 47 times, a sampling frequency is 4.6 Hz. The algorithm improves signal-toACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 23

noise ratio from 9.8 to 22.3 times, when the sampling frequency is 25 Hz. The algorithm is intended for use of the detector data processing in a portable or autonomous capillary electrophoresis systems.

Acknowledgement This research was funded by a grant (Nr. 09.3.3-LMT-K-712-02-0202) from the Research Council of Lithuania. The authors also want to express the appreciation to: (i) Agilent Technologies (Dr. Gerard Rozing) for donating capillary electrophoresis system HP3DCE. We would like to express our gratitude to Shadeh Ferris-Francis and Nathan Kovarik (CalTech, NASA, JPL, Microdevices and Sensor Systems, Advanced optical and electro-mechanical microsystems). The authors dedicate this work to the Centennial of the restored Baltic States: Estonia, Latvia and Lithuania.

ACS Paragon Plus Environment

Page 15 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

References (1)

Yang, L. L.; Lee, H. Y.; Wang, M. K.; Lin, X. Y.; Hsu, K. H.; Chang, Y. R.; Fann, W.; White, J. D. Real-Time Data Acquisition Incorporating High-Speed Software Correlator for Single-Molecule Spectroscopy. J. Microsc. 2009, 234 (3), 302–310.

(2)

Wojenski, A.; Pozniak, K.; Kasprowicz, G.; Zabolotny, W.; Byszuk, A.; Zienkiewicz, P.; Chernyshova, M.; Czarski, T. Fast Data Acquisition Measurement System for Plasma Diagnostics Using GEM Detectors. In First EPs Conference on Plasma Diagnostics - 1st ECPD; PoS - Proceedings of Science: Rome, 2015; pp 1–10.

(3)

AnalogDevices. LTC2500-32 - 32-Bit Oversampling ADC with Configurable Digital Filter. 2018, pp 1–54.

(4)

Texas Instruments. ADS126x 32-Bit, Precision, 38-kSPS, Analog-to-Digital Converrter (ADC) with Programmable Gain Amplifier (PGA) and Voltage Reference. 2015, pp 1–133.

(5)

Zheng, W.; Liu, R.; Zhang, M.; Zhuang, G.; Yuan, T. Design of FPGA Based High-Speed Data Acquisition and Real-Time Data Processing System on J-TEXT Tokamak. In Fusion Engineering and Design; North-Holland, 2014; Vol. 89, pp 698–701.

(6)

Matos, J. T. V; Duarte, R. M. B. O.; Duarte, A. C. Trends in Data Processing of Comprehensive Two-Dimensional Chromatography: State of the Art. J. Chromatogr. B Anal. Technol. Biomed. Life Sci. 2012, 910, 31–45.

(7)

Elkin, K. R. Portable, Fully Autonomous, Ion Chromatography System for on-Site Analyses. J. Chromatogr. A 2014, 1352, 38–45.

(8)

Sharma, S.; Tolley, L. T.; Tolley, H. D.; Plistil, A.; Stearns, S. D.; Lee, M. L. Hand-Portable Liquid Chromatographic Instrumentation. J. Chromatogr. A 2015, 1421, 38–47.

(9)

Ríos, Á.; Zougagh, M.; Avila, M. Miniaturization through Lab-on-a-Chip: Utopia or Reality for Routine Laboratories? A Review. Anal. Chim. Acta 2012, 740, 1–11.

(10)

Koczka, P. I.; Bodoki, E.; Gáspár, A. Application of Capacitively Coupled Contactless Conductivity as an External Detector for Zone Electrophoresis in Poly(dimethylsiloxane) Chips. Electrophoresis 2016, 37 (3), 398–405.

(11)

Campos, C. D. M.; da Silva, J. a. F. Applications of Autonomous Microfluidic Systems in Environmental Monitoring. RSC Adv. 2013, 3 (40), 18216.

(12)

Kaljurand, M. Analytical Separations in Harsh Environments. Trends Environ. Anal. Chem. 2014, 1, e2–e7.

(13)

Gorbatsova, J.; Jaanus, M.; Vaher, M.; Kaljurand, M. Digital Microfluidics Platform for Interfacing Solid-Liquid Extraction Column with Portable Capillary Electropherograph for Analysis of Soil Amino Acids. Electrophoresis 2016, 37 (3), 472–475.

(14)

Morbioli, G. G.; Mazzu-Nascimento, T.; Aquino, A.; Cervantes, C.; Carrilho, E. Recombinant Drugs-on-a-Chip: The Usage of Capillary Electrophoresis and Trends in Miniaturized Systems – A Review. Anal. Chim. Acta 2016, 935, 44–57.

(15)

Kidd, R. D.; Noell, A.; Kazarians, G.; Aubrey, A. D.; Scianmarello, N.; Tai, Y.-C. Ion Chromatography-on-a-Chip for Water Quality Analysis. Int. Conf. Environ. Syst. 2015, 45 (July), 1–8.

(16)

Kiplagat, I. K.; Kubáň, P.; Pelcová, P.; Kubáň, V. Portable, Lightweight, Low Power, Ion Chromatographic System with Open Tubular Capillary Columns. J. Chromatogr. A 2010, 1217 (31), 5116–5123.

(17)

Nazario, C. E. D.; Silva, M. R.; Franco, M. S.; Lanças, F. M. Evolution in Miniaturized Column Liquid Chromatography Instrumentation and Applications: An Overview. J. Chromatogr. A 2015, 1421, 18–37. ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 23

(18)

Özer, N.; Aksoy, Y.; Ögüs, I. H. New Sample Preparation Method for the Capillary Electrophoretic Determination of Adenylate Energy Charge in Human Erythrocytes. J. Biochem. Biophys. Methods 2000, 45 (2), 141–146.

(19)

Drevinskas, T.; Bartkuvienė, V.; Maruška, A. Electrophoretic Separation of Pressurized Hot Water Extracts of Hibiscus Sabdariffa L . Using Vitamin C as Background Electrolyte and Contactless Conductivity. Chemija 2014, 25 (4), 206–212.

(20)

Galievsky, V. A.; Stasheuski, A. S.; Krylov, S. N. “Getting the Best Sensitivity from on-Capillary Fluorescence Detection in Capillary Electrophoresis” – A Tutorial. Anal. Chim. Acta 2016, 935, 58–81.

(21)

Bizzotto, C. S.; Meinhart, A. D.; Ballus, C. A.; Ghiselli, G.; Godoy, H. T. Comparison of Capillary Electrophoresis and High Performance Liquid Chromatography Methods for Caffeine Determination in Decaffeinated Coffee. Food Sci. Technol. 2013, 33 (1), 186–191.

(22)

Hempel, G. Strategies to Improve the Sensitivity in Capillary Electrophoresis for the Analysis of Drugs in Biological Fluids. Electrophoresis 2000, 21 (4), 691–698.

(23)

Locatelli, M.; Melucci, D.; Carlucci, G.; Locatelli, C. Recent HPLC Strategies to Improve Sensitivity and Selectivity for the Analysis of Complex Matrices. Instrum. Sci. Technol. 2012, 40 (2–3), 122–137.

(24)

Savitzky, A.; Golay, M. J. E. Smoothing and Differentiation of Data by Simplified Least Squares Procedures. Anal. Chem. 1964, 36 (8), 1627–1639.

(25)

Gorry, P. A. General Least-Squares Smoothing and Differentiation by the Convolution (Savitzky-Golay) Method. Anal. Chem. 1990, 62 (6), 570–573.

(26)

Strauss, W. Digital Signal Processing, 2nd ed.; California Technical Publishing, 2000; Vol. 17.

(27)

Jackson, K. D.; Walton, S. J.; Campbell, D. Signal Enhancement Techniques for Chromatography Detection Systems. J. Automat. Chem. 1997, 19 (5), 145–152.

(28)

Kubáň, P.; Hauser, P. C. Contactless Conductivity Detection for Analytical Techniques— Developments from 2014 to 2016. Electrophoresis 2017, 38 (1), 95–114.

(29)

Kubáň, P.; Hauser, P. C. Ten Years of Axial Capacitively Coupled Contactless Conductivity Detection for CZE - A Review. Electrophoresis 2009, 30 (1), 176–188.

(30)

Wenig, P.; Odermatt, J. OpenChrom: A Cross-Platform Open Source Software for the Mass Spectrometric Analysis of Chromatographic Data. BMC Bioinformatics 2010, 11 (1), 405.

(31)

Drevinskas, T.; Kaljurand, M.; Maruška, A. Capacitance-to-Digital: A Single Chip Detector for Capillary Electrophoresis. Electrophoresis 2014, 35 (16), 2401–2407.

(32)

Drevinskas, T.; Maruška, A.; Briedis, V. Capacitance-to-Digital: The Upgrades of Single Chip Detector. Electrophoresis 2015, 36 (2), 292–297.

(33)

AnalogDevices. AD7745/AD7746 24-Bit Capacitance-to-Digital Converter with Temperature Sensor. Data Sheet. Analog Devices: Norwood, MA 2007, pp 1–28.

(34)

Kaukėnas, J.; Navickas, G.; Telksnys, L. Human-Computer Audiovisual Interface. Inf. Technol. Control 2006, 35 (2), 87–93.

(35)

Whatley, H. Clinical and Forensic Applications of Capillary Electrophoresis. Clin. Forensic Appl. Capill. Electrophor. 2001, 21–58. ACS Paragon Plus Environment

Page 17 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

(36)

Huang, X.; Coleman, W. F.; Zare, R. N. Analysis of Factors Causing Peak Broadening in Capillary Zone Electrophoresis. J. Chromatogr. A 1989, 480 (C), 95–110.

(37)

Fellinger, A. Data Analysis and Signal Processing in Chromatography; Elsevier, 1992.

(38)

Grimmett, G.; Stirzaker, D. Probability and Random Processes; Oxford University Press, 1989.

(39)

National Instruments. Data Acquisition ( DAQ ) Fundamentals; 1999.

(40)

Texas Instruments (Inc., B.-B. Principles of Data Acquisition and Conversion; 1994.

(41)

Wahab, M. F.; Dasgupta, P. K.; Kadjo, A. F.; Armstrong, D. W. Sampling Frequency, Response Times and Embedded Signal Filtration in Fast, High Efficiency Liquid Chromatography: A Tutorial. Anal. Chim. Acta 2016, 907, 31–44.

(42)

Yang, B.; Zhang, M.; Kanyanee, T.; Stamos, B. N.; Dasgupta, P. K. An Open Tubular Ion Chromatograph. Anal. Chem. 2014, 86 (23), 11554–11561.

(43)

Submitted_For_Publication_Drevinskas, T.; Telksnys, L.; Maruška, A.; Gorbatsova, J.; Kaljurand, M. Compensation of the Baseline Temperature Fluctuations for Autonomous CE-C4D Instrument Working in Harsh Environments. Electrophoresis 2018, 1–8.

(44)

Drevinskas, T.; Mickiene, R.; Maruška, A.; Stankevičius, M.; Tiso, N.; Mikašauskaite, J.; Ragažinskiene, O.; Levišauskas, D.; Bartkuviene, V.; Snieškiene, V.; et al. Downscaling the in Vitro Test of Fungal Bioremediation of Polycyclic Aromatic Hydrocarbons: Methodological Approach. Anal. Bioanal. Chem. 2016, 408 (4), 1043–1053.

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

For TOC only

ACS Paragon Plus Environment

Page 18 of 23

Page 19 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Figure 1. Representation of signal analysis using STADIA3 software (A) imported electropherogram signal into STADIA3 software: x(j) = xr(j) (j=1,…,J; r=1,…,R). Sampling period ∆ = 0.216 s34. r – the number of the local peak of electrophoresis signal. xr(j) (j=rb, rb+1,…,re); rb is the starting point of xr(j). re – the end point of the xr(j). (B) Calculated estimation of spectral density function Sr(f) of the rth peak and Frcutoff – cutoff frequency of the rth peak 70x28mm (300 x 300 DPI)

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2. Determination of cutoff frequency (A) and averaging window size (B) for peaks 59x21mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 20 of 23

Page 21 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

Figure 3. Smoothing of simulated signals. A, B and C – repetitive step unit functions (baseline at 0 and step at 1.0), window size 10 s, sampling frequency 4.6 Hz. Step position at (A) 30 s, (B) 220.5 s, (C) 602.5 s. D, E and F – Gaussian functions (baseline at 0) a = 1.0, c = 1. Step position at (D) 30 s, (E) 220.5 s, (F) 602.5 s. (G) White noise smoothing (sampling frequency 4.6 Hz), (H) Gaussian noise (sampling frequency 4.6 Hz) smoothing, (I) smoothing of simulated spikes in the electropherogram 108x69mm (300 x 300 DPI)

ACS Paragon Plus Environment

Analytical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 4. Visual comparison of original and differently smoothed signals. Electropherograms: (A) collected at a sampling frequency of 4.6 Hz; (B) at 13 Hz; (C) at 25 Hz. Peaks: 1 – K+, 2 – Ca2+, 3 – Na+, 4 – AMP and 5 – TRIS. CE Conditions: 13 kV, 40 °C, Ltot 59 cm, Leff 47 cm, I.D. 50 µm, O.D. 365 µm. Detection: contactless conductivity at 3.3 V square wave of 32 kHz. K+ peak smoothing (D, E, and F) compared to different filters: proposed adaptive averaging, conventional averaging with a window size of 7 and 31, FFT filter with 7 and 31 points and Savitzky-Golay filter with 7 (4 + 3) and 31 (16 + 15) points 101x59mm (300 x 300 DPI)

ACS Paragon Plus Environment

Page 22 of 23

Page 23 of 23 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Analytical Chemistry

TOC 42x22mm (600 x 600 DPI)

ACS Paragon Plus Environment