While the COVID-19 pandemic continues to worsen, effective medicines that target the life cycle of SARS-CoV-2 are still under development. As more highly infective and dangerous variants of the coronavirus emerge, the protective power of vaccines will decrease or vanish. Thus, the development of drugs, which are free of drug resistance is direly needed. The aim of this study is to identify allosteric binding modulators from a large compound library to inhibit the binding between the Spike protein of the SARS-CoV-2 virus and human angiotensin-converting enzyme 2 (hACE2). The binding of the Spike protein to hACE2 is the first step of the infection of host cells by the coronavirus. We first built a compound library containing 77 448 antiviral compounds. Molecular docking was then conducted to preliminarily screen compounds which can potently bind to the Spike protein at two allosteric binding sites. Next, molecular dynamics simulations were performed to accurately calculate the binding affinity between the spike protein and an identified compound from docking screening and to investigate whether the compound can interfere with the binding between the Spike protein and hACE2. We successfully identified two possible drug binding sites on the Spike protein and discovered a series of antiviral compounds which can weaken the interaction between the Spike protein and hACE2 receptor through conformational changes of the key Spike residues at the Spike–hACE2 binding interface induced by the binding of the ligand at the allosteric binding site. We also applied our screening protocol to another compound library which consists of 3407 compounds for which the inhibitory activities of Spike/hACE2 binding were measured. Encouragingly, in vitro data supports that the identified compounds can inhibit the Spike–ACE2 binding. Thus, we developed a promising computational protocol to discover allosteric inhibitors of the binding of the Spike protein of SARS-CoV-2 to the hACE2 receptor, and several promising allosteric modulators were discovered.

Structure-based virtual screening utilizes molecular docking to explore and analyze ligand–macromolecule interactions, crucial for identifying and developing potential drug candidates. Although there is availability of several widely used docking programs, the accurate prediction of binding affinity and binding mode still presents challenges. In this study, we introduced a novel protocol that combines our in-house geometry optimization algorithm, the conjugate gradient with backtracking line search (CG-BS), which is capable of restraining and constraining rotatable torsional angles and other geometric parameters with a highly accurate machine learning potential, ANI-2x, renowned for its precise molecular energy predictions reassembling the wB97X/6-31G(d) model. By integrating this protocol with binding pose prediction using the Glide, we conducted additional structural optimization and potential energy prediction on 11 small molecule–macromolecule and 12 peptide–macromolecule systems. We observed that ANI-2x/CG-BS greatly improved the docking power, not only optimizing binding poses more effectively, particularly when the RMSD of the predicted binding pose by Glide exceeded around 5 Å, but also achieving a 26% higher success rate in identifying those native-like binding poses at the top rank compared to Glide docking. As for the scoring and ranking powers, ANI-2x/CG-BS demonstrated an enhanced performance in predicting and ranking hundreds or thousands of ligands over Glide docking. For example, Pearson’s and Spearman’s correlation coefficients remarkedly increased from 0.24 and 0.14 with Glide docking to 0.85 and 0.69, respectively, with the addition of ANI-2x/CG-BS for optimizing and ranking small molecules binding to the bacterial ribosomal aminoacyl-tRNA receptor. These results suggest that ANI-2x/CG-BS holds considerable potential for being integrated into virtual screening pipelines due to its enhanced docking performance.

The “Long-COVID syndrome” has posed significant challenges due to a lack of validated therapeutic options. We developed a novel multi-step virtual screening strategy to reliably identify inhibitors against 3-chymotrypsin-like protease of SARS-CoV-2 from abundant flavonoids, which represents a promising source of antiviral and immune-boosting nutrients. We identified 57 interacting residues as contributors to the protein-ligand binding pocket. Their energy interaction profiles constituted the input features for Machine Learning (ML) models. The consensus of 25 classifiers trained using various ML algorithms attained 93.9% accuracy and a 6.4% false-positive-rate. The consensus of 10 regression models for binding energy prediction also achieved a low root-mean-square error of 1.18 kcal/mol. We screened out 120 flavonoid hits first and retained 50 drug-like hits after predefined ADMET filtering to ensure bioavailability and safety profiles. Furthermore, molecular dynamics simulations prioritized nine bioactive flavonoids as promising anti-SARS-CoV-2 agents exhibiting both high structural stability (root-mean-square deviation < 5 Å for 218 ns) and low MM/PBSA binding free energy (<−6 kcal/mol). Among them, KB-2 (PubChem-CID, 14630497) and 9-O-Methylglyceofuran (PubChem-CID, 44257401) displayed excellent binding affinity and desirable pharmacokinetic capabilities. These compounds have great potential to serve as oral nutraceuticals with therapeutic and prophylactic properties as care strategies for patients with long-COVID syndrome.

Malaria is a severe parasite infectious disease with high fatality. As one of the approved treatments of this disease, hydroxychloroquine (HCQ) lacks clinical administration guidelines for patients with special health conditions and co-morbidities. This may result in improper dosing for different populations and lead them to suffer from severe side effects. One of the most important toxicities of HCQ overdose is cardiotoxicity. In this study, we built and validated a physiologically based pharmacokinetic modeling (PBPK) model for HCQ. With the full-PBPK model, we predicted the pharmacokinetic (PK) profile for malaria patients without other co-morbidities under the HCQ dosing regimen suggested by Food and Drug Administration (FDA) guidance. The PK profiles for different special populations were also predicted and compared to the normal population. Moreover, we proposed a series of adjusted dosing regimens for different populations with special health conditions and predicted the concentration-time (C-T) curve of the drug plasma concentration in these populations which include the pregnant population, elderly population, RA patients, and renal impairment populations. The recommended special population-dependent dosage regimens can maintain the similar drug levels observed in the virtual healthy population under the original dosing regimen provided by FDA. Last, we developed mathematic formulas for predicting dosage based on a patient’s body measurements and two indexes of renal function (glomerular filtration rate and serum creatine level) for the pediatric and morbidly obese populations. Those formulas can facilitate personalized treatment of this disease. We hope to provide some advice to clinical practice when taking HCQ as a treatment for malaria patients with special health conditions or co-morbidities so that they will not suffer from severe side effects due to higher drug plasma concentration, especially cardiotoxicity.

Pharmacokinetic (PK) data analysis in drug discovery is challenged by the inherent variability of experimental and clinical study designs, which hinders data integration and predictive modeling. To address this, we have curated a comprehensive, human-derived, and machine learning (ML)-ready PK dataset from authoritative, multiedition sources. This novel resource represents a systematically curated, human-derived PK dataset that integrates compound structures, clinical study design information, and experimental variability annotations, providing a structured foundation for data-driven analysis and modeling of human pharmacokinetics. We first implemented a rigorous standardization and filtering protocol to prepare the ML-ready dataset, and we demonstrated its utility through chemoinformatic analysis and ML classification model evaluations at across multiple classification systems. Fragment analysis showed a clear association between molecular structure and PK behavior, with hydrophilic fragments correlating with low distribution and high excretion, while lipophilic fragments were linked to enhanced absorption and plasma protein binding. By leveraging consensus predictions from an ensemble of classification models trained on calculated properties and molecular descriptors for each PK parameter, we achieved the most accurate predictions: ternary classifiers excelled in total clearance, whereas binary classifiers performed better for the others. In conclusion, this study provides a solid foundation for PK parameter classification and predictive modeling using a well curated human PK dataset. Integrating comprehensive data curation with ML presents a powerful strategy for accelerating drug design and enhancing rational therapeutic decision making. SIGNIFICANCE STATEMENT: Accurate prediction of absorption, distribution, metabolism, and excretion and pharmacokinetics (PK) properties is crucial for drug discovery and development. Data curation and availability is invaluable to the progress of the development of robust machine learning models to predict drug PK behavior. Standardized labeling and/or classification of human PK parameters would provide better interpretability and predictability. Integrating comprehensive data curation with machine learning presents a powerful strategy for accelerating drug design and enhancing rational therapeutic decision making.

Intranasal diamorphine (IND), approved for managing breakthrough pain in the UK, has been identified as an acceptable alternative offering effective, expedient, and less traumatic analgesia for children. However, the current dose regimen in pediatric populations relies on clinical expertise while the pharmacokinetics properties are poorly understood. This study aimed to develop diamorphine population pharmacokinetics (pop‐PK) models and simulate the IND dosing in virtual pediatric subjects. An integrated four‐compartment pop‐PK model with first‐order absorption and elimination provided an appropriate fit and characterized publicly available 385 concentration measurements of diamorphine, 6‐monoacetylmorphine, and morphine collected from adults. Body weight allometry and renal function maturation (age) were incorporated into the final model, serving as two covariates. The estimated IND relative bioavailability was around 52% compared with intramuscularly injected diamorphine. Using this final model, the morphine plasma concentrations, as the active metabolite for pain relief, were simulated in virtual subjects. The utility of model extrapolation was supported by external verification with acceptable average fold errors of 1.06 ± 0.30 and 0.83 ± 0.07 for morphine maximum concentration and exposures. Meanwhile, the simulated morphine concentration–time profiles could recover the PK profiles observed in children after a single dose of IND. The model‐based dosing simulations were therefore assessed in four children age groups to match the therapeutic window of morphine concentrations in steady state (10–20 μg/L). Our study demonstrates that the dose regimen of 0.3 mg/kg loading dose plus 0.1 mg/kg hourly maintenance dose is generally appropriate for multiple pediatric populations with breakthrough pain, in the view of PK.