1
|
Nguyen TH, Thai QM, Pham MQ, Minh PTH, Phung HTT. Machine learning combines atomistic simulations to predict SARS-CoV-2 Mpro inhibitors from natural compounds. Mol Divers 2024; 28:553-561. [PMID: 36823394 PMCID: PMC9950021 DOI: 10.1007/s11030-023-10601-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2022] [Accepted: 01/04/2023] [Indexed: 02/25/2023]
Abstract
To date, the COVID-19 pandemic has still been infectious around the world, continuously causing social and economic damage on a global scale. One of the most important therapeutic targets for the treatment of COVID-19 is the main protease (Mpro) of SARS-CoV-2. In this study, we combined machine-learning (ML) model with atomistic simulations to computationally search for highly promising SARS-CoV-2 Mpro inhibitors from the representative natural compounds of the National Cancer Institute (NCI) Database. First, the trained ML model was used to scan the library quickly and reliably for possible Mpro inhibitors. The ML output was then confirmed using atomistic simulations integrating molecular docking and molecular dynamic simulations with the linear interaction energy scheme. The results turned out to show that there was evidently good agreement between ML and atomistic simulations. Ten substances were proposed to be able to inhibit SARS-CoV-2 Mpro. Seven of them have high-nanomolar affinity and are very potential inhibitors. The strategy has been proven to be reliable and appropriate for fast prediction of SARS-CoV-2 Mpro inhibitors, benefiting for new emerging SARS-CoV-2 variants in the future accordingly.
Collapse
Affiliation(s)
- Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Quynh Mai Thai
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Minh Quan Pham
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Pham Thi Hong Minh
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Huong Thi Thu Phung
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City, Vietnam
| |
Collapse
|
2
|
Zuo K, Kranjc A, Capelli R, Rossetti G, Nechushtai R, Carloni P. Metadynamics simulations of ligands binding to protein surfaces: a novel tool for rational drug design. Phys Chem Chem Phys 2023; 25:13819-13824. [PMID: 37184538 DOI: 10.1039/d3cp01388j] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Structure-based drug design protocols may encounter difficulties to investigate poses when the biomolecular targets do not exhibit typical binding pockets. In this study, by providing two concrete examples from our labs, we suggest that the combination of metadynamics free energy methods (validated against affinity measurements), along with experimental structural information (by X-ray crystallography and NMR), can help to identify the poses of ligands on protein surfaces. The simulation workflow proposed here was implemented in a widely used code, namely GROMACS, and it could straightforwardly be applied to various drug-design campaigns targeting ligands' binding to protein surfaces.
Collapse
Affiliation(s)
- Ke Zuo
- Computational Biomedicine, Institute of Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52425, Germany.
- Department of Physics, RWTH Aachen University, Aachen 52074, Germany
- The Alexander Silberman Institute of Life Science, The Hebrew University of Jerusalem, Edmond J. Safra Campus at Givat Ram, Jerusalem 91904, Israel
- Department of Physics, Università degli Studi di Ferrara, Ferrara 44121, Italy
| | - Agata Kranjc
- Computational Biomedicine, Institute of Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52425, Germany.
| | - Riccardo Capelli
- Department of Biosciences, Università degli Studi di Milano, Via Celoria 26, Milan 20133, Italy
| | - Giulia Rossetti
- Computational Biomedicine, Institute of Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52425, Germany.
- Jülich Supercomputing Center (JSC), Forschungszentrum Jülich GmbH, Jülich 52425, Germany
- Department of Neurology, Faculty of Medicine, RWTH Aachen University, Aachen 52074, Germany
| | - Rachel Nechushtai
- The Alexander Silberman Institute of Life Science, The Hebrew University of Jerusalem, Edmond J. Safra Campus at Givat Ram, Jerusalem 91904, Israel
| | - Paolo Carloni
- Computational Biomedicine, Institute of Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich 52425, Germany.
- Department of Physics, RWTH Aachen University, Aachen 52074, Germany
- JARA Institute: Molecular Neuroscience and Imaging, Institute of Neuroscience and Medicine INM-11, Forschungszentrum Jülich GmbH, Jülich 52425, Germany
| |
Collapse
|
3
|
Ngo ST, Nguyen TH, Tung NT, Vu VV, Pham MQ, Mai BK. Characterizing the ligand-binding affinity toward SARS-CoV-2 Mpro via physics- and knowledge-based approaches. Phys Chem Chem Phys 2022; 24:29266-29278. [PMID: 36449268 DOI: 10.1039/d2cp04476e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Computational approaches, including physics- and knowledge-based methods, have commonly been used to determine the ligand-binding affinity toward SARS-CoV-2 main protease (Mpro or 3CLpro). Strong binding ligands can thus be suggested as potential inhibitors for blocking the biological activity of the protease. In this context, this paper aims to provide a short review of computational approaches that have recently been applied in the search for inhibitor candidates of Mpro. In particular, molecular docking and molecular dynamics (MD) simulations are usually combined to predict the binding affinity of thousands of compounds. Quantitative structure-activity relationship (QSAR) is the least computationally demanding and therefore can be used for large chemical collections of ligands. However, its accuracy may not be high. Moreover, the quantum mechanics/molecular mechanics (QM/MM) method is most commonly used for covalently binding inhibitors, which also play an important role in inhibiting the activity of SARS-CoV-2. Furthermore, machine learning (ML) models can significantly increase the searching space of ligands with high accuracy for binding affinity prediction. Physical insights into the binding process can then be confirmed via physics-based calculations. Integration of ML models into computational chemistry provides many more benefits and can lead to new therapies sooner.
Collapse
Affiliation(s)
- Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Vietnam. .,Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Advanced Institute of Materials Science, Ton Duc Thang University, Ho Chi Minh City, Vietnam. .,Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Nguyen Thanh Tung
- Institute of Materials Science, Vietnam Academy of Science and Technology, Hanoi, Vietnam. .,Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Van V Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City, Vietnam
| | - Minh Quan Pham
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam.,Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Binh Khanh Mai
- Department of Chemistry, University of Pittsburgh, Pittsburgh, PA, USA
| |
Collapse
|
4
|
Aida H, Shigeta Y, Harada R. Ligand Binding Path Sampling Based on Parallel Cascade Selection Molecular Dynamics: LB-PaCS-MD. MATERIALS 2022; 15:ma15041490. [PMID: 35208030 PMCID: PMC8878848 DOI: 10.3390/ma15041490] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/22/2021] [Revised: 02/08/2022] [Accepted: 02/10/2022] [Indexed: 01/09/2023]
Abstract
Parallel cascade selection molecular dynamics (PaCS-MD) is a rare-event sampling method that generates transition pathways between a reactant and product. To sample the transition pathways, PaCS-MD repeats short-time MD simulations from important configurations as conformational resampling cycles. In this study, PaCS-MD was extended to sample ligand binding pathways toward a target protein, which is referred to as LB-PaCS-MD. In a ligand-concentrated environment, where multiple ligand copies are randomly arranged around the target protein, LB-PaCS-MD allows for the frequent sampling of ligand binding pathways. To select the important configurations, we specified the center of mass (COM) distance between each ligand and the relevant binding site of the target protein, where snapshots generated by the short-time MD simulations were ranked by their COM distance values. From each cycle, snapshots with smaller COM distance values were selected as the important configurations to be resampled using the short-time MD simulations. By repeating conformational resampling cycles, the COM distance values gradually decreased and converged to constants, meaning that a set of ligand binding pathways toward the target protein was sampled by LB-PaCS-MD. To demonstrate relative efficiency, LB-PaCS-MD was applied to several proteins, and their ligand binding pathways were sampled more frequently than conventional MD simulations.
Collapse
Affiliation(s)
- Hayato Aida
- Graduate School of Science and Technology, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan;
| | - Yasuteru Shigeta
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan;
| | - Ryuhei Harada
- Center for Computational Sciences, University of Tsukuba, 1-1-1 Tennodai, Tsukuba 305-8577, Japan;
- Correspondence:
| |
Collapse
|
5
|
In silico screening of potential β-secretase (BACE1) inhibitors from VIETHERB database. J Mol Model 2022; 28:60. [DOI: 10.1007/s00894-022-05051-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2021] [Accepted: 02/07/2022] [Indexed: 11/25/2022]
|
6
|
Tam NM, Nguyen TH, Ngan VT, Tung NT, Ngo ST. Unbinding ligands from SARS-CoV-2 Mpro via umbrella sampling simulations. ROYAL SOCIETY OPEN SCIENCE 2022; 9:211480. [PMID: 35116157 PMCID: PMC8790385 DOI: 10.1098/rsos.211480] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/25/2021] [Accepted: 12/20/2021] [Indexed: 05/03/2023]
Abstract
The umbrella sampling (US) simulation is demonstrated to be an efficient approach for determining the unbinding pathway and binding affinity to the SARS-CoV-2 Mpro of small molecule inhibitors. The accuracy of US is in the same range as the linear interaction energy (LIE) and fast pulling of ligand (FPL) methods. In detail, the correlation coefficient between US and experiments does not differ from FPL and is slightly smaller than LIE. The root mean square error of US simulations is smaller than that of LIE. Moreover, US is better than FPL and poorer than LIE in classifying SARS-CoV-2 Mpro inhibitors owing to the reciever operating characteristic-area under the curve analysis. Furthermore, the US simulations also provide detailed insights on unbinding pathways of ligands from the binding cleft of SARS-CoV-2 Mpro. The residues Cys44, Thr45, Ser46, Leu141, Asn142, Gly143, Glu166, Leu167, Pro168, Ala191, Gln192 and Ala193 probably play an important role in the ligand dissociation. Therefore, substitutions at these points may change the mechanism of binding of inhibitors to SARS-CoV-2 Mpro.
Collapse
Affiliation(s)
- Nguyen Minh Tam
- Computational Chemistry Research Group, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Trung Hai Nguyen
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Vu Thi Ngan
- Laboratory of Computational Chemistry and Modelling, Department of Chemistry, Quy Nhon University, Quy Nhon, Vietnam
| | - Nguyen Thanh Tung
- Institute of Materials Science, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Son Tung Ngo
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| |
Collapse
|
7
|
Pham TNH, Nguyen TH, Tam NM, Y Vu T, Pham NT, Huy NT, Mai BK, Tung NT, Pham MQ, V Vu V, Ngo ST. Improving ligand-ranking of AutoDock Vina by changing the empirical parameters. J Comput Chem 2021; 43:160-169. [PMID: 34716930 DOI: 10.1002/jcc.26779] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Revised: 10/10/2021] [Accepted: 10/14/2021] [Indexed: 01/09/2023]
Abstract
AutoDock Vina (Vina) achieved a very high docking-success rate, p ^ , but give a rather low correlation coefficient, R , for binding affinity with respect to experiments. This low correlation can be an obstacle for ranking of ligand-binding affinity, which is the main objective of docking simulations. In this context, we evaluated the dependence of Vina R coefficient upon its empirical parameters. R is affected more by changing the gauss2 and rotation than other terms. The docking-success rate p ^ is sensitive to the alterations of the gauss1, gauss2, repulsion, and hydrogen bond parameters. Based on our benchmarks, the parameter set1 has been suggested to be the most optimal. The testing study over 800 complexes indicated that the modified Vina provided higher correlation with experiment R set 1 = 0.556 ± 0.025 compared with R Default = 0.493 ± 0.028 obtained by the original Vina and R Vina 1.2 = 0.503 ± 0.029 by Vina version 1.2. Besides, the modified Vina can be also applied more widely, giving R ≥ 0.500 for 32/48 targets, compared with the default package, giving R ≥ 0.500 for 31/48 targets. In addition, validation calculations for 1036 complexes obtained from version 2019 of PDBbind refined structures showed that the set1 of parameters gave higher correlation coefficient ( R set 1 = 0.617 ± 0.017 ) than the default package ( R Default = 0.543 ± 0.020 ) and Vina version 1.2 ( R Vina 1.2 = 0.540 ± 0.020 ). The version of Vina with set1 of parameters can be downloaded at https://github.com/sontungngo/mvina. The outcomes would enhance the ranking of ligand-binding affinity using Autodock Vina.
Collapse
Affiliation(s)
- T Ngoc Han Pham
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam.,Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Nguyen Minh Tam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam.,Computational Chemistry Research Group, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Thien Y Vu
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Nhat Truong Pham
- Faculty of Electrical and Electronics Engineering, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Nguyen Truong Huy
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Binh Khanh Mai
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania, USA
| | - Nguyen Thanh Tung
- Institute of Materials Science, Vietnam Academy of Science and Technology, Hanoi, Vietnam.,Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Minh Quan Pham
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam.,Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Van V Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City, Vietnam
| | - Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam.,Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| |
Collapse
|
8
|
Ngo ST, Vu KB, Pham MQ, Tam NM, Tran PT. Marine derivatives prevent wMUS81 in silico studies. ROYAL SOCIETY OPEN SCIENCE 2021; 8:210974. [PMID: 34527278 PMCID: PMC8424343 DOI: 10.1098/rsos.210974] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Accepted: 08/11/2021] [Indexed: 05/15/2023]
Abstract
The winged-helix domain of the methyl methanesulfonate and ultraviolet-sensitive 81 (wMUS81) is a potential cancer drug target. In this context, marine fungi compounds were indicated to be able to prevent wMUS81 structure via atomistic simulations. Eight compounds such as D197 (Tryptoquivaline U), D220 (Epiremisporine B), D67 (Aspergiolide A), D153 (Preussomerin G), D547 (12,13-dihydroxyfumitremorgin C), D152 (Preussomerin K), D20 (Marinopyrrole B) and D559 (Fumuquinazoline K) were indicated that they are able to prevent the conformation of wMUS81 via forming a strong binding affinity to the enzyme via perturbation approach. The electrostatic interaction is the dominant factor in the binding process of ligands to wMUS81. The residues Trp55, Arg59, Leu62, His63 and Arg69 were found to frequently form non-bonded contacts and hydrogen bonds to inhibitors. Moreover, the influence of the ligand D197, which formed the lowest binding free energy to wMUS81, on the structural change of enzyme was investigated using replica exchange molecular dynamics simulations. The obtained results indicated that D197, which forms a strong binding affinity, can modify the structure of wMUS81. Overall, the marine compounds probably inhibit wMUS81 due to forming a strong binding affinity to the enzyme as well as altering the enzymic conformation.
Collapse
Affiliation(s)
- Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Khanh B. Vu
- Department of Chemical Engineering, International University, Ho Chi Minh City, Vietnam
- Vietnam National University, Ho Chi Minh City, Vietnam
| | - Minh Quan Pham
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Nguyen Minh Tam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Phuong-Thao Tran
- Department of Pharmaceutical Chemistry, Hanoi University of Pharmacy, Hanoi, Vietnam
| |
Collapse
|
9
|
Cao DT, Huong Doan TM, Pham VC, Minh Le TH, Chae JW, Yun HY, Na MK, Kim YH, Pham MQ, Nguyen VH. Molecular design of anticancer drugs from marine fungi derivatives. RSC Adv 2021; 11:20173-20179. [PMID: 35479875 PMCID: PMC9033662 DOI: 10.1039/d1ra01855h] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/09/2021] [Accepted: 05/24/2021] [Indexed: 12/21/2022] Open
Abstract
Heat shock protein 90 (Hsp90) is one of the most potential targets in cancer therapy. We have demonstrated using a combination of molecular docking and fast pulling of ligand (FPL) simulations that marine fungi derivatives can be possible inhibitors, preventing the biological activity of Hsp90. The computational approaches were validated and compared with previous experiments. Based on the benchmark of available inhibitors of Hsp90, the GOLD docking package using the ChemPLP scoring function was found to be superior over both Autodock Vina and Autodock4 in the preliminary estimation of the ligand-binding affinity and binding pose with the Pearson correlation, R = -0.62. Moreover, FPL calculations were also indicated as a suitable approach to refine docking simulations with a correlation coefficient with the experimental data of R = -0.81. Therefore, the binding affinity of marine fungi derivatives to Hsp90 was evaluated. Docking and FPL calculations suggest that five compounds including 23, 40, 46, 48, and 52 are highly potent inhibitors for Hsp90. The obtained results enhance cancer therapy research.
Collapse
Affiliation(s)
- Duc Tuan Cao
- Hai Phong University of Medicine and Pharmacy Haiphong Vietnam
| | - Thi Mai Huong Doan
- Institute of Marine Biochemistry, Vietnam Academy of Science and Technology Hanoi Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology Hanoi Vietnam
| | - Van Cuong Pham
- Institute of Marine Biochemistry, Vietnam Academy of Science and Technology Hanoi Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology Hanoi Vietnam
| | - Thi Hong Minh Le
- Institute of Marine Biochemistry, Vietnam Academy of Science and Technology Hanoi Vietnam
| | - Jung-Woo Chae
- College of Pharmacy, Chungnam National University Daejeon Republic of Korea
| | - Hwi-Yeol Yun
- College of Pharmacy, Chungnam National University Daejeon Republic of Korea
| | - Min-Kyun Na
- College of Pharmacy, Chungnam National University Daejeon Republic of Korea
| | - Young-Ho Kim
- College of Pharmacy, Chungnam National University Daejeon Republic of Korea
| | - Minh Quan Pham
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology Hanoi Vietnam
- Institute of Natural Products Chemistry, Vietnam Academy of Science and Technology Hanoi Vietnam
| | - Van Hung Nguyen
- Hai Phong University of Medicine and Pharmacy Haiphong Vietnam
| |
Collapse
|
10
|
Ngo ST, Nguyen TH, Pham DH, Tung NT, Nam PC. Thermodynamics and kinetics in antibody resistance of the 501Y.V2 SARS-CoV-2 variant. RSC Adv 2021; 11:33438-33446. [PMID: 35497518 PMCID: PMC9042284 DOI: 10.1039/d1ra04134g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2021] [Accepted: 10/06/2021] [Indexed: 02/01/2023] Open
Abstract
Understanding the thermodynamics and kinetics of the binding process of an antibody to the SARS-CoV-2 receptor-binding domain (RBD) of the spike protein is very important for the development of COVID-19 vaccines. In particular, it is essential to understand how the binding mechanism may change under the effects of RBD mutations. In this context, we have demonstrated that the South African variant (B1.351 or 501Y.V2) can resist the neutralizing antibody (NAb). Three substitutions in the RBD including K417N, E484K, and N501Y alter the free energy landscape, binding pose, binding free energy, binding kinetics, hydrogen bonding, nonbonded contacts, and unbinding pathway of RBD + NAb complexes. The low binding affinity of NAb to 501Y.V2 RBD confirms the antibody resistance of the South African variant. Moreover, the fragment of NAb + RBD can be used as an affordable model to investigate changes in the binding process between the mutated RBD and antibodies. Increasing FEL minima of 501Y.V2 RBD + antibody in comparison with the WT RBD systems imply that the complex 501Y.V2 RBD + antibody is more unstable than the WT one.![]()
Collapse
Affiliation(s)
- Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Trung Hai Nguyen
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City, Vietnam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City, Vietnam
| | - Duc-Hung Pham
- Division of Immunobiology, Cincinnati Children's Hospital Medical Center, Cincinnati 45229, OH, USA
| | - Nguyen Thanh Tung
- Institute of Materials Science, Vietnam Academy of Science and Technology, Hanoi, Vietnam
- Graduate University of Science and Technology, Vietnam Academy of Science and Technology, Hanoi, Vietnam
| | - Pham Cam Nam
- Department of Chemical Engineering, The University of Da Nang, University of Science and Technology, Da Nang City, Vietnam
| |
Collapse
|