The successful development of nonviral delivery systems for nucleic acids has been reported extensively over the past years. Increasingly employed to improve the delivery efficiency and therapeutic efficacy of RNA are lipid nanoparticles (LNPs). Many of the various critical formulation parameters can affect the quality attributes and effectiveness of these nano-formulations. Therefore, the systematic drug development approach (QbD) and multivariate design and statistical analysis (DOE) can be very helpful and recommended for the optimization of the composition and production of RNA–LNPs. This review addresses the concepts and applications of QbD and/or DOE for the development of lipid nanoparticles for the delivery of different types of RNA, reporting examples published in the ten recent years presenting the latest trends and regulatory requirements as well as the modern mathematical and statistical design methods. As the topic explored in this review is a novel approach, the full QbD has been described in only a few papers, and a few refer only to some aspects of QbD. In contrast, the DOE approach has been used in most of the optimization works. Different approaches and innovations in DOE have been observed. Traditional statistical tests and modeling (ANOVA, regression analysis) are slowly being replaced by artificial intelligence and machine learning methods.

TrkB is a tyrosine kinase receptor that is activated upon binding to brain-derived neurotrophic factor (BDNF). To date, the search for low-molecular-weight molecules mimicking BDNF’s action has been unsuccessful. Several molecules exerting antidepressive effects in vivo, such as 7,8-DHF, have been suggested to be TrkB agonists. However, more recent publications question this hypothesis. In this study, we developed a set of experimental procedures including the evaluation of direct interactions, dimerization, downstream signaling, and cytoprotection in parallel with physicochemical and ADME methods to verify the pharmacology of 7,8-DHF and other potential reference compounds, and perform screening for novel TrkB agonists. 7,8 DHF bound to TrkB with K_d = 1.3 μM; however, we were not able to observe any other activity against the TrkB receptor in SN56 T48 and differentiated SH-SY5Y cell lines. Moreover, the pharmacokinetic and pharmacodynamic effects of 7,8-DHF at doses of 1 and 50 mg/kg were examined in mice after i.v and oral administration, respectively. The poor pharmacokinetic properties and lack of observed activation of TrkB-dependent signaling in the brain confirmed that 7,8-DHF is not a relevant tool for studying TrkB activation in vivo. The binding profile for 133 molecular targets revealed a significant lack of selectivity of 7,8-DHF, suggesting a distinct functional profile independent of interaction with TrkB. Additionally, a compound library was screened in search of novel low-molecular-weight orthosteric TrkB agonists; however, we were not able to identify reliable drug candidates. Our results suggest that published reference compounds including 7,8-DHF do not activate TrkB, consistent with canonical dogma, which indicates that the reported pharmacological activity of these compounds should be interpreted carefully in a broad functional context.

The purpose of this work was to demonstrate the use of the AQbD with the DOE approach to the methodical step-by-step development of a UHPLC method for the quantitative determination of the impurity profile of new CPL409116 substance (JAK/ROCK inhibitor) on the preclinical and clinical step of drug discovery studies. The critical method parameters (CMPs) have been tested extensively: the kind of stationary phase (8 different columns), pH of the aqueous mobile phase (2.6, 3.2, 4.0, 6.8), and start (20–25%) and stop (85–90%) percentage of organic mobile phase (ACN). The critical method attributes (CMAs) are the resolution between the peaks (≥2.0) and peak symmetry of analytes (≥0.8 and ≤1.8). In the screening step, the effects of different levels of CMPs on the CMAs were evaluated based on a full fractional design 2². The robustness tests were established from the knowledge space of the screening step and performed by application fractional factorial design 2⁽⁴⁻¹⁾. Method operable design region (MODR) was generated. The probability of meeting the specifications for the CMAs was calculated by Monte-Carlo simulations. In relation to literature such a complete AQbD approach including screening, optimization, and validation steps for the development of a new method for the quantitative determination of the full profile of nine impurities of an innovative pharmaceutical substance with the structure-based pre-development pointed out the novelty of our work. The final working conditions were as follows: column Zorbax Eclipse Plus C18, aqueous mobile phase 10 mM ± 1 mM aqueous solution of HCOOH, pH 2.6, 20% ± 1% of ACN at the start and 85% ± 1% of ACN at the end of the gradient, and column temperature 30 °C ± 2 °C. The method was validated in compliance with ICH guideline Q2(R1). The optimized method is specified, linear, precise, and robust. LOQ is on the reporting threshold level of 0.05% and LOD at 0.02% for all impurities.

Phosphoinositide 3-kinase (PI3K) is the family of lipid kinases participating in vital cellular processes such as cell proliferation, growth, migration, or cytokines production. Due to the high expression of these proteins in many human cells and their involvement in metabolism regulation, normal embryogenesis, or maintaining glucose homeostasis, the inhibition of PI3K (especially the first class which contains four subunits: α, β, γ, δ) is considered to be a promising therapeutic strategy for the treatment of inflammatory and autoimmune diseases such as systemic lupus erythematosus (SLE) or multiple sclerosis. In this work, we synthesized a library of benzimidazole derivatives of pyrazolo[1,5-a]pyrimidine representing a collection of new, potent, active, and selective inhibitors of PI3Kδ, displaying IC₅₀ values ranging from 1.892 to 0.018 μM. Among all compounds obtained, CPL302415 (6) showed the highest activity (IC₅₀ value of 18 nM for PI3Kδ), good selectivity (for PI3Kδ relative to other PI3K isoforms: PI3Kα/δ = 79; PI3Kβ/δ = 1415; PI3Kγ/δ = 939), and promising physicochemical properties. As a lead compound synthesized on a relatively large scale, this structure is considered a potential future candidate for clinical trials in SLE treatment.

Herein, we describe the development of a green, scalable flow Pd-catalyzed aerobic oxidation for the key step in the synthesis of CPL302415, which is a new PI3Kδ inhibitor. Applying this environmental-friendly, sustainable catalytic oxidation we significantly increased product yield (up to 84%) and by eliminating of workup step, we improved the waste index and E factor (up to 0.13) in comparison with the stoichiometric synthesis. The process was optimized by using the DoE approach.

Phosphoinositide 3-kinase δ (PI3Kδ), a member of the class I PI3K family, is an essential signaling biomolecule that regulates the differentiation, proliferation, migration, and survival of immune cells. The overactivity of this protein causes cellular dysfunctions in many human disorders, for example, inflammatory and autoimmune diseases, including asthma or chronic obstructive pulmonary disease (COPD). In this work, we designed and synthesized a new library of small-molecule inhibitors based on indol-4-yl-pyrazolo[1,5-a]pyrimidine with IC₅₀ values in the low nanomolar range and high selectivity against the PI3Kδ isoform. CPL302253 (54), the most potent compound of all the structures obtained, with IC₅₀ = 2.8 nM, is a potential future candidate for clinical development as an inhaled drug to prevent asthma.