Abstract
Docking simulation, a key technique in virtual screening, typically treats proteins as rigid bodies. However, proteins are inherently flexible, and ligand binding can induce significant conformational changes, affecting prediction accuracy. This study proposes a new approach to identify protein binding pockets that exhibit substantial conformational changes upon lig-and binding, potentially impacting docking simulation accuracy. In this research, we developed a prediction model using graph neural network to identify protein pockets with large conformational changes. To train the model, we constructed a dataset by calculating conformational changes in ligand-binding sites between multiple holo structures corresponding to the apo structure. We evaluated the performance of the prediction model and the results demonstrated that our model could iden-tify proteins with significant conformational changes, although the prediction accuracy remains low, with an AUC of 0.58 on the test data. This study highlights the potential of deep learning approaches in addressing the challenges of protein flexi-bility in virtual screening and docking simulations.
Supplementary materials
Title
supporting infomation
Description
This document contains details on the method of constructing the dataset, the impact of data augmentation techniques on learning, and a description of the CSV file for the dataset constructed this study.
Actions
Title
Dataset of Pocket RMSD between apo structure and holo structure
Description
We provide a CSV file containing the root-mean-square deviation (RMSD) values of the pockets between apo and holo structures of proteins. The file includes columns for the PDB ID of the apo protein (apo_name), the chain name of the apo protein (apo_chain), the PDB ID of the holo protein (holo_name), the chain name of the holo protein (holo_chain), the RMSD value of the pocket between the apo and holo proteins (pocket_rmsd), an identifier for the protein (protein_id), an identifier for the group of proteins with sequence similarity of 50% or more (family50_id), the name of the ligand (ligand), the number of atoms in the ligand (ligand_atom_count), and the number of residues in the pocket of the holo protein (pocket_residue_count). The ligand column is blank for ligands formed by polymers, and the ligand_atom_count is zero for these ligands, as they cannot be specified.
Actions