1
|
Apostolopoulos V, Bojarska J, Chai TT, Elnagdy S, Kaczmarek K, Matsoukas J, New R, Parang K, Lopez OP, Parhiz H, Perera CO, Pickholz M, Remko M, Saviano M, Skwarczynski M, Tang Y, Wolf WM, Yoshiya T, Zabrocki J, Zielenkiewicz P, AlKhazindar M, Barriga V, Kelaidonis K, Sarasia EM, Toth I. A Global Review on Short Peptides: Frontiers and Perspectives. Molecules 2021; 26:E430. [PMID: 33467522 PMCID: PMC7830668 DOI: 10.3390/molecules26020430] [Citation(s) in RCA: 168] [Impact Index Per Article: 56.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2020] [Revised: 12/23/2020] [Accepted: 01/09/2021] [Indexed: 12/13/2022] Open
Abstract
Peptides are fragments of proteins that carry out biological functions. They act as signaling entities via all domains of life and interfere with protein-protein interactions, which are indispensable in bio-processes. Short peptides include fundamental molecular information for a prelude to the symphony of life. They have aroused considerable interest due to their unique features and great promise in innovative bio-therapies. This work focusing on the current state-of-the-art short peptide-based therapeutical developments is the first global review written by researchers from all continents, as a celebration of 100 years of peptide therapeutics since the commencement of insulin therapy in the 1920s. Peptide "drugs" initially played only the role of hormone analogs to balance disorders. Nowadays, they achieve numerous biomedical tasks, can cross membranes, or reach intracellular targets. The role of peptides in bio-processes can hardly be mimicked by other chemical substances. The article is divided into independent sections, which are related to either the progress in short peptide-based theranostics or the problems posing challenge to bio-medicine. In particular, the SWOT analysis of short peptides, their relevance in therapies of diverse diseases, improvements in (bio)synthesis platforms, advanced nano-supramolecular technologies, aptamers, altered peptide ligands and in silico methodologies to overcome peptide limitations, modern smart bio-functional materials, vaccines, and drug/gene-targeted delivery systems are discussed.
Collapse
Affiliation(s)
- Vasso Apostolopoulos
- Institute for Health and Sport, Victoria University, Melbourne, VIC 3030, Australia; (V.A.); (J.M.); (V.B.)
| | - Joanna Bojarska
- Institute of General and Ecological Chemistry, Faculty of Chemistry, Lodz University of Technology, Żeromskiego 116, 90-924 Lodz, Poland
| | - Tsun-Thai Chai
- Department of Chemical Science, Faculty of Science, Universiti Tunku Abdul Rahman, Kampar 31900, Malaysia;
| | - Sherif Elnagdy
- Botany and Microbiology Department, Faculty of Science, Cairo University, Gamaa St., Giza 12613, Egypt; (S.E.); (M.A.)
| | - Krzysztof Kaczmarek
- Institute of Organic Chemistry, Faculty of Chemistry, Lodz University of Technology, Żeromskiego 116, 90-924 Lodz, Poland; (K.K.); (J.Z.)
| | - John Matsoukas
- Institute for Health and Sport, Victoria University, Melbourne, VIC 3030, Australia; (V.A.); (J.M.); (V.B.)
- NewDrug, Patras Science Park, 26500 Patras, Greece;
- Department of Physiology and Pharmacology, Cumming School of Medicine, University of Calgary, Calgary, AB T2N 4N1, Canada
| | - Roger New
- Vaxcine (UK) Ltd., c/o London Bioscience Innovation Centre, London NW1 0NH, UK;
- Faculty of Science & Technology, Middlesex University, The Burroughs, London NW4 4BT, UK;
| | - Keykavous Parang
- Center for Targeted Drug Delivery, Department of Biomedical and Pharmaceutical Sciences, Chapman University School of Pharmacy, Harry and Diane Rinker Health Science Campus, Irvine, CA 92618, USA;
| | - Octavio Paredes Lopez
- Centro de Investigación y de Estudios Avanzados del IPN, Departamento de Biotecnología y Bioquímica, Irapuato 36824, Guanajuato, Mexico;
| | - Hamideh Parhiz
- Infectious Disease Division, Department of Medicine, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104-6073, USA;
| | - Conrad O. Perera
- School of Chemical Sciences, The University of Auckland, Private Bag 92019, Auckland 1142, New Zealand;
| | - Monica Pickholz
- Departamento de Física, Facultad de Ciencias Exactas y Naturales, Universidad de Buenos Aires, Buenos Aires 1428, Argentina;
- Instituto de Física de Buenos Aires (IFIBA, UBA-CONICET), Argentina, Buenos Aires 1428, Argentina
| | - Milan Remko
- Remedika, Luzna 9, 85104 Bratislava, Slovakia;
| | - Michele Saviano
- Institute of Crystallography (CNR), Via Amendola 122/o, 70126 Bari, Italy;
| | - Mariusz Skwarczynski
- School of Chemistry & Molecular Biosciences, The University of Queensland, St Lucia, QLD 4072, Australia; (M.S.); (I.T.)
| | - Yefeng Tang
- Key Laboratory of Bioorganic Phosphorus Chemistry & Chemical Biology (MOE), School of Pharma Ceutical Sciences, Tsinghua University, Beijing 100084, China;
| | - Wojciech M. Wolf
- Institute of General and Ecological Chemistry, Faculty of Chemistry, Lodz University of Technology, Żeromskiego 116, 90-924 Lodz, Poland
| | | | - Janusz Zabrocki
- Institute of Organic Chemistry, Faculty of Chemistry, Lodz University of Technology, Żeromskiego 116, 90-924 Lodz, Poland; (K.K.); (J.Z.)
| | - Piotr Zielenkiewicz
- Institute of Biochemistry and Biophysics, Polish Academy of Sciences, Pawinskiego 5a, 02-106 Warsaw, Poland;
- Department of Systems Biology, Institute of Experimental Plant Biology and Biotechnology, University of Warsaw, Miecznikowa 1, 02-096 Warsaw, Poland
| | - Maha AlKhazindar
- Botany and Microbiology Department, Faculty of Science, Cairo University, Gamaa St., Giza 12613, Egypt; (S.E.); (M.A.)
| | - Vanessa Barriga
- Institute for Health and Sport, Victoria University, Melbourne, VIC 3030, Australia; (V.A.); (J.M.); (V.B.)
| | | | | | - Istvan Toth
- School of Chemistry & Molecular Biosciences, The University of Queensland, St Lucia, QLD 4072, Australia; (M.S.); (I.T.)
- Institute for Molecular Bioscience, The University of Queensland, St Lucia, QLD 4072, Australia
- School of Pharmacy, The University of Queensland, Woolloongabba, QLD 4102, Australia
| |
Collapse
|
2
|
Saleh N, Ibrahim P, Saladino G, Gervasio FL, Clark T. An Efficient Metadynamics-Based Protocol To Model the Binding Affinity and the Transition State Ensemble of G-Protein-Coupled Receptor Ligands. J Chem Inf Model 2017; 57:1210-1217. [PMID: 28453271 DOI: 10.1021/acs.jcim.6b00772] [Citation(s) in RCA: 91] [Impact Index Per Article: 13.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
A generally applicable metadynamics scheme for predicting the free energy profile of ligand binding to G-protein-coupled receptors (GPCRs) is described. A common and effective collective variable (CV) has been defined using the ideally placed and highly conserved Trp6.48 as a reference point for ligand-GPCR distance measurement and the common orientation of GPCRs in the cell membrane. Using this single CV together with well-tempered multiple-walker metadynamics with a funnel-like boundary allows an efficient exploration of the entire ligand binding path from the extracellular medium to the orthosteric binding site, including vestibule and intermediate sites. The protocol can be used with X-ray structures or high-quality homology models (based on a high-quality template and after thorough refinement) for the receptor and is universally applicable to agonists, antagonists, and partial and reverse agonists. The root-mean-square error (RMSE) in predicted binding free energies for 12 diverse ligands in five receptors (a total of 23 data points) is surprisingly small (less than 1 kcal mol-1). The RMSEs for simulations that use receptor X-ray structures and homology models are very similar.
Collapse
Affiliation(s)
- Noureldin Saleh
- Computer-Chemie-Centrum and Interdisciplinary Center for Molecular Materials Friedrich-Alexander-Universität Erlangen-Nürnberg , Nägelsbachstraße 25, 91052 Erlangen, Germany
| | - Passainte Ibrahim
- Computer-Chemie-Centrum and Interdisciplinary Center for Molecular Materials Friedrich-Alexander-Universität Erlangen-Nürnberg , Nägelsbachstraße 25, 91052 Erlangen, Germany
| | - Giorgio Saladino
- Department of Chemistry, University College London , London WC1H 0AJ, United Kingdom
| | - Francesco Luigi Gervasio
- Department of Chemistry, University College London , London WC1H 0AJ, United Kingdom.,Institute of Structural and Molecular Biology, University College London , London WC1E 6BT, United Kingdom
| | - Timothy Clark
- Computer-Chemie-Centrum and Interdisciplinary Center for Molecular Materials Friedrich-Alexander-Universität Erlangen-Nürnberg , Nägelsbachstraße 25, 91052 Erlangen, Germany
| |
Collapse
|
3
|
Kuzmanic A, Sutto L, Saladino G, Nebreda AR, Gervasio FL, Orozco M. Changes in the free-energy landscape of p38α MAP kinase through its canonical activation and binding events as studied by enhanced molecular dynamics simulations. eLife 2017; 6. [PMID: 28445123 PMCID: PMC5406204 DOI: 10.7554/elife.22175] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2016] [Accepted: 04/06/2017] [Indexed: 01/03/2023] Open
Abstract
p38α is a Ser/Thr protein kinase involved in a variety of cellular processes and pathological conditions, which makes it a promising pharmacological target. Although the activity of the enzyme is highly regulated, its molecular mechanism of activation remains largely unexplained, even after decades of research. By using state-of-the-art molecular dynamics simulations, we decipher the key elements of the complex molecular mechanism refined by evolution to allow for a fine tuning of p38α kinase activity. Our study describes for the first time the molecular effects of different regulators of the enzymatic activity, and provides an integrative picture of the activation mechanism that explains the seemingly contradictory X-ray and NMR data. DOI:http://dx.doi.org/10.7554/eLife.22175.001
Collapse
Affiliation(s)
- Antonija Kuzmanic
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Ludovico Sutto
- Department of Chemistry, University College London, London, United Kingdom
| | - Giorgio Saladino
- Department of Chemistry, University College London, London, United Kingdom
| | - Angel R Nebreda
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Barcelona, Spain.,Catalan Institution for Research and Advanced Studies (ICREA), Barcelona, Spain
| | | | - Modesto Orozco
- Institute for Research in Biomedicine (IRB Barcelona), The Barcelona Institute of Science and Technology, Barcelona, Spain.,Joint BSC-CRG-IRB Program in Computational Biology, Barcelona, Spain.,Department of Biochemistry, University of Barcelona, Barcelona, Spain
| |
Collapse
|
4
|
Meloni R, Camilloni C, Tiana G. Properties of low-dimensional collective variables in the molecular dynamics of biopolymers. Phys Rev E 2016; 94:052406. [PMID: 27967023 DOI: 10.1103/physreve.94.052406] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2016] [Indexed: 11/07/2022]
Abstract
The description of the dynamics of a complex, high-dimensional system in terms of a low-dimensional set of collective variables Y can be fruitful if the low-dimensional representation satisfies a Langevin equation with drift and diffusion coefficients that depend only on Y. We present a computational scheme to evaluate whether a given collective variable provides a faithful low-dimensional representation of the dynamics of a high-dimensional system. The scheme is based on the framework of a finite-difference Langevin equation, similar to that used for molecular-dynamics simulations. This allows one to calculate the drift and diffusion coefficients in any point of the full-dimensional system. The width of the distribution of drift and diffusion coefficients in an ensemble of microscopic points at the same value of Y indicates to what extent the dynamics of Y is described by a simple Langevin equation. Using a simple protein model, we show that collective variables often used to describe biopolymers display a non-negligible width both in the drift and in the diffusion coefficients. We also show that the associated effective force is compatible with the equilibrium free energy calculated from a microscopic sampling, but it results in markedly different dynamical properties.
Collapse
Affiliation(s)
- Roberto Meloni
- Department of Physics, Università degli Studi di Milano, and INFN, via Celoria 16, 20133 Milano, Italy
| | - Carlo Camilloni
- Department of Chemistry and Institute for Advanced Study, Technische Universität München, Lichtenbergstrasse 4, 85747 Garching, Germany
| | - Guido Tiana
- Center for Complexity and Biosystems and Department of Physics, Università degli Studi di Milano, and INFN, via Celoria 16, 20133 Milano, Italy
| |
Collapse
|
5
|
Mahanti M, Bhakat S, Nilsson UJ, Söderhjelm P. Flap Dynamics in Aspartic Proteases: A Computational Perspective. Chem Biol Drug Des 2016; 88:159-77. [PMID: 26872937 DOI: 10.1111/cbdd.12745] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Recent advances in biochemistry and drug design have placed proteases as one of the critical target groups for developing novel small-molecule inhibitors. Among all proteases, aspartic proteases have gained significant attention due to their role in HIV/AIDS, malaria, Alzheimer's disease, etc. The binding cleft is covered by one or two β-hairpins (flaps) which need to be opened before a ligand can bind. After binding, the flaps close to retain the ligand in the active site. Development of computational tools has improved our understanding of flap dynamics and its role in ligand recognition. In the past decade, several computational approaches, for example molecular dynamics (MD) simulations, coarse-grained simulations, replica-exchange molecular dynamics (REMD) and metadynamics, have been used to understand flap dynamics and conformational motions associated with flap movements. This review is intended to summarize the computational progress towards understanding the flap dynamics of proteases and to be a reference for future studies in this field.
Collapse
Affiliation(s)
- Mukul Mahanti
- Centre for Analysis and Synthesis, Department of Chemistry, Lund University, PO Box 124, SE-221 00, Lund, Sweden
| | - Soumendranath Bhakat
- Division of Biophysical Chemistry, Department of Chemistry, Lund University, PO Box 124, SE-22100, Lund, Sweden
| | - Ulf J Nilsson
- Centre for Analysis and Synthesis, Department of Chemistry, Lund University, PO Box 124, SE-221 00, Lund, Sweden
| | - Pär Söderhjelm
- Division of Biophysical Chemistry, Department of Chemistry, Lund University, PO Box 124, SE-22100, Lund, Sweden
| |
Collapse
|
6
|
De Vivo M, Masetti M, Bottegoni G, Cavalli A. Role of Molecular Dynamics and Related Methods in Drug Discovery. J Med Chem 2016; 59:4035-61. [DOI: 10.1021/acs.jmedchem.5b01684] [Citation(s) in RCA: 538] [Impact Index Per Article: 67.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Affiliation(s)
- Marco De Vivo
- Laboratory
of Molecular Modeling and Drug Discovery, Istituto Italiano di Tecnologia, Via Morego 30, 16163 Genova, Italy
- IAS-5/INM-9 Computational
Biomedicine Forschungszentrum Jülich, Wilhelm-Johnen-Straße, 52428 Jülich, Germany
| | - Matteo Masetti
- Department
of Pharmacy and Biotechnology, University of Bologna, Via Belmeloro
6, I-40126 Bologna, Italy
| | - Giovanni Bottegoni
- CompuNet, Istituto
Italiano di Tecnologia, Via Morego
30, 16163 Genova, Italy
- BiKi Technologies
srl, Via XX Settembre 33/10, 16121 Genova, Italy
| | - Andrea Cavalli
- Department
of Pharmacy and Biotechnology, University of Bologna, Via Belmeloro
6, I-40126 Bologna, Italy
- CompuNet, Istituto
Italiano di Tecnologia, Via Morego
30, 16163 Genova, Italy
| |
Collapse
|
7
|
Zheng S, Pfaendtner J. Enhanced sampling of chemical and biochemical reactions with metadynamics. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.923574] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
8
|
Mata-Munguía C, Escoto-Delgadillo M, Torres-Mendoza B, Flores-Soto M, Vázquez-Torres M, Gálvez-Gastelum F, Viniegra-Osorio A, Castillero-Manzano M, Vázquez-Valls E. Natural polymorphisms and unusual mutations in HIV-1 protease with potential antiretroviral resistance: a bioinformatic analysis. BMC Bioinformatics 2014; 15:72. [PMID: 24629078 PMCID: PMC4003850 DOI: 10.1186/1471-2105-15-72] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2012] [Accepted: 03/05/2014] [Indexed: 11/22/2022] Open
Abstract
Background The correlations of genotypic and phenotypic tests with treatment, clinical history and the significance of mutations in viruses of HIV-infected patients are used to establish resistance mutations to protease inhibitors (PIs). Emerging mutations in human immunodeficiency virus type 1 (HIV-1) protease confer resistance to PIs by inducing structural changes at the ligand interaction site. The aim of this study was to establish an in silico structural relationship between natural HIV-1 polymorphisms and unusual HIV-1 mutations that confer resistance to PIs. Results Protease sequences isolated from 151 Mexican HIV-1 patients that were naïve to, or subjected to antiretroviral therapy, were examined. We identified 41 unrelated resistance mutations with a prevalence greater than 1%. Among these mutations, nine exhibited positive selection, three were natural polymorphisms (L63S/V/H) in a codon associated with drug resistance, and six were unusual mutations (L5F, D29V, L63R/G, P79L and T91V). The D29V mutation, with a prevalence of 1.32% in the studied population, was only found in patients treated with antiretroviral drugs. Using in silico modelling, we observed that D29V formed unstable protease complexes when were docked with lopinavir, saquinavir, darunavir, tipranavir, indinavir and atazanavir. Conclusions The structural correlation of natural polymorphisms and unusual mutations with drug resistance is useful for the identification of HIV-1 variants with potential resistance to PIs. The D29V mutation likely confers a selection advantage in viruses; however, in silico, presence of this mutation results in unstable enzyme/PI complexes, that possibly induce resistance to PIs.
Collapse
Affiliation(s)
| | | | | | | | | | | | | | | | - Eduardo Vázquez-Valls
- Laboratorio de Inmunodeficiencias y Retrovirus Humanos, Centro de Investigación Biomédica de Occidente, CMNO, IMSS, Guadalajara 44340, México.
| |
Collapse
|
9
|
Tiana G, Camilloni C. Ratcheted molecular-dynamics simulations identify efficiently the transition state of protein folding. J Chem Phys 2012; 137:235101. [DOI: 10.1063/1.4769085] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
10
|
Using metadynamics and path collective variables to study ligand binding and induced conformational transitions. Methods Mol Biol 2012; 819:501-13. [PMID: 22183554 DOI: 10.1007/978-1-61779-465-0_29] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Large-scale conformational transitions represent both a challenge and an opportunity for computational drug design. Exploring the conformational space of a druggable target with sufficient detail is computationally demanding. However, if it were possible to fully account for target flexibility, one could exploit this knowledge to rationally design more potent and more selective drug candidates. Here, we discuss how molecular dynamics together with free energy algorithms based on Metadynamics and Path Collective Variables can be used to study both large-scale conformational transitions and ligand binding to flexible targets. We show real-life examples of how these methods have been applied in the case of cyclin-dependent kinases, a family of flexible targets that shows promise in cancer therapy.
Collapse
|
11
|
Barducci A, Bonomi M, Parrinello M. Metadynamics. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2011. [DOI: 10.1002/wcms.31] [Citation(s) in RCA: 712] [Impact Index Per Article: 54.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Alessandro Barducci
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI Campus, Lugano, Switzerland
| | - Massimiliano Bonomi
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI Campus, Lugano, Switzerland
| | - Michele Parrinello
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, USI Campus, Lugano, Switzerland
| |
Collapse
|
12
|
Bonomi M, Barducci A, Gervasio FL, Parrinello M. Multiple routes and milestones in the folding of HIV-1 protease monomer. PLoS One 2010; 5:e13208. [PMID: 20967249 PMCID: PMC2954147 DOI: 10.1371/journal.pone.0013208] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2010] [Accepted: 09/11/2010] [Indexed: 11/25/2022] Open
Abstract
Proteins fold on a time scale incompatible with a mechanism of random search in conformational space thus indicating that somehow they are guided to the native state through a funneled energetic landscape. At the same time the heterogeneous kinetics suggests the existence of several different folding routes. Here we propose a scenario for the folding mechanism of the monomer of HIV–1 protease in which multiple pathways and milestone events coexist. A variety of computational approaches supports this picture. These include very long all-atom molecular dynamics simulations in explicit solvent, an analysis of the network of clusters found in multiple high-temperature unfolding simulations and a complete characterization of free-energy surfaces carried out using a structure-based potential at atomistic resolution and a combination of metadynamics and parallel tempering. Our results confirm that the monomer in solution is stable toward unfolding and show that at least two unfolding pathways exist. In our scenario, the formation of a hydrophobic core is a milestone in the folding process which must occur along all the routes that lead this protein towards its native state. Furthermore, the ensemble of folding pathways proposed here substantiates a rational drug design strategy based on inhibiting the folding of HIV–1 protease.
Collapse
Affiliation(s)
- Massimiliano Bonomi
- Computational Science, Department of Chemistry and Applied Biosciences, ETH Zurich, Lugano, Switzerland.
| | | | | | | |
Collapse
|
13
|
Linking well-tempered metadynamics simulations with experiments. Biophys J 2010; 98:L44-6. [PMID: 20441734 DOI: 10.1016/j.bpj.2010.01.033] [Citation(s) in RCA: 52] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2009] [Revised: 01/27/2010] [Accepted: 01/28/2010] [Indexed: 11/23/2022] Open
Abstract
Linking experiments with the atomistic resolution provided by molecular dynamics simulations can shed light on the structure and dynamics of protein-disordered states. The sampling limitations of classical molecular dynamics can be overcome using metadynamics, which is based on the introduction of a history-dependent bias on a small number of suitably chosen collective variables. Even if such bias distorts the probability distribution of the other degrees of freedom, the equilibrium Boltzmann distribution can be reconstructed using a recently developed reweighting algorithm. Quantitative comparison with experimental data is thus possible. Here we show the potential of this combined approach by characterizing the conformational ensemble explored by a 13-residue helix-forming peptide by means of a well-tempered metadynamics/parallel tempering approach and comparing the reconstructed nuclear magnetic resonance scalar couplings with experimental data.
Collapse
|
14
|
Antoniou X, Borsello T. Cell Permeable Peptides: A Promising Tool to Deliver Neuroprotective Agents in the Brain. Pharmaceuticals (Basel) 2010; 3:379-392. [PMID: 27713257 PMCID: PMC4033915 DOI: 10.3390/ph3020379] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2008] [Revised: 01/20/2010] [Accepted: 01/28/2010] [Indexed: 12/03/2022] Open
Abstract
The inability of most drugs to cross the blood-brain barrier and/or plasma membrane limits their use for biomedical applications in the brain. Cell Permeable Peptides (CPPs) overcome this problem and are effective in vivo, crossing the plasma membrane and the blood-brain barrier. CPPs deliver a wide variety of compounds intracellularly in an active form. In fact, many bioactive cargoes have neuroprotective properties, and due to their ability to block protein-protein interactions, offer exciting perspectives in the clinical setting. In this review we give an overview of the Cell Permeable Peptides strategy to deliver neuroprotectants against neurodegeneration in the CNS.
Collapse
Affiliation(s)
- Xanthi Antoniou
- Istituto di Ricerche Farmacologiche "Mario Negri", Via La Masa 19, 20157 Milano, Italy.
| | - Tiziana Borsello
- Istituto di Ricerche Farmacologiche "Mario Negri", Via La Masa 19, 20157 Milano, Italy.
| |
Collapse
|
15
|
Approaches to the design of HIV protease inhibitors with improved resistance profiles. Curr Opin HIV AIDS 2009; 3:633-41. [PMID: 19373035 DOI: 10.1097/coh.0b013e328313911d] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
PURPOSE OF REVIEW This review describes current approaches to HIV protease inhibitor design, with a focus on improving their profile against drug-resistant mutants. Potential explanations for the flat resistance profile of some potent protease inhibitors and discrepancies between the apparent fold change of potency at the enzyme level and in cell-based assays are discussed. RECENT FINDINGS Despite new ideas and a clear rationale for designing inhibitors that bind outside the enzyme active site, all current protease inhibitors with potent antiviral activity target this site. Several bis-tetrahydrofuran-containing compounds including darunavir, brecanavir, GS-8374, and Sequoia protease inhibitors exhibit excellent potency against mutant HIV strains that are resistant to clinically used protease inhibitors. The apparently flat resistance profiles of these and some other protease inhibitors may, at least in part, be explained by their high potency against wild-type enzyme. The substrate envelope and solvent-anchoring hypotheses have been used to design and/or rationalize improved resistance profiles. Traditional approaches yielded a lysine sulfonamide PL-100 with a unique resistance profile. SUMMARY Several theories on how to design HIV protease inhibitors with improved resistance profiles have been proposed during the review period. The general concepts that are incorporated into most design strategies include maximizing the interactions with the backbone and conserved side chains of the enzyme while minimizing inhibitor size and maintaining conformational flexibility to allow for modified binding modes.
Collapse
|
16
|
Spiwok V, Tvaroska I. Conformational free energy surface of alpha-N-acetylneuraminic acid: an interplay between hydrogen bonding and solvation. J Phys Chem B 2009; 113:9589-94. [PMID: 19374421 DOI: 10.1021/jp8113495] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The conformational free energy surface of alpha-N-acetylneuraminic acid (Neu5Ac, sialic acid) in the space of ring-puckering coordinates was calculated using the metadynamics method. Free energy surfaces in vacuum and with an explicit solvent were calculated in GLYCAM 06 force field. In vacuum three structures are almost equivalently populated, namely, the (2)C(5) chair and the B(3,6)/(2)S(6) and (O)S(3) boat/skew-boat conformations. The B(3,6)/(2)S(6) structure is stabilized by an ionic hydrogen bond between the amide N-H bond and the carboxylic group. However, this structure is unfavorable in a water environment in which the experimentally observed (2)C(5) chair conformation is predicted to be more stable than the other structures. These results indicate that environment significantly influences conformation of Neu5Ac and that Neu5Ac-processing enzymes might modify a conformation of their substrates solely by a changing polarity of the environment. The structure of Neu5Ac bound in influenza neuraminidase ((4)S(2)/B(2,5)) belongs to conformations preferred in a water environment. The free energy penalty of this conformational change was calculated (relative to (2)C(5)) as 10.2 +/- 2.0 and 17.3 +/- 2.0 kJ/mol for (4,O)B/(O)S(3) and (4)S(2), respectively. This result indicates that mimicking of the enzyme-bound conformation is likely to be a viable strategy for the design of neuraminidase inhibitors.
Collapse
Affiliation(s)
- Vojtech Spiwok
- Department of Structure and Function of Saccharides, Centre for Glycomics, Slovak Academy of Sciences, 84538 Bratislava, Slovak Republic.
| | | |
Collapse
|
17
|
Rubinstein M, Niv MY. Peptidic modulators of protein-protein interactions: progress and challenges in computational design. Biopolymers 2009; 91:505-13. [PMID: 19226619 DOI: 10.1002/bip.21164] [Citation(s) in RCA: 76] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
With the decline in productivity of drug-development efforts, novel approaches to rational drug design are being introduced and developed. Naturally occurring and synthetic peptides are emerging as novel promising compounds that can specifically and efficiently modulate signaling pathways in vitro and in vivo. We describe sequence-based approaches that use peptides to mimic proteins in order to inhibit the interaction of the mimicked protein with its partners. We then discuss a structure-based approach, in which protein-peptide complex structures are used to rationally design and optimize peptidic inhibitors. We survey flexible peptide docking techniques and discuss current challenges and future directions in the rational design of peptidic inhibitors.
Collapse
Affiliation(s)
- Mor Rubinstein
- The Institute of Biochemistry, Food Science and Nutrition, The Hebrew University of Jerusalem, Rehovot 76100, Israel
| | | |
Collapse
|
18
|
Michel C, Laio A, Milet A. Tracing the Entropy along a Reactive Pathway: The Energy As a Generalized Reaction Coordinate. J Chem Theory Comput 2009; 5:2193-6. [DOI: 10.1021/ct900177h] [Citation(s) in RCA: 14] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Carine Michel
- Département de Chimie Moléculaire, Chimie Théorique, UMR-5250, ICMG FR-2607, CNRS, Université Joseph Fourier Grenoble I, DU BP 53 38041 Grenoble Cedex 09 France, SISSA, ISAS, Via Beirut 2-4, I-34014 Trieste, Italy, and Institut Universitaire de France, France
| | - Alessandro Laio
- Département de Chimie Moléculaire, Chimie Théorique, UMR-5250, ICMG FR-2607, CNRS, Université Joseph Fourier Grenoble I, DU BP 53 38041 Grenoble Cedex 09 France, SISSA, ISAS, Via Beirut 2-4, I-34014 Trieste, Italy, and Institut Universitaire de France, France
| | - Anne Milet
- Département de Chimie Moléculaire, Chimie Théorique, UMR-5250, ICMG FR-2607, CNRS, Université Joseph Fourier Grenoble I, DU BP 53 38041 Grenoble Cedex 09 France, SISSA, ISAS, Via Beirut 2-4, I-34014 Trieste, Italy, and Institut Universitaire de France, France
| |
Collapse
|
19
|
Spiwok V, Tvaroska I. Metadynamics modelling of the solvent effect on primary hydroxyl rotamer equilibria in hexopyranosides. Carbohydr Res 2009; 344:1575-81. [PMID: 19539272 DOI: 10.1016/j.carres.2009.05.019] [Citation(s) in RCA: 24] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2009] [Revised: 04/29/2009] [Accepted: 05/13/2009] [Indexed: 11/16/2022]
Abstract
Accurate modelling of rotamer equilibria for the primary hydroxyl groups of monosaccharides continues to be a great challenge of computational glycochemistry. The metadynamics technique was applied to study the conformational free energy surfaces of methyl alpha-D-glucopyranoside and methyl alpha-D-galactopyranoside, employing the GLYCAM06 force field. For both molecules, seven to eight conformational free-energy minima, differing in the omega (O-5-C-5-C-6-O-6) and chi (C-3-C-4-O-4-HO-4) dihedral angles, were identified in vacuum or in a water environment. The calculated rotamer equilibrium of the primary hydroxyl group is significantly different in vacuum than in water. The major effect of a water environment is the destabilisation of a hydrogen bond between O-4-HO-4 and O-6-HO-6 groups. It was possible to calculate the free-energy differences of individual rotamers with an accuracy of better than 2 kJ/mol. The calculated gg, gt and tg rotamer populations in water are in close agreement with experimental measurements, and therefore support the theoretical background of metadynamics.
Collapse
Affiliation(s)
- Vojtech Spiwok
- Department of Biochemistry and Microbiology, Institute of Chemical Technology in Prague, Technická 3, 166 28 Prague 6, Czech Republic.
| | | |
Collapse
|
20
|
Xiu P, Zhou B, Qi W, Lu H, Tu Y, Fang H. Manipulating Biomolecules with Aqueous Liquids Confined within Single-Walled Nanotubes. J Am Chem Soc 2009; 131:2840-5. [DOI: 10.1021/ja804586w] [Citation(s) in RCA: 62] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Affiliation(s)
- Peng Xiu
- School of Physics, Shandong University, Jinan, 250100, China, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, P.O. Box 800-204, Shanghai 201800, China, Graduate School of the Chinese Academy of Sciences, Beijing 100080, China, Department of Physics, Zhejiang Normal University, 321004, Jinhua, China, and Theoretical Physics Center for Science Facilities (TPCSF), CAS, 19(B) Yuquan Road, Beijing 100049, China
| | - Bo Zhou
- School of Physics, Shandong University, Jinan, 250100, China, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, P.O. Box 800-204, Shanghai 201800, China, Graduate School of the Chinese Academy of Sciences, Beijing 100080, China, Department of Physics, Zhejiang Normal University, 321004, Jinhua, China, and Theoretical Physics Center for Science Facilities (TPCSF), CAS, 19(B) Yuquan Road, Beijing 100049, China
| | - Wenpeng Qi
- School of Physics, Shandong University, Jinan, 250100, China, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, P.O. Box 800-204, Shanghai 201800, China, Graduate School of the Chinese Academy of Sciences, Beijing 100080, China, Department of Physics, Zhejiang Normal University, 321004, Jinhua, China, and Theoretical Physics Center for Science Facilities (TPCSF), CAS, 19(B) Yuquan Road, Beijing 100049, China
| | - Hangjun Lu
- School of Physics, Shandong University, Jinan, 250100, China, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, P.O. Box 800-204, Shanghai 201800, China, Graduate School of the Chinese Academy of Sciences, Beijing 100080, China, Department of Physics, Zhejiang Normal University, 321004, Jinhua, China, and Theoretical Physics Center for Science Facilities (TPCSF), CAS, 19(B) Yuquan Road, Beijing 100049, China
| | - Yusong Tu
- School of Physics, Shandong University, Jinan, 250100, China, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, P.O. Box 800-204, Shanghai 201800, China, Graduate School of the Chinese Academy of Sciences, Beijing 100080, China, Department of Physics, Zhejiang Normal University, 321004, Jinhua, China, and Theoretical Physics Center for Science Facilities (TPCSF), CAS, 19(B) Yuquan Road, Beijing 100049, China
| | - Haiping Fang
- School of Physics, Shandong University, Jinan, 250100, China, Shanghai Institute of Applied Physics, Chinese Academy of Sciences, P.O. Box 800-204, Shanghai 201800, China, Graduate School of the Chinese Academy of Sciences, Beijing 100080, China, Department of Physics, Zhejiang Normal University, 321004, Jinhua, China, and Theoretical Physics Center for Science Facilities (TPCSF), CAS, 19(B) Yuquan Road, Beijing 100049, China
| |
Collapse
|
21
|
Caldarini M, Vasile F, Provasi D, Longhi R, Tiana G, Broglia RA. Identification and characterization of folding inhibitors of hen egg lysozyme: an example of a new paradigm of drug design. Proteins 2009; 74:390-9. [PMID: 18623063 DOI: 10.1002/prot.22161] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Abstract
Studies of protein folding indicate the presence of native contacts in the denatured state, giving rise to folding elements which contribute to the accomplishment of the native state. The possibility of finding molecules which can interact with specific folding elements of a target protein preventing it from reaching its native state, and hence from becoming biologically active, is particularly attractive. The notion that folding elements not only provide molecular recognition directing the folding process, but also have conserved sequence, implies that targeting such elements will make protein folding inhibitors less susceptible to mutations which, in many cases, abrogate drug effects. The folding-inhibition strategy can lead to a truly novel and rational approach to drug design, aside from providing new insight into folding. This is illustrated in the case of hen egg lysozyme.
Collapse
Affiliation(s)
- M Caldarini
- Department of Physics, University of Milano and INFN, via Celoria 16, 20133 Milano, Italy
| | | | | | | | | | | |
Collapse
|
22
|
Verkhivker G. Computational proteomics analysis of binding mechanisms and molecular signatures of the HIV-1 protease drugs. Artif Intell Med 2009; 45:197-206. [DOI: 10.1016/j.artmed.2008.08.011] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2007] [Revised: 08/12/2008] [Accepted: 08/19/2008] [Indexed: 11/25/2022]
|
23
|
Moebius K, Eichler J. HIV-derived peptide mimics. DRUG DISCOVERY TODAY. TECHNOLOGIES 2009; 6:e1-e40. [PMID: 24128988 DOI: 10.1016/j.ddtec.2009.09.001] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/02/2023]
|
24
|
Abstract
A molecular dynamics investigation and coarse-grained analysis of inactivated actin-related protein (Arp) 2/3 complex is presented. It was found that the nucleotide binding site within Arp3 remained in a closed position with bound ATP or ADP, but opened when simulation with no nucleotide was performed. In contrast, simulation of the isolated Arp3 subunit with bound ATP, showed a fast opening of the nucleotide binding cleft. A homology model for the missing subdomains 1 and 2 of Arp2 was constructed, and it was also found that the Arp2 binding cleft remained closed with bound nucleotide. Within the nucleotide binding cleft a distinct opening and closing period of 10 ns was observed in many of the simulations of Arp2/3 as well as isolated Arp3. Substitution studies were employed, and several alanine substitutions were found to induce a partial opening of the ATP binding cleft in Arp3 and Arp2, whereas only a single substitution was found to induce opening of the ADP binding cleft. It was also found that the nucleotide type did not cause a substantial change on interfacial contacts between Arp3 and the ArpC2, ArpC3 and ArpC4 subunits. Nucleotide-free Arp3 had generally less stable contacts, but the overall contact architecture was constant. Finally, nucleotide-dependent coarse-grained models for Arp3 are developed that serve to further highlight the structural differences induced in Arp3 by nucleotide hydrolysis.
Collapse
|
25
|
Quigley D, Rodger PM. Metadynamics simulations of ice nucleation and growth. J Chem Phys 2008; 128:154518. [DOI: 10.1063/1.2888999] [Citation(s) in RCA: 112] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
Abstract
Biochemical experiments have recently revealed that the p-S8 peptide, with an amino-acid sequence identical to the conserved fragment 83-93 (S8) of the HIV-1 protease, can inhibit catalytic activity of the enzyme by interfering with protease folding and dimerization. In this study, we introduce a hierarchical modeling approach for understanding the molecular basis of the HIV-1 protease folding inhibition. Coarse-grained molecular docking simulations of the flexible p-S8 peptide with the ensembles of HIV-1 protease monomers have revealed structurally different complexes of the p-S8 peptide, which can be formed by targeting the conserved segment 24-34 (S2) of the folding nucleus (folding inhibition) and by interacting with the antiparallel termini beta-sheet region (dimerization inhibition). All-atom molecular dynamics simulations of the inhibitor complexes with the HIV-1 PR monomer have been independently carried out for the predicted folding and dimerization binding modes of the p-S8 peptide, confirming the thermodynamic stability of these complexes. Binding free-energy calculations of the p-S8 peptide and its active analogs are then performed using molecular dynamics trajectories of the peptide complexes with the HIV-1 PR monomers. The results of this study have provided a plausible molecular model for the inhibitor intervention with the HIV-1 PR folding and dimerization and have accurately reproduced the experimental inhibition profiles of the active folding inhibitors.
Collapse
|
27
|
Broglia RA, Levy Y, Tiana G. HIV-1 protease folding and the design of drugs which do not create resistance. Curr Opin Struct Biol 2008; 18:60-6. [DOI: 10.1016/j.sbi.2007.10.004] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2007] [Accepted: 10/29/2007] [Indexed: 10/22/2022]
|
28
|
Camilloni C, Provasi D, Tiana G, Broglia RA. Exploring the protein G helix free-energy surface by solute tempering metadynamics. Proteins 2007; 71:1647-54. [DOI: 10.1002/prot.21852] [Citation(s) in RCA: 58] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
|