1
|
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.3] [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
|
2
|
Fan F, Zhao Y, Cao Z. Insight into the delivery channel and selectivity of multiple binding sites in bovine serum albumin towards naphthalimide-polyamine derivatives. Phys Chem Chem Phys 2019; 21:7429-7439. [PMID: 30892331 DOI: 10.1039/c9cp00527g] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Naphthalimide derivatives are types of small-molecule anticancer drug candidates; however, their negative factors and potential side effects make their application limited. The pharmacophores select a direct access into the tumor cells as the first choice; this can reduce the side effect of the anti-cancer drugs on the normal cells. Herein, the delivery and binding of the naphthalimide-polyamine complex assisted by the bovine serum albumin (BSA) protein have been studied by combining several molecular dynamic simulations. The plausible transportation channels and the most favorable pathways for the delivery of the naphthalimide-polyamine complex to two drug sites (DSI and DSII), their thermodynamic and dynamic properties and the mechanisms have been discussed in detail. The residues His287 and Phe394 acted as guards in the DSI and DSII, respectively, which played a gating-switch role by flipping the ring from open to close during the compound delivery. The binding mode, binding energy and substituent effects have been also identified. The two drug sites have different preferences towards the compound with the electron-withdrawing and electron-donating substituents, and their strong interactions are more sensitive to the number of the substituent groups. The naphthalimide-polyamine complexes are more likely to choose DSI, both thermodynamically and dynamically, as compared to DSII. This selective specificity of these two drug sites manipulated by the electron-withdrawing and electron-donating substituents is quite promising for the design of new naphthalimide drugs.
Collapse
Affiliation(s)
- Fangfang Fan
- State Key Laboratory of Physical Chemistry of Solid Surfaces and Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, College of Chemistry and Chemical Engineering, Xiamen University, Xiamen, China.
| | | | | |
Collapse
|