Bacteriophages (phages), as natural antibacterial agents, are being rediscovered because of the growing threat of multi- and pan-drug-resistant bacterial pathogens globally. However, with an estimated 10³¹ phages on the planet, finding the right phage to recognize a specific bacterial host is like looking for a needle in a trillion haystacks. The host range of a phage is primarily determined by phage tail fibers (or spikes), which initially mediate reversible and specific recognition and adsorption by susceptible bacteria. Recent significant advances at single-molecule and atomic levels have begun to unravel the structural organization of tail fibers and underlying mechanisms of phage–host interactions. Here, we discuss the molecular mechanisms and models of the tail fibers of the well-characterized T4 phage’s interaction with host surface receptors. Structure–function knowledge of tail fibers will pave the way for reprogramming phage host range and will bring future benefits through more-effective phage therapy in medicine. Furthermore, the design strategies of tail fiber engineering are briefly summarized, including machine-learning-assisted engineering inspired by the increasingly enormous amount of phage genetic information.

Endocrine-disrupting chemicals (EDCs) are notable for their persistence, bioaccumulation, and associations with cancer. Human nuclear receptors (hNRs) are primary targets disrupted by these persistent EDCs, resulting in alterations to xenobiotic metabolism, lipid homeostasis, and endocrine function, which can lead to carcinogenic effects. Despite their hazardous effects, comprehensive studies on EDC interactions and their impacts on hNRs remain limited. Here, we profiled the interactions of persistent EDCs, including PFAS, plastic additives, bisphenols, polybrominated diphenyl ethers, and phthalates, with key hNRs such as PXR, CAR, PPARα, PPARγ, PPARδ, AR, and RORγt. Through controlled molecular docking simulations, we observed strong binding of the EDCs to these receptors. Further analysis showed that EDCs exhibit strong binding activity towards hNRs by preferentially interacting with hydrophobic amino acids, namely leucine, isoleucine, methionine, and phenylalanine. PFAS demonstrated the highest binding affinity, characterized by a combination of complementary hydrophobic interactions from their fluorinated carbon chains and polar interactions from their functional groups (e.g., carboxylate, sulfonate) across all receptors. Distinct polycyclic and hydrophobic trends, contributing to strong NR binding, were evident in non-PFAS and nonplastic EDCs. The hNR activity assay in HepG2 cells revealed agonistic effects of dicyclohexyl phthalate (DCHP) and di-2-ethylhexyl phthalate (DEHP) on most receptors, except for PPARα. The hNR transcription factor pathway assay in HepG2 cells demonstrated increased gene expression of VDRE and PXR, suggesting potential chronic effects on xenobiotic metabolism and calcium homeostasis. Overall, our findings demonstrate the need for further research into the endocrine disruption and carcinogenic effects of these persistent EDCs.

Accurately predicting protein–ligand binding affinity is an important step in the drug discovery process. Deep learning (DL) methods have improved binding affinity prediction by using diverse categories of molecular data. However, many models rely heavily on interaction and sequence data, which impedes proper learning and limits performance in de novo applications. To address these limitations, we developed a novel graph neural network model, called StructureNet (structure-based graph neural network), to predict protein–ligand binding affinity. StructureNet improves existing DL methods by focusing entirely on structural descriptors to mitigate data memorization issues introduced by sequence and interaction data. StructureNet represents the protein and ligand structures as graphs, which are processed using a GNN-based ensemble deep learning model. StructureNet achieved a PCC of 0.68 and an AUC of 0.75 on the PDBBind v.2020 Refined Set, outperforming similar structure-based models. External validation on the DUDE-Z dataset showed that StructureNet can effectively distinguish between active and decoy ligands. Further testing on a small subset of well-known drugs indicates that StructureNet has high potential for rapid virtual screening applications. We also hybridized StructureNet with interaction- and sequence-based models to investigate their impact on testing accuracy and found minimal difference (0.01 PCC) between merged models and StructureNet as a standalone model. An ablation study found that geometric descriptors were the key drivers of model performance, with their removal leading to a PCC decrease of over 15.7%. Lastly, we tested StructureNet on ensembles of binding complex conformers generated using molecular dynamics (MD) simulations and found that incorporating multiple conformations of the same complex often improves model accuracy by capturing binding site flexibility. Overall, the results show that structural data alone are sufficient for binding affinity predictions and can address pattern recognition challenges introduced by sequence and interaction features. Additionally, structural representations of protein–ligand complexes can be considerably improved using geometric and topological descriptors. We made StructureNet GUI interface freely available online.

Per- and polyfluoroalkyl substances (PFAS) are environmentally persistent organofluorines linked to cancer, organ dysfunction, and other health problems. This study used quantitative structure–property relationship (QSPR) and quantitative structure–activity relationship (QSAR) modeling to examine the binding of PFAS to estrogen receptor alpha (ERα) and beta (ERβ). Molecular docking of 14,591 PFAS compounds was performed, and docking scores were used as a measure of receptor affinity. QSPR models were built for two datasets: the ERα and ERβ top binders (TBs), and a set of commonly exposed (CE) PFAS. These models quantified how chemical descriptors influence binding affinity. Across the models, higher density and electrophilicity indicated positive correlations with affinity, while surface tension indicated negative correlations. Electrostatic descriptors, including HOMO energy and positive Fukui index (F⁺ max), were part of the models but showed inconsistent trends. The CE QSPR models displayed correlations that conflicted with those of the TB models. Following QSPR analysis, 66 QSAR models were developed using a mix of top binders and experimental data. These models achieved strong performance, with R² values averaging 0.95 for training sets and 0.78 for test sets, that indicated reliable predictive ability. To improve generalizability, large-set QSAR models were created for each receptor. After outlier removal, these models reached R² values of 0.68–0.71, which supports their use in screening structurally diverse PFAS. Overall, QSPR and QSAR analyses reveal key chemical features that influence PFAS–ER binding. This predictive approach provides a scalable framework to assess the binding interactions of structurally diverse PFAS to ERs and other nuclear receptors. All the codes, data, and the GUI visualization of the results are freely available at sivaGU/QSPR-QSAR-Molecular-Visualization-Tool.

The current epitope selection methods for peptide vaccines often rely on epitope binding affinity predictions, prompting the need for the development of more sophisticated in silico methods to determine immunologically relevant epitopes. Here, we developed AutoPepVax to expedite and improve the in silico epitope selection for peptide vaccine design. AutoPepVax is a novel program that automatically identifies non-toxic and non-allergenic epitopes capable of inducing tumor-infiltrating lymphocytes by considering various epitope characteristics. AutoPepVax employs random forest classification and linear regression machine-learning-based models, which are trained with datasets derived from tumor samples. AutoPepVax, along with documentation on how to run the program, is freely available on GitHub. We used AutoPepVax to design a pan-cancer peptide vaccine targeting epidermal growth factor receptor (EGFR) missense mutations commonly found in lung adenocarcinoma (LUAD), colorectal adenocarcinoma (CRAD), glioblastoma multiforme (GBM), and head and neck squamous cell carcinoma (HNSCC). These mutations have been previously targeted in clinical trials for EGFR-specific peptide vaccines in GBM and LUAD, and they show promise but lack demonstrated clinical efficacy. Using AutoPepVax, our analysis of 96 EGFR mutations identified 368 potential MHC-I-restricted epitope–HLA pairs from 49,113 candidates and 430 potential MHC-II-restricted pairs from 168,669 candidates. Notably, 19 mutations presented viable epitopes for MHC I and II restrictions. To evaluate the potential impact of a pan-cancer vaccine composed of these epitopes, we used our program, PCOptim, to curate a minimal list of epitopes with optimal population coverage. The world population coverage of our list ranged from 81.8% to 98.5% for MHC Class II and Class I epitopes, respectively. From our list of epitopes, we constructed 3D epitope–MHC models for six MHC-I-restricted and four MHC-II-restricted epitopes, demonstrating their epitope binding potential and interaction with T-cell receptors. AutoPepVax’s comprehensive approach to in silico epitope selection addresses vaccine safety, efficacy, and broad applicability. Future studies aim to validate the AutoPepVax-designed vaccines with murine tumor models that harbor the studied mutations.

JEI is a scientific journal for middle and high school scientists