Skip to main content
  • Research article
  • Open access
  • Published:

Kinetic modeling of nicotine in mainstream cigarette smoking



The attempt to understand the kinetic behavior of nicotine in tobacco will provide a basis for unraveling its energetics in tobacco burning and the formation of free radicals considered harmful to the cigarette smoking community. To the best of our knowledge, the high temperature destruction kinetic characteristics of nicotine have not been investigated before; hence this study is necessary especially at a time addiction science and tobacco research in general is gaining intense attention.


The pyrolysis of tobacco under conditions simulating cigarette smoking in the temperature region 200–700 °C has been investigated for the evolution of nicotine and pyridine from two commercial cigarettes coded ES1 and SM1 using gas chromatography hyphenated to a mass selective detector (MSD). Moreover, a kinetic model on the thermal destruction of nicotine within a temperature window of 673 and 973 K is proposed using pseudo-first order reaction kinetics. A reaction time of 2.0 s was employed in line with the average puff time in cigarette smoking. Nonetheless, various reaction times were considered for the formation kinetics of nicotine.


GC–MS results showed the amount of nicotine evolved decreased with increase in the puff time. This observation was remarkably consistent with UV–Vis data reported in this investigation. Generally, the temperature dependent rate constants for the destruction of nicotine were found to be \(k = 2.1\; \times \;10^{6} T^{n} \; \times \;e^{{ - \,\frac{108.85}{RT}}}\) s−1 and \(k = 3.0\; \times \;10^{7} T^{n} \; \times \;e^{{ - \,\frac{136.52}{RT}}}\) s−1 for ES1 and SM1 cigarettes respectively. In addition, the amount of nicotine evolved by ES1 cigarette was ~10 times more than the amount of nicotine released by SM1 cigarette.


The suggested mechanistic model for the formation of pyridine from the thermal degradation of nicotine in tobacco has been found to be agreement with the kinetic model proposed in this investigation. Consequently, the concentration of radical intermediates of tobacco smoke such as pyridinyl radical can be determined indirectly from a set of integrated rate laws. This study has also shown that different cigarettes can yield varying amounts of nicotine and pyridine depending on the type of cigarette primarily because of potential different growing conditions and additives introduced during tobacco processing. The activation energy of nicotine articulated in this work is consistent with that reported in literature.

The anatomy of tobacco cigarette and the major chemistry involved during combustion (pyrolysis, GC–MS analysis, and kinetic modeling)


Tobacco smoke is a highly dynamic and very complex matrix consisting of over 6000 compounds which makes a cigarette behave like a chemical reactor where several complex chemical processes take place during pyrolysis [16]. Pyrolysis can be described as the direct decomposition of an organic matrix to obtain a range of reaction products in limited oxygen [710]. Accordingly, the thermal degradation reaction mechanisms are complex and therefore it is necessary to simplify input parameters and physical properties in order to simulate the largest possible influence on the overall kinetic characteristics of biomass pyrolysis including tobacco [8, 9]. A kinetic scheme of biomass pyrolysis must therefore involve the solution of a high-dimensional system of differential equations [1113].

The thermal destruction of nicotine in this investigation was conducted within a temperature window of 673 and 973 K at an average reaction time of 2.0 s as reported in literature [1416]. For simplicity, a consecutive first order reaction with rate constants k 1 and k 2 has been considered in which a global kinetic model [1720] was employed to obtain the kinetic parameters for the thermal destruction of nicotine in mainstream cigarette smoking. Accordingly, pseudo-unimolecular reactions were applied in which the empirical rate of decomposition of the initial product is first order and expressed by Eq. 1.

$$C = C_{o} e^{ - kt}$$

where Co and C are respective concentrations of the reactant at time, t = 0, and time, t = 2.0 s, while k is the pseudo-unimolecular rate constant in the Arrhenius expression (cf. Eq. 2).

$$k = Ae^{{ - \;\frac{Ea}{RT}}}$$

A is the pre-exponential factor (s−1), Ea is the activation energy (kJmol−1), R is the universal gas constant (8.314 JK−1mol−1), and T is the temperature in K. Despite all the criticisms against the Arrhenius rate law, it remains the only kinetic expression that can satisfactorily account for the temperature-dependent behavior of even the most unconventional reactions including biomass pyrolysis [9]. The integrated form of the first order rate law (cf. Eq. 3) was used to calculate the rate constant for the pyrolysis behavior of tobacco at a reaction time of 2.0 s.

$$k = \ln \left( {\frac{Co}{C}} \right)\frac{1}{t}$$

The activation energy was determined from the Arrhenius plots (ln k vs. \(1/T\)) which establishes a linear relationship between the pre-exponential factor A and the rate constant k as given by Eq. 4, where ln A is the y-intercept and \(- \frac{Ea}{RT}\) is the slope.

$$\ln k = \ln A - \frac{Ea}{RT}$$

To the best of our knowledge, there is no known destruction kinetic modeling of nicotine reported in literature. Consequently, this is perhaps the first such study on the destruction kinetics of nicotine. Although, the results obtained in this study are estimated from experimental data and may require further tests, we believe this an important step in the study of kinetics of reaction products in complex biomass materials such as plant matter. In this work, we have used GC-Area counts to determine the destruction rate constants because according to the first order reaction kinetics (Eq. 3, vide infra) the ratio of concentrations at various temperatures is a constant. Therefore, calibration of nicotine will still achieve similar results.

The primary focus of this study is to give a general kinetic account of the destruction kinetics of nicotine and demonstrate how the concentration of intermediates, in this case, pyridinyl radical can be determined indirectly and estimate the kinetic parameters of nicotine in ES1 and SM1 cigarette. The kinetics of nicotine destruction is based on high temperature regimes characteristic of cigarette burning [16, 21]. The results reported in this investigation are no doubt different from the kinetics of nicotine inhaled into the blood system which is beyond the scope of this study. Therefore, this work considers only the gas-phase kinetics of nicotine deemed fundamental towards understanding the inhalation kinetics of mainstream cigarette smoke. Furthermore, attempts have been made to identify and describe kinetically the intermediate radicals produced by the thermal degradation of nicotine from two different commercial cigarette samples (ES1 and SM1). Radicals such as pyridinyl radical which is the focus of this work have been known to cause serious health impacts because they are highly reactive towards biological tissues such as DNA, lipids, and microphages [2225]. Free radicals such as pyridinyl radical has the ability to generate reactive oxygen species when it reacts with biological tissues and thus accelerating the growth of tumours, cancer cells, cell injury and oxidative stress [2527].

From a quantum chemical perspective, the scission of the phenyl C–C linkage in nicotine has been explored using the density functional theory (DFT) in order to determine the energetics for the formation of pyridinyl radical from pure nicotine (in absence of other tobacco components). Although this is critical in understanding the mechanistic formation of pyridine from nicotine, it will only be discussed briefly.

Experimental protocol


The heater (muffle furnace) was purchased from Thermo Scientific Inc., USA while the quartz reactor was locally fabricated in our laboratory by a glass-blower. Commercial cigarettes coded SM1 and ES1 (for confidential reasons, cannot be revealed) were purchased from retail outlets and used without further treatment. Methanol (purity >>99 %) used to dissolve cigarette pyrolysate was purchased from Sigma Aldrich Inc. (USA). All experiments in this work were conducted under ISO conditions reported in Reference [16].

Sample preparation

Processed tobacco (from ES1 and SM1) of 50 ± 0.2 mg was weight and packed in a quartz reactor of dimensions: i.e. 1 cm × 2 cm (volume ≈ 1.6 cm3). The tobacco sample in the quartz reactor was placed in an electrical heater furnace whose maximum heating temperature is 1000 °C. The tobacco sample was heated in flowing nitrogen (pyrolysis gas) and the smoke effluent was allowed to pass through a transfer column and collected in 10 mL methanol in a conical flask for a total pyrolysis time of 2 min and sampled into a 2 mL crimp top amber vials for GC–MS analysis. The pyrolysis gas flow rate was designed to maintain a constant residence time of 2.0 s representative of cigarette smoking [1416, 28]. The goal of many studies, however; is to establish the relationship between tobacco constituents and smoke products under conditions that simulate actual human smoking, but this desire remains a challenge because of the large number of processes occurring inside a burning cigarette involving varying temperatures and changes in oxygen concentration [3, 4]. It turns out that the burning conditions in a cigarette change significantly from the way the cigarette burns from the oxygen rich peripheral surface towards the interior of the cigarette where oxygen is either low or generally absent [28]. This combustion experiment was conducted under conventional pyrolysis described in literature [29] and the evolution of nicotine and pyridine were monitored between 200 and 700 °C as shown in Fig. 5.

GC–MS determination of nicotine and pyridine from ES1 and SM1 tobacco

Analysis of nicotine and pyridine was carried out using an Agilent Technologies 7890A GC system connected to an Agilent Technologies 5975C inert XL Electron Ionization/Chemical Ionization (EI/CI) with a triple axis mass selective detector, using HP-5MS 5 % phenyl methyl siloxane column (30 m × 250 µm × 0.25 µm). The temperature of the injector port was set at 200 °C to vaporize the organic components for GC–MS analysis. The carrier gas was ultra-high pure (UHP) helium (99.999 %) and the flow rate was 3.3 mL min−1. Temperature programming was applied at a heating rate of 15 °C for 10 min, holding for 1 min at 200 °C, followed by a heating rate of 25 °C for 4 min, and holding for 10 min at 300 °C. Electron Impact ionization energy of 70 eV was used. To ensure that the right compounds were detected, standards were run through the GC–MS system and the peak shapes as well as retention times were compared with those of nicotine and pyridine. The data was run through the NIST and the Agilent Chemstation library databases—MS-fragmentation patterns, as additional tools to confirm the identity of the compounds (nicotine and pyridine) [29]. The MS-fragmentation patterns for these compounds are presented in the support information (Additional file 1): S1(MS-Fragmentation pattern of nicotine) and S2(MS-Fragmentation pattern of pyridine). Experimental results were averaged replicates of two or more data points.

GC–MS and UV–Vis analysis of nicotine in ES1 cigarette

The rate of formation of nicotine from ES1 cigarette was determined experimentally at modest puff times (2, 5, and 10 s) using laboratory designed apparatus (Fig. 1). For every puff time, the concentration of nicotine was determined using a GC–MS hyphenated to a mass selective detector as discussed in the section above. To qualify the characteristic kinetics for the formation of nicotine at various puff times, the absorbance measurements of nicotine were taken and absorbance curves plotted. The results remarkably were similar to the GC–MS data. Maximum absorbance of nicotine in UV–Vis occurred at 220 nm. The absorbance was confirmed by running nicotine standard through the UV–Vis instrument. Methanol was used as a blank in UV–Vis analysis. The model of the instrument used for UV–Vis analysis was SHIMADZU, UV 1800.

Fig. 1
figure 1

Apparatus set up for trapping cigarette smoke from cigarette burning

The kinetic model

During the kinetic modeling of nicotine from the thermal degradation of tobacco biomass, decent assumptions were considered (Fig. 2): (1) the rate of formation of nicotine prevails the rate of destruction, (2) at the peak of the curve, the rates of formation and destruction are approximately the same, and (3) as the temperature is increased, the rate of destruction overwhelms the rate of formation. These assumptions are made based on the fact that pyrolysis of tobacco leads to the formation of nicotine, one of the major tobacco alkaloids as articulated in literature [6, 24, 30, 31]. This is consistent with our experiments which show that the pyrolysis of tobacco yields significant amounts of nicotine (Fig. 4). Therefore, from these assumptions, it is possible to determine the apparent kinetic parameters for the destruction of nicotine from the temperature dependence of its yields. A simple single step reaction mechanism during the thermal degradation of nicotine as presented in Eq. (5) is considered. Although tobacco pyrolysis is very complex, we believe some understanding on the kinetic behavior of certain reaction products from basic kinetic equations can be deduced. Therefore, modeling does not necessarily need to be complex to describe complex reactions systems. In essence, even simple models based on relevant assumptions may yield reasonable results as presented in this work, and compared with literature data.

Fig. 2
figure 2

The relationship between the rates of formation of the intermediate product (Rf) vs. the rate of destruction (Rd). Co is taken as the maximum concentration of the reaction product

$$Nic \mathop{\longrightarrow}\limits_{}^{k_{1}} I \mathop{\longrightarrow}\limits_{}^{k_{2}} Product$$

Conventionally, the differential rate laws for each species Nic (nicotine), I (intermediate), and the final product are given by Eqs. 6, 7, and 8 respectively.

$$\frac{d\,[Nic]}{dt} = - k_{1} [Nic]$$
$$\frac{d\,[I]}{dt} = k_{1} [Nic] - k_{2} [I]$$
$$\frac{d\,[Product]}{dt} = k_{2} [I]$$

If these equations are solved analytically, then the integrated rate laws are as given by Eqs. 9 and 10.

$$[Nic] = [Nic]_{0} e^{{ - k_{1} t}}$$

Equations 10 and 11 give the respective concentrations of the intermediate I and the product at any time t.

$$[I] = \frac{{k_{1} [Nic]_{0} }}{{k_{2} - k_{1} }}\left( {e^{{ - k_{1} t}} - e^{{ - k_{2} t}} } \right)$$
$$[Product] = [Nic]_{0} \left[ {1 + \frac{{k_{1} }}{{k_{1} - k_{2} }}\left( {k_{2} e^{{ - k_{1} t}} - k_{1} e^{{ - k_{2} t}} } \right)} \right]$$

In order to simplify Eq. 11 further, we will assume that step two (Eq. 5) is the rate determining step so that k 2 << k 1 and thus the term \(e^{{ - k_{1} t}}\) decays more rapidly than the term \(e^{{ - k_{2} t}}\) [32]. Therefore Eq. 11 reduces to Eq. 12. This assumption is valid based on previous studies documented in literature [11, 32, 33].

$$[Product] = [Nic]_{0} \left( {1 - e^{{ - k_{2} t}} } \right)$$

Results and discussion

To mimic actual cigarette smoking conditions, smoking apparatus were designed according to ISO 3402:1999 standards [16]. Whereas the destruction kinetics of nicotine was explored for both ES1 and SM1 cigarettes, only ES1 cigarette was investigated for nicotine formation. For formation kinetics, smoking residence times usually representative of real world cigarette smoking conditions (2, 5, and 10 s) were explored. Consequently, a plot of ln k as a function of puff (smoking) time yielded a straight line with a slope of −0.1323 (Fig. 3) from which the formation rate constant of nicotine (0.13 s−1) was calculated. The plot, although an estimation from restricted smoking times is consistent with first order reaction kinetics. The original amount of nicotine in ES1 cigarette was estimated from the y-intercept and established to be 9.1 × 108 GC-Area counts. This value is remarkably close to that obtained from experimental modeling of tobacco burning from ES1, ~8.0 × 108 GC-Area counts.

Fig. 3
figure 3

Formation kinetics of nicotine (A) and absorbance of nicotine at various puff times (B) in ES1 cigarette

Interesting data have been reported in this work concerning the decrease of nicotine with smoking times (Fig. 3A, B). This suggests that longer residence times may lead to possible side reactions which result in the conversion of nicotine to other by-products. It is well known in literature that shorter residence times minimize secondary reactions but longer residence times may lead to radical formation, recombination, and pyrosynthesis of new by-products [29, 34]. Thus, these processes reduce the yield of the parent compound, in this case, nicotine. The UV–Vis data was basically qualitative but remarkably corroborates GC–MS data. Therefore, the longer the smoking times the lower the concentration of nicotine reaching the lungs of the cigarette smoker. Longer puff times may be beneficial to the smoking community based on the results obtained from this work.

Molecular distribution of nicotine and pyridine

The product distribution of nicotine in the temperature region 200–700 °C is presented in Fig. 4. Clearly, ES1 cigarette yielded high levels of nicotine and pyridine in comparison to SM1 cigarette. The nicotine levels from the two commercial cigarettes peaked at different pyrolysis temperatures. For instance, nicotine from ESI peaked at 400 °C while nicotine from SM1 peaked at about 500 °C. Interestingly, pyridine from the two cigarettes reached a maximum at about 500 °C. The two cigarettes, based on this data are significantly different. This result may be attributed to possible different growing conditions and additives during the processing of the two cigarettes. Interestingly, the total nicotine content in the entire pyrolysis range in ES1 tobacco was ~10 times the amount of nicotine released by SM1 tobacco in the same pyrolysis temperature region (200–700 °C). This may imply that SM1 cigarette is much safer than ES1 cigarette based on nicotine and pyridine data alone presented in this study. Accordingly, a close examination of the curves in Fig. 4 indicates that nicotine from the pyrolysis of tobacco is formed even at lower temperatures than the lowest temperature selected in this study (200 °C). This behaviour is explained in literature [6]. Accordingly, Forster et al. [6] proposes that the concentration of nicotine should increase with increase in the pyrolysis temperature hence the shift in nicotine yields at 200 °C as presented in Fig. 4.

Fig. 4
figure 4

Evolution of nicotine and pyridine from ES1 and SM1 cigarette tobacco

The overlay chromatograms showing the formation of nicotine and pyridine at two pyrolysis temperatures (300–400 °C) is presented in Fig. 5. Clearly, from Fig. 5, nicotine has a high intensity at 400 °C in agreement with predictions made by Forster et al. [6]. The intensity of pyridine also increases with increase in temperature. Nonetheless, like other reaction products of tobacco and other biomass pyrolysis, nicotine peaks between 300 and 500 °C before decreasing significantly with increase in temperature [29, 30, 35] (Fig. 4). The region where the concentration of nicotine begins to decrease with increase in temperature as illustrated in Fig. 2 forms the basis for modeling the destruction kinetics of nicotine which is the main subject of this investigation.

Fig. 5
figure 5

Overlay chromatograms showing the peaks for pyridine and nicotine for the pyrolysis of ES1 tobacco at 300 °C (red line) and 400 °C (blue line)

Destruction kinetics of nicotine

The destruction kinetics revealed that nicotine from ES1 has activation energy of 108.85 kJmol−1 while SM1 has activation energy of 136.52 kJmol−1 (Table 1). This implies that the two cigarettes may have different matrix composition. Thus the activation energies of nicotine in the two cigarettes may not necessarily be the same considering the fact that additives of varying composition introduced during cigarettes processing may act as catalysts and ultimately reduce the activation energy of a given compound in a complex biomass material such as tobacco. Remarkably, the activation energy determined from this study is comparable to that documented in literature in which the average activation energy of nicotine was found to be 120 kJmol−1 [6]. Moreover, the activation energies determined from this work are similar to the results from the kinetic modeling of the pyrolysis of other biomass materials such as cellulose [11]. Arrhenius plots for the destruction of nicotine from the cigarettes under study are presented in Fig. 6. Nonetheless, in modeling the destruction kinetics of nicotine, we are aware that the kinetic characteristics of a given heterogeneous system such as plant matter may change during the process of pyrolysis and so it is possible that the complete reaction mechanism cannot be represented adequately by a specific kinetic model [9, 36]. Although we have assumed a linear relationship between ln k and \(1/T\) we note that not all reactions will necessarily obey this relation. Therefore in order to estimate the Arrhenius dependent rate constants consistent with experimental rate constants, the modified Arrhenius rate expression is applied.

Table 1 The Arrhenius parameters for the destruction of nicotine from the pyrolysis of ES1 and SM1 cigarette tobacco
Fig. 6
figure 6

Arrhenius plots for the destruction kinetics of nicotine in ES1 and SM1 cigarette tobacco

$$k = AT^{n} e^{{ - \;\frac{Ea}{RT}}}$$

For a given temperature, since the rate constant has been determine experimentally and all the other parameters are known, the value of n can be determined from Eq. 13. For instance, the value of n at 673 K was determined and found to be 0.55 and 1.05 for the destruction of nicotine in ES1 and SM1 cigarettes respectively. Equation 13 can be used to calculate the value of n at any particular temperature, since the rate constants are temperature dependent.

The destruction rate constant k 1 at 673 K for ES1 was 0.31 s−1 while that of SM1 at the same temperature was estimated as 0.74 s−1. At the highest pyrolysis temperature (973 K), the respective rate constants were 2.12–1.0 s−1. Accordingly, the average destruction rate constant for ES1 was found to be 1.11 s−1. Table 1 presents the Arrhenius parameters from the destruction kinetics of nicotine (Activation energies and Arrhenius factors). Whereas the activation energies are comparably close, the pre-exponential factors for the two cigarettes under study differ by a whole magnitude.

If wish to calculate the rate constant k 2 for the formation of the product, for instance pyridine (a by-product of nicotine pyrolysis), then we will need to use the differential rate law provided in Eq. 12. To be able to do this, serious assumptions have to be taken into account. For instance, one of the major by-products from the destruction of nicotine pyrolysis must be pyridine [30, 37]. This assumption is valid if we take into consideration the reactive nature of the H radical relative to the methyl radical which may yield 3-methylpyridine (a minor product) [20, 29]. Furthermore it has been proven experimentally that one of the major by-products from the thermal destruction of nicotine is pyridine [4, 30]. These findings corroborate our kinetic model on the thermal destruction of nicotine at high temperature smoking regimes.

Therefore, by substituting the original concentration of nicotine for ES1 (8.0 × 108 GC-Area counts) and the maximum concentration of the product, in this case, pyridine (4.4 × 108 GC-Area counts) into Eq. 12, vide supra, the value of k 2 was computed and found to be 0.13 s−1. This shows that the value of k 2 is less than the value of k 1 by 1 magnitude. Secondly, since the rate constants k 1 and k 2 have been estimated, and the original value of nicotine is known, then the concentration of the intermediate, pyridinyl radical, can be calculated from Eq. 10. Accordingly, the concentration of pyridinyl radical was determined as 6.1 × 108 GC-Area counts. Similar calculations were conducted for the kinetics of nicotine in SM1 cigarette and the value of k 2 was estimated as 0.67 s−1 while its pyridinyl radical intermediate had a concentration of 3.31 × 108 GC-Area counts. From these data, the concentration of pyridinyl radical in ES1 is ~2 times the concentration of pyridinyl radical in SM1.

Evidently, the sum of the concentrations of the intermediate and the proposed final product (pyridine) for each cigarette was greater than the original concentration of nicotine evolved by each cigarette. This is expected because in the pyrolysis of a complex matrix such as plant matter, various heterogeneous reactions occur. Thus the thermal degradation of nicotine may not be the only route for pyridine formation. This argument is acceptable if we consider experimentally that both nicotine and pyridine are evolved simultaneously during pyrolysis (Fig. 5). Nevertheless, nicotine destruction is suggested as the major route for the formation of pyridine [30, 37]. The ratio of original nicotine to the sum of concentrations of the intermediate (pyridinyl radical) and pyridine for ES1 and SM1 cigarettes were respectively 0.76 and 0.60. On the other hand, the ratio of pyridine (presumed the major by-product of nicotine destruction) to the original nicotine was determined as 0.55 and 0.52 for ES1 and SM1 respectively. These findings indicate that it might be possible that ~45 % of nicotine in ES1 and ~48 % in SM1 may have been transferred intact into the smoker. Schmeltz et al. [30] puts this figure at <41 %. This discrepancy may be attributed to a number of factors; the type of tobacco and the pyrolysis conditions. In our study, we have used an inert atmosphere to simulate cigarette smoking which implies extensive fragmentation may occur during the thermal degradation of tobacco resulting in high yields of pyridine as reported in literature [4].

Mechanistic description for the formation of pyridine from nicotine

It is possible by inspection to envisage that the scission of the C–C phenyl bond in nicotine should result in the formation of pyridine despite the complex nature of pyrolytic processes taking place in plant matter such as tobacco. In order to appreciate this assumption, we have designed a mechanistic model for the formation of pyridine from nicotine as presented in Scheme 1 to support our kinetic model. Rearrangement and dehydrogenation reactions that may yield compounds such as β-nicotyrine from nicotine may not be thermodynamically feasible. This is in agreement with our experimental results in which insignificant yields of β-nicotyrine were detected. The other assumption is 1-methylpyrrolidine is a minor product. From an experimental perspective, this assumption is true because no 1-methylpyrrolidine was detected in the entire range of tobacco pyrolysis whereas significant amounts of pyridine was detected, Fig. 5, vide supra. Although, pyridine may not be the only by-product of nicotine decomposition owing to the complex processes occurring during tobacco pyrolysis, it is definitely one of the major products [6, 30]. Nonetheless, its yields depends entirely on the growing conditions of tobacco, additives introduced during tobacco processing, and the pyrolysis atmosphere in tobacco burning. This observation is clear based on the results of the two cigarettes reported in this study.

Scheme 1
scheme 1

Mechanistic destruction of nicotine to radical intermediates and possible by-products

The bond dissociation energy via the rate constant k 1 and the bond formation energy via rate constant k 2 (scheme 1) were estimated using the density functional theory framework at the B3LYP energy functional in conjunction with 6-31G basis set. Nonetheless, the bond energies will not be discussed further because they are the subject of critical discussions in our next article. The scheme, however; proposes a plausible mechanistic pathway for the thermal degradation of nicotine to the intermediate (pyridinyl radical) and ultimately to pyridine.

Toxicological impacts of nicotine, pyridine, and pyridinyl radical

Animal studies support biological evidence for accelerated motor activity, neurobehavioral, learning and memory deficits, and alteration of neurotransmitter function due to exposure to nicotine [38, 39]. Nicotine also affects the cardiovascular system in many ways that is by activating the sympathetic nervous system; nicotine induces increased heart rate and myocardial contraction, vasoconstriction in the skin and adrenal, reproductive problems and neural release of catecholamine [4042]. Nicotine can also affect lipid metabolism [43], accelerate the development of atherosclerosis [44], induce endothelial dysfunction [45], and has been suspected as a carcinogen [42]. After a puff, high levels of nicotine reach the brain in 10–20 s, faster than with intravenous administration, producing rapid behavioural reinforcement [46]. On the other hand, pyridine has been implicated in the inhibition of the growth of chick chorioallantoic membrane and reproductive health issues [37, 47, 48]. In this study, the radicals including pyridinyl and 1-methylpyrrolidinyl radicals are good candidates for cell injury and oxidative stress during cigarette smoking. The molecular structure of nicotine and other alkaloid related compounds investigated in this work may covalently bond to the DNA, lipids, nuclei acids, and body cells before metabolizing into harmful by-products that are potential risks to the human health [23, 26, 27, 42]. In addition, pyridinyl radical can react with biological molecules to enhance the production of reactive oxygen species which can cause oxidative stress, tumourogens, and cancer [23, 4952].


The temperature dependent destruction kinetics of nicotine has been presented for the first time in this investigation. A mechanistic model showing the formation of pyridine from the thermal destruction of nicotine has been proposed and found to be in agreement with the kinetic model reported in this study. We therefore believe the results presented in this investigation will form the basis of further research towards understanding the fate of nicotine during cigarette smoking. The two cigarettes investigated in this work coded ES1 and SM1 have exhibited various kinetic characteristics possibly because of their different biomass composition attributed mainly to their growing conditions and additives during tobacco processing. Moreover, this study has established that the activation energy of nicotine is remarkably consistent with that reported in literature. The concentration of the intermediate (pyridinyl radical) has been estimated from kinetic modeling of nicotine. This is remarkable since the concentrations of intermediates in complex reaction systems such as biomass are usually tedious to determine experimentally.


  1. Adam T et al (2010) Influence of filter ventilation on the chemical composition of cigarette mainstream smoke. Anal Chim Acta 657(1):36–44

    Article  CAS  Google Scholar 

  2. Adam T et al (2009) Simultaneous on-line size and chemical analysis of gas phase and particulate phase of cigarette mainstream smoke. Anal Bioanal Chem 394(4):1193–1203

    Article  CAS  Google Scholar 

  3. Czegeny Z et al (2009) Formation of selected toxicants from tobacco under different pyrolysis conditions. J Anal Appl Pyrol 85(1–2):47–53

    Article  CAS  Google Scholar 

  4. Baker RR, Bishop LJ (2004) The pyrolysis of tobacco ingredients. J Anal Appl Pyrol 71(1):223–311

    Article  CAS  Google Scholar 

  5. Baker RR (2006) Smoke generation inside a burning cigarette: modifying combustion to develop cigarettes that may be less hazardous to health. Prog Energy Combust Sci 32(4):373–385

    Article  CAS  Google Scholar 

  6. Forster M et al (2015) An experimental method to study emissions from heated tobacco between 100–200 °C. Chem Cent J 9(1):1–10

    Article  CAS  Google Scholar 

  7. Babu BV (2008) Biomass pyrolysis: a state-of-the-art review. Biofuels Bioprod Biorefin 2(5):393–414

    Article  CAS  Google Scholar 

  8. Huang L et al (2004) Rapid nitric acid digestion of plant material with an open-vessel microwave system. Commun Soil Sci Plant Anal 35(3–4):427–440

    Article  CAS  Google Scholar 

  9. White JE, Catallo WJ, Legendre BL (2011) Biomass pyrolysis kinetics: a comparative critical review with relevant agricultural residue case studies. J Anal Appl Pyrol 91(1):1–33

    Article  CAS  Google Scholar 

  10. Demirbas A (2009) Biorefineries: current activities and future developments. Energy Convers Manag 50(11):2782–2801

    Article  CAS  Google Scholar 

  11. Zhang XL, Yang WH, Blasiak W (2012) Kinetics of levoglucosan and formaldehyde formation during cellulose pyrolysis process. Fuel 96(1):383–391

    Article  CAS  Google Scholar 

  12. Zhang XL, Yang WH, Blasiak W (2011) Modeling study of woody biomass: interactions of cellulose, hemicellulose, and lignin. Energy Fuels 25(10):4786–4795

    Article  CAS  Google Scholar 

  13. Chen HX, Liu NA, Fan WC (2006) Two-step consecutive reaction model and kinetic parameters relevant to the decomposition of Chinese forest fuels. J Appl Polym Sci 102(1):571–576

    Article  CAS  Google Scholar 

  14. Cardoso CR et al (2011) Determination of kinetic parameters and analytical pyrolysis of tobacco waste and sorghum bagasse. J Anal Appl Pyrol 92(2):392–400

    Article  CAS  Google Scholar 

  15. Seeman JI et al (1999) The form of nicotine in tobacco: thermal transfer of nicotine and nicotinic acid salts to nicotine in the gas phase. J Agric Food Chem 47(12):5133–5145

    Article  CAS  Google Scholar 

  16. Busch C et al (2012) Pyrolysis and combustion of tobacco in a cigarette smoking simulator under air and nitrogen atmosphere. Anal Bioanal Chem 403(2):419–430

    Article  CAS  Google Scholar 

  17. Dellinger B et al (2010) Theoretical estimation of incinerability of halons and hydrochlorofluorocarbons. Environ Eng Sci 27(7):587–591

    Article  CAS  Google Scholar 

  18. Taylor PH, Tirey DA, Dellinger B (1996) A detailed kinetic model of high-temperature pyrolysis of tetrachloroethene (vol 104, pg 260, 1996). Combust Flame 107(1–2):193–195

    Google Scholar 

  19. Taylor PH, Dellinger B, Lee CC (1990) Development of a thermal-stability based ranking of hazardous organic-compound incinerability. Environ Sci Technol 24(3):316–328

    Article  CAS  Google Scholar 

  20. Khachatryan L, Asatryan R, Dellinger B (2003) Development of expanded and core kinetic models for the gas phase formation of dioxins from chlorinated phenols. Chemosphere 52(4):695–708

    Article  CAS  Google Scholar 

  21. Zhou S et al (2011) The pyrolysis of cigarette paper under the conditions that simulate cigarette smouldering and puffing. J Therm Anal Calorim 104(3):1097–1106

    Article  CAS  Google Scholar 

  22. Dellinger B et al (2001) Role of free radicals in the toxicity of airborne fine particulate matter. Chem Res Toxicol 14(10):1371–1377

    Article  CAS  Google Scholar 

  23. Dellinger B et al (2000) The role of combustion-generated radicals in the toxicity of PM2.5. Proc Combust Inst 28:2675–2681

    Article  CAS  Google Scholar 

  24. Borgerdinga M, Klus K (2005) Analysis of complex mixtures—cigarette smoke. Exp Toxicol Pathol 57:43–73

    Article  Google Scholar 

  25. Kehrer JP et al (1988) Contemporary issues in toxicology. Toxicol Appl Pharmacol 95:349–362

    Article  CAS  Google Scholar 

  26. Bolton JL et al (2000) Role of quinones in toxicology. Chem Res Toxicol 13(3):135–160

    Article  CAS  Google Scholar 

  27. Sidney DN, Pearson PG (1990) Covalent and non covalent interactions in acute lethal cell injury caused by chemicals. Annu Rev Pharmacol Toxicol 30:169–195

    Article  Google Scholar 

  28. Thielen A, Klus H, Muller L (2008) Tobacco smoke: unraveling a controversial subject. Exp Toxicol Pathol 60(2–3):141–156

    Article  CAS  Google Scholar 

  29. Kibet J, Khachatryan L, Dellinger B (2012) Molecular products and radicals from pyrolysis of lignin. Environ Sci Technol 46(23):12994–13001

    Article  CAS  Google Scholar 

  30. Schmeltz I et al (1979) Chemical studies in tobacco smoke. 63. The fate of nicotine during pyrolysis and in burning a cigarette. J Agric Food Chem 27(3):602–608

    Article  CAS  Google Scholar 

  31. Armitage AK et al (2004) The effect of inhalation volume and breath-hold duration on the retention of nicotine and solanesol in the human respiratory tract and on subsequent plasma nicotine concentrations during cigarette smoking. Beitr Tabakforsch Int 21:240–249

    CAS  Google Scholar 

  32. McQuarrie DA, Simon DJ (1997) Physical chemistry: a molecular approach. In: McGuire A (ed) University Science Books, vol 1. Edwards Brothers Inc, Sausalito, pp 1–1360

    Google Scholar 

  33. Zhang XL et al (2011) formation mechanism of levoglucosan and formaldehyde during cellulose pyrolysis. Energy Fuels 25(8):3739–3746

    Article  CAS  Google Scholar 

  34. Kibet J, Khachatryan L, Dellinger B (2013) Molecular products from the pyrolysis and oxidative pyrolysis of tyrosine. Chemosphere 91(7):1026–1034

    Article  CAS  Google Scholar 

  35. Sharma RK et al (2004) Characterization of chars from pyrolysis of lignin. Fuel 83(11–12):1469–1482

    Article  CAS  Google Scholar 

  36. Wang SR et al (2012) Mechanism research on cellulose pyrolysis by Py-GC/MS and subsequent density functional theory studies. Bioresour Technol 104:722–728

    Article  CAS  Google Scholar 

  37. Kulshreshtha NP, Moldoveanu SC (1981) Analysis of pyridines in mainstream cigarette smoke. J Chromatogr A 985:303–312

    Article  Google Scholar 

  38. Zhao G et al (2012) Trends in health-related behavioral risk factors among pregnant women in the United States: 2001–2009. J Women’s Health 21(3):255–263

    Article  Google Scholar 

  39. Cho K et al (2013) Prenatal exposure to nicotine and impaired reading performance. J Pediatr 162(4):713–718

    Article  CAS  Google Scholar 

  40. Benowitz N (1996) Pharmacology of nicotine: addiction and therapeutics. Annu Rev Pharmacol Toxicol 36:597–613

    Article  CAS  Google Scholar 

  41. Cho K, Frijters JC, Zhang H, Miller LL, Gruen JR (2013) Prenatal exposure to nicotine and impaired reading performance. J Pediatr 162:713–718

    Article  CAS  Google Scholar 

  42. Elshal MF et al (2009) Sperm head defects and disturbances in spermatozoal chromatin and DNA integrities in idiopathic infertile subjects: association with cigarette smoking. Clin Biochem 42:589–594

    Article  CAS  Google Scholar 

  43. Cluette-Brown J et al (1986) Oral nicotine induces an atherogenic lipoprotein profile. Proc Soc Exp Biol Med 182:409–413

    Article  CAS  Google Scholar 

  44. Strohschneider T et al (1994) Effect of chronic nicotine delivery on the proliferation rate of endothelial and smooth muscle cells in experimentally induced vascular wall plaques. J Clin Investig 72:908–912

    CAS  Google Scholar 

  45. Chalon S et al (2000) Nicotine impairs endothelium-dependent dilatation in human veins in vivo. Clin Pharmacol Ther 67:391–397

    Article  CAS  Google Scholar 

  46. Benowitz N (1990) Clinical pharmacology of inhaled drugs of abuse: implications in understanding nicotine dependence. NIDA Res Monogr J 99:12–29

    CAS  Google Scholar 

  47. Ji L et al (2002) Identification of pyridine compounds in cigarette smoke solution that inhibit growth of chick chorioallantoic membrane. Toxicol Sci 69:217–225

    Article  CAS  Google Scholar 

  48. Riveles K et al (2003) Pyridines in cigarette smoke inhibit hamster oviductal functioning in picomolar doses. Reprod Toxicol 17:191–202

    Article  CAS  Google Scholar 

  49. Cormier SA et al (2006) Origin and health impacts of emissions of toxic by-products and fine particles from combustion and thermal treatment of hazardous wastes and materials. Environ Health Perspect 114(6):810–817

    Article  CAS  Google Scholar 

  50. Pryor WA, Prier DG, Church DF (1983) ESR study of mainstream and sidestream cigarette smoke: nature of free radicals in gas-phase smoke and in cigarette tar. Environ Health Perspect 47:345–355

    Article  CAS  Google Scholar 

  51. Kapur S et al (2000) Toxicology of benzyl alcohols: a QSAR analysis. Chemosphere 41(10):1643–1649

    Article  CAS  Google Scholar 

  52. Selassie D et al (1999) On the toxicity of phenols to fast growing cells. A QSAR model for a radical-based toxicity. J Chem Soc Perkin Trans 2(12):2729–2733

    Article  Google Scholar 

Download references

Authors’ contributions

CK prepared tobacco and cigarette samples, and conducted experimental analysis of nicotine and pyridine using GC–MS. CK also conducted UV–Vis analysis of nicotine under the supervision of JK and SM, and wrote the first draft of the manuscript. JK offered technical support during data interpretation, calculations, and compiling of the manuscript. SM and NK with technical advice from JK designed the mechanistic pathway for the conversion of nicotine to pyridinyl radical intermediate, and ultimately to pyridine. NK conducted computational calculations reported in this investigation. JB assisted CK during sample preparation and helped proof read the manuscript before it was submitted to JK for critical review and final editing. All authors read and approved the final manuscript.


The authors appreciate partial funding from the Directorate of Research & Extension (R&E) at Egerton University (Njoro).

Competing interests

The authors declare that there are no competing interests regarding the publication of this article.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Joshua Kibet.

Additional file


Additional file 1. This has beeb corrected accordingly under the section GC-MS determination of nicotine and pyridine in ES1 and SM1 tobacco.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kibet, J., Kurgat, C., Limo, S. et al. Kinetic modeling of nicotine in mainstream cigarette smoking. Chemistry Central Journal 10, 60 (2016).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: