Employing Hexahydroquinolines as PfCDPK4 Inhibitors to Combat Malaria Transmission: An Advanced Computational Approach

Background

Existing antimalarial drugs primarily target blood-stage parasites, but there is a need for transmission-blocking drugs to combat malaria effectively. Plasmodium falciparum Calcium-dependent Protein Kinase 4 (CDPK4) is a promising target for such drugs. This study employed advanced in silico analyses of hexahydroquinolines (HHQ) derivatives to identify PfCDPK4 inhibitors capable of disrupting malaria transmission. Structure-based virtual screening (SBVS) was employed to discover HHQ derivatives with the highest binding affinities against the 3D structure of PfCDPK4 (PDB 1D: 4QOX).

Methods

Interaction analysis of protein-ligand complexes utilized Discovery Studio Client, while druglikeness and ADMET properties were assessed using SwissADME and pkCSM web servers, respectively. Quantum mechanical calculations of the top hits were conducted using density functional theory (DFT), and GROMACS was employed to perform the molecular dynamics (MD) simulations. Binding free energy was predicted using the MMPBSA.py tool from the AMBER package.

Results

SBVS identified ten best hits possessing docking scores within the range of -11.2 kcal/mol and -10.6 kcal/mol, surpassing the known inhibitor, BKI-1294 (-9.9 kcal/mol). Among these, 4-[4-(Furan-2-carbonyl)piperazin-1-yl]-1-(naphthalen-2-ylmethyl)-2-oxo-4a,5,6,7,8,8a-hexahydroquinoline-3-carbonitrile (PubChem ID: 145784778) exhibited the highest binding affinity (-11.2 kcal/mol) against PfCDPK4.

Conclusion

Comparative analysis of this compound with BKI-1294 using advanced computational approaches demonstrated competitive potential. These findings suggest the potential of 4-[4-(Furan-2-carbonyl)piperazin-1-yl]-1-(naphthalen-2-ylmethyl)-2-oxo-4a,5,6,7,8,8a-hexahydroquinoline-3-carbonitrile as a promising PfCDPK4 inhibitor for disrupting malaria transmission. However, further experimental studies are warranted to validate its efficacy and safety profile.

Antiangiogenic potential of phytochemicals from Clerodendrum inerme (L.) Gaertn investigated through in silico and quantum computational methods

Suppressing vascular endothelial growth factor (VEGF), its receptor (VEGFR2), and the VEGF/VEGFR2 signaling cascade system to inhibit angiogenesis has emerged as a possible cancer therapeutic target. The present work was designed to discover and evaluate bioactive phytochemicals from the Clerodendrum inerme (L.) Gaertn plant for their anti-angiogenic potential. Molecular docking of twenty-one phytochemicals against the VEGFR-2 (PDB ID: 3VHE) protein was performed, followed by ADMET profiling and molecular docking simulations. These investigations unveiled two hit compounds, cirsimaritin (- 12.29 kcal/mol) and salvigenin (- 12.14 kcal/mol), with the highest binding energy values when compared to the reference drug, Sorafenib (- 15.14 kcal/mol). Furthermore, only nine phytochemicals (cirsimaritin and salvigenin included) obeyed Lipinski's rule of five and passed ADMET filters. Molecular dynamics simulations run over 100 ns revealed that the protein-ligand complexes remained stable with minimal backbone fluctuations. The binding free energy values of cirsimaritin (- 52.35 kcal/mol) and salvigenin (- 55.89 kcal/mol), deciphered by MM-GBSA analyses, further corroborated the docking interactions. The HOMO-LUMO band energy gap (ΔE) was calculated using density-functional theory (DFT) and substantiated using density of state (DOS) spectra. The chemical reactivity analyses revealed that salvigenin exhibited the highest chemical softness value (6.384 eV), the lowest hardness value (0.07831 eV), and the lowest ΔE value (0.1566 eV), which implies salvigenin was less stable and chemically more reactive than cirsimaritin and sorafenib. These findings provide further evidence that cirsimaritin and salvigenin have the ability to prevent angiogenesis and the development of cancer. Nevertheless, more in vitro and in vivo confirmation is necessary.

Computational Approach to Drug Penetration across the Blood-Brain and Blood-Milk Barrier Using Chromatographic Descriptors

Drug penetration through biological barriers is an important aspect of pharmacokinetics. Although the structure of the blood-brain and blood-milk barriers is different, a connection can be found in the literature between drugs entering the central nervous system (CNS) and breast milk. This study was created to reveal such a relationship with the use of statistical modelling. The basic physicochemical properties of 37 active pharmaceutical compounds (APIs) and their chromatographic retention data (TLC and HPLC) were incorporated into calculations as molecular descriptors (MDs). Chromatography was performed in a thin layer format (TLC), where the plates were impregnated with bovine serum albumin to mimic plasma protein binding. Two columns were used in high performance liquid chromatography (HPLC): one with immobilized human serum albumin (HSA), and the other containing an immobilized artificial membrane (IAM). Statistical methods including multiple linear regression (MLR), cluster analysis (CA) and random forest regression (RF) were performed with satisfactory results: the MLR model explains 83% of the independent variable variability related to CNS bioavailability; while the RF model explains up to 87%. In both cases, the parameter related to breast milk penetration was included in the created models. A significant share of reversed-phase TLC retention values was also noticed in the RF model.

RP-18 TLC Chromatographic and Computational Study of Skin Permeability of Steroids

The skin permeability of steroids, as investigated in this study, is important because some of these compounds are, or could, be used in preparations applied topically. Several models of skin permeability, involving thin layer chromatographic and calculated descriptors, were generated and validated using K_p reference values obtained in silico and then tested on a group of solutes whose experimental K_p values could be found (log K_p^exp). The study established that the most applicable log K_p model is based on RP-18 thin layer chromatographic data (R_M) and the calculated descriptors V_M (molar volume) and PSA (polar surface area). Two less efficient, yet simple, equations based on PSA or V_M combined with HD (H-donor count) can be used with caution for rapid, rough estimations of compounds’ skin permeability prior to their chemical synthesis.

Assessment of Lipophilicity Descriptors of Selected NSAIDs Obtained at Different TLC Stationary Phases

Lipophilicity study of selected NSAIDs, the group of the bioactive compounds usually used in humans and animals medicine, with the use of experimental and calculation methods was evaluated. LogP values are proposed and compared as descriptors of the lipophilicity of eleven compounds (from oxicams and coxibs). Obtained data were designated by thin-layer chromatography (TLC) in various chromatographic conditions, with stationary phases with different properties. The mobile phase systems were prepared by mixing the respective amounts of water and organic modifier, methanol and acetone, in the range of 30 to 80% (v/v) in 5% increments. Retention parameters (RF, RM and RM0) were calculated and statistically evaluated to establish correlations. All experimentally determined RM0 values were compared with partition coefficients obtained by computational methods using linear regression analysis. Moreover, in order to extract information about the lipophilicity of compounds from large retention datasets, two chemometric approaches, namely principal component analysis (PCA) and cluster analysis (CA) were carried out. Established models of lipophilicity may have the potential to predict the biological activity of a number of drugs. The presented knowledge may also be of use during drug discovery processes, broadening the knowledge of potential ways to modify the physicochemical properties of chemical compounds.

Application of RP-18 TLC Retention Data to the Prediction of the Transdermal Absorption of Drugs

Several chromatographic parameters (R_M⁰ and S obtained from RP-18 TLC with methanol—pH 7.4 phosphate buffer mobile phases by extrapolation to zero concentration of methanol; R_f and R_M obtained from RP-18 TLC with acetonitrile—pH 7.4 phosphate buffer 70:30 v/v as a mobile phase) and calculated molecular descriptors (molecular weight—M_W; molar volume—V_M; polar surface area—PSA; total count of nitrogen and oxygen atoms—(N+O); H-bond donor count—HD; H-bond acceptor count—HA; distribution coefficient—log D; total energy—E_T; binding energy—E_b; hydration energy—E_h; energy of the highest occupied molecular orbital—E_HOMO; energy of the lowest unoccupied orbital—E_LUMO; electronic energy—E_e; surface area—S_a; octanol-water partition coefficient—log P; dipole moment—DM; refractivity—R, polarizability—α) and their combinations (R_f/PSA, R_M/M_W, R_M/V_M) were tested in order to generate useful models of solutes’ skin permeability coefficient log K_p. It was established that neither R_M⁰ nor S obtained in the conditions used in this study is a good predictor of the skin permeability coefficient. The chromatographic parameters R_f and R_f/PSA were also unsuitable for this purpose. A simple and potentially useful, purely computational model based on (N+O), log D and HD as independent variables and accounting for ca. 83% of total variability was obtained. The evaluation of parameters derived from R_M (R_M, R_M/M_W, R_M/V_M) as independent variables in log K_p models proved that R_M/V_M is the most suitable descriptor belonging to this group. In a search for a reliable log K_p model based on this descriptor two possibilities were considered: a relatively simple model based on 5 independent variables: (N+O), log D, R_M/V_M, E_T and E_h and a more complex one, involving also E_b, M_W and PSA.

Towards Deep Neural Network Models for the Prediction of the Blood-Brain Barrier Permeability for Diverse Organic Compounds

Permeation through the blood–brain barrier (BBB) is among the most important processes controlling the pharmacokinetic properties of drugs and other bioactive compounds. Using the fragmental (substructural) descriptors representing the occurrence number of various substructures, as well as the artificial neural network approach and the double cross-validation procedure, we have developed a predictive in silico LogBB model based on an extensive and verified dataset (529 compounds), which is applicable to diverse drugs and drug-like compounds. The model has good predictivity parameters (Q2=0.815, RMSEcv=0.318) that are similar to or better than those of the most reliable models available in the literature. Larger datasets, and perhaps more sophisticated network architectures, are required to realize the full potential of deep neural networks. The analysis of fragment contributions reveals patterns of influence consistent with the known concepts of structural characteristics that affect the BBB permeability of organic compounds. The external validation of the model confirms good agreement between the predicted and experimental LogBB values for most of the compounds. The model enables the evaluation and optimization of the BBB permeability of potential neuroactive agents and other drug compounds.

Comparison of statistical methods for predicting penetration capacity of drugs into human breast milk using physicochemical, pharmacokinetic and chromatographic descriptors

In silico methods are often used for predicting pharmacokinetic properties of drugs due to their simplicity and cost-effectiveness. This study evaluates the penetration of 83 active pharmaceutical ingredients into human breast milk with an experimental milk-to-plasma ratio (M/P) obtained from the literature. Multiple linear regression (MLR), partial least squares (PLS) and random forest (RF) regression methods were compared to uncover the relationship between physicochemical, pharmacokinetic and membrane crossing properties of active pharmaceutical ingredients (APIs) using their rapid reference measurement value (R_f values), thin-layer chromatography (TLC) data from albumin-impregnated plates. Molecular descriptors of APIs proven to be important for their crossing into breast milk, including protein binding, ionisation state and lipophilicity and TLC data, have been included in the development of the prediction models. The best regression results were achieved by MLR (r ² = 0.83 and r ² = 0.86, n = 28) and RF (r ² = 0.85, n = 58). In addition, the discriminant function analysis (DFA) was performed on acidic, basic and neutral drugs separately and showed a prediction accuracy of 93%, with M/P included as the discriminating variable.

Application of separation methods for in vitro prediction of blood-brain barrier permeability-The state of the art

Despite many efforts, drug discovery pipeline is still a highly inefficient process. Nowadays, when combinatorial chemistry enables to synthesize hundreds of new drugs candidates, methods for rapid assessment of biopharmaceutical parameters of new compounds are highly desired. Over one-third of drugs candidates is rejected because of unsatisfactory pharmacokinetic properties. In the drug discovery process, the blood-brain barrier (BBB) permeability plays a critical role for central nervous system active drugs candidates as well as non-central nervous system active drugs. For this reason, knowledge on the BBB permeability of compounds is essential in the development of new medicines. The review was focused on the application of different separation methods for BBB permeability assessment. Both chromatographic and electrophoretic methods were described. In the article, the advantages and limitations of well-established chromatographic methods like immobilized artificial membrane chromatography or micellar liquid chromatography, and less common techniques were discussed. Special attention was devoted to methods were microemulsion is used as mobile or pseudostationary phases.

Design, synthesis and pharmacological evaluation of CVIB, a codrug of carvacrol and ibuprofen as a novel anti-inflammatory agent

The search for new drugs with anti-inflammatory properties remains a challenge for modern medicine. Among the various strategies for drug discovery, deriving new chemical entities from known bioactive natural and/or synthetic compounds remains a promising approach. Here, we designed and synthesized CVIB, a codrug developed by association of carvacrol (a phenolic monoterpene) with ibuprofen (a non-steroidal anti-inflammatory drug). In silico pharmacokinetic and physicochemical properties evaluation indicated low aqueous solubility (LogP ≥5.0). Nevertheless, the hybrid presented excellent oral bioavailability, gastrointestinal tract absorption, and low toxicity. CVIB did not present cytotoxicity in peripheral blood mononuclear cells (PBMCs), and promoted a significant reduction in IL-2, IL-10, IL-17, and IFN-γ cytokine levels in vitro. The LD₅₀ was estimated to be approximately 5000 mg/kg. CVIB was stable and detectable in human plasma after 24 h. In vivo anti-inflammatory evaluations revealed that CVIB at 10 and 50 mg/kg i.p. caused a significant decrease in total leukocyte count (p < 0.01) and provoked a significant reduction in IL-1β (p < 0.01). CVIB at 10 mg/kg i.p. efficiently decreased inflammatory parameters better than the physical mixture (carvacrol + ibuprofen 10 mg/kg i.p.). The results suggest that the codrug approach is a good option for drug design and development, creating the possibility of combining NSAIDs with natural products in order to obtain new hybrid drugs may be useful for treatment of inflammatory diseases.

RP-18 HPLC Analysis of Drugs’ Ability to Cross the Blood-Brain Barrier

One hundred ten compounds of diverse structures (actives and excipients used in pharmaceutical preparations) were studied by RP-18 HPLC with acetonitrile-pH 7.4 phosphate buffer 1 : 1 (v/v) as the mobile phase. The relationships between the BBB permeation coefficients and the chromatographic parameters log k and (log k)/PSA were compared to those between the blood-brain barrier (BBB) permeation parameters and the RP-18 TLC descriptors Rf and Rf/PSA known from our earlier studies. It was found that the correlations between the BBB permeability and the HPLC data are slightly worse than those achieved for the thin-layer chromatographic data. MLR analysis based upon the physicochemical data confirmed the value of the molecular descriptors, related to the CNS bioavailability. These variables, combined with the HPLC data, made it possible to generate computational models, explaining 70–96% of the total variance of the CNS bioavailability. Contrary to TLC Rf, the advantage of the modification of HPLC log k with PSA (polar surface area) has not been confirmed and the results obtained with log k are superior to those obtained after a novel (log k)/PSA parameter has been introduced. Establishing a firm threshold limit of (log k)/PSA, log k, or even k and k/PSA to distinguish between the CNS+ and CNS− compounds was impossible. On the other hand, discriminant function analyses involving log k and (log k)/PSA as discriminating variables separated the CNS+ and CNS− compounds with the success rate ca. 90%. On the basis of these results, it was concluded that the RP-18 HPLC analytical models are entirely successful in studies and predictions of the BBB permeability.