1
|
Ferreiro D, Khalil R, Sousa SF, Arenas M. Substitution Models of Protein Evolution with Selection on Enzymatic Activity. Mol Biol Evol 2024; 41:msae026. [PMID: 38314876 PMCID: PMC10873502 DOI: 10.1093/molbev/msae026] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2023] [Revised: 01/25/2024] [Accepted: 01/31/2024] [Indexed: 02/07/2024] Open
Abstract
Substitution models of evolution are necessary for diverse evolutionary analyses including phylogenetic tree and ancestral sequence reconstructions. At the protein level, empirical substitution models are traditionally used due to their simplicity, but they ignore the variability of substitution patterns among protein sites. Next, in order to improve the realism of the modeling of protein evolution, a series of structurally constrained substitution models were presented, but still they usually ignore constraints on the protein activity. Here, we present a substitution model of protein evolution with selection on both protein structure and enzymatic activity, and that can be applied to phylogenetics. In particular, the model considers the binding affinity of the enzyme-substrate complex as well as structural constraints that include the flexibility of structural flaps, hydrogen bonds, amino acids backbone radius of gyration, and solvent-accessible surface area that are quantified through molecular dynamics simulations. We applied the model to the HIV-1 protease and evaluated it by phylogenetic likelihood in comparison with the best-fitting empirical substitution model and a structurally constrained substitution model that ignores the enzymatic activity. We found that accounting for selection on the protein activity improves the fitting of the modeled functional regions with the real observations, especially in data with high molecular identity, which recommends considering constraints on the protein activity in the development of substitution models of evolution.
Collapse
Affiliation(s)
- David Ferreiro
- CINBIO, Universidade de Vigo, 36310 Vigo, Spain
- Department of Biochemistry, Genetics and Immunology, Universidade de Vigo, 36310 Vigo, Spain
| | - Ruqaiya Khalil
- CINBIO, Universidade de Vigo, 36310 Vigo, Spain
- Department of Biochemistry, Genetics and Immunology, Universidade de Vigo, 36310 Vigo, Spain
| | - Sergio F Sousa
- UCIBIO/REQUIMTE, BioSIM, Departamento de Biomedicina, Faculdade de Medicina da Universidade do Porto, 4200-319 Porto, Portugal
| | - Miguel Arenas
- CINBIO, Universidade de Vigo, 36310 Vigo, Spain
- Department of Biochemistry, Genetics and Immunology, Universidade de Vigo, 36310 Vigo, Spain
| |
Collapse
|
2
|
Dou B, Zhu Z, Merkurjev E, Ke L, Chen L, Jiang J, Zhu Y, Liu J, Zhang B, Wei GW. Machine Learning Methods for Small Data Challenges in Molecular Science. Chem Rev 2023; 123:8736-8780. [PMID: 37384816 PMCID: PMC10999174 DOI: 10.1021/acs.chemrev.3c00189] [Citation(s) in RCA: 21] [Impact Index Per Article: 21.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/01/2023]
Abstract
Small data are often used in scientific and engineering research due to the presence of various constraints, such as time, cost, ethics, privacy, security, and technical limitations in data acquisition. However, big data have been the focus for the past decade, small data and their challenges have received little attention, even though they are technically more severe in machine learning (ML) and deep learning (DL) studies. Overall, the small data challenge is often compounded by issues, such as data diversity, imputation, noise, imbalance, and high-dimensionality. Fortunately, the current big data era is characterized by technological breakthroughs in ML, DL, and artificial intelligence (AI), which enable data-driven scientific discovery, and many advanced ML and DL technologies developed for big data have inadvertently provided solutions for small data problems. As a result, significant progress has been made in ML and DL for small data challenges in the past decade. In this review, we summarize and analyze several emerging potential solutions to small data challenges in molecular science, including chemical and biological sciences. We review both basic machine learning algorithms, such as linear regression, logistic regression (LR), k-nearest neighbor (KNN), support vector machine (SVM), kernel learning (KL), random forest (RF), and gradient boosting trees (GBT), and more advanced techniques, including artificial neural network (ANN), convolutional neural network (CNN), U-Net, graph neural network (GNN), Generative Adversarial Network (GAN), long short-term memory (LSTM), autoencoder, transformer, transfer learning, active learning, graph-based semi-supervised learning, combining deep learning with traditional machine learning, and physical model-based data augmentation. We also briefly discuss the latest advances in these methods. Finally, we conclude the survey with a discussion of promising trends in small data challenges in molecular science.
Collapse
Affiliation(s)
- Bozheng Dou
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Zailiang Zhu
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Ekaterina Merkurjev
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Lu Ke
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Long Chen
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Jian Jiang
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
| | - Yueying Zhu
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Jie Liu
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Bengong Zhang
- Research Center of Nonlinear Science, School of Mathematical and Physical Sciences,Wuhan Textile University, Wuhan 430200, P, R. China
| | - Guo-Wei Wei
- Department of Mathematics, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Electrical and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
3
|
Li LP, Li HX, Zhou H, Li WY, Wang RL, Zhang YC, Ma Y. Exploring the mechanism of C473D mutation on CDC25B causing weak binding affinity with CDK2/CyclinA by molecular dynamics study. J Biomol Struct Dyn 2023; 41:12552-12564. [PMID: 36655391 DOI: 10.1080/07391102.2023.2166995] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2022] [Accepted: 01/05/2023] [Indexed: 01/20/2023]
Abstract
CDC25B belongs to the CDC25 family, and it plays an important part in regulating the activity of CDK/CyclinA. Studies have shown that CDC25B is closely related to cancer development. When CYS473 on CDC25B is mutated into ASP, the affinity between CDC25B and CDK2/CyclinA weakens, and their dissociation speed is greatly improved. However, the mechanism by which the CDC25BC473D mutant weakens its binding to CDK2/CyclinA is unclear. In order to study the effect of CDC25BC473D mutants on CDK2/CyclinA substrates, we constructed and verified the rationality of the CDC25BWT:CDK2/CyclinA system and CDC25BC473D:CDK2/CyclinA system and conducted molecular dynamics (MD) simulation analysis. In the post-analysis, the fluctuations of residues ARG488-SER499, LYS541-TRP550 on CDC25B and residues ASP206-ASP210 on CDK2 were massive in the mutant CDC25BC473D:CDK2/CyclinA system. And the interactions between residue ARG492 and residue GLU208, residue ARG544 and residue GLU42, residue ARG544 and TRP550 were weakened in the mutant CDC25BC473D:CDK2/CyclinA system. The results showed that when CYS473 on CDC25B was mutated into ASP473, the mutant CDC25BC473D:CDK2/CyclinA system was less stable than the wild-type CDC25BWT:CDK2/CyclinA system. Finally, active site CYS473 of CDC25B was speculated to be the key residue, which had great effects on the binding between CDC25BCYS473 and CDK2 in the CDC25BC473D:CDK2/CyclinA system. Consequently, overall analyses appeared in this study ultimately provided a useful understanding of the weak interactions between CDC25BCYS473D and CDK2/CyclinA.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Li-Peng Li
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, People's Republic of China
| | - Hao-Xin Li
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, People's Republic of China
| | - Hui Zhou
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, People's Republic of China
| | - Wei-Ya Li
- China Department of Pharmacy, Tianjin Medical University, The First Hospital of Hebei Medical University, Shijiazhuang, China
| | - Run-Ling Wang
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, People's Republic of China
| | - Ying-Chi Zhang
- State Key Laboratory of Experimental Hematology, National Clinical Research Center for Blood Diseases, Institute of Hematology & Blood Diseases Hospital, Chinese Academy of Medical Sciences & Peking Union Medical College, Tianjin, China
| | - Ying Ma
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, People's Republic of China
| |
Collapse
|
4
|
Yu YX, Wang W, Sun HB, Zhang LL, Wang LF, Yin YY. Decoding drug resistant mechanism of V32I, I50V and I84V mutations of HIV-1 protease on amprenavir binding by using molecular dynamics simulations and MM-GBSA calculations. SAR AND QSAR IN ENVIRONMENTAL RESEARCH 2022; 33:805-831. [PMID: 36322686 DOI: 10.1080/1062936x.2022.2140708] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/08/2022] [Accepted: 10/20/2022] [Indexed: 06/16/2023]
Abstract
Mutations V32I, I50V and I84V in the HIV-1 protease (PR) induce drug resistance towards drug amprenavir (APV). Multiple short molecular dynamics (MSMD) simulations and molecular mechanics generalized Born surface area (MM-GBSA) method were utilized to investigate drug-resistant mechanism of V32I, I50V and I84V towards APV. Dynamic information arising from MSMD simulations suggest that V32I, I50V and I84V highly affect structural flexibility, motion modes and conformational behaviours of two flaps in the PR. Binding free energies calculated by MM-GBSA method suggest that the decrease in binding enthalpy and the increase in binding entropy induced by mutations V32I, I50V and I84V are responsible for drug resistance of the mutated PRs on APV. The energetic contributions of separate residues on binding of APV to the PR show that V32I, I50V and I84V highly disturb the interactions of two flaps with APV and mostly drive the decrease in binding ability of APV to the PR. Thus, the conformational changes of two flaps in the PR caused by V32I, I50V and I84V play key roles in drug resistance of three mutated PR towards APV. This study can provide useful dynamics information for the design of potent inhibitors relieving drug resistance.
Collapse
Affiliation(s)
- Y X Yu
- School of Science, Shandong Jiaotong University, Jinan, China
| | - W Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - H B Sun
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L L Zhang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - L F Wang
- School of Science, Shandong Jiaotong University, Jinan, China
| | - Y Y Yin
- School of Science, Shandong Jiaotong University, Jinan, China
| |
Collapse
|
5
|
Baltrukevich H, Podlewska S. From Data to Knowledge: Systematic Review of Tools for Automatic Analysis of Molecular Dynamics Output. Front Pharmacol 2022; 13:844293. [PMID: 35359865 PMCID: PMC8960308 DOI: 10.3389/fphar.2022.844293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2021] [Accepted: 01/26/2022] [Indexed: 12/02/2022] Open
Abstract
An increasing number of crystal structures available on one side, and the boost of computational power available for computer-aided drug design tasks on the other, have caused that the structure-based drug design tools are intensively used in the drug development pipelines. Docking and molecular dynamics simulations, key representatives of the structure-based approaches, provide detailed information about the potential interaction of a ligand with a target receptor. However, at the same time, they require a three-dimensional structure of a protein and a relatively high amount of computational resources. Nowadays, as both docking and molecular dynamics are much more extensively used, the amount of data output from these procedures is also growing. Therefore, there are also more and more approaches that facilitate the analysis and interpretation of the results of structure-based tools. In this review, we will comprehensively summarize approaches for handling molecular dynamics simulations output. It will cover both statistical and machine-learning-based tools, as well as various forms of depiction of molecular dynamics output.
Collapse
Affiliation(s)
- Hanna Baltrukevich
- Maj Institute of Pharmacology, Polish Academy of Sciences, Kraków, Poland
- Faculty of Pharmacy, Chair of Technology and Biotechnology of Medical Remedies, Jagiellonian University Medical College in Krakow, Kraków, Poland
| | - Sabina Podlewska
- Maj Institute of Pharmacology, Polish Academy of Sciences, Kraków, Poland
| |
Collapse
|
6
|
Ritonavir and xk263 Binding-Unbinding with HIV-1 Protease: Pathways, Energy and Comparison. LIFE (BASEL, SWITZERLAND) 2022; 12:life12010116. [PMID: 35054509 PMCID: PMC8779838 DOI: 10.3390/life12010116] [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/10/2021] [Revised: 01/04/2022] [Accepted: 01/10/2022] [Indexed: 01/22/2023]
Abstract
Understanding non-covalent biomolecular recognition, which includes drug-protein bound states and their binding/unbinding processes, is of fundamental importance in chemistry, biology, and medicine. Fully revealing the factors that govern the binding/unbinding processes can further assist in designing drugs with desired binding kinetics. HIV protease (HIVp) plays an integral role in the HIV life cycle, so it is a prime target for drug therapy. HIVp has flexible flaps, and the binding pocket can be accessible by a ligand via various pathways. Comparing ligand association and dissociation pathways can help elucidate the ligand-protein interactions such as key residues directly involved in the interaction or specific protein conformations that determine the binding of a ligand under certain pathway(s). Here, we investigated the ligand unbinding process for a slow binder, ritonavir, and a fast binder, xk263, by using unbiased all-atom accelerated molecular dynamics (aMD) simulation with a re-seeding approach and an explicit solvent model. Using ritonavir-HIVp and xk263-HIVp ligand-protein systems as cases, we sampled multiple unbinding pathways for each ligand and observed that the two ligands preferred the same unbinding route. However, ritonavir required a greater HIVp motion to dissociate as compared with xk263, which can leave the binding pocket with little conformational change of HIVp. We also observed that ritonavir unbinding pathways involved residues which are associated with drug resistance and are distal from catalytic site. Analyzing HIVp conformations sampled during both ligand-protein binding and unbinding processes revealed significantly more overlapping HIVp conformations for ritonavir-HIVp rather than xk263-HIVp. However, many HIVp conformations are unique in xk263-HIVp unbinding processes. The findings are consistent with previous findings that xk263 prefers an induced-fit model for binding and unbinding, whereas ritonavir favors a conformation selection model. This study deepens our understanding of the dynamic process of ligand unbinding and provides insights into ligand-protein recognition mechanisms and drug discovery.
Collapse
|
7
|
Li HX, Yang WY, Li LP, Zhou H, Li WY, Ma Y, Wang RL. Molecular dynamics study of CDC25B R492L mutant causing the activity decrease of CDC25B. J Mol Graph Model 2021; 109:108030. [PMID: 34509094 DOI: 10.1016/j.jmgm.2021.108030] [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/09/2021] [Revised: 09/01/2021] [Accepted: 09/05/2021] [Indexed: 11/25/2022]
Abstract
Cell division cycle 25B (CDC25B) was responsible for regulating the various stages of cell division in the cell cycle. R492L was one of the common types of CDC25B mutants. Researches showed that compared to CDC25BWT, CDC25BR492L mutant had a ∼100-fold reduction in the rate constant for forming phosphatase intermediate (k2). However, the molecular basis of how the CDC25BR492L mutant influenced the process of binding between CDC25B and CDK2/CyclinA was not yet known. Therefore, the optimizations of three-dimensional structure of the CDC25BWT-CDK2/CyclinA system and the CDC25BR492L-CDK2/CyclinA system were constructed by ZDOCK and RDOCK, and five methods were employed to verify the reasonability of the docking structure. Then the molecular dynamics simulations on the two systems were performed to explore the reason why CDC25BR492L mutant caused the weak interactions between CDC25BR492L and CDK2/CyclinA, respectively. The remote docking site (Arg488-Tyr497) and the second active site (Lys538-Arg544) of CDC25B were observed to have high fluctuations in the CDC25BR492L-CDK2/CyclinA system with post-analysis, where the high fluctuation of these two regions resulted in weak interactions between CD25B and CDK2. In addition, Asp38-Glu42 and Asp206-Asp210 of CDK2 showed the slightly descending fluctuation, and CDK2 revealed an enhanced the self-interaction, which made CDK2 keep a relatively stable state in the CDC25BR492L-CDK2/CyclinA system. Finally, Leu492 of CDC25B was speculated to be the key residue, which had great effects on the binding between CDC25BR492L and CDK2 in the CDC25BR492L-CDK2/CyclinA system. Consequently, overall analyses appeared in this study ultimately offered a helpful understanding of the weak interactions between CDC25BR492L and CDK2.
Collapse
Affiliation(s)
- Hao-Xin Li
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, 300070, People's Republic of China
| | - Wen-Yu Yang
- State Key Laboratory of Experimental Hematology, National Clinical Research Center for Blood Diseases, Institute of Hematology & Blood Diseases Hospital, Chinese Academy of Medical Sciences & Peking Union Medical College, Tianjin, 300020, People's Republic of China
| | - Li-Peng Li
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, 300070, People's Republic of China
| | - Hui Zhou
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, 300070, People's Republic of China
| | - Wei-Ya Li
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, 300070, People's Republic of China
| | - Ying Ma
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, 300070, People's Republic of China.
| | - Run-Ling Wang
- Tianjin Key Laboratory of Technologies Enabling Development of Clinical Therapeutics and Diagnostics, School of Pharmacy, Tianjin Medical University, Tianjin, 300070, People's Republic of China.
| |
Collapse
|
8
|
Vora J, Athar M, Sinha S, Jha PC, Shrivastava N. Binding Insight of Anti-HIV Phytocompounds with Prime Targets of HIV: A Molecular Dynamics Simulation Analysis. Curr HIV Res 2021; 18:132-141. [PMID: 31995010 DOI: 10.2174/1570162x18666200129112509] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2019] [Revised: 12/27/2019] [Accepted: 01/09/2020] [Indexed: 12/29/2022]
Abstract
BACKGROUND Despite intense efforts, AIDS is difficult to tackle by current anti-retroviral therapy (ART) due to its side effects; therefore, there is an urgent need to discover potential, multitarget and low-cost anti-HIV compounds. OBJECTIVE We have shown that few phytocompounds can potentially inhibit the prime targets of HIV namely GP120 envelope protein, reverse transcriptase, protease, integrase and ribonulcease. In this study, top ranked prioritized compounds were subjected to Molecular Dynamics (MD) simulation in order to study the conformational dynamics and integrity of crucial interaction in the receptor sites. METHODS The system was built for selected protein-ligand complex using TIP3P water model and OPLS_2005 force field. Trajectories were recorded up to 20 ns simulation time in Desmond module of Schrödinger software. RESULTS As a result of a comprehensive analysis of molecular properties and dynamics of the complexes, it has been concluded that Chebulic acid, Curcumin and Mulberroside C could be developed as envelope glycoprotein GP120 inhibitor, reverse transcriptase inhibitor and protease inhibitor respectively. However, the fluctuation of Chebulic acid with respect to integrase and ribonuclease protein was higher during the simulation. CONCLUSION These findings can aid in the designing of the structural properties for more effective anti-HIV compounds against the given targets.
Collapse
Affiliation(s)
- Jaykant Vora
- B.V. Patel Pharmaceutical Education and Research Development (PERD) Centre, Ahmedabad, Gujarat, India.,Department of Life Science, Gujarat University, Ahmedabad, Gujarat; India
| | - Mohd Athar
- School of Chemical Sciences, Central University of Gujarat, Gandhinagar, India
| | - Sonam Sinha
- B.V. Patel Pharmaceutical Education and Research Development (PERD) Centre, Ahmedabad, Gujarat, India.,Department of Life Science, Gujarat University, Ahmedabad, Gujarat; India
| | - Prakash C Jha
- Centre for Applied Chemistry, Central University of Gujarat, Gandhinagar, India
| | - Neeta Shrivastava
- B.V. Patel Pharmaceutical Education and Research Development (PERD) Centre, Ahmedabad, Gujarat, India
| |
Collapse
|
9
|
Eche S, Kumar A, Sonela N, Gordon ML. Acquired HIV-1 Protease Conformational Flexibility Associated with Lopinavir Failure May Shape the Outcome of Darunavir Therapy after Antiretroviral Therapy Switch. Biomolecules 2021; 11:489. [PMID: 33805099 PMCID: PMC8064090 DOI: 10.3390/biom11040489] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2021] [Revised: 03/08/2021] [Accepted: 03/12/2021] [Indexed: 11/16/2022] Open
Abstract
Understanding the underlying molecular interaction during a therapy switch from lopinavir (LPV) to darunavir (DRV) is essential to achieve long-term virological suppression. We investigated the kinetic and structural characteristics of multidrug-resistant South African HIV-1 subtype C protease (HIV-1 PR) during therapy switch from LPV to DRV using enzyme activity and inhibition assay, fluorescence spectroscopy, and molecular dynamic simulation. The HIV-1 protease variants were from clinical isolates with a combination of drug resistance mutations; MUT-1 (M46I, I54V, V82A, and L10F), MUT-2 (M46I, I54V, L76V, V82A, L10F, and L33F), and MUT-3 (M46I, I54V, L76V, V82A, L90M, and F53L). Enzyme kinetics analysis shows an association between increased relative resistance to LPV and DRV with the progressive decrease in the mutant HIV-1 PR variants' catalytic efficiency. A direct relationship between high-level resistance to LPV and intermediate resistance to DRV with intrinsic changes in the three-dimensional structure of the mutant HIV-1 PR as a function of the multidrug-resistance mutation was observed. In silico analysis attributed these structural adjustments to the multidrug-resistance mutations affecting the LPV and DRV binding landscape. Though DRV showed superiority to LPV, as a lower concentration was needed to inhibit the HIV-1 PR variants, the inherent structural changes resulting from mutations selected during LPV therapy may dynamically shape the DRV treatment outcome after the therapy switch.
Collapse
Affiliation(s)
- Simeon Eche
- Discipline of Virology, School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban 4001, South Africa;
| | - Ajit Kumar
- Discipline of Microbiology, School of Life Sciences, University of KwaZulu-Natal (Westville Campus), Durban 4000, South Africa;
| | - Nelson Sonela
- School of Medicine, Physical and Natural Sciences, University of Rome Tor Vegata, 1-00133 Rome, Italy;
- Chantal Biya International Reference Center for Research on the Management and Prevention of HIV/AIDS (CIRCB), Yaoundé P.O. Box 3077, Cameroon
| | - Michelle L. Gordon
- Discipline of Virology, School of Laboratory Medicine and Medical Sciences, University of KwaZulu-Natal, Durban 4001, South Africa;
| |
Collapse
|
10
|
Ngo ST, Hong ND, Quynh Anh LH, Hiep DM, Tung NT. Effective estimation of the inhibitor affinity of HIV-1 protease via a modified LIE approach. RSC Adv 2020; 10:7732-7739. [PMID: 35492181 PMCID: PMC9049864 DOI: 10.1039/c9ra09583g] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2019] [Accepted: 02/06/2020] [Indexed: 01/07/2023] Open
Abstract
The inhibition of the Human Immunodeficiency Virus Type 1 Protease (HIV-1 PR) can prevent the synthesis of new viruses. Computer-aided drug design (CADD) would enhance the discovery of new therapies, through which the estimation of ligand-binding affinity is critical to predict the most efficient inhibitor. A time-consuming binding free energy method would reduce the usefulness of CADD. The modified linear interaction energy (LIE) approach emerges as an appropriate protocol that performs this task. In particular, the polar interaction free energy, which is obtained via numerically resolving the linear Poisson-Boltzmann equation, plays as an important role in driving the binding mechanism of the HIV-1 PR + inhibitor complex. The electrostatic interaction energy contributes to the attraction between two molecules, but the vdW interaction acts as a repulsive factor between the ligand and the HIV-1 PR. Moreover, the ligands were found to adopt a very strong hydrophobic interaction with the HIV-1 PR. Furthermore, the results obtained corroborate the high accuracy and precision of computational studies with a large correlation coefficient value R = 0.83 and a small RMSE δ RMSE = 1.25 kcal mol-1. This method is less time-consuming than the other end-point methods, such as the molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) and free energy perturbation (FEP) approaches. Overall, the modified LIE approach would provide ligand-binding affinity with HIV-1 PR accurately, precisely, and rapidly, resulting in a more efficient design of new inhibitors.
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
| | - Nam Dao Hong
- University of Medicine and Pharmacy at Ho Chi Minh City Ho Chi Minh City Vietnam
| | - Le Huu Quynh Anh
- Department of Climate Change and Renewable Energy, Ho Chi Minh City University of Natural Resources and Environment Ho Chi Minh City Vietnam
| | | | - Nguyen Thanh Tung
- Institute of Materials Science & Graduate University of Science and Technology, Vietnam Academy of Science and Technology Hanoi Vietnam
| |
Collapse
|
11
|
Triki D, Kermarrec M, Visseaux B, Descamps D, Flatters D, Camproux AC, Regad L. Exploration of the effects of sequence variations between HIV-1 and HIV-2 proteases on their three-dimensional structures. J Biomol Struct Dyn 2019; 38:5014-5026. [PMID: 31830870 DOI: 10.1080/07391102.2019.1704877] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
HIV protease inhibitors (PIs) approved by the FDA (US Food and Drug Administration) are a major class of antiretroviral. HIV-2 protease (PR2) is naturally resistant to most of them as PIs were designed for HIV-1 protease (PR1). In this study, we explored the impact of amino-acid substitutions between PR1 and PR2 on the structure of protease (PR) by comparing the structural variability of 13 regions using 24 PR1 and PR2 structures complexed with diverse ligands. Our analyses confirmed structural rigidity of the catalytic region and highlighted the important role of three regions in the conservation of the catalytic region conformation. Surprisingly, we showed that the flap region, corresponding to a flexible region, exhibits similar conformations in PR1 and PR2. Furthermore, we identified regions exhibiting different conformations in PR1 and PR2, which could be explained by the intrinsic flexibility of these regions, by crystal packing, or by PR1 and PR2 substitutions. Some substitutions induce structural changes in the R2 and R4 regions that could have an impact on the properties of PI-binding site and could thus modify PI binding mode. Substitutions involved in structural changes in the elbow region could alter the flexibility of the PR2 flap regions relative to PR1, and thus play a role in the transition from the semi-open form to the closed form, and have an impact on ligand binding. These results improve the understanding of the impact of sequence variations between PR1 and PR2 on the natural resistance of HIV-2 to commercially available PIs.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Dhoha Triki
- Department of Chemistry, Bioinformatics Research Center, College of Science, North Carolina State University, Raleigh, North Carolina, USA.,Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, Paris, France
| | - Maxime Kermarrec
- Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, Paris, France
| | - Benoît Visseaux
- Université de Paris, IAME, UMR 1137, INSERM, AP-HP, Laboratoire de Virologie, Hôpital Bichat-Claude Bernard, Paris, France
| | - Diane Descamps
- Université de Paris, IAME, UMR 1137, INSERM, AP-HP, Laboratoire de Virologie, Hôpital Bichat-Claude Bernard, Paris, France
| | - Delphine Flatters
- Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, Paris, France
| | | | - Leslie Regad
- Université de Paris, BFA, UMR 8251, CNRS, ERL U1133, Inserm, Paris, France
| |
Collapse
|