1
|
Arı E, Şahin N, Üstün E, Dündar M, Karcı H, Özdemir İ, Koç A, Gürbüz N, Özdemir İ. Synthesis, antimicrobial activity and molecular docking study of benzyl functionalized benzimidazole silver(I) complexes. J Biol Inorg Chem 2023; 28:725-736. [PMID: 37934281 DOI: 10.1007/s00775-023-02024-y] [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: 05/17/2023] [Accepted: 09/10/2023] [Indexed: 11/08/2023]
Abstract
In this study, a series of N-functionalized benzimidazole silver(I) complexes were prepared and characterized by FT-IR, 1H, 13C{1H} NMR spectroscopy, and elemental analysis. Synthesized N-benzylbenzimidazole silver(I) complexes were evaluated for their antimicrobial activities against bacteria Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus, and the fungal strains Candida albicans and Candida glabrata. The results indicated that N-alkylbenzimidazole silver(I) complexes exhibited good antimicrobial activity compared to N-alkylbenzimidazole derivatives. Especially, complex 2e presented perfect antimicrobial activity than the other complexes. The characterized molecules were optimized by DFT-based calculation methods and the optimized molecules were analyzed in detail by molecular docking methods against bacterial DNA-gyrase and CYP51. The amino acid residues detected for both target molecules are consistent with expectations, and the calculated binding affinities and inhibition constants are promising for further studies. A series of N-alkylbenzimidazole silver(I) complexes were synthesized and fully characterized by means of 1H NMR, 13C NMR, and FT-IR spectroscopies. Synthesized N-alkylbenzimidazole silver(I) complexes were investigated for their antimicrobial activities against bacteria Escherichia coli, Pseudomonas aeruginosa, Staphylococcus aureus, and the fungal strains Candida albicans and Candida glabrata. All complexes showed better activity according to Ampicilin against Pseudomonas aeruginosa. The molecules which were firstly optimized by DFT-based calculation methods were also analyzed by molecular docking methods against DNA gyrase of E. Coli and CYP51. 338 × 190 mm (96 × 96 DPI).
Collapse
Affiliation(s)
- Erkan Arı
- Department of Chemistry, Faculty of Science and Art, İnönü University, 44280, Malatya, Turkey
| | - Neslihan Şahin
- Department of Science Education, Faculty of Education, Cumhuriyet University, 58040, Sivas, Turkey
| | - Elvan Üstün
- Department of Chemistry, Faculty of Art and Science, Ordu University, 52200, Ordu, Turkey
| | - Muhammed Dündar
- Department of Molecular Biology and Genetics, Faculty of Science and Art, İnönü University, Malatya, Turkey
- Drug Application and Research Center, İnönü University, 44280, Malatya, Turkey
| | - Hüseyin Karcı
- Department of Chemistry, Faculty of Science and Art, İnönü University, 44280, Malatya, Turkey
- Drug Application and Research Center, İnönü University, 44280, Malatya, Turkey
| | - İlknur Özdemir
- Department of Chemistry, Faculty of Science and Art, İnönü University, 44280, Malatya, Turkey.
- Drug Application and Research Center, İnönü University, 44280, Malatya, Turkey.
| | - Ahmet Koç
- Department of Genetics, Faculty of Medicine, İnönü University, Malatya, Turkey
| | - Nevin Gürbüz
- Department of Chemistry, Faculty of Science and Art, İnönü University, 44280, Malatya, Turkey
- Drug Application and Research Center, İnönü University, 44280, Malatya, Turkey
| | - İsmail Özdemir
- Department of Chemistry, Faculty of Science and Art, İnönü University, 44280, Malatya, Turkey
- Drug Application and Research Center, İnönü University, 44280, Malatya, Turkey
| |
Collapse
|
2
|
Singh K, Patil RB, Patel V, Remenyik J, Hegedűs T, Goda K. Synergistic Inhibitory Effect of Quercetin and Cyanidin-3O-Sophoroside on ABCB1. Int J Mol Sci 2023; 24:11341. [PMID: 37511101 PMCID: PMC10379049 DOI: 10.3390/ijms241411341] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/06/2023] [Revised: 06/30/2023] [Accepted: 07/05/2023] [Indexed: 07/30/2023] Open
Abstract
The human ABCB1 (P-glycoprotein, Pgp) protein is an active exporter expressed in the plasma membrane of cells forming biological barriers. In accordance with its broad substrate spectrum and tissue expression pattern, it affects the pharmacokinetics of numerous chemotherapeutic drugs and it is involved in unwanted drug-drug interactions leading to side effects or toxicities. When expressed in tumor tissues, it contributes to the development of chemotherapy resistance in malignancies. Therefore, the understanding of the molecular details of the ligand-ABCB1 interactions is of crucial importance. In a previous study, we found that quercetin (QUR) hampers both the transport and ATPase activity of ABCB1, while cyandin-3O-sophroside (C3S) stimulates the ATPase activity and causes only a weak inhibition of substrate transport. In the current study, when QUR and C3S were applied together, both a stronger ATPase inhibition and a robust decrease in substrate transport were observed, supporting their synergistic ABCB1 inhibitory effect. Similar to cyclosporine A, a potent ABCB1 inhibitor, co-treatment with QUR and C3S shifted the conformational equilibrium to the "inward-facing" conformer of ABCB1, as it was detected by the conformation-selective UIC2 mAb. To gain deeper insight into the molecular details of ligand-ABCB1 interactions, molecular docking experiments and MD simulations were also carried out. Our in silico studies support that QUR and C3S can bind simultaneously to ABCB1. The most favourable ligand-ABCB1 interaction is obtained when C3S binds to the central substrate binding site and QUR occupies the "access tunnel". Our results also highlight that the strong ABCB1 inhibitory effect of the combined treatment with QUR and C3S may be exploited in chemotherapy protocols for the treatment of multidrug-resistant tumors or for improving drug delivery through pharmacological barriers.
Collapse
Affiliation(s)
- Kuljeet Singh
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, 4032 Debrecen, Hungary
- Doctoral School of Molecular Cell and Immune Biology, University of Debrecen, 4032 Debrecen, Hungary
| | - Rajesh B Patil
- Department of Pharmaceutical Chemistry, Sinhgad Technical Education Society's Sinhgad College of Pharmacy, OffSinhgad Road, Vadgaon (Bk), Pune 411041, Maharashtra, India
| | - Vikas Patel
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, 4032 Debrecen, Hungary
| | - Judit Remenyik
- Institute of Food Technology, Faculty of Agricultural and Food Sciences and Environmental Management, University of Debrecen, 4032 Debrecen, Hungary
| | - Tamás Hegedűs
- Department of Biophysics and Radiation Biology, Semmelweis University, 1085 Budapest, Hungary
- ELKH-SE Biophysical Virology Research Group, Eötvös Loránd Research Network, 1052 Budapest, Hungary
| | - Katalin Goda
- Department of Biophysics and Cell Biology, Faculty of Medicine, University of Debrecen, 4032 Debrecen, Hungary
| |
Collapse
|
3
|
Mishra S, Rout M, Panda S, Singh SK, Sinha R, Dehury B, Pati S. An immunoinformatic approach towards development of a potent and effective multi-epitope vaccine against monkeypox virus (MPXV). J Biomol Struct Dyn 2023; 41:11714-11727. [PMID: 36591724 DOI: 10.1080/07391102.2022.2163426] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Accepted: 12/22/2022] [Indexed: 01/03/2023]
Abstract
Monkeypox is a viral zoonotic disease, often transmitted to humans from animals. While the whole world is haggling with the COVID-19 pandemic, the emergence of the monkeypox virus (MPXV) arose as a new challenge to mankind. Till date, numerous cases related to the MPXV have been reported in several countries across the globe, but, its momentary distribution in the current time has left everyone in fright with increasing mortality and limited clinically approved treatments. Therefore, it is of immense importance to develop a potent and highly effective vaccine capable of inducing desired immunogenic responses against the highly contagious MPXV. Herein, using various immunoinformatic and computational biology tools, we made an attempt to develop a multi-epitope vaccine construct against the MPXV which is antigenic, non-allergen and non-toxic in nature and capable of exhibiting immunogenic behavior. The sequence of vaccine construct was designed using the proposed 4 MHC-I, 3 MHC-II and 4 B-cell epitopes linked with suitable adjuvant and linkers. The modeled structure of the vaccine construct was used to assess its interaction with the Toll-like Receptor 4 (TLR4) using ClusPro and HADDOCK. All-atoms molecular dynamics simulation of the MPXV vaccine construct-TLR4 complex followed by a high level of gene expression of the construct within the bacterial system affirmed its stability along with induction of immunogenic response within the host cell. Altogether, our immunoinformatic approach aid in the development of a stable chimeric vaccine construct against MPXV and needs further experimental validation for its immunological relevance and usefulness as a vaccine candidate.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Sarbani Mishra
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| | - Madhusmita Rout
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| | - Sunita Panda
- Mycology Laboratory, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| | - Susheel Kumar Singh
- Vaccine and Diagnostic Laboratory, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| | - Rohan Sinha
- Department of Computer Science and Engineering, National Institute of Technology Patna, Patna, Bihar, India
| | - Budheswar Dehury
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| | - Sanghamitra Pati
- Bioinformatics Division, ICMR-Regional Medical Research Centre, Bhubaneswar, Odisha, India
| |
Collapse
|
4
|
Höing A, Kirupakaran A, Beuck C, Pörschke M, Niemeyer FC, Seiler T, Hartmann L, Bayer P, Schrader T, Knauer SK. Recognition of a Flexible Protein Loop in Taspase 1 by Multivalent Supramolecular Tweezers. Biomacromolecules 2022; 23:4504-4518. [PMID: 36200481 DOI: 10.1021/acs.biomac.2c00652] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Many natural proteins contain flexible loops utilizing well-defined complementary surface regions of their interacting partners and usually undergo major structural rearrangements to allow perfect binding. The molecular recognition of such flexible structures is still highly challenging due to the inherent conformational dynamics. Notably, protein-protein interactions are on the other hand characterized by a multivalent display of complementary binding partners to enhance molecular affinity and specificity. Imitating this natural concept, we here report the rational design of advanced multivalent supramolecular tweezers that allow addressing two lysine and arginine clusters on a flexible protein surface loop. The protease Taspase 1, which is involved in cancer development, carries a basic bipartite nuclear localization signal (NLS) and thus interacts with Importin α, a prerequisite for proteolytic activation. Newly established synthesis routes enabled us to covalently fuse several tweezer molecules into multivalent NLS ligands. The resulting bi- up to pentavalent constructs were then systematically compared in comprehensive biochemical assays. In this series, the stepwise increase in valency was robustly reflected by the ligands' gradually enhanced potency to disrupt the interaction of Taspase 1 with Importin α, correlated with both higher binding affinity and inhibition of proteolytic activity.
Collapse
Affiliation(s)
- Alexander Höing
- Molecular Biology II, Center of Medical Biotechnology (ZMB), University of Duisburg-Essen, Universitätsstrasse 5, 45141 Essen, Germany
| | - Abbna Kirupakaran
- Institute of Organic Chemistry I, University of Duisburg-Essen, Universitätsstrasse 7, 45141 Essen, Germany
| | - Christine Beuck
- Structural and Medicinal Biochemistry, Center of Medical Biotechnology (ZMB), University of Duisburg-Essen, Universitätsstrasse 5, 45141 Essen, Germany
| | - Marius Pörschke
- Structural and Medicinal Biochemistry, Center of Medical Biotechnology (ZMB), University of Duisburg-Essen, Universitätsstrasse 5, 45141 Essen, Germany
| | - Felix C Niemeyer
- Institute of Organic Chemistry I, University of Duisburg-Essen, Universitätsstrasse 7, 45141 Essen, Germany
| | - Theresa Seiler
- Organic Chemistry and Macromolecular Chemistry, Heinrich Heine University Düsseldorf, Universitätsstrasse 1, 40225 Düsseldorf, Germany
| | - Laura Hartmann
- Organic Chemistry and Macromolecular Chemistry, Heinrich Heine University Düsseldorf, Universitätsstrasse 1, 40225 Düsseldorf, Germany
| | - Peter Bayer
- Structural and Medicinal Biochemistry, Center of Medical Biotechnology (ZMB), University of Duisburg-Essen, Universitätsstrasse 5, 45141 Essen, Germany
| | - Thomas Schrader
- Institute of Organic Chemistry I, University of Duisburg-Essen, Universitätsstrasse 7, 45141 Essen, Germany
| | - Shirley K Knauer
- Molecular Biology II, Center of Medical Biotechnology (ZMB), University of Duisburg-Essen, Universitätsstrasse 5, 45141 Essen, Germany
| |
Collapse
|
5
|
Chaudhari AM, Joshi M, Kumar D, Patel A, Lokhande KB, Krishnan A, Hanack K, Filipek S, Liepmann D, Renugopalakrishnan V, Paulmurugan R, Joshi C. Evaluation of immune evasion in SARS-CoV-2 Delta and Omicron variants. Comput Struct Biotechnol J 2022; 20:4501-4516. [PMID: 35965661 PMCID: PMC9359593 DOI: 10.1016/j.csbj.2022.08.010] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2022] [Revised: 08/04/2022] [Accepted: 08/04/2022] [Indexed: 12/18/2022] Open
Abstract
Emerging SARS-CoV-2 variants with higher transmissibility and immune escape remain a persistent threat across the globe. This is evident from the recent outbreaks of the Delta (B.1.617.2) and Omicron variants. These variants have originated from different continents and spread across the globe. In this study, we explored the genomic and structural basis of these variants for their lineage defining mutations of the spike protein through computational analysis, protein modeling, and molecular dynamic (MD) simulations. We further experimentally validated the importance of these deletion mutants for their immune escape using a pseudovirus-based neutralization assay, and an antibody (4A8) that binds directly to the spike protein's NTD. Delta variant with the deletion and mutations in the NTD revealed a better rigidity and reduced flexibility as compared to the wild-type spike protein (Wuhan isolate). Furthermore, computational studies of 4A8 monoclonal antibody (mAb) revealed a reduced binding of Delta variant compared to the wild-type strain. Similarly, the MD simulation data and virus neutralization assays revealed that the Omicron also exhibits immune escape, as antigenic beta-sheets appear to be disrupted. The results of the present study demonstrate the higher possibility of immune escape and thereby achieved better fitness advantages by the Delta and Omicron variants, which warrants further demonstrations through experimental evidences. Our study, based on in-silico computational modelling, simulations, and pseudovirus-based neutralization assay, highlighted and identified the probable mechanism through which the Delta and Omicron variants are more pathogenically evolved with higher transmissibility as compared to the wild-type strain.
Collapse
Affiliation(s)
- Armi M Chaudhari
- Gujarat Biotechnology Research Centre (GBRC), Department of Science and Technology, Government of Gujarat, Gandhinagar 382011, India
| | - Madhvi Joshi
- Gujarat Biotechnology Research Centre (GBRC), Department of Science and Technology, Government of Gujarat, Gandhinagar 382011, India
| | - Dinesh Kumar
- Gujarat Biotechnology Research Centre (GBRC), Department of Science and Technology, Government of Gujarat, Gandhinagar 382011, India
| | - Amrutlal Patel
- Gujarat Biotechnology Research Centre (GBRC), Department of Science and Technology, Government of Gujarat, Gandhinagar 382011, India
| | - Kiran Bharat Lokhande
- Gujarat Biotechnology Research Centre (GBRC), Department of Science and Technology, Government of Gujarat, Gandhinagar 382011, India
| | - Anandi Krishnan
- Cellular Pathway Imaging Laboratory (CPIL), Department of Radiology, Stanford University School of Medicine, Palo Alto, CA 94304, United States
- Department of Pathology, Stanford University School of Medicine, Palo Alto, CA 94304, United States
| | - Katja Hanack
- Immunotechnology Group, Department of Biochemistry and Biology, University of Potsdam, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany
| | - Slawomir Filipek
- Faculty of Chemistry & Biological and Chemical Research, Centre, University of Warsaw, ul, Pasteura 1, 02-093 Warsaw, Poland
| | - Dorian Liepmann
- Department of Bioengineering, University of California, Berkeley, Berkeley, CA 94720, United States
| | - Venkatesan Renugopalakrishnan
- Department of Chemistry, Northeastern University, Boston Children's Hospital, Harvard Medical School, Boston, MGB Center for COVID Innovation, MA 02115, United States
| | - Ramasamy Paulmurugan
- Cellular Pathway Imaging Laboratory (CPIL), Department of Radiology, Stanford University School of Medicine, Palo Alto, CA 94304, United States
| | - Chaitanya Joshi
- Gujarat Biotechnology Research Centre (GBRC), Department of Science and Technology, Government of Gujarat, Gandhinagar 382011, India
| |
Collapse
|
6
|
Wang Y, Tang H, Gao C, Ge M, Li Z, Dong Z, Zhao L. Flexibility-aware graph model for accurate epitope identification. Comput Biol Med 2022; 149:106064. [DOI: 10.1016/j.compbiomed.2022.106064] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2022] [Revised: 08/05/2022] [Accepted: 08/27/2022] [Indexed: 11/25/2022]
|
7
|
Santa-Coloma TA. Overlapping synthetic peptides as a tool to map protein-protein interactions ̶ FSH as a model system of nonadditive interactions. Biochim Biophys Acta Gen Subj 2022; 1866:130153. [DOI: 10.1016/j.bbagen.2022.130153] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2022] [Revised: 04/06/2022] [Accepted: 04/12/2022] [Indexed: 10/18/2022]
|
8
|
Nguyen PH, Derreumaux P. Computer Simulations Aimed at Exploring Protein Aggregation and Dissociation. Methods Mol Biol 2022; 2340:175-196. [PMID: 35167075 DOI: 10.1007/978-1-0716-1546-1_9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Protein aggregation can lead to well-defined structures that are functional, but is also the cause of the death of neuron cells in many neurodegenerative diseases. The complexity of the molecular events involved in the aggregation kinetics of amyloid proteins and the transient and heterogeneous characters of all oligomers prevent high-resolution structural experiments. As a result, computer simulations have been used to determine the atomic structures of amyloid proteins at different association stages as well as to understand fibril dissociation. In this chapter, we first review the current computer simulation methods used for aggregation with some atomistic and coarse-grained results aimed at better characterizing the early formed oligomers and amyloid fibril formation. Then we present the applications of non-equilibrium molecular dynamics simulations to comprehend the dissociation of protein assemblies.
Collapse
Affiliation(s)
- Phuong H Nguyen
- Laboratoire de Biochimie Théorique, UPR 9080, CNRS, Université de Paris, Paris, France
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, UPR 9080, CNRS, Université de Paris, Paris, France.
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris, France.
| |
Collapse
|
9
|
Siddiqui SA, Dubey KD. Can the local electric field be a descriptor of catalytic activity? A case study on chorismate mutase. Phys Chem Chem Phys 2021; 24:1974-1981. [PMID: 34757367 DOI: 10.1039/d1cp03978d] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
The current theoretical perception of enzymatic activity is highly reliant on the determination of the activation energy of the reactions, which is often calculated using computationally demanding quantum mechanical calculations. With the ever-increasing use of bioengineering techniques that produce too many variants of the same enzyme, a fast and accurate way to study the relative efficiency of enzymes is currently in high demand. Here, we propose the local electric field (LEF) of the enzyme along the reaction axis as a descriptor for the enzymatic activity using the example of chorismate mutase in its native form and several variants (R90A, R90G, and R90K/C88S). The study shows a direct correlation between the calculated enzymatic EF and the enzymatic activity for all the complexes. MD simulations of the Michaelis complex and the transition state analog (TSA) show a stabilizing force on the TSA due to the enzymatic EF. QM/MM and QM-only DFT calculations in the presence of an external electric field (EEF) oriented along the reaction axis show that the electric field can interact with the dipole moment of the TS, thereby stabilizing it and thus lowering the activation energy.
Collapse
Affiliation(s)
- Shakir Ali Siddiqui
- Department of Chemistry, School of Natural Sciences, Shiv Nadar University, Delhi-NCR, 201314, India.
| | - Kshatresh Dutta Dubey
- Department of Chemistry, School of Natural Sciences, Shiv Nadar University, Delhi-NCR, 201314, India. .,Center for Informatics, Department of Chemistry, School of Natural Sciences, Shiv Nadar University, Delhi-NCR, 201314, India
| |
Collapse
|
10
|
Rants'o TA, Johan van der Westhuizen C, van Zyl RL. Optimization of covalent docking for organophosphates interaction with Anopheles acetylcholinesterase. J Mol Graph Model 2021; 110:108054. [PMID: 34688161 DOI: 10.1016/j.jmgm.2021.108054] [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: 08/20/2021] [Revised: 10/13/2021] [Accepted: 10/14/2021] [Indexed: 11/29/2022]
Abstract
Organophosphates (OPs) used as potent insecticides for malaria vector control, covalently phosphorylate the catalytic serine residue of Anopheles gambiae AChE (AgAChE) in a reaction that liberates their leaving groups. In the recent 10-year insecticide use assessment, OPs were the most frequently used World Health Organization prequalified insecticides. Molecular modelling programs are best suited to display molecular interactions between ligands and the target proteins. The docking modes that generate ligand poses closer to the binding site show high accuracy in predicting the ligand binding mode. The implicit solvation approach such as molecular mechanics-generalized born surface area (MM-GBSA) is a more reliable method to predict ligand onformations and binding affinities. Apart from covalent docking studies being scarce, current molecular docking programs do not adequately possess the covalent docking reaction algorithm to display the molecular mechanism of OPs at the AgAChE catalytic site. This results into OP docking studies commonly being conducted through noncovalent pannels. The aim of this study was to establish the optimim covalent docking system for OPs through manual customization of Schrödinger's Glide covalent docking reaction algorithm. To achieve this, a newly customized covalent reaction algorithm was assessed on a set of ligands covering aromatic, non-aromatic and hydrophobic OPs and compared to the noncovalent docking results in terms of reliability based on the reported X-ray diffraction molecular interactions and crystal poses. The study established that by virtue of omitting the well-known OP hydrolysis, noncovalent mode suggested molecular interactions that were further from the catalytic triad and could not otherwise occur when the molecule is hydrolyzed as in the customized covalent docking mode. Moreover, the MM-GBSA concurred with the optimized covalent docking in eliminating such inaccurate molecular interactions. Additionally, the covalent docking mode confined the interactions and ligand poses to the catalytic site indicating relatively high accuracy and reliability. This study reports the optimized covalent docking panel that effectively confirmed the molecular mechanisms of OPs, as well as indentifying the corresponding amino acid residues required to stabilize the aromatic, non-aromatic and hydrophobic OPs at the AgAChE catalytic site in line with the reported X-ray diffraction studies. As such, the proposed manual customization of the Schrödinger's Glide covalent docking platform can be used to reliably predict molecular interactions between OPs and AgAChE target.
Collapse
Affiliation(s)
- Thankhoe A Rants'o
- Pharmacology Division, Department of Pharmacy and Pharmacology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, 2193, South Africa; WITS Research Institute for Malaria (WRIM), Faculty of Health Sciences, University of Witwatersrand, Johannesburg, 2193, South Africa.
| | - C Johan van der Westhuizen
- Council for Scientific and Industrial Research (CSIR), Future Production: Chemicals Cluster, Meiring Naude Road, Brummeria, Pretoria, 0001, South Africa
| | - Robyn L van Zyl
- Pharmacology Division, Department of Pharmacy and Pharmacology, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, 2193, South Africa; WITS Research Institute for Malaria (WRIM), Faculty of Health Sciences, University of Witwatersrand, Johannesburg, 2193, South Africa
| |
Collapse
|
11
|
Mutations of SARS-CoV-2 RBD May Alter Its Molecular Structure to Improve Its Infection Efficiency. Biomolecules 2021; 11:biom11091273. [PMID: 34572486 PMCID: PMC8466379 DOI: 10.3390/biom11091273] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2021] [Revised: 07/23/2021] [Accepted: 08/14/2021] [Indexed: 01/08/2023] Open
Abstract
The receptor-binding domain (RBD) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mediates the viral-host interaction and is a target for most neutralizing antibodies. Nevertheless, SARS-CoV-2 RBD mutations pose a threat due to their role in host cell entry via the human angiotensin-converting enzyme 2 receptor that might strengthen SARS-CoV-2 infectivity, viral load, or resistance against neutralizing antibodies. To understand the molecular structural link between RBD mutations and infectivity, the top five mutant RBDs (i.e., N501Y, E484K L452R, S477N, and N439K) were selected based on their recorded case numbers. These mutants along with wild-type (WT) RBD were studied through all-atom molecular dynamics (MD) simulations of 100 ns. The principal component analysis and the free energy landscape were used too. Interestingly, N501Y, N439K, and E484K mutations were observed to increase the rigidity in some RBD regions while increasing the flexibility of the receptor-binding motif (RBM) region, suggesting a compensation of the entropy penalty. However, S477N and L452R RBDs were observed to increase the flexibility of the RBM region while maintaining similar flexibility in other RBD regions in comparison to WT RBD. Therefore, both mutations (especially S477N) might destabilize the RBD structure, as loose conformation compactness was observed. The destabilizing effect of S477N RBD was consistent with previous work on S477N mutation. Finally, the free energy landscape results showed that mutations changed WT RBD conformation while local minima were maintained for all mutant RBDs. In conclusion, RBD mutations definitely impact the WT RBD structure and conformation as well as increase the binding affinity to angiotensin-converting enzyme receptor.
Collapse
|
12
|
Pitard I, Monet D, Goossens PL, Blondel A, Malliavin TE. Analyzing In Silico the Relationship Between the Activation of the Edema Factor and Its Interaction With Calmodulin. Front Mol Biosci 2020; 7:586544. [PMID: 33344505 PMCID: PMC7746812 DOI: 10.3389/fmolb.2020.586544] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Accepted: 11/02/2020] [Indexed: 11/25/2022] Open
Abstract
Molecular dynamics (MD) simulations have been recorded on the complex between the edema factor (EF) of Bacilllus anthracis and calmodulin (CaM), starting from a structure with the orthosteric inhibitor adefovir bound in the EF catalytic site. The starting structure has been destabilized by alternately suppressing different co-factors, such as adefovir ligand or ions, revealing several long-distance correlations between the conformation of CaM, the geometry of the CaM/EF interface, the enzymatic site and the overall organization of the complex. An allosteric communication between CaM/EF interface and the EF catalytic site, highlighted by these correlations, was confirmed by several bioinformatics approaches from the literature. A network of hydrogen bonds and stacking interactions extending from the helix V of of CaM, and the residues of the switches A, B and C, and connecting to catalytic site residues, is a plausible candidate for the mediation of allosteric communication. The greatest variability in volume between the different MD conditions was also found for cavities present at the EF/CaM interface and in the EF catalytic site. The similarity between the predictions from literature and the volume variability might introduce the volume variability as new descriptor of allostery.
Collapse
Affiliation(s)
- Irène Pitard
- Unité de Bioinformatique Structurale, Institut Pasteur and CNRS UMR 3528, Paris, France.,Center de Bioinformatique, Biostatistique et Biologie Intégrative, Institut Pasteur and CNRS USR 3756, Paris, France.,Ecole Doctorale Université Paris Sorbonne, Paris, France
| | - Damien Monet
- Unité de Bioinformatique Structurale, Institut Pasteur and CNRS UMR 3528, Paris, France.,Center de Bioinformatique, Biostatistique et Biologie Intégrative, Institut Pasteur and CNRS USR 3756, Paris, France.,Ecole Doctorale Université Paris Sorbonne, Paris, France
| | | | - Arnaud Blondel
- Unité de Bioinformatique Structurale, Institut Pasteur and CNRS UMR 3528, Paris, France.,Center de Bioinformatique, Biostatistique et Biologie Intégrative, Institut Pasteur and CNRS USR 3756, Paris, France
| | - Thérèse E Malliavin
- Unité de Bioinformatique Structurale, Institut Pasteur and CNRS UMR 3528, Paris, France.,Center de Bioinformatique, Biostatistique et Biologie Intégrative, Institut Pasteur and CNRS USR 3756, Paris, France
| |
Collapse
|
13
|
Decherchi S, Cavalli A. Thermodynamics and Kinetics of Drug-Target Binding by Molecular Simulation. Chem Rev 2020; 120:12788-12833. [PMID: 33006893 PMCID: PMC8011912 DOI: 10.1021/acs.chemrev.0c00534] [Citation(s) in RCA: 95] [Impact Index Per Article: 23.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2020] [Indexed: 12/19/2022]
Abstract
Computational studies play an increasingly important role in chemistry and biophysics, mainly thanks to improvements in hardware and algorithms. In drug discovery and development, computational studies can reduce the costs and risks of bringing a new medicine to market. Computational simulations are mainly used to optimize promising new compounds by estimating their binding affinity to proteins. This is challenging due to the complexity of the simulated system. To assess the present and future value of simulation for drug discovery, we review key applications of advanced methods for sampling complex free-energy landscapes at near nonergodicity conditions and for estimating the rate coefficients of very slow processes of pharmacological interest. We outline the statistical mechanics and computational background behind this research, including methods such as steered molecular dynamics and metadynamics. We review recent applications to pharmacology and drug discovery and discuss possible guidelines for the practitioner. Recent trends in machine learning are also briefly discussed. Thanks to the rapid development of methods for characterizing and quantifying rare events, simulation's role in drug discovery is likely to expand, making it a valuable complement to experimental and clinical approaches.
Collapse
Affiliation(s)
- Sergio Decherchi
- Computational
and Chemical Biology, Fondazione Istituto
Italiano di Tecnologia, 16163 Genoa, Italy
| | - Andrea Cavalli
- Computational
and Chemical Biology, Fondazione Istituto
Italiano di Tecnologia, 16163 Genoa, Italy
- Department
of Pharmacy and Biotechnology, University
of Bologna, 40126 Bologna, Italy
| |
Collapse
|
14
|
Burevschi E, Alonso ER, Sanz ME. Binding Site Switch by Dispersion Interactions: Rotational Signatures of Fenchone-Phenol and Fenchone-Benzene Complexes. Chemistry 2020; 26:11327-11333. [PMID: 32428270 PMCID: PMC7497235 DOI: 10.1002/chem.202001713] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2020] [Indexed: 12/15/2022]
Abstract
Non-covalent interactions between molecules determine molecular recognition and the outcome of chemical and biological processes. Characterising how non-covalent interactions influence binding preferences is of crucial importance in advancing our understanding of these events. Here, we analyse the interactions involved in smell and specifically the effect of changing the balance between hydrogen-bonding and dispersion interactions by examining the complexes of the common odorant fenchone with phenol and benzene, mimics of tyrosine and phenylalanine residues, respectively. Using rotational spectroscopy and quantum chemistry, two isomers of each complex have been identified. Our results show that the increased weight of dispersion interactions in these complexes changes the preferred binding site in fenchone and sets the basis for a better understanding of the effect of different residues in molecular recognition and binding events.
Collapse
|
15
|
Dehury B, Raina V, Misra N, Suar M. Effect of mutation on structure, function and dynamics of receptor binding domain of human SARS-CoV-2 with host cell receptor ACE2: a molecular dynamics simulations study. J Biomol Struct Dyn 2020; 39:7231-7245. [PMID: 32762417 PMCID: PMC7484587 DOI: 10.1080/07391102.2020.1802348] [Citation(s) in RCA: 31] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Recent studies have pointed the role of angiotensin-converting enzyme-II (ACE2) in mediating the entry of SARS-CoV-2 to the host cell by binding to the receptor-binding domain (RBD) of viral spike protein, and successive priming by cellular proteases initiates the infection. SARS-CoV replication rate and disease severity is controlled by the binding affinity of RBD with ACE2. To understand, how mutations in the conserved residues of RBD affect the molecular interaction with ACE2, we generated five alanine mutants i.e. Y449A, N487A, Y489A, N501A and Y505A in the receptor binding motif (RBM) of the ACE2-RBD SARS-CoV-2 complex (PDB: 6M0J). Computational site directed mutagenesis induced dynamics in wild-type and mutant complexes were extensively studied through all-atoms molecular dynamics (MD) simulations of 150 ns. In silico mutational analysis revealed loss of important intermolecular hydrogen bonds and other non-bonded contacts, critical for molecular recognition of SARS-CoV-2 RBD to ACE2, which is well supported by saturation mutagenesis study of binding interface residues. MD simulations results showed that RBM motif is flexible, where mutant residues are relatively more mobile than corresponding wild-type residues. Global motion analysis through principal component studies revealed that RBD exhibits protuberant in-ward motion towards the human ACE2 binding interface which may be crucial for molecular interaction. Conclusively, the present finding are in congruence with previous experimental reports and provides detailed information on the structural basis of receptor binding by human SARS-CoV-2, which will crucial for the development of novel inhibitors or drugs to combat against SARS-CoV-2. Communicated by Ramaswamy H. Sarma
Collapse
Affiliation(s)
- Budheswar Dehury
- School of Biotechnology, Kalinga Institute of Industrial Technology (KIIT), Deemed to be University, Bhubaneswar, India
| | - Vishakha Raina
- School of Biotechnology, Kalinga Institute of Industrial Technology (KIIT), Deemed to be University, Bhubaneswar, India
| | - Namrata Misra
- School of Biotechnology, Kalinga Institute of Industrial Technology (KIIT), Deemed to be University, Bhubaneswar, India.,KIIT-Technology Business Incubator (KIIT-TBI), Kalinga Institute of Industrial Technology (KIIT), Deemed to be University, Bhubaneswar, India
| | - Mrutyunjay Suar
- School of Biotechnology, Kalinga Institute of Industrial Technology (KIIT), Deemed to be University, Bhubaneswar, India.,KIIT-Technology Business Incubator (KIIT-TBI), Kalinga Institute of Industrial Technology (KIIT), Deemed to be University, Bhubaneswar, India
| |
Collapse
|
16
|
Structures of the intrinsically disordered Aβ, tau and α-synuclein proteins in aqueous solution from computer simulations. Biophys Chem 2020; 264:106421. [PMID: 32623047 DOI: 10.1016/j.bpc.2020.106421] [Citation(s) in RCA: 82] [Impact Index Per Article: 20.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2020] [Revised: 06/17/2020] [Accepted: 06/18/2020] [Indexed: 12/21/2022]
Abstract
Intrinsically disordered proteins (IDPs) play many biological roles in the human proteome ranging from vesicular transport, signal transduction to neurodegenerative diseases. The Aβ and tau proteins, and the α-synuclein protein, key players in Alzheimer's and Parkinson's diseases, respectively are fully disordered at the monomer level. The structural heterogeneity of the monomeric and oligomeric states and the high self-assembly propensity of these three IDPs have precluded experimental structural determination. Simulations have been used to determine the atomic structures of these IDPs. In this article, we review recent computer models to capture the equilibrium ensemble of Aβ, tau and α-synuclein proteins at different association steps in aqueous solution and present new results of the PEP-FOLD framework on α-synuclein monomer.
Collapse
|
17
|
Derreumaux P, Man VH, Wang J, Nguyen PH. Tau R3-R4 Domain Dimer of the Wild Type and Phosphorylated Ser356 Sequences. I. In Solution by Atomistic Simulations. J Phys Chem B 2020; 124:2975-2983. [PMID: 32216358 DOI: 10.1021/acs.jpcb.0c00574] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
In Alzheimer's disease, neurofibrillary lesions correlate with cognitive deficits and consist of inclusions of tau protein with cross-β structure. A stable dimeric form of soluble tau has been evidenced in the cells, but its high-resolution structure is missing in solution. We know, however, that cryo-electron microscopy (c-EM) of full-length tau in the brain of an individual with AD displays a core of eight β-sheets with a C-shaped architecture spanning the R3-R4 repeat domain, while the rest of the protein is very flexible. To address the conformational ensemble of the dimer, we performed atomistic replica exchange molecular dynamics simulations on the tau R3-R4 domain starting from the c-EM configuration. We find that the wild type tau R3-R4 dimer explores elongated, U-shaped, V-shaped, and globular forms rather than the C-shape. Phosphorylation of Ser356, pSer356, is known to block the interaction between the tau protein and the amyloid-β42 peptide. Standard molecular dynamics simulations of this phosphorylated sequence for a total of 5 μs compared to its wild type counterpart show a modulation of the population of β-helices and accessible topologies and a decrease of intermediates near the fibril-like conformers.
Collapse
Affiliation(s)
- Philippe Derreumaux
- Laboratory of Theoretical Chemistry, Ton Duc Thang University, 33000, Ho Chi Minh City, Vietnam.,Faculty of Pharmacy, Ton Duc Thang University, 33000, Ho Chi Minh City, Vietnam
| | - Viet Hoang Man
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Junmei Wang
- Department of Pharmaceutical Sciences, School of Pharmacy, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Phuong H Nguyen
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, F-75005, Paris, France.,Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, 75000, Paris, France
| |
Collapse
|
18
|
Minh Hung H, Nguyen MT, Tran PT, Truong VK, Chapman J, Quynh Anh LH, Derreumaux P, Vu VV, Ngo ST. Impact of the Astaxanthin, Betanin, and EGCG Compounds on Small Oligomers of Amyloid Aβ 40 Peptide. J Chem Inf Model 2020; 60:1399-1408. [PMID: 32105466 DOI: 10.1021/acs.jcim.9b01074] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
There is experimental evidence that the astaxanthin, betanin, and epigallocatechin-3-gallate (EGCG) compounds slow down the aggregation kinetics and the toxicity of the amyloid-β (Aβ) peptide. How these inhibitors affect the self-assembly at the atomic level remains elusive. To address this issue, we have performed for each ligand atomistic replica exchange molecular dynamic (REMD) simulations in an explicit solvent of the Aβ11-40 trimer from the U-shape conformation and MD simulations starting from Aβ1-40 dimer and tetramer structures characterized by different intra- and interpeptide conformations. We find that the three ligands have similar binding free energies on small Aβ40 oligomers but very distinct transient binding sites that will affect the aggregation of larger assemblies and fibril elongation of the Aβ40 peptide.
Collapse
Affiliation(s)
- Huynh Minh Hung
- Department of Chemistry, KU Leuven, Celestijnenlaan 200F, B-3001 Leuven, Belgium
| | - Minh Tho Nguyen
- Computational Chemistry Research Group, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
| | - Phuong-Thao Tran
- Department of Pharmaceutical Chemistry, Hanoi University of Pharmacy, Hanoi 100000, Vietnam
| | - Vi Khanh Truong
- School of Science, RMIT University, GPO Box 2476, Melbourne 3001, Australia
| | - James Chapman
- School of Science, RMIT University, GPO Box 2476, Melbourne 3001, Australia
| | - 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 700000, Vietnam
| | - Philippe Derreumaux
- Laboratory of Theoretical Chemistry, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam.,Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam.,Laboratoire de Biochimie Théorique, UPR9080, CNRS, Université de Paris, 13 rue Pierre et Marie Curie, F-75005 Paris, France.,Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, 75005 Paris, France
| | - Van V Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University, Ho Chi Minh City 700000, Vietnam
| | - Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam.,Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City 700000, Vietnam
| |
Collapse
|
19
|
Ngo ST, Nguyen PH, Derreumaux P. Impact of A2T and D23N Mutations on Tetrameric Aβ42 Barrel within a Dipalmitoylphosphatidylcholine Lipid Bilayer Membrane by Replica Exchange Molecular Dynamics. J Phys Chem B 2020; 124:1175-1182. [DOI: 10.1021/acs.jpcb.9b11881] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Son Tung Ngo
- Laboratory of Theoretical and Computational Biophysics, Ton Duc Thang University, Ho Chi Minh City 33000, Vietnam
- Faculty of Applied Sciences, Ton Duc Thang University, Ho Chi Minh City 33000, Vietnam
| | - Phuong H. Nguyen
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, 13 rue Pierre et Marie Curie, F-75005 Paris, France
- Institut de Biologie Physico-Chimique, Fondation Edmond de Rothschild, PSL Research University, Paris 75005, France
| | - Philippe Derreumaux
- Laboratory of Theoretical Chemistry, Ton Duc Thang University, Ho Chi Minh City 33000, Vietnam
- Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City 33000, Vietnam
| |
Collapse
|
20
|
Ngo ST, Nguyen PH, Derreumaux P. Stability of Aβ11-40 Trimers with Parallel and Antiparallel β-Sheet Organizations in a Membrane-Mimicking Environment by Replica Exchange Molecular Dynamics Simulation. J Phys Chem B 2020; 124:617-626. [PMID: 31931566 DOI: 10.1021/acs.jpcb.9b10982] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
The aggregation of the amyloid (Aβ) peptide of 39-43 amino acids into plaques is observed in the brain of Alzheimer's disease (AD) patients, but the mechanisms underlying the neurotoxicity of Aβ oligomers are still elusive. One suggested initial mechanism is related to the implications of amyloid membrane interactions, but characterization of these assemblies is challenging by experimental means. In this study, we have explored the stability of a trimer of Aβ11-40 in parallel and antiparallel β-sheet structures for the wild-type sequence and its F20W mutant in a dipalmitoylphosphatidylcholine membrane using atomistic replica exchange molecular dynamic simulations. We show that both the U-shape organization and the assembly of β-hairpins are maintained in the membrane and are resistant to the mutation F20W. In contrast the models are destabilized by the F19P mutation. Overall, our results indicate that these two assemblies represent minimal seeds or nuclei for the formation of either amyloid fibrils, a variety of β-barrel pores, or various aggregates for many Aβ sequences in a membrane-mimicking environment.
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
| | - Phuong H Nguyen
- Laboratoire de Biochimie Théorique , UPR 9080, CNRS, Université de Paris , 13 rue Pierre et Marie Curie , 75005 , Paris , France.,Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University , 75005 Paris , France
| | - Philippe Derreumaux
- Laboratory of Theoretical Chemistry , Ton Duc Thang University , Ho Chi Minh City , Vietnam.,Faculty of Pharmacy , Ton Duc Thang University , Ho Chi Minh City , Vietnam
| |
Collapse
|
21
|
Nguyen PH, Sterpone F, Derreumaux P. Aggregation of disease-related peptides. PROGRESS IN MOLECULAR BIOLOGY AND TRANSLATIONAL SCIENCE 2020; 170:435-460. [PMID: 32145950 DOI: 10.1016/bs.pmbts.2019.12.002] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
Abstract
Protein misfolding and aggregation of amyloid proteins is the fundamental cause of more than 20 diseases. Molecular mechanisms of the self-assembly and the formation of the toxic aggregates are still elusive. Computer simulations have been intensively used to study the aggregation of amyloid peptides of various amino acid lengths related to neurodegenerative diseases. We review atomistic and coarse-grained simulations of short amyloid peptides aimed at determining their transient oligomeric structures and the early and late aggregation steps.
Collapse
Affiliation(s)
- Phuong H Nguyen
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Fabio Sterpone
- CNRS, Université de Paris, UPR 9080, Laboratoire de Biochimie Théorique, Paris, France; Institut de Biologie Physico-Chimique-Fondation Edmond de Rothschild, PSL Research University, Paris, France
| | - Philippe Derreumaux
- Laboratory of Theoretical Chemistry, Ton Duc Thang University, Ho Chi Minh City, Vietnam; Faculty of Pharmacy, Ton Duc Thang University, Ho Chi Minh City, Vietnam.
| |
Collapse
|
22
|
Naderi M, Lemoine JM, Govindaraj RG, Kana OZ, Feinstein WP, Brylinski M. Binding site matching in rational drug design: algorithms and applications. Brief Bioinform 2019; 20:2167-2184. [PMID: 30169563 PMCID: PMC6954434 DOI: 10.1093/bib/bby078] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2018] [Revised: 07/18/2018] [Accepted: 07/29/2018] [Indexed: 01/06/2023] Open
Abstract
Interactions between proteins and small molecules are critical for biological functions. These interactions often occur in small cavities within protein structures, known as ligand-binding pockets. Understanding the physicochemical qualities of binding pockets is essential to improve not only our basic knowledge of biological systems, but also drug development procedures. In order to quantify similarities among pockets in terms of their geometries and chemical properties, either bound ligands can be compared to one another or binding sites can be matched directly. Both perspectives routinely take advantage of computational methods including various techniques to represent and compare small molecules as well as local protein structures. In this review, we survey 12 tools widely used to match pockets. These methods are divided into five categories based on the algorithm implemented to construct binding-site alignments. In addition to the comprehensive analysis of their algorithms, test sets and the performance of each method are described. We also discuss general pharmacological applications of computational pocket matching in drug repurposing, polypharmacology and side effects. Reflecting on the importance of these techniques in drug discovery, in the end, we elaborate on the development of more accurate meta-predictors, the incorporation of protein flexibility and the integration of powerful artificial intelligence technologies such as deep learning.
Collapse
Affiliation(s)
- Misagh Naderi
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Jeffrey Mitchell Lemoine
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, USA
- Division of Computer Science and Engineering, Louisiana State University, Baton Rouge, LA 70803, USA
| | | | - Omar Zade Kana
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Wei Pan Feinstein
- High-Performance Computing, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Michal Brylinski
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, USA
- Center for Computation & Technology, Louisiana State University, Baton Rouge, LA 70803, USA
| |
Collapse
|
23
|
Ghosh S, T D, Baul U, Vemparala S. Aggregation dynamics of charged peptides in water: Effect of salt concentration. J Chem Phys 2019; 151:074901. [DOI: 10.1063/1.5100890] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Affiliation(s)
- Susmita Ghosh
- The Institute of Mathematical Sciences, C.I.T. Campus, Taramani, Chennai 600113, India
- Homi Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
| | - Devanand T
- The Institute of Mathematical Sciences, C.I.T. Campus, Taramani, Chennai 600113, India
- Homi Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
| | - Upayan Baul
- Institue of Physics, Albert-Ludwigs-University of Freiburg, Hermann-Herder-Strasse 3, 79104 Freiburg, Germany
| | - Satyavani Vemparala
- The Institute of Mathematical Sciences, C.I.T. Campus, Taramani, Chennai 600113, India
- Homi Bhabha National Institute, Training School Complex, Anushakti Nagar, Mumbai 400094, India
| |
Collapse
|
24
|
Ngo ST, Mai BK, Derreumaux P, Vu VV. Adequate prediction for inhibitor affinity of Aβ 40 protofibril using the linear interaction energy method. RSC Adv 2019; 9:12455-12461. [PMID: 35515829 PMCID: PMC9063661 DOI: 10.1039/c9ra01177c] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Accepted: 04/11/2019] [Indexed: 11/21/2022] Open
Abstract
The search for efficient inhibitors targeting Aβ oligomers and fibrils is an important issue in Alzheimer's disease treatment. As a consequence, an accurate and computationally cheap approach to estimate the binding affinity for many ligands interacting with Aβ peptides is very important. Here, the calculated binding free energies of 30 ligands interacting with 12Aβ11-40 peptides using the linear interaction energy (LIE) approach are found to be in good correlation with experimental data (R = 0.79). The binding affinities of these complexes are also calculated by using free energy perturbation (FEP) and molecular mechanic/Poisson-Boltzmann surface area (MM/PBSA) methods. The time-consuming FEP method provides results with similar correlation (R = 0.72), whereas MM/PBSA calculations show very low correlation with experimental data (R = 0.27). In all complexes, van der Waals interactions contribute much more than electrostatic interactions. The LIE model, which is much less time-consuming than both the FEP and MM/PBSA methods, opens the door to accurate and rapid affinity prediction of ligands with Aβ peptides and the design of new ligands.
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
| | - Binh Khanh Mai
- Institute for Computational Science and Technology (ICST), Quang Trung Software City Ho Chi Minh City Vietnam
| | - Philippe Derreumaux
- Laboratoire de Biochimie Théorique, UPR 9080 CNRS, IBPC, Université Paris Diderot 13 rue Pierre et Marie Curie 75005 Paris France
- Laboratory of Theoretical Chemistry, Ton Duc Thang University Ho Chi Minh City Vietnam
- Faculty of Pharmacy, Ton Duc Thang University Ho Chi Minh City Vietnam
| | - Van V Vu
- NTT Hi-Tech Institute, Nguyen Tat Thanh University Ho Chi Minh City Vietnam
| |
Collapse
|
25
|
Matsuzaka Y, Uesawa Y. Optimization of a Deep-Learning Method Based on the Classification of Images Generated by Parameterized Deep Snap a Novel Molecular-Image-Input Technique for Quantitative Structure-Activity Relationship (QSAR) Analysis. Front Bioeng Biotechnol 2019; 7:65. [PMID: 30984753 PMCID: PMC6447703 DOI: 10.3389/fbioe.2019.00065] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2018] [Accepted: 03/07/2019] [Indexed: 12/22/2022] Open
Abstract
Numerous chemical compounds are distributed around the world and may affect the homeostasis of the endocrine system by disrupting the normal functions of hormone receptors. Although the risks associated with these compounds have been evaluated by acute toxicity testing in mammalian models, the chronic toxicity of many chemicals remains due to high cost of the compounds and the testing, etc. However, computational approaches may be promising alternatives and reduce these evaluations. Recently, deep learning (DL) has been shown to be promising prediction models with high accuracy for recognition of images, speech, signals, and videos since it greatly benefits from large datasets. Recently, a novel DL-based technique called DeepSnap was developed to conduct QSAR analysis using three-dimensional images of chemical structures. It can be used to predict the potential toxicity of many different chemicals to various receptors without extraction of descriptors. DeepSnap has been shown to have a very high capacity in tests using Tox21 quantitative qHTP datasets. Numerous parameters must be adjusted to use the DeepSnap method but they have not been optimized. In this study, the effects of these parameters on the performance of the DL prediction model were evaluated in terms of the loss in validation as an indicator for evaluating the performance of the DL using the toxicity information in the Tox21 qHTP database. The relations of the parameters of DeepSnap such as (1) number of molecules per SDF split into (2) zoom factor percentage, (3) atom size for van der waals percentage, (4) bond radius, (5) minimum bond distance, and (6) bond tolerance, with the validation loss following quadratic function curves, which suggests that optimal thresholds exist to attain the best performance with these prediction models. Using the parameter values set with the best performance, the prediction model of chemical compounds for CAR agonist was built using 64 images, at 105° angle, with AUC of 0.791. Thus, based on these parameters, the proposed DeepSnap-DL approach will be highly reliable and beneficial to establish models to assess the risk associated with various chemicals.
Collapse
Affiliation(s)
| | - Yoshihiro Uesawa
- Department of Medical Molecular Informatics, Meiji Pharmaceutical University, Tokyo, Japan
| |
Collapse
|
26
|
Borah K, Sharma S, Silla Y. Structural bioinformatics-based identification of putative plant based lead compounds for Alzheimer Disease Therapy. Comput Biol Chem 2019; 78:359-366. [DOI: 10.1016/j.compbiolchem.2018.12.012] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2018] [Accepted: 12/25/2018] [Indexed: 12/19/2022]
|
27
|
Hassan SA. Self-adaptive multiscaling algorithm for efficient simulations of many-protein systems in crowded conditions. Phys Chem Chem Phys 2018; 20:28544-28557. [PMID: 30421760 PMCID: PMC6752035 DOI: 10.1039/c8cp05517c] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A method is described for the efficient simulation of multiprotein systems in crowded environments. It is based on an adaptive, reversible structural coarsening algorithm that preserves relevant physical features of the proteins across scales. Water is treated implicitly whereas all the other components of the aqueous solution, such as ions, cosolutes, or osmolytes, are treated in atomic detail. The focus is on the analytical adaptation of the solvent model to different levels of molecular resolutions, which allows continuous, on-the-fly transitions between scales. This permits the analytical calculation of forces during dynamics and preserves detailed balance in Monte Carlo simulations. A major computational speedup can be achieved in systems containing hundreds of proteins without cutting off the long-range interactions. The method can be combined with a self-adaptive configurational-bias sampling technique described previously, designed to detect strong, weak, or ultra-weak protein associations and shown to improve sampling efficiency and convergence. The implementation aims to simulate early stages of multimeric complexation, aggregation, or self-assembly. The method can be adopted as the basis for a more general algorithm to identify vertices, edges, and hubs in protein interaction networks or to predict critical steps in signal transduction pathways.
Collapse
Affiliation(s)
- Sergio A Hassan
- Center for Molecular Modeling, OIR/CIT, National Institutes of Health, U.S. DHHS, USA.
| |
Collapse
|
28
|
Zheng Z, Pei J, Bansal N, Liu H, Song LF, Merz KM. Generation of Pairwise Potentials Using Multidimensional Data Mining. J Chem Theory Comput 2018; 14:5045-5067. [PMID: 30183299 DOI: 10.1021/acs.jctc.8b00516] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The rapid development of molecular structural databases provides the chemistry community access to an enormous array of experimental data that can be used to build and validate computational models. Using radial distribution functions collected from experimentally available X-ray and NMR structures, a number of so-called statistical potentials have been developed over the years using the structural data mining strategy. These potentials have been developed within the context of the two-particle Kirkwood equation by extending its original use for isotropic monatomic systems to anisotropic biomolecular systems. However, the accuracy and the unclear physical meaning of statistical potentials have long formed the central arguments against such methods. In this work, we present a new approach to generate molecular energy functions using structural data mining. Instead of employing the Kirkwood equation and introducing the "reference state" approximation, we model the multidimensional probability distributions of the molecular system using graphical models and generate the target pairwise Boltzmann probabilities using the Bayesian field theory. Different from the current statistical potentials that mimic the "knowledge-based" PMF based on the 2-particle Kirkwood equation, the graphical-model-based structure-derived potential developed in this study focuses on the generation of lower-dimensional Boltzmann distributions of atoms through reduction of dimensionality. We have named this new scoring function GARF, and in this work we focus on the mathematical derivation of our novel approach followed by validation studies on its ability to predict protein-ligand interactions.
Collapse
Affiliation(s)
- Zheng Zheng
- Department of Chemistry , Michigan State University , 578 South Shaw Lane , East Lansing , Michigan 48824 , United States
| | - Jun Pei
- Department of Chemistry , Michigan State University , 578 South Shaw Lane , East Lansing , Michigan 48824 , United States
| | - Nupur Bansal
- Department of Chemistry , Michigan State University , 578 South Shaw Lane , East Lansing , Michigan 48824 , United States
| | - Hao Liu
- Department of Chemistry , Michigan State University , 578 South Shaw Lane , East Lansing , Michigan 48824 , United States
| | - Lin Frank Song
- Department of Chemistry , Michigan State University , 578 South Shaw Lane , East Lansing , Michigan 48824 , United States
| | - Kenneth M Merz
- Department of Chemistry , Michigan State University , 578 South Shaw Lane , East Lansing , Michigan 48824 , United States
| |
Collapse
|
29
|
Guedes IA, Pereira FSS, Dardenne LE. Empirical Scoring Functions for Structure-Based Virtual Screening: Applications, Critical Aspects, and Challenges. Front Pharmacol 2018; 9:1089. [PMID: 30319422 PMCID: PMC6165880 DOI: 10.3389/fphar.2018.01089] [Citation(s) in RCA: 144] [Impact Index Per Article: 24.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2018] [Accepted: 09/07/2018] [Indexed: 12/19/2022] Open
Abstract
Structure-based virtual screening (VS) is a widely used approach that employs the knowledge of the three-dimensional structure of the target of interest in the design of new lead compounds from large-scale molecular docking experiments. Through the prediction of the binding mode and affinity of a small molecule within the binding site of the target of interest, it is possible to understand important properties related to the binding process. Empirical scoring functions are widely used for pose and affinity prediction. Although pose prediction is performed with satisfactory accuracy, the correct prediction of binding affinity is still a challenging task and crucial for the success of structure-based VS experiments. There are several efforts in distinct fronts to develop even more sophisticated and accurate models for filtering and ranking large libraries of compounds. This paper will cover some recent successful applications and methodological advances, including strategies to explore the ligand entropy and solvent effects, training with sophisticated machine-learning techniques, and the use of quantum mechanics. Particular emphasis will be given to the discussion of critical aspects and further directions for the development of more accurate empirical scoring functions.
Collapse
Affiliation(s)
- Isabella A Guedes
- Grupo de Modelagem Molecular em Sistemas Biológicos, Laboratório Nacional de Computação Científica, Petrópolis, Brazil
| | - Felipe S S Pereira
- Grupo de Modelagem Molecular em Sistemas Biológicos, Laboratório Nacional de Computação Científica, Petrópolis, Brazil
| | - Laurent E Dardenne
- Grupo de Modelagem Molecular em Sistemas Biológicos, Laboratório Nacional de Computação Científica, Petrópolis, Brazil
| |
Collapse
|
30
|
Phosphorylation promotes binding affinity of Rap-Raf complex by allosteric modulation of switch loop dynamics. Sci Rep 2018; 8:12976. [PMID: 30154518 PMCID: PMC6113251 DOI: 10.1038/s41598-018-31234-7] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2018] [Accepted: 08/13/2018] [Indexed: 12/14/2022] Open
Abstract
The effects of phosphorylation of a serine residue on the structural and dynamic properties of Ras-like protein, Rap, and its interactions with effector protein Ras binding domain (RBD) of Raf kinase, in the presence of GTP, are investigated via molecular dynamics simulations. The simulations show that phosphorylation significantly effects the dynamics of functional loops of Rap which participate in the stability of the complex with effector proteins. The effects of phosphorylation on Rap are significant and detailed conformational analysis suggest that the Rap protein, when phosphorylated and with GTP ligand, samples different conformational space as compared to non-phosphorylated protein. In addition, phosphorylation of SER11 opens up a new cavity in the Rap protein which can be further explored for possible drug interactions. Residue network analysis shows that the phosphorylation of Rap results in a community spanning both Rap and RBD and strongly suggests transmission of allosteric effects of local alterations in Rap to distal regions of RBD, potentially affecting the downstream signalling. Binding free energy calculations suggest that phosphorylation of SER11 residue increases the binding between Rap and Raf corroborating the network analysis results. The increased binding of the Rap-Raf complex can have cascading effects along the signalling pathways where availability of Raf can influence the oncogenic effects of Ras proteins. These simulations underscore the importance of post translational modifications like phosphorylation on the functional dynamics in proteins and can be an alternative to drug-targeting, especially in notoriously undruggable oncoproteins belonging to Ras-like GTPase family.
Collapse
|
31
|
Anti-inflammatory activity of dimethyl octenol and oleanene tetrol isolated from Trianthema decandra L. Mol Biol Rep 2018; 45:547-561. [PMID: 29797173 DOI: 10.1007/s11033-018-4193-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2018] [Accepted: 05/14/2018] [Indexed: 10/16/2022]
Abstract
Dimethyl octenol from chloroform extract and oleanene tetrol from water extract of Trianthema decandra (TD) were isolated and characterized by using HPLC, UV, FT-IR, NMR, LC-MS and CHNS, their structure were elucidated from their respective spectral data. The anti-inflammatory activity of chloroform extract, water extract, dimethyl octenol and oleanene tetrol of T. decandra were studied and underlying cellular and molecular mechanisms of action were investigated in vitro and in vivo using macrophage-like cell line (RAW264.7 cells) and type II collagen induced arthritis mice models. Nitric oxide production was inhibited and TNF-α secretion was supressed in stimulated RAW cells treated with the chloroform extract and dimethyl octenol of T. decandra. Further, the chloroform and water extract, dimethyl octenol and oleanene tetrol inhibited protein denaturation and stabilized HRBC membranes in vitro. Reduction in inflammation as a measure of paw diameter was recorded in all the treated animals when compared to control animals. Catalase, peroxidase and glutathione peroxidase levels significantly increased in the joint tissue of treated groups. The possible mechanism of action of these compounds was studied using in silico molecular docking methods with phospholipase A2 (PLA2), cycloxygenase-1 (COX-1) and cycloxygenase-2 (COX-2) as targets. Among the three target proteins, the inhibition of the inflammatory protein PLA2 and COX-2 towards dimethyl octenol and oleanene tetrol respectively. Our results contribute towards confirmation of the traditional use of TD and its compounds for the therapy of rheumatoid arthritis and other inflammatory joint disorders.
Collapse
|
32
|
Al-Sha'er MA, Taha MO. Ligand-based modeling of Akt3 lead to potent dual Akt1/Akt3 inhibitor. J Mol Graph Model 2018; 83:153-166. [PMID: 29456101 DOI: 10.1016/j.jmgm.2018.02.001] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2017] [Revised: 01/01/2018] [Accepted: 02/02/2018] [Indexed: 11/26/2022]
Abstract
Akt1 and Akt3 are important serine/threonine-specific protein kinases involved in G2 phase required by cancer cells to maintain cell cycle and to prevent cell death. Accordingly, inhibitors of these kinases should have potent anti-cancer properties. This prompted us to use pharmacophore/QSAR modeling to identify optimal binding models and physicochemical descriptors that explain bioactivity variation within a set of 74 diverse Akt3 inhibitors. Two successful orthogonal pharmacophores were identified and further validated using receiver operating characteristic (ROC) curve analyses. The pharmacophoric models and associated QSAR equation were applied to screen the national cancer institute (NCI) list of compounds for new Akt3 inhibitors. Six hits showed significant experimental anti-Akt3 IC50 values, out of which one compound exhibited dual low micromolar anti-Akt1 and anti-Akt3 inhibitory profiles.
Collapse
Affiliation(s)
| | - Mutasem O Taha
- Department of Pharmaceutical Sciences, Faculty of Pharmacy, The University of Jordan, Amman, Jordan.
| |
Collapse
|
33
|
Muthusamy K, Nagamani S. Vitamin D receptor (VDR) non-synonymous single nucleotide polymorphisms (nsSNPs) affect the calcitriol drug response - A theoretical insight. J Mol Graph Model 2018; 81:14-24. [PMID: 29476931 DOI: 10.1016/j.jmgm.2018.02.004] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2017] [Revised: 12/30/2017] [Accepted: 02/05/2018] [Indexed: 11/19/2022]
Abstract
Pharmacogenetics and pharmacogenomics have become presumptive with advancements in next-generation sequencing technology. In complex diseases, distinguishing the feasibility of pathogenic and neutral disease-causing variants is a time consuming and expensive process. Recent drug research and development processes mainly rely on the relationship between the genotype and phenotype through Single nucleotide polymorphisms (SNPs). The SNPs play an indispensable role in elucidating the individual's vulnerability to disease and drug response. The understanding of the interplay between these leads to the establishment of personalized medicine. In order to address this issue, we developed a computational pipeline of vitamin D receptor (VDR) for SNP centered study by application of elegant molecular docking and molecular dynamics simulation approaches. In a few SNPs the volume of the binding cavities has increased in mutant structures when compared to the wild type, indicating a weakening in interaction (699.1 Å3 in wild type Vs. 738.8 in Leu230Val, 820.7 Å3 in Arg247Leu). This also differently reflected in the H-bond interactions and binding free energies -169.93 kcal/mol (wild type) Vs -156.43 kcal/mol (R154W), -105.49 kcal/mol (R274L) in Leu230Val and Arg247Leu respectively. Although we could not find noteworthy changes in the binding free energies and binding pocket in the remaining mutations, the H-bond interactions made these SNPs deleterious. Thus, we further analyzed the H-bond interactions and distances using molecular dynamics (MD) simulation studies.
Collapse
Affiliation(s)
| | - Selvaraman Nagamani
- Department of Bioinformatics, Alagappa University, Karaikudi, 630 004, India
| |
Collapse
|
34
|
Mutations at protein-protein interfaces: Small changes over big surfaces have large impacts on human health. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2017; 128:3-13. [DOI: 10.1016/j.pbiomolbio.2016.10.002] [Citation(s) in RCA: 107] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/30/2016] [Revised: 10/15/2016] [Accepted: 10/19/2016] [Indexed: 12/22/2022]
|
35
|
Nguyen DD, Xiao T, Wang M, Wei GW. Rigidity Strengthening: A Mechanism for Protein–Ligand Binding. J Chem Inf Model 2017; 57:1715-1721. [DOI: 10.1021/acs.jcim.7b00226] [Citation(s) in RCA: 56] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Affiliation(s)
- Duc D. Nguyen
- Department of Mathematics, ‡Department of Biochemistry and Molecular Biology, and §Department of Electrical
and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, United States
| | - Tian Xiao
- Department of Mathematics, ‡Department of Biochemistry and Molecular Biology, and §Department of Electrical
and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, United States
| | - Menglun Wang
- Department of Mathematics, ‡Department of Biochemistry and Molecular Biology, and §Department of Electrical
and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, United States
| | - Guo-Wei Wei
- Department of Mathematics, ‡Department of Biochemistry and Molecular Biology, and §Department of Electrical
and Computer Engineering, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
36
|
Hamzeh-Mivehroud M, Sokouti B, Dastmalchi S. Molecular Docking at a Glance. Oncology 2017. [DOI: 10.4018/978-1-5225-0549-5.ch030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
The current chapter introduces different aspects of molecular docking technique in order to give an overview to the readers about the topics which will be dealt with throughout this volume. Like many other fields of science, molecular docking studies has experienced a lagging period of slow and steady increase in terms of acquiring attention of scientific community as well as its frequency of application, followed by a pronounced era of exponential expansion in theory, methodology, areas of application and performance due to developments in related technologies such as computational resources and theoretical as well as experimental biophysical methods. In the following sections the evolution of molecular docking will be reviewed and its different components including methods, search algorithms, scoring functions, validation of the methods, and area of applications plus few case studies will be touched briefly.
Collapse
Affiliation(s)
| | | | - Siavoush Dastmalchi
- Biotechnology Research Center, Tabriz University of Medical Sciences, Iran & School of Pharmacy, Tabriz University of Medical Sciences, Iran
| |
Collapse
|
37
|
Riley TP, Ayres CM, Hellman LM, Singh NK, Cosiano M, Cimons JM, Anderson MJ, Piepenbrink KH, Pierce BG, Weng Z, Baker BM. A generalized framework for computational design and mutational scanning of T-cell receptor binding interfaces. Protein Eng Des Sel 2016; 29:595-606. [PMID: 27624308 PMCID: PMC5181382 DOI: 10.1093/protein/gzw050] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2016] [Revised: 08/19/2016] [Accepted: 08/23/2016] [Indexed: 11/13/2022] Open
Abstract
T-cell receptors (TCRs) have emerged as a new class of therapeutics, most prominently for cancer where they are the key components of new cellular therapies as well as soluble biologics. Many studies have generated high affinity TCRs in order to enhance sensitivity. Recent outcomes, however, have suggested that fine manipulation of TCR binding, with an emphasis on specificity may be more valuable than large affinity increments. Structure-guided design is ideally suited for this role, and here we studied the generality of structure-guided design as applied to TCRs. We found that a previous approach, which successfully optimized the binding of a therapeutic TCR, had poor accuracy when applied to a broader set of TCR interfaces. We thus sought to develop a more general purpose TCR design framework. After assembling a large dataset of experimental data spanning multiple interfaces, we trained a new scoring function that accounted for unique features of each interface. Together with other improvements, such as explicit inclusion of molecular flexibility, this permitted the design new affinity-enhancing mutations in multiple TCRs, including those not used in training. Our approach also captured the impacts of mutations and substitutions in the peptide/MHC ligand, and recapitulated recent findings regarding TCR specificity, indicating utility in more general mutational scanning of TCR-pMHC interfaces.
Collapse
Affiliation(s)
- Timothy P Riley
- Department of Chemistry & Biochemistry and the Harper Cancer Research Institute, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, IN 46556, USA
| | - Cory M Ayres
- Department of Chemistry & Biochemistry and the Harper Cancer Research Institute, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, IN 46556, USA
| | - Lance M Hellman
- Department of Chemistry & Biochemistry and the Harper Cancer Research Institute, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, IN 46556, USA
| | - Nishant K Singh
- Department of Chemistry & Biochemistry and the Harper Cancer Research Institute, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, IN 46556, USA
| | - Michael Cosiano
- Department of Chemistry & Biochemistry and the Harper Cancer Research Institute, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, IN 46556, USA
| | - Jennifer M Cimons
- Department of Chemistry & Biochemistry and the Harper Cancer Research Institute, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, IN 46556, USA
| | - Michael J Anderson
- Department of Chemistry & Biochemistry and the Harper Cancer Research Institute, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, IN 46556, USA
| | - Kurt H Piepenbrink
- Department of Chemistry & Biochemistry and the Harper Cancer Research Institute, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, IN 46556, USA
| | - Brian G Pierce
- Institute for Bioscience and Biotechnology Research, University of Maryland, 9600 Gudelsky Drive, Rockville, MD 20850, USA
| | - Zhiping Weng
- Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, MA 01605, USA
| | - Brian M Baker
- Department of Chemistry & Biochemistry and the Harper Cancer Research Institute, University of Notre Dame, 251 Nieuwland Science Hall, Notre Dame, IN 46556, USA
| |
Collapse
|
38
|
Hubbard TA, Brown AJ, Bell IAW, Cockroft SL. The Limit of Intramolecular H-Bonding. J Am Chem Soc 2016; 138:15114-15117. [DOI: 10.1021/jacs.6b09130] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Affiliation(s)
| | - Alisdair J. Brown
- Afton Chemical Limited, London Road, Bracknell, Berkshire RG12 2UW, U.K
| | - Ian A. W. Bell
- Afton Chemical Limited, London Road, Bracknell, Berkshire RG12 2UW, U.K
| | | |
Collapse
|
39
|
In silico structural characterization of protein targets for drug development against Trypanosoma cruzi. J Mol Model 2016; 22:244. [DOI: 10.1007/s00894-016-3115-9] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2016] [Accepted: 09/02/2016] [Indexed: 10/21/2022]
|
40
|
Kmiecik S, Gront D, Kolinski M, Wieteska L, Dawid AE, Kolinski A. Coarse-Grained Protein Models and Their Applications. Chem Rev 2016; 116:7898-936. [DOI: 10.1021/acs.chemrev.6b00163] [Citation(s) in RCA: 555] [Impact Index Per Article: 69.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Sebastian Kmiecik
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Dominik Gront
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| | - Michal Kolinski
- Bioinformatics
Laboratory, Mossakowski Medical Research Center of the Polish Academy of Sciences, Pawinskiego 5, 02-106 Warsaw, Poland
| | - Lukasz Wieteska
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
- Department
of Medical Biochemistry, Medical University of Lodz, Mazowiecka 6/8, 92-215 Lodz, Poland
| | | | - Andrzej Kolinski
- Faculty
of Chemistry, University of Warsaw, Pasteura 1, 02-093 Warsaw, Poland
| |
Collapse
|
41
|
Sessa L, Di Biasi L, Concilio S, Cattaneo G, De Santis A, Iannelli P, Piotto S. A New Flexible Protocol for Docking Studies. COMMUNICATIONS IN COMPUTER AND INFORMATION SCIENCE 2016. [DOI: 10.1007/978-3-319-32695-5_11] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
|
42
|
Abstract
Drug discovery utilizes chemical biology and computational drug design approaches for the efficient identification and optimization of lead compounds. Chemical biology is mostly involved in the elucidation of the biological function of a target and the mechanism of action of a chemical modulator. On the other hand, computer-aided drug design makes use of the structural knowledge of either the target (structure-based) or known ligands with bioactivity (ligand-based) to facilitate the determination of promising candidate drugs. Various virtual screening techniques are now being used by both pharmaceutical companies and academic research groups to reduce the cost and time required for the discovery of a potent drug. Despite the rapid advances in these methods, continuous improvements are critical for future drug discovery tools. Advantages presented by structure-based and ligand-based drug design suggest that their complementary use, as well as their integration with experimental routines, has a powerful impact on rational drug design. In this article, we give an overview of the current computational drug design and their application in integrated rational drug development to aid in the progress of drug discovery research.
Collapse
Affiliation(s)
- Stephani Joy Y Macalino
- National Leading Research Laboratory of Molecular Modeling and Drug Design, College of Pharmacy and Graduate School of Pharmaceutical Sciences, and Global Top 5 Research Program, Ewha Womans University, Seoul, 120-750, Korea
| | - Vijayakumar Gosu
- National Leading Research Laboratory of Molecular Modeling and Drug Design, College of Pharmacy and Graduate School of Pharmaceutical Sciences, and Global Top 5 Research Program, Ewha Womans University, Seoul, 120-750, Korea
| | - Sunhye Hong
- National Leading Research Laboratory of Molecular Modeling and Drug Design, College of Pharmacy and Graduate School of Pharmaceutical Sciences, and Global Top 5 Research Program, Ewha Womans University, Seoul, 120-750, Korea
| | - Sun Choi
- National Leading Research Laboratory of Molecular Modeling and Drug Design, College of Pharmacy and Graduate School of Pharmaceutical Sciences, and Global Top 5 Research Program, Ewha Womans University, Seoul, 120-750, Korea.
| |
Collapse
|
43
|
Cardone A, Bornstein A, Pant HC, Brady M, Sriram R, Hassan SA. Detection and characterization of nonspecific, sparsely populated binding modes in the early stages of complexation. J Comput Chem 2015; 36:983-95. [PMID: 25782918 DOI: 10.1002/jcc.23883] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2014] [Revised: 02/02/2015] [Accepted: 02/08/2015] [Indexed: 12/11/2022]
Abstract
A method is proposed to study protein-ligand binding in a system governed by specific and nonspecific interactions. Strong associations lead to narrow distributions in the proteins configuration space; weak and ultraweak associations lead instead to broader distributions, a manifestation of nonspecific, sparsely populated binding modes with multiple interfaces. The method is based on the notion that a discrete set of preferential first-encounter modes are metastable states from which stable (prerelaxation) complexes at equilibrium evolve. The method can be used to explore alternative pathways of complexation with statistical significance and can be integrated into a general algorithm to study protein interaction networks. The method is applied to a peptide-protein complex. The peptide adopts several low-population conformers and binds in a variety of modes with a broad range of affinities. The system is thus well suited to analyze general features of binding, including conformational selection, multiplicity of binding modes, and nonspecific interactions, and to illustrate how the method can be applied to study these problems systematically. The equilibrium distributions can be used to generate biasing functions for simulations of multiprotein systems from which bulk thermodynamic quantities can be calculated.
Collapse
Affiliation(s)
- Antonio Cardone
- Software and System Division, National Institute of Standards and Technology, Gaithersburg, Maryland, 20899; Institute for Advanced Computer Studies, University of Maryland, College Park, Maryland, 20742
| | | | | | | | | | | |
Collapse
|
44
|
Haxton TK, Mannige RV, Zuckermann RN, Whitelam S. Modeling Sequence-Specific Polymers Using Anisotropic Coarse-Grained Sites Allows Quantitative Comparison with Experiment. J Chem Theory Comput 2014; 11:303-15. [DOI: 10.1021/ct5010559] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Thomas K. Haxton
- Molecular
Foundry, Lawrence Berkeley
National Laboratory, Berkeley, California 94720, United States
| | - Ranjan V. Mannige
- Molecular
Foundry, Lawrence Berkeley
National Laboratory, Berkeley, California 94720, United States
| | - Ronald N. Zuckermann
- Molecular
Foundry, Lawrence Berkeley
National Laboratory, Berkeley, California 94720, United States
| | - Stephen Whitelam
- Molecular
Foundry, Lawrence Berkeley
National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
45
|
Smith LJ, Gunsteren WFV, Allison JR. Multiple binding modes for palmitate to barley lipid transfer protein facilitated by the presence of proline 12. Protein Sci 2014; 22:56-64. [PMID: 23139016 DOI: 10.1002/pro.2184] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/04/2012] [Revised: 10/17/2012] [Accepted: 10/19/2012] [Indexed: 11/09/2022]
Abstract
Molecular dynamics simulations have been used to characterise the binding of the fatty acid ligand palmitate in the barley lipid transfer protein 1 (LTP) internal cavity. Two different palmitate binding modes (1 and 2), with similar protein-ligand interaction energies, have been identified using a variety of simulation strategies. These strategies include applying experimental protein-ligand atom-atom distance restraints during the simulation, or protonating the palmitate ligand, or using the vacuum GROMOS 54B7 force-field parameter set for the ligand during the initial stages of the simulations. In both the binding modes identified the palmitate carboxylate head group hydrogen bonds with main chain amide groups in helix A, residues 4 to 19, of the protein. In binding mode 1 the hydrogen bonds are to Lys 11, Cys 13, and Leu 14 and in binding mode 2 to Thr 15, Tyr 16, Val 17, Ser 24 and also to the OH of Thr 15. In both cases palmitate binding exploits irregularity of the intrahelical hydrogen-bonding pattern in helix A of barley LTP due to the presence of Pro 12. Simulations of two variants of barley LTP, namely the single mutant Pro12Val and the double mutant Pro12Val Pro70Val, show that Pro 12 is required for persistent palmitate binding in the LTP cavity. Overall, the work identifies key MD simulation approaches for characterizing the details of protein-ligand interactions in complexes where NMR data provide insufficient restraints.
Collapse
Affiliation(s)
- Lorna J Smith
- Department of Chemistry, University of Oxford, Inorganic Chemistry Laboratory, Oxford OX1 3QR, United Kingdom.
| | | | | |
Collapse
|
46
|
Chemical specificity and conformational flexibility in proteinase-inhibitor interaction: scaffolds for promiscuous binding. PROGRESS IN BIOPHYSICS AND MOLECULAR BIOLOGY 2014; 116:151-7. [PMID: 25151636 DOI: 10.1016/j.pbiomolbio.2014.08.003] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/16/2014] [Revised: 07/26/2014] [Accepted: 08/02/2014] [Indexed: 11/24/2022]
Abstract
One of the most important roles of proteins in cellular milieu is recognition of other biomolecules including other proteins. Protein-protein complexes are involved in many essential cellular processes. Interfaces of protein-protein complexes are traditionally known to be conserved in evolution and less flexible than other solvent interacting tertiary structural surface. But many examples are emerging where these features do not hold good. An understanding of inter-play between flexibility and sequence conservation is emerging, providing a fresh dimension to the paradigm of sequence-structure-function relationship. The functional manifestation of the inter-relation between sequence conservation and flexibility of interface is exemplified in this review using proteinase-inhibitor protein complexes.
Collapse
|
47
|
Molecular Dynamics Simulations of Bromodomains Reveal Binding-Site Flexibility and Multiple Binding Modes of the Natural Ligand Acetyl-Lysine. Isr J Chem 2014. [DOI: 10.1002/ijch.201400009] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
|
48
|
Brut M, Estève A, Landa G, Djafari Rouhani M. Toward in silico biomolecular manipulation through static modes: atomic scale characterization of HIV-1 protease flexibility. J Phys Chem B 2014; 118:2821-30. [PMID: 24568689 DOI: 10.1021/jp4113156] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Probing biomolecular flexibility with atomic-scale resolution is a challenging task in current computational biology for fundamental understanding and prediction of biomolecular interactions and associated functions. This paper makes use of the static mode method to study HIV-1 protease considered as a model system to investigate the full biomolecular flexibility at the atomic scale, the screening of active site biomechanical properties, the blind prediction of allosteric sites, and the design of multisite strategies to target deformations of interest. Relying on this single calculation run of static modes, we demonstrate that in silico predictive design of an infinite set of complex excitation fields is reachable, thanks to the storage of the static modes in a data bank that can be used to mimic single or multiatom contact and efficiently anticipate conformational changes arising from external stimuli. All along this article, we compare our results to data previously published and propose a guideline for efficient, predictive, and custom in silico experiments.
Collapse
Affiliation(s)
- Marie Brut
- CNRS , LAAS, 7 avenue du colonel Roche, F-31400 Toulouse, France
| | | | | | | |
Collapse
|
49
|
The conformational musings of a medicinal chemist. Drug Discov Today 2014; 19:320-5. [DOI: 10.1016/j.drudis.2013.10.016] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2013] [Revised: 09/04/2013] [Accepted: 10/22/2013] [Indexed: 11/22/2022]
|
50
|
Structure-based design of small-molecule protein–protein interaction modulators: the story so far. Future Med Chem 2014; 6:343-57. [DOI: 10.4155/fmc.13.204] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
As the pivotal role of protein–protein interactions in cell growth, transcriptional activity, intracellular trafficking, signal transduction and pathological conditions has been assessed, experimental and in silico strategies have been developed to design protein–protein interaction modulators. State-of-the-art structure-based design methods, mainly pharmacophore modeling and docking, which have succeeded in the identification of enzyme inhibitors, receptor agonists and antagonists, and new tools specifically conceived to target protein–protein interfaces (e.g., hot-spot and druggable pocket prediction methods) have been applied in the search for small-molecule protein–protein interaction modulators. Many successful applications of structure-based design approaches that, despite the challenge of targeting protein–protein interfaces with small molecules, have led to the identification of micromolar and submicromolar hits are reviewed here.
Collapse
|