Research article Special Issues

ToMI-FBA: A genome-scale metabolic flux based algorithm to select optimum hosts and media formulations for expressing pathways of interest

  • The Total Membrane Influx constrained Flux Balance Analysis (ToMI-FBA) algorithm was developed in this research as a new tool to help researchers decide which microbial host and medium formulation are optimal for expressing a new metabolic pathway. ToMI-FBA relies on genome-scale metabolic flux modeling and a novel in silico cell membrane influx constraint that specifies the flux of atoms (not molecules) into the cell through all possible membrane transporters. The ToMI constraint is constructed through the addition of an extra row and column to the stoichiometric matrix of a genome-scale metabolic flux model. In this research, the mathematical formulation of the ToMI constraint is given along with four case studies that demonstrate its usefulness. In Case Study 1, ToMI-FBA returned an optimal culture medium formulation for the production of isobutanol from Bacillus subtilis. Significant levels of L-valine were recommended to optimize production, and this result has been observed experimentally. In Case Study 2, it is demonstrated how the carbon to nitrogen uptake ratio can be specified as an additional ToMI-FBA constraint. This was investigated for maximizing medium chain length polyhydroxyalkanoates (mcl-PHA) production from Pseudomonas putida KT2440. In Case Study 3, ToMI-FBA revealed a strategy of adding cellobiose as a means to increase ethanol selectivity during the stationary growth phase of Clostridium acetobutylicum ATCC 824. This strategy was also validated experimentally. Finally, in Case Study 4, B. subtilis was identified as a superior host to Escherichia coli, Saccharomyces cerevisiae, and Synechocystis PCC6803 for the production of artemisinate.

    Citation: Hadi Nazem-Bokaee, Ryan S. Senger. ToMI-FBA: A genome-scale metabolic flux based algorithm to select optimum hosts and media formulations for expressing pathways of interest[J]. AIMS Bioengineering, 2015, 2(4): 335-374. doi: 10.3934/bioeng.2015.4.335

    Related Papers:

    [1] Ruijin Hong, Jialin Ji, Chunxian Tao, Daohua Zhang, Dawei Zhang . Fabrication of Au/graphene oxide/Ag sandwich structure thin film and its tunable energetics and tailorable optical properties. AIMS Materials Science, 2017, 4(1): 223-230. doi: 10.3934/matersci.2017.1.223
    [2] Muhamed Shajudheen V P, Saravana Kumar S, Senthil Kumar V, Uma Maheswari A, Sivakumar M, Sreedevi R Mohan . Enhancement of anticorrosion properties of stainless steel 304L using nanostructured ZnO thin films. AIMS Materials Science, 2018, 5(5): 932-944. doi: 10.3934/matersci.2018.5.932
    [3] Ahmed Z. Abdullah, Adawiya J. Haider, Allaa A. Jabbar . Pure TiO2/PSi and TiO2@Ag/PSi structures as controllable sensor for toxic gases. AIMS Materials Science, 2022, 9(4): 522-533. doi: 10.3934/matersci.2022031
    [4] G. A. El-Awadi . Review of effective techniques for surface engineering material modification for a variety of applications. AIMS Materials Science, 2023, 10(4): 652-692. doi: 10.3934/matersci.2023037
    [5] Johana Gamez, Luis Reyes-Osorio, Oscar Zapata, Roberto Cabriales, Luis Lopez, Miguel Delgado-Pamanes . Study of protective hard coatings of SiO2-TiO2 on aluminum substrates. AIMS Materials Science, 2024, 11(2): 200-215. doi: 10.3934/matersci.2024011
    [6] Haridas D. Dhaygude, Surendra K. Shinde, Ninad B. Velhal, G.M. Lohar, Vijay J. Fulari . Synthesis and characterization of ZnO thin film by low cost modified SILAR technique. AIMS Materials Science, 2016, 3(2): 349-356. doi: 10.3934/matersci.2016.2.349
    [7] Neerajkumar Wayzode, Vinod Suryawanshi . Mechanical properties of graphene nanoplatelets reinforced glass/epoxy composites manufactured using resin film infusion process. AIMS Materials Science, 2023, 10(4): 693-709. doi: 10.3934/matersci.2023038
    [8] Lyudmyla V. Borkovska, Olexander Gudymenko, Olexander Stroyuk, Alexandra Raevskaya, Olena Fesenko, Tetyana Kryshtab . Photoinduced transformations of optical properties of CdSe and Ag-In-S nanocrystals embedded in the films of polyvinyl alcohol. AIMS Materials Science, 2016, 3(2): 658-668. doi: 10.3934/matersci.2016.2.658
    [9] Anshu Mli Gaur, Dinesh Singh Rana . Investigation of structural, dielectric and sensing properties of MgCl2/PVDF composite films prepared via solution casting technique. AIMS Materials Science, 2016, 3(3): 1117-1124. doi: 10.3934/matersci.2016.3.1117
    [10] Alain Mauger, Christian M. Julien . V2O5 thin films for energy storage and conversion. AIMS Materials Science, 2018, 5(3): 349-401. doi: 10.3934/matersci.2018.3.349
  • The Total Membrane Influx constrained Flux Balance Analysis (ToMI-FBA) algorithm was developed in this research as a new tool to help researchers decide which microbial host and medium formulation are optimal for expressing a new metabolic pathway. ToMI-FBA relies on genome-scale metabolic flux modeling and a novel in silico cell membrane influx constraint that specifies the flux of atoms (not molecules) into the cell through all possible membrane transporters. The ToMI constraint is constructed through the addition of an extra row and column to the stoichiometric matrix of a genome-scale metabolic flux model. In this research, the mathematical formulation of the ToMI constraint is given along with four case studies that demonstrate its usefulness. In Case Study 1, ToMI-FBA returned an optimal culture medium formulation for the production of isobutanol from Bacillus subtilis. Significant levels of L-valine were recommended to optimize production, and this result has been observed experimentally. In Case Study 2, it is demonstrated how the carbon to nitrogen uptake ratio can be specified as an additional ToMI-FBA constraint. This was investigated for maximizing medium chain length polyhydroxyalkanoates (mcl-PHA) production from Pseudomonas putida KT2440. In Case Study 3, ToMI-FBA revealed a strategy of adding cellobiose as a means to increase ethanol selectivity during the stationary growth phase of Clostridium acetobutylicum ATCC 824. This strategy was also validated experimentally. Finally, in Case Study 4, B. subtilis was identified as a superior host to Escherichia coli, Saccharomyces cerevisiae, and Synechocystis PCC6803 for the production of artemisinate.


    1. Introduction

    One of the challenges in the chemistry of materials is to achieve control over the structuration of molecules to enhance devices performance or driven to new applications [1,2,3]. In this context, supramolecular organization plays an important role in sensors performance. For instance, in the work of Volpati et al. [4] the role of the ultrathin film architecture composed by phthalocyanines in the sensing unit performance was exploited. That work revealed that a sensor array could be built using sensing units with identical or similar materials changing the architectures of the ultrathin films.

    Sensors based on nanostructured films are assembled from many different fabrication techniques on the nanoscale. Techniques known to produce films with nanostructured architectures in solid substrates are Langmuir, Langmuir-Blodgett, Langmuir-Schaefer, Layer-by-Layer, thermal evaporation between other [5,6,7,8]. These techniques allow the formation of thin films of organic molecules including extended π -conjugated ring systems, as phthalocyanines [1,9,10], and perylene derivatives [9,11,12,13], forming organized structures and superstructures driven by π-π-stacking interactions [1,9]. Both molecules, phthalocyanines, and perylenes are interesting to applications in electric and electrochemical sensors due to the noncovalent assemblies of π-conjugated chromophores owing to π-π interaction impacting the percolated pathways for charge transport [2].

    Based on the dependence of sensors performance in function of molecular organization and nanostructuration of films composed by phthalocyanine and perylene, we explore the possibility of generating films from iron phthalocyanine (FePc) and n-butylimidoethylenamine perylene (PTCD-NH2) with distinct electrical and electrochemical properties. The formation of the films was performed using Langmuir-Blodgett (LB) and physical vapor deposition (PVD) deposition techniques, while the properties and sensing application were evaluated using electrical impedance under interdigitated electrodes (IDE) (electrical properties), and differential pulse voltammetry (DPV) under ITO glass electrode (electrochemical properties). The sensorial performance was carried out in the presence of standard catechol solution, and the real application was evaluated in tea sample to polyphenol quantification.


    2. Materials and Methods


    2.1. PVD and LB films deposition

    The molecular structure of both materials used here, iron phthalocyanine (FePc) (Kodak) and n-butylimidoethylenamine perylene (PTCD-NH2) (provided by Dr. J. Duff from the Xerox Research Center of Canada) are shown in Figure 1. The PVD films of FePc and PTCD-NH2 were fabricated in a Boc Edwards vacuum system model Auto 306. The powders were placed into a metallic Ta boat. The depositions were carried out by an electrical current at 10−6 torr, being the thickness controlled by a quartz crystal microbalance.

    The LB films were grown in a KSV Langmuir trough model 2000. The Langmuir films of FePc were formed using a chloroform solution 0.5 mg/mL and deposited at a surface pressure of 25 mN/m, and PTCD-NH2 were formed using a trifluoroacetic/chloroform (1:10) solution with 0.15 mg/mL and deposited at a surface pressure of 30 mN/m. Both films were obtained by the symmetric compression of the barriers at 10 mm/min with ultrapure water in the subphase.

    Figure 1. representation of the molecular structure of iron phthalocyanine (FePc) and n-butylimidoethylenamine perylene (PTCD-NH2).

    2.2. Films characterization

    The role of film composition and nanostructuration on the impedance performance was evaluated using an assembly composed of four Pt interdigitated electrodes (IDE) covered with films produced with both materials, FePc and PTCD-NH2 using LB (5 layers) and PVD (10 nm) techniques. The approximate thickness of 5 LB layers is 9 nm for FePc and 15 nm for PTCD-NH2. The Pt IDE contains 50 pairs of digits, each digit presenting 10 μm in width, 5 mm in length and 100 nm in height, spaced 10 μm each other. The impedance spectroscopy was carried out using an impedance analyzer (Solartron 1260A). The frequency range considered was 80 Hz to 1 MHz and, applied signal with 50 mV of amplitude. The IDE modified units were immersed in ultrapure water and also in 1 μM catechol standard solution (Sigma -Aldrich Brazil).

    The role of film composition and nanostructuration on the electrochemical behavior was also evaluated using ITO electrode covered with films produced with FePc and PTCD-NH2 using LB (5 layers) and PVD (10 nm) techniques. The electrochemical behavior was carried out using a potentiostat/galvanostat μ-autolab type Ⅲ (EcoChimie). The electrochemical measurements were performed using the differential pulse voltammetry (DPV) applying a potential range of −1.0 to + 1.0 V vs. Ag/AgCl, 50 mV of pulse amplitude, and 10 mV/s of scan rate. The measurements were performed using a voltammetric cell of three electrodes, with Ag/AgCl saturated 3M KCl as a reference electrode, a platinum wire as a counter electrode and ITO covered with the PVD or LB films as working electrode. The voltammograms were obtained in supporting electrolyte (0.1 mol/L KCl) and also in 120 μM catechol standard solution.

    The roughness for all films deposited onto quartz substrate was estimated using atomic force microscopy (AFM). The measurements were performed out in a Nanosurf microscope model Easyscan2 in tapping mode and resolution of 1024 dots/line.


    2.3. Analytical curve

    The PVD PTCD-NH2 film was used to obtain the analytical curve for catechol detection. The analytical curves were obtained by impedance spectroscopy and DPV. The analytical curve from impedance spectroscopy was obtained through the addition of 1.0 × 10−9, 3.0 × 10−9, 6.0 × 10−9, 1.0 × 10−8, 3.0 × 10−8, 6.0 × 10−8, 1.0 × 10−7, 5.0 × 10−7 and 1.0 × 10−6 mol/L catechol standard solution, while from DPV was obtained through addition of 2 × 10−6, 4 × 10−6, 8 × 10−6, 1.2 × 10−5, 2 × 10−5, 4 × 10−5, 6 × 10−5, 8 × 10−5, 1.2 × 10−4, and 1.6 × 10−4 mol/L catechol standard solution. The electrical impedance technique (absence of supporting electrolyte) is more sensitive than electrochemical techniques. Thus, a high concentration of the analyte produces a saturation of the impedance-based sensor response (interdigitated electrodes). On the other hand, the electrochemical measurements show a response to catechol oxidation only above of 1 µmol/L.


    2.4. Real sample

    The PVD PTCD-NH2 film was applied as a sensor to detect the polyphenol content in mate tea sample. The tea sample was prepared by weighing 1.5 g of herb and adding 100 mL of ultrapure water, heating at 90 °C for 15 minutes with constant stirring (300 rpm). The mate tea sample was filtrated (hot filtration) and the volume adjusted to 100 mL with ultrapure water. The polyphenol determination by impedance spectroscopy and DPV were performed using a 0.002% and 1% (v/v) of the tea extract, respectively.


    2.5. Fittings parameters

    Equivalent electric circuits can model impedance curves. Z-View software package (Scribner Associates, Southern Pines, NC, USA) was employed for the equivalent circuit fitting. The models of equivalent circuit models for IDEs sensors illustrated in Figure 2 were applied for fitting the impedance spectra. The first equivalent circuit presented in Figure 2(a) was used for spectra fitting when only one semicircle appears in a spectrum at high frequencies, and the low frequency part of the Nyquist plot is quite linear. When the IDEs systems showed spectra with two semicircles, the second equivalent circuit was used to fit them, as presented in Figure 2(b). The parameter Rs represents the bulk solution resistance according to Randles Model. The parallel combination of the Cg-Rct represents the geometrical capacitance of the IDE sensor and the charge transfer resistance of the films layer. In series with the first RC parallel circuit, the second RC composed by a constant phase element, CPE, and R2 present the capacitance and resistance at the interface between the film and the solution [14]. The CPE was modelled as a non-ideal capacitor given by CPE = −1/(Ciω)α; where C is the capacitance which describes the charge separation at the double layer interface, ω is the frequency (rad/s) and α exponent associated with the heterogeneity of the electrode surface (0.5 < α < 1). In this circuit, the addition of CPE on the circuit is necessary, which can be ascribed to blocking of the diffusion process, once the impedance measurements are made in ultrapure water (absence of electrolyte). Thus, the charge separation occurs, and the interface acts as a non-ideal capacitor (CPE) [15].

    Fit precision is assessed through the chi-square parameter which is the sum of squares of the ratio of the standard deviation between the original data and the calculated spectrum. For all spectra showed in this work, the chi-square parameter was smaller than 5.6 × 10−4.

    Figure 2. electrical equivalent circuit fitted to IDE bare (a), and IDE covered with films (b). (c) Scheme illustrating each physic-chemical process included in the model.

    3. Results and Discussion


    3.1. Electrical impedance characterization

    The effect of the thin film composition and nanostructuration on the performance of sensors based on impedance spectroscopy was first evaluated in ultrapure water. Figure 3 shows the impedance results of the LB and PVD films for both materials and bare IDE immersed in ultrapure water regarding Bode (|Z| vs. frequency) and Nyquist (Z' vs. Z”) plot.

    The impedance spectra profiles of both materials are similar for PVD and LB films. Furthermore, the profile found here is similar to that found by Yang working with the detection of Salmonella cell suspensions in deionized water using gold interdigitated microelectrodes [16]. This impedance spectrum profile (Bode plot) is typically found in systems where the polarization is due to a combination of kinetic and diffusion processes [16]. Comparing the Bode plot for all films studied, the major differences are observed to the PTCD-NH2 film. The different impedance values for the same material in the form of PVD and LB films are related to the molecular organization of these materials when processed in different ways. For instance, as discussed in Volpati et al. [4], FePc presents an orientation almost parallel to the substrate (between 0° and 45°) in LB film and a preferred orientation perpendicular to the substrate (45° to 90°) in PVD film. For PTCD-NH2, the molecules are oriented at an angle of approximately 45° with the substrate in PVD films and rather more inclined (> 45o) when in LB films (unpublished results). The important point is that the conductivity is greatly affected by materials that present anisotropy in the molecular orientation. For instance, in the work of Dey and Pal [5], the grown of LbL films of NiTsPc under a magnetic field were studied. Their result shows that the magnetic field induces a preferential molecular orientation on the LbL films of NiTsPc and this change in orientation was responsible for the observed changes in conductivity values and dielectric constant films. According to Dey and Pal [5], films in which the molecules were immobilized with the macrocycle orientation parallel to the substrate had lower values of dielectric constant and conductivity. Note that in the case of Dey and Pal [5], the electrical devices have been made in Sandwich-Al ITO. This effect is explained by the authors due to the very close packing of the molecules to each other, and the process of electronic conduction occurring in the molecule to molecule is improved. The authors also suggest that when the molecules are oriented perpendicular to the substrate surface, and the macrocycle in a particular direction, the electronic conduction by hopping still prevails compared with films in which the molecules have no specific orientation or organization [5].

    Besides, all PVD and LB films promote an increase in the impedance for lower frequencies (~50 kHz), as seen in the Bode plot (Figure 3(a)). The order in the impedance values is LB/PTCD-NH2 > PVD/PTCD-NH2 > LB/FePc > PVD/FePc. In the low frequencies range (<100Hz), impedance related to the double layer (Cd) appear in the spectra and contributes to the impedance value in a frequency-dependent manner: Z = −j/√2πfC, where C is capacitance and f is frequency. In this way, the major Cd is presented by the FePc films. The double layer capacitance is attributed to the reversible ion adsorption onto the film surface, which gives rise to an electrolyte/film interface, where no faradaic (redox) processes are involved. This capacitance can be described according to C = εA/d, where ε is the electrolyte dielectric constant, A is the film surface area accessible to ions, and d is the distance between the center of the ion and film surface. Thus, once in our case, the same electrolyte was used, higher values of Cd could be related to two mainly factors: i) larger surface areas promoted by porous (or aggregated) films and ii) shorter distances between ion-surface promoted by the material characteristic.

    Figure 3. Impedance spectra of IDE bare and IDE covered with PVD and LB films immersed in ultrapure water using Bode plot impedance (a) and phase (b), Nyquist plot (c) and the magnification of the high-frequency range (d) from the plot (c). The solid lines represent fitted data using the equivalent circuits.

    All impedance spectra show a semicircle in the high frequency range (kinetic control of the charge-transfer process) and a quite linear range at a lower frequency (diffusion control) in Nyquist plots. The impedance spectrum of bare IDE electrode was modeled using the equivalent circuit in Figure 2(a). The impedance spectra of the IDEs coated with films were modeled using the equivalent circuit in Figure 2(b). In the presence of film on IDE surface, a resistance parameter was added in comparison with the equivalent circuit from the bare. The first RC parallel can be associated with the charge transfer equilibrium at IDE surface, once that the circuit parameters show no significant changes between the films (Table 1). However, the second RC parallel showed significant changes, suggesting an influence of the films on diffusion process. These differences in circuit parameters for the same material in the form of PVD and LB films are related to the molecular organization of these materials when processed in different ways. Moreover, in the FePc case, the roughness decreases from 10.27 nm to 1.4 nm from the LB to PVD films. This decrease of roughness can be ascribed to FePc characteristic in aggregates itself through π-π interaction in solution [17], which are transferred to Langmuir films. Thus, the PVD technique which uses solid compound produces more homogeneous surface during the film deposition. Based on molecular organization and surface roughness the PVD technique allows a molecular orientation that favors the insertion of water molecules on the FePc films, increasing the CPE values.

    Table 1. Parameters obtained from equivalent circuit fitting from impedance measurements in ultrapure water and roughness of the IDE bare and IDE covered with PVD and LB films.
    Film R /
    C/µF sα-1 Rct/
    CPE/µF sα-1 α R2/
    Roughness/ nm
    Bare 301.7 1.30 × 10−10 25122 1.15 0.65 -
    LB/FePc 117.4 1.35 × 10−10 26080 0.0828 0.70 4.19 × 105 10.27
    PVD/FePc 175.1 1.28 × 10−10 26141 0.267 0.68 1.00 × 106 1.4
    LB/PTCD-NH2 131.9 1.35 × 10−10 26471 0.0169 0.82 6.27 × 105 4.10
    PVD/PTCD-NH2 188.7 1.33 × 10−10 22363 0.0164 0.84 1.01 × 106 2.00
     | Show Table
    DownLoad: CSV

    The PTCD-NH2 films showed small variation not only on the molecular organization (between 45° and 90°) but also on roughness. The latter being from 4.10 nm to 2.00 nm for LB and PVD films, respectively. In agreement with the impedance results obtained to FePc films, the PTCD-NH2 showed the similar tendency, being the PDV films more homogeneous than LB films. Thus we can conclude that molecular organization of thin films on substrate surface influences directly on impedance properties in ultrapure water.

    After impedance analysis in ultrapure water, the nanostructured films measurements in standard catechol solution were done to investigate the influence of the different composition and molecular orientation of the thin films as impedance sensor. The IDEs modified with LB and PVD films were immersed in 1.0 × 10−6 mol/L catechol standard solution. All films showed decreases of a semicircle, ascribed to increases of charge transfer on the electrode surface (Figure 4). The results indicate the LB and PVD films from FePc and PTCD-NH2 showed interaction with catechol molecules. The equivalent circuit modeled in this situation was the same that applied to films immersed in ultrapure water. However, the Rct decreases as expected, which can be ascribed to the increase of charge transfer between the film and catechol molecules, increasing the conductivity of the system. The greatest decrease of the semicircle diameter (Rct) in catechol solution was observed to PVD PTCD-NH2 films, indicating the best sensitivity. This sensitivity is an interesting property to the sensor application. Thus, the PVD PTCD-NH2 sensor was used to evaluate the sensing performance in low catechol concentrations.

    Table 2. Parameters obtained from equivalent circuit fitting from impedance measurements in catechol solution of the IDE bare and IDE covered with PVD and LB films.
    Film R /
    C/µF sα-1 Rct/
    CPE/µF sα-1 α R2/
    LB/FePc 124.5 1.35 × 10−10 22322 0.0734 0.71 4.66 × 105
    PVD/FePc 183.5 1.27 × 10−10 23068 0.284 0.68 1.10 × 106
    LB/PTCD-NH2 136.3 1.34 × 10−10 22429 0.0162 0.83 5.87 × 105
    PVD/PTCD-NH2 189.4 1.33 × 10−10 19763 0.0168 0.84 8.99 × 105
     | Show Table
    DownLoad: CSV
    Figure 4. Nyquist plot of IDE with PVD and LB films immersed in ultrapure water (square curve) and 1.0 × 10−6 mol/L catechol solution (circle curve). Lines indicates equivalent circuit fitting. Inset: magnification of the high-frequency range.

    The impedance measurements were performed using catechol concentrations in ultrapure water: 1.0 × 10−9, 3.0 × 10−9, 6.0 × 10−9, 1.0 × 10−8, 3.0 × 10−8, 6.0 × 10−8, 1.0 × 10−7, 5.0 × 10−7 and 1.0 × 10−6 mol/L. The increases of catechol concentration promote increases of charge transfer observed through the decreases semicircle diameter (decreases of Rct). The curve of 1/Rct in function of catechol concentration showed an exponential dependence, as showed at Figure 5. Linearity can be observed from 1.0 × 10−9 to 1.0 × 10−6 mol/L, with linear regression of 1/Rct (1/Ω) = 2.09255 × 10−5 + 2.09002 × 10−6 [catechol].

    Figure 5. response of the PVD PTCD-NH2 impedance sensor in ultrapure water containing different catechol concentrations. Measurements performed to three different electrodes (average from measurements in triplicate).

    3.2. Differential pulse voltammetry

    The results from voltammetric measurements collaborate with the results observed from impedance in ultrapure water. The films formed by PVD technique showed better stability and voltammetric signal in supporting electrolyte (0.1 mol/L KCl). The redox potential observed in negative range (Figure 6) can be ascribed to ring oxidation of perylene or phthalocyanine ligand [18,19]. The electrodes immersed in catechol solution (Figure 6) showed a peak potential ascribed to catechol oxidation [20], which showed a shift of oxidation potential to lower potential in comparison with ITO electrode uncovered (catechol oxidation at 0.231 V vs. Ag/AgCl). The small catechol oxidation potential was observed to PVD films, being 0.046 and 0.076 V vs. Ag/AgCl for PTCD-NH2 and FePc films, respectively. The potentials for catechol oxidation are summarized in Table 3. The PDV films are more homogeneous than LB films, suggesting like as in impedance behavior, the homogeneity can influence on the electrocatalytic effect of this films. The increase of homogeneity and the decreases of roughness observed from PVD to LB films decrease the onset potential of catechol oxidation.

    Table 3. Potential of catechol oxidation using PTCD-NH2 and FePc films
    Film Eox (V) vs. Ag/AgCl
    Bare 0.231
    PVD/PTCD-NH2 0.046
    LB/PTCD-NH2 0.167
    PVD/FePc 0.076
    LB/FePc 0.195
     | Show Table
    DownLoad: CSV
    Figure 6. differentialpulse voltammograms of PVD and LB PTCD-NH2 films (a), and PVD and LB FePc films (b) immersed in supporting electrolyte (0.1 mol/L KCl) and standard catechol solution (1.2 × 10−4 mol/L). The voltammogram of ITO bare in catechol solution also is shown.

    Also, the thickness of PVD films showed influence on film oxidation but has no influence on catechol peak oxidation (supporting information). This behavior reinforces the hypothesis that the interaction of catechol occurs on the film surface, which promotes a fast saturation of films surface/signal principally observed on impedance measurements.

    The voltammetry in agreement with impedance results collaborates with the hypothesis that molecular organization of thin films on substrate surface influences directly on electrochemical and electrical properties of the films. Thus, the results from DPV and impedance measurements suggesting PTCD-NH2 films showed the better potential (i.e. smaller potential to catechol oxidation) for the development of sensor devices to the determination of catechol and/or polyphenols derivatives.

    The PVD PTCD-NH2 films were also used to evaluate the electrochemical sensing performance. The DPV measurements were performed using 2 × 10−6, 4 × 10−6, 8 × 10−6, 1.2 × 10−5, 2 × 10−5, 4 × 10−5, 6 × 10−5, 8 × 10−5, 1.2 × 10−4 and 1.6 × 10−4 mol/L catechol standard solution. The anodic peak current at 0.046 V vs. Ag/AgCl (ascribed to catechol oxidation) increase linearly with the catechol concentration from 2.0 × 10−6 to 1.6 × 10−4 mol/L (Figure 7). The linear regression equation can be represented by Ipa (A) = 1.956 × 10−7 + 0.2516 [catechol mol/L] (r = 0.998), with limit of detection (LOD) of 1.62 × 10−7 mol/L (0.162 μmol/L) calculated according to 3 x SD/slope criteria.

    Figure 7. response of the PVD PTCD-NH2 voltammetric sensor in KCl solution for different catechol concentrations. Measurements performed in triplicate using the same electrode, and three different electrodes as well (average of nine addition measurements).

    The analytical performance of PVD PTCD-NH2 was similar or better than other sensors and biosensor as summarized in Table 4. The PVD PTCD-NH2 films show a satisfactory limit of detection to the application as catechol sensor.

    Table 4. Comparison of analytical performance of sensors for the detection of catechol
    Film Detectiontechnique Linear range
    (μmol/L)
    LOD (μmol/L) Ref.
    Tyr/AA/LuPc2 Cyclic voltammetry 1.67-21.46 1.71 [20]
    (PAH/FeTsPc+ DPPG/AgNp)5 Cyclic voltammetry 2.0-100 0.87 [21]
    PEDOT/CPE DPV 1.0-250 0.50 [22]
    Polypyrrol/AuNp Cyclic voltammetry 0.1-1.0 3-88 [23]
    3-APBA-PTCA-CNTs/GC DPV 0.5-30 0.10 [24]
    PVD PTCD-NH2 DPV 2-160 0.16 This work
    Abbreviation: Tyr = tyrosine enzyme; AA = arachidic acid; LuPc2 = lutetium phthalocyanine; PAH = poly (allylamine) hydrochlorate; FeTsPc = iron tetrasulfonated phthalocya-nine; DPPG = phospholipid 1, 2-dipalmitoyl-sn-3-glycero-(phosphor-rac-(1-glycerol); AgNp = silver nanoparticles; PEDOT = poly(3, 4-ethylenedioxythiophene); CPE = carbon past electrode; AuNp = gold nanoparticles; APBA-PTCA = 3-aminophenylboronic acid-3, 4, 9, 10-perylene tetracarboxylic acid; CNTs = carbon nanotubes, GC = glassy carbon.
     | Show Table
    DownLoad: CSV

    3.3. Analyses in mate tea samples

    PVD PTCD-NH2 films were applied as a sensor to determine the catechol content in mate tea samples using electrical impedance and DPV methods. The catechol concentration in mate tea sample was determined using standard addition in both analytical methods. The catechol concentrations in mate tea samples determined by DPV measurements was 171.3 mg/L (±15.6), while the value determined using impedance measurements was 0.09 mg/L (± 0.13). Comparing with the value found using Folin-Cioalteu method (297 mg/L ± 8.5), a standard method for polyphenol quantification, DPV method showed higher efficiency than impedance method. The results for polyphenols quantification are summarized in Table 4.

    The impedance measurements are very sensitive to the variation on interface IDE surface/solution. Thus, the complex tea sample matrix showed a higher influence on charge transfer with leading to an error on catechol quantification. Besides, an adsorption of polyphenols and the other compounds present in tea matrix can occur. The poisoning or "memory effect" can be observed through comparison of the impedance response in the initial ultrapure water (before measurements in catechol) and final ultrapure water (after measurements in catechol) in Figure 8. As is shown below for PTCD-NH2, the impedance profile changes due to the catechol adsorption. This effect was not observed in the electrochemical measurements. Thus, in this case, the electrical impedance cannot be used as an analytical method using the PVD PTCD-NH2 films.

    Figure 8. Nyquist plot of IDE with PTCD-NH2 PVD film immersed in ultrapure water before (initial water) and after (final water) catechol measurements.

    However, using DPV method the concentration found was close to that determined by the standard method. The difference observed are assigned to oxidation potential used in DPV the measurements. The maximum peak potential for the mate tea sample using PDV PTCD-NH2 films was observed at 0.205 V vs. Ag/AgCl (1% of tea), and to catechol oxidation at 0.046 V vs. Ag/AgCl. The oxidation peak observed in the tea sample can be ascribed to oxidation of epicatechin flavan-3-ol derivatives [25]. The mate tea sample showed higher levels of EC than other flavan-3-ol. Thus, the values found in this work include only polyphenols oxidized around 0.046 V.

    Table 4. Polyphenol quantification in mate tea sample
    Method* Amount of polyphenol (mg/L)
    Differential pulse voltammetry 171.3 (±15.6)
    Electrical impedance 0.090 (±0.13)
    Folin-Ciocalteu 297.0 (±8.50)
    *the differential pulse voltammetry and electrical impedance were applied using PDV PTCD-NH2 films.
     | Show Table
    DownLoad: CSV

    4. Conclusion

    The Langmuir-Blodgett and physical vapor deposition techniques allow the formation of thin films from FePc and PTCD-NH2 compounds with different molecular organization and nanostructuration. Molecular organization and nanostructuration directly influenced the electrical and electrochemical properties of these films. The impedance results in ultrapure water showed the films onto IDE modified the electrical response. In general, equivalent circuits reveal an increase in the resistance of the system in the low frequency range with the films deposition, associated with the diffusional process. Comparing LB and PVD films, the PDV deposition promotes more homogeneous films and with lower roughness. These characteristics favor the insertion of water molecules, increasing the CPE values. The same tendency was observed in catechol solution, once PDV films better interact with catechol molecules. The latter was analyzed by the decreasing of resistance values in high and low frequency range. The electrocatalytic properties were also evaluated by using differential pulse voltammetry for films deposited onto ITO electrode. Both PVD films (FePc and PTCD-NH2) showed smaller onset potential for catechol oxidation in comparison with LB films, which can be directly associated with the molecular organization, homogeneity, and roughness. However, the thicknesses not show significant influence.

    These observations agree with impedance results, which confirm nanostructuration play an important role on electrical and electrochemical properties of this films (principally in catechol solution). Among FePc and PTCD-NH2 films (LB and PVD), the PVD PTCD-NH2 film showed interesting electrical and electrochemical properties for the application as catechol sensor. The analytical curve for catechol detection was obtained using PVD PTCD-NH2 film through impedance and DPV methods. The impedance signal showed an exponential dependence with low catechol concentration, indicating a sensitive interaction on IDE interface. However, the voltammetric signal was linear with catechol concentration with a limit of detection of 1.62x10-7 mol/L. The PVD PTCD-NH2 film was tested as a sensor for quantification of catechol contents in mate tea samples. Using DPV method the concentration of catechol found (171.3 mg/L) was close to that determined by the Folin-Cioalteu standard method.


    Acknowledgments

    This work was supported by FAPESP, CNPq, and CAPES.


    Conflict of Interest

    All authors declare no conflicts of interest in this paper.


    [1] Hanai T, Atsumi S, Liao JC (2007) Engineered synthetic pathway for isopropanol production in Escherichia coli. Appl Environ Microbiol 73: 7814-7818.
    [2] Atsumi S, Cann AF, Connor MR, et al. (2008) Metabolic engineering of Escherichia coli for 1-butanol production. Metabolic Eng 10: 305-311. doi: 10.1016/j.ymben.2007.08.003
    [3] Connor MR, Liao JC (2008) Engineering of an Escherichia coli strain for the production of 3-methyl-1-butanol. Appl Environ Microbiol 74: 5769-5775. doi: 10.1128/AEM.00468-08
    [4] McKenna R, Nielsen DR (2011) Styrene biosynthesis from glucose by engineered E. coli. Metabolic Eng 13: 544-554. doi: 10.1016/j.ymben.2011.06.005
    [5] McKenna R, Pugh S, Thompson B, et al. (2013) Microbial production of the aromatic building-blocks (S)-styrene oxide and (R)-1,2-phenylethanediol from renewable resources. Biotechnol J 8: 1465-1475.
    [6] Moon TS, Dueber JE, Shiue E, et al. (2010) Use of modular, synthetic scaffolds for improved production of glucaric acid in engineered E. coli. Metabolic Eng 12: 298-305. doi: 10.1016/j.ymben.2010.01.003
    [7] Tseng HC, Martin CH, Nielsen DR, et al. (2009) Metabolic engineering of Escherichia coli for enhanced production of (R)- and (S)-3-hydroxybutyrate. Appl Environ Microbiol 75: 3137-3145. doi: 10.1128/AEM.02667-08
    [8] Ajikumar PK, Xiao WH, Tyo KE, et al. (2010) Isoprenoid pathway optimization for Taxol precursor overproduction in Escherichia coli. Science 330: 70-74. doi: 10.1126/science.1191652
    [9] Engels B, Dahm P, Jennewein S (2008) Metabolic engineering of taxadiene biosynthesis in yeast as a first step towards Taxol (Paclitaxel) production. Metabolic Eng 10: 201-206. doi: 10.1016/j.ymben.2008.03.001
    [10] Scaife MA, Burja AM, Wright PC (2009) Characterization of cyanobacterial beta-carotene ketolase and hydroxylase genes in Escherichia coli, and their application for astaxanthin biosynthesis. Biotechnol Bioeng 103: 944-955. doi: 10.1002/bit.22330
    [11] Lemuth K, Steuer K, Albermann C (2011) Engineering of a plasmid-free Escherichia coli strain for improved in vivo biosynthesis of astaxanthin. Microb Cell Fact10: 29.
    [12] Tsuruta H, Paddon CJ, Eng D, et al. (2009) High-level production of amorpha-4,11-diene, a precursor of the antimalarial agent artemisinin, in Escherichia coli. PloS One 4: e4489. doi: 10.1371/journal.pone.0004489
    [13] Becker J, Wittmann C (2012) Bio-based production of chemicals, materials and fuels -Corynebacterium glutamicum as versatile cell factory. Curr Opin Biotechnol 23: 631-640.
    [14] Lindberg P, Park S, Melis A (2010) Engineering a platform for photosynthetic isoprene production in cyanobacteria, using Synechocystis as the model organism. Metabolic Eng 12: 70-79. doi: 10.1016/j.ymben.2009.10.001
    [15] Hatzimanikatis V, Li C, Ionita JA, et al. (2005) Exploring the diversity of complex metabolic networks. Bioinformatics 21: 1603-1609. doi: 10.1093/bioinformatics/bti213
    [16] Finley SD, Broadbelt LJ, Hatzimanikatis V (2009) Computational framework for predictive biodegradation. Biotechnol Bioeng 104: 1086-1097. doi: 10.1002/bit.22489
    [17] Henry CS, Broadbelt LJ, Hatzimanikatis V (2010) Discovery and analysis of novel metabolic pathways for the biosynthesis of industrial chemicals: 3-hydroxypropanoate. Biotechnol Bioeng 106: 462-473.
    [18] Fisher AK, Freedman BG, Bevan DR, et al. (2014) A review of metabolic and enzymatic engineering strategies for designing and optimizing performance of microbial cell factories. Comput Struct Biotechnol J 11: 91-99. doi: 10.1016/j.csbj.2014.08.010
    [19] Lee SK, Chou H, Ham TS, et al. (2008) Metabolic engineering of microorganisms for biofuels production: from bugs to synthetic biology to fuels. Curr Opin Biotechnol 19: 556-563.
    [20] McEwen JT, Atsumi S (2012) Alternative biofuel production in non-natural hosts. Curr Opin Biotechnol 23: 744-750.
    [21] Jang YS, Park JM, Choi S, et al. (2012) Engineering of microorganisms for the production of biofuels and perspectives based on systems metabolic engineering approaches. Biotechnol Adv 30: 989-1000. doi: 10.1016/j.biotechadv.2011.08.015
    [22] Atsumi S, Hanai T, Liao JC (2008) Non-fermentative pathways for synthesis of branched-chain higher alcohols as biofuels. Nature 451: 86-89. doi: 10.1038/nature06450
    [23] Atsumi S, Higashide W, Liao JC (2009) Direct photosynthetic recycling of carbon dioxide to isobutyraldehyde. Nat Biotechnol 27: 1177-1180. doi: 10.1038/nbt.1586
    [24] Smith KM, Cho KM, Liao JC (2010) Engineering Corynebacterium glutamicum for isobutanol production. Appl Microbiol Biotechnol 87: 1045-1055. doi: 10.1007/s00253-010-2522-6
    [25] Jia X, Li S, Xie S, et al. (2011) Engineering a metabolic pathway for isobutanol biosynthesis in Bacillus subtilis. Appl Biochem Biotechnol 168: 1-9.
    [26] Chen X, Nielsen KF, Borodina I, et al. (2011) Increased isobutanol production in Saccharomyces cerevisiae by overexpression of genes in valine metabolism. Biotechnol Biofuels 4: 21. doi: 10.1186/1754-6834-4-21
    [27] Higashide W, Li Y, Yang Y, et al. (2011) Metabolic engineering of Clostridium cellulolyticum for production of isobutanol from cellulose. Appl Environ Microbiol 77: 2727-2733.
    [28] Orth JD, Thiele I, Palsson BO (2010) What is flux balance analysis? Nat Biotechnol 28: 245-248. doi: 10.1038/nbt.1614
    [29] Varma A, Palsson BO (1994) Metabolic flux balancing - basic concepts, scientific and practical use. Bio-Technol 12: 994-998. doi: 10.1038/nbt1094-994
    [30] Reed JL, Palsson BO (2004) Genome-scale in silico models of E. coli have multiple equivalent phenotypic states: assessment of correlated reaction subsets that comprise network states. Genome Res 14: 1797-1805.
    [31] Henry CS, DeJongh M, Best AA, et al. (2010) High-throughput generation, optimization and analysis of genome-scale metabolic models. Nat Biotechnol 28: 977-982. doi: 10.1038/nbt.1672
    [32] Devoid S, Overbeek R, DeJongh M, et al. (2013) Automated genome annotation and metabolic model reconstruction in the SEED and Model SEED. Methods Mol Biol 985: 17-45. doi: 10.1007/978-1-62703-299-5_2
    [33] DeJongh M, Formsma K, Boillot P, et al. (2007) Toward the automated generation of genome-scale metabolic networks in the SEED. BMC Bioinformatics 8: 139. doi: 10.1186/1471-2105-8-139
    [34] Zakrzewski P, Medema MH, Gevorgyan A, et al. (2012) MultiMetEval: comparative and multi-objective analysis of genome-scale metabolic models. PloS One 7: e51511. doi: 10.1371/journal.pone.0051511
    [35] Hucka M, Finney A, Sauro HM, et al. (2003) The systems biology markup language (SBML): a medium for representation and exchange of biochemical network models. Bioinformatics 19: 524-531. doi: 10.1093/bioinformatics/btg015
    [36] Schellenberger J, Que R, Fleming RM, et al. (2011) Quantitative prediction of cellular metabolism with constraint-based models: the COBRA Toolbox v2.0. Nat Protoc 6: 1290-1307. doi: 10.1038/nprot.2011.308
    [37] McAnulty MJ, Yen JY, Freedman BG, et al. (2012) Genome-scale modeling using flux ratio constraints to enable metabolic engineering of clostridial metabolism in silico. BMC Syst Biol 6: 42. doi: 10.1186/1752-0509-6-42
    [38] Yen JY, Nazem-Bokaee H, Freedman BG, et al. (2013) Deriving metabolic engineering strategies from genome-scale modeling with flux ratio constraints. Biotechnol J 8: 581-594.
    [39] de Oliveira Dal'Molin CG, Quek LE, Palfreyman RW, et al. (2010) AraGEM, a genome-scale reconstruction of the primary metabolic network in Arabidopsis. Plant Physiol 152: 579-589. doi: 10.1104/pp.109.148817
    [40] Duarte NC, Herrgard MJ, Palsson BO (2004) Reconstruction and validation of Saccharomyces cerevisiae iND750, a fully compartmentalized genome-scale metabolic model. Genome Res 14: 1298-1309. doi: 10.1101/gr.2250904
    [41] Nogales J, Gudmundsson S, Knight EM, et al. (2012) Detailing the optimality of photosynthesis in cyanobacteria through systems biology analysis. Proc Natl Acad Sci U S A 109: 2678-2683. doi: 10.1073/pnas.1117907109
    [42] Feist AM, Henry CS, Reed JL, et al. (2007) A genome-scale metabolic reconstruction for Escherichia coli K-12 MG1655 that accounts for 1260 ORFs and thermodynamic information. Mol Syst Biol 3: 121.
    [43] Jones SW, Paredes CJ, Tracy B, et al. (2008) The transcriptional program underlying the physiology of clostridial sporulation. Genome Biol 9: R114. doi: 10.1186/gb-2008-9-7-r114
    [44] Lewis NE, Hixson KK, Conrad TM, et al. (2010) Omic data from evolved E. coli are consistent with computed optimal growth from genome-scale models. Mol Syst Biol 6: 390.
    [45] Keating SM, Bornstein BJ, Finney A, et al. (2006) SBMLToolbox: an SBML toolbox for MATLAB users. Bioinformatics 22: 1275-1277. doi: 10.1093/bioinformatics/btl111
    [46] Henry CS, Zinner JF, Cohoon MP, et al. (2009) iBsu1103: a new genome-scale metabolic model of Bacillus subtilis based on SEED annotations. Genome Biol 10: R69. doi: 10.1186/gb-2009-10-6-r69
    [47] Ro DK, Paradise EM, Ouellet M, et al. (2006) Production of the antimalarial drug precursor artemisinic acid in engineered yeast. Nature 440: 940-943. doi: 10.1038/nature04640
    [48] Senger RS, Papoutsakis ET (2008) Genome-scale model for Clostridium acetobutylicum: Part II. Development of specific proton flux states and numerically determined sub-systems. Biotechnol Bioeng 101: 1053-1071.
    [49] Baez A, Cho KM, Liao JC (2011) High-flux isobutanol production using engineered Escherichia coli: a bioreactor study with in situ product removal. Appl Microbiol Biotechnol 90: 1681-1690. doi: 10.1007/s00253-011-3173-y
    [50] Blombach B, Riester T, Wieschalka S, et al. (2011) Corynebacterium glutamicum tailored for efficient isobutanol production. Appl Environ Microbiol 77: 3300-3310. doi: 10.1128/AEM.02972-10
    [51] Li S, Wen J, Jia X (2011) Engineering Bacillus subtilis for isobutanol production by heterologous Ehrlich pathway construction and the biosynthetic 2-ketoisovalerate precursor pathway overexpression. Appl Microbiol Biotechnol 91: 577-589. doi: 10.1007/s00253-011-3280-9
    [52] Follonier S, Panke S, Zinn M (2011) A reduction in growth rate of Pseudomonas putida KT2442 counteracts productivity advances in medium-chain-length polyhydroxyalkanoate production from gluconate. Microb Cell Fact 10: 25.
    [53] Tracy BP, Jones SW, Fast AG, et al. (2012) Clostridia: the importance of their exceptional substrate and metabolite diversity for biofuel and biorefinery applications. Curr Opin Biotechnol 23: 364-381. doi: 10.1016/j.copbio.2011.10.008
    [54] Senger RS, Papoutsakis ET (2008) Genome-scale model for Clostridium acetobutylicum: Part I. Metabolic network resolution and analysis. Biotechnol Bioeng 101: 1036-1052.
    [55] Martin VJ, Pitera DJ, Withers ST, et al. (2003) Engineering a mevalonate pathway in Escherichia coli for production of terpenoids. Nat Biotechnol 21: 796-802. doi: 10.1038/nbt833
  • This article has been cited by:

    1. Shaowei Ding, Suprem R. Das, Benjamin J. Brownlee, Kshama Parate, Taylor M. Davis, Loreen R. Stromberg, Edward K.L. Chan, Joseph Katz, Brian D. Iverson, Jonathan C. Claussen, CIP2A immunosensor comprised of vertically-aligned carbon nanotube interdigitated electrodes towards point-of-care oral cancer screening, 2018, 117, 09565663, 68, 10.1016/j.bios.2018.04.016
    2. J. Alvarez-Ramirez, J. Vazquez-Arenas, M. González, Y. Carrera-Tarela, Non-Linear First-Harmonic Balance to Compute the Electrochemical Impedance of Butler-Volmer Equation, 2018, 165, 0013-4651, H517, 10.1149/2.1091809jes
    3. Cibely S. Martin, Priscila Alessio, 2020, 9780128166796, 359, 10.1016/B978-0-12-816679-6.00011-5
    4. Cristina Stefanov, Jacobus (Koos) Frederick van Staden, Raluca-Ioana Stefan-van Staden, Review—Enzymatic and Non-Enzymatic (bio)sensors Based on Phthalocyanines. A Minireview, 2020, 9, 2162-8777, 051012, 10.1149/2162-8777/ab9dc8
    5. Georgi B. Hadjichristov, Yordan G. Marinov, Todor E. Vlakhov, Nicola Scaramuzza, 2021, 34, 9780323914994, 129, 10.1016/bs.abl.2021.11.005
  • Reader Comments
  • © 2015 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(6507) PDF downloads(1166) Cited by(5)

Article outline

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog