Exploring the computational methods for protein-ligand binding site prediction

Artificial intelligence for RNA-ligand interaction prediction: advances and prospects

Accurate prediction of RNA-ligand interactions is vital for understanding biological processes and advancing RNA-targeted drug discovery. Given their complexity, artificial intelligence (AI) is revolutionizing the study of RNA-ligand interactions, offering insights into the complex dynamics and therapeutic potential of RNA. In this review, we highlight advances in AI-driven RNA-ligand binding site identification, structure modeling, binding mode and binding affinity prediction, and virtual screening (VS). We also discuss key challenges, such as data set scarcity and modeling RNA flexibility. Future directions emphasize integrating cutting-edge AI techniques with physics-based models and expanding experimental data sets to enhance RNA-ligand interaction predictions.

MVSF-AB: accurate antibody-antigen binding affinity prediction via multi-view sequence feature learning

MOTIVATION

Predicting the binding affinity between antigens and antibodies accurately is crucial for assessing therapeutic antibody effectiveness and enhancing antibody engineering and vaccine design. Traditional machine learning methods have been widely used for this purpose, relying on interfacial amino acids' structural information. Nevertheless, due to technological limitations and high costs of acquiring structural data, the structures of most antigens and antibodies are unknown, and sequence-based methods have gained attention. Existing sequence-based approaches designed for protein-protein affinity prediction exhibit a significant drop in performance when applied directly to antibody-antigen affinity prediction due to imbalanced training data and lacking design in the model framework specifically for antibody-antigen, hindering the learning of key features of antibodies and antigens. Therefore, we propose MVSF-AB, a Multi-View Sequence Feature learning for accurate Antibody-antigen Binding affinity prediction.

RESULTS

MVSF-AB designs a multi-view method that fuses semantic features and residue features to fully utilize the sequence information of antibody-antigen and predicts the binding affinity. Experimental results demonstrate that MVSF-AB outperforms existing approaches in predicting unobserved natural antibody-antigen affinity and maintains its effectiveness when faced with mutant strains of antibodies.

AVAILABILITY AND IMPLEMENTATION

Datasets we used and source code are available on our public GitHub repository https://github.com/TAI-Medical-Lab/MVSF-AB.

Structural Bioinformatics Applied to Acetylcholinesterase Enzyme Inhibition

Acetylcholinesterase (AChE) is a critical enzyme involved in neurotransmission by hydrolyzing acetylcholine at the synaptic cleft, making it a key target for drug discovery, particularly in the treatment of neurodegenerative disorders such as Alzheimer's disease. Computational approaches, particularly molecular docking and molecular dynamics (MD) simulations, have become indispensable tools for identifying and optimizing AChE inhibitors by predicting ligand-binding affinities, interaction mechanisms, and conformational dynamics. This review serves as a comprehensive guide for future research on AChE using molecular docking and MD simulations. It compiles and analyzes studies conducted over the past five years, providing a critical evaluation of the most widely used computational tools, including AutoDock, AutoDock Vina, and GROMACS, which have significantly contributed to the advancement of AChE inhibitor screening. Furthermore, we identify PDB ID: 4EY7, the most frequently used AChE crystal structure in docking studies, and highlight Donepezil, a well-established reference molecule widely employed as a control in computational screening for novel inhibitors. By examining these key aspects, this review aims to enhance the accuracy and reliability of virtual screening approaches and guide researchers in selecting the most appropriate computational methodologies. The integration of docking and MD simulations not only improves hit identification and lead optimization but also provides deeper mechanistic insights into AChE-ligand interactions, contributing to the rational design of more effective AChE inhibitors.

SELFprot: Effective and Efficient Multitask Finetuning Methods for Protein Parameter Prediction

Accurately predicting protein-ligand interactions and enzymatic kinetics remains a challenge for computational biology. Here, we present SELFprot, a suite of modular transformer-based machine learning architectures that leverage the ESM2-35M model architecture for protein sequence and small molecule embeddings to improve predictions of complex biochemical interactions. SELFprot employs multitask learning and parameter-efficient finetuning through low-rank adaptation, allowing for adaptive, data-driven model refinement. Furthermore, ensemble learning techniques are used to enhance the robustness and reduce the prediction variance. Evaluated on the BindingDB and CatPred-DB data sets, SELFprot achieves competitive performance with notable improvements in parameter-efficient prediction of kcat, Km, Ki, Kd, IC50, and EC50 values as well as the classification of functional site residues. With comparable accuracy to existing models and an order of magnitude fewer parameters, SELFprot demonstrates versatility and efficiency, making it a valuable tool for protein-ligand interaction studies in bioengineering.

Computer aided aptamer selection for fabrication of electrochemical sensor to detect Aflatoxin B1

Aflatoxin B1 (AFB1) is a naturally occurring toxin produced by Aspergillus flavus and Aspergillus parasiticus. The AFB1 is classified as a potent carcinogen and poses significant health risks both to humans and animals. Early detection of the toxin in post-harvest agricultural products will save lives and promote healthy food production. In this study, stratified docking approach was utilized to screen and identify potential aptamers that can bind to AFB1. ssDNA sequences were acquired from the Mendeley dataset, secondary and tertiary structures were predicted through a series of bioinformatics pipelines. Further, the final DNA tertiary structures were minimized and SiteMap algorithm was used to probe and locate binding cavities. According to the final XP docking result, a 34 nt sequence (5'-ATCCTGTGAGGAATGCTCATGCATAGCAAGGGCT-3') aptamer with a docking score of -5.959 kcal/mol was considered for 200 ns MD Simulation. Finally, the screened DNA-aptamer was immobilized over the gold surface based on Au-S chemistry and utilized for the detection of AFB1. The fabricated DNA-aptamer electrode demonstrated a good analytical performance including wide linear range (1.0 to 1000 ng L-1), detection limit (1.0 ng L-1), high stability, and reproducibility.

Screening of phytochemicals from Clerodendrum inerme (L.) Gaertn as potential anti-breast cancer compounds targeting EGFR: an in-silico approach

Breast cancer (BC) is the most prevalent malignancy among women around the world. The epidermal growth factor receptor (EGFR) is a tyrosine kinase receptor (RTK) of the ErbB/HER family. It is essential for triggering the cellular signaling cascades that control cell growth and survival. However, perturbations in EGFR signaling lead to cancer development and progression. Hence, EGFR is regarded as a prominent therapeutic target for breast cancer. Therefore, in the current investigation, EGFR was targeted with phytochemicals from Clerodendrum inerme (L.) Gaertn (C. inerme). A total of 121 phytochemicals identified by gas chromatography-mass spectrometry (GC-MS) analysis were screened against EGFR through molecular docking, ADMET analyses (Absorption, Distribution, Metabolism, Excretion, and Toxicity), PASS predictions, and molecular dynamics simulation, which revealed three potential hit compounds with CIDs 10586 [i.e. alpha-bisabolol (-6.4 kcal/mol)], 550281 [i.e. 2,(4,4-Trimethyl-3-hydroxymethyl-5a-(3-methyl-but-2-enyl)-cyclohexene) (-6.5 kcal/mol)], and 161271 [i.e. salvigenin (-7.4 kcal/mol)]. The FDA-approved drug gefitinib was used to compare the inhibitory effects of the phytochemicals. The top selected compounds exhibited good ADMET properties and obeyed Lipinski's rule of five (ROF). The molecular docking analysis showed that salvigenin was the best among the three compounds and formed bonds with the key residue Met 793. Furthermore, the molecular mechanics generalized born surface area (MMGBSA) calculations, molecular dynamics simulation, and normal mode analysis validated the binding affinity of the compounds and also revealed the strong stability and compactness of phytochemicals at the docked site. Additionally, DFT and DOS analyses were done to study the reactivity of the compounds and to further validate the selected phytochemicals. These results suggest that the identified phytochemicals possess high inhibitory potential against the target EGFR and can treat breast cancer. However, further in vitro and in vivo investigations are warranted towards the development of these constituents into novel anti-cancer drugs.

Machine learning approaches for predicting protein-ligand binding sites from sequence data

Proteins, composed of amino acids, are crucial for a wide range of biological functions. Proteins have various interaction sites, one of which is the protein-ligand binding site, essential for molecular interactions and biochemical reactions. These sites enable proteins to bind with other molecules, facilitating key biological functions. Accurate prediction of these binding sites is pivotal in computational drug discovery, helping to identify therapeutic targets and facilitate treatment development. Machine learning has made significant contributions to this field by improving the prediction of protein-ligand interactions. This paper reviews studies that use machine learning to predict protein-ligand binding sites from sequence data, focusing on recent advancements. The review examines various embedding methods and machine learning architectures, addressing current challenges and the ongoing debates in the field. Additionally, research gaps in the existing literature are highlighted, and potential future directions for advancing the field are discussed. This study provides a thorough overview of sequence-based approaches for predicting protein-ligand binding sites, offering insights into the current state of research and future possibilities.

Prediction of Ca2+ Binding Site in Proteins With a Fast and Accurate Method Based on Statistical Mechanics and Analysis of Crystal Structures

Predicting the precise locations of metal binding sites within metalloproteins is a crucial challenge in biophysics. A fast, accurate, and interpretable computational prediction method can complement the experimental studies. In the current work, we have developed a method to predict the location of Ca2+ ions in calcium-binding proteins using a physics-based method with an all-atom description of the proteins, which is substantially faster than the molecular dynamics simulation-based methods with accuracy as good as data-driven approaches. Our methodology uses the three-dimensional reference interaction site model (3D-RISM), a statistical mechanical theory, to calculate Ca2+ ion density around protein structures, and the locations of the Ca2+ ions are obtained from the density. We have taken previously used datasets to assess the efficacy of our method as compared to previous works. Our accuracy is 88%, comparable with the FEATURE program, one of the well-known data-driven methods. Moreover, our method is physical, and the reasons for failures can be ascertained in most cases. We have thoroughly examined the failed cases using different structural and crystallographic measures, such as B-factor, R-factor, electron density map, and geometry at the binding site. It has been found that x-ray structures have issues in many of the failed cases, such as geometric irregularities and dubious assignment of ion positions. Our algorithm, along with the checks for structural accuracy, is a major step in predicting calcium ion positions in metalloproteins.

Quantitative Prediction of Protein-Polyelectrolyte Binding Thermodynamics: Adsorption of Heparin-Analog Polysulfates to the SARS-CoV-2 Spike Protein RBD

Interactions of polyelectrolytes (PEs) with proteins play a crucial role in numerous biological processes, such as the internalization of virus particles into host cells. Although docking, machine learning methods, and molecular dynamics (MD) simulations are utilized to estimate binding poses and binding free energies of small-molecule drugs to proteins, quantitative prediction of the binding thermodynamics of PE-based drugs presents a significant obstacle in computer-aided drug design. This is due to the sluggish dynamics of PEs caused by their size and strong charge-charge correlations. In this paper, we introduce advanced sampling methods based on a force-spectroscopy setup and theoretical modeling to overcome this barrier. We exemplify our method with explicit solvent all-atom MD simulations of the interactions between anionic PEs that show antiviral properties, namely heparin and linear polyglycerol sulfate (LPGS), and the SARS-CoV-2 spike protein receptor binding domain (RBD). Our prediction for the binding free-energy of LPGS to the wild-type RBD matches experimentally measured dissociation constants within thermal energy, k B T, and correctly reproduces the experimental PE-length dependence. We find that LPGS binds to the Delta-variant RBD with an additional free-energy gain of 2.4 k B T, compared to the wild-type RBD, due to the additional presence of two mutated cationic residues contributing to the electrostatic energy gain. We show that the LPGS-RBD binding is solvent dominated and enthalpy driven, though with a large entropy-enthalpy compensation. Our method is applicable to general polymer adsorption phenomena and predicts precise binding free energies and reconfigurational friction as needed for drug and drug-delivery design.

CryptoBench: cryptic protein-ligand binding sites dataset and benchmark

MOTIVATION

Structure-based methods for detecting protein-ligand binding sites play a crucial role in various domains, from fundamental research to biomedical applications. However, current prediction methodologies often rely on holo (ligand-bound) protein conformations for training and evaluation, overlooking the significance of the apo (ligand-free) states. This oversight is particularly problematic in the case of cryptic binding sites (CBSs) where holo-based assessment yields unrealistic performance expectations.

RESULTS

To advance the development in this domain, we introduce CryptoBench, a benchmark dataset tailored for training and evaluating novel CBS prediction methodologies. CryptoBench is constructed upon a large collection of apo-holo protein pairs, grouped by UniProtID, clustered by sequence identity, and filtered to contain only structures with substantial structural change in the binding site. CryptoBench comprises 1107 structures with predefined cross-validation splits, making it the most extensive CBS dataset to date. To establish a performance baseline, we measured the predictive power of sequence- and structure-based CBS residue prediction methods using the benchmark. We selected PocketMiner as the state-of-the-art representative of the structure-based methods for CBS detection, and P2Rank, a widely-used structure-based method for general binding site prediction that is not specifically tailored for cryptic sites. For sequence-based approaches, we trained a neural network to classify binding residues using protein language model embeddings. Our sequence-based approach outperformed PocketMiner and P2Rank across key metrics, including area under the curve, area under the precision-recall curve, Matthew's correlation coefficient, and F1 scores. These results provide baseline benchmark results for future CBS and potentially also non-CBS prediction endeavors, leveraging CryptoBench as the foundational platform for further advancements in the field.

AVAILABILITY AND IMPLEMENTATION

The CryptoBench dataset, including the benchmark model, is available on Open Science Framework-https://osf.io/pz4a9/. The code and tutorial are available at the GitHub repository-https://github.com/skrhakv/CryptoBench/.

Neuroprotective Potential of Indole-Based Compounds: A Biochemical Study on Antioxidant Properties and Amyloid Disaggregation in Neuroblastoma Cells

Based on the established neuroprotective properties of indole-based compounds and their significant potential as multi-targeted therapeutic agents, a series of synthetic indole-phenolic compounds was evaluated as multifunctional neuroprotectors. Each compound demonstrated metal-chelating properties, particularly in sequestering copper ions, with quantitative analysis revealing approximately 40% chelating activity across all the compounds. In cellular models, these hybrid compounds exhibited strong antioxidant and cytoprotective effects, countering reactive oxygen species (ROS) generated by the Aβ(25-35) peptide and its oxidative byproduct, hydrogen peroxide, as demonstrated by quantitative analysis showing on average a 25% increase in cell viability and a reduction in ROS levels to basal states. Further analysis using thioflavin T fluorescence assays, circular dichroism, and computational studies indicated that the synthesized derivatives effectively promoted the self-disaggregation of the Aβ(25-35) fragment. Taken together, these findings suggest a unique profile of neuroprotective actions for indole-phenolic derivatives, combining chelating, antioxidant, and anti-aggregation properties, which position them as promising compounds for the development of multifunctional agents in Alzheimer's disease therapy. The methods used provide reliable in vitro data, although further in vivo validation and assessment of blood-brain barrier penetration are needed to confirm therapeutic efficacy and safety.

In Silico Tools to Score and Predict Cholesterol-Protein Interactions

Cholesterol is structurally distinct from other lipids, which confers it with singular roles in membrane organization and protein function. As a signaling molecule, cholesterol engages in discrete interactions with transmembrane, peripheral, and certain soluble proteins to control cellular responses. Accordingly, the cholesterol-protein interface is central to cholesterol-related diseases and is an essential consideration in drug design. However, cholesterol's hydrophobic, un-drug-like nature presents a unique challenge to traditional in silico analyses. In this Perspective, we survey a collection of tools designed to predict and evaluate cholesterol binding sites in proteins, including classical sequence motifs, molecular docking, template-based strategies, molecular dynamics simulations, and recent artificial intelligence approaches. We then comment on contemporary tools to evaluate ligand-protein interactions, their applicability to cholesterol, and the yet-untapped potential of cholesterol-protein interactions in human health and disease.

Structure-based protein and small molecule generation using EGNN and diffusion models: A comprehensive review

Recent breakthroughs in deep learning have revolutionized protein sequence and structure prediction. These advancements are built on decades of protein design efforts, and are overcoming traditional time and cost limitations. Diffusion models, at the forefront of these innovations, significantly enhance design efficiency by automating knowledge acquisition. In the field of de novo protein design, the goal is to create entirely novel proteins with predetermined structures. Given the arbitrary positions of proteins in 3-D space, graph representations and their properties are widely used in protein generation studies. A critical requirement in protein modelling is maintaining spatial relationships under transformations (rotations, translations, and reflections). This property, known as equivariance, ensures that predicted protein characteristics adapt seamlessly to changes in orientation or position. Equivariant graph neural networks offer a solution to this challenge. By incorporating equivariant graph neural networks to learn the score of the probability density function in diffusion models, one can generate proteins with robust 3-D structural representations. This review examines the latest deep learning advancements, specifically focusing on frameworks that combine diffusion models with equivariant graph neural networks for protein generation.

Databases of ligand-binding pockets and protein-ligand interactions

Many research groups and institutions have created a variety of databases curating experimental and predicted data related to protein-ligand binding. The landscape of available databases is dynamic, with new databases emerging and established databases becoming defunct. Here, we review the current state of databases that contain binding pockets and protein-ligand binding interactions. We have compiled a list of such databases, fifty-three of which are currently available for use. We discuss variation in how binding pockets are defined and summarize pocket-finding methods. We organize the fifty-three databases into subgroups based on goals and contents, and describe standard use cases. We also illustrate that pockets within the same protein are characterized differently across different databases. Finally, we assess critical issues of sustainability, accessibility and redundancy.

Revealing Missing Protein-Ligand Interactions Using AlphaFold Predictions

Protein-ligand interactions represent an essential step to understand the bases of molecular recognition, an intense field of research in many scientific areas. Structural biology has played a central role in unveiling protein-ligand interactions, but current techniques are still not able to reliably describe the interactions of ligands with highly flexible regions. In this work, we explored the capacity of AlphaFold2 (AF2) to estimate the presence of interactions between ligands and residues belonging to disordered regions. As these interactions are missing in the crystallographic-derived structures, we called them "ghost interactions". Using a set of protein structures experimentally obtained after AF2 was trained, we found that the obtained models are good predictors of regions associated with order-disorder transitions. Additionally, we found that AF2 predicts residues making ghost interactions with ligands, which are mostly buried and show differential evolutionary conservation with the rest of the residues located in the flexible region. Our findings could fuel current areas of research that consider, given their biological relevance and their involvement in diseases, intrinsically disordered proteins as potentially valuable targets for drug development.

Indian Medicinal Plant-Derived Phytochemicals as Potential Antidotes for Snakebite: A Pharmacoinformatic Study of Atrolysin Inhibitors

Snakebite envenoming is a significant health threat, particularly in tropical regions, causing substantial morbidity and mortality. Traditional treatments, including antivenom therapy, have limitations and associated risks. This research aims to discover novel phytochemical antidotes for snakebites, specifically targeting the western diamondback rattlesnake (Crotalus atrox) venom metalloproteinase Atrolysin. Utilizing pharmacoinformatic techniques such as molecular docking, high-throughput ligand screening, pharmacophore mapping, pharmacokinetic profiling, and molecular dynamics (MD) simulations, we analyzed phytochemicals from the Indian Medicinal Plants, Phytochemistry And Therapeutics (IMPPAT) database alongside well-known nine metalloproteinase inhibitors from the PubChem database. From an initial set of 17,967 compounds, 4708 unique compounds were identified for further study. These compounds were evaluated based on drug likeness, molecular descriptors, ADME properties, and toxicity profiles. Binding site predictions and molecular docking identified key interacting residues and binding energies, highlighting several promising compounds. Density functional theory (DFT) analysis provided insights into these compounds' electronic properties and stability. MD simulations assessed the dynamic stability of protein-ligand complexes using parameters such as RMSD, RMSF, the radius of gyration, and hydrogen bond interactions. This study identified top candidates, including CID5291, IMPHY001495, IMPHY014737, IMPHY008983, IMPHY008176, and IMPHY003833, based on their favorable binding energies, interaction forces, and structural stability. These findings suggest that the selected phytochemicals have the potential to serve as effective alternatives to traditional antivenom treatments, offering a promising avenue for further research and development in snakebite management.

Thermal decomposition and atmospheric pressure chemical ionization of alanine using ion mobility spectrometry and computational study

This study investigates the impact of thermal decomposition on the ion mobility spectrum of L-alanine using ion mobility spectrometry (IMS) and computational methods. By employing a post-injection delay system, we examined the evolution of ion peaks corresponding to thermal decomposition products and their interaction with protonated alanine. Experimental results revealed that the observed ion mobility spectra predominantly feature protonated isomers and adduct ions. Computational analysis using Density Functional Theory (DFT) predicted the thermodynamically favored structures and stabilities of these products. Findings indicate that protonation at the nitrogen site in alanine is more stable than at the oxygen site, and observed peaks correspond to protonated isomers and adducts formed with ammonium ions. Further investigations showed that thermal decomposition of alanine generates ammonia, contributing to the formation of new adduct ions. This research provides new insights into the behavior of amino acids under thermal conditions with implications for analytical chemistry and biochemistry.

Combined Physics- and Machine-Learning-Based Method to Identify Druggable Binding Sites Using SILCS-Hotspots

Identifying druggable binding sites on proteins is an important and challenging problem, particularly for cryptic, allosteric binding sites that may not be obvious from X-ray, cryo-EM, or predicted structures. The Site-Identification by Ligand Competitive Saturation (SILCS) method accounts for the flexibility of the target protein using all-atom molecular simulations that include various small molecule solutes in aqueous solution. During the simulations, the combination of protein flexibility and comprehensive sampling of the water and solute spatial distributions can identify buried binding pockets absent in experimentally determined structures. Previously, we reported a method for leveraging the information in the SILCS sampling to identify binding sites (termed Hotspots) of small mono- or bicyclic compounds, a subset of which coincide with known binding sites of drug-like molecules. Here, we build on that physics-based approach and present a ML model for ranking the Hotspots according to the likelihood they can accommodate drug-like molecules (e.g., molecular weight >200 Da). In the independent validation set, which includes various enzymes and receptors, our model recalls 67% and 89% of experimentally validated ligand binding sites in the top 10 and 20 ranked Hotspots, respectively. Furthermore, we show that the model's output Decision Function is a useful metric to predict binding sites and their potential druggability in new targets. Given the utility the SILCS method for ligand discovery and optimization, the tools presented represent an important advancement in the identification of orthosteric and allosteric binding sites and the discovery of drug-like molecules targeting those sites.

Turbocharging protein binding site prediction with geometric attention, inter-resolution transfer learning, and homology-based augmentation

BACKGROUND

Locating small molecule binding sites in target proteins, in the resolution of either pocket or residue, is critical in many drug-discovery scenarios. Since it is not always easy to find such binding sites using conventional methods, different deep learning methods to predict binding sites out of protein structures have been developed in recent years. The existing deep learning based methods have several limitations, including (1) the inefficiency of the CNN-only architecture, (2) loss of information due to excessive post-processing, and (3) the under-utilization of available data sources.

METHODS

We present a new model architecture and training method that resolves the aforementioned problems. First, by layering geometric self-attention units on top of residue-level 3D CNN outputs, our model overcomes the problems of CNN-only architectures. Second, by configuring the fundamental units of computation as residues and pockets instead of voxels, our method reduced the information loss from post-processing. Lastly, by employing inter-resolution transfer learning and homology-based augmentation, our method maximizes the utilization of available data sources to a significant extent.

RESULTS

The proposed method significantly outperformed all state-of-the-art baselines regarding both resolutions-pocket and residue. An ablation study demonstrated the indispensability of our proposed architecture, as well as transfer learning and homology-based augmentation, for achieving optimal performance. We further scrutinized our model's performance through a case study involving human serum albumin, which demonstrated our model's superior capability in identifying multiple binding sites of the protein, outperforming the existing methods.

CONCLUSIONS

We believe that our contribution to the literature is twofold. Firstly, we introduce a novel computational method for binding site prediction with practical applications, substantiated by its strong performance across diverse benchmarks and case studies. Secondly, the innovative aspects in our method- specifically, the design of the model architecture, inter-resolution transfer learning, and homology-based augmentation-would serve as useful components for future work.

Plant-based seafood alternatives: Current insights on the nutrition, protein-flavour interactions, and the processing of these foods

Fish are an important food source; however, the sustainability of current seafood supplies is a major concern for key stakeholders. The development of plant-based seafood alternatives may be suitable products to alleviate some of the pressures on aquatic ecosystems and help support environmental sustainability. However, the wide-spread adoption of these products weighs heavily on the ingredients used in the formulations which should not only satisfy nutritional and sustainability targets but must also meet consumer approval and functionality. In this review, we highlight recent advances in our understanding of the nutritional quality and sensory challenges in particular flavour (which includes taste and aroma), that have so far proven difficult to overcome in the development of plant-based seafood alternatives. Protein interactions that contribute to flavour development in plant-based seafood alternatives and the factors that impact these interactions are also discussed. We also review the recent advances in the innovative technologies used to improve the texture of products in this emerging food category. Finally, we highlight key areas for targeted research to advance the development of this growing segment of food products.

Precise redesign for improving enzyme robustness based on coevolutionary analysis and multidimensional virtual screening

Natural enzymes are able to function effectively under optimal physiological conditions, but the intrinsic performance often fails to meet the demands of industrial production. Existing strategies are based mainly on the evaluation and subsequent combination of single-point mutations; however, this approach often suffers from a limited number of designable residues and from low accuracy. Here, we propose a strategy (Co-MdVS) based on coevolutionary analysis and multidimensional virtual screening for precise design to improve enzyme robustness, employing nattokinase as a model. Using this strategy, we efficiently screened 8 dual mutants with enhanced thermostability from a virtual mutation library containing 7980 mutants. After further iterative combination, the optimal mutant M6 exhibited a 31-fold increase in half-life at 55 °C, significantly enhanced acid resistance, and improved catalytic efficiency with different substrates. Molecular dynamics simulations indicated that the reduced flexibility of thermal and acid-sensitive regions resulted in a significantly increased robustness of M6. Furthermore, the potential of multidimensional virtual screening in enhancing design precision has been validated on l-rhamnose isomerase and PETase. Therefore, the Co-MdVS strategy introduced in this research may offer a viable approach for developing enzymes with enhanced robustness.

Elucidation of the CadA Protein 3D Structure and Affinity for Metals

The mitigation of cadmium (Cd) pollution, a significant ecological threat, is of paramount importance. Pseudomonas aeruginosa harbors 2 Cd resistance genes, namely, cadR and cadA. Presently, our focus is on the identification and characterization of the cation-transporting P-type ATPase (cadA) in Pseudomonas aeruginosa BC15 through in silico methods. The CadA protein and its binding capacities remain poorly understood, with no available structural elucidation. The presence of the cadA gene in P aeruginosa was confirmed, showing a striking 99% sequence similarity with both P aeruginosa and P putida. Phylogenetic analysis unveiled the evolutionary relationship between CadA protein sequences from various Pseudomonas species. Physicochemical analysis demonstrated the stability of CadA, revealing a composition of 690 amino acids, a molecular weight of 73 352.85, and a predicted isoelectric point (PI) of 5.39. Swiss-Model homology modelling unveiled a 33.73% sequence homology with CopA (3J09), and the projected structure indicated that 89.3% of amino acid residues were situated favourably within the Ramachandran plot, signifying energetic stability. Notably, the study identified metal-binding sites in CadA, namely, H3, C30, C32, C35, H48, C89, and C106. Docking studies revealed a higher efficiency of Cd binding with CadA compared to other heavy metals. This underscores the crucial role of N-terminal cysteine residues in Cd removal. It is evident that CadA of P aeruginosa BC15 plays a crucial role in Cd tolerance, rendering it a potential microorganism for Cd toxicity bioremediation. The structural and functional elucidation of CadA, facilitated by this study, holds promise for advancing cost-effective strategies in the remediation of cadmium-contaminated environments.

Protein ligand binding site prediction using graph transformer neural network

Ligand binding site prediction is a crucial initial step in structure-based drug discovery. Although several methods have been proposed previously, including those using geometry based and machine learning techniques, their accuracy is considered to be still insufficient. In this study, we introduce an approach that leverages a graph transformer neural network to rank the results of a geometry-based pocket detection method. We also created a larger training dataset compared to the conventionally used sc-PDB and investigated the correlation between the dataset size and prediction performance. Our findings indicate that utilizing a graph transformer-based method alongside a larger training dataset could enhance the performance of ligand binding site prediction.

LVPocket: integrated 3D global-local information to protein binding pockets prediction with transfer learning of protein structure classification

BACKGROUND

Previous deep learning methods for predicting protein binding pockets mainly employed 3D convolution, yet an abundance of convolution operations may lead the model to excessively prioritize local information, thus overlooking global information. Moreover, it is essential for us to account for the influence of diverse protein folding structural classes. Because proteins classified differently structurally exhibit varying biological functions, whereas those within the same structural class share similar functional attributes.

RESULTS

We proposed LVPocket, a novel method that synergistically captures both local and global information of protein structure through the integration of Transformer encoders, which help the model achieve better performance in binding pockets prediction. And then we tailored prediction models for data of four distinct structural classes of proteins using the transfer learning. The four fine-tuned models were trained on the baseline LVPocket model which was trained on the sc-PDB dataset. LVPocket exhibits superior performance on three independent datasets compared to current state-of-the-art methods. Additionally, the fine-tuned model outperforms the baseline model in terms of performance.

SCIENTIFIC CONTRIBUTION

We present a novel model structure for predicting protein binding pockets that provides a solution for relying on extensive convolutional computation while neglecting global information about protein structures. Furthermore, we tackle the impact of different protein folding structures on binding pocket prediction tasks through the application of transfer learning methods.

CpG Island Definition and Methylation Mapping of the T2T-YAO Genome

Precisely defining and mapping all cytosine (C) positions and their clusters, known as CpG islands (CGIs), as well as their methylation status, are pivotal for genome-wide epigenetic studies, especially when population-centric reference genomes are ready for timely application. Here, we first align the two high-quality reference genomes, T2T-YAO and T2T-CHM13, from different ethnic backgrounds in a base-by-base fashion and compute their genome-wide density-defined and position-defined CGIs. Second, by mapping some representative genome-wide methylation data from selected organs onto the two genomes, we find that there are about 4.7%-5.8% sequence divergency of variable categories depending on quality cutoffs. Genes among the divergent sequences are mostly associated with neurological functions. Moreover, CGIs associated with the divergent sequences are significantly different with respect to CpG density and observed CpG/expected CpG (O/E) ratio between the two genomes. Finally, we find that the T2T-YAO genome not only has a greater CpG coverage than that of the T2T-CHM13 genome when whole-genome bisulfite sequencing (WGBS) data from the European and American populations are mapped to each reference, but also shows more hyper-methylated CpG sites as compared to the T2T-CHM13 genome. Our study suggests that future genome-wide epigenetic studies of the Chinese populations rely on both acquisition of high-quality methylation data and subsequent precision CGI mapping based on the Chinese T2T reference.

Highly accurate carbohydrate-binding site prediction with DeepGlycanSite

As the most abundant organic substances in nature, carbohydrates are essential for life. Understanding how carbohydrates regulate proteins in the physiological and pathological processes presents opportunities to address crucial biological problems and develop new therapeutics. However, the diversity and complexity of carbohydrates pose a challenge in experimentally identifying the sites where carbohydrates bind to and act on proteins. Here, we introduce a deep learning model, DeepGlycanSite, capable of accurately predicting carbohydrate-binding sites on a given protein structure. Incorporating geometric and evolutionary features of proteins into a deep equivariant graph neural network with the transformer architecture, DeepGlycanSite remarkably outperforms previous state-of-the-art methods and effectively predicts binding sites for diverse carbohydrates. Integrating with a mutagenesis study, DeepGlycanSite reveals the guanosine-5'-diphosphate-sugar-recognition site of an important G-protein coupled receptor. These findings demonstrate DeepGlycanSite is invaluable for carbohydrate-binding site prediction and could provide insights into molecular mechanisms underlying carbohydrate-regulation of therapeutically important proteins.

OLB-AC: toward optimizing ligand bioactivities through deep graph learning and activity cliffs

MOTIVATION

Deep graph learning (DGL) has been widely employed in the realm of ligand-based virtual screening. Within this field, a key hurdle is the existence of activity cliffs (ACs), where minor chemical alterations can lead to significant changes in bioactivity. In response, several DGL models have been developed to enhance ligand bioactivity prediction in the presence of ACs. Yet, there remains a largely unexplored opportunity within ACs for optimizing ligand bioactivity, making it an area ripe for further investigation.

RESULTS

We present a novel approach to simultaneously predict and optimize ligand bioactivities through DGL and ACs (OLB-AC). OLB-AC possesses the capability to optimize ligand molecules located near ACs, providing a direct reference for optimizing ligand bioactivities with the matching of original ligands. To accomplish this, a novel attentive graph reconstruction neural network and ligand optimization scheme are proposed. Attentive graph reconstruction neural network reconstructs original ligands and optimizes them through adversarial representations derived from their bioactivity prediction process. Experimental results on nine drug targets reveal that out of the 667 molecules generated through OLB-AC optimization on datasets comprising 974 low-activity, noninhibitor, or highly toxic ligands, 49 are recognized as known highly active, inhibitor, or nontoxic ligands beyond the datasets' scope. The 27 out of 49 matched molecular pairs generated by OLB-AC reveal novel transformations not present in their training sets. The adversarial representations employed for ligand optimization originate from the gradients of bioactivity predictions. Therefore, we also assess OLB-AC's prediction accuracy across 33 different bioactivity datasets. Results show that OLB-AC achieves the best Pearson correlation coefficient (r2) on 27/33 datasets, with an average improvement of 7.2%-22.9% against the state-of-the-art bioactivity prediction methods.

AVAILABILITY AND IMPLEMENTATION

The code and dataset developed in this work are available at github.com/Yueming-Yin/OLB-AC.

DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors for a wide range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed "DrugMap," an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NF-κB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NF-κB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription-factor activity.

Learnt representations of proteins can be used for accurate prediction of small molecule binding sites on experimentally determined and predicted protein structures

Protein-ligand binding site prediction is a useful tool for understanding the functional behaviour and potential drug-target interactions of a novel protein of interest. However, most binding site prediction methods are tested by providing crystallised ligand-bound (holo) structures as input. This testing regime is insufficient to understand the performance on novel protein targets where experimental structures are not available. An alternative option is to provide computationally predicted protein structures, but this is not commonly tested. However, due to the training data used, computationally-predicted protein structures tend to be extremely accurate, and are often biased toward a holo conformation. In this study we describe and benchmark IF-SitePred, a protein-ligand binding site prediction method which is based on the labelling of ESM-IF1 protein language model embeddings combined with point cloud annotation and clustering. We show that not only is IF-SitePred competitive with state-of-the-art methods when predicting binding sites on experimental structures, but it performs better on proxies for novel proteins where low accuracy has been simulated by molecular dynamics. Finally, IF-SitePred outperforms other methods if ensembles of predicted protein structures are generated.

Antiviral Drug Target Identification and Ligand Discovery

This chapter intends to provide a general overview of web-based resources available for antiviral drug discovery studies. First, we explain how the structure for a potential viral protein target can be obtained and then highlight some of the main considerations in preparing for the application of receptor-based molecular docking techniques. Thereafter, we discuss the resources to search for potential drug candidates (ligands) against this target protein receptor, how to screen them, and preparing their analogue library. We make specific reference to free, online, open-source tools and resources which can be applied for antiviral drug discovery studies.

Computational approaches for modeling and structural design of biological systems: A comprehensive review

The convergence of biology and computational science has ushered in a revolutionary era, revolutionizing our understanding of biological systems and providing novel solutions to global problems. The field of genetic engineering has facilitated the manipulation of genetic codes, thus providing opportunities for the advancement of innovative disease therapies and environmental enhancements. The emergence of bio-molecular simulation represents a significant advancement in this particular field, as it offers the ability to gain microscopic insights into molecular-level biological processes over extended periods. Biomolecular simulation plays a crucial role in advancing our comprehension of organismal mechanisms by establishing connections between molecular structures, interactions, and biological functions. The field of computational biology has demonstrated its significance in deciphering intricate biological enigmas through the utilization of mathematical models and algorithms. The process of decoding the human genome has resulted in the advancement of therapies for a wide range of genetic disorders, while the simulation of biological systems contributes to the identification of novel pharmaceutical compounds. The potential of biomolecular simulation and computational biology is vast and limitless. As the exploration of the underlying principles that govern living organisms progresses, the potential impact of this understanding on cancer treatment, environmental restoration, and other domains is anticipated to be transformative. This review examines the notable advancements achieved in the field of computational biology, emphasizing its potential to revolutionize the comprehension and enhancement of biological systems.

zPoseScore model for accurate and robust protein-ligand docking pose scoring in CASP15

We introduce a deep learning-based ligand pose scoring model called zPoseScore for predicting protein-ligand complexes in the 15th Critical Assessment of Protein Structure Prediction (CASP15). Our contributions are threefold: first, we generate six training and evaluation data sets by employing advanced data augmentation and sampling methods. Second, we redesign the "zFormer" module, inspired by AlphaFold2's Evoformer, to efficiently describe protein-ligand interactions. This module enables the extraction of protein-ligand paired features that lead to accurate predictions. Finally, we develop the zPoseScore framework with zFormer for scoring and ranking ligand poses, allowing for atomic-level protein-ligand feature encoding and fusion to output refined ligand poses and ligand per-atom deviations. Our results demonstrate excellent performance on various testing data sets, achieving Pearson's correlation R  = 0.783 and 0.659 for ranking docking decoys generated based on experimental and predicted protein structures of CASF-2016 protein-ligand complexes. Additionally, we obtain an averaged local distance difference test (lDDT pli = 0.558) of AIchemy LIG2 in CASP15 for de novo protein-ligand complex structure predictions. Detailed analysis shows that accurate ligand binding site prediction and side-chain orientation are crucial for achieving better prediction performance. Our proposed model is one of the most accurate protein-ligand pose prediction models and could serve as a valuable tool in small molecule drug discovery.

Enzyme function and evolution through the lens of bioinformatics

Enzymes have been shaped by evolution over billions of years to catalyse the chemical reactions that support life on earth. Dispersed in the literature, or organised in online databases, knowledge about enzymes can be structured in distinct dimensions, either related to their quality as biological macromolecules, such as their sequence and structure, or related to their chemical functions, such as the catalytic site, kinetics, mechanism, and overall reaction. The evolution of enzymes can only be understood when each of these dimensions is considered. In addition, many of the properties of enzymes only make sense in the light of evolution. We start this review by outlining the main paradigms of enzyme evolution, including gene duplication and divergence, convergent evolution, and evolution by recombination of domains. In the second part, we overview the current collective knowledge about enzymes, as organised by different types of data and collected in several databases. We also highlight some increasingly powerful computational tools that can be used to close gaps in understanding, in particular for types of data that require laborious experimental protocols. We believe that recent advances in protein structure prediction will be a powerful catalyst for the prediction of binding, mechanism, and ultimately, chemical reactions. A comprehensive mapping of enzyme function and evolution may be attainable in the near future.

DrugMap: A quantitative pan-cancer analysis of cysteine ligandability

Cysteine-focused chemical proteomic platforms have accelerated the clinical development of covalent inhibitors of a wide-range of targets in cancer. However, how different oncogenic contexts influence cysteine targeting remains unknown. To address this question, we have developed DrugMap , an atlas of cysteine ligandability compiled across 416 cancer cell lines. We unexpectedly find that cysteine ligandability varies across cancer cell lines, and we attribute this to differences in cellular redox states, protein conformational changes, and genetic mutations. Leveraging these findings, we identify actionable cysteines in NFκB1 and SOX10 and develop corresponding covalent ligands that block the activity of these transcription factors. We demonstrate that the NFκB1 probe blocks DNA binding, whereas the SOX10 ligand increases SOX10-SOX10 interactions and disrupts melanoma transcriptional signaling. Our findings reveal heterogeneity in cysteine ligandability across cancers, pinpoint cell-intrinsic features driving cysteine targeting, and illustrate the use of covalent probes to disrupt oncogenic transcription factor activity.

Prediction of the Peptide-TIM3 Binding Site in Inhibiting TIM3-Galectin 9 Binding Pathways

T-cell immunoglobulin and mucin domain-containing protein-3 (TIM3) is an important receptor protein that modulates the immune system. The binding of TIM3 with Galectin 9 (GAL9) triggers immune system suppression, but the TIM3-GAL9 binding can be inhibited by binding of the peptide P26 to TIM3. A fast and accurate prediction of the P26-TIM3 binding site is crucial and a prerequisite for the investigation of P26-TIM3 interactions and TIM3-GAL9 binding pathways. Here, we present a machine learning approach, which considers protein conformational changes, to quickly identify the ligand-binding site on TIM3. Our results show that the P26 binding site is located near the C″-D loop of TIM3. Further simulations show that the binding pose is stabilized by a variety of electrostatic and hydrophobic interactions. Binding of P26 can alter the conformations of nearby glycan side chains on TIM3, providing possible mechanisms of how P26 inhibits TIM3-GAL9 binding pathways. The insights from this work will facilitate the identification of other peptides or antibodies that may also inhibit the TIM3-GAL9 pathways and eventually lead to improved attempts in the modulation of the TIM3-GAL9 immunosuppression pathways. The strategies and machine learning method can be generalized to study ligand-receptor binding when the conformational changes during the binding are important.

Protein-DNA Binding Residues Prediction Using a Deep Learning Model With Hierarchical Feature Extraction

Biologically important effects occur when proteins bind to other substances, of which binding to DNA is a crucial one. Therefore, accurate identification of protein-DNA binding residues is important for further understanding of the protein-DNA interaction mechanism. Although wet-lab methods can accurately obtain the location of bound residues, it requires significant human, financial and time costs. There is thus an urgent need to develop efficient computational-based methods. Most current state-of-the-art methods are two-step approaches: the first step uses a sliding window technique to extract residue features; the second step uses each residue as an input to the model for prediction. This has a negative impact on the efficiency of prediction and ease of use. In this study, we propose a sequence-to-sequence (seq2seq) model that can input the entire protein sequence of variable length and use two modules, Transformer Encoder Block and Feature Extracting Block, for hierarchical feature extraction, where Transformer Encoder Block is used to extract global features, and then Feature Extracting Block is used to extract local features to further improve the recognition capability of the model. The comparison results on two benchmark datasets, namely PDNA-543 and PDNA-41, prove the effectiveness of our method in identifying protein-DNA binding residues.

Modified Curcuminoid-Rich Extract Liposomal CRE-SDInhibits Osteoclastogenesis via the Canonical NF-κB Signaling Pathway

Curcuminoids, namely curcumin, demethoxycurcumin, and bisdemethoxycurcumin, are the major active compounds found in Curcuma longa L. (turmeric). Although their suppressive effects on bone resorption have been demonstrated, their pharmacokinetic disadvantages remain a concern. Herein, we utilized solid dispersion of a curcuminoid-rich extract (CRE), comprising such curcuminoids, to prepare CRE-SD; subsequently, we performed liposome encapsulation of the CRE-SD to yield liposomal CRE-SD. In vitro release assessment revealed that a lower cumulative mass percentage of CRE-SD was released from liposomal CRE-SD than from CRE-SD samples. After culture of murine RANKL-stimulated RAW 264.7 macrophages, our in vitro examinations confirmed that liposomal CRE-SD may impede osteoclastogenesis by suppressing p65 and IκBα phosphorylation, together with nuclear translocation and transcriptional activity of phosphorylated p65. Blind docking simulations showed the high binding affinity between curcuminoids and the IκBα/p50/p65 protein complex, along with many intermolecular interactions, which corroborated our in vitro findings. Therefore, liposomal CRE-SD can inhibit osteoclastogenesis via the canonical NF-κB signaling pathway, suggesting its pharmacological potential for treating bone diseases with excessive osteoclastogenesis.

Active components and molecular mechanisms of Sagacious Confucius' Pillow Elixir to treat cognitive impairment based on systems pharmacology

BACKGROUND

Sagacious Confucius' Pillow Elixir (SCPE) is a common clinical prescription to treat cognitive impairment (CI) in East Asia.

OBJECTIVE

To predict the active components of SCPE, identify the associated signaling pathway, and explore the molecular mechanism using systems pharmacology and an animal study.

METHODS

Systems pharmacology and Python programming language-based molecular docking were used to select and analyze the active components and targets. Senescence-accelerated prone 8 mice were used as a CI model. The molecular mechanism was evaluated using the water maze test, neuropathological observation, cerebrospinal fluid microdialysis, and Western blotting.

RESULTS

Thirty active components were revealed by screening relevant databases and performing topological analysis. Additionally, 376 differentially expressed genes for CI were identified. Pathway enrichment analysis, protein-protein interaction (PPI) network analysis and molecular docking indicated that SCPE played a crucial role in modulating the PI3K/Akt/mTOR signaling pathway, and 23 SCPE components interacted with it. In the CI model, SCPE improved cognitive function, increased the levels of the neurotransmitter 5-hydroxytryptamine (5-HT) and metabolite 5-hydroxyindole acetic acid (5-HIAA), ameliorated pathological damage and regulated the PI3K/AKT/mTOR signaling pathway. SCPE increased the LC3-II/LC3-I, p-PI3K p85/PI3K p85, p-AKT/AKT, and p-mTOR/mTOR protein expression ratios and inhibited P62 expression in the hippocampal tissue of the CI model.

CONCLUSION

Our study revealed that 23 active SCPE components improve CI by increasing the levels of the neurotransmitter 5-HT and metabolite 5-HIAA, suppressing pathological injury and regulating the PI3K/Akt/mTOR signaling pathway to improve cognitive function.

Deep Learning Approach for the Discovery of Tumor-Targeting Small Organic Ligands from DNA-Encoded Chemical Libraries

DNA-Encoded Chemical Libraries (DELs) have emerged as efficient and cost-effective ligand discovery tools, which enable the generation of protein-ligand interaction data of unprecedented size. In this article, we present an approach that combines DEL screening and instance-level deep learning modeling to identify tumor-targeting ligands against carbonic anhydrase IX (CAIX), a clinically validated marker of hypoxia and clear cell renal cell carcinoma. We present a new ligand identification and hit-to-lead strategy driven by machine learning models trained on DELs, which expand the scope of DEL-derived chemical motifs. CAIX-screening datasets obtained from three different DELs were used to train machine learning models for generating novel hits, dissimilar to elements present in the original DELs. Out of the 152 novel potential hits that were identified with our approach and screened in an in vitro enzymatic inhibition assay, 70% displayed submicromolar activities (IC₅₀ < 1 μM). To generate lead compounds that are functionalized with anticancer payloads, analogues of top hits were prioritized for synthesis based on the predicted CAIX affinity and synthetic feasibility. Three lead candidates showed accumulation on the surface of CAIX-expressing tumor cells in cellular binding assays. The best compound displayed an in vitro K_D of 5.7 nM and selectively targeted tumors in mice bearing human renal cell carcinoma lesions. Our results demonstrate the synergy between DEL and machine learning for the identification of novel hits and for the successful translation of lead candidates for in vivo targeting applications.

Analyzing the Geometry and Dynamics of Viral Structures: A Review of Computational Approaches Based on Alpha Shape Theory, Normal Mode Analysis, and Poisson-Boltzmann Theories

The current SARS-CoV-2 pandemic highlights our fragility when we are exposed to emergent viruses either directly or through zoonotic diseases. Fortunately, our knowledge of the biology of those viruses is improving. In particular, we have more and more structural information on virions, i.e., the infective form of a virus that includes its genomic material and surrounding protective capsid, and on their gene products. It is important to have methods that enable the analyses of structural information on such large macromolecular systems. We review some of those methods in this paper. We focus on understanding the geometry of virions and viral structural proteins, their dynamics, and their energetics, with the ambition that this understanding can help design antiviral agents. We discuss those methods in light of the specificities of those structures, mainly that they are huge. We focus on three of our own methods based on the alpha shape theory for computing geometry, normal mode analyses to study dynamics, and modified Poisson-Boltzmann theories to study the organization of ions and co-solvent and solvent molecules around biomacromolecules. The corresponding software has computing times that are compatible with the use of regular desktop computers. We show examples of their applications on some outer shells and structural proteins of the West Nile Virus.

Updates on drug designing approach through computational strategies: a review

The drug discovery and development (DDD) process in pursuit of novel drug candidates is a challenging procedure requiring lots of time and resources. Therefore, computer-aided drug design (CADD) methodologies are used extensively to promote proficiency in drug development in a systematic and time-effective manner. The point in reference is SARS-CoV-2 which has emerged as a global pandemic. In the absence of any confirmed drug moiety to treat the infection, the science fraternity adopted hit and trial methods to come up with a lead drug compound. This article is an overview of the virtual methodologies, which assist in finding novel hits and help in the progression of drug development in a short period with a specific medicinal solution.

Cell line-specific features of 3D chromatin organization in hepatocellular carcinoma

Liver cancer, particularly hepatocellular carcinoma (HCC), poses a significant global threat to human lives. To advance the development of innovative diagnostic and treatment approaches, it is essential to examine the hidden features of HCC, particularly its 3D genome architecture, which is not well understood. In this study, we investigated the 3D genome organization of four HCC cell lines-Hep3B, Huh1, Huh7, and SNU449-using in situ Hi-C and assay for transposase-accessible chromatin sequencing. Our findings revealed that HCC cell lines had more long-range interactions, both intra-and interchromosomal, compared to human mammary epithelial cells (HMECs). Unexpectedly, HCC cell lines displayed cell line-specific compartmental modifications at the megabase (Mb) scale, which could potentially be leveraged in determining HCC subtypes. At the sub-Mb scale, we observed decreases in intra-TAD (topologically associated domain) interactions and chromatin loops in HCC cell lines compared to HMECs. Lastly, we discovered a correlation between gene expression and the 3D chromatin architecture of SLC8A1, which encodes a sodium-calcium antiporter whose modulation is known to induce apoptosis by comparison between HCC cell lines and HMECs. Our findings suggest that HCC cell lines have a distinct 3D genome organization that is different from those of normal and other cancer cells based on the analysis of compartments, TADs, and chromatin loops. Overall, we take this as evidence that genome organization plays a crucial role in cancer phenotype determination. Further exploration of epigenetics in HCC will help us to better understand specific gene regulation mechanisms and uncover novel targets for cancer treatment.

Molecular Structural Analysis of Porcine CMAH-Native Ligand Complex and High Throughput Virtual Screening to Identify Novel Inhibitors

Porcine meat is the most consumed red meat worldwide. Pigs are also vital tools in biological and medical research. However, xenoreactivity between porcine's N-glycolylneuraminic acid (Neu5Gc) and human anti-Neu5Gc antibodies poses a significant challenge. On the one hand, dietary Neu5Gc intake has been connected to particular human disorders. On the other hand, some pathogens connected to pig diseases have a preference for Neu5Gc. The Cytidine monophospho-N-acetylneuraminic acid hydroxylase (CMAH) catalyses the conversion of N-acetylneuraminic acid (Neu5Ac) to Neu5Gc. In this study, we predicted the tertiary structure of CMAH, performed molecular docking, and analysed the protein-native ligand complex. We performed a virtual screening from a drug library of 5M compounds and selected the two top inhibitors with Vina scores of -9.9 kcal/mol for inhibitor 1 and -9.4 kcal/mol for inhibitor 2. We further analysed their pharmacokinetic and pharmacophoric properties. We conducted stability analyses of the complexes with molecular dynamic simulations of 200 ns and binding free energy calculations. The overall analyses revealed the inhibitors' stable binding, which was further validated by the MMGBSA studies. In conclusion, this result may pave the way for future studies to determine how to inhibit CMAH activities. Further in vitro studies can provide in-depth insight into these compounds' therapeutic potential.

A bioinformatics approach to characterize a hypothetical protein Q6S8D9_SARS of SARS-CoV

Characterization as well as prediction of the secondary and tertiary structure of hypothetical proteins from their amino acid sequences uploaded in databases by in silico approach are the critical issues in computational biology. Severe acute respiratory syndrome-associated coronavirus (SARS-CoV), which is responsible for pneumonia alike diseases, possesses a wide range of proteins of which many are still uncharacterized. The current study was conducted to reveal the physicochemical characteristics and structures of an uncharacterized protein Q6S8D9_SARS of SARS-CoV. Following the common flowchart of characterizing a hypothetical protein, several sophisticated computerized tools e.g., ExPASy Protparam, CD Search, SOPMA, PSIPRED, HHpred, etc. were employed to discover the functions and structures of Q6S8D9_SARS. After delineating the secondary and tertiary structures of the protein, some quality evaluating tools e.g., PROCHECK, ProSA-web etc. were performed to assess the structures and later the active site was identified also by CASTp v.3.0. The protein contains more negatively charged residues than positively charged residues and a high aliphatic index value which make the protein more stable. The 2D and 3D structures modeled by several bioinformatics tools ensured that the proteins had domain in it which indicated it was functional protein having the ability to trouble host antiviral inflammatory cytokine and interferon production pathways. Moreover, active site was found in the protein where ligand could bind. The study was aimed to unveil the features and structures of an uncharacterized protein of SARS-CoV which can be a therapeutic target for development of vaccines against the virus. Further research are needed to accomplish the task.

Advancing Targeted Protein Degradation via Multiomics Profiling and Artificial Intelligence

Only around 20% of the human proteome is considered to be druggable with small-molecule antagonists. This leaves some of the most compelling therapeutic targets outside the reach of ligand discovery. The concept of targeted protein degradation (TPD) promises to overcome some of these limitations. In brief, TPD is dependent on small molecules that induce the proximity between a protein of interest (POI) and an E3 ubiquitin ligase, causing ubiquitination and degradation of the POI. In this perspective, we want to reflect on current challenges in the field, and discuss how advances in multiomics profiling, artificial intelligence, and machine learning (AI/ML) will be vital in overcoming them. The presented roadmap is discussed in the context of small-molecule degraders but is equally applicable for other emerging proximity-inducing modalities.

CAPLA: improved prediction of protein-ligand binding affinity by a deep learning approach based on a cross-attention mechanism

MOTIVATION

Accurate and rapid prediction of protein-ligand binding affinity is a great challenge currently encountered in drug discovery. Recent advances have manifested a promising alternative in applying deep learning-based computational approaches for accurately quantifying binding affinity. The structure complementarity between protein-binding pocket and ligand has a great effect on the binding strength between a protein and a ligand, but most of existing deep learning approaches usually extracted the features of pocket and ligand by these two detached modules.

RESULTS

In this work, a new deep learning approach based on the cross-attention mechanism named CAPLA was developed for improved prediction of protein-ligand binding affinity by learning features from sequence-level information of both protein and ligand. Specifically, CAPLA employs the cross-attention mechanism to capture the mutual effect of protein-binding pocket and ligand. We evaluated the performance of our proposed CAPLA on comprehensive benchmarking experiments on binding affinity prediction, demonstrating the superior performance of CAPLA over state-of-the-art baseline approaches. Moreover, we provided the interpretability for CAPLA to uncover critical functional residues that contribute most to the binding affinity through the analysis of the attention scores generated by the cross-attention mechanism. Consequently, these results indicate that CAPLA is an effective approach for binding affinity prediction and may contribute to useful help for further consequent applications.

AVAILABILITY AND IMPLEMENTATION

The source code of the method along with trained models is freely available at https://github.com/lennylv/CAPLA.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Application of machine learning on understanding biomolecule interactions in cellular machinery

Machine learning (ML) applications have become ubiquitous in all fields of research including protein science and engineering. Apart from protein structure and mutation prediction, scientists are focusing on knowledge gaps with respect to the molecular mechanisms involved in protein binding and interactions with other components in the experimental setups or the human body. Researchers are working on several wet-lab techniques and generating data for a better understanding of concepts and mechanics involved. The information like biomolecular structure, binding affinities, structure fluctuations and movements are enormous which can be handled and analyzed by ML. Therefore, this review highlights the significance of ML in understanding the biomolecular interactions while assisting in various fields of research such as drug discovery, nanomedicine, nanotoxicity and material science. Hence, the way ahead would be to force hand-in hand of laboratory work and computational techniques.

COLLAPSE: A representation learning framework for identification and characterization of protein structural sites

The identification and characterization of the structural sites which contribute to protein function are crucial for understanding biological mechanisms, evaluating disease risk, and developing targeted therapies. However, the quantity of known protein structures is rapidly outpacing our ability to functionally annotate them. Existing methods for function prediction either do not operate on local sites, suffer from high false positive or false negative rates, or require large site-specific training datasets, necessitating the development of new computational methods for annotating functional sites at scale. We present COLLAPSE (Compressed Latents Learned from Aligned Protein Structural Environments), a framework for learning deep representations of protein sites. COLLAPSE operates directly on the 3D positions of atoms surrounding a site and uses evolutionary relationships between homologous proteins as a self-supervision signal, enabling learned embeddings to implicitly capture structure-function relationships within each site. Our representations generalize across disparate tasks in a transfer learning context, achieving state-of-the-art performance on standardized benchmarks (protein-protein interactions and mutation stability) and on the prediction of functional sites from the Prosite database. We use COLLAPSE to search for similar sites across large protein datasets and to annotate proteins based on a database of known functional sites. These methods demonstrate that COLLAPSE is computationally efficient, tunable, and interpretable, providing a general-purpose platform for computational protein analysis.

In Silico Structural and Functional Analyses of NLRP3 Inflammasomes to Provide Insights for Treating Neurodegenerative Diseases

Inflammasomes are cytoplasmic intracellular multiprotein complexes that control the innate immune system's activation of inflammation in response to derived chemicals. Recent advancements increased our molecular knowledge of activation of NLRP3 inflammasomes. Although several studies have been done to investigate the role of inflammasomes in innate immunity and other diseases, structural, functional, and evolutionary investigations are needed to further understand the clinical consequences of NLRP3 gene. The purpose of this study is to investigate the structural and functional impact of the NLRP3 protein by using a computational analysis to uncover putative protein sites involved in the stabilization of the protein-ligand complexes with inhibitors. This will allow for a deeper understanding of the molecular mechanism underlying these interactions. It was found that human NLRP3 gene coexpresses with PYCARD, NLRC4, CASP1, MAVS, and CTSB based on observed coexpression of homologs in other species. The NACHT, LRR, and PYD domain-containing protein 3 is a key player in innate immunity and inflammation as the sensor subunit of the NLRP3 inflammasome. The inflammasome polymeric complex, consisting of NLRP3, PYCARD, and CASP1, is formed in response to pathogens and other damage-associated signals (and possibly CASP4 and CASP5). Comprehensive structural and functional analyses of NLRP3 inflammasome components offer a fresh approach to the development of new treatments for a wide variety of human disorders.