1
|
Kariyawasam NL, Ploetz EA, Swint-Kruse L, Smith PE. Simulated pressure changes in LacI suggest a link between hydration and functional conformational changes. Biophys Chem 2024; 304:107126. [PMID: 37924711 PMCID: PMC10842697 DOI: 10.1016/j.bpc.2023.107126] [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: 08/28/2023] [Revised: 10/18/2023] [Accepted: 10/25/2023] [Indexed: 11/06/2023]
Abstract
The functions of many proteins are associated with interconversions among conformational substates. However, these substates can be difficult to measure experimentally, and determining contributions from hydration changes can be especially difficult. Here, we assessed the use of pressure perturbations to sample the substates accessible to the Escherichia coli lactose repressor protein (LacI) in various liganded forms. In the presence of DNA, the regulatory domain of LacI adopts an Open conformation that, in the absence of DNA, changes to a Closed conformation. Increasing the simulation pressure prevented the transition from an Open to a Closed conformation, in a similar manner to the binding of DNA and anti-inducer, ONPF. The results suggest the hydration of specific residues play a significant role in determining the population of different LacI substates and that simulating pressure perturbation could be useful for assessing the role of hydration changes that accompany functionally-relevant amino acid substitutions.
Collapse
Affiliation(s)
- Nilusha L Kariyawasam
- Department of Chemistry, 213 CBC Building, 1212 Mid-Campus Dr. North, Kansas State University, Manhattan, KS 66506, USA
| | - Elizabeth A Ploetz
- Department of Chemistry, 213 CBC Building, 1212 Mid-Campus Dr. North, Kansas State University, Manhattan, KS 66506, USA
| | - Liskin Swint-Kruse
- Department of Biochemistry and Molecular Biology, MSN 3030, 3901 Rainbow Blvd, The University of Kansas Medical Center, Kansas City, KS 66160, USA
| | - Paul E Smith
- Department of Chemistry, 213 CBC Building, 1212 Mid-Campus Dr. North, Kansas State University, Manhattan, KS 66506, USA.
| |
Collapse
|
2
|
Gracia Carmona O, Gillhofer M, Tomasiak L, De Ruiter A, Oostenbrink C. Accelerated Enveloping Distribution Sampling to Probe the Presence of Water Molecules. J Chem Theory Comput 2023. [PMID: 37167545 DOI: 10.1021/acs.jctc.3c00109] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Determining the presence of water molecules at protein-ligand interfaces is still a challenging task in free-energy calculations. The inappropriate placement of water molecules results in the stabilization of wrong conformational orientations of the ligand. With classical alchemical perturbation methods, such as thermodynamic integration (TI), it is essential to know the amount of water molecules in the active site of the respective ligands. However, the resolution of the crystal structure and the correct assignment of the electron density do not always lead to a clear placement of water molecules. In this work, we apply the one-step perturbation method named accelerated enveloping distribution sampling (AEDS) to determine the presence of water molecules in the active site by probing them in a fast and straightforward way. Based on these results, we combined the AEDS method with standard TI to calculate accurate binding free energies in the presence of buried water molecules. The main idea is to perturb the water molecules with AEDS such that they are allowed to alternate between regular water molecules and non-interacting dummy particles while treating the ligand with TI over an alchemical pathway. We demonstrate the use of AEDS to probe the presence of water molecules for six different test systems. For one of these, previous calculations showed difficulties to reproduce the experimental binding free energies, and here, we use the combined TI-AEDS approach to tackle these issues.
Collapse
Affiliation(s)
- Oriol Gracia Carmona
- Institute for Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Michael Gillhofer
- Institute for Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Lisa Tomasiak
- Institute for Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Anita De Ruiter
- Institute for Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190 Vienna, Austria
| | - Chris Oostenbrink
- Institute for Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190 Vienna, Austria
- Christian Doppler Laboratory for Molecular Informatics in the Biosciences, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190 Vienna, Austria
| |
Collapse
|
3
|
Omori A, Sasaki S, Kikukawa T, Shimono K, Miyauchi S. Elucidation of a Thermodynamical Feature Attributed to Substrate Binding to the Prokaryotic H +/Oligopeptide Cotransporter YdgR with Calorimetric Analysis: The Substrate Binding Driven by the Change in Entropy Implies the Release of Bound Water Molecules from the Binding Pocket. Biochemistry 2023. [PMID: 37163674 DOI: 10.1021/acs.biochem.2c00673] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/12/2023]
Abstract
Here, we have elucidated the substrate recognition mechanism by a prokaryotic H+/oligopeptide cotransporter, YdgR, using isothermal titration calorimetry. Under acidic conditions (pH 6.0), the binding of a dipeptide, Val-Ala, to YdgR elicited endothermic enthalpy, which compensated for the increase in entropy due to dipeptide binding. A series of dipeptides were used in the binding titration. The dipeptides represent Val-X and X-Val, where X is Ala, Ser, Val, Tyr, or Phe. Most dipeptides revealed endothermic enthalpy, which was completely compensated by the increase in entropy due to dipeptide binding. The change in enthalpy due to binding correlated well with the change in entropy, whereas the Gibbs free energy involved in the binding of the dipeptide to YdgR remained unchanged irrespective of dipeptide sequences, implying that the binding reaction was driven by entropy, that is, the release of bound water molecules in the binding pocket. It is also important to clarify that, based on the prediction of water molecules in the ligand-binding pocket of YdgR, the release of three bound water molecules in the putative substrate binding pocket occurred through binding to YdgR. In the comparison of Val-X and X-Val dipeptides, the N-terminal region of the binding pocket might contain more bound water molecules than the C-terminal region. In light of these findings, we suggest that bound water molecules might play an important role in substrate recognition and binding by YdgR.
Collapse
Affiliation(s)
- Akiko Omori
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
| | - Shotaro Sasaki
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
| | - Takashi Kikukawa
- Faculty of Advanced Life Science, Hokkaido University, Kita 10, Nishi 8, Kita-ku, Sapporo 060-0810, Japan
| | - Kazumi Shimono
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
- Faculty of Pharmaceutical Sciences, Sojo University, 4-22-1 Ikeda Nishi-ku, Kumamoto 860-0082, Japan
| | - Seiji Miyauchi
- Faculty of Pharmaceutical Sciences, Toho University, 2-2-1 Miyama, Funabashi, Chiba 274-8510, Japan
| |
Collapse
|
4
|
Eberhardt J, Forli S. WaterKit: Thermodynamic Profiling of Protein Hydration Sites. J Chem Theory Comput 2023; 19:2535-2556. [PMID: 37094087 PMCID: PMC10732097 DOI: 10.1021/acs.jctc.2c01087] [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] [Indexed: 04/26/2023]
Abstract
Water desolvation is one of the key components of the free energy of binding of small molecules to their receptors. Thus, understanding the energetic balance of solvation and desolvation resulting from individual water molecules can be crucial when estimating ligand binding, especially when evaluating different molecules and poses as done in High-Throughput Virtual Screening (HTVS). Over the most recent decades, several methods were developed to tackle this problem, ranging from fast approximate methods (usually empirical functions using either discrete atom-atom pairwise interactions or continuum solvent models) to more computationally expensive and accurate ones, mostly based on Molecular Dynamics (MD) simulations, such as Grid Inhomogeneous Solvation Theory (GIST) or Double Decoupling. On one hand, MD-based methods are prohibitive to use in HTVS to estimate the role of waters on the fly for each ligand. On the other hand, fast and approximate methods show an unsatisfactory level of accuracy, with low agreement with results obtained with the more expensive methods. Here we introduce WaterKit, a new grid-based sampling method with explicit water molecules to calculate thermodynamic properties using the GIST method. Our results show that the discrete placement of water molecules is successful in reproducing the position of crystallographic waters with very high accuracy, as well as providing thermodynamic estimates with accuracy comparable to more expensive MD simulations. Unlike these methods, WaterKit can be used to analyze specific regions on the protein surface, (such as the binding site of a receptor), without having to hydrate and simulate the whole receptor structure. The results show the feasibility of a general and fast method to compute thermodynamic properties of water molecules, making it well-suited to be integrated in high-throughput pipelines such as molecular docking.
Collapse
Affiliation(s)
- Jerome Eberhardt
- Department of Integrative Structural and Computational Biology, Scripps Research, La Jolla, California 92037, United States
| | - Stefano Forli
- Department of Integrative Structural and Computational Biology, Scripps Research, La Jolla, California 92037, United States
| |
Collapse
|
5
|
Ge Y, Melling OJ, Dong W, Essex JW, Mobley DL. Enhancing sampling of water rehydration upon ligand binding using variants of grand canonical Monte Carlo. J Comput Aided Mol Des 2022; 36:767-779. [PMID: 36198874 PMCID: PMC9869699 DOI: 10.1007/s10822-022-00479-w] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2022] [Accepted: 09/15/2022] [Indexed: 01/26/2023]
Abstract
Water plays an important role in mediating protein-ligand interactions. Water rearrangement upon a ligand binding or modification can be very slow and beyond typical timescales used in molecular dynamics (MD) simulations. Thus, inadequate sampling of slow water motions in MD simulations often impairs the accuracy of the accuracy of ligand binding free energy calculations. Previous studies suggest grand canonical Monte Carlo (GCMC) outperforms normal MD simulations for water sampling, thus GCMC has been applied to help improve the accuracy of ligand binding free energy calculations. However, in prior work we observed protein and/or ligand motions impaired how well GCMC performs at water rehydration, suggesting more work is needed to improve this method to handle water sampling. In this work, we applied GCMC in 21 protein-ligand systems to assess the performance of GCMC for rehydrating buried water sites. While our results show that GCMC can rapidly rehydrate all selected water sites for most systems, it fails in five systems. In most failed systems, we observe protein/ligand motions, which occur in the absence of water, combine to close water sites and block instantaneous GCMC water insertion moves. For these five failed systems, we both extended our GCMC simulations and tested a new technique named grand canonical nonequilibrium candidate Monte Carlo (GCNCMC). GCNCMC combines GCMC with the nonequilibrium candidate Monte Carlo (NCMC) sampling technique to improve the probability of a successful water insertion/deletion. Our results show that GCNCMC and extended GCMC can rehydrate all target water sites for three of the five problematic systems and GCNCMC is more efficient than GCMC in two out of the three systems. In one system, only GCNCMC can rehydrate all target water sites, while GCMC fails. Both GCNCMC and GCMC fail in one system. This work suggests this new GCNCMC method is promising for water rehydration especially when protein/ligand motions may block water insertion/removal.
Collapse
Affiliation(s)
- Yunhui Ge
- Department of Pharmaceutical Sciences, University of California, Irvine, Irvine, CA, 92697, USA
| | - Oliver J Melling
- School of Chemistry, University of Southampton, Southampton, SO17 1BJ, UK
| | - Weiming Dong
- Department of Pharmaceutical Sciences, University of California, Irvine, Irvine, CA, 92697, USA
| | - Jonathan W Essex
- School of Chemistry, University of Southampton, Southampton, SO17 1BJ, UK
| | - David L Mobley
- Department of Pharmaceutical Sciences, University of California, Irvine, Irvine, CA, 92697, USA.
- Department of Chemistry, University of California,Irvine, Irvine, CA, 92697, USA.
| |
Collapse
|
6
|
Chai HH, Ham JS, Kim TH, Lim D. Identifying ligand-binding specificity of the oligopeptide receptor OppA from Bifidobacterium longum KACC91563 by structure-based molecular modeling. ARAB J CHEM 2022. [DOI: 10.1016/j.arabjc.2022.104198] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
|
7
|
Ge Y, Wych DC, Samways ML, Wall ME, Essex JW, Mobley DL. Enhancing Sampling of Water Rehydration on Ligand Binding: A Comparison of Techniques. J Chem Theory Comput 2022; 18:1359-1381. [PMID: 35148093 PMCID: PMC9241631 DOI: 10.1021/acs.jctc.1c00590] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Water often plays a key role in protein structure, molecular recognition, and mediating protein-ligand interactions. Thus, free energy calculations must adequately sample water motions, which often proves challenging in typical MD simulation time scales. Thus, the accuracy of methods relying on MD simulations ends up limited by slow water sampling. Particularly, as a ligand is removed or modified, bulk water may not have time to fill or rearrange in the binding site. In this work, we focus on several molecular dynamics (MD) simulation-based methods attempting to help rehydrate buried water sites: BLUES, using nonequilibrium candidate Monte Carlo (NCMC); grand, using grand canonical Monte Carlo (GCMC); and normal MD. We assess the accuracy and efficiency of these methods in rehydrating target water sites. We selected a range of systems with varying numbers of waters in the binding site, as well as those where water occupancy is coupled to the identity or binding mode of the ligand. We analyzed the rehydration of buried water sites in binding pockets using both clustering of trajectories and direct analysis of electron density maps. Our results suggest both BLUES and grand enhance water sampling relative to normal MD and grand is more robust than BLUES, but also that water sampling remains a major challenge for all of the methods tested. The lessons we learned for these methods and systems are discussed.
Collapse
Affiliation(s)
- Yunhui Ge
- Department of Pharmaceutical Sciences, University of California, Irvine, California 92697, United States
| | - David C Wych
- Department of Pharmaceutical Sciences, University of California, Irvine, California 92697, United States
- Computer, Computational, and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Marley L Samways
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, United Kingdom
| | - Michael E Wall
- Computer, Computational, and Statistical Sciences Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Jonathan W Essex
- School of Chemistry, University of Southampton, Southampton SO17 1BJ, United Kingdom
| | - David L Mobley
- Department of Pharmaceutical Sciences, University of California, Irvine, California 92697, United States
- Department of Chemistry, University of California, Irvine, California 92697, United States
| |
Collapse
|
8
|
Öhlknecht C, Katz S, Kröß C, Sprenger B, Engele P, Schneider R, Oostenbrink C. Efficient In Silico Saturation Mutagenesis of a Member of the Caspase Protease Family. J Chem Inf Model 2021; 61:1193-1203. [PMID: 33570387 PMCID: PMC8023567 DOI: 10.1021/acs.jcim.0c01216] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2020] [Indexed: 12/28/2022]
Abstract
Rational-design methods have proven to be a valuable toolkit in the field of protein design. Numerical approaches such as free-energy calculations or QM/MM methods are fit to widen the understanding of a protein-sequence space but require large amounts of computational time and power. Here, we apply an efficient method for free-energy calculations that combines the one-step perturbation (OSP) with the third-power-fitting (TPF) approach. It is fit to calculate full free energies of binding from three different end states only. The nonpolar contribution to the free energies are calculated for a set of chosen amino acids from a single simulation of a judiciously chosen reference state. The electrostatic contributions, on the other hand, are predicted from simulations of the neutral and charged end states of the individual amino acids. We used this method to perform in silico saturation mutagenesis of two sites in human Caspase-2. We calculated relative binding free energies toward two different substrates that differ in their P1' site and in their affinity toward the unmutated protease. Although being approximate, our approach showed very good agreement upon validation against experimental data. 76% of the predicted relative free energies of amino acid mutations was found to be true positives or true negatives. We observed that this method is fit to discriminate amino acid mutations because the rate of false negatives is very low (<1.5%). The approach works better for a substrate with medium/low affinity with a Matthews correlation coefficient (MCC) of 0.63, whereas for a substrate with very low affinity, the MCC was 0.38. In all cases, the combined TPF + OSP approach outperformed the linear interaction energy method.
Collapse
Affiliation(s)
- Christoph Öhlknecht
- Institute
of Molecular Modeling and Simulation, University
of Natural Resources and Life Sciences, Vienna A-1190, Austria
- Austrian
Centre of Industrial Biotechnology, Petersgasse 14, Graz 8041, Austria
| | - Sonja Katz
- Institute
of Molecular Modeling and Simulation, University
of Natural Resources and Life Sciences, Vienna A-1190, Austria
- Austrian
Centre of Industrial Biotechnology, Petersgasse 14, Graz 8041, Austria
| | - Christina Kröß
- Austrian
Centre of Industrial Biotechnology, Petersgasse 14, Graz 8041, Austria
- Institute
of Biochemistry and Center of Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck 6020, Austria
| | - Bernhard Sprenger
- Austrian
Centre of Industrial Biotechnology, Petersgasse 14, Graz 8041, Austria
- Institute
of Biochemistry and Center of Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck 6020, Austria
| | - Petra Engele
- Austrian
Centre of Industrial Biotechnology, Petersgasse 14, Graz 8041, Austria
- Institute
of Biochemistry and Center of Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck 6020, Austria
| | - Rainer Schneider
- Institute
of Biochemistry and Center of Molecular Biosciences Innsbruck, University of Innsbruck, Innsbruck 6020, Austria
| | - Chris Oostenbrink
- Institute
of Molecular Modeling and Simulation, University
of Natural Resources and Life Sciences, Vienna A-1190, Austria
| |
Collapse
|
9
|
Bergazin TD, Ben-Shalom IY, Lim NM, Gill SC, Gilson MK, Mobley DL. Enhancing water sampling of buried binding sites using nonequilibrium candidate Monte Carlo. J Comput Aided Mol Des 2021; 35:167-177. [PMID: 32968887 PMCID: PMC7904576 DOI: 10.1007/s10822-020-00344-8] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Accepted: 09/16/2020] [Indexed: 11/26/2022]
Abstract
Water molecules can be found interacting with the surface and within cavities in proteins. However, water exchange between bulk and buried hydration sites can be slow compared to simulation timescales, thus leading to the inefficient sampling of the locations of water. This can pose problems for free energy calculations for computer-aided drug design. Here, we apply a hybrid method that combines nonequilibrium candidate Monte Carlo (NCMC) simulations and molecular dynamics (MD) to enhance sampling of water in specific areas of a system, such as the binding site of a protein. Our approach uses NCMC to gradually remove interactions between a selected water molecule and its environment, then translates the water to a new region, before turning the interactions back on. This approach of gradual removal of interactions, followed by a move and then reintroduction of interactions, allows the environment to relax in response to the proposed water translation, improving acceptance of moves and thereby accelerating water exchange and sampling. We validate this approach on several test systems including the ligand-bound MUP-1 and HSP90 proteins with buried crystallographic waters removed. We show that our BLUES (NCMC/MD) method enhances water sampling relative to normal MD when applied to these systems. Thus, this approach provides a strategy to improve water sampling in molecular simulations which may be useful in practical applications in drug discovery and biomolecular design.
Collapse
Affiliation(s)
| | - Ido Y Ben-Shalom
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, La Jolla, CA, 92093, USA
| | - Nathan M Lim
- Department of Pharmaceutical Sciences, University of California, Irvine, CA, 92697, USA
| | - Sam C Gill
- Department of Chemistry, University of California, Irvine, Irvine, CA, 92697, USA
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, La Jolla, CA, 92093, USA
| | - David L Mobley
- Department of Pharmaceutical Sciences, University of California, Irvine, CA, 92697, USA.
- Department of Chemistry, University of California, Irvine, Irvine, CA, 92697, USA.
| |
Collapse
|
10
|
Ben-Shalom IY, Lin Z, Radak BK, Lin C, Sherman W, Gilson MK. Accounting for the Central Role of Interfacial Water in Protein-Ligand Binding Free Energy Calculations. J Chem Theory Comput 2020; 16:7883-7894. [PMID: 33206520 PMCID: PMC7725968 DOI: 10.1021/acs.jctc.0c00785] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Rigorous binding free energy methods in drug discovery are growing in popularity because of a combination of methodological advances, improvements in computer hardware, and workflow automation. These calculations typically use molecular dynamics (MD) to sample from the Boltzmann distribution of conformational states. However, when part or all of the binding sites is inaccessible to the bulk solvent, the time needed for water molecules to equilibrate between bulk solvent and the binding site can be well beyond what is practical with standard MD. This sampling limitation is problematic in relative binding free energy calculations, which compute the reversible work of converting ligand 1 to ligand 2 within the binding site. Thus, if ligand 1 is smaller and/or more polar than ligand 2, the perturbation may allow additional water molecules to occupy a region of the binding site. However, this change in hydration may not be captured by standard MD simulations and may therefore lead to errors in the computed free energy. We recently developed a hybrid Monte Carlo/MD (MC/MD) method, which speeds up the equilibration of water between bulk solvent and buried cavities, while sampling from the intended distribution of states. Here, we report on the use of this approach in the context of alchemical binding free energy calculations. We find that using MC/MD markedly improves the accuracy of the calculations and also reduces hysteresis between the forward and reverse perturbations, relative to matched calculations using only MD with or without the crystallographic water molecules. The present method is available for use in AMBER simulation software.
Collapse
Affiliation(s)
- Ido Y Ben-Shalom
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, 92093 La Jolla, California, United States
| | - Zhixiong Lin
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Brian K Radak
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Charles Lin
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Woody Sherman
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, 92093 La Jolla, California, United States
| |
Collapse
|
11
|
Ben-Shalom IY, Lin Z, Radak BK, Lin C, Sherman W, Gilson MK. Accounting for the Central Role of Interfacial Water in Protein-Ligand Binding Free Energy Calculations. J Chem Theory Comput 2020. [PMID: 33206520 DOI: 10.26434/chemrxiv.12668816.v1] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2023]
Abstract
Rigorous binding free energy methods in drug discovery are growing in popularity because of a combination of methodological advances, improvements in computer hardware, and workflow automation. These calculations typically use molecular dynamics (MD) to sample from the Boltzmann distribution of conformational states. However, when part or all of the binding sites is inaccessible to the bulk solvent, the time needed for water molecules to equilibrate between bulk solvent and the binding site can be well beyond what is practical with standard MD. This sampling limitation is problematic in relative binding free energy calculations, which compute the reversible work of converting ligand 1 to ligand 2 within the binding site. Thus, if ligand 1 is smaller and/or more polar than ligand 2, the perturbation may allow additional water molecules to occupy a region of the binding site. However, this change in hydration may not be captured by standard MD simulations and may therefore lead to errors in the computed free energy. We recently developed a hybrid Monte Carlo/MD (MC/MD) method, which speeds up the equilibration of water between bulk solvent and buried cavities, while sampling from the intended distribution of states. Here, we report on the use of this approach in the context of alchemical binding free energy calculations. We find that using MC/MD markedly improves the accuracy of the calculations and also reduces hysteresis between the forward and reverse perturbations, relative to matched calculations using only MD with or without the crystallographic water molecules. The present method is available for use in AMBER simulation software.
Collapse
Affiliation(s)
- Ido Y Ben-Shalom
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, 92093 La Jolla, California, United States
| | - Zhixiong Lin
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Brian K Radak
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Charles Lin
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Woody Sherman
- Silicon Therapeutics LLC, Boston, Massachusetts 02110, United States
| | - Michael K Gilson
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California, San Diego, 92093 La Jolla, California, United States
| |
Collapse
|
12
|
Hahn DF, König G, Hünenberger PH. Overcoming Orthogonal Barriers in Alchemical Free Energy Calculations: On the Relative Merits of λ-Variations, λ-Extrapolations, and Biasing. J Chem Theory Comput 2020; 16:1630-1645. [DOI: 10.1021/acs.jctc.9b00853] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- David F. Hahn
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Gerhard König
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Philippe H. Hünenberger
- Laboratory of Physical Chemistry, Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| |
Collapse
|
13
|
Zou J, Simmerling C, Raleigh DP. Dissecting the Energetics of Intrinsically Disordered Proteins via a Hybrid Experimental and Computational Approach. J Phys Chem B 2019; 123:10394-10402. [PMID: 31702919 PMCID: PMC7291390 DOI: 10.1021/acs.jpcb.9b08323] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Intrinsically disordered proteins (IDPs) play important roles in biology, but little is known about the energetics of their inter-residue interactions. Methods that have been successfully applied to analyze the energetics of globular proteins are not applicable to the fluctuating partially ordered ensembles populated by IDPs. A combined computational experimental strategy is introduced for analyzing the energetic role of individual residues in the free state of IDPs. The approach combines experimental measurements of the binding of wild-type and mutant IDPs to their partners with alchemical free energy calculations of the structured complexes. These data allow quantitative information to be deduced about the free state via a thermodynamic cycle. The approach is validated by the analysis of the effects of mutations upon the binding free energy of the ovomucoid inhibitor third binding domain to its partners and is applied to the C-terminal domain of the measles virus nucleoprotein, a 125-residue IDP involved in the RNA transcription and replication of measles virus. The analysis reveals significant inter-residue interactions in the unbound IDP and suggests a biological role for them. The work demonstrates that advances in force fields and computational hardware have now led to the point where it is possible to develop methods, which integrate experimental and computational techniques to reveal insights that cannot be studied using either technique alone.
Collapse
Affiliation(s)
- Junjie Zou
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794-3400, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-3400, United S tates
| | - Carlos Simmerling
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794-3400, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-3400, United S tates
| | - Daniel P. Raleigh
- Department of Chemistry, Stony Brook University, Stony Brook, New York 11794-3400, United States
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, New York 11794-3400, United S tates
| |
Collapse
|
14
|
Maurer M, Oostenbrink C. Water in protein hydration and ligand recognition. J Mol Recognit 2019; 32:e2810. [PMID: 31456282 PMCID: PMC6899928 DOI: 10.1002/jmr.2810] [Citation(s) in RCA: 72] [Impact Index Per Article: 14.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2019] [Revised: 07/31/2019] [Accepted: 08/01/2019] [Indexed: 12/16/2022]
Abstract
This review describes selected basics of water in biomolecular recognition. We focus on a qualitative understanding of the most important physical aspects, how these change in magnitude between bulk water and protein environment, and how the roles that water plays for proteins arise from them. These roles include mechanical support, thermal coupling, dielectric screening, mass and charge transport, and the competition with a ligand for the occupation of a binding site. The presence or absence of water has ramifications that range from the thermodynamic binding signature of a single ligand up to cellular survival. The large inhomogeneity in water density, polarity and mobility around a solute is hard to assess in experiment. This is a source of many difficulties in the solvation of protein models and computational studies that attempt to elucidate or predict ligand recognition. The influence of water in a protein binding site on the experimental enthalpic and entropic signature of ligand binding is still a point of much debate. The strong water‐water interaction in enthalpic terms is counteracted by a water molecule's high mobility in entropic terms. The complete arrest of a water molecule's mobility sets a limit on the entropic contribution of a water displacement process, while the solvent environment sets limits on ligand reactivity.
Collapse
Affiliation(s)
- Manuela Maurer
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| | - Chris Oostenbrink
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Austria
| |
Collapse
|
15
|
Wahl J, Smieško M. Assessing the Predictive Power of Relative Binding Free Energy Calculations for Test Cases Involving Displacement of Binding Site Water Molecules. J Chem Inf Model 2019; 59:754-765. [DOI: 10.1021/acs.jcim.8b00826] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Joel Wahl
- Molecular Modeling, Department of Pharmaceutical Sciences, University of Basel, Klingelbergstrasse 50, CH-4056 Basel, Switzerland
| | - Martin Smieško
- Molecular Modeling, Department of Pharmaceutical Sciences, University of Basel, Klingelbergstrasse 50, CH-4056 Basel, Switzerland
| |
Collapse
|
16
|
Insights into the Molecular Mechanisms of Protein-Ligand Interactions by Molecular Docking and Molecular Dynamics Simulation: A Case of Oligopeptide Binding Protein. COMPUTATIONAL AND MATHEMATICAL METHODS IN MEDICINE 2018; 2018:3502514. [PMID: 30627209 PMCID: PMC6305025 DOI: 10.1155/2018/3502514] [Citation(s) in RCA: 79] [Impact Index Per Article: 13.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/29/2018] [Accepted: 10/16/2018] [Indexed: 01/14/2023]
Abstract
Protein-ligand interactions are a necessary prerequisite for signal transduction, immunoreaction, and gene regulation. Protein-ligand interaction studies are important for understanding the mechanisms of biological regulation, and they provide a theoretical basis for the design and discovery of new drug targets. In this study, we analyzed the molecular interactions of protein-ligand which was docked by AutoDock 4.2 software. In AutoDock 4.2 software, we used a new search algorithm, hybrid algorithm of random drift particle swarm optimization and local search (LRDPSO), and the classical Lamarckian genetic algorithm (LGA) as energy optimization algorithms. The best conformations of each docking algorithm were subjected to molecular dynamic (MD) simulations to further analyze the molecular mechanisms of protein-ligand interactions. Here, we analyze the binding energy between protein receptors and ligands, the interactions of salt bridges and hydrogen bonds in the docking region, and the structural changes during complex unfolding. Our comparison of these complexes highlights differences in the protein-ligand interactions between the two docking methods. It also shows that salt bridge and hydrogen bond interactions play a crucial role in protein-ligand stability. The present work focuses on extracting the deterministic characteristics of docking interactions from their dynamic properties, which is important for understanding biological functions and determining which amino acid residues are crucial to docking interactions.
Collapse
|
17
|
Maurer M, Hansen N, Oostenbrink C. Comparison of free-energy methods using a tripeptide-water model system. J Comput Chem 2018; 39:2226-2242. [PMID: 30280398 PMCID: PMC6220940 DOI: 10.1002/jcc.25537] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2018] [Revised: 06/29/2018] [Accepted: 06/30/2018] [Indexed: 11/24/2022]
Abstract
We investigate the ability of several free-energy calculation methods to combine two alchemical changes. We use Bennett acceptance ratio (BAR), thermodynamic integration (TI), extended TI (X-TI), and enveloping distribution sampling (EDS) to perturb a water molecule, which is restrained to an amino acid that is also being perturbed. In addition to these pairwise methods, we present two two-dimensional approaches, EDS-TI and two-dimensional TI (2D-TI). We compare feasibility, efficiency and usability of these methods in regard to our simple model system, which mimics the displacement of a water molecule in the active site of a protein on residue mutation. The correct treatment of structural water has been shown to greatly aid binding affinity calculations in some cases that remained elusive otherwise. This is of broad interest in, for example, drug design, and we conclude that thus far, the pairwise method BAR and also the newer X-TI remain the most suitable methods to treat this problem as long as few end states are involved. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Manuela Maurer
- Department of Material Sciences and Process EngineeringInstitute of Molecular Modeling and Simulation, University of Natural Resources and Life SciencesMuthgasse 18, A‐1190, ViennaAustria
| | - Niels Hansen
- Department of Energy‐, Process‐ and Bio‐Engineering, University of StuttgartInstitute of Thermodynamics and Thermal Process EngineeringPfaffenwaldring 9, 70569, StuttgartGermany
| | - Chris Oostenbrink
- Department of Material Sciences and Process EngineeringInstitute of Molecular Modeling and Simulation, University of Natural Resources and Life SciencesMuthgasse 18, A‐1190, ViennaAustria
| |
Collapse
|
18
|
Wahl J, Smieško M. Thermodynamic Insight into the Effects of Water Displacement and Rearrangement upon Ligand Modifications using Molecular Dynamics Simulations. ChemMedChem 2018; 13:1325-1335. [DOI: 10.1002/cmdc.201800093] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2018] [Revised: 04/07/2018] [Indexed: 01/11/2023]
Affiliation(s)
- Joel Wahl
- Molecular Modeling, Department of Pharmaceutical Sciences; University of Basel; Klingelbergstrasse 50 4056 Basel Switzerland
| | - Martin Smieško
- Molecular Modeling, Department of Pharmaceutical Sciences; University of Basel; Klingelbergstrasse 50 4056 Basel Switzerland
| |
Collapse
|
19
|
Rudling A, Orro A, Carlsson J. Prediction of Ordered Water Molecules in Protein Binding Sites from Molecular Dynamics Simulations: The Impact of Ligand Binding on Hydration Networks. J Chem Inf Model 2018; 58:350-361. [PMID: 29308882 PMCID: PMC6716772 DOI: 10.1021/acs.jcim.7b00520] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
![]()
Water
plays a major role in ligand binding and is attracting increasing
attention in structure-based drug design. Water molecules can make
large contributions to binding affinity by bridging protein–ligand
interactions or by being displaced upon complex formation, but these
phenomena are challenging to model at the molecular level. Herein,
networks of ordered water molecules in protein binding sites were
analyzed by clustering of molecular dynamics (MD) simulation trajectories.
Locations of ordered waters (hydration sites) were first identified
from simulations of high resolution crystal structures of 13 protein–ligand
complexes. The MD-derived hydration sites reproduced 73% of the binding
site water molecules observed in the crystal structures. If the simulations
were repeated without the cocrystallized ligands, a majority (58%)
of the crystal waters in the binding sites were still predicted. In
addition, comparison of the hydration sites obtained from simulations
carried out in the absence of ligands to those identified for the
complexes revealed that the networks of ordered water molecules were
preserved to a large extent, suggesting that the locations of waters
in a protein–ligand interface are mainly dictated by the protein.
Analysis of >1000 crystal structures showed that hydration sites
bridged
protein–ligand interactions in complexes with different ligands,
and those with high MD-derived occupancies were more likely to correspond
to experimentally observed ordered water molecules. The results demonstrate
that ordered water molecules relevant for modeling of protein–ligand
complexes can be identified from MD simulations. Our findings could
contribute to development of improved methods for structure-based
virtual screening and lead optimization.
Collapse
Affiliation(s)
- Axel Rudling
- Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University , SE-106 91 Stockholm, Sweden
| | - Adolfo Orro
- Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University , SE-106 91 Stockholm, Sweden
| | - Jens Carlsson
- Science for Life Laboratory, Department of Cell and Molecular Biology, Uppsala University, BMC , Box 596, SE-751 24 Uppsala, Sweden
| |
Collapse
|
20
|
Jandova Z, Fast D, Setz M, Pechlaner M, Oostenbrink C. Saturation Mutagenesis by Efficient Free-Energy Calculation. J Chem Theory Comput 2018; 14:894-904. [PMID: 29262673 PMCID: PMC5813279 DOI: 10.1021/acs.jctc.7b01099] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
Single-point mutations
in proteins can greatly influence protein
stability, binding affinity, protein function or its expression per
se. Here, we present accurate and efficient predictions of the free
energy of mutation of amino acids. We divided the complete mutational
free energy into an uncharging step, which we approximate by a third-power
fitting (TPF) approach, and an annihilation step, which we approximate
using the one-step perturbation (OSP) method. As a diverse set of
test systems, we computed the solvation free energy of all amino acid
side chain analogues and obtained an excellent agreement with thermodynamic
integration (TI) data. Moreover, we calculated mutational free energies
in model tripeptides and established an efficient protocol involving
a single reference state. Again, the approximate methods agreed excellently
with the TI references, with a root-mean-square error of only 3.6
kJ/mol over 17 mutations. Our combined TPF+OSP approach does show
not only a very good agreement but also a 2-fold higher efficiency
than full blown TI calculations.
Collapse
Affiliation(s)
- Zuzana Jandova
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| | - Daniel Fast
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| | - Martina Setz
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| | - Maria Pechlaner
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| | - Chris Oostenbrink
- Institute of Molecular Modeling and Simulation, University of Natural Resources and Life Sciences , Vienna A-1190, Austria
| |
Collapse
|
21
|
Chen D, Li Y, Zhao M, Tan W, Li X, Savidge T, Guo W, Fan X. Effective lead optimization targeting the displacement of bridging receptor–ligand water molecules. Phys Chem Chem Phys 2018; 20:24399-24407. [DOI: 10.1039/c8cp04118k] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
Enhancing the binding affinities of ligands by means of lead modifications that displace bridging water molecules at protein–ligand interfaces is an important and widely studied lead optimization strategy.
Collapse
Affiliation(s)
- Deliang Chen
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Yibao Li
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Mingming Zhao
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Wen Tan
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Xun Li
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Tor Savidge
- Department of Pathology & Immunology
- Baylor College of Medicine
- Houston
- USA
- Texas Children's Microbiome Center
| | - Wei Guo
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| | - Xiaolin Fan
- Jiangxi Key Laboratory of Organo-Pharmaceutical Chemistry
- Chemistry and Chemical Engineering College
- Gannan Normal University
- Ganzhou
- P. R. China
| |
Collapse
|
22
|
Williams-Noonan BJ, Yuriev E, Chalmers DK. Free Energy Methods in Drug Design: Prospects of “Alchemical Perturbation” in Medicinal Chemistry. J Med Chem 2017; 61:638-649. [DOI: 10.1021/acs.jmedchem.7b00681] [Citation(s) in RCA: 93] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Billy J. Williams-Noonan
- Medicinal Chemistry, Monash
Institute of Pharmaceutical Sciences, Monash University, 381 Royal
Parade, Parkville, Victoria 3052, Australia
| | - Elizabeth Yuriev
- Medicinal Chemistry, Monash
Institute of Pharmaceutical Sciences, Monash University, 381 Royal
Parade, Parkville, Victoria 3052, Australia
| | - David K. Chalmers
- Medicinal Chemistry, Monash
Institute of Pharmaceutical Sciences, Monash University, 381 Royal
Parade, Parkville, Victoria 3052, Australia
| |
Collapse
|
23
|
Graf MMH, Maurer M, Oostenbrink C. Free-energy calculations of residue mutations in a tripeptide using various methods to overcome inefficient sampling. J Comput Chem 2016; 37:2597-605. [PMID: 27634475 PMCID: PMC5082540 DOI: 10.1002/jcc.24488] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2016] [Revised: 08/22/2016] [Accepted: 08/24/2016] [Indexed: 01/19/2023]
Abstract
Previous free-energy calculations have shown that the seemingly simple transformation of the tripeptide KXK to KGK in water holds some unobvious challenges concerning the convergence of the forward and backward thermodynamic integration processes (i.e., hysteresis). In the current study, the central residue X was either alanine, serine, glutamic acid, lysine, phenylalanine, or tyrosine. Interestingly, the transformation from alanine to glycine yielded the highest hysteresis in relation to the extent of the chemical change of the side chain. The reason for that could be attributed to poor sampling of φ2 /ψ2 dihedral angles along the transformation. Altering the nature of alanine's Cβ atom drastically improved the sampling and at the same time led to the identification of high energy barriers as cause for it. Consequently, simple strategies to overcome these barriers are to increase simulation time (computationally expensive) or to use enhanced sampling techniques such as Hamiltonian replica exchange molecular dynamics and one-step perturbation. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Michael M H Graf
- Institute of Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, BOKU, University of Natural Resources and Life Sciences, Muthgasse 18, AT-1190, Vienna, Austria
| | - Manuela Maurer
- Institute of Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, BOKU, University of Natural Resources and Life Sciences, Muthgasse 18, AT-1190, Vienna, Austria
| | - Chris Oostenbrink
- Institute of Molecular Modeling and Simulation, Department of Material Sciences and Process Engineering, BOKU, University of Natural Resources and Life Sciences, Muthgasse 18, AT-1190, Vienna, Austria.
| |
Collapse
|