Betasciencepress.com | STM Publishing Corporation                                                                                                                                             

Research Article

Development and Validation of a Bioanalytical Method for the Simultaneous Determination of 14 Antiretroviral Drugs using Liquid Chromatography-Tandem Mass Spectrometry

Alper Daskapan1, Kai van Hateren1, Ymkje Stienstra2, Jos Kosterink1,3, Tjip van der Werf2, Daan Touw1,4 Jan-Willem Alffenaar1,*

1University of Groningen, University Medical Center Groningen, Department of Clinical Pharmacy and Pharmacology, Groningen, The Netherlands.
2University of Groningen, University Medical Center Groningen, Department of Internal Medicine – Infectious Diseases, Groningen, The Netherlands.
3University of Groningen, Groningen Research Institute of Pharmacy, PharmacoTherapy,- Epidemiology and-Economy, Groningen, The Netherlands.
4University of Groningen, Groningen Research Institute of Pharmacy, Unit Pharmacokinetics, Toxicology and Targeting, Groningen, The Netherlands.
Vol.4, No. 2, Pages 37-50, doi: 10.17145/jab.18.007. (ISSN 2405-710X).   Download PDF  

Correspondence

University Medical Center Groningen, Department of Clinical Pharmacy and Pharmacology. P.O. Box 30.001, 9700 RB Groningen, the Netherlands. Phone: +31 503614071; Fax: +31 503614087. Email: This email address is being protected from spambots. You need JavaScript enabled to view it.

Funding/Manuscript writing assistance

The authors have no financial support or funding to report and they also declare that no writing assistance was utilized in the production of this article.

Competing interest

The authors have declared that no competing interest exist.

Article history

Received: 29 November 2017, Revised 17 January 2018, Accepted 25 January 2018.

 

Abstract

OBJECTIVES: The aim was to develop and validate a LC-MS/MS assay to determine antiretrovirals in human plasma for routine therapeutic drug monitoring.
METHODS: The selectivity, sensitivity, linearity, accuracy, precision, recovery, matrix effect, stability and dilution integrity and carry-over were validated according to EMA and FDA standards.
RESULTS: For accuracy and precision, the highest overall bias was 11.3% at LLOQ of both lopinavir and saquinavir. The highest overall CV was 15.6% at the LLOQ of darunavir. Storage stability at 4°C, 20–25°C and 10°C had a maximum CV of 13.2% at low QC level (0.2 mg/L) for saquinavir. Freeze-thaw stability had a maximum overall bias of 7.4% at low QC level (0.8 mg/L) for tipranavir. Selectivity and specificity showed no interfering peaks of more than 20% of the LLOQ.
CONCLUSIONS: The bioanalytical method is suitable for both TDM in standard care and clinical studies.

Keywords: lc-ms/ms, antiretrovirals, hiv, TDM.

Introduction
The treatment of Human Immunodeficiency Virus (HIV) has improved significantly over the last three decades. This improvement was entirely due to the broad introduction of antiretroviral (ARV) drugs in a combined treatment regimen (cART) [1]. With the advent of cART the mortality and morbidity associated with HIV-1 infection and acquired immunodeficiency syndrome (AIDS) sharply decreased [2]. ARV drugs are designed to intervene in the HIV replicative cycle, which eventually results in an HIV ribonucleic acid (RNA) load decrease and subsequently, the recovery of the host immune system. Currently, 25 ARV drugs have been approved for the treatment of HIV by the Food and Drug Administration (FDA) and the European Medicines Agency (EMA) [3]. These drugs are divided in six classes, each class representing a target in the HIV life cycle: Nucleoside Reverse Transcriptase Inhibitors (NRTI), Nonnucleoside Reverse Transcriptase Inhibitors (NNRTI), Protease Inhibitors (PI), Fusion Inhibitors (FI), Entry Inhibitors and integrase strand transfer inhibitors (INI) [3].
In order to optimize treatment outcome and to prevent drug resistance, therapeutic drug monitoring (TDM) is recommended in cases of: drug-drug interaction; renal or hepatic morbidity; pregnancy; administration of drug doses not commonly used; virologic failure; suspicion of nonadherence and adverse events [4]. In order to be eligible for TDM a correlation should exist between drug concentrations and effect or adverse effects; or, a drug should have large inter-individual pharmacokinetic variability; or a narrow therapeutic index [5]. In daily practice at this point NNRTIs, PIs and INIs qualify for TDM [5-7].
Previously, several assays have been published describing the simultaneous determination of ARV drugs for TDM [8-11]. The analytical methods described often used UV detection or included an extensive sample preparation such as liquid-liquid extraction. These methods are time-consuming and incur high costs. Lliquid chromatography (LC) tandem mass spectrometry (MS) is considered to be the most important analytical technique which is superior to UV detection due to its specificity [12]. Further, using a stable isotope internal standard (IS) for LC tandem MS is highly recommended as it corrects for injection and ionization variability [13]. A stable isotope IS will compensate for these deviations and ensures a robust, high-throughput assay [13].
The aim of this study was to develop and validate a bioanalytical assay to determine both older and new PIs and NNRTIs, and two INIs with a simple sample preparation and using stable isotope IS followed by LC-MS/MS analysis in human plasma as this approach is more cost-effective, less time-consuming than most methods using UV detection or extensive sample preparation.

Materials and Methods
Analysis
The following ARVs were included in the method: amprenavir, atazanavir, darunavir, dolutegravir, etravirine, indinavir, lopinavir, nelfinavir, nelfinavir-M8, nevirapine, raltegravir, ritonavir, saquinavir and tipranavir. The corresponding stable isotope IS were used: [2H4]-amprenavir, [2H6]-atazanavir, [13C6]-darunavir, [13C,2H5]-dolutegravir, [13C6]-etravirine, [13C4,15N]-indinavir, [2H8]-lopinavir, [2H3]-nelfinavir, [2H4]-nevirapine, [2H6]-raltegravir, [13C, 2H3]-ritonavir, [2H9]-saquinavir and [2H4]-tipranavir. The ARVs and IS used (Figure 1), were purchased from Alsachim (Illkirch, Graffenstaden, France), except for nelfinavir-M8 which was purchased from Pfizer Inc. (New York, Unites States of America). Acetonitrile for LC-MS/MS was purchased from BioSolve (Valkenswaard, The Netherlands). The chemicals used, including methanol and trifluoroacetic acid, were of high-pressure liquid chromatography (HPLC) or analytical grade and were purchased from VWR (Amsterdam, The Netherlands). Purified water was obtained from a Milli-Q water purifying system (Millipore Corporation, Billerica, MA, USA). Both pooled human serum samples and pooled human plasma samples with ethylene diamine tetraacetic acid (EDTA) as anticoagulant (Becton Dickinson Vacutainer®, K2E 7.2 mg) were made available in accordance with the standard operating procedures of our medical center.
One milligram of each stable isotope IS was added to methanol to a concentration of 1 g/L and stored at -80°C. The precipitation reagent consisted of a mixture of methanol Lichrosolv, acetonitrile (4:21, v/v) and IS. Twenty-five microliters per IS was standard spiked directly in 250 mL precipitation reagent. To 10 µL of each sample, a volume of 750 µL of precipitation reagents combined with the IS was added to a vial. The samples were vortexed for 1 min. The vials were centrifuged for 5 min at 9500g. Ten microliters of the upper layer was injected into the LC-MS/MS.
The analysis was performed on a triple-quadrupole LC-MS/MS (Thermo Scientific, San Jose, CA, USA) with an MS Pump Plus (Finnigan Surveyor) and autosampler Plus (Finnigan Surveyor). The mass spectrometer was a triple-stage quadrupole Quantum Access Max mass spectrometer. The autosampler temperature was set at 10°C. Liquid chromatographic separation was performed on a HyPURITY C18 analytical column (50 by 2.1 mm, 3 µm particle size; Thermo Scientific, Interscience, Breda, The Netherlands), and the temperature was set at 20°C. The mobile phase had a flow rate of 0.5 mL/min and consisted of purified water, acetonitrile, and an aqueous buffer (containing ammonium acetate [5 g/L], acetic acid [35 mL/L of water], and trifluoroacetic acid [2 mL/L of water]. The method had a run time of 2.9 min and chromatographic separation was performed by means of an elution gradient shown in Table 1. Peak area ratios of the ARVs and their IS were used to calculate concentrations.
The MS was operating in positive electrospray ionization mode and selected reaction monitoring (SRM) mode with a positive spray voltage of 5000V, a capillary temperature of 350°C, a vaporizer temperature of 150°C and a sheath gas pressure and auxiliary pressure of 60 and 10 arbitrary units, respectively. High-purity argon was used for both the source and collision gas flows and high-purity nitrogen for the sheath gas. The mass transitions used in ARV drug quantification and the detail mass spectrometry conditions are demonstrated in Table 2. A scan width of 0.5 m/z was used for all components. Peak area integration for all components was calculated by Xcalibur software version 2.0.7 (Thermo Fisher, San Jose, CA, USA).

Method validation
The analytical method was validated in accordance with the guidance for Industry of the FDA and the EMA guidelines [14,5]. For the validation of the assay the selectivity, sensitivity, linearity, accuracy, precision, recovery, matrix effect, stability and dilution integrity were determined.
For all components separate stock solutions were made in a concentration of 10 g/L, with the exception of tipranavir; for this drug a stock solution of 20 g/L was made. All stock solutions were diluted to 500 mg/L, except for the stock solution of tipranavir that was diluted to 2 g/L. Dimethylsulfoxide (DMSO) was used as solvent for dilution. Subsequently 1, 2, 5, 10, 25, 50, 100, 150 and 200 µL of each diluted stock solution was spiked to 5 mL EDTA plasma, resulting in a final composition of 4% organic solvent in the calibration standards. The analysis for linearity was conducted in 5 replicates per concentration. The selection of the assay range was based on the utilization for TDM in daily practice. The calibration standards, blanks, and quality control (QC) samples were fully thawed at room temperature. The criteria used for the selection of the QC concentration levels were based on FDA and EMA guidelines [14,15], LOW was 2 or 3 times the LLOQ; MED was at 40% and HIGH was at 80%. QC samples and calibration standards were stored at -20°C. QC samples with 4 different concentrations of each ARV were used. The concentrations and number of samples used for the calibration curves in combination with the concentrations used for the QC samples per component are shown in Table 3.
For selectivity, 6 pooled human plasma samples were examined for interference and their responses were compared with those of the LLOQ samples. Over 3 days, each day a single calibration curve in plasma was analysed and accuracy was measured by evaluation of five determinations per QC sample on three consecutive days. Precision was divided into within-run and between-run values using the same method as used for the accuracy. The coefficient of variation (CV) for the LLOQ was maximized at 20% deviation and the CV for the other QC samples should not exceed 15%. The recovery was determined on three levels (LOW, MED, and HIGH) and was done in five replicates. The stabilities of the ARVs were tested for storage stability and freeze-thaw stability. Storage stability of the ARVs was examined by storing QC samples at room temperature (20°C to 25°C) in a refrigerator at 4°C and after sample preparation in the autosampler at 10°C, all for 7 days. Stability was also tested using five freeze-thaw cycles at -20°C. All stability tests were done using two different QC levels (LOW and HIGH) in five determinations per concentration. The FDA does not describe a maximum CV requirement for stability, therefore a maximum CV of 15% was used for the current method in accordance with the EMA guidelines [14]. In order to determine potential differences between the analysis of ARVs in human plasma and in human serum, a matrix comparison was performed. Since protein precipitation was the single required sample preparation step in this method, relative recovery was measured by comparing the ratios of integrated peak area of the ARVs and the corresponding IS of the QC samples processed with the average peak area of the recovery samples. Recovery samples (low, med, and high) were post-extraction blank samples spiked at the same concentrations as the QC samples. To determine the dilution integrity, on three consecutive days, for each ARV drug a sample was diluted 10 times and then prepared in five replicates. To determine the carry-over a blank sample was analysed after analysing the highest calibration standard per component.

Clinical application
After the method was validated and approved it was used for TDM in standard care. For use in daily practice blood was collected in a EDTA blood collection tube. Whole blood was centrifuged with 9500 x g for 5 min and a minimum of 0.500 mL plasma was transferred to a vial and stored until analysis in a – 80°C refrigerator. For analysis 10 µL of plasma was transferred to a vial and 750 µL of precipitation reagent was added. The sample was vortexed for 1 min and subsequently centrifuged for 5 min at 9500 x g. Ten microliters of the upper layer was injected into the LC-MS/MS. One calibration curve was used containing all 14 components. The performance of the assay is routinely evaluated by participation in the international quality control program of ARV drugs of the Association for Quality Assessment in Therapeutic Drug Monitoring and Clinical Toxicology (KKGT) [7,16].

Statistics
Results were analysed using one-way analysis of variance (ANOVA) in validated Excel sheets (Microsoft, Redmond, WA).

Results
Method development
For this bioanalytical method several points of consideration were identified regarding mobile- and stationary phase, sample pre-treatment, MS parameters and usability in daily clinical practice. Earlier studies showed adequate liquid-liquid extraction [8,10,17], however simple sample precipitation was preferred for the current method since it is less time-consuming and less expensive [18]. The precipitation reagent has been used for many other bioanalytical methods in the current composition, among others for the determination of azole antifungals and ertapenem [19,20], and therefore was found to be suitable for this bioanalytical method. In order to obtain adequate retention of the highly polar ARVs an analytical column with high stability under aqueous mobile phase conditions was required. The HyPURITY C18 analytical column has a good polar retention capacity and showed adequate peak shape, retention time and has been extensively used in our daily routine making it a sufficient column for bioanalytical determinations in daily practice in combination with the mobile phase used in a stepwise gradient as demonstrated in Table 1. Since the particle diameter of the HyPURITY C18 analytical column was 3 µm the flow rate of the mobile phase was set at 0.5 mL/min in order to ensure optimized resolution and retention time (Figure 2). When the flow rate of the mobile phase was set at 0.5 mL/min it was important to set the detector capillary temperature at 350°C. The adjustments of these parameters was necessary to maximize the ionization and subsequently the sensitivity of the detector.

Method validation results
The validation results for linearity, accuracy and precision are demonstrated in  Table 4. For accuracy and precision, the highest overall bias found during the validation was 11.3% at LLOQ of both lopinavir and saquinavir. The highest overall CV was 15.6% at the LLOQ of darunavir. The results for the storage- and freeze-thaw stability are displayed in Table 5. Storage stability in a refrigerator (4°C), at room temperature (20 – 25°C) and in the autosampler (10°C) met the requirements with a maximum CV of 13.2% at low QC level (0.2 mg/L) for saquinavir. Furthermore, freeze-thaw stability was within FDA and EMA margins, with a maximum overall bias of 7.4% at low QC level (0.8 mg/L) for tipranavir. Selectivity and specificity showed no interfering peaks of more than 20% of the LLOQ. For the matrix comparison the integrated peak height ratio of a component and its corresponding IS in plasma was compared to that in serum and no statistically significant difference between the two matrices was found for each component. Recovery samples (LOW, MED, HIGH) were assayed and recovery of all components was higher than 85%, with the most deviating recovery of 93.4% at low QC level for etravirine. The dilution integrity was proven with a maximum overall bias of 5% for raltegravir. The carry-over testing showed that there was no carry-over. All results of the method validation were within the EMA and FDA guidelines.

Clinical application results
This method has successfully been used in clinical practice for TDM in cases of drug-drug interactions, renal or hepatic morbidities, pregnancy, virologic failure, suspicion of nonadherence and adverse events and it was also used in a clinical study on darunavir pharmacokinetics [21]. In the period from January 2016 until May 2017 403 plasma samples were analysed using this bioanalytical method (Table 6). Performance of the assay was within the accepted margins (accuracy and precision <20%) of the international quality control program of ARV drugs of the KKGT in our laboratory.

Figure 1

 Figure 1n

Figure 2

 Figure2 jab18007

Figure 3

 

Table 1

 Table 1

Table 2

 

 Table2 jab18007

Table 3

 Table3 jab18007

 

Table 4

 Table 4 jab18007

 

Table 4

Table 5 jab18007

 

Table 5

Table 5n jab18007

 

Table 6

 Table 6 jab18007

Discussion and conclusion
The bioanalytical method described here is rapid and simple and provides an efficient tool for TDM of ARV drugs. This bioanalytical method contains both new ARVs as older ARVs such as tipranavir, saquinavir , indinavir and nelfinavir, making it also suitable for resource limited settings were these drugs are currently commonly used. We were struck by the good bias and precision data illustrating the robustness of our analysis.
Previously, other comparable assays have been published describing the simultaneous determination of PIs, NNRTs and INIs [8,10,17,22]. Three of these studies required a time consuming liquid-liquid extraction and did not include stable isotope IS [8,10,17], what could result in a poorer compensation for inefficiencies in extraction and sample preparation steps, as well as for any matrix effect. One of these bioanalytical methods contained 17 ARVs, but required liquid-liquid extraction and did not include stable isotope IS [8]. All four of these studies had a relatively long run time varying from 10 – 30 min. For TDM purposes a relatively small sample volume and simple sample preparation is desirable. Therefore, the method described by Kromdijk et al. appears to be the best applicable in daily practice compared to the other three methods published [22]. The current described bioanalytical method has a shorter run time and includes more ARVs which can be determined simultaneously than the study of Kromdijk et al. Due to the simple sample preparation, which only included protein precipitation using a precipitation reagent, no time-consuming and expensive liquid-liquid extractions and solid-phase extractions were required. The short run time of 2.9 min, the use of small volumes of plasma and the possibility of the simultaneous determination of 14 ARV drugs makes the method we described highly suitable for routine service. The currently developed method and the applicability for routine service is supported by a recent study by Baldelli et al. [18] where a bioanalytical method for elvitegravir was developed similar to the current one. Baldelli et al. utilized the quality by design approach for the method development which entails a systematic method for quality standards by looking at the entire development system and product life cycle [18].
In the current bioanalytical method overlapping retention times of some components were demonstrated. For the development of this method we achieved a high resolution on the one hand and a fast run-time on the other hand. However, due to its desired applicability for routine service to clinicians, and the high sensitivity and discriminatory power of LC-MS/MS, a faster run-time prevailed over higher resolution.
The bioanalytical method was simple, specific, robust, reproducible and demonstrated a high sensitivity for all components and better cost-effectiveness than the commonly used methods. Our assay platform is suitable for TDM in standard care and in clinical studies for old and new ARVs.

References
1. Hull MW, Montaner J. Antiretroviral therapy: a key component of a comprehensive HIV prevention strategy. Curr HIV/AIDS Rep 8(2), 85-93 (2011). [CrossRef]
2. Arts EJ, Hazuda DJ. HIV-1 antiretroviral drug therapy. Cold Spring Harb Perspect Med 2(4), a007161 (2012). [CrossRef]
3. U.S. Department of Health and Human Services. U.S. Food and Drug Administration: Antiretroviral drugs used in the treatment of HIV infection. 2017(9/19) (2016). https://www.fda.gov/forpatients/illness/hivaids/treatment/ucm118915.htm
4. Panel on Antiretroviral Guidelines for Adults and Adolescents. Guidelines for the use of antiretroviral agents in HIV-1 infected adults and adolscents. Departmet of Health and Human Services. 2017(9/19) (2017). https://aidsinfo.nih.gov/guidelines.
5. Punyawudho B, Singkham N, Thammajaruk N et al. Therapeutic drug monitoring of antiretroviral drugs in HIV-infected patients. Expert Rev Clin Pharmacol 9(12), 1583-1595 (2016). [CrossRef]
6. Acosta EP, Gerber JG, Adult Pharmacology Committee of the AIDS Clinical Trials Group. Position paper on therapeutic drug monitoring of antiretroviral agents. AIDS Res Hum Retroviruses 18(12), 825-834 (2002). [CrossRef]
7. Burger D, Krens S, Robijns K, Aarnoutse R, Bruggemann R, Touw D. Poor performance of laboratories assaying newly developed antiretroviral agents: results for darunavir, etravirine, and raltegravir from the international quality control program for therapeutic drug monitoring of antiretroviral drugs in human plasma/serum. Ther Drug Monit 36(6), 824-827 (2014). [CrossRef]
8. Jung BH, Rezk NL, Bridges AS, Corbett AH, Kashuba AD. Simultaneous determination of 17 antiretroviral drugs in human plasma for quantitative analysis with liquid chromatography-tandem mass spectrometry. Biomed Chromatogr 21(10), 1095-1104 (2007). [CrossRef]
9. Proust V, Toth K, Hulin A, Taburet AM, Gimenez F, Singlas E. Simultaneous high-performance liquid chromatographic determination of the antiretroviral agents amprenavir, nelfinavir, ritonavir, saquinavir, delavirdine and efavirenz in human plasma. J Chromatogr B Biomed Sci Appl 742(2), 453-458 (2000). [CrossRef]
10. Turner ML, Reed-Walker K, King JR, Acosta EP. Simultaneous determination of nine antiretroviral compounds in human plasma using liquid chromatography. J Chromatogr B Analyt Technol Biomed Life Sci 784(2), 331-341 (2003). [CrossRef]
11. Villani P, Feroggio M, Gianelli L et al. Antiretrovirals: simultaneous determination of five protease inhibitors and three nonnucleoside transcriptase inhibitors in human plasma by a rapid high-performance liquid chromatography--mass spectrometry assay. Ther Drug Monit 23(4), 380-388 (2001). [CrossRef]
12. Veringa A, Sturkenboom MG, Dekkers BG et al. LC-MS/MS for therapeutic drug monitoring of anti-infective drugs. TrAC Trends in Analytical Chemistry. 84, 34-40 (2016). [CrossRef]
13. Lowes S, Jersey J, Shoup R et al. Recommendations on: internal standard criteria, stability, incurred sample reanalysis and recent 483s by the Global CRO Council for Bioanalysis. Bioanalysis 3(12), 1323-1332 (2011). [CrossRef]
14. European Medicine Agency (EMA), Committe for Medicinal Product for Human Use. Guideline on bioanalytical method validation. 2017 (2012). http://www.ema.europa.eu/docs/en_GB/document_library/Scientific_guideline/2011/08/WC500109686.pdf
15. U.S. Department of Health and Human Services Food and Drug Administration. Guidance for industry, bioanalytical method validation. Center for Drug Evaluation and Research, Rockville, MD. 2017 (2001). https://www.fda.gov/downloads/Drugs/Guidance/ucm070107.pdf
16. Burger D, Teulen M, Eerland J, Harteveld A, Aarnoutse R, Touw D. The International Interlaboratory Quality Control Program for Measurement of Antiretroviral Drugs in Plasma: a global proficiency testing program. Ther Drug Monit 33(2), 239-243 (2011). [CrossRef]
17. Fayet A, Beguin A, Zanolari B et al. A LC-tandem MS assay for the simultaneous measurement of new antiretroviral agents: Raltegravir, maraviroc, darunavir, and etravirine. J Chromatogr B Analyt Technol Biomed Life Sci 877(11-12), 1057-1069 (2009). [CrossRef]
18. Baldelli S, Marrubini G, Cattaneo D, Clementi E, Cerea M. Application of Quality by Design Approach to Bioanalysis: Development of a Method for Elvitegravir Quantification in Human Plasma. Ther Drug Monit 39(5), 531-542 (2017). [CrossRef]
19. Alffenaar JW, Wessels AM, van Hateren K, Greijdanus B, Kosterink JG, Uges DR. Method for therapeutic drug monitoring of azole antifungal drugs in human serum using LC/MS/MS. J Chromatogr B Analyt Technol Biomed Life Sci 878(1), 39-44 (2010). [CrossRef]
20. van Rijn SP, Wessels AM, Greijdanus B, Touw DJ, Alffenaar JW. Quantification and validation of ertapenem using a liquid chromatography-tandem mass spectrometry method. Antimicrob Agents Chemother 58(6), 3481-3484 (2014). [CrossRef]
21. Daskapan A, Dijkema D, de Weerd DA et al. Food intake and darunavir plasma concentrations in people living with HIV in an outpatient setting. Br J Clin Pharmacol 83(10), 2325-2329 (2017). [CrossRef]
22. Kromdijk W, Pereira SA, Rosing H, Mulder JW, Beijnen JH, Huitema AD. Development and validation of an assay for the simultaneous determination of zidovudine, abacavir, emtricitabine, lamivudine, tenofovir and ribavirin in human plasma using liquid chromatography-tandem mass spectrometry. J Chromatogr B Analyt Technol Biomed Life Sci 919-920, 43-51 (2013). [CrossRef]

 

All site content, except where otherwise noted, is licensed under a Creative Commons Attribution 4.0 License.

  CC by

Quantification of isoniazid, pyrazinamide and ethambutol in serum using liquid chromatography-tandem mass spectrometry

Click on the badge to see more information about this article.

 

Article Tools

 Open Access Article

 Peer-reviewed Article

cloud download 2 Free download (PDF)

 Creative Commons Attribution 4.0 License 

Go to top

Betasciencepress Publishing © 2014-2018