1
|
Gilson MK, Kurtzman T. Free Energy Density of a Fluid and Its Role in Solvation and Binding. J Chem Theory Comput 2024; 20:2871-2887. [PMID: 38536144 PMCID: PMC11197885 DOI: 10.1021/acs.jctc.3c01173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/10/2024]
Abstract
The concept that a fluid has a position-dependent free energy density appears in the literature but has not been fully developed or accepted. We set this concept on an unambiguous theoretical footing via the following strategy. First, we set forth four desiderata that should be satisfied by any definition of the position-dependent free energy density, f(R), in a system comprising only a fluid and a rigid solute: its volume integral, plus the fixed internal energy of the solute, should be the system free energy; it deviates from its bulk value, fbulk, near a solute but should asymptotically approach fbulk with increasing distance from the solute; it should go to zero where the solvent density goes to zero; and it should be well-defined in the most general case of a fluid made up of flexible molecules with an arbitrary interaction potential. Second, we use statistical thermodynamics to formulate a definition of the free energy density that satisfies these desiderata. Third, we show how any free energy density satisfying the desiderata may be used to analyze molecular processes in solution. In particular, because the spatial integral of f(R) equals the free energy of the system, it can be used to compute free energy changes that result from the rearrangement of solutes as well as the forces exerted on the solutes by the solvent. This enables the use of a thermodynamic analysis of water in protein binding sites to inform ligand design. Finally, we discuss related literature and address published concerns regarding the thermodynamic plausibility of a position-dependent free energy density. The theory presented here has applications in theoretical and computational chemistry and may be further generalizable beyond fluids, such as to solids and macromolecules.
Collapse
Affiliation(s)
- Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, and Department of Chemistry and Biochemistry, UC San Diego, La Jolla, CA, 92093, USA
| | - Tom Kurtzman
- PhD Programs in Chemistry, Biochemistry, and Biology, The Graduate Center of the City University of New York, New York, 10016, USA; Department of Chemistry, Lehman College, The City University of New York, Bronx, New York, 10468, USA
| |
Collapse
|
2
|
Carney DW, Leffler AE, Bell JA, Chandrasinghe AS, Cheng C, Chang E, Dornford A, Dougan DR, Frye LL, Grimes ME, Knehans T, Knight JL, Komandla M, Lane W, Li H, Newman SR, Phimister K, Saikatendu KS, Silverstein H, Vafaei S. Exploiting high-energy hydration sites for the discovery of potent peptide aldehyde inhibitors of the SARS-CoV-2 main protease with cellular antiviral activity. Bioorg Med Chem 2024; 103:117577. [PMID: 38518735 DOI: 10.1016/j.bmc.2023.117577] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Revised: 12/21/2023] [Accepted: 12/22/2023] [Indexed: 03/24/2024]
Abstract
Small-molecule antivirals that prevent the replication of the SARS-CoV-2 virus by blocking the enzymatic activity of its main protease (Mpro) are and will be a tenet of pandemic preparedness. However, the peptidic nature of such compounds often precludes the design of compounds within favorable physical property ranges, limiting cellular activity. Here we describe the discovery of peptide aldehyde Mpro inhibitors with potent enzymatic and cellular antiviral activity. This structure-activity relationship (SAR) exploration was guided by the use of calculated hydration site thermodynamic maps (WaterMap) to drive potency via displacement of waters from high-energy sites. Thousands of diverse compounds were designed to target these high-energy hydration sites and then prioritized for synthesis by physics- and structure-based Free-Energy Perturbation (FEP+) simulations, which accurately predicted biochemical potencies. This approach ultimately led to the rapid discovery of lead compounds with unique SAR that exhibited potent enzymatic and cellular activity with excellent pan-coronavirus coverage.
Collapse
Affiliation(s)
- Daniel W Carney
- Takeda Development Center Americas, Inc, 9625 Towne Centre Drive, San Diego, CA 92121, United States.
| | - Abba E Leffler
- Schrödinger, Inc, 1540 Broadway, New York, NY 10036, United States.
| | - Jeffrey A Bell
- Schrödinger, Inc, 1540 Broadway, New York, NY 10036, United States
| | | | - Cecilia Cheng
- Schrödinger, Inc, 9868 Scranton Road, Suite 3200, San Diego, CA 92121, United States
| | - Edcon Chang
- Takeda Development Center Americas, Inc, 9625 Towne Centre Drive, San Diego, CA 92121, United States
| | - Adam Dornford
- Schrödinger, Inc, 1 Main St, 11th Floor, Cambridge, MA 02142, United States
| | - Douglas R Dougan
- Takeda Development Center Americas, Inc, 9625 Towne Centre Drive, San Diego, CA 92121, United States
| | - Leah L Frye
- Schrödinger, Inc, 101 SW Main Street, Suite 1300, Portland, OR 97204, United States
| | - Mary E Grimes
- Schrödinger, Inc, 101 SW Main Street, Suite 1300, Portland, OR 97204, United States
| | - Tim Knehans
- Schrödinger GmbH, Glücksteinallee 25, 68163 Mannheim, Germany
| | | | - Mallareddy Komandla
- Takeda Development Center Americas, Inc, 9625 Towne Centre Drive, San Diego, CA 92121, United States
| | - Weston Lane
- Treeline Biosciences, 500 Arsenal Way, Watertown, MA 02472, United States
| | - Hubert Li
- Schrödinger, Inc, 9868 Scranton Road, Suite 3200, San Diego, CA 92121, United States
| | - Sophia R Newman
- Schrödinger, Inc, 101 SW Main Street, Suite 1300, Portland, OR 97204, United States
| | - Katalin Phimister
- Schrödinger Technologies Limited, 1st Floor West, Davidson House, Forbury Square, Reading RG1 3EU, United Kingdom
| | - Kumar S Saikatendu
- Takeda Development Center Americas, Inc, 9625 Towne Centre Drive, San Diego, CA 92121, United States
| | - Hercules Silverstein
- Schrödinger, Inc, 101 SW Main Street, Suite 1300, Portland, OR 97204, United States
| | | |
Collapse
|
3
|
Fedotova MV, Chuev GN. The Three-Dimensional Reference Interaction Site Model Approach as a Promising Tool for Studying Hydrated Viruses and Their Complexes with Ligands. Int J Mol Sci 2024; 25:3697. [PMID: 38612508 PMCID: PMC11011341 DOI: 10.3390/ijms25073697] [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: 03/09/2024] [Revised: 03/22/2024] [Accepted: 03/25/2024] [Indexed: 04/14/2024] Open
Abstract
Viruses are the most numerous biological form living in any ecosystem. Viral diseases affect not only people but also representatives of fauna and flora. The latest pandemic has shown how important it is for the scientific community to respond quickly to the challenge, including critically assessing the viral threat and developing appropriate measures to counter this threat. Scientists around the world are making enormous efforts to solve these problems. In silico methods, which allow quite rapid obtention of, in many cases, accurate information in this field, are effective tools for the description of various aspects of virus activity, including virus-host cell interactions, and, thus, can provide a molecular insight into the mechanism of virus functioning. The three-dimensional reference interaction site model (3D-RISM) seems to be one of the most effective and inexpensive methods to compute hydrated viruses, since the method allows us to provide efficient calculations of hydrated viruses, remaining all molecular details of the liquid environment and virus structure. The pandemic challenge has resulted in a fast increase in the number of 3D-RISM calculations devoted to hydrated viruses. To provide readers with a summary of this literature, we present a systematic overview of the 3D-RISM calculations, covering the period since 2010. We discuss various biophysical aspects of the 3D-RISM results and demonstrate capabilities, limitations, achievements, and prospects of the method using examples of viruses such as influenza, hepatitis, and SARS-CoV-2 viruses.
Collapse
Affiliation(s)
- Marina V. Fedotova
- G.A. Krestov Institute of Solution Chemistry, The Russian Academy of Sciences, Akademicheskaya St., 1, 153045 Ivanovo, Russia
| | - Gennady N. Chuev
- Institute of Theoretical and Experimental Biophysics, The Russian Academy of Sciences, Institutskaya St., 142290 Pushchino, Russia
| |
Collapse
|
4
|
The Development of Pharmacophore Models for the Search of New Natural Inhibitors of SARS-CoV-2 Spike RBD-ACE2 Binding Interface. Molecules 2022; 27:molecules27248938. [PMID: 36558067 PMCID: PMC9788546 DOI: 10.3390/molecules27248938] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/14/2022] [Revised: 12/11/2022] [Accepted: 12/12/2022] [Indexed: 12/23/2022] Open
Abstract
To date, some succeeding variants of SARS-CoV-2 have become more contagious. This virus is known to enter human cells by binding the receptor-binding domain (RBD) of spike protein with the angiotensin-converting enzyme 2 (ACE2), the latter being a membrane protein that regulates the renin-angiotensin system. Since the host cell receptor plays a critical role in viral entry, inhibition of the RBD-ACE2 complex is a promising strategy for preventing COVID-19 infection. In the present communication, we propose and utilize an approach based on the generation of a complex of pharmacophore models and subsequent Induced Fit Docking (IFD) to identify potential inhibitors of the main binding sites of the Omicron SARS-CoV-2 RBD(S1)-ACE2 complex (PDB ID: 7T9L) among a number of natural products of various types and origins. Several natural compounds have been found to provide a high affinity for the receptor of interest. It is expected that the present results will stimulate further research aimed at the development of specialized drugs against this virus.
Collapse
|
5
|
Ebrahim A, Riley BT, Kumaran D, Andi B, Fuchs MR, McSweeney S, Keedy DA. The temperature-dependent conformational ensemble of SARS-CoV-2 main protease (M pro). IUCRJ 2022; 9:682-694. [PMID: 36071812 PMCID: PMC9438506 DOI: 10.1107/s2052252522007497] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/12/2022] [Accepted: 07/21/2022] [Indexed: 05/12/2023]
Abstract
The COVID-19 pandemic, instigated by the SARS-CoV-2 coronavirus, continues to plague the globe. The SARS-CoV-2 main protease, or Mpro, is a promising target for the development of novel antiviral therapeutics. Previous X-ray crystal structures of Mpro were obtained at cryogenic tem-per-ature or room tem-per-ature only. Here we report a series of high-resolution crystal structures of unliganded Mpro across multiple tem-per-atures from cryogenic to physiological, and another at high humidity. We inter-rogate these data sets with parsimonious multiconformer models, multi-copy ensemble models, and isomorphous difference density maps. Our analysis reveals a perturbation-dependent conformational landscape for Mpro, including a mobile zinc ion inter-leaved between the catalytic dyad, mercurial conformational heterogeneity at various sites including a key substrate-binding loop, and a far-reaching intra-molecular network bridging the active site and dimer inter-face. Our results may inspire new strategies for antiviral drug development to aid preparation for future coronavirus pandemics.
Collapse
Affiliation(s)
- Ali Ebrahim
- Diamond Light Source, Harwell Science and Innovation Campus, Didcot, OX11 0DE, England, United Kingdom
- Structural Biology Initiative, CUNY Advanced Science Research Center, New York, NY 10031, USA
| | - Blake T. Riley
- Structural Biology Initiative, CUNY Advanced Science Research Center, New York, NY 10031, USA
| | - Desigan Kumaran
- Biology Department, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - Babak Andi
- Center for BioMolecular Structure, NSLS-II, Brookhaven National Laboratory, Upton, NY 11973, USA
- National Virtual Biotechnology Laboratory (NVBL), US Department of Energy, Washington, DC, USA
| | - Martin R. Fuchs
- Center for BioMolecular Structure, NSLS-II, Brookhaven National Laboratory, Upton, NY 11973, USA
| | - Sean McSweeney
- Center for BioMolecular Structure, NSLS-II, Brookhaven National Laboratory, Upton, NY 11973, USA
- National Virtual Biotechnology Laboratory (NVBL), US Department of Energy, Washington, DC, USA
| | - Daniel A. Keedy
- Structural Biology Initiative, CUNY Advanced Science Research Center, New York, NY 10031, USA
- Department of Chemistry and Biochemistry, City College of New York, New York, NY 10031, USA
- PhD Programs in Biochemistry, Biology, and Chemistry, The Graduate Center–City University of New York, New York, NY 10016, USA
| |
Collapse
|
6
|
Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex. Molecules 2022; 27:molecules27030799. [PMID: 35164065 PMCID: PMC8839245 DOI: 10.3390/molecules27030799] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2021] [Revised: 01/20/2022] [Accepted: 01/21/2022] [Indexed: 12/02/2022] Open
Abstract
The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.
Collapse
|
7
|
Waibl F, Kraml J, Fernández-Quintero ML, Loeffler JR, Liedl KR. Explicit solvation thermodynamics in ionic solution: extending grid inhomogeneous solvation theory to solvation free energy of salt-water mixtures. J Comput Aided Mol Des 2022; 36:101-116. [PMID: 35031880 PMCID: PMC8907097 DOI: 10.1007/s10822-021-00429-y] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Accepted: 10/28/2021] [Indexed: 12/03/2022]
Abstract
Hydration thermodynamics play a fundamental role in fields ranging from the pharmaceutical industry to environmental research. Numerous methods exist to predict solvation thermodynamics of compounds ranging from small molecules to large biomolecules. Arguably the most precise methods are those based on molecular dynamics (MD) simulations in explicit solvent. One theory that has seen increased use is inhomogeneous solvation theory (IST). However, while many applications require accurate description of salt-water mixtures, no implementation of IST is currently able to estimate solvation properties involving more than one solvent species. Here, we present an extension to grid inhomogeneous solvation theory (GIST) that can take salt contributions into account. At the example of carbazole in 1 M NaCl solution, we compute the solvation energy as well as first and second order entropies. While the effect of the first order ion entropy is small, both the water-water and water-ion entropies contribute strongly. We show that the water-ion entropies are efficiently approximated using the Kirkwood superposition approximation. However, this approach cannot be applied to the water-water entropy. Furthermore, we test the quantitative validity of our method by computing salting-out coefficients and comparing them to experimental data. We find a good correlation to experimental salting-out constants, while the absolute values are overpredicted due to the approximate second order entropy. Since ions are frequently used in MD, either to neutralize the system or as a part of the investigated process, our method greatly extends the applicability of GIST. The use-cases range from biopharmaceuticals, where many assays require high salt concentrations, to environmental research, where solubility in sea water is important to model the fate of organic substances.
Collapse
Affiliation(s)
- Franz Waibl
- Department of General, Inorganic, and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, 6020, Innsbruck, Austria
| | - Johannes Kraml
- Department of General, Inorganic, and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, 6020, Innsbruck, Austria
| | - Monica L Fernández-Quintero
- Department of General, Inorganic, and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, 6020, Innsbruck, Austria
| | - Johannes R Loeffler
- Department of General, Inorganic, and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, 6020, Innsbruck, Austria
| | - Klaus R Liedl
- Department of General, Inorganic, and Theoretical Chemistry, University of Innsbruck, Innrain 80/82, 6020, Innsbruck, Austria.
| |
Collapse
|
8
|
Ebrahim A, Riley BT, Kumaran D, Andi B, Fuchs MR, McSweeney S, Keedy DA. The temperature-dependent conformational ensemble of SARS-CoV-2 main protease (M pro). BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2021:2021.05.03.437411. [PMID: 33972941 PMCID: PMC8109201 DOI: 10.1101/2021.05.03.437411] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
Abstract
The COVID-19 pandemic, instigated by the SARS-CoV-2 coronavirus, continues to plague the globe. The SARS-CoV-2 main protease, or Mpro, is a promising target for development of novel antiviral therapeutics. Previous X-ray crystal structures of Mpro were obtained at cryogenic temperature or room temperature only. Here we report a series of high-resolution crystal structures of unliganded Mpro across multiple temperatures from cryogenic to physiological, and another at high humidity. We interrogate these datasets with parsimonious multiconformer models, multi-copy ensemble models, and isomorphous difference density maps. Our analysis reveals a temperature-dependent conformational landscape for Mpro, including mobile solvent interleaved between the catalytic dyad, mercurial conformational heterogeneity in a key substrate-binding loop, and a far-reaching intramolecular network bridging the active site and dimer interface. Our results may inspire new strategies for antiviral drug development to counter-punch COVID-19 and combat future coronavirus pandemics.
Collapse
Affiliation(s)
- Ali Ebrahim
- Diamond Light Source, Harwell Science and Innovation Campus, Didcot OX11 0DE, England
- Structural Biology Initiative, CUNY Advanced Science Research Center, New York, NY 10031
| | - Blake T. Riley
- Structural Biology Initiative, CUNY Advanced Science Research Center, New York, NY 10031
| | - Desigan Kumaran
- Biology Department, Brookhaven National Laboratory, Upton, NY 11973
| | - Babak Andi
- National Synchrotron Light Source II, Brookhaven National Laboratory, Upton, NY 11973
- National Virtual Biotechnology Laboratory (NVBL), US Department of Energy, Washington, DC, United States
| | - Martin R. Fuchs
- National Synchrotron Light Source II, Brookhaven National Laboratory, Upton, NY 11973
| | - Sean McSweeney
- National Synchrotron Light Source II, Brookhaven National Laboratory, Upton, NY 11973
- National Virtual Biotechnology Laboratory (NVBL), US Department of Energy, Washington, DC, United States
| | - Daniel A. Keedy
- Structural Biology Initiative, CUNY Advanced Science Research Center, New York, NY 10031
- Department of Chemistry and Biochemistry, City College of New York, New York, NY 10031
- Ph.D. Programs in Biochemistry, Biology, and Chemistry, The Graduate Center – City University of New York, New York, NY 10016
| |
Collapse
|
9
|
Kobryn AE, Maruyama Y, Velázquez-Martínez CA, Yoshida N, Gusarov S. Modeling the interaction of SARS-CoV-2 binding to the ACE2 receptor via molecular theory of solvation. NEW J CHEM 2021. [DOI: 10.1039/d1nj02015c] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The angiotensin-converting enzyme 2 (ACE2) protein is a cell gate receptor for the SARS-CoV-2 virus, responsible for the development of symptoms associated with the Covid-19 disease.
Collapse
Affiliation(s)
- Alexander E. Kobryn
- Nanotechnology Research Centre, National Research Council Canada, 11421 Saskatchewan Drive NW, Edmonton, Alberta, T6G 2M9, Canada
| | - Yutaka Maruyama
- Architecture Development Team, FLAGSHIP 2020 Project, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - Carlos A. Velázquez-Martínez
- 2142-L Katz Group Centre for Research, University of Alberta, 11315-87 Avenue NW, Edmonton, Alberta, T6G 2H5, Canada
| | - Norio Yoshida
- Department of Chemistry, Graduate School of Science, Kyushu University, 744, Motooka, Nishi-ku, Fukuoka, 819-0395, Japan
| | - Sergey Gusarov
- Nanotechnology Research Centre, National Research Council Canada, 11421 Saskatchewan Drive NW, Edmonton, Alberta, T6G 2M9, Canada
| |
Collapse
|