Next Article in Journal
Analysis of Corticosteroid-Induced Glaucoma Using the Japanese Adverse Drug Event Reporting Database
Next Article in Special Issue
Development and Validation of a Novel UHPLC-MS/MS Method for the Quantification of Plinabulin in Plasma and Its Application in a Pharmacokinetic Study with Leukopenic Rats
Previous Article in Journal
Novel Thiazolidine-2,4-dione-trimethoxybenzene-thiazole Hybrids as Human Topoisomerases Inhibitors
Previous Article in Special Issue
Determination of Oxaliplatin by a UHPLC-MS/MS Method: Application to Pharmacokinetics and Tongue Tissue Distribution Studies in Rats
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of the Analytical Procedure Lifecycle Concept to a Quantitative 1H NMR Method for Total Dammarane-Type Saponins

1
Pharmaceutical Informatics Institute, College of Pharmaceutical Sciences, Zhejiang University, Hangzhou 310058, China
2
State Key Laboratory of Component-Based Chinese Medicine, Innovation Center, Zhejiang University, Hangzhou 310058, China
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Pharmaceuticals 2023, 16(7), 947; https://0-doi-org.brum.beds.ac.uk/10.3390/ph16070947
Submission received: 29 May 2023 / Revised: 27 June 2023 / Accepted: 27 June 2023 / Published: 29 June 2023
(This article belongs to the Special Issue Analytical Techniques in the Pharmaceutical Sciences 2023)

Abstract

:
Dammarane-type saponins (DTSs) exist in various medicinal plants, which are a class of active ingredients with effects on improving myocardial ischemia and immunomodulation. In this study, a quantitative 1H NMR method of total DTSs in herbal medicines was developed based on the analytical procedure lifecycle. In the first stage (analytical procedure design), the Ishikawa diagram and failure mode effects and criticality analysis were used to conduct risk identification and risk ranking. Plackett–Burman design and central composite design were used to screen and optimize critical analytical procedure parameter. Then, the method operable design region was obtained through modeling. In the second stage (analytical procedure performance qualification), the performance of methodological indexes was investigated based on analytical quality by design. As examples of continued procedure performance verification, the method was successfully applied to determine the total DTSs in herbal pharmaceutical preparations and botanical extracts. As a general analytical method to quantify total DTSs in medicinal plants or pharmaceutical preparations, the developed method provides a new quality control strategy for various products containing dammarane-type saponin.

1. Introduction

Dammarane-type saponins (DTSs), classified as tetracyclic triterpene saponins, are widely distributed in Araliaceae and Cucurbitaceae herbs, such as Ginseng, American Ginseng, Panax notoginseng, and Gynostemma. Thus far, more than 760 DTSs have been reported from more than 130 plant species [1]. Studies have revealed that DTSs possess a diverse range of effects including antitumor, anti-inflammatory, myocardial ischemic, immunomodulatory, hypoglycemic, anti-shock, hepatoprotective, sedative, and tranquilizing properties [2], so DTSs are considered as crucial active pharmaceutical ingredients (APIs) in herbal preparations. In the standards of relevant preparations, the content of total saponins is specified as a critical quality evaluation indicator. The glacial acetic acid–perchloric acid–UV chromogenic method is a classical approach for determining total saponins content. However, this method is laborious, necessitates the use of highly corrosive reagents, and can only provide fluctuating estimated outcomes, rather than absolute molar concentrations [3]. Several studies have reported the simultaneous determination of up to 10 different DTSs using HPLC or UPLC methods [4,5,6]. Nonetheless, it has been shown that some preparations contain more than 50 kinds of DTSs [7]. Therefore, it is imperative to make further advancements in methods development for the analysis of DTSs that are more efficient, accurate, and comprehensive.
Li et al. [3] developed a proton quantitative nuclear magnetic resonance (1H qNMR) method for the determination of total ginsenosides. Although this method has limitations in its scope, as it cannot be directly applied to other classes of DTSs or other herbal species, it still demonstrates the potential of qNMR methods for the analysis of DTSs. 1H qNMR spectroscopy is a non-destructive and rapid analytical technique with quantitative ability of high accuracy and good reproducibility. Due to the broad application prospect of qNMR in the analysis of complex systems, this study aims to develop a 1H qNMR method for the determination of DTSs in herbal medicines. In recent years, several standards and guidelines have recommended the application of lifecycle management to the analytical methodology development procedure in order to obtain a more robust parameter space [8,9].
The concept of lifecycle management was mentioned in 1975 for the first time in the scientific literature [10] and has been introduced into the field of chemical analysis in the 21st century [11]. In 2021, the United States Pharmacopeia (USP-2021) included the General Chapter <1220> “The Analytical Procedure Lifecycle (APLC)” formally, which addresses the application of lifecycle management to the analytical procedure [8]. APLC integrates and extends the steps of analytical method development, validation, verification, transfer, and maintenance in traditional procedures, which are regarded as a continuous dynamic cycle to ensure valid data generation throughout its lifecycle. Analytical target profile (ATP) is a fundamental component of the lifecycle approach to analytical procedures, which determines the criteria for the application purpose and the expected performance characteristics of the analytical measurements [8]. There are three stages of the lifecycle approach to analytical procedures, namely analytical procedure design (APD), procedure performance qualification (APPQ), and continued procedure performance verification (CPPV). In the first stage (APD), quality risk management (QRM) means are applied to carry out risk analysis and risk ranking of parameters. Then, the design of experiments (DoE) methods are applied to optimize critical analytical procedure parameters (CAPPs) and obtain the method operable design region (MODR) based on the concept of analytical quality-by-design (AQbD). In the second stage (APPQ), verification experiments were conducted to confirm that the method will operate (in routine use) as intended and meets the previously defined ATP criteria [11]. In the last stage (CPPV), the established method is continuously monitored and validated in its daily application based on a set of analytical control strategies.
The attempts to apply APLC In pharmaceutical analysis-related studies have gradually increased in the past decade, but reports involving the analysis of herbal medicines are still relatively few [12]. Most of the analytical methods for herbal medicines developed by applying APLC and AQbD are chromatographic methods aiming at the separation and quantification of active ingredients of herbal medicines. For example, there are studies about the development of HPLC-UV/PDA methods for the analysis of flavonoids [13,14], saponins [15,16], terpenoids [17], quinones [9], and phenolic [18]; HPLC-MS methods to determine terpene lactones [19] and polyphenols [20]; HPLC-ELSD methods to quantify sugars and their derivatives [21], etc. In addition, other chromatographic methods have also been developed based on AQbD, such as capillary electrophoresis [22,23], supercritical fluid chromatography [24,25], gas chromatography [26,27], etc. However, the practice of applying APLC and AQbD to developing qNMR methods for herbal medicines has not been reported.
In this study, APLC and AQbD approaches were applied to develop a quantitative 1H NMR (1H qNMR) method for total DTSs in herbal medicines. Failure mode effects and criticality analysis (FMECA) was used for risk assessment. Plackett–Burman design (PBD) [28] and central composite design (CCD) were applied to screen and optimize CAPPs. The developed method has been successfully applied to quantify total DTSs in a liquid preparation, a solid preparation, and a production process intermediate, and has performed well in the validation of all the methodological indexes.

2. Results and Discussion

2.1. Risk Identification and Assessment Results Based on the Ishikawa Diagram and FMECA

In this study, risk factors that may affect the results of the 1H qNMR method were analyzed as potential failure factors of FMECA from five aspects, including sample preparation, spectra acquisition, data processing, environment, and instrument, which were summarized in an Ishikawa diagram (Figure 1A) [12]. Risk factors identified in the Ishikawa diagram were evaluated by FMECA. The detailed process of FMECA was shown in Section S4. The definition of failure mode for the 1H qNMR method as well as S, O, D, and RPN scores of each factor were listed in Tables S2 and S3. As shown in Figure 1B, the maximum value of RPN was 125, the minimum value is 1, and the standard deviation (SD) was 42.88. Considering the 95% confidence interval, the statistical mean of RPN was 20.91~55.55. Based on the RPN score, risk factors could be divided into three levels: low, medium, and high. The high-risk factors for the method should be investigated and optimized in the method development process.
From the RPN calculation results, it was clear that sample concentration (M), number of dummy scans (DS), number of scans (NS), sampling temperature (Temp), size of fid (TD), receiver gain (RG), relaxation delay (D1), and integration method were the high-risk factors that need to be further investigated and strictly controlled. Among these factors, only the integration method was a data processing parameter, while the others were experimental parameters. Since there was no interaction between experimental and data processing parameters, the integration method would be examined separately from the other high-risk factors. Based on literature research [29] and single-factor experimental comparisons, the “line fitting” method in the MestReNova software was selected as the integration method. The remaining seven parameters will be further screened and optimized by DoE as potential CAPPs.

2.2. DoE-Based Analytical Process Design (APD)

2.2.1. CAPPs Screening Results Based on PBD

PBD Screening Experiments

The total of 15 runs for screening experiments were obtained by PBD, which were used to investigate the effects of M (x1), DS (x2), NS (x3), TD (x4), RG (x5), D1 (x6), and Temp (x7) on the 1H qNMR method. Considering that the most important expected requirements in the ATPs are good accuracy and precision, and precision is mainly related to signal resolution and signal intensity, accuracy, SNR, and resolution were used as evaluation indexes of CAPPs (y). Parameters of screening experiments in PBD and the coding values of each factor were shown in Table 1.
The samples used for PBD experiments were obtained by mixing representative DTSs ginsenoside Rb1, ginsenoside Rg1, ginsenoside Re, notoginsenoside R1, notoginsenoside Fe, gypenoside XLIX, gypenoside XVII, and gypenoside XLVI. The above mixing DTSs were dissolved in CD3OD containing 1.46 mmol/L DMT to prepare samples with total DTSs concentrations of 5, 2.5, and 0.5 mmol/L, respectively. The 1H NMR spectra of all samples were recorded with a Bruker AVANCE III 500 MHz NMR spectrometer (Bruker Technologies GmbH, Germany, 5 mm BBO probe with Topspin workstation). A classical water suppression pulse ZGPR was used to optimize the parameter of transmitter frequency offset (O1) at different temperatures, and the optimization results were 2436.22, 2428.67, and 2416.84 Hz at 296, 298, and 300 K, respectively. Then, NOESYGPPR1D pulse was used to acquire 1H NMR spectra, and the parameters were set according to the experimental conditions in Table 1. The parameters not specified in the table were set as follows: pulse width (SW) was 12.016 ppm, mixing time (D8) was 0.05 s, and the 90° pulse width (P1) was 14.75 μs. Finally, a total of fifteen 1H NMR spectra were acquired, and the data were processed according to the method in Section 3.2.5.

CAPPs Screening Results

The absolute value of the percentage of relative deviation between the measured value of 1H qNMR and the actual weighing value of the DTSs content in the samples was used as the accuracy index. The peak width factor of the quantitative peaks (calculated by Mestrenova software (version 14.0.0)) was used as the resolution index; the ratio of the internal standard peak height to the noise level (calculated by Mestrenova software (version 14.0.0)) was used as the SNR index. The results of the indexes for the 15 PBD experiments are shown in Table 2. The analytical parameters with significant effects on the four indexes were selected as CAPPs by using the weighted multiple linear regression coefficient method. The multiple linear regression model between the analytical parameters (xi) and the evaluation indexes (y) is represented by Equation (1), where a 0 is a constant term and x i and a i are the level and partial regression coefficients of the corresponding analytical parameters, respectively.
y = a 0 + i = 1 7 a i x i
After fitting, the regression coefficients corresponding to the three evaluation indexes are shown in Figure 2, and the determination coefficients (r2) are 0.644, 0.705, and 0.787, respectively. The error line represented the 95% confidence interval for the results of the 15 experiments. If the range of this confidence interval contains 0, the parameter was not significant for the model. As shown in Figure 2, the significant parameter in the SNR model was NS, the significant parameters in the accuracy model were D1 and TD, and the significant parameter in the resolution model was Temp. The principles of significant parameters affecting the three evaluation indexes were further analyzed. A 1H NMR spectrum was usually superimposed by multiple acquisitions, and a larger NS led to higher signal intensity and SNR. D1 affected the time interval between two acquisitions. A too small D1 resulted in incomplete relaxation of the signals, which affected accuracy. TD reflected the inherent resolution of the spectra, and a sufficient amount of data points (TD) could accurately describe the shape of the signal peaks. Temperature affected the velocity of molecular motion, thus affecting the characteristic displacement of the molecule. Strengthening the temperature control means could ensure a good signal resolution, so there is no need to further optimize the temperature setting. Finally, the three remaining significant parameters NS, TD, and D1 were selected as CAPPs for this study.

2.2.2. The Results of CCD-based CAPP Optimization

CCD Experimental Design Results

After determining the CAPPs (NS, TD, and D1), this study used CCD experiments to find the combination of parameters that ensure robust operation. Before that, a parameter closely related to TD and D1 is considered: acquisition time (AQ). AQ is not a directly set parameter but is determined by SW and TD, as shown in Equation (2):
1 A Q = 2 S W T D
The time lag (Tl) between two acquisitions is the sum of AQ and D1. D1 affects the accuracy by influencing Tl, but the change of TD also affects Tl. In order to avoid the interaction between TD and D1 involving the modeling, the three factors in CCD were set as NS (X1), TD (X2), and Tl (X3 = AQ + D1). The parameter levels of NS and TD were determined with reference to the literature and previous experimental experience, and the level of Tl was set according to the longest (longitudinal relaxation time) T1 of the internal standard signal and quantitative signals. The final design yielded 17 optimized experiments with 3 factors and 5 levels. The CCD experiments table and the coding values of each factor are shown in Table 3.
The evaluation indexes of CCD experiments were determined based on ATPs and the results of PBD experiments. The accuracy and SNR indexes are the same as those of the PBD experiments. Since the NS, TD, and D1 were not significant in the resolution model, resolution was no longer used as an evaluation index. Considering the requirement of ATPs on the analysis time, the total acquisition time (Tq) was included in the evaluation index of the optimization experiment.
The sample used for the CCD experiments was the sample with total DTSs concentrations of 5 mmol/L in Section 2.2.1. The 1H NMR spectrum was acquired with NOESYGPPR1D pulse. The acquisition parameters were set according to Table 3, and the parameters not specified in the table were set as follows: SW at 12.016 ppm, Temp at 298 K, O1 at 2428.67 Hz, P1 at 14.75 μs, DS at 4, and RG at 71.8.

MODR Calculation Results Based on the Optimization Model

The results of indexes for the 17 CCD experiments were shown in Table 3. MODDE 13 software (Umetrics Inc., Sartorius AG, Göttingen, Germany.) was used to establish the PLS regression models between the parameters and the evaluation indexes, as shown in Equation (3). Where b 0 were constants, b i , b i i , and b i j were regression coefficients, Y was the evaluation index of the CCD optimization experiments, X i and X j were the coded values of each CCD factor. The determination coefficients, r2 and adjusted r2 (r2adj), are shown in Table 4. The r2 was greater than 0.85 for all three models, which proved that the models fitted well.
Y = b 0 + i = 1 3 b i X i + i = 1 3 b i i X i 2 + i = 1 2 i = 1 3 b i j X i X j
Based on the above models, the Monte Carlo method was applied to calculate the probability-based design space. According to the ATPs, the upper limit of accuracy was set at 3%, the lower limit of SNR was set at 1000, and the upper limit of analysis time was set at 25 min. The calculation parameters were set as follows: the resolution of the simulation set point was 64, the number of simulations calculated for each point was 10,000 times, and the failure probability was 1%. The calculation results are shown in Figure 3. The gradient of color from green to red in the figure indicated that the probability of failure of each index under the corresponding parameter conditions (i.e., the probability of failure) gradually increases from 0.5%. The contour lines and the numbers on the lines in the figure represented the probability of failure of each evaluation index in the corresponding area. In this study, the region with a probability of failure less than or equal to 1% was used as the design space (the green region in Figure 3), and the region that was convenient for practical operation was selected as MODR. MODR was set as follows: TD was 65,536 or 98,304; NS was 32, 40, or 48; and Tl was 20~25 s.

2.3. Procedure Performance Qualification (APPQ)

2.3.1. Determination of Acquisition Parameters

According to ATPs, the analysis time needed to be as short as possible while ensuring the robustness of the results. Therefore, the two parameters related to the analysis time, NS and Tl, should be as small as possible. The final selection of parameter combination is shown as the star point in Figure 3, and the total duration of acquisition was 13 min. The acquisition parameters were determined as follows: 1H NMR spectrum was acquired with a NOESYGPPR1D pulse, with a solvent CD3OD locking field, SW of 12.016 ppm, Temp of 298 K, DS of 4, NS of 32, D1 of 18.73 s, TD of 65,536, O1 of 2428.67 Hz, D8 of 0.05 s, P1 of 14.75 μs, RG of 14.75 μs, and RG of 0.05 s.

2.3.2. APPQ Index Examination

MODR Robustness Examination

To verify the robustness of MODR, six experimental sites were selected in MODR to conduct validation experiments. The CAPPs (NS, TD, D1) settings for the six experiments were (32, 98,304, 19 s), (48, 98,304, 16 s), (40, 98,304, 17 s), (48, 65,536, 20 s), (32, 65,536, 18 s), and (40, 65,536, 19 s). The other parameters are set according to Section 2.3.1. The 1H NMR spectra were collected with the above parameters and analyzed for the indexes specified in Section 2.2.2. The results showed that the accuracy of the six experiments was within 1.78%, the SNR was greater than 1245.87, and the analysis time was within 19.76 min, which demonstrated a robust MODR.

Linearity, Accuracy Investigation, and Calculation of LOD and LOQ

Seven samples with total DTSs concentrations of 9, 5, 2.5, 1.5, 1.0, 0.5, and 0.25 mmol/L were obtained by gradient dilution. The seven samples were analyzed according to the method described in Section 2.3.1. A linear regression was performed with the actual concentration prepared as “y” and the measured concentration as “x”. The regression equation was obtained as y = 0.9968x + 0.0472, with a coefficient (r2) of 0.9999, demonstrating the excellent linearity of the method. The slope of the equation was close to 1.00, indicating good accuracy. Based on the equation, the LOD and LOQ were calculated as described in Section 3.6 and the results were LOD = 0.0794 mmol/L and LOQ = 0.2647 mmol/L. The recovery rates of the six samples with total DTSs content above LOQ were calculated. The results showed that the recoveries of all six samples were in the range of 98.87~102.89%, which met the ATP requirements.

Precision and Sample Stability Investigation

Intra-day and inter-day precision were investigated by repeated experiments, and relative standard deviations (RSDs) of six parallel experiments were 0.67% and 1.14%, respectively. A sample was placed at room temperature and analyzed three times at 0 h and 48 h, respectively, to examine the stability within 48 h. The results show that the RSD of the 6 experiments was 1.81%. All these results met the ATP requirements and proved that the intra-day precision, inter-day precision, and 48 h stability of the samples were good.

Signal Specificity Investigation

According to the structural information provided by the NMR spectrum, the quantitative signals and the internal standard signals were highly specific in the absence of overlapping and interfering impurity signals. Since the samples in method development were prepared with standards, impurity interference could be excluded. Signal overlapping may exist when the method is applied to actual TCM samples. Therefore, quantitative signal specificity needs to be examined at the time of method transfer for each type of sample.

Measurement Uncertainty

The combined uncertainty (CU) and the expanded uncertainty (EU) for the total DTSs 1H qNMR method were calculated from Equations (S2) and (S3). The results were 1.21% and 2.42%, respectively, which met the ATP requirements. CU was the combination of uncertainties from the molecular weight of the internal standard and DTSs, the purity of the internal standard, the weighing measurement process, the volume measurement process, and the integral operation. The contributions of the uncertainties from different sources are shown in Figure 4. The most significant contribution was the integration result error, followed by the volume measurement error.

2.4. CPPV Example: Method Transfer and Application

Traditional method transfer strategies include comparative testing, method co-validation, method re-validation, and transfer exemptions. In the APLC concept, the transfer of analytical methods can be considered part of the method performance validation. In this study, the transfer and application of the 1H qNMR method were performed on three samples of herbal medicines containing DTSs. The three samples included a liquid preparation (Shenmai Injection), a solid preparation (Xuesaitong Injection), and a process intermediate (Gynostemma column chromatography eluting intermediates). For the 1H qNMR method in this study, the conditions of reagents and instruments were unchanged. When the method was transferred to the determination of different kinds of samples, the most significant issue was signal overlapping of impurity signals, which affected the method’s accuracy. Therefore, the transfer of the method required revalidation of the specificity of the quantitative signals and the accuracy of the results, while other methodological indexes did not need to be revalidated. Given that the signals corresponding to both H-26 and H-27 protons can independently serve as quantitative peaks, the method remains valid even if one signal experiences interference while the other remains unaffected. The specificity of the internal standard signals of the above three samples has been verified by 2D NMR experiments. 1H-13C HSQC and 1H-1H COSY results both showed that the internal standard signals were independent. The quantitative signals were the overlapping signals of multiple saponins so their specificity could not be verified simply by 2D experiments, which were confirmed by negative control experiments in this study.

2.4.1. 1H qNMR Analysis of Total Dammarane-Type Ginsenosides in the Shenmai Injection

Shenmai Injection is an herbal preparation with ginsenoside as API, prepared from red ginseng and Ophiopogon japonicus. In total, 600 µL of Shenmai Injection was measured and lyophilized to volatilize the solvent completely. Then, 600 µL of CD3OD solution containing DMT was added to dissolve the remaining solid, and the supernatant was centrifuged at 10,000 rpm for 10 min to obtain the Shenmai Injection NMR sample.
The D101 type macroporous resin was activated by soaking in ethanol for 24 h, mounted on a column (1.5 cm × 12 cm), and washed with deionized water until there was no alcohol smell. In total, 1 mL Shenmai Injection was added to the column and eluted with 25 mL of deionized water at a flow rate of 0.5 mL/min. The eluate was collected, at which time the DTSs were retained in the macroporous resin column. After completely evaporating the solvent of the eluent, 600 µL of CD3OD containing DMT was added to re-dissolve the remaining solid. Then, the mixture was centrifuged at 10,000 rpm for 10 min to get supernatant as the negative control sample A of Shenmai Injection.
In addition, to exclude the influence of Ophiopogon japonicus-related components on the quantitative signals, the Ophiopogon japonicus-related process intermediate “Ophiopogon japonicus aqueous supernatant” from the production of Shenmai Injection was taken and the sample was prepared in the same way as that of Shenmai Injection to obtain the negative control sample B.
The 1H NMR spectra of Shenmai Injection, negative control samples A and B were collected by the method described in Section 2.3.1, and the results are shown in Figure 5. The two negative control samples showed no signal response at 1.68 ppm, which proved that the quantitative peak specificity at 1.68 ppm was good.
After that, the accuracy of the method was investigated by recovery experiments. Three groups of samples were prepared with standard addition amounts of 80%, 100%, and 120%, respectively. The spectra of the above samples were collected to calculate the recovery rates. The average recovery rates were 98.95%, 99.58%, and 100.84% for the low, medium, and high concentration groups, respectively, indicating that 1H qNMR could accurately quantify the total DTSs in the Shenmai Injection.

2.4.2. 1H qNMR Analysis of Total Notoginsenosides in the Xuesaitong Injection

In total, 5.29 mg of Xuesaitong Injection was weighed precisely in a centrifuge tube, and 600 µL of CD3OD containing DMT was added and dissolved to obtain the Xuesaitong NMR sample. Xuesaitong Injection is purified notoginsenoside, the DTSs of which accounted for more than 90% of the total solids. Therefore, the 1H qNMR spectrum Xuesaitong Injection was similar to the samples described in Section 2.2.1, and the interference of impurity signals could be excluded. The signals showed high specificity in the spectrum.
The accuracy of the method was also investigated by recovery experiments. The recovery experiments procedure was similar to that in Section 2.4.1. The average recovery rates of the three concentration groups were 100.62%, 98.14%, and 99.63%, respectively, which demonstrated that 1H qNMR could accurately quantify the total notoginsenosides in the Xuesaitong Injection.

2.4.3. 1H qNMR Analysis of Total Gypenosides in the Gynostemma Process Intermediates

The 70% ethanol eluting intermediate of gynostemma column chromatography was precisely pipetted into a centrifuge tube and concentrated by centrifugation until the water evaporated, after which 600 µL of CD3OD containing DMT was added to dissolve the remaining solid, and the mixture was centrifuged at 10,000× g rpm for 10 min to obtain the supernatant as gynostemma NMR sample. The negative sample was prepared in the same way on water eluting intermediate of gynostemma column chromatography. The NMR spectra of the gynostemma NMR sample and the negative sample were obtained by the method described in Section 2.3.1. The results are shown in Figure 6, and the negative sample showed no signal response at 1.68 ppm and 1.62 ppm, which proved the excellent specificity of the quantitative signals at both locations.
The accuracy was verified by comparison with an independent method of known authenticity. Six batches of gynostemma column chromatography eluting intermediates were analyzed by applying the classical glacial acetic acid-perchloric acid UV chromogenic method and 1H qNMR method respectively to determine total gypenosides. Table 5 presents the outcomes obtained from the utilization of both methods, indicating that the relative deviations of two methods fall within the permissible range of 5%. Paired t-tests were performed on the results for the two groups of measurements. The 95% confidence interval for the deviation of the results between the two methods was (−1.690, 0.409) with a test p-value of 0.178 (>0.05), which proved that there was no significant difference between the two groups of measurements. Both the results of relative deviations and paired t-test demonstrated that the 1H qNMR method was in high agreement with the analytical results of independent methods of known authenticity, i.e., 1H qNMR can accurately quantify the content of total gypenosides in gynostemma process intermediates. Notwithstanding, the 1H NMR technique exhibits a superior level of operational ease, environmental sustainability, and efficiency in contrast to traditional methodologies.

3. Materials and Methods

3.1. Reagents and Materials

Deuterated methanol containing 0.03% v/v tetramethylsilane (TMS) (CD3OD, 99.8% D) was purchased from Cambridge Isotope Laboratories, Inc. (Tewksbury, Massachusetts, USA); dimethyl terephthalate (DMT, Lot No. BCBF6171V) was purchased from Sigma-Aldrich (Darmstadt, Germany); standard substances Ginsenoside Rb1 (Lot No. 200725), Ginsenoside Rg1 (Lot No. 200709), Ginsenoside Re (Lot No. 200603), Panaxoside R1 (Lot No. 190813), Panaxoside Fe (Lot No. 190803), Gynostemma saponin XLIX (Lot No. 210826), Gynostemma saponin XVII (Lot No. 210903), and Gynostemma saponin XLVI (Lot No. 210730) were purchased from Shanghai Winherb Medical Technology Co., Ltd. (Shanghai, China); Shenmai Injection was provided by Chiatai Qingchunbao Pharmaceutical Co., Ltd. (China); Xuesaitong Injection was provided by Heilongjiang ZBD Pharmaceutical Group Co., Ltd. (Hangzhou, China); and gynostemma column chromatography eluting intermediate was provided by Wanbangde Pharmaceutical Group Co., Ltd. (Taizhou, China).

3.2. Primary Analysis Conditions

In this study, the internal standard method was used to quantify total DTSs. Before screening and optimization of CAPPs, some primary analytical conditions need to be initially confirmed, including deuterated reagents, pulse sequences, internal standard substances, and signals for quantification.

3.2.1. Selection of Deuterated Reagents

In NMR-related studies of DTSs, commonly used deuterated reagents include deuterated pyridine (C5D5N) and deuterated methanol (CD3OD). Both solvents can better solubilize DTSs, but the signal resolution in CD3OD is much higher than that in C5D5N. Therefore, CD3OD was chosen as the deuterated reagent in this study.

3.2.2. Selection of Pulse Sequences

The 1H NMR spectrum acquisition with CD3OD as the solvent can be completed directly with a standard hydrogen pulse ZG30. However, the water signal in the spectrum acquired by ZG30 is high, which may affect the sensitivity of signals for quantification, as shown in Figure S1. Therefore, the pre-saturated water suppression pulse NOESYGPPR1D was considered. This study compared the average signal-to-noise ratio (SNR) of the same signals in the spectra acquired under the NOESYGPPR1D and ZG30 pulses. The SNR of the signal acquired under the NOESYGPPR1D pulse was significantly higher than that of the signal acquired under the ZG30 pulse, and the higher percentage was more than 50%, indicating that the use of NOESYGPPR1D pulse has a more remarkable improvement in sensitivity, so this pulse was chosen in this study. The specific experimental results are shown in Section S1.

3.2.3. Selection of Internal Standard

The requirements for the internal standard are high solubility, and good stability without signal overlapping in the sample. DMT was chosen as the internal standard in this study because the spin-lattice relaxation time (T1) of the protons of DMT in CD3OD is relatively short, which can support a shorter analysis time. In the 1H NMR spectrum, DMT has two signals at 3.94 ppm and 8.11 ppm, respectively. The singlet of -O-CH3 at 3.94 ppm tends to overlap with the other signals in the sample, so the singlet of -CH at 8.11 ppm was used for quantification, and the number of equivalent protons for this signal was four.

3.2.4. Selection of Signals for Quantification

The structural skeleton of DTSs is shown in Figure 7A. The H-26 and H-27 in the side chain of the methyl structure are in almost all DTSs. Although some side chain modifications of DTSs, such as cyclization, double bond transfer, etc., have been reported in recent years, the saponins with the above conditions belong to trace components, accounting for less than 1% of the total content of DTSs [2]. The effect of these saponins is negligible for the macroscopic index of “total DTSs” content. The H-26 and H-27 protons signals in the 1H NMR spectrum are at 1.62 ppm and 1.68 ppm, which are the signal I and II labeled in Figure 7B.
In a mixture of multiple DTSs, H-26 and H-27 signals of different DTSs overlap and stack at the positions of signal I and II, respectively, so the areas of these two signals contain the information of total DTSs content and the number of equivalent protons for both sets of peaks is three. In summary, signal I and II can be used individually or in combination as the signals for quantification of total DTSs; and their areas can be used to calculate the exact molar concentration of total DTSs in the sample, as shown in Equation (4):
M T S = N D M T · A T S N T S · A D M T · M D M T
where M T S and M D M T are the molar concentrations of total DTSs and the internal standard (DMT), respectively, A T S and A D M T are the integrated areas of signals of total DTSs and DMT, respectively, and N T S and N D M T are the number of equivalent protons for signals of DTSs and DMT, respectively.

3.2.5. Data Processing

After spectra acquisition, Fourier transform and manual phase correction of the free induction decay (FID) signals were performed in Topspin software (version 3.2.6), and the chemical shift correction was performed with the TMS peak (0.00 ppm) as reference. The window function was selected in exponential mode and the line width (LB) was set as 0.3 Hz. Afterward, the signals of DMT and total DTSs were integrated by “line fitting” in MestReNova software (version 14.0.0), and the molar concentrations of total DTSs were calculated according to Equation (1).

3.3. Determination of ATPs

The prospective requirements of the qNMR method in this study were well defined in the ATPs, which included all elements such as target sample, target APIs, sample preparation, required analytical technique, method requirement, target application, reportable quality attributes, and critical analytical characteristics. Further information regarding the aforementioned details can be found in Section S2.

3.4. Risk Assessment Methodology

The main tools used to perform risk assessment in this study were the Ishikawa diagram and FMECA. Ishikawa diagram was used for risk identification, but it cannot determine the correlation between parameters. Therefore, FMECA was introduced as a risk-ranking method in this study. The risk priority number (RPN) is determined by the severity (S), probability of occurrence (O), and ease of detection (D) of the failure factor. The higher the RPN score, the higher the influence of the factor on the analysis results, which needs to be focused on in the process of method development [12].

3.5. DoE Methods

The purpose of DoE is to establish an experimental program to obtain “maximum information” through statistical methods with “minimum experimentation”. The statistical models are used in DoE to help researchers understand the impact of APPs on method performance and provide reliable estimates of analytical method detection. This study used PBD to further screen the high-risk parameters. The screened CAPPs were optimized by CCD, which is suitable for fitting multivariate non-linear models, ensuring the model’s fit without increasing the number of runs. Based on the optimized model established by CCD, the design space was calculated according to the ATPs requirements and MODR was selected. Design of the experiments, statistical modeling, and design space calculation in this study were performed on MODDE 13 software (Umetrics Inc., Sartorius AG, Göttingen, Germany).

3.6. Indexes of APPQ

Based on the requirements in the ATPs, accuracy, precision, specificity, limit of detection (LOD), limit of quantification (LOQ), linearity, robustness, and measurement uncertainty were determined with reference to ICH Q2 (R1) guideline, the specific validation process is shown in Section S3.

4. Conclusions

In this study, a 1H qNMR method for the determination of total DTSs in herbal medicines was developed based on APLC. In the initial stage of APLC (APD), the Ishikawa diagram and FMECA were applied to identify the potential risk factors and select seven parameters as potential CAPPs. DoE methods were used to screen and optimize the CAPPs, and a robust MODR was obtained through modeling. In the subsequent stage (APPQ), all the methodological indexes performed well; the combined uncertainty and the expanded uncertainty for the 1H qNMR method were 1.21% and 2.42%, respectively. Finally, the developed 1H qNMR method was transferred and successfully applied to the quantitative analysis of total DTSs in three herbal medicines, providing an example of CPPV. It should be noted that the APLC has yet to be completed, and routine monitoring of the method should be carried out during subsequent applications. Timely investigation as well as corrective and preventive measures should be taken if there is any indication that the analytical method is out of control.
In comparison to the conventional UV colorimetric method, the 1H qNMR method developed in this study offers several advantages, including a smaller required sampling amount, simplified sample pretreatment procedures, and avoidance of corrosive reagents. Significantly, the total DTSs content determined via this methodology represents the true value as measured by molar quantity, in contrast to the estimated value derived from the traditional approach. Therefore, the 1H qNMR method provides a higher degree of accuracy. The method development process applied QRM and AQbD concepts, allowing for a comprehensive examination of various factors. It effectively minimized the potential for errors and omissions while optimizing parameter efficiency. In addition, the utilization of models facilitated a quantitative assessment of multiparameter effects and their interactions.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/ph16070947/s1, Section S1: Determination of the acquisition pulse; Section S2: Determination of ATPs; Section S3: Process performance validation; Section S4: Risk identification and assessment results based on Ishikawa diagram and FMECA; Figure S1: Results of SNR paired t-test for signals collected under two pulses.

Author Contributions

Conceptualization, W.L. and H.Q.; methodology, W.L. and J.Y.; validation, W.L., J.Y. and F.Z.; formal analysis, F.Z. and X.X.; investigation, J.Y.; resources, H.Q.; data curation, W.L.; writing—original draft preparation, W.L.; writing—review and editing, J.Y., X.X. and H.Q.; visualization, W.L.; supervision, J.P.; project administration, H.Q.; funding acquisition, H.Q. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National S&T Major Project of China (2018ZX09201011-002).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data underlying this article are available in the article and in its online supplementary material.

Conflicts of Interest

The authors have no conflict of interest to declare.

Abbreviations

DTSs: dammarane-type saponins; 1H NMR: proton nuclear magnetic resonance; 1H qNMR: proton quantitative nuclear magnetic resonance; APIs: active pharmaceutical ingredients; APLC: analytical procedure lifecycle; ATP: analytical target profile; APD: analytical procedure design; APPQ: procedure performance qualification; CPPV: continued procedure performance verification; QRM: quality risk management; DoE: design of experiments; CAPPs: critical analytical procedure parameters; MODR: method operable design region; AQbD: analytical quality-by-design; FMECA: failure mode effects and criticality analysis; PBD: Plackett–Burman design; CCD: central composite design; TMS: tetramethylsilane; DMT: dimethyl terephthalate; SNR: signal-to-noise ratio; RPN: risk priority number; LOD: limit of detection; LOQ: limit of quantification; DS: number of dummy scans; NS: number of scans; TD size of fid; RG: receiver gain; D1: relaxation delay; SW: pulse width; P1: 90° pulse width; Tl: time lag; CU: combined uncertainty; EU: expanded uncertainty.

References

  1. Ruan, J.; Zheng, C.; Qu, L.; Liu, Y.; Han, L.; Yu, H.; Zhang, Y.; Wang, T. Plant resources, 13C-NMR spectral characteristic and pharmacological activities of dammarane-type triterpenoids. Molecules 2016, 21, E1047. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Cao, J.; Zhang, X.; Qu, F.; Guo, Z.; Zhao, Y. Dammarane triterpenoids for pharmaceutical use: A patent review (2005–2014). Expert Opin. Ther. Pat. 2015, 25, 805–817. [Google Scholar] [CrossRef] [PubMed]
  3. Li, W.; Yang, J.; Zhao, F.; Pan, J.; Qu, H. 1H-qNMR quantification of total ginsenosides in Shenmai Injection. Chin. J. Chin. Mater. Med. 2022, 47, 587–592. [Google Scholar] [CrossRef]
  4. Ynag, L.; Chen, J.; Xu, X.; Gu, L.; Zhang, S.; You, Z.; Zheng, G.; Chen, J.; Xin, Y. Simultaneous determination of eleven components in Shenmai injection by UHPLC. Chin. J. Pharm. Anal. 2019, 39, 1660–1665. [Google Scholar] [CrossRef]
  5. Stavrianidi, A.; Stekolshchikova, E.; Porotova, A.; Rodin, I.; Shpigun, O. Combination of HPLC–MS and QAMS as a new analytical approach for determination of saponins in ginseng containing products. J. Pharm. Biomed. Anal. 2017, 132, 87–92. [Google Scholar] [CrossRef]
  6. Abashev, M.; Stekolshchikova, E.; Stavrianidi, A. Quantitative aspects of the hydrolysis of ginseng saponins: Application in HPLC-MS analysis of herbal products. J. Ginseng Res. 2021, 45, 246–253. [Google Scholar] [CrossRef]
  7. Wang, R.; Wang, H.; Fan, X.; Duan, M.; Li, H.; Wang, S. Analysis of chemical constituents in Shenmai Injection by LC-Q-TOF-MS and LC-IT-MS. Chin. J. Chin. Mater. Med. 2020, 45, 105–114. [Google Scholar] [CrossRef]
  8. Verch, T.; Campa, C.; Chéry, C.; Frenkel, R.; Graul, T.; Jaya, N.; Nakhle, B.; Springall, J.; Starkey, J.; Wypych, J.; et al. Analytical Quality by Design, Life Cycle Management, and Method Control. AAPS J. 2022, 24, 34. [Google Scholar] [CrossRef]
  9. Tiwari, R.; Kumar, A.; Solanki, P.; Dhobi, M.; Sundaresan, V.; Kalaiselvan, V.; Raghuvanshi, R.S. Analytical quality-by-design (AQbD) guided development of a robust HPLC method for the quantification of plumbagin from Plumbago species. J. Liq. Chromatogr. Relat. Technol. 2021, 44, 529–537. [Google Scholar] [CrossRef]
  10. Baglow, R.L. The reliability parameter and its importance for life cycle management. Microelectron. Reliab. 1975, 14, 91–104. [Google Scholar] [CrossRef]
  11. Parr, M.K.; Schmidt, A.H. Life cycle management of analytical methods. J. Pharm. Biomed. Anal. 2018, 147, 506–517. [Google Scholar] [CrossRef]
  12. Park, G.; Kim, M.K.; Go, S.H.; Choi, M.; Jang, Y.P. Analytical quality by design (AQbD) approach to the development of analytical procedures for medicinal plants. Plants 2022, 11, 2960. [Google Scholar] [CrossRef]
  13. Kim, M.K.; Park, S.C.; Park, G.; Choi, E.; Ji, Y.; Jang, Y.P. Analytical quality by design methodology for botanical raw material analysis: A case study of flavonoids in Genkwa Flos. Sci. Rep. 2021, 11, 11936. [Google Scholar] [CrossRef]
  14. Zhang, H.; Wang, J.; Chen, Y.; Shen, X.; Jiang, H.; Gong, X.; Yan, J. Establishing the chromatographic fingerprint of traditional Chinese medicine standard decoction based on quality by design approach: A case study of Licorice. J. Sep. Sci. 2019, 42, 1144–1154. [Google Scholar] [CrossRef]
  15. Dai, S.; Xu, B.; Zhang, Y.; Sun, F.; Li, J.; Shi, X.; Qiao, Y. Robust design space development for HPLC analysis of five chemical components in Panax notoginseng saponins. J. Liq. Chromatogr. Relat. Technol. 2016, 39, 504–512. [Google Scholar] [CrossRef]
  16. Gong, X.; Zhang, Y.; Chen, H.; Chen, T.; Pan, J.; Wang, X.; Qu, H. Development of an analytical method by defining a design space: A case study of saponin determination for Panax notoginseng extracts. Anal. Methods 2016, 8, 2282–2289. [Google Scholar] [CrossRef]
  17. Kim, M.K.; Park, G.; Hong, S.P.; Jang, Y.P. Analytical quality by design methodology approach for simultaneous quantitation of paeoniflorin and decursin in herbal medicine by RP-HPLC analysis. Nat. Prod. Sci. 2021, 27, 264–273. [Google Scholar] [CrossRef]
  18. Parab Gaonkar, V.; Hullatti, K. Quality assessment and RP-HPLC method development for estimation of curcuminoids in Curcuma longa: A quality by design approach. J. Liq. Chromatogr. Relat. Technol. 2021, 44, 95–102. [Google Scholar] [CrossRef]
  19. Zhang, S.; Wu, L.; Wang, X.; Gong, X.; Qu, H. Development of an HPLC-MS method for the determination of four terpene trilactones in Ginkgo biloba leaf extract via quality by design. Biomed. Chromatogr. 2021, 35, e5170. [Google Scholar] [CrossRef]
  20. Ancillotti, C.; Orlandini, S.; Ciofi, L.; Pasquini, B.; Caprini, C.; Droandi, C.; Furlanetto, S.; Del Bubba, M. Quality by design compliant strategy for the development of a liquid chromatography–tandem mass spectrometry method for the determination of selected polyphenols in Diospyros kaki. J. Chromatogr. A 2018, 1569, 79–90. [Google Scholar] [CrossRef]
  21. Silva, P.; Silva, C.L.; Perestrelo, R.; Nunes, F.M.; Câmara, J.S. Application of quality-by-design approach in the analytical method development for quantification of sugars in sugarcane honey by reversed-phase liquid chromatography. Food Anal. Methods 2020, 13, 1634–1649. [Google Scholar] [CrossRef]
  22. Van Tricht, E.; Geurink, L.; Backus, H.; Germano, M.; Somsen, G.W.; Sänger-van de Griend, C.E. One single, fast and robust capillary electrophoresis method for the direct quantification of intact adenovirus particles in upstream and downstream processing samples. Talanta 2017, 166, 8–14. [Google Scholar] [CrossRef] [PubMed]
  23. Perovani, I.S.; Serpellone, C.O.; de Oliveira, A.R.M. An appraisal of experimental designs: Application to enantioselective capillary electromigration techniques. Electrophoresis 2021, 42, 1726–1743. [Google Scholar] [CrossRef] [PubMed]
  24. Aboushady, D.; Hanafi, R.S.; Parr, M.K. Quality by design approach for enantioseparation of terbutaline and its sulfate conjugate metabolite for bioanalytical application using supercritical fluid chromatography. J. Chromatogr. A 2022, 1676, 463285. [Google Scholar] [CrossRef] [PubMed]
  25. Muteki, K.; Morgado, J.E.; Reid, G.L.; Wang, J.; Xue, G.; Riley, F.W.; Harwood, J.W.; Fortin, D.T.; Miller, I.J. Quantitative Structure retention relationship models in an analytical quality by design framework: Simultaneously accounting for compound properties, mobile-phase conditions, and stationary-phase properties. Ind. Eng. Chem. Res. 2013, 52, 12269–12284. [Google Scholar] [CrossRef]
  26. Freitas, J.; Silva, P.; Vaz-Pires, P.; Câmara, J.S. A systematic AQbD approach for optimization of the most influential experimental parameters on analysis of fish spoilage-related volatile amines. Foods 2020, 9, E1321. [Google Scholar] [CrossRef]
  27. Robu, S.; Romila, A.; Buzia, O.D.; Spac, A.F.; Diaconu, C.; Tutunaru, D.; Lisa, E.; Nechita, A. Contribution to the optimization of a gas chromatographic method by QbD approach used for analysis of essential oils from Salvia officinalis. Rev. Chim. 2019, 70, 2015–2020. [Google Scholar] [CrossRef]
  28. Plackett, R.L.; Burman, J.P. The design of optimum multifactorial experiments. Biometrika 1946, 33, 305–325. [Google Scholar] [CrossRef]
  29. Gödecke, T.; Napolitano, J.G.; Rodríguez-Brasco, M.F.; Chen, S.-N.; Jaki, B.U.; Lankin, D.C.; Pauli, G.F. Validation of a generic quantitative 1H NMR method for natural products analysis. Phytochem. Anal. 2013, 24, 581–597. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (A) Risk identification for the 1H qNMR method based on Ishikawa diagram; (B) Risk ranking of the 1H qNMR method based on FMECA.
Figure 1. (A) Risk identification for the 1H qNMR method based on Ishikawa diagram; (B) Risk ranking of the 1H qNMR method based on FMECA.
Pharmaceuticals 16 00947 g001
Figure 2. Regression coefficient of multiple linear models of (A) Accuracy, (B) SNR, and (C) resolution.
Figure 2. Regression coefficient of multiple linear models of (A) Accuracy, (B) SNR, and (C) resolution.
Pharmaceuticals 16 00947 g002
Figure 3. Design space of the total DTSs 1H qNMR method.
Figure 3. Design space of the total DTSs 1H qNMR method.
Pharmaceuticals 16 00947 g003
Figure 4. Contribution of different sources of errors to the combined uncertainty.
Figure 4. Contribution of different sources of errors to the combined uncertainty.
Pharmaceuticals 16 00947 g004
Figure 5. 1H NMR spectra of (A) negative control sample A, (B) negative control sample B, and (C) Shenmai Injection sample. The red box indicates the location of a quantitative peaks.
Figure 5. 1H NMR spectra of (A) negative control sample A, (B) negative control sample B, and (C) Shenmai Injection sample. The red box indicates the location of a quantitative peaks.
Pharmaceuticals 16 00947 g005
Figure 6. 1H NMR spectra of (A) the negative control sample and (B) the positive control sample of gypenosides. The red box indicates the location of a quantitative peaks.
Figure 6. 1H NMR spectra of (A) the negative control sample and (B) the positive control sample of gypenosides. The red box indicates the location of a quantitative peaks.
Pharmaceuticals 16 00947 g006
Figure 7. (A) Structure of dammarane-type tetracyclic triterpenes; (B) Signals for quantification of total DTSs.
Figure 7. (A) Structure of dammarane-type tetracyclic triterpenes; (B) Signals for quantification of total DTSs.
Pharmaceuticals 16 00947 g007
Table 1. Experimental design table for PBD screening and coding values of each factor.
Table 1. Experimental design table for PBD screening and coding values of each factor.
NO.x1 (mmol/L)x2 (Times)x3 (Times)x4 (pcs)x5x6 (Seconds)x7 (K)
15 (+1)0 (−1)128 (+1)16,384 (−1)40.3 (−1)2 (−1)300 (+1)
258 (+1)16 (−1)65,536 (+1)40.32296 (−1)
30.5 (−1)812816,384161 (+1)2296
45012865,53640.330 (+1)296
5581665,5361612300
65812816,38416130296
70.5812865,53640.330300
80.5012865,5361612300
90.501665,53616130296
10501616,38416130300
110.581616,38440.330300
120.501616,38440.32296
132 (0)4 (0)64 (0)32,768 (0)90.5 (0)10 (0)298 (0)
14246432,76890.510298
15246432,76890.510298
Table 2. Results of evaluation indexes in PBD experiments.
Table 2. Results of evaluation indexes in PBD experiments.
NO.Accuracy (%)SNRResolution
177.072446.80.71
217.651316.00.62
355.751601.01.35
40.2914174.40.68
527.691343.40.61
60.7092771.81.29
70.4113407.80.96
823.073404.60.75
90.7471660.50.72
109.5721589.50.72
114.133965.61.36
1251.68658.61.15
136.2023288.70.63
146.0553451.90.60
156.2633318.10.60
Table 3. Parameters setting and evaluation index results of the CCD experiments.
Table 3. Parameters setting and evaluation index results of the CCD experiments.
NO.Parameter SettingEvaluation Indexes
X1 (Times)X2 (pcs)X3 (Seconds)Accuracy (%)SNRTq (Minutes)
132 (−1)32,768 (−1)20 (−1)1.7421314.4913.68
23298,304 (+1)20.001.3271270.5913.68
380 (+1)32,76820.001.5671904.4331.93
48098,30420.000.9471908.1231.93
53232,76830.00 (+1)0.7851225.8819.68
63298,30430.000.3961217.2319.68
78032,76830.000.1862006.7445.93
88098,30430.000.0372010.8845.93
95616,384 (−α)25.002.1611621.7427.80
1056131,072 (+α)25.000.0661678.8827.80
1116 (−α)65,53625.001.238843.209.27
1296 (+α)65,53625.000.1302215.8946.35
135665,53616.59 (−α)1.8301689.5719.40
145665,53633.41 (+α)0.0491609.0036.22
1556 (0)65,536 (0)25.00 (0)0.7231577.3427.80
165665,53625.000.6231637.7727.80
175665,53625.000.7181652.3327.80
Table 4. Regression coefficient, r2, and r2adj of PLS models of evaluation index for CCD experiments.
Table 4. Regression coefficient, r2, and r2adj of PLS models of evaluation index for CCD experiments.
Parameter ItemsRegression Coefficient
Accuracy (%)SNRTq (Minutes)
b00.8551633.05027.802
X1−0.228345.69610.240
X2−0.338--
X3−0.485−4.900114.619
X12-−29.7650.002
X32--0.002
X X13-36.9011.701
r20.8530.9881.000
r2adj0.8190.9831.000
Note: Parameters that are not significant in all three models are not listed.
Table 5. Comparison of the ultraviolet colorimetric method and 1H qNMR method.
Table 5. Comparison of the ultraviolet colorimetric method and 1H qNMR method.
No.UV
(mmol/L)
1H qNMR
(mmol/L)
Relative Deviation
(%)
136.4536.53−0.21
239.5641.12−3.78
343.0742.261.92
440.8641.80−2.25
538.9440.80−4.55
636.4536.68−0.61
Mean values39.2239.86−1.61
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, W.; Yang, J.; Zhao, F.; Xie, X.; Pan, J.; Qu, H. Application of the Analytical Procedure Lifecycle Concept to a Quantitative 1H NMR Method for Total Dammarane-Type Saponins. Pharmaceuticals 2023, 16, 947. https://0-doi-org.brum.beds.ac.uk/10.3390/ph16070947

AMA Style

Li W, Yang J, Zhao F, Xie X, Pan J, Qu H. Application of the Analytical Procedure Lifecycle Concept to a Quantitative 1H NMR Method for Total Dammarane-Type Saponins. Pharmaceuticals. 2023; 16(7):947. https://0-doi-org.brum.beds.ac.uk/10.3390/ph16070947

Chicago/Turabian Style

Li, Wenzhu, Jiayu Yang, Fang Zhao, Xinyuan Xie, Jianyang Pan, and Haibin Qu. 2023. "Application of the Analytical Procedure Lifecycle Concept to a Quantitative 1H NMR Method for Total Dammarane-Type Saponins" Pharmaceuticals 16, no. 7: 947. https://0-doi-org.brum.beds.ac.uk/10.3390/ph16070947

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop