Abstract
Background: Computational molecular docking plays an important role in determining the precise receptor-ligand conformation, which becomes a powerful tool for drug discovery. In the past 30 years, most computational docking methods have treated the receptor structure as a rigid body, although flexible docking often yields higher accuracy. The main disadvantage of flexible docking is its significantly higher computational cost. Due to the fact that different protein pocket residues exhibit different degrees of flexibility, semi-flexible docking methods, balancing rigid docking and flexible docking, have demonstrated success in predicting highly accurate conformations with a relatively low computational cost.
Methods: In our study, the number of flexible pocket residues was assessed by quantitative analysis, and a novel adaptive residue flexibility prediction method, named A-RFP, was proposed to improve the docking performance. Based on the homologous information, a joint strategy is used to predict the pocket residue flexibility by combining RMSD, the distance between the residue sidechain and the ligand, and the sidechain orientation. For each receptor-ligand pair, A-RFP provides a docking conformation with the optimal affinity.
Results: By analyzing the docking affinities of 3507 target-ligand pairs in 5 different values ranging from 0 to 10, we found there is a general trend that the larger number of flexible residues inevitably improves the docking results by using Autodock Vina. However, a certain number of counterexamples still exist. To validate the effectiveness of A-RFP, the experimental assessment was tested in a small-scale virtual screening on 5 proteins, which confirmed that A-RFP could enhance the docking performance. And the flexible-receptor virtual screening on a low-similarity dataset with 85 receptors validates the accuracy of residue flexibility comprehensive evaluation. Moreover, we studied three receptors with FDA-approved drugs, which further proved A-RFP can play a suitable role in ligand discovery.
Conclusion: Our analysis confirms that the screening performance of the various numbers of flexible residues varies wildly across receptors. It suggests that a fine-grained docking method would offset the aforementioned deficiency. Thus, we presented A-RFP, an adaptive pocket residue flexibility prediction method based on homologous information. Without considering computational resources and time costs, A-RFP provides the optimal docking result.
[http://dx.doi.org/10.1371/journal.pone.0117955] [PMID: 25822509]
[http://dx.doi.org/10.1038/nrd3478] [PMID: 21804595]
[PMID: 33663346]
[http://dx.doi.org/10.1016/j.bmc.2014.05.068] [PMID: 24961874]
[http://dx.doi.org/10.1007/978-1-4939-6952-4_20]
[http://dx.doi.org/10.1016/S0958-1669(99)00051-8] [PMID: 10679349]
[http://dx.doi.org/10.1186/s13321-018-0320-9] [PMID: 30552524]
[http://dx.doi.org/10.1002/jcc.23108] [PMID: 22961860]
[http://dx.doi.org/10.1142/S0219720015410073] [PMID: 25800162]
[http://dx.doi.org/10.1002/jcc.21256] [PMID: 19399780]
[http://dx.doi.org/10.1002/jcc.21334] [PMID: 19499576]
[http://dx.doi.org/10.1371/journal.pcbi.1004586] [PMID: 26629955]
[http://dx.doi.org/10.1111/cbdd.13764] [PMID: 32679606]
[http://dx.doi.org/10.1038/nchem.1954] [PMID: 24950326]
[http://dx.doi.org/10.1007/s10822-011-9523-0] [PMID: 22246345]
[http://dx.doi.org/10.1016/j.jmb.2013.03.023] [PMID: 23524133]
[http://dx.doi.org/10.1016/j.jmb.2003.11.041] [PMID: 14729347]
[http://dx.doi.org/10.1021/ci300390h] [PMID: 23005250]
[http://dx.doi.org/10.1093/bioinformatics/btl158] [PMID: 16731699]
[http://dx.doi.org/10.1002/jcc.20084] [PMID: 15264254]
[http://dx.doi.org/10.1002/(SICI)1097-0134(20000515)39:3261::AID-PROT903.0.CO;2-4] [PMID: 10737948]
[http://dx.doi.org/10.1093/nar/gky356] [PMID: 29762700]
[http://dx.doi.org/10.1039/C6CP01555G] [PMID: 27108770]
[http://dx.doi.org/10.1021/acs.jcim.8b00312] [PMID: 29989806]
[http://dx.doi.org/10.1021/jm0306430] [PMID: 15027865]
[http://dx.doi.org/10.1093/nar/gkj120] [PMID: 16381872]