1
|
Liu X, Zheng L, Qin C, Li Y, Zhang JZH, Sun Z. Screening Power of End-Point Free-Energy Calculations in Cucurbituril Host-Guest Systems. J Chem Inf Model 2023; 63:6938-6946. [PMID: 37908066 DOI: 10.1021/acs.jcim.3c01356] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
End-point free-energy methods as an indispensable component in virtual screening are commonly recognized as a tool with a certain level of screening power in pharmaceutical research. While a huge number of records could be found for end-point applications in protein-ligand, protein-protein, and protein-DNA complexes from academic and industrial reports, up to now, there is no large-scale benchmark in host-guest complexes supporting the screening power of end-point free-energy techniques. A good benchmark requires a data set of sufficient coverage of pharmaceutically relevant chemical space, a long-time sampling length supporting the trajectory approximation of the ensemble average, and a sufficient sample size of receptor-acceptor pairs to stabilize the performance statistics. In this work, selecting a popular family of macrocyclic hosts named cucurbiturils, we construct a large data set containing 154 host-guest pairs, perform extensive end-point sampling of several hundred nanosecond lengths for each system, and extract the free-energy estimates with a variety of end-point free-energy techniques, including the advanced three-trajectory dielectric-constant-variable regime proposed in our recent work. The best-performing end-point protocol employs GAFF2 for solute descriptions, the three-trajectory end-point sampling regime, and the MM/GBSA Hamiltonian in free-energy extraction, achieving a high ranking metrics of Kendall τ > 0.6, a Pearlman predictive index of ∼0.8, and a high scoring power of Pearson r > 0.8. The current project as the first large-scale systematic benchmark of end-point methods in host-guest complexes in academic publications provides solid evidence of the applicability of end-point techniques and direct guidance of computational setups in practical host-guest systems.
Collapse
Affiliation(s)
- Xiao Liu
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai 201620, China
| | - Lei Zheng
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
| | - Chu Qin
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai 201620, China
| | - Yang Li
- College of Information Science and Engineering, Shandong Agricultural University, Tai'an 271018, China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Department of Chemistry, New York University, New York, New York 10003, United States
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- Faculty of Synthetic Biology and Institute of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen Guangdong 518055, China
| | - Zhaoxi Sun
- Changping Laboratory, Beijing 102206, China
| |
Collapse
|
2
|
Liu X, Zheng L, Cong Y, Gong Z, Yin Z, Zhang JZH, Liu Z, Sun Z. Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: II. regression and dielectric constant. J Comput Aided Mol Des 2022; 36:879-894. [PMID: 36394776 DOI: 10.1007/s10822-022-00487-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/26/2022] [Accepted: 10/29/2022] [Indexed: 11/18/2022]
Abstract
End-point free energy calculations as a powerful tool have been widely applied in protein-ligand and protein-protein interactions. It is often recognized that these end-point techniques serve as an option of intermediate accuracy and computational cost compared with more rigorous statistical mechanic models (e.g., alchemical transformation) and coarser molecular docking. However, it is observed that this intermediate level of accuracy does not hold in relatively simple and prototypical host-guest systems. Specifically, in our previous work investigating a set of carboxylated-pillar[6]arene host-guest complexes, end-point methods provide free energy estimates deviating significantly from the experimental reference, and the rank of binding affinities is also incorrectly computed. These observations suggest the unsuitability and inapplicability of standard end-point free energy techniques in host-guest systems, and alteration and development are required to make them practically usable. In this work, we consider two ways to improve the performance of end-point techniques. The first one is the PBSA_E regression that varies the weights of different free energy terms in the end-point calculation procedure, while the second one is considering the interior dielectric constant as an additional variable in the end-point equation. By detailed investigation of the calculation procedure and the simulation outcome, we prove that these two treatments (i.e., regression and dielectric constant) are manipulating the end-point equation in a somehow similar way, i.e., weakening the electrostatic contribution and strengthening the non-polar terms, although there are still many detailed differences between these two methods. With the trained end-point scheme, the RMSE of the computed affinities is improved from the standard ~ 12 kcal/mol to ~ 2.4 kcal/mol, which is comparable to another altered end-point method (ELIE) trained with system-specific data. By tuning PBSA_E weighting factors with the host-specific data, it is possible to further decrease the prediction error to ~ 2.1 kcal/mol. These observations along with the extremely efficient optimized-structure computation procedure suggest the regression (i.e., PBSA_E as well as its GBSA_E extension) as a practically applicable solution that brings end-point methods back into the library of usable tools for host-guest binding. However, the dielectric-constant-variable scheme cannot effectively minimize the experiment-calculation discrepancy for absolute binding affinities, but is able to improve the calculation of affinity ranks. This phenomenon is somehow different from the protein-ligand case and suggests the difference between host-guest and biomacromolecular (protein-ligand and protein-protein) systems. Therefore, the spectrum of tools usable for protein-ligand complexes could be unsuitable for host-guest binding, and numerical validations are necessary to screen out really workable solutions in these 'prototypical' situations.
Collapse
Affiliation(s)
- Xiao Liu
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China.
| | - Lei Zheng
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| | - Yalong Cong
- School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Zhihao Gong
- School of Micro-Nano Electronics, Zhejiang University, Hangzhou, 310027, China.,Hangzhou Global Scientific and Technological Innovation Center, Zhejiang University, Hangzhou, 310027, China
| | - Zhixiang Yin
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China. .,School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China. .,Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong, China. .,Department of Chemistry, New York University, NY, NY, 10003, USA.
| | - Zhirong Liu
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China
| | - Zhaoxi Sun
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China.
| |
Collapse
|
3
|
Liu X, Zheng L, Qin C, Zhang JZH, Sun Z. Comprehensive evaluation of end-point free energy techniques in carboxylated-pillar[6]arene host-guest binding: I. Standard procedure. J Comput Aided Mol Des 2022; 36:735-752. [PMID: 36136209 DOI: 10.1007/s10822-022-00475-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2022] [Accepted: 09/06/2022] [Indexed: 10/14/2022]
Abstract
Despite the massive application of end-point free energy methods in protein-ligand and protein-protein interactions, computational understandings about their performance in relatively simple and prototypical host-guest systems are limited. In this work, we present a comprehensive benchmark calculation with standard end-point free energy techniques in a recent host-guest dataset containing 13 host-guest pairs involving the carboxylated-pillar[6]arene host. We first assess the charge schemes for solutes by comparing the charge-produced electrostatics with many ab initio references, in order to obtain a preliminary albeit detailed view of the charge quality. Then, we focus on four modelling details of end-point free energy calculations, including the docking procedure for the generation of initial condition, the charge scheme for host and guest molecules, the water model used in explicit-solvent sampling, and the end-point methods for free energy estimation. The binding thermodynamics obtained with different modelling schemes are compared with experimental references, and some practical guidelines on maximizing the performance of end-point methods in practical host-guest systems are summarized. Further, we compare our simulation outcome with predictions in the grand challenge and discuss further developments to improve the prediction quality of end-point free energy methods. Overall, unlike the widely acknowledged applicability in protein-ligand binding, the standard end-point calculations cannot produce useful outcomes in host-guest binding and thus are not recommended unless alterations are performed.
Collapse
Affiliation(s)
- Xiao Liu
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China.
| | - Lei Zheng
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China
| | - Chu Qin
- School of Mathematics, Physics and Statistics, Shanghai University of Engineering Science, Shanghai, 201620, China
| | - John Z H Zhang
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China.,School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China.,Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong, China.,Department of Chemistry, New York University, New York, NY, 10003, USA
| | - Zhaoxi Sun
- College of Chemistry and Molecular Engineering, Peking University, Beijing, 100871, China.
| |
Collapse
|
4
|
Messiad FA, Ammouchi N, Belhocine Y, Alhussain H, Ghoniem MG, Said RB, Ali FAM, Rahali S. In Search of Preferential Macrocyclic Hosts for Sulfur Mustard Sensing and Recognition: A Computational Investigation through the New Composite Method r2SCAN-3c of the Key Factors Influencing the Host-Guest Interactions. NANOMATERIALS 2022; 12:nano12152517. [PMID: 35893486 PMCID: PMC9329917 DOI: 10.3390/nano12152517] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/29/2022] [Revised: 07/18/2022] [Accepted: 07/19/2022] [Indexed: 11/29/2022]
Abstract
Sulfur mustard (SM) is a harmful warfare agent that poses a serious threat to human health and the environment. Thus, the design of porous materials capable of sensing and/or capturing SM is of utmost importance. In this paper, the interactions of SM and its derivatives with ethylpillar[5]arene (EtP[5]) and the interactions between SM and a variety of host macrocycles were investigated through molecular docking calculations and non-covalent interaction (NCI) analysis. The electronic quantum parameters were computed to assess the chemical sensing properties of the studied hosts toward SM. It was found that dispersion interactions contributed significantly to the overall complexation energy, leading to the stabilization of the investigated systems. DFT energy computations showed that SM was more efficiently complexed with DCMP[5] than the other hosts studied here. Furthermore, the studied macrocyclic containers could be used as host-based chemical sensors or receptors for SM. These findings could motivate experimenters to design efficient sensing and capturing materials for the detection of SM and its derivatives.
Collapse
Affiliation(s)
- Fatine Ali Messiad
- Department of Process Engineering, Faculty of Technology, Université 20 Août 1955, El Hadaik Road, Skikda 21000, Algeria;
- LRPCSI-Laboratoire de Recherche sur la Physico-Chimie des Surfaces et Interfaces, Université 20 Août 1955, Skikda 21000, Algeria
| | - Nesrine Ammouchi
- LRPCSI-Laboratoire de Recherche sur la Physico-Chimie des Surfaces et Interfaces, Université 20 Août 1955, Skikda 21000, Algeria
- Département de Technologie, Faculté de Technologie, Université 20 Août 1955, B.P. 26, Route d’El Hadaiek, Skikda 21000, Algeria
- Correspondence: (N.A.); (Y.B.); (S.R.)
| | - Youghourta Belhocine
- Department of Process Engineering, Faculty of Technology, Université 20 Août 1955, El Hadaik Road, Skikda 21000, Algeria;
- Correspondence: (N.A.); (Y.B.); (S.R.)
| | - Hanan Alhussain
- Department of Chemistry, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia; (H.A.); (M.G.G.); (F.A.M.A.)
| | - Monira Galal Ghoniem
- Department of Chemistry, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia; (H.A.); (M.G.G.); (F.A.M.A.)
| | - Ridha Ben Said
- Department of Chemistry, College of Science and Arts, Qassim University, P.O. 53, Ar Rass 51921, Saudi Arabia;
- Laboratoire de Caractérisations, Applications et Modélisations des Matériaux, Faculté des Sciences de Tunis, Université Tunis El Manar, Tunis 2092, Tunisia
| | - Fatima Adam Mohamed Ali
- Department of Chemistry, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia; (H.A.); (M.G.G.); (F.A.M.A.)
| | - Seyfeddine Rahali
- Department of Chemistry, College of Science and Arts, Qassim University, P.O. 53, Ar Rass 51921, Saudi Arabia;
- Correspondence: (N.A.); (Y.B.); (S.R.)
| |
Collapse
|
5
|
Unveiling the nature of interactions in the supramolecular complex of Cucurbit [7-8] uril with ascorbic acid, dopamine and uric acid for biosensing applications: a computational study. Struct Chem 2022. [DOI: 10.1007/s11224-022-01945-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
6
|
Adam FA, Ghoniem MG, Diawara M, Rahali S, Abdulkhair BY, Elamin MR, Ben Aissa MA, Seydou M. Enhanced adsorptive removal of indigo carmine dye by bismuth oxide doped MgO based adsorbents from aqueous solution: equilibrium, kinetic and computational studies. RSC Adv 2022; 12:24786-24803. [PMID: 36199891 PMCID: PMC9431457 DOI: 10.1039/d2ra02636h] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/25/2022] [Accepted: 08/24/2022] [Indexed: 01/03/2023] Open
Abstract
The IC adsorption mechanism on the Bi2O3 doped MgO nanosorbents occurred through the chemisorption process.
Collapse
Affiliation(s)
- Fatima A. Adam
- Department of Chemistry, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
| | - M. G. Ghoniem
- Department of Chemistry, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
| | - Moussa Diawara
- Laboratoire de Centre de Calcul de Modélisation et de Simulation (CCMS), DER de Physique de La Faculté des Sciences et Techniques (FST), Université des Sciences des Techniques et des Technologies de Bamako (USTTB-Mali), Bamako, Mali
| | - Seyfeddine Rahali
- Department of Chemistry, College of Science and Arts, Qassim University, Ar Rass, Saudi Arabia
| | - Babiker Y. Abdulkhair
- Department of Chemistry, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
| | - M. R. Elamin
- Department of Chemistry, College of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
| | - Mohamed Ali Ben Aissa
- Department of Chemistry, College of Science and Arts, Qassim University, Ar Rass, Saudi Arabia
| | - Mahamadou Seydou
- Université de Paris, ITODYS, CNRS, UMR 7086, 15 Rue J-A de Baïf, 75013, Paris, France
| |
Collapse
|