Open Access

A simple method for HPLC retention time prediction: linear calibration using two reference substances

  • Lei Sun1, 2,
  • Hong-yu Jin1,
  • Run-tao Tian3,
  • Ming-juan Wang1,
  • Li-na Liu1,
  • Liu-ping Ye4,
  • Tian-tian Zuo1 and
  • Shuang-cheng Ma1Email author
Chinese Medicine201712:16

DOI: 10.1186/s13020-017-0137-x

Received: 20 February 2017

Accepted: 30 May 2017

Published: 19 June 2017



Analysis of related substances in pharmaceutical chemicals and multi-components in traditional Chinese medicines needs bulk of reference substances to identify the chromatographic peaks accurately. But the reference substances are costly. Thus, the relative retention (RR) method has been widely adopted in pharmacopoeias and literatures for characterizing HPLC behaviors of those reference substances unavailable. The problem is it is difficult to reproduce the RR on different columns due to the error between measured retention time (tR) and predicted tR in some cases. Therefore, it is useful to develop an alternative and simple method for prediction of tR accurately.


In the present study, based on the thermodynamic theory of HPLC, a method named linear calibration using two reference substances (LCTRS) was proposed. The method includes three steps, procedure of two points prediction, procedure of validation by multiple points regression and sequential matching. The tR of compounds on a HPLC column can be calculated by standard retention time and linear relationship.


The method was validated in two medicines on 30 columns.


It was demonstrated that, LCTRS method is simple, but more accurate and more robust on different HPLC columns than RR method. Hence quality standards using LCTRS method are easy to reproduce in different laboratories with lower cost of reference substances.


RP-HPLC Retention time Relative retention Linear calibration using two reference substances Multi-component analysis Traditional Chinese medicines


Multi-components analysis is an effective strategy for quality control of traditional Chinese medicines (TCMs), which have complex chemical profiles. But the classic external standard method was severely confined in its application due to the high cost of reference substances. As a consequence, substitute reference substance methods such as extractive reference substance (ERS) method and single standard to determine multi-components (SSDMC) method for overall quality control of TCMs have emerged, and widely used in Chinese pharmacopoeia 2015 edition, the United States Pharmacopoeia (USP39-NF34) and literatures [110]. In general, ERS method provides only one reference chromatogram in the pharmacopoeias, instructions of ERS and literatures. But there are hundreds of brands of C18 columns in the market. It means that the reference chromatogram may be different from the actual chromatogram. Due to the column types and other various factors, the error between measured retention time (tR) and predicted tR by the relative retention (RR) method cannot be ignored sometimes.

In order to improve the reproducibility of chromatographic separation and RR, the method of classification of C18 columns has been proposed [1115]. The columns were divided into three types: A, B and EP. Although the same type of columns was used to repeat the analytical methods, the differences in the performance and the separation effects were still large. And then the methods for selecting columns with equivalent selectivity, such as the USP approach [16], the PQRI approach [17, 18] and Katholieke Universiteit Leuven column classification system [1921] were proposed. Take PQRI approach [17, 18] as an example, hydrophobicity (H), steric interaction (S), hydrogen-bond acidity (A), hydrogen-bond basicity (B) and ion-exchange capacity (C), were used to describe the performance of the column. And the similarity between a column and the reference column was calculated by these five parameters. When the similarity was less than three, the two columns were regarded to be equivalent. Using the equivalent column, the reproducibility of separation and RR could be improved to some extent. However, in addition to column, many other factors also have great influences on the chromatograms, such as the dead volume of chromatographic system, the different structure of analytes, the complexity of the chromatographic conditions, and so on. Therefore, it is necessary to develop a method that takes all aforementioned factors into account to reduce the prediction error of the tR.

According to the thermodynamic theory of liquid chromatography, there is a linear relationship between the tR of the compounds on two different HPLC systems (including chromatographs and columns) [22]. For better understanding, the pdf of reference 22 (in Chinese, Additional file 1) and the English version of reference 22 (only the section of theory was translated, Additional file 2) are provided. Combined with the above principle and previous studies [2325], a novel method using two reference substances for predicting HPLC tR has been proposed (linear calibration by two reference substances, LCTRS). The StR (arithmetic average of tR for the same compound on different HPLC system under the same chromatographic conditions) is used as the reference value, and the linear regression is used as the basic algorithm for tR prediction. In this study, the method was validated in two medicines on 30 C18 columns. Compared with the RR method, LCTRS method is proved to be more accurate, and more robust on different HPLC columns. Hence, it provides a good prospective application in quantification of multi-components in TCMs as well as related substances in pharmaceutical chemicals.


The Minimum Standards of Reporting Checklist contains details of the experimental design, and statistics, and resources used in this study (Additional file 3).

Instruments and reagents

Waters e2695 HPLC (2998PDA detector), Agilent 1260 HPLC (DAD detector), and Shimadzu LC-2010A HT HPLC (UV–Vis detector) were used. Matlab software was provided by Math Works Inc. USA. 30 C18 columns (shown in Table 1), from 13 manufacturers, included A, B, and EP types were used. And most columns belong to type B according to the previous study [1115]. According to the PQRI approach [17, 18] and using the data from the USP website (, the similarity of columns were calculated using col1 as the reference column. The similarity (0–13.18) showed that the differences among the columns were large, which indicated that the selected columns are in a wide range and have good representative trait.
Table 1

tR (min) of four saponins in Paridis on different columns



Chonglou saponin VII

Chonglou saponin VI

Chonglou saponin II

Chonglou saponin I


Discovery C 18 a

21.234 ± 0.021

23.001 ± 0.006

32.773 ± 0.007

35.118 ± 0.004


Discovery C 18 b

16.555 ± 0.010

17.989 ± 0.006

27.161 ± 0.011

29.154 ± 0.002


Xbridge C 18 a

21.101 ± 0.004

23.070 ± 0.027

32.483 ± 0.021

34.963 ± 0.006


BDS HypersilC 18 a

21.014 ± 0.011

22.898 ± 0.018

32.679 ± 0.032

35.170 ± 0.036


Inertsil ODS-2a

22.132 ± 0.009

24.502 ± 0.004

33.176 ± 0.017

35.936 ± 0.003


Kromasil C 18 a

21.276 ± 0.012

23.693 ± 0.016

32.618 ± 0.006

35.929 ± 0.007


Luna C18(2)a

20.760 ± 0.018

23.362 ± 0.008

30.941 ± 0.011

33.813 ± 0.003


Luna C18(2)b

16.551 ± 0.008

18.865 ± 0.023

25.923 ± 0.016

28.310 ± 0.011


Inertsil ODS-3b

17.225 ± 0.006

19.640 ± 0.005

26.754 ± 0.013

29.471 ± 0.005


Alltima C 18 a

20.856 ± 0.008

23.752 ± 0.021

31.687 ± 0.011

34.872 ± 0.004


Symmetry C 18 a

21.016 ± 0.015

22.076 ± 0.016

32.470 ± 0.023

35.476 ± 0.010


Gemini C 18 a

21.300 ± 0.014

23.756 ± 0.007

31.599 ± 0.003

34.337 ± 0.015



21.076 ± 0.012

23.828 ± 0.004

31.627 ± 0.006

34.695 ± 0.001


Zorbax Extend-C 18 a

17.201 ± 0.022

19.731 ± 0.003

27.525 ± 0.009

30.504 ± 0.005


Sunfire C 18 a

21.652 ± 0.013

24.065 ± 0.024

32.375 ± 0.010

35.197 ± 0.006


Sunfire C 18 b

17.501 ± 0.006

19.571 ± 0.018

27.401 ± 0.004

29.826 ± 0.010


Nucleosil C18 HDa

21.452 ± 0.013

23.735 ± 0.006

32.747 ± 0.019

35.513 ± 0.011


ODS Hypersila

19.669 ± 0.010

21.583 ± 0.022

30.361 ± 0.007

32.657 ± 0.005



20.092 ± 0.013

22.560 ± 0.010

29.568 ± 0.006

32.218 ± 0.005


Spherisorb ODS2a

18.684 ± 0.006

21.334 ± 0.021

28.580 ± 0.008

31.385 ± 0.004


Zorbax SB-C 18 a

18.438 ± 0.015

20.986 ± 0.010

28.085 ± 0.008

30.849 ± 0.007


DiamonsilC 18 a

22.082 ± 0.021

25.078 ± 0.005

32.475 ± 0.020

35.721 ± 0.023


DiamonsilC 18 b

16.769 ± 0.007

19.154 ± 0.008

26.240 ± 0.011

28.935 ± 0.004


Diamonsil C18(2)a

19.296 ± 0.016

22.819 ± 0.004

31.192 ± 0.019

34.710 ± 0.009


Kromasil Eternity C 18 a

19.732 ± 0.017

22.099 ± 0.018

29.535 ± 0.004

32.153 ± 0.021


Shim-pack VP-ODSb

19.015 ± 0.007

21.059 ± 0.014

29.401 ± 0.002

31.845 ± 0.001


Agilent HC-C 18 a

23.492 ± 0.247

25.839 ± 0.250

34.884 ± 0.344

37.637 ± 0.330


Agilent TC-C 18 a

20.609 ± 0.003

22.574 ± 0.005

30.806 ± 0.003

32.783 ± 0.002


Venusil MP C 18 b

17.970 ± 0.007

20.468 ± 0.011

27.287 ± 0.031

29.985 ± 0.007


Nucleosil C18 ABa

18.348 ± 0.012

20.209 ± 0.017

29.233 ± 0.011

31.893 ± 0.014


StR (Average tR)





a4.6 mm × 250 mm × 5 μm

b4.6 mm × 150 mm × 5 μm

Reference substances of psoralen, isopsoralen, Chonglou saponin I, Chonglou saponin II, Chonglou saponin VI, Chonglou saponin VII, ethinylestradiol, and herbal reference substances including Psoraleae Fructus (Psoraleae) and Paridis Rhizome (Paridis), were supplied by the National Institutes for Food and Drug Control, China. Methanol, acetonitrile, and phosphoric acid were HPLC graded and supplied by the Fischer Company, USA. Ammonium nitrate (analytical grade) was supplied by Beijing Chemical Works. Water was prepared by Milli-Q system, Millipore Company, USA.

Preparation of sample solution

Psoraleae [26]: weigh 0.5 g of the powder and place it in a 50-mL stopper conical flask, then add 25 mL of ethanol. Sonicate the mixture for 30 min and centrifuge for 5 min. Filter the supernatant through a 0.45-μm PTFE filter.

Paridis [27]: weigh 0.5 g of the powder to a stopper conical flask, add 25 mL of ethanol, heat under reflux on a water bath for 30 min, cool and filter the supernatant through a 0.45-μm PTFE filter.

Chromatographic conditions

Psoraleae [26]: mobile phase A was water and mobile phase B was methanol. The elution procedure was shown below: 0–20 min, 50%B→70%B; 20–45 min, 70%B→85%B; 45–50 min, 85%B→90%B; 50–60 min, 90%B. Detection wavelength was 308 nm. Paridis [27]: mobile phase A was water and mobile phase B was acetonitrile. The elution procedure was shown below: 0–40 min, 30%B→60%B; 40–50 min, 60%B→30%B; 50–60 min, 30%B. Detection wavelength was 203 nm. Column temperatures were both set at 30 °C and flow rate was 1.0 mL/min.


HPLC chromatogram of samples

The typical chromatograms of Psoraleae and Paridis were shown in Fig. 1. The peaks were mainly identified by the reference substances. For those peaks without reference substances, UV–Vis spectrum and mass spectrum were used for identification.
Fig. 1

HPLC chromatograms of samples

Standard retention time (StR)

Under the same chromatographic conditions, measured retention time (tRmea) of the four saponins in Paridis on different chromatographic systems (which includes HPLC instruments and columns, hereinafter referred to as columns due to the differences of tR mainly caused by columns) were shown in Table 1. The arithmetic average of tR for the same compound on different columns is called StR, formula (1). Just like RR, StR is the reference value for calculating the predicted retention time (tRpre) of analyte in the samples. Theoretically, under the same chromatographic condition, the RR calculated by different columns is constants, but StR is not. It will be discussed in Section "Minimum number of columns for StR calculation" that the advantages of using StR was better than tR of any single column. In this paper, the deviation (ΔtR) of tRmea and tRpre (formula 2) was used to evaluate the merits and defects of RR method and LCTRS method.
$$ St_{R} = \mathop \sum \limits_{i = 1}^{n} t_{Ri} /n \quad (n \ge 1) $$
$$ \Delta t_{R} = \left| {t_{R} mea - t_{R} pre} \right| $$

Linear principle of LCTRS

According to the chromatographic thermodynamic theory, Wang et al. proved that there was a linear relationship between the tR of the same compounds on different HPLC system (mainly considered as columns) under the same chromatographic conditions [22], as expressed in formula (3) and Fig. 2a, b.
$$ t_{R} coli = a \times t_{R} colj + b $$
Fig. 2

Linear fitting results of Psoraleae and Paridis, code No. is the same as that in Fig. 1

Since formulas (1) and (3) are both linear, thus there is a linear relationship between tR and StR for each compound, as shown in formula (4) and Fig. 2c, d. It is noteworthy that the correlation coefficient of the linear regression is higher than that shown in Fig. 2a, b.
$$ t_{R} coli = a \times St_{R} + b $$

Minimum number of columns for StR calculation

Theoretically, tR on any column can be used as reference value for linear fitting. But the ΔtR calculated with random column were instable. Thus, the reasonable number of columns for StR calculation was thoroughly investigated by random sampling. StR was calculated based on 1, 5, 10, 15, 20, 25, and 30 columns combined with non-replicate random sampling times of 30, 100, 100, 100, 100, 100, and 1, respectively. The value of StR with tRmea on 30 columns was used to fit multiple point linear equation. The averages of ΔtR (average ± standard deviation) were calculated, as shown in Fig. 3. For both medicines, the prediction deviation was reduced with increasing number of columns. However, the prediction accuracy will not be significantly improved when the number of columns reaches five, which is considered as a low-cost and reasonable limit. It is recommended to choose five–fifteen columns for StR calculation.
Fig. 3

Prediction results of different number of columns for StR calculation

Even for the columns with same type of packing material, there are still some differences among the column stationary phase, packing techniques and errors in the process of chromatographic analysis. Those differences will cause deviation of tR. The physical explanation of StR calculation was to evenly mix and refill the stationary phase of the columns selected. Because of reducing the random and system errors, the prediction result was accurate and robust.

Procedure of two points prediction

For RR method, only one compound was chosen as reference compound (reference substance required), and RR of all other compounds were used as reference value for calculating tRpre. For LCTRS method, two compounds were chosen as reference compounds (reference substances required), and StR of all other compounds were used as reference value for calculating tRpre. The reference compounds, the value of RR and StR were shown in Tables 2 and 3.
Table 2

RR and reference compound for RR method


Peak 1

Peak 2

Peak 3

Peak 4

Peak 5

Peak 6

Peak 7

Peak 8

Peak 9

Peak 10

Peak 11

Psoraleae (n = 23)












Paridis (n = 30)






aReference compound

Table 3

StR (min) and reference compound for LCTRS method


Peak 1

Peak 2

Peak 3

Peak 4

Peak 5

Peak 6

Peak 7

Peak 8

Peak 9

Peak 10

Peak 11

Psoraleae (n = 23)












Paridis (n = 30)






aReference compound

Take Paridis as an example. First of all, reference substances solution of two reference compounds (Chonglou saponin VII and Chonglou saponin I) and sample solution were performed on a C18 column (col4: BDS Hypersil C18). The tRmea (21.014 and 35.170 min) of two reference compounds in the sample solution were obtained by the reference substances solution (Fig. 4a). Then two points, Chonglou saponin VII (19.803, 21.014) and Chonglou saponin I (33.035, 35.170), could be determined in the coordinate using StR as abscissa and tRmea as ordinate. Based on the two points, the following linear equation was given: y = 1.0698x − 0.1719 (Fig. 4b). Taking StR of analytes (Chonglou saponin VI and Chonglou saponin II) into equation, the tRpre of Chonglou saponin VI (23.481 min) and Chonglou saponin II (32.263 min) were attained, respectively. Finally, in the chromatogram of the sample solution, the corresponding peaks of Chonglou saponin VI and Chonglou saponin II can be found within the range of tRpre ± tRW (tRW is abbreviation of tR window, in this case is 0.6 min), as shown in Fig. 4c. It can be seen that ΔtR of analytes calculated by prediction of two points were 0.583 min and 0.416 min (The tRmea of analytes were 22.898 min and 32.679 min).
Fig. 4

Flow chart of LCTRS (Paridis, code No. is the same as that in Fig. 1)

Procedure of multiple points regression

After assignment of the peaks of analytes in the sample solution by prediction of two points regression, the tRmea of those peaks should be validated by multiple points regression. In this procedure, tRmea and StR of reference compounds and analytes were used to fit a multiple points linear regression: Y = 1.1038x − 1.1075 (Fig. 4d). Taking StR of analytes (Chonglou saponin VI and Chonglou saponin II) into new equation, the new tRpre of Chonglou saponin VI (23.297 min) and Chonglou saponin II (32.358 min) were calculated. If ΔtR of all analytes were less than the given tRL (tRL: tR limit, in this case is 0.5 min), the prediction was success, otherwise failure (Fig. 4e). In this case, ΔtR were 0.399 and 0.321 min, respectively. The step of validation by multiple points was based on the principle of stepwise linear regression, which can further improve the prediction accuracy.

The purpose of setting that the tRW is larger than tRL is to increase the amount of suitable columns and to improve the accuracy of prediction. Generally, the recommended ranges of tRW and tRL are 0.8–2.0 and 0.5–1.5 min, respectively. If necessary, the values can be adjusted in accordance with different samples under different chromatographic conditions. If the ΔtR of some compounds are large, their tRW and tRL can be set individually.

Sequential matching rule

If the tR of two peaks are too close, e.g. less than 2 min, there would be a mistake for peak matching by the least ΔtR rule. Take Psoraleae for example, as shown in Fig. 5a, peak #6 was assigned to peak A in the sample solution on col6 (Kromasil C18) with a small ΔtR of 0.515 min. However, peak #5 was not found and peak B in the sample solution was not matched. When tRW was set as 1.2 min, tRmea of peak A was within the window of tRpre of peak #5. tRmea of the peaks A and B would both fall into the window of tRpre of peak #6. Because of the existence of one common peak (peak A), peaks #5 and 6 should be treated as peak series for sequential matching. That is, the earlier tRpre will be matched to the peak with the earlier tRmea. Although ΔtR of peak #6 increased to 1.036 min, the match results were correct, as shown in Fig. 5b. This rule can be further applied to multiple-peak series, which has a close tR.
Fig. 5

Advantage of sequential matching (Psoraleae, code No. is the same as that in Fig. 1)

Comparison between LCTRS method and RR method

The comparison among unadjusted RR method, adjusted RR method (dead time was measured by ammonium nitrate as probe compound), prediction by two points, and validation by multiple points was summarized in Tables 4 and 5. The results showed that the unadjusted RR method and adjusted RR method were similar, their prediction accuracy were bigger and suitable for less positive columns. But the prediction deviation was reduced and the number of positive columns was increased by LCTRS method. The best was validation by multiple points which was based on the prediction by two points.
Table 4

Comparison result by four methods (Psoraleae)


Maximum of ΔtR/min

Average of ΔtR/min

Number of positive columnsa

Unadjusted RR




Adjusted RR




Prediction by two points




Validation by multiple points




aThe columns which meet the following requirements are called positive column, (1) the resolution of peaks meets the requirements; (2) ΔtR of all pending test compounds are no more than tRL (for Psoraleae is 1.2 min)

Table 5

Comparison of prediction result in four methods (Paridis)


Maximum of ΔtR/min

Average of ΔtR/min

Number of positive columnsa

Unadjusted RR




Adjusted RR




Prediction by two points




Validation by multiple points




atRL = 0.5 min

Exclusion of column and compound by linear fitting

Nonlinear shift of tR for a compound on different columns could be caused either by different column packing materials and use of other packing techniques, or by the different compound structure. In order to exclude the columns and compounds with relatively large nonlinear shift, linear fitting of tRmea and StR were performed. The following rules were used to identify the outlier column and compound. (1) In a regression scatter plot, the compounds obviously deviated from a regression line (the correlation coefficient is usually less than 0.99). (2) ΔtR was usually larger than 1–2 min. The excluded columns and compounds would not be used for StR calculating.

For Psoraleae: no obvious nonlinear deviation was observed of all 11 compounds. 23 columns met the requirements (the average of correlation coefficient was 0.9989). The outlier columns were col2, 8, 12 (Fig. 6a), 16, 19, 20 and 30. For Paridis: no obvious nonlinear deviation was observed of four saponins. All 30 columns met the requirement with average correlation coefficient of 0.9993.
Fig. 6

Outlier column (a) and Outlier compounds (b), code No. is the same as that in Fig. 1

In order to simulate tR of compound with large structural difference, reference substances solution of ethinylestradiol mixing with four Chonglou saponins were used to measure tR of those five compounds on 30 columns. Nonlinear shift of ethinylestradiol was observed on col1 (Fig. 6b), 2–6, 8, 11, 15–18, 26–28 and 30. It appears that the HPLC retention behaviors of ethinylestradiol and four Chonglou saponins were significantly different on this chromatographic condition. It further indicated that the classification and similarity evaluation of columns should be based on the characteristics of columns as well as analytes.

If the outlier compounds cannot be excluded. The following approaches could be used: (1) specify one or more suitable columns; (2) provide reference substances for those compounds; (3) use UV–Vis spectrum and/or mass spectrum for assistant peak identification.

Selection of two reference compounds

Ideally there should be no difference in selecting any of the two compounds as reference compounds. However, because of the difference of HPLC instruments, columns, compounds structure, complexity of elution condition, and accidental error of analysis, different selection of reference compound pairs will make differences. In order to find out the rule for reference compounds selection, each combination of possible reference compound pairs for the two medicines was studied. The average of ΔtR corresponding to each reference compound pair were calculated and shown in Fig. 7. It can be seen that, for the two medicines, the ΔtR of prediction by two points step would be decreased with increasing coverage of tR [as shown in formula (5). The coverage of tR is a reflex of the relative position of the two reference compounds. The first compound is at one end (with smaller tR), the last compound is at the other end (with bigger tR). If the coverage is high, the two reference compounds are near both ends, otherwise they are in the middle or near the same end]. The coverage corresponding to the smallest ΔtR was 80–100%.The results of Psoraleae (Fig. 7a) showed the advantage of choosing reference compounds with smaller linear deviation, when the coverages of tR were similar. Therefore, the optimized reference compounds for Psoraleae were peak #2 and peak #8, rather than peak #1 and peak #11 which had a maximum coverage of tR but with more deviation from linearity. The selection rule decreases the randomness of choosing reference compounds and the amount of calculation (or the ΔtR of all possible reference compound pairs will be calculated every time). The accurate and simple selection procedures were as follows. Firstly, Select two reference compound pairs with large tR coverage (80–100%). Secondly, exclude compounds with large linear deviation based on the linear fitting results. Thirdly, calculate the ΔtR of the rest reference compound pairs and select reference compound pairs with the smallest ΔtR.
$$ Coverage \; of \; t_{R} = \frac{{t_{R2} - t_{R1} }}{{t_{Rlast} - t_{Rfirst} }} $$
tR2 is tR (or StR) of second reference compound; tR1 is tR (or StR) of first reference compound; tRlast is tR (or StR) of last compound; tRfirst is tR (or StR) of first compound.
Fig. 7

Selection of two reference compound (abbreviated as RC)

In summary, the establishment procedures of LCTRS were as follows. Firstly, select five–fifteen different brands of C18 columns, and record the HPLC chromatograms of reference substances and sample on all columns. Secondly, calculate initial StR by using all columns, and perform linear fitting of tR on each column with StR. Exclude outlier columns and compounds, and recalculate final StR using remaining columns. Finally, select two reference compounds with large tR coverage and low linear deviation.


Advantages of LCTRS method

According to the study of Wang et al. [22], tR of the compounds on different HPLC system follows the linearity principle. The RR method can be regarded as external standard one point method, which means the regression line is forced to pass origin. However, most of the linear equations have intercepts, which is why the deviation of unadjusted RR method was large. For considering the dead time, adjusted RR method should be better than unadjusted RR method in theory. But the probe compound for dead time measurement would be interacted with mobile phases and stationary phases of the columns. The interaction would increase the error in dead time measurement. So the prediction accuracy of this method was not improved in practice. For prediction by two points and validation by multiple points, dead time, gradient delay, volume exclusion effect of stationary phase, retention behavior of homologous compounds and so forth, were fully considered. Thus, the prediction accuracy was significantly improved. Stepwise linear regression was used in the validation by multiple points step, which further improved the prediction effect.

Compatibility of LCTRS method and RR method

Both LCTRS method and RR method are equivalent in mathematics. Formulas can be expressed in the same form. In the LCTRS, calibrated retention (CR) is defined as the ratio of StR of analytes to reference compounds, as shown in formula (6). Different from RR, CR is based on statistics of StR. Thus, its prediction accuracy was equal to LCTRS (only equal to prediction by two points).
$$ CR = RR = \frac{{t_{Ri} - t_{R1} }}{{t_{R2} - t_{R1} }} $$
where tRi is StR of analytes in CR, or tR of analytes in RR; tR1 is StR of the first reference compound in CR, or dead time in adjusted RR, or zero in unadjusted RR; tR2 is StR of the second reference compound in CR, or tR of reference compound in RR.


A new method for tR prediction of HPLC chromatographic peaks was proposed. 16 compounds in two medicines under isocratic or gradient elution conditions were tested through three brands of HPLC instruments with 30 different brands of C18 columns. It is demonstrated that the method is simple, accurate, and robust for more HPLC columns. Furthermore, the calculation approach of StR and the selection rule of the two reference compounds were discussed.

Both multi-components analysis in TCMs and determination of related substances in pharmaceutical chemicals need lots of reference substances for peak identification. But it may be not affordable for routine analysis and research using all reference substances. LCTRS is a simple and low-cost alternative method for peak identification. Compared with RR method, it need one more reference substance but is more accurate and suitable for more HPLC columns. LCTRS method provides a good prospective application for overall quality evaluation of TCMs and impurities analysis in pharmaceutical chemicals.



traditional Chinese medicines


relative retention


linear calibration using two reference substances


retention time


single standard to determine multi-components


arithmetic average of tR


deviation of tR


predicted retention time


measured retention time


calibrated retention


Authors’ contributions

LS and SCM conceived the study. LS and HYJ designed the study. LS and LPY conducted the experiments. LS, RTT, MJW, LNL and TTZ analyzed the data. LS wrote the manuscript. RTT and SCM revised the manuscript. All authors read and approved the final manuscript.


We would like to thank Ms. Yu Pang from Center for Drug Reevaluation and Mr. Chao Zhang from Suzhou Institute for Drug Control for the useful suggestions on this paper.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

All data generated or analysed during this study are included in this published article [and its supplementary information files].


This project was financially supported by Youth Development Research Foundation of NIFDC (2013WA8), National Natural Foundation of China (81303214) and National Major Scientific and Technological Special Project for “Significant New Drugs Development” (2014ZX09304307-002).

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis 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.

Authors’ Affiliations

National Institutes for Food and Drug Control
Xinjiang Institute for Food and Drug Control
Chemmind Technologies Co., Ltd.
Huainan Food and Drug Inspection Center


  1. Gao XY, Jiang Y, Lu JQ, Tu PF. One single standard substance for the determination of multiple anthraquinone derivatives in rhubarb using high-performance liquid chromatography-diode array detection. J Chromatogr A. 2009;1216(11):2118–23.View ArticlePubMedGoogle Scholar
  2. Hou JJ, Wu WY, Da J, Yao S, Long HL, Yang Z, Cai LY, Yang M, Liu X, Jiang BH, Guo DA. Ruggedness and robustness of conversion factors in method of simultaneous determination of multi-components with single reference standard. J Chromatogr A. 2011;1218(33):5618–27.View ArticlePubMedGoogle Scholar
  3. Hou JJ, Wu WY, Liang J, Yang Z, Long HL, Cai LY, Fang L, Wang DD, Yao S, Liu X, Jiang BH, Guo DA. A single, multi-faceted, enhanced strategy to quantify the chromatographically diverse constituents in the roots of Euphorbia kansui. J Pharm Biomed Anal. 2014;88(88C):321–30.View ArticlePubMedGoogle Scholar
  4. Wu X, He J, Xu H, Bi K, Li Q. Quality assessment of Cinnamomi Ramulus by the simultaneous analysis of multiple active components using high-performance thin-layer chromatography and high-performance liquid chromatography. J Sep Sci. 2014;37(18):2490–8.View ArticlePubMedGoogle Scholar
  5. Liao H, Li Q, Liu R, Liu J, Bi K. Fingerprint analysis and multi-ingredient determination using a single reference standard for Saposhnikoviae Radix. Anal Sci. 2014;30(12):1157–63.View ArticlePubMedGoogle Scholar
  6. Zhong JS, Wan JZ, Liu YP, Ding WJ, Wu XF, Xie ZY. Simultaneous HPLC quantification of seven chromones in Aloe barbadensis Miller using a single reference standard. Anal Methods. 2014;6(12):4388–95.View ArticleGoogle Scholar
  7. Du YY, Li Q, Liu JJ, Yin YD, Bi KS. Combinative method using multi-components quantitation by single reference standard and HPLC fingerprint for comprehensive evaluation of Rhodiola crenulatav H. Ohba. Anal Methods. 2014;6(15):5891–8.View ArticleGoogle Scholar
  8. Da J, Cheng CR, Yao S, Long HL, Wang YH, Khan LA, Li YF, Wang QR, Cai LY, Jiang BH, Liu X, Wu WY, Guo DA. A reproducible analytical system based on the multi-component analysis of triterpene acids in Ganoderma lucidum. Phytochemistry. 2015;114:146–54.View ArticlePubMedGoogle Scholar
  9. Lu WY, Liu YG, Yang HS, Sheng Y, Shi HM, Yu LL. Simultaneous HPLC quantification of five major triterpene alcohol and sterol ferulates in rice bran oil using a single reference standard. Food Chem. 2014;148(4):329–34.View ArticlePubMedGoogle Scholar
  10. Zhang YW, Li Q, Lv CX, Liu XJ, Chen XH, Bi KS. Simultaneous determination of four active components in Alisma orientale (Sam.) Juz. by HPLC–DAD using a single reference standard. J Pharmaceut Anal. 2015;5(2):85–92.Google Scholar
  11. Gilroy JJ, Dolan JW, Snyder LR. Column selectivity in reversed-phase liquid chromatography: IV type-B alkyl-silica columns. J Chromatogr A. 2003;1000:757–78.View ArticlePubMedGoogle Scholar
  12. Gilroy JJ, Dolan JW, Carr PW, Snyder LR. Column selectivity in reversed-phase liquid chromatography: V. Higher metal content (type-A) alkyl-silica columns. J Chromatogr A. 1026;2004(1026):77–89.Google Scholar
  13. Wilson NS, Gilroy J, Dolan JW, Snyder LR. Column selectivity in reversed-phase liquid chromatography: VI. Columns with embedded or end-capping polar groups. J Chromatogr A. 2004;1026(1–2):91–100.View ArticlePubMedGoogle Scholar
  14. Visky D, Vander Heyden Y, Iványi T, Baten P, De Beer J, Kovács Z, Noszál B, Dehouck P, Roets E, Massart DL, Hoogmartens J. Characterisation of reversed-phase liquid chromatographic columns by chromatographic tests: rational column classification by a minimal number of column test parameters. J Chromatogr A. 2003;1012(1):11–29.View ArticlePubMedGoogle Scholar
  15. Dehouck P, Visky D, Vander Heyden Y, Adams E, Kovács Z, Noszál B, Massart DL, Hoogmartens J. Characterisation of reversed-phase liquid-chromatographic columns by chromatographic tests: comparing column classification based on chromatographic parameters and column performance for the separation of acetylsalicylic acid and related compounds. J Chromatogr A. 2004;1025(2):189–200.View ArticlePubMedGoogle Scholar
  16. Bidlingmeyer B, Chan CC, Fastino P, Henry R, Koerner P, Maule AT, Marques MRC, Neue U, Ng L, Pappa H, Sander L, Santasania C, Snyder L, Wozniak T. HPLC column classification. Pharm Forum. 2005;31(2):637–45.Google Scholar
  17. Dolan JW, Maule A, Bingley D, Wrisley L, Chan CC, Angod M, Lunte C, Krisko R, Winston JM, Homeier BA, McCalley DV, Snyder LR. Choosing an equivalent replacement column for a reversed-phase liquid chromatographic assay procedure. J Chromatogr A. 2004;1057(1–2):59–74.View ArticlePubMedGoogle Scholar
  18. Snyder LR, Dolan JW, Carr PW. The hydrophobic-subtraction model of reversed-phase column selectivity. J Chromatogr A. 2004;1060(1–2):77–116.View ArticlePubMedGoogle Scholar
  19. Visky D, Haghedooren E, Dehouck P, Kovács Z, Kóczián K, Noszál B, Hoogmartens J, Adams E. Facilitated column selection in pharmaceutical analyses using a simple column classification system. J Chromatogr A. 2006;1101(1–2):103–14.View ArticlePubMedGoogle Scholar
  20. Kóczián K, Haghedooren E, Dragovic S, Noszál B, Hoogmartens J, Adams E. Column selection for pharmaceutical analyses based on a column classification using four test parameters. J Pharm Biomed Anal. 2007;44(4):894–905.View ArticlePubMedGoogle Scholar
  21. Dragovic S, Haghedooren E, Németh T, Palabiyik IM, Hoogmartens J, Adams E. Evaluation of two approaches to characterise liquid chromatographic columns using pharmaceutical separations. J Chromatogr A. 2009;1216(15):3210–6.View ArticlePubMedGoogle Scholar
  22. Wang LX, Xiao HB, Liang XM. A new method to improve the reproducibility of retention time on reversed phase C18 columns in different laboratories. Chin J Anal Chem. 2003;31(10):1232–6.Google Scholar
  23. Sun L, Jin HY, Pang Y, Ma SC. Two reference substances for determination of multiple components (1): linear calibration using two reference substances for identification of chromatographic peaks. Chin J Pharm Anal. 2013;33(8):1424–30.Google Scholar
  24. Sun L, Pang Y, Jin HY, Liu LN, Ma SC. Two reference substances for determination of multiple components II influence of detection wavelength selection on quantitative analysis. Chin J Pharm Anal. 2013;33(9):1578–86.Google Scholar
  25. He B, Yang SY, Zhang Y. A new method of calibration and positioning in quantitative analysis of multi components by single marker[J]. Acta Pharm Sin. 2012;12:1653–9.Google Scholar
  26. Department of Health Hong Kong Special Administrative Region The People’s Republic of China. Hong Kong Chinese Materia Medica Standards, vol. 3, Government Logistics Department. Beijing; 2012. p. 101.Google Scholar
  27. National Commission of Chinese Pharmacopoeia. Pharmacopoeia of People’s Republic of China, vol. I, China Medical Science Press. Beijing; 2015. p. 260.Google Scholar


© The Author(s) 2017