Skip to main content

ORIGINAL RESEARCH article

Front. Med., 22 July 2024
Sec. Translational Medicine
This article is part of the Research Topic Recent Advances in Pharmaceutical Analysis: Applications and New Challenges for the Quality of Medicines View all 4 articles

Advanced modeling of pharmaceutical solubility in solvents using artificial intelligence techniques: assessment of drug candidate for nanonization processing

Turki Al Hagbani
Turki Al Hagbani1*Sameer AlshehriSameer Alshehri2Sami BawazeerSami Bawazeer3
  • 1Department of Pharmaceutics, College of Pharmacy, University of Hail, Hail, Saudi Arabia
  • 2Department of Pharmaceutics and Industrial Pharmacy, College of Pharmacy, Taif University, Taif, Saudi Arabia
  • 3Department of Pharmaceutical Sciences, Faculty of Pharmacy, Umm Al-Qura University, Makkah, Saudi Arabia

This research is an analysis of multiple regression models developed for predicting ketoprofen solubility in supercritical carbon dioxide under different levels of T(K) and P(bar) as input features. Solubility of the drug was correlated to pressure and temperature as major operational variables. Selected models for this study are Piecewise Polynomial Regression (PPR), Kernel Ridge Regression (KRR), and Tweedie Regression (TDR). In order to improve the performance of the models, hyperparameter tuning is executed utilizing the Water Cycle Algorithm (WCA). Among, the PPR model obtained the best performance, with an R2 score of 0.97111, alongside an MSE of 1.6867E-09 and an MAE of 3.01040E-05. Following closely, the KRR model demonstrated a good performance with an R2 score of 0.95044, an MSE of 2.5499E-09, and an MAE of 3.49707E-05. In contrast, the TDR model produces a lower R2 score of 0.84413 together with an MSE of 7.4249E-09 and an MAE of 5.69159E-05.

1 Introduction

Simulation and modeling of pharmaceutical processes are great tools for development of pharmaceutical manufacturing and would help for shifting from batch toward continuous manufacture mode. Process analytical technology (PAT) and process modeling are important elements of process understanding for design of continuous manufacturing in pharmaceutical area (13) where they can be exploited to enhance the efficiency of process and products quality. Process modeling can be performed by finding the relationships between the process parameters and critical quality attributes of finished products. Once the relationship has been established, one can implement Quality-by-Design (QbD) for improvement of process and products (4, 5). Thus, development of robust and rigorous models for pharmaceutical processing is a major challenge which should be addressed.

There are different processing routes for manufacture of solid-dosage oral products which can be optimized via process modeling and theoretical computations. For instance, the method of wet granulation can be simulated via population balance model (PBM) to predict the granule size distribution. As the granules’ properties can affect the tablet characteristics, building relationship between granule size and tablet properties can be done via process modeling which would help for process understanding (68).

One of the main problems in pharmaceutical area is the poor solubility of APIs (Active Pharmaceutical Ingredients) in aqueous media which makes patients to take more dosage of drugs to obtain the therapeutics effectiveness. Taking more dosages would consequently result in side effects for patients. Therefore, some techniques have been developed to enhance the solubility of APIs in aqueous media such as production of nanomedicines (911). The processing of drugs with supercritical carbon dioxide has been reported as one of the methods for creation of nanomedicines. This method has attracted much attention and can be considered as a green method for preparation of nanosized APIs particles. Some computational models have been developed for description of this process, where the majority of studies focus on correlation of drug solubility dataset (1215). Machine learning models are among the most commonly used methods for correlation of solubility data which can be used for a given dataset (16, 17).

Modern data analysis now heavily depends on machine learning (ML), which provides strong instruments and methods for deriving important insights from large datasets. ML includes a range of algorithms that let computers learn from data and, in the absence of explicit programming, make predictions or decisions. These algorithms can identify patterns, relationships, and trends within data, making ML particularly useful for tasks such as regression, classification, clustering, and anomaly detection (18, 19). The ML models have been used recently for correlation of pharmaceutical solubility in supercritical solvents such as CO2. Abouzied et al. (20) investigated the drug solubility in supercritical CO2 using multi-layer perceptron, k-nearest neighbors, and GPR methods. A great fitting accuracy was obtained with R2 more than 0.99 which confirmed the validity of ML models in estimating drug solubility. Support vector machine (SVM) has been one of the major methods used for evaluation of drug solubility in supercritical solvent which is useful in this area with great fitting accuracy (21).

In this paper, we focus on the application of ML to model the solubility of ketoprofen in supercritical carbon dioxide, a critical factor in pharmaceutical manufacturing processes. Accurate modeling of solubility is essential for optimizing production efficiency and ensuring the quality of pharmaceutical products. To achieve this, we employed three regression techniques:

1. Piecewise Polynomial Regression (PPR): PPR partitions the input data range into segments and applies polynomial functions to each segment, enabling a flexible and localized estimation of the regression function.

2. Kernel Ridge Regression (KRR): KRR combines ridge regression with kernel functions, enabling it to model non-linear relationships in the data by projecting it into higher-dimensional spaces.

3. Tweedie Regression (TDR): TDR is a generalized linear model capable of accommodating various distribution types, rendering it well-suited for modeling continuous, non-negative data with variance scaling with a power of the mean.

To optimize the performance of these regression models, we utilized the Water Cycle Algorithm (WCA), a nature-inspired optimization method that simulates the water cycle process to find optimal solutions. WCA has proven effective in navigating complex search spaces and identifying optimal parameter settings for various ML models.

The selected models are highly suitable for small datasets because of their adaptability and resilience in capturing complex relationships. PPR enables localized polynomial fits that can adapt to specific segments of the data, KRR effectively handles non-linearities even with limited data points using kernel functions, and TDR is designed for modeling continuous non-negative data with scalable variance. The careful application of these models, along with the Water Cycle Algorithm (WCA) for optimal parameter tuning, ensures that they are capable of delivering accurate and reliable solubility predictions despite the small dataset size.

Multiple contributions are made by this paper. First, we evaluate the three regression models for ketoprofen solubility prediction, highlighting their pros and cons. Second, we demonstrate how the Water Cycle Algorithm improves model performance by tuning hyperparameters. Finally, we present detailed dataset visualizations and statistical analyses to reveal T(K), P(bar), and solubility relationships. This study uses advanced ML techniques and robust optimization strategies to predict supercritical fluid solubility and advance pharmaceutical process optimization. Indeed, Piecewise Polynomial Regression (PPR), Kernel Ridge Regression (KRR), and Tweedie Regression (TDR) were used for the first time with the Water Cycle Algorithm (WCA) optimizer to improve the prediction accuracy of models for solubility of ketoprofen in supercritical CO2. The models are then used for evaluation of effect of temperature and pressure on the solubility variations.

2 Data set description

The dataset includes solubility measurements at temperatures spanning from 308.15 K to 338.15 K and pressures ranging from 160 bar to 400 bar, sourced from (22). The entire data points of the dataset are shown in Table 1. So, T and P are taken as inputs of ML models, and the drug solubility has been considered as the single output for all models in this study.

Table 1
www.frontiersin.org

Table 1. Complete dataset of ketoprofen solubility in supercritical carbon dioxide (22).

Figure 1 depicts distribution plots for temperature (T), pressure (P), and ketoprofen solubility in supercritical carbon dioxide. These plots use kernel density estimates (KDE) to represent the probability density functions of the variables. The solubility distribution is right-skewed, indicating that higher solubility values occur less frequently in the dataset.

Figure 1
www.frontiersin.org

Figure 1. Distribution plots of temperature (T), pressure (P), and solubility of ketoprofen in supercritical CO2.

Furthermore, the violin plots in Figure 2 depict the temperature (T), pressure (P), and solubility of ketoprofen in supercritical carbon dioxide. Each plot combines a boxplot with a kernel density estimate (KDE). The KDE component shows the probability density of the data at various values, whereas the boxplot component within the violin plot displays the median, interquartile range, and potential outliers. The temperature and pressure data have relatively uniform distributions, whereas the solubility data has a more skewed distribution, indicating that solubility values vary across experimental conditions.

Figure 2
www.frontiersin.org

Figure 2. Violin plots illustrating the distributions of temperature (T), pressure (P), and the solubility of ketoprofen in supercritical CO2.

3 Methodology

3.1 Water cycle algorithm

Inspired by the natural water cycle, the Water Cycle Algorithm (WCA) is a population-based optimization algorithm. The method relies on the water cycle, which includes evaporation, cloud formation, and precipitation. The WCA replicates this cycle to find the best solutions. Initialization, evaporation, precipitation, and river formation are key algorithm steps (23, 24).

During the initialization phase, a set of potential solutions is generated in a random manner. In this context, every solution is defined by a set of hyper-parameter values. For instance, in a function optimization scenario, these parameters may represent the input variable values (25).

The evaluation of the solutions’ fitness values takes place in the evaporation phase. The fitness of a solution reflects its quality, with higher fitness values corresponding to superior solutions. These fitness values play a role in computing the evaporation rate, dictating the amount of water that evaporates from each solution (25, 26).

In the precipitation stage, the evaporated water is converted into clouds, which are subsequently dispersed randomly among the solution population. Each cloud symbolizes a potential enhancement to a solution. The fitness values of the clouds are assessed, and the most superior cloud is selected (27). The introduced steps of WCA are reiterated until a termination condition is met. The stopping criterion could involve attaining an specific fitness value, reaching a maximum number of iterations, or consuming a maximum computational time (24, 28).

Figure 3 shows the basic workflow of WCA algorithm. One of the strengths of the WCA is its capability to address multiple objectives. Multi-objective optimization finds the best solutions for conflicting goals (accuracy and generality in ML tasks). The dominance principle is used to extend the WCA to multiple goals. A solution dominates another solution if it excels in at least one objective without being inferior in any. Utilizing this idea, the WCA can identify a collection of solutions that are independent of one another (28). Natural-inspired processes give the Water Cycle Algorithm (WCA) benefits over Particle Swarm Optimization (PSO) and Differential Evolution (DE). WCA’s iterative evaporation, cloud formation, and river construction balance exploration and exploitation to avoid local optima and increase convergence. Its self-adjusting evaporation rates and cloud dispersal improve performance without parameter adjustment. The robust diversification mechanism of WCA provides constant exploration of new regions, while its structured navigation of complex search spaces via river formation efficiently avoids suboptimal regions. These properties make WCA more adaptable and effective than PSO and DE.

Figure 3
www.frontiersin.org

Figure 3. Water cycle algorithm (WCA).

This algorithm is used for model optimization (hyper-parameter tuning) in this study. We used deterministic optimization here while there exists methods for optimization under uncertainty (29, 30). Here, each solution consists of a combination of hyper-parameter values and one of the objective functions is the RMSE error rate of the model build on each solution which should be maximized.

Also, by selecting the architecture with the least Akaike Information Criterion (AIC) value, the models are filtered to prevent overfitting, thus promoting generalization and robustness in the forecasting models. This method already shown promising results in avoiding overfitted models (31).

3.2 Piecewise polynomial regression

Piecewise Polynomial Regression (PRR) is the process of estimating a regression function by fitting multiple polynomial functions to different segments of the dataset. In this regression model, several polynomial functions are used to approximate the regression function in specific data segments (32).

PPR divides the dataset into segments and uses polynomial functions to approximate the regression function. The foundation of this segmentation is input space partitioning. Every segment is equivalent to a polynomial function that denotes the regression relationship inside that certain interval (32, 33).

To determine the optimal piecewise polynomial estimator, the research paper suggests considering various models defined by partitions and polynomial degrees. It utilizes a penalized least squares criterion to identify the model with an estimator closely resembling the best one in terms of quadratic risk. Additionally, the study establishes a non-asymptotic risk bound to ensure the selected model’s performance.

Extending the methodology to tree-structured partitions akin to those in the CART (34) algorithm offers a novel approach to constructing piecewise polynomial estimators for regression functions. This extension involves iteratively optimizing the selection of polynomial functions within each segment to best represent the underlying regression function across the entire data range.

3.3 Kernel ridge regression

The ideas of linear ridge regression are extended in the robust non-linear regression method known as kernel ridge regression (KRR). Utilizing the kernel trick, KRR adeptly captures the non-linear patterns inherent in data, rendering it applicable across diverse scenarios. The regularization mechanisms in KRR, including the ridge penalty, play a pivotal role in shaping the optimized kernel function, thereby mitigating overfitting concerns and bolstering predictive precision (35, 36).

Assume we have a dataset x i y i i = 1 N consisting of N data rows that have been sampled from a distribution ℙ over the Cartesian product of X and the real numbers (ℝ). The objective is to identify the best function that, with the expectation computed collectively over (X,Y) pairs, minimizes the Mean Squared Error (MSE) of the data f x y 2 .

The conditional mean f x := E Y | X = x is widely regarded as the most suitable function (37). Using a squared Hilbert norm penalty and the M-estimator with the lowest squares loss on the dataset is a feasible approach to forecast the unknown function f (37).

f ̂ : = argmin f H 1 N i = 1 N f x i y i 2 + f H 2 ,

Here, H denotes a reproducing kernel Hilbert space and ⋋ > 0 represents a regularization parameter. The kernel ridge regression estimate (referred to as KRR) serves as the estimator employed in the equation above (38).

3.4 Tweedie regression

The Tweedie regression model is a versatile tool for analyzing data that is non-negative, right-skewed, and has a high probability of being zero (39). It provides a single model that can effectively handle various types of continuous data automatically, making model selection easier during fitting.

The Tweedie regression model is predicated upon the assumption that the response variable Y follows a Tweedie distribution, denoted as Y T w p u ϕ p , where u represents the mean, ϕ denotes the dispersion parameter, and p signifies the power parameter. This model postulates a relationship between the mean and variance of Y, where the variance is proportional to the mean raised to the power parameter p, expressed as V a r Y = ϕ · u p . The regression parameters b are linked to the mean of Y through a specified link function g · , such that E Y = g X b , where X represents the design matrix. The likelihood function for the Tweedie regression model involves the probability density function of the Tweedie distribution, which entails complexity due to the presence of the power parameter p, necessitating numerical methods for its evaluation. Estimation of the regression parameters b and the dispersion parameter ϕ typically relies on maximum likelihood, quasi-likelihood, or pseudo-likelihood methods, leveraging second-moment assumptions to facilitate efficient and adaptable modeling of continuous data (39, 40).

4 Results and discussion

In this study, we implemented regression models to predict the solubility of ketoprofen in supercritical CO₂ using the provided dataset. The dataset includes measurements of temperature (T in Kelvin) and pressure (P in bar) as inputs, and solubility as the output. The regression models employed are PPR, KRR, and TDR. To optimize the hyperparameters of these models, we utilized the Water Cycle Algorithm (WCA), a robust optimization method known for its efficiency in finding optimal solutions in complex search spaces. The dataset was partitioned into training and testing subsets (80% training and 20% test) randomly, and the effectiveness of the models was gaged using a diverse set of metrics. The models are implemented on a machine with core-i7 CPU and 8Gb RAM which takes very small time (near Realtime) for each model to be executed. The effectiveness of regression models was appraised based on the following criteria (41):

1. R2 Score (Coefficient of Determination):

R 2 = 1 i = 1 n Y observed , i Y predicted , i 2 i = 1 n Y observed , i Y observed ¯ 2

These metric measures how well model predictions match data. An R2 score near 1 indicates a highly accurate model.

1. Mean squared error

M S E = 1 n i = 1 n Y observed , i Y predicted , i 2

The average squared difference between the outcomes as seen in reality and those the model predicts is known as MSE. Lower MSE values indicate better model performance.

1. Mean absolute error

M A E = 1 n i = 1 n Y observed , i Y predicted , i

Without taking into account the direction of the errors, MAE determines the average magnitude of errors in a set of predictions. Like MSE, lower MAE values indicate better model performance.

Table 2 summarizes the numerical findings for every regression model. The results indicate that Piecewise Polynomial Regression (PPR) outperforms both Kernel Ridge Regression (KRR) and Tweedie Regression (TDR) across all evaluation metrics. PPR achieved the highest R2 score of 0.97111, indicating that it explains approximately 97% of the variance in solubility. Additionally, PPR had the lowest MSE (1.6867E-09) and MAE (3.01040E-05), further demonstrating its superior predictive accuracy. Figure 4 shows a comparison of predicted solubility values and their corresponding actual values using the PPR model.

Table 2
www.frontiersin.org

Table 2. Performance metrics for the regression models predicting ketoprofen solubility.

Figure 4
www.frontiersin.org

Figure 4. Actual and predicted values comparison (PPR model).

The results obtained in this study showed superior performance compared to the previous machine learning models developed for prediction of drug solubility in supercritical CO2. The accuracy reported for machine learning modeling of Hyoscine solubility in supercritical CO2 was reported to be 2.1680E-01 for the best model in terms of RMSE which is higher than our model developed in this study (20). A RMSE value of 2.74912E-01 was reported by Almehizia et al. (42) for prediction of multiple drugs solubility in supercritical CO2. The best model was reported to be HS-PR (Harmony Search-Polynomial Regression).

Kernel Ridge Regression (KRR) also performed well, with an R2 score of 0.95044, an MSE of 2.5499E-09, and an MAE of 3.49707E-05. Although KRR’s performance was slightly inferior to PPR, it still showed strong predictive capabilities, making it a viable alternative for modeling ketoprofen solubility. Figure 5 compares the KRR model-predicted and corresponding actual solubility values.

Figure 5
www.frontiersin.org

Figure 5. Actual and predicted values comparison (KRR model).

Tweedie Regression (TDR) exhibited the lowest performance among the three models, with R2 score of 0.84413, an MSE of 7.4249E-09, and an MAE of 5.69159E-05. Despite its relatively lower performance, TDR provided reasonable predictions, but it is clear that more complex models like PPR and KRR better capture the relationship between input features and solubility. The comparison between predicted and actual solubility values using the TDR model is depicted in Figure 6. Although some deviations can be observed in the testing step for this model, the overall fitting accuracy is acceptable considering the complexity of the process and dataset. Moreover, the models have been optimized in a way to minimize the risk of overprediction for the test dataset.

Figure 6
www.frontiersin.org

Figure 6. Actual and predicted values comparison (TDR model).

Ultimately, the Piecewise Polynomial Regression (PPR) model shows excellent accuracy and dependability in predicting the solubility of ketoprofen in supercritical CO₂. Comparing this model to the other models assessed in this work, it performs better because it can partition the data space and fit polynomial functions within each segment to efficiently capture the underlying patterns in the data. Figures 7, 8 illustrate, with this model, how inputs affect the solubility values. Furthermore, shown in Figure 9 is the solubility in a contour plot and three-dimensional manner as a function of T(K) and P(bar). The trend for temperature shows exponential increase of solubility with temperature rise. On the other hand, a linear trend was observed for the influence of pressure on the drug solubility (see Figure 8). Thus, the maximum amount of ketoprofen solubility is determined at the highest values of T and P based on the ML models. Indeed, there is no optimum point for the solubility, and the optimum conditions should be determined from the process cost and economic evaluations.

Figure 7
www.frontiersin.org

Figure 7. Impact of temperature on the solubility on different pressure levels.

Figure 8
www.frontiersin.org

Figure 8. Impact of pressure on the solubility on different temperature levels.

Figure 9
www.frontiersin.org

Figure 9. Final PPR model: (A) the 3D representation of predicted solubility values (B) contour plot of predicted solubility values.

5 Conclusion

In this study, we successfully developed and evaluated three machine learning regression models–PPR, KRR, and TDR–to predict the solubility of ketoprofen in supercritical CO2. By employing the Water Cycle Algorithm for hyperparameter tuning, we optimized each model’s performance, demonstrating the effectiveness of this approach in enhancing predictive accuracy. Our comparative analysis revealed that PPR outperformed the other models, providing the most accurate predictions with an R2 score of 0.97111, a MSE of 1.6867 × 10−9, and a MAE of 3.01040 × 10−5. This paper makes multiple contributions. Initially, we assess the three regression models used for predicting ketoprofen solubility, emphasizing their advantages and disadvantages. Furthermore, we illustrate the enhancement of model performance through the optimization of hyperparameters using the Water Cycle Algorithm. Ultimately, we provide comprehensive visual representations and statistical examinations to uncover the connections between temperature (T), pressure (P), and solubility. This study employs sophisticated machine learning techniques and resilient optimization strategies to forecast the solubility of supercritical fluids and enhance the optimization of pharmaceutical processes.

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.

Author contributions

TH: Conceptualization, Validation, Formal analysis, Writing – Review & Editing. SA: Supervision, Validation, Investigation, Resources, Writing – Review & Editing. SB: Conceptualization, Formal analysis, Investigation, Methodology, Resources, Validation, Visualization, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by Taif University, Saudi Arabia, Project No. (TU-DSPP-2024-61).

Acknowledgments

The authors extend their appreciation to Taif University, Saudi Arabia, for supporting this work through project number (TU-DSPP-2024-61).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

1. Hyer, A, Gregory, D, Kay, K, le, Q, Turnage, J, Gupton, F, et al. Continuous manufacturing of active pharmaceutical ingredients: current trends and perspectives. Adv Synth Catal. (2024) 366:357–89. doi: 10.1002/adsc.202301137

Crossref Full Text | Google Scholar

2. Liu, P, Jin, H, Chen, Y, Wang, D, Yan, H, Wu, M, et al. Process analytical technologies and self-optimization algorithms in automated pharmaceutical continuous manufacturing. Chin Chem Lett. (2024) 35:108877. doi: 10.1016/j.cclet.2023.108877

Crossref Full Text | Google Scholar

3. Sundarkumar, V, Wang, W, Mills, M, Oh, SW, Nagy, Z, and Reklaitis, G. Developing a modular continuous drug product manufacturing system with real time quality Assurance for Producing Pharmaceutical Mini-Tablets. J Pharm Sci. (2024) 113:937–47. doi: 10.1016/j.xphs.2023.09.024

PubMed Abstract | Crossref Full Text | Google Scholar

4. Camacho Vieira, C, Peltonen, L, Karttunen, AP, and Ribeiro, AJ. Is it advantageous to use quality by design (QbD) to develop nanoparticle-based dosage forms for parenteral drug administration? Int J Pharm. (2024) 657:124163. doi: 10.1016/j.ijpharm.2024.124163

PubMed Abstract | Crossref Full Text | Google Scholar

5. Simão, J, Chaudhary, SA, and Ribeiro, AJ. Implementation of quality by design (QbD) for development of bilayer tablets. Eur J Pharm Sci. (2023) 184:106412. doi: 10.1016/j.ejps.2023.106412

PubMed Abstract | Crossref Full Text | Google Scholar

6. Barrera Jiménez, AA, Matsunami, K, van Hauwermeiren, D, Peeters, M, Stauffer, F, dos Santos Schultz, E, et al. Partial least squares regression to calculate population balance model parameters from material properties in continuous twin-screw wet granulation. Int J Pharm. (2023) 640:123040. doi: 10.1016/j.ijpharm.2023.123040

PubMed Abstract | Crossref Full Text | Google Scholar

7. Bellinghausen, S, Gavi, E, Jerke, L, Barrasso, D, Salman, AD, and Litster, JD. Model-driven design using population balance modelling for high-shear wet granulation. Powder Technol. (2022) 396:578–95. doi: 10.1016/j.powtec.2021.10.028

Crossref Full Text | Google Scholar

8. Muthancheri, I, Chaturbedi, A, Bétard, A, and Ramachandran, R. A compartment based population balance model for the prediction of steady and induction granule growth behavior in high shear wet granulation. Adv Powder Technol. (2021) 32:2085–96. doi: 10.1016/j.apt.2021.04.021

Crossref Full Text | Google Scholar

9. Tong, F, Wang, Y, and Gao, H. Progress and challenges in the translation of cancer nanomedicines. Curr Opin Biotechnol. (2024) 85:103045. doi: 10.1016/j.copbio.2023.103045

Crossref Full Text | Google Scholar

10. Xiao, Y, Zhong, L, Liu, J, Chen, L, Wu, Y, and Li, G. Progress and application of intelligent nanomedicine in urinary system tumors. J. Pharm. Anal. (2024):100964. doi: 10.1016/j.jpha.2024.100964

Crossref Full Text | Google Scholar

11. Zhang, S, Li, R, Jiang, T, Gao, Y, Zhong, K, Cheng, H, et al. Inhalable nanomedicine for lung cancer treatment. Smart Mater. Med. (2024) 5:261–80. doi: 10.1016/j.smaim.2024.04.001

Crossref Full Text | Google Scholar

12. Abdelbasset, WK, Elkholi, SM, Ahmed Ismail, K, Alalwani, TAA, Hachem, K, Mohamed, A, et al. Modeling and computational study on prediction of pharmaceutical solubility in supercritical CO2 for manufacture of nanomedicine for enhanced bioavailability. J Mol Liq. (2022) 359:119306. doi: 10.1016/j.molliq.2022.119306

Crossref Full Text | Google Scholar

13. Abourehab, MAS, Salah al-Shati, A, Venkatesan, K, Alshehri, S, Alzhrani, RM, Alsubaiyel, AM, et al. Theoretical investigations on the manufacture of drug nanoparticles using green supercritical processing: estimation and prediction of drug solubility in the solvent using advanced methods. J Mol Liq. (2022) 368:120559. doi: 10.1016/j.molliq.2022.120559

Crossref Full Text | Google Scholar

14. Bagheri, H, Notej, B, Shahsavari, S, and Hashemipour, H. Supercritical carbon dioxide utilization in drug delivery: experimental study and modeling of paracetamol solubility. Eur J Pharm Sci. (2022) 177:106273. doi: 10.1016/j.ejps.2022.106273

PubMed Abstract | Crossref Full Text | Google Scholar

15. Banchero, M, and Manna, L. Solubility of fenamate drugs in supercritical carbon dioxide by using a semi-flow apparatus with a continuous solvent-washing step in the depressurization line. J Supercrit Fluids. (2016) 107:400–7. doi: 10.1016/j.supflu.2015.10.008

Crossref Full Text | Google Scholar

16. Chinh Nguyen, H, Alamray, F, Kamal, M, Diana, T, Mohamed, A, Algarni, M, et al. Computational prediction of drug solubility in supercritical carbon dioxide: thermodynamic and artificial intelligence modeling. J Mol Liq. (2022) 354:118888. doi: 10.1016/j.molliq.2022.118888

Crossref Full Text | Google Scholar

17. Kostyrin, EV, Ponkratov, VV, and Salah Al-Shati, A. Development of machine learning model and analysis study of drug solubility in supercritical solvent for green technology development. Arab J Chem. (2022) 15:104346. doi: 10.1016/j.arabjc.2022.104346

Crossref Full Text | Google Scholar

18. Zhou, ZH . Machine learning. Berlin: Springer Nature (2021).

Google Scholar

19. Zhong, S, Zhang, K, Bagheri, M, Burken, JG, Gu, A, Li, B, et al. Machine learning: new ideas and tools in environmental science and engineering. Environ Sci Technol. (2021) 55:12741–54. doi: 10.1021/acs.est.1c01339

PubMed Abstract | Crossref Full Text | Google Scholar

20. Abouzied, AS, Alshahrani, SM, Hani, U, Obaidullah, AJ, al Awadh, AA, Lahiq, AA, et al. Assessment of solid-dosage drug nanonization by theoretical advanced models: modeling of solubility variations using hybrid machine learning models. Case Stud. Therm. Eng. (2023) 47:103101. doi: 10.1016/j.csite.2023.103101

Crossref Full Text | Google Scholar

21. Sadeghi, A, Su, CH, Khan, A, Lutfor Rahman, M, Sani Sarjadi, M, and Sarkar, SM. Machine learning simulation of pharmaceutical solubility in supercritical carbon dioxide: prediction and experimental validation for busulfan drug. Arab J Chem. (2022) 15:103502. doi: 10.1016/j.arabjc.2021.103502

Crossref Full Text | Google Scholar

22. Sabegh, MA, Rajaei, H, Esmaeilzadeh, F, and Lashkarbolooki, M. Solubility of ketoprofen in supercritical carbon dioxide. J Supercrit Fluids. (2012) 72:191–7. doi: 10.1016/j.supflu.2012.08.008

Crossref Full Text | Google Scholar

23. Abou El-Ela, AA, El-Sehiemy, RA, and Abbas, AS. Optimal placement and sizing of distributed generation and capacitor banks in distribution systems using water cycle algorithm. IEEE Syst J. (2018) 12:3629–36. doi: 10.1109/JSYST.2018.2796847

Crossref Full Text | Google Scholar

24. Ghazwani, M, and Begum, MY. Computational intelligence modeling of hyoscine drug solubility and solvent density in supercritical processing: gradient boosting, extra trees, and random forest models. Sci Rep. (2023) 13:10046. doi: 10.1038/s41598-023-37232-8

PubMed Abstract | Crossref Full Text | Google Scholar

25. Eskandar, H, Sadollah, A, Bahreininejad, A, and Hamdi, M. Water cycle algorithm–A novel metaheuristic optimization method for solving constrained engineering optimization problems. Comput Struct. (2012) 110-111:151–66. doi: 10.1016/j.compstruc.2012.07.010

Crossref Full Text | Google Scholar

26. Sadollah, A, Eskandar, H, and Kim, JH. Water cycle algorithm for solving constrained multi-objective optimization problems. Appl Soft Comput. (2015) 27:279–98. doi: 10.1016/j.asoc.2014.10.042

Crossref Full Text | Google Scholar

27. Razmjooy, N, Khalilpour, M, and Ramezani, M. A new meta-heuristic optimization algorithm inspired by FIFA world cup competitions: theory and its application in PID designing for AVR system. J Control Autom Electr Syst. (2016) 27:419–40. doi: 10.1007/s40313-016-0242-6

Crossref Full Text | Google Scholar

28. Jafar, R. M. S., Geng, S., Ahmad, W., Hussain, S., and Wang, H. A comprehensive evaluation: water cycle algorithm and its applications. In: Bio-inspired computing: Theories and applications: 13th international conference, BIC-TA 2018, Beijing, China, November 2–4, 2018, proceedings, part II 13. Berlin: Springer. (2018).

Google Scholar

29. Sharma, S, Mahadevan, J, Giri, L, and Mitra, K. Identification of optimal flow rate for culture media, cell density, and oxygen toward maximization of virus production in a fed-batch baculovirus-insect cell system. Biotechnol Bioeng. (2023) 120:3529–42. doi: 10.1002/bit.28558

PubMed Abstract | Crossref Full Text | Google Scholar

30. Sharma, S, Keerthi, PN, Giri, L, and Mitra, K. Toward performance improvement of a Baculovirus–insect cell system under uncertain environment: A robust multiobjective dynamic optimization approach for Semibatch suspension culture. Ind Eng Chem Res. (2022) 62:111–25. doi: 10.1021/acs.iecr.2c03355

Crossref Full Text | Google Scholar

31. Pujari, KN, Miriyala, SS, Mittal, P, and Mitra, K. Better wind forecasting using evolutionary neural architecture search driven green deep learning. Expert Syst Appl. (2023) 214:119063. doi: 10.1016/j.eswa.2022.119063

Crossref Full Text | Google Scholar

32. Sauve, M . Piecewise polynomial estimation of a regression function. IEEE Trans Inf Theory. (2009) 56:597–613. doi: 10.1109/TIT.2009.2027481

Crossref Full Text | Google Scholar

33. Chaudhuri, P, Huang, MC, Loh, WY, Yao, R, et al. Piecewise-polynomial regression trees. Stat Sin. (1994) 4:143–67.

Google Scholar

34. Breiman, L, Friedman, JH, Olshen, RA, and Stone, CJ. Classification and regression trees. Abingdon: Routledge (2017).

Google Scholar

35. Tandon, R, Si, S, Ravikumar, P, and Dhillon, I. Kernel ridge regression via partitioning. arXiv (2016).

Google Scholar

36. Welling, M . Kernel ridge regression. Max Welling’s classnotes in machine learning, pp. 1–3. (2013).

Google Scholar

37. Byrne, E, and Schniter, P. Sparse multinomial logistic regression via approximate message passing. IEEE Trans Signal Process. (2016) 64:5485–98. doi: 10.1109/TSP.2016.2593691

Crossref Full Text | Google Scholar

38. Zhang, Y, Duchi, J, and Wainwright, M. Divide and conquer kernel ridge regression. In conference on learning theory. Proc Machine Learn Res. 206:6245–62.

PubMed Abstract | Google Scholar

39. Bonat, WH, and Kokonendji, CC. Flexible Tweedie regression models for continuous data. J Stat Comput Simul. (2017) 87:2138–52. doi: 10.1080/00949655.2017.1318876

Crossref Full Text | Google Scholar

40. Hassine, A, Masmoudi, A, and Ghribi, A. Tweedie regression model: a proposed statistical approach for modelling indoor signal path loss. Int J Numer Modell. (2017) 30:e2243. doi: 10.1002/jnm.2243

Crossref Full Text | Google Scholar

41. Botchkarev, A . Performance metrics (error measures) in machine learning regression, forecasting and prognostics: Properties and typology. arXiv, (2018).

Google Scholar

42. Almehizia, AA, Naglah, AM, Alkahtani, HM, Hani, U, and Ghazwani, M. Numerical optimization of drug solubility inside the supercritical carbon dioxide system using different machine learning models. J Mol Liq. (2023) 392:123466. doi: 10.1016/j.molliq.2023.123466

Crossref Full Text | Google Scholar

Keywords: drug development, solubility prediction, optimization, machine learning, modeling

Citation: Al Hagbani T, Alshehri S and Bawazeer S (2024) Advanced modeling of pharmaceutical solubility in solvents using artificial intelligence techniques: assessment of drug candidate for nanonization processing. Front. Med. 11:1435675. doi: 10.3389/fmed.2024.1435675

Received: 20 May 2024; Accepted: 26 June 2024;
Published: 22 July 2024.

Edited by:

Ovidiu Constantin Baltatu, Anhembi Morumbi University, Brazil

Reviewed by:

Kishalay Mitra, Indian Institute of Technology Hyderabad, India
A.M. Elsawah, Beijing Normal University–Hong Kong Baptist University United International College, China
Mahboubeh Pishnamazi, Duy Tan University, Vietnam

Copyright © 2024 Al Hagbani, Alshehri and Bawazeer. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Turki Al Hagbani, T.alhagbani@uoh.edu.sa

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.