1
|
Zanetti-Polzi L, Daidone I, Iacobucci C, Amadei A. Thermodynamic Evolution of a Metamorphic Protein: A Theoretical-Computational Study of Human Lymphotactin. Protein J 2023:10.1007/s10930-023-10123-7. [PMID: 37233895 DOI: 10.1007/s10930-023-10123-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/03/2023] [Indexed: 05/27/2023]
Abstract
Metamorphic, or fold-switching, proteins feature different folds that are physiologically relevant. The human chemokine XCL1 (or Lymphotactin) is a metamorphic protein that features two native states, an [Formula: see text] and an all[Formula: see text] fold, which have similar stability at physiological condition. Here, extended molecular dynamics (MD) simulations, principal component analysis of atomic fluctuations and thermodynamic modeling based on both the configurational volume and free energy landscape, are used to obtain a detailed characterization of the conformational thermodynamics of human Lymphotactin and of one of its ancestors (as was previously obtained by genetic reconstruction). Comparison of our computational results with the available experimental data show that the MD-based thermodynamics can explain the experimentally observed variation of the conformational equilibrium between the two proteins. In particular, our computational data provide an interpretation of the thermodynamic evolution in this protein, revealing the relevance of the configurational entropy and of the shape of the free energy landscape within the essential space (i.e., the space defined by the generalized internal coordinates providing the largest, typically non-Gaussian, structural fluctuations).
Collapse
Affiliation(s)
- Laura Zanetti-Polzi
- Center S3, CNR-Institute of Nanoscience, Via Campi 213/A, 100190, Modena, Italy
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, Via Vetoio (Coppito 1), 67010, L'Aquila, Italy
| | - Claudio Iacobucci
- Department of Physical and Chemical Sciences, University of L'Aquila, Via Vetoio (Coppito 1), 67010, L'Aquila, Italy
| | - Andrea Amadei
- Department of Chemical Science and Technology, University of Rome "Tor Vergata", Via Della Ricerca Scientifica 1, 00185, Rome, Italy.
| |
Collapse
|
2
|
Capone M, Zanetti-Polzi L, Leonzi I, Spreti N, Daidone I. Evidence for a high pK a of an aspartic acid residue in the active site of CALB by a fully atomistic multiscale approach. J Biomol Struct Dyn 2022:1-8. [PMID: 35593533 DOI: 10.1080/07391102.2022.2077834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Abstract
Candida antarctica Lipase B (CALB) is a paradigm for the family of lipases. At pH 7, the optimal pH for catalysis, the protonation state of an aspartic acid of the active site (Asp134) could not be conclusively assigned. In fact, the pKa estimate provided by a widely used computational tool, namely PropKa, that predicts pKa values of ionizable groups in proteins based on the crystallographic structure, is only slightly above 7 (pKa = 7.25). This, along with the lack of an experimental evaluation, makes the assignment of its protonation state at neutral pH challenging. Here, we calculate the pKa of Asp134 by means of a fully atomistic multiscale computational approach based on classical molecular dynamics (MD) simulation and the perturbed matrix method (PMM), namely the MD-PMM approach. MD-PMM is able to take into account the dynamics of the system and, at the same time, to treat the deprotonation step at the quantum level. The calculations provide a pKa value of 8.9 ± 1.1, hence suggesting that Asp134 in CALB should be protonated at neutral, and even at slightly basic, pH.
Collapse
Affiliation(s)
- Matteo Capone
- Department of Physical and Chemical Sciences, University of L'Aquila, L'Aquila, Italy
| | | | - Ilenia Leonzi
- Department of Physical and Chemical Sciences, University of L'Aquila, L'Aquila, Italy
| | - Nicoletta Spreti
- Department of Physical and Chemical Sciences, University of L'Aquila, L'Aquila, Italy
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, L'Aquila, Italy
| |
Collapse
|
3
|
Pinto SMV, Tasinato N, Barone V, Zanetti-Polzi L, Daidone I. A computational insight into the relationship between side chain IR line shapes and local environment in fibril-like structures. J Chem Phys 2021; 154:084105. [PMID: 33639764 DOI: 10.1063/5.0038913] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Infrared spectroscopy is a widely used technique to characterize protein structures and protein mediated processes. While the amide I band provides information on proteins' secondary structure, amino acid side chains are used as infrared probes for the investigation of protein reactions and local properties. In this paper, we use a hybrid quantum mechanical/classical molecular dynamical approach based on the perturbed matrix method to compute the infrared band due to the C=O stretching mode of amide-containing side chains. We calculate, at first, the infrared band of zwitterionic glutamine in water and obtain results in very good agreement with the experimental data. Then, we compute the signal arising from glutamine side chains in a microcrystal of the yeast prion Sup35-derived peptide, GNNQQNY, with a fibrillar structure. The infrared bands obtained by selective isotopic labeling of the two glutamine residues, Q4 and Q5, of each peptide were experimentally used to investigate the local hydration in the fibrillar microcrystal. The experimental spectra of the two glutamine residues, which experience different hydration environments, feature different spectral signals that are well reproduced by the corresponding calculated spectra. In addition, the analysis of the simulated spectra clarifies the molecular origin of the experimentally observed spectroscopic differences that arise from the different local electric field experienced by the two glutamine residues, which is, in turn, determined by a different hydrogen bonding pattern.
Collapse
Affiliation(s)
- Sandra M V Pinto
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Nicola Tasinato
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore, Piazza dei Cavalieri 7, I-56126 Pisa, Italy
| | | | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, Via Vetoio, I-67100 L'Aquila, Italy
| |
Collapse
|
4
|
Vieira Pinto SM, Tasinato N, Barone V, Amadei A, Zanetti-Polzi L, Daidone I. Modeling amino-acid side chain infrared spectra: the case of carboxylic residues. Phys Chem Chem Phys 2020; 22:3008-3016. [DOI: 10.1039/c9cp04774c] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Infrared (IR) spectroscopy is commonly utilized for the investigation of protein structures and protein-mediated processes.
Collapse
Affiliation(s)
- Sandra Mónica Vieira Pinto
- Scuola Normale Superiore
- I-56126 Pisa
- Italy
- Department of Physical and Chemical Sciences
- University of L'Aquila
| | | | | | - Andrea Amadei
- Department of Chemical and Technological Sciences
- University of Rome “Tor Vergata
- I-00185 Rome
- Italy
| | - Laura Zanetti-Polzi
- Department of Physical and Chemical Sciences
- University of L'Aquila
- I-67010 L'Aquila
- Italy
- CNR Institute of Nanoscience
| | - Isabella Daidone
- Department of Physical and Chemical Sciences
- University of L'Aquila
- I-67010 L'Aquila
- Italy
| |
Collapse
|
5
|
Amadei A, Aschi M. Theoretical-computational modeling of charge transfer and intersystem crossing reactions in complex chemical systems. RSC Adv 2018; 8:27900-27918. [PMID: 35542751 PMCID: PMC9083445 DOI: 10.1039/c8ra03900c] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2018] [Accepted: 07/23/2018] [Indexed: 11/21/2022] Open
Abstract
In this paper we present a theoretical-computational methodology specifically aimed at describing processes involving internal conversion or intersystem crossing, from atomistic (semiclassical) simulations and, hence, very suitable for treating complex atomic-molecular systems. The core of the presented approach is the evaluation of the diabatic perturbed energy surfaces of a portion of the whole system, treated at the quantum level and therefore preventively selected, in semi-classical interaction with the atomic-molecular environment. Subsequently, the estimation of the coupling between the diabatic surfaces and the inclusion of the obtained observables within a properly designed kinetic model allows the reconstruction of the whole phenomenology directly comparable to the experimental (typically kinetic) data. Application to two systems has demonstrated that the proposed approach can represent a valuable tool, somewhat complementary to other methods based on explicit quantum-dynamical approaches, for the theoretical-computational investigations of large and complex atomic-molecular systems.
Collapse
Affiliation(s)
- Andrea Amadei
- Dipartimento di Scienze e Tecnologie Chimiche, Universita' di Roma 'Tor Vergata' Roma Italy +390672594905
| | - Massimiliano Aschi
- Dipartimento di Scienze Fisiche e Chimiche, Universita' di L'Aquila L'Aquila Italy +390862433775
| |
Collapse
|
6
|
Daidone I, Amadei A, Aschi M, Zanetti-Polzi L. On the nature of solvatochromic effect: The riboflavin absorption spectrum as a case study. SPECTROCHIMICA ACTA. PART A, MOLECULAR AND BIOMOLECULAR SPECTROSCOPY 2018; 192:451-457. [PMID: 29212059 DOI: 10.1016/j.saa.2017.11.031] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2017] [Revised: 10/27/2017] [Accepted: 11/14/2017] [Indexed: 06/07/2023]
Abstract
We present here the calculation of the absorption spectrum of riboflavin in acetonitrile and dimethyl sulfoxide using a hybrid quantum/classical approach, namely the perturbed matrix method, based on quantum mechanical calculations and molecular dynamics simulations. The calculated spectra are compared to the absorption spectrum of riboflavin previously calculated in water and to the experimental spectra obtained in all three solvents. The experimentally observed variations in the absorption spectra upon change of the solvent environment are well reproduced by the calculated spectra. In addition, the nature of the excited states of riboflavin interacting with different solvents is investigated, showing that environment effects determine a recombination of the gas-phase electronic states and that such a recombination is strongly affected by the polarity of the solvent inducing significant changes in the absorption spectra.
Collapse
Affiliation(s)
- Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy
| | - Andrea Amadei
- Department of Chemical and Technological Sciences, University of Rome "Tor Vergata", Via della Ricerca Scientifica, 00185 Rome, Italy
| | - Massimiliano Aschi
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy
| | - Laura Zanetti-Polzi
- Department of Physical and Chemical Sciences, University of L'Aquila, via Vetoio (Coppito 1), L'Aquila 67010, Italy.
| |
Collapse
|
7
|
Carrillo-Parramon O, Del Galdo S, Aschi M, Mancini G, Amadei A, Barone V. Flexible and Comprehensive Implementation of MD-PMM Approach in a General and Robust Code. J Chem Theory Comput 2017; 13:5506-5514. [DOI: 10.1021/acs.jctc.7b00341] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Affiliation(s)
| | - Sara Del Galdo
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7 I-56126, Pisa, Italy
| | - Massimiliano Aschi
- Dipartimento
di Scienze Fisiche e Chimiche, Universitá di L’Aquila, Via
Vetoio s.n.c.67100, L’Aquila, Italy
| | - Giordano Mancini
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7 I-56126, Pisa, Italy
- Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| | - Andrea Amadei
- Dipartimento
di Scienze e Tecnologie Chimiche, Universitá di Roma ’Tor Vergata’, Via Della Ricerca Scientifica, 00100 Roma, Italy
| | - Vincenzo Barone
- Scuola Normale Superiore di Pisa, Piazza dei Cavalieri 7 I-56126, Pisa, Italy
- Istituto Nazionale di Fisica Nucleare (INFN) sezione di Pisa, Largo Bruno Pontecorvo 3, 56127 Pisa, Italy
| |
Collapse
|
8
|
Zanetti-Polzi L, Davis CM, Gruebele M, Dyer RB, Amadei A, Daidone I. Parallel folding pathways of Fip35 WW domain explained by infrared spectra and their computer simulation. FEBS Lett 2017; 591:3265-3275. [PMID: 28881468 DOI: 10.1002/1873-3468.12836] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2017] [Revised: 08/21/2017] [Accepted: 08/31/2017] [Indexed: 11/06/2022]
Abstract
We present a calculation of the amide I' infrared (IR) spectra of the folded, unfolded, and intermediate states of the WW domain Fip35, a model system for β-sheet folding. Using an all-atom molecular dynamics simulation in which multiple folding and unfolding events take place we identify six conformational states and then apply perturbed matrix method quantum-mechanical calculations to determine their amide I' IR spectra. Our analysis focuses on two states previously identified as Fip35 folding intermediates and suggests that a three-stranded core similar to the folded state core is the main source of the spectroscopic differences between the two intermediates. In particular, we propose a hypothesis for why folding via one of these intermediates was not experimentally observed by IR T-jump.
Collapse
Affiliation(s)
| | - Caitlin M Davis
- Departments of Chemistry and Physics, University of Illinois at Urbana-Champaign, IL, USA
| | - Martin Gruebele
- Departments of Chemistry and Physics, University of Illinois at Urbana-Champaign, IL, USA.,Center for Biophysics and Quantitative Biology, University of Illinois at Urbana-Champaign, IL, USA
| | - R Brian Dyer
- Department of Chemistry, Emory University, Atlanta, GA, USA
| | - Andrea Amadei
- Department of Chemical and Technological Sciences, University of Rome "Tor Vergata", Italy
| | - Isabella Daidone
- Department of Physical and Chemical Sciences, University of L'Aquila, Italy
| |
Collapse
|
9
|
Influence of Glu/Arg, Asp/Arg, and Glu/Lys Salt Bridges on α-Helical Stability and Folding Kinetics. Biophys J 2017; 110:2328-2341. [PMID: 27276251 DOI: 10.1016/j.bpj.2016.04.015] [Citation(s) in RCA: 52] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2015] [Revised: 03/17/2016] [Accepted: 04/08/2016] [Indexed: 01/24/2023] Open
Abstract
Using a combination of ultraviolet circular dichroism, temperature-jump transient-infrared spectroscopy, and molecular dynamics simulations, we investigate the effect of salt bridges between different types of charged amino-acid residue pairs on α-helix folding. We determine the stability and the folding and unfolding rates of 12 alanine-based α-helical peptides, each of which has a nearly identical composition containing three pairs of positively and negatively charged residues (either Glu(-)/Arg(+), Asp(-)/Arg(+), or Glu(-)/Lys(+)). Within each set of peptides, the distance and order of the oppositely charged residues in the peptide sequence differ, such that they have different capabilities of forming salt bridges. Our results indicate that stabilizing salt bridges (in which the interacting residues are spaced and ordered such that they favor helix formation) speed up α-helix formation by up to 50% and slow down the unfolding of the α-helix, whereas salt bridges with an unfavorable geometry have the opposite effect. Comparing the peptides with different types of charge pairs, we observe that salt bridges between side chains of Glu(-) and Arg(+) are most favorable for the speed of folding, probably because of the larger conformational space of the salt-bridging Glu(-)/Arg(+) rotamer pairs compared to Asp(-)/Arg(+) and Glu(-)/Lys(+). We speculate that the observed impact of salt bridges on the folding kinetics might explain why some proteins contain salt bridges that do not stabilize the final, folded conformation.
Collapse
|
10
|
Computational evidence support the hypothesis of neuroglobin also acting as an electron transfer species. J Biol Inorg Chem 2017; 22:615-623. [DOI: 10.1007/s00775-017-1455-2] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2016] [Accepted: 03/28/2017] [Indexed: 12/31/2022]
|
11
|
Kitzig S, Thilemann M, Cordes T, Rück-Braun K. Light-Switchable Peptides with a Hemithioindigo Unit: Peptide Design, Photochromism, and Optical Spectroscopy. Chemphyschem 2016; 17:1252-63. [DOI: 10.1002/cphc.201501050] [Citation(s) in RCA: 64] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2015] [Indexed: 12/21/2022]
Affiliation(s)
- S. Kitzig
- Institut für Chemie; Technische Universität Berlin; Str. des 17. Juni 135 10623 Berlin Germany
| | - M. Thilemann
- Institut für Chemie; Technische Universität Berlin; Str. des 17. Juni 135 10623 Berlin Germany
| | - T. Cordes
- Molecular Microscopy Research Group; Zernike Institute for Advanced Materials; University of Groningen; Nijenborgh 4 9747 AG Groningen The Netherlands
| | - Karola Rück-Braun
- Institut für Chemie; Technische Universität Berlin; Str. des 17. Juni 135 10623 Berlin Germany
| |
Collapse
|
12
|
Daidone I, Thukral L, Smith JC, Amadei A. Monitoring the Folding Kinetics of a β-Hairpin by Time-Resolved IR Spectroscopy in Silico. J Phys Chem B 2015; 119:4849-56. [DOI: 10.1021/acs.jpcb.5b01477] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Isabella Daidone
- Department
of Physical and Chemical Sciences, University of L’Aquila, via
Vetoio (Coppito 1), 67010 L’Aquila, Italy
| | - Lipi Thukral
- CSIR-Institute of Genomics and Integrative Biology, South Campus, Mathura Road, New Delhi 110020, India
| | - Jeremy C. Smith
- University of Tennessee/Oak Ridge National Laboratory, Center for Molecular Biophysics, P.O. Box 2008, Oak Ridge, Tennessee 37831-6309, United States
- Department
of Biochemistry and Cellular and Molecular Biology, University of Tennessee, M407 Walters Life Sciences, 1414 Cumberland Avenue, Knoxville, Tennessee 37996, United States
| | - Andrea Amadei
- Department
of Chemical Sciences and Technologies, University of Rome “Tor Vergata”, via della Ricerca Scientifica 1, 00133 Rome, Italy
| |
Collapse
|
13
|
Schmitt T, de Groot FMF, Rubensson JE. Prospects of high-resolution resonant X-ray inelastic scattering studies on solid materials, liquids and gases at diffraction-limited storage rings. JOURNAL OF SYNCHROTRON RADIATION 2014; 21:1065-76. [PMID: 25177995 PMCID: PMC4151682 DOI: 10.1107/s1600577514017123] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/24/2014] [Accepted: 07/24/2014] [Indexed: 05/29/2023]
Abstract
The spectroscopic technique of resonant inelastic X-ray scattering (RIXS) will particularly profit from immensely improved brilliance of diffraction-limited storage rings (DLSRs). In RIXS one measures the intensities of excitations as a function of energy and momentum transfer. DLSRs will allow for pushing the achievable energy resolution, signal intensity and the sampled spot size to new limits. With RIXS one nowadays probes a broad range of electronic systems reaching from simple molecules to complex materials displaying phenomena like peculiar magnetism, two-dimensional electron gases, superconductivity, photovoltaic energy conversion and heterogeneous catalysis. In this article the types of improved RIXS studies that will become possible with X-ray beams from DLSRs are envisioned.
Collapse
Affiliation(s)
- Thorsten Schmitt
- Research Department Synchrotron Radiation and Nanotechnology, Paul Scherrer Institut, Swiss Light Source, WSLA/123, 5232 Villigen PSI, Switzerland
| | - Frank M. F. de Groot
- Inorganic Chemistry and Catalysis, Debye Institute for Nanomaterials Science, Utrecht University, Universiteitsweg 99, Utrecht 3584 CG, The Netherlands
| | - Jan-Erik Rubensson
- Department of Physics and Astronomy, Uppsala University, Box 516, Uppsala 751 20, Sweden
| |
Collapse
|