1
|
Dey P, Biswas P. Relaxation dynamics measure the aggregation propensity of amyloid-β and its mutants. J Chem Phys 2023; 158:105101. [PMID: 36922119 DOI: 10.1063/5.0138189] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/17/2023] Open
Abstract
Atomistic molecular dynamics simulations are employed to investigate the global and segmental relaxation dynamics of the amyloid-β protein and its causative and protective mutants. Amyloid-β exhibits significant global/local dynamics that span a broad range of length and time scales due to its intrinsically disordered nature. The relaxation dynamics of the amyloid-β protein and its mutants is quantitatively correlated with its experimentally measured aggregation propensity. The protective mutant has slower relaxation dynamics, whereas the causative mutants exhibit faster global dynamics compared with that of the wild-type amyloid-β. The local dynamics of the amyloid-β protein or its mutants is governed by a complex interplay of the charge, hydrophobicity, and change in the molecular mass of the mutated residue.
Collapse
Affiliation(s)
- Priya Dey
- Department of Chemistry, University of Delhi, Delhi 110007, India
| | - Parbati Biswas
- Department of Chemistry, University of Delhi, Delhi 110007, India
| |
Collapse
|
2
|
Oliva R, Winter R. Harnessing Pressure-Axis Experiments to Explore Volume Fluctuations, Conformational Substates, and Solvation of Biomolecular Systems. J Phys Chem Lett 2022; 13:12099-12115. [PMID: 36546666 DOI: 10.1021/acs.jpclett.2c03186] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
Intrinsic thermodynamic fluctuations within biomolecules are crucial for their function, and flexibility is one of the strategies that evolution has developed to adapt to extreme environments. In this regard, pressure perturbation is an important tool for mechanistically exploring the causes and effects of volume fluctuations in biomolecules and biomolecular assemblies, their role in biomolecular interactions and reactions, and how they are affected by the solvent properties. High hydrostatic pressure is also a key parameter in the context of deep-sea and subsurface biology and the study of the origin and physical limits of life. We discuss the role of pressure-axis experiments in revealing intrinsic structural fluctuations as well as high-energy conformational substates of proteins and other biomolecular systems that are important for their function and provide some illustrative examples. We show that the structural and dynamic information obtained from such pressure-axis studies improves our understanding of biomolecular function, disease, biological evolution, and adaptation.
Collapse
Affiliation(s)
- Rosario Oliva
- Department of Chemistry and Chemical Biology, Physical Chemistry I, Biophysical Chemistry, TU Dortmund University, Otto-Hahn-Strasse 6, Dortmund44227, Germany
- Department of Chemical Sciences, University of Naples Federico II, Via Cintia 4, 80126Naples, Italy
| | - Roland Winter
- Department of Chemistry and Chemical Biology, Physical Chemistry I, Biophysical Chemistry, TU Dortmund University, Otto-Hahn-Strasse 6, Dortmund44227, Germany
| |
Collapse
|
3
|
Alvarado YJ, Olivarez Y, Lossada C, Vera-Villalobos J, Paz JL, Vera E, Loroño M, Vivas A, Torres FJ, Jeffreys LN, Hurtado-León ML, González-Paz L. Interaction of the new inhibitor paxlovid (PF-07321332) and ivermectin with the monomer of the main protease SARS-CoV-2: A volumetric study based on molecular dynamics, elastic networks, classical thermodynamics and SPT. Comput Biol Chem 2022; 99:107692. [PMID: 35640480 PMCID: PMC9107165 DOI: 10.1016/j.compbiolchem.2022.107692] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Revised: 04/28/2022] [Accepted: 05/02/2022] [Indexed: 02/04/2023]
Abstract
The COVID-19 pandemic has accelerated the study of drugs, most notably ivermectin and more recently Paxlovid (PF-07321332) which is in phase III clinical trials with experimental data showing covalent binding to the viral protease Mpro. Theoretical developments of catalytic site-directed docking support thermodynamically feasible non-covalent binding to Mpro. Here we show that Paxlovid binds non-covalently at regions other than the catalytic sites with energies stronger than reported and at the same binding site as the ivermectin B1a homologue, all through theoretical methodologies, including blind docking. We volumetrically characterize the non-covalent interaction of the ivermectin homologues (avermectins B1a and B1b) and Paxlovid with the mMpro monomer, through molecular dynamics and scaled particle theory (SPT). Using the fluctuation-dissipation theorem (FDT), we estimated the electric dipole moment fluctuations at the surface of each of complex involved in this study, with similar trends to that observed in the interaction volume. Using fluctuations of the intrinsic volume and the number of flexible fragments of proteins using anisotropic and Gaussian elastic networks (ANM+GNM) suggests the complexes with ivermectin are more dynamic and flexible than the unbound monomer. In contrast, the binding of Paxlovid to mMpro shows that the mMpro-PF complex is the least structurally dynamic of all the species measured in this investigation. The results support a differential molecular mechanism of the ivermectin and PF homologues in the mMpro monomer. Finally, the results showed that Paxlovid despite beingbound in different sites through covalent or non-covalent forms behaves similarly in terms of its structural flexibility and volumetric behaviour.
Collapse
Affiliation(s)
- Ysaias José Alvarado
- Instituto Venezolano de Investigaciones Científicas (IVIC), Centro de Investigación y Tecnología de Materiales (CITeMA), Laboratorio de Caracterización Molecular y Biomolecular, 4001 Maracaibo, Bolivarian Republic of Venezuela.
| | - Yosmari Olivarez
- Universidad del Zulia (LUZ). Facultad Experimental de Ciencias (FEC), Departamento de Quimica, Laboratorio de Electronica Molecular, 4001 Maracaibo, Bolivarian Republic of Venezuela
| | - Carla Lossada
- Instituto Venezolano de Investigaciones Científicas (IVIC), Centro de Investigación y Tecnología de Materiales (CITeMA), Laboratorio de Caracterización Molecular y Biomolecular, 4001 Maracaibo, Bolivarian Republic of Venezuela
| | - Joan Vera-Villalobos
- Facultad de Ciencias Naturales y Matemáticas, Departamento de Química y Ciencias Ambientales, Laboratorio de Análisis Químico Instrumental (LAQUINS), Escuela Superior Politécnica del Litoral, Guayaquil, Ecuador
| | - José Luis Paz
- Departamento Académico de Química Inorgánica, Facultad de Química e Ingeniería Química, Universidad Nacional Mayor de San Marcos, Lima, Peru
| | - Eddy Vera
- Universidad del Zulia (LUZ). Facultad Experimental de Ciencias (FEC), Departamento de Quimica, Laboratorio de Electronica Molecular, 4001 Maracaibo, Bolivarian Republic of Venezuela
| | - Marcos Loroño
- Departamento Académico de Química Analítica e Instrumental, Facultad de Química e Ingeniería Química, Universidad Nacional Mayor de San Marcos, Lima, Peru
| | - Alejandro Vivas
- Universidad del Zulia (LUZ). Facultad Experimental de Ciencias (FEC), Departamento de Quimica, Laboratorio de Electronica Molecular, 4001 Maracaibo, Bolivarian Republic of Venezuela
| | - Fernando Javier Torres
- Grupo de Química Computacional y Teórica (QCT-UR), Facultad de Ciencias Naturales, Universidad del Rosario, Bogotá, Colombia; Grupo de Química Computacional y Teórica (QCT-USFQ), Instituto de Simulación Computacional (ISC-USFQ), Departamento de Ingeniería Química, Universidad San Francisco de Quito (USFQ), Quito, Ecuador
| | - Laura N Jeffreys
- Centre for Drugs and Diagnostics, Department of Tropical Disease Biology, Liverpool School of Tropical Medicine, Pembroke Place, Liverpool L3 5QA, UK
| | - María Laura Hurtado-León
- Universidad del Zulia (LUZ), Facultad Experimental de Ciencias (FEC), Departamento de Biología, Laboratorio de Genética y Biología Molecular (LGBM), Maracaibo 4001, Zulia, Bolivarian Republic of Venezuela
| | - Lenin González-Paz
- Universidad del Zulia (LUZ), Facultad Experimental de Ciencias (FEC), Departamento de Biología, Laboratorio de Genética y Biología Molecular (LGBM), Maracaibo 4001, Zulia, Bolivarian Republic of Venezuela; Instituto Venezolano de Investigaciones Científicas (IVIC), Centro de Estudios Botanicos y Agroforestales, (CEBA), Laboratorio de Proteccion Vegetal, 4001 Maracaibo, Bolivarian Republic of Venezuela.
| |
Collapse
|
4
|
Aggarwal L, Biswas P. Hydration Thermodynamics of Familial Parkinson's Disease-Linked Mutants of α-Synuclein. J Chem Inf Model 2021; 61:1850-1858. [PMID: 33749266 DOI: 10.1021/acs.jcim.1c00034] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
The hydration thermodynamics of different mutants of α-synuclein (α-syn) related to familial Parkinson's disease (PD) is explored using a computational approach that combines both molecular dynamics simulations in water and integral equation theory of molecular liquids. This analysis focuses on the change in conformational entropy, hydration free energy (HFE), and partial molar volume of α-syn upon mutation. The results show that A53T, A30P, E46K, and H50Q mutants aggregate more readily and display increased HFE and less negative interaction volume than the wild-type α-syn. In contrast, an opposite trend is observed for the G51D mutant with a lower experimental aggregation rate. The residuewise decomposition analysis of the HFE highlights that the dehydration/hydration of the hydrophilic residue-rich N- and C-termini of α-syn majorly contributes to the change upon mutation. The hydration shell contributions of different residues to the interaction volume are consistent with its increase/decrease upon mutation. This work shows that both HFE and interaction volume determine the aggregation kinetics of α-syn upon mutation and may serve as an appropriate benchmark for the treatment of PD.
Collapse
Affiliation(s)
- Leena Aggarwal
- Department of Chemistry, University of Delhi, Delhi 110007, India
| | - Parbati Biswas
- Department of Chemistry, University of Delhi, Delhi 110007, India
| |
Collapse
|
5
|
Aggarwal L, Biswas P. Hydration Thermodynamics of the N-Terminal FAD Mutants of Amyloid-β. J Chem Inf Model 2021; 61:298-310. [PMID: 33440932 DOI: 10.1021/acs.jcim.0c01286] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The hydration thermodynamics of amyloid-β (Aβ) and its pathogenic familial Alzheimer's disease (FAD) mutants such as A2V, Taiwan (D7H), Tottori (D7N), and English (H6R) and the protective A2T mutant is investigated by a combination of all-atom, explicit water molecular dynamics (MD) simulations and the three-dimensional reference interaction site model (3D-RISM) theory. The change in the hydration free energy on mutation is decomposed into the energetic and entropic components, which comprise electrostatic and nonelectrostatic contributions. An increase in the hydration free energy is observed for A2V, D7H, D7N, and H6R mutations that increase the aggregation propensity of Aβ and lead to an early onset of Alzheimer's disease, while a reverse trend is noted for the protective A2T mutation. An antiphase correlation is found between the change in the hydration energy and the internal energy of Aβ upon mutation. A residue-wise decomposition analysis shows that the change in the hydration free energy of Aβ on mutation is primarily due to the hydration/dehydration of the side-chain atoms of the negatively charged residues. The decrease in the hydration of the negatively charged residues on mutation may decrease the solubility of the mutant, which increases the observed aggregation propensity of the FAD mutants. Results obtained from the theory show an excellent match with the experimentally reported data.
Collapse
Affiliation(s)
- Leena Aggarwal
- Department of Chemistry, University of Delhi, Delhi 110007, India
| | - Parbati Biswas
- Department of Chemistry, University of Delhi, Delhi 110007, India
| |
Collapse
|