1
|
Calculation of Crystal-Solution Dissociation Constants. Biomolecules 2022; 12:biom12020147. [PMID: 35204648 PMCID: PMC8961641 DOI: 10.3390/biom12020147] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2021] [Revised: 01/08/2022] [Accepted: 01/10/2022] [Indexed: 11/17/2022] Open
Abstract
The calculation of dissociation constants is an important problem in molecular biophysics. For such a calculation, it is important to correctly calculate both terms of the binding free energy; that is, the enthalpy and entropy of binding. Both these terms can be computed using molecular dynamics simulations, but this approach is very computationally expensive, and entropy calculations are especially slow. We develop an alternative very fast method of calculating the binding entropy and dissociation constants. The main part of our approach is based on the evaluation of movement ranges of molecules in the bound state. Then, the range of molecular movements in the bound state (here, in molecular crystals) is used for the calculation of the binding entropies and, then (using, in addition, the experimentally measured sublimation enthalpies), the crystal-to-vapor dissociation constants. Previously, we considered the process of the reversible sublimation of small organic molecules from crystals to vapor. In this work, we extend our approach by considering the dissolution of molecules, in addition to their sublimation. Similar to the sublimation case, our method shows a good correlation with experimentally measured dissociation constants at the dissolution of crystals.
Collapse
|
2
|
Borgis D, Luukkonen S, Belloni L, Jeanmairet G. Accurate prediction of hydration free energies and solvation structures using molecular density functional theory with a simple bridge functional. J Chem Phys 2021; 155:024117. [PMID: 34266282 DOI: 10.1063/5.0057506] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper assesses the ability of molecular density functional theory to predict efficiently and accurately the hydration free energies of molecular solutes and the surrounding microscopic water structure. A wide range of solutes were investigated, including hydrophobes, water as a solute, and the FreeSolv database containing 642 drug-like molecules having a variety of shapes and sizes. The usual second-order approximation of the theory is corrected by a third-order, angular-independent bridge functional. The overall functional is parameter-free in the sense that the only inputs are bulk water properties, independent of the solutes considered. These inputs are the direct correlation function, compressibility, liquid-gas surface tension, and excess chemical potential of the solvent. Compared to molecular simulations with the same force field and the same fixed solute geometries, the present theory is shown to describe accurately the solvation free energy and structure of both hydrophobic and hydrophilic solutes. Overall, the method yields a precision of order 0.5 kBT for the hydration free energies of the FreeSolv database, with a computer speedup of 3 orders of magnitude. The theory remains to be improved for a better description of the H-bonding structure and the hydration free energy of charged solutes.
Collapse
Affiliation(s)
- Daniel Borgis
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Sohvi Luukkonen
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Luc Belloni
- Universié Paris-Saclay, CEA, CNRS, NIMBE, 91191 Gif-sur-Yvette, France
| | - Guillaume Jeanmairet
- Sorbonne Université, CNRS, Physico-Chimie des Électrolytes et Nanosystèmes Interfaciaux, PHENIX, F-75005 Paris, France
| |
Collapse
|
3
|
Gobbo D, Ballone P, Garabato BD. Coarse-Grained Model of Entropy-Driven Demixing. J Phys Chem B 2020; 124:9267-9274. [DOI: 10.1021/acs.jpcb.0c07575] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- D. Gobbo
- Computational and Chemical Biology, Fondazione Istituto Italiano di Tecnologia, Genova 16163, Italy
| | - P. Ballone
- School of Physics, University College, Dublin 4, Ireland
- Conway Institute for Biomolecular and Biomedical Research, University College, Dublin 4, Ireland
| | - B. D. Garabato
- Computational and Chemical Biology, Fondazione Istituto Italiano di Tecnologia, Genova 16163, Italy
| |
Collapse
|
4
|
Luukkonen S, Levesque M, Belloni L, Borgis D. Hydration free energies and solvation structures with molecular density functional theory in the hypernetted chain approximation. J Chem Phys 2020; 152:064110. [DOI: 10.1063/1.5142651] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Sohvi Luukkonen
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Maximilien Levesque
- PASTEUR, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| | - Luc Belloni
- LIONS, NIMBE, CEA, CNRS, Université Paris-Saclay, 91191 Gif-sur-Yvette, France
| | - Daniel Borgis
- Maison de la Simulation, USR 3441 CNRS-CEA-Université Paris-Saclay, 91191 Gif-sur-Yvette, France
- PASTEUR, Département de Chimie, École Normale Supérieure, PSL University, Sorbonne Université, CNRS, 75005 Paris, France
| |
Collapse
|
5
|
Sabe VT, Tolufashe GF, Ibeji CU, Maseko SB, Govender T, Maguire GEM, Lamichhane G, Honarparvar B, Kruger HG. Identification of potent L,D-transpeptidase 5 inhibitors for Mycobacterium tuberculosis as potential anti-TB leads: virtual screening and molecular dynamics simulations. J Mol Model 2019; 25:328. [PMID: 31656981 DOI: 10.1007/s00894-019-4196-z] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2018] [Accepted: 08/28/2019] [Indexed: 11/26/2022]
Abstract
Virtual screening is a useful in silico approach to identify potential leads against various targets. It is known that carbapenems (doripenem and faropenem) do not show any reasonable inhibitory activities against L,D-transpeptidase 5 (LdtMt5) and also an adduct of meropenem exhibited slow acylation. Since these drugs are active against L,D-transpeptidase 2 (LdtMt2), understanding the differences between these two enzymes is essential. In this study, a ligand-based virtual screening of 12,766 compounds followed by molecular dynamics (MD) simulations was applied to identify potential leads against LdtMt5. To further validate the obtained virtual screening ranking for LdtMt5, we screened the same libraries of compounds against LdtMt2 which had more experimetal and calculated binding energies reported. The observed consistency between the binding affinities of LdtMt2 validates the obtained virtual screening binding scores for LdtMt5. We subjected 37 compounds with docking scores ranging from - 7.2 to - 9.9 kcal mol-1 obtained from virtual screening for further MD analysis. A set of compounds (n = 12) from four antibiotic classes with ≤ - 30 kcal mol-1 molecular mechanics/generalized born surface area (MM-GBSA) binding free energies (ΔGbind) was characterized. A final set of that, all β-lactams (n = 4), was considered. The outcome of this study provides insight into the design of potential novel leads for LdtMt5. Graphical abstract.
Collapse
Affiliation(s)
- Victor T Sabe
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Gideon F Tolufashe
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Collins U Ibeji
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Sibusiso B Maseko
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Thavendran Govender
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Glenn E M Maguire
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa
- School of Chemistry and Physics, University of KwaZulu-Natal, Durban, 4001, South Africa
| | - Gyanu Lamichhane
- Center for Tuberculosis Research, Division of Infectious Diseases, School of Medicine, Johns Hopkins University, Baltimore, MD, 21205, USA
| | - Bahareh Honarparvar
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.
| | - Hendrik G Kruger
- Catalysis and Peptide Research Unit, School of Health Sciences, University of KwaZulu-Natal, Durban, 4001, South Africa.
| |
Collapse
|
6
|
Hayashi T, Matsuda T, Nagata T, Katahira M, Kinoshita M. Mechanism of protein-RNA recognition: analysis based on the statistical mechanics of hydration. Phys Chem Chem Phys 2019; 20:9167-9180. [PMID: 29560998 DOI: 10.1039/c8cp00155c] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
We investigate the RBD1-r(GUAGU) binding as a case study using all-atom models for the biomolecules, molecular models for water, and the currently most reliable statistical-mechanical method. RBD1 is one of the RNA-binding domains of mammalian Musashi1 (Msi1), and r(GUAGU) contains the minimum recognition sequence for Msi1, r(GUAG). We show that the binding is driven by a large gain of configurational entropy of water in the entire system. It is larger than the sum of conformational-entropy losses for RBD1 and r(GUAGU). The decrease in RBD1-r(GUAGU) interaction energy upon binding is largely cancelled out by the increase in the sum of RBD1-water, r(GUAGU)-water, and water-water interaction energies. We refer to this increase as "energetic dehydration". The decrease is larger than the increase for the van der Waals component, whereas the opposite is true for the electrostatic component. We give a novel reason for the empirically known fact that protein residues possessing side chains with positive charges and with flat moieties frequently appear within protein-RNA binding interfaces. A physical picture of the general protein-RNA binding mechanism is then presented. To achieve a sufficiently large water-entropy gain, shape complementarity at the atomic level needs to be constructed by utilizing the stacking and sandwiching of flat moieties (aromatic rings of the protein and nucleobases of RNA) as fundamental motifs. To compensate for electrostatic energetic dehydration, charge complementarity becomes crucial within the binding interface. We argue the reason why the RNA recognition motif (RRM) is the most ubiquitous RNA binding domain.
Collapse
Affiliation(s)
- Tomohiko Hayashi
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan.
| | - Tomoaki Matsuda
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan.
| | - Takashi Nagata
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan.
| | - Masato Katahira
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan.
| | - Masahiro Kinoshita
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan.
| |
Collapse
|
7
|
Garbuzynskiy SO, Finkelstein AV. Evaluation of the Accuracy of Calculation of the Standard Binding Entropy of Molecules from their Average Mobility in Molecular Crystals. Mol Biol 2018. [DOI: 10.1134/s0026893318010053] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
8
|
Garbuzynskiy SO, Finkelstein AV. Sublimation Entropy and Dissociation Constants Prediction by Quantitative Evaluation of Molecular Mobility in Crystals. J Phys Chem Lett 2017; 8:2758-2763. [PMID: 28558247 DOI: 10.1021/acs.jpclett.7b00915] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Prediction of binding free energies (or dissociation constants) is a crucial challenge for computational biochemistry. One of the main problems here consists in fast and accurate evaluation of binding entropy, which is far more time-consuming than evaluation of binding enthalpy. Here, we offer a fast and rather accurate approach to evaluate the sublimation entropy (i.e., entropy of binding of a vapor molecule to a crystal, taken with the opposite sign) from the average range of molecular movements in the solid state. To estimate this range (and the corresponding amplitude), we considered equilibrium sublimation of small organic molecules from molecular crystals. The calculations were based on experimental data on the sublimation enthalpy, pressure of saturated vapor, and structural characteristics of the molecule in question. The resulting average amplitude (0.17 ± 0.01 Å) of molecular movements in crystals was used to predict sublimation entropies and dissociation constants for sublimation of 28 molecular crystals. The results of these predictions are in close agreement with experimental values.
Collapse
Affiliation(s)
- Sergiy O Garbuzynskiy
- Laboratory of Protein Physics, Institute of Protein Research, Russian Academy of Sciences , 4 Institutskaya Street, 142290 Pushchino, Moscow Region, Russia
| | - Alexei V Finkelstein
- Laboratory of Protein Physics, Institute of Protein Research, Russian Academy of Sciences , 4 Institutskaya Street, 142290 Pushchino, Moscow Region, Russia
| |
Collapse
|
9
|
Oshima H, Hayashi T, Kinoshita M. Statistical Thermodynamics for Actin-Myosin Binding: The Crucial Importance of Hydration Effects. Biophys J 2017; 110:2496-2506. [PMID: 27276267 DOI: 10.1016/j.bpj.2016.05.006] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/19/2016] [Revised: 05/03/2016] [Accepted: 05/04/2016] [Indexed: 01/01/2023] Open
Abstract
Actomyosin is an important molecular motor, and the binding of actin and myosin is an essential research target in biophysics. Nevertheless, the physical factors driving or opposing the binding are still unclear. Here, we investigate the role of water in actin-myosin binding using the most reliable statistical-mechanical method currently available for assessing biomolecules immersed in water. This method is characterized as follows: water is treated not as a dielectric continuum but as an ensemble of molecules; the polyatomic structures of proteins are taken into consideration; and the binding free energy is decomposed into physically insightful entropic and energetic components by accounting for the hydration effect to its full extent. We find that the actin-myosin binding brings large gains of electrostatic and Lennard-Jones attractive interactions. However, these gains are accompanied by even larger losses of actin-water and myosin-water electrostatic and LJ attractive interactions. Although roughly half of the energy increase due to the losses is cancelled out by the energy decrease arising from structural reorganization of the water released upon binding, the remaining energy increase is still larger than the energy decrease brought by the gains mentioned above. Hence, the net change in system energy is positive, which opposes binding. Importantly, the binding is driven by a large gain of configurational entropy of water, which surpasses the positive change in system energy and the conformational entropy loss occurring for actin and myosin. The principal physical origin of the large water-entropy gain is as follows: the actin-myosin interface is closely packed with the achievement of high shape complementarity on the atomic level, leading to a large increase in the total volume available to the translational displacement of water molecules in the system and a resultant reduction of water crowding (i.e., entropic correlations among water molecules).
Collapse
Affiliation(s)
- Hiraku Oshima
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto, Japan
| | - Tomohiko Hayashi
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto, Japan
| | | |
Collapse
|
10
|
Sorin EJ, Alvarado W, Cao S, Radcliffe A, La P, An Y. Ensemble Molecular Dynamics of a Protein-Ligand Complex: Residual Inhibitor Entropy Enhances Drug Potency in Butyrylcholinesterase. ACTA ACUST UNITED AC 2017; 6. [PMID: 28944107 DOI: 10.4172/2167-7662.1000145] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Abstract
Butyrylcholinesterase is a key enzyme that catalyzes the hydrolysis of the neurotransmitter acetylcholine and shows an increased activity in patients suffering from Alzheimer's disease (AD), making this enzyme a primary target in treating AD. Central to this problem, and to similar scenarios involving biomolecular recognition, is our understanding of the nature of the protein-ligand complex. The butyrylcholinesterase enzyme was studied via all-atom, explicit solvent, ensemble molecular dynamics simulations sans inhibitor and in the presence of three dialkyl phenyl phosphate inhibitors of known potency to a cumulative sampling of over 40 μs. Following the relaxation of these ensembles to conformational equilibria, binding modes for each inhibitor were identified. While classical models, which assume significant reduction in protein and ligand conformational entropies, continue to be favored in contemporary studies, our observations contradict those assumptions: bound ligands occupy many conformational states, thereby stabilizing the complex, while also promoting protein flexibility.
Collapse
Affiliation(s)
- Eric J Sorin
- Department of Chemistry and Biochemistry, California State University Long Beach, California, USA
| | - Walter Alvarado
- Department of Physics and Astronomy, California State University Long Beach, California, USA
| | - Samantha Cao
- Department of Chemistry and Biochemistry, California State University Long Beach, California, USA
| | - Amethyst Radcliffe
- Department of Physics and Astronomy, California State University Long Beach, California, USA
| | - Phuc La
- Department of Chemistry and Biochemistry, California State University Long Beach, California, USA
| | - Yi An
- Department of Chemistry and Biochemistry, California State University Long Beach, California, USA
| |
Collapse
|
11
|
Yan Z, Wang J. Scoring Functions of Protein-Ligand Interactions. Oncology 2017. [DOI: 10.4018/978-1-5225-0549-5.ch036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Scoring function of protein-ligand interactions is used to recognize the “native” binding pose of a ligand on the protein and to predict the binding affinity, so that the active small molecules can be discriminated from the non-active ones. Scoring function is widely used in computationally molecular docking and structure-based drug discovery. The development and improvement of scoring functions have broad implications in pharmaceutical industry and academic research. During the past three decades, much progress have been made in methodology and accuracy for scoring functions, and many successful cases have be witnessed in virtual database screening. In this chapter, the authors introduced the basic types of scoring functions and their derivations, the commonly-used evaluation methods and benchmarks, as well as the underlying challenges and current solutions. Finally, the authors discussed the promising directions to improve and develop scoring functions for future molecular docking-based drug discovery.
Collapse
|
12
|
Garbuzynskiy SO, Finkelstein AV. Calculation of mobility and entropy of the binding of molecules by crystals. Mol Biol 2016. [DOI: 10.1134/s0026893316020060] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
13
|
Hayashi T, Oshima H, Yasuda S, Kinoshita M. Mechanism of One-to-Many Molecular Recognition Accompanying Target-Dependent Structure Formation: For the Tumor Suppressor p53 Protein as an Example. J Phys Chem B 2015; 119:14120-9. [PMID: 26421917 DOI: 10.1021/acs.jpcb.5b08513] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The new type of molecular recognition, in which an intrinsically disordered region (IDR) of a protein binds to many different target proteins with target-dependent structure formation, is indispensable to the expression of life phenomena and also implicated in a number of diseases. According to the prevailing view, the physicochemical factors responsible for the binding are also target dependent. Here we consider an IDR of the tumor suppressor p53 protein, p53CTD, as an important example related to carcinogenesis and analyze its binding to four targets accompanying the formation of target-dependent structures (i.e., helix, sheet, and two different coils) using our statistical-mechanical method combined with molecular models for water. We find that all of the seemingly different binding processes are driven by a large gain of the translational, configurational entropy of water in the system. The gain originates from sufficiently high shape complementarity on the atomic level within the p53CTD-target interface. It is also required that the electrostatic complementarity be ensured as much as possible to compensate for the dehydration. Such complementarities are achieved in harmony with the portion of the target to which p53CTD binds, leading to a large diversity of structures of p53CTD formed upon binding: If they are not achievable, the binding does not occur. This finding is made possible only by calculating the changes in thermodynamic quantities upon binding and decomposing them into physically insightful components.
Collapse
Affiliation(s)
- Tomohiko Hayashi
- Institute of Advanced Energy, Kyoto University , Uji, Kyoto 611-0011, Japan
| | - Hiraku Oshima
- Institute of Advanced Energy, Kyoto University , Uji, Kyoto 611-0011, Japan
| | - Satoshi Yasuda
- Institute of Advanced Energy, Kyoto University , Uji, Kyoto 611-0011, Japan
| | - Masahiro Kinoshita
- Institute of Advanced Energy, Kyoto University , Uji, Kyoto 611-0011, Japan
| |
Collapse
|
14
|
Hayashi T, Chiba S, Kaneta Y, Furuta T, Sakurai M. ATP-induced conformational changes of nucleotide-binding domains in an ABC transporter. Importance of the water-mediated entropic force. J Phys Chem B 2014; 118:12612-20. [PMID: 25302667 DOI: 10.1021/jp507930e] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
ATP binding cassette (ABC) proteins belong to a superfamily of active transporters. Recent experimental and computational studies have shown that binding of ATP to the nucleotide binding domains (NBDs) of ABC proteins drives the dimerization of NBDs, which, in turn, causes large conformational changes within the transmembrane domains (TMDs). To elucidate the active substrate transport mechanism of ABC proteins, it is first necessary to understand how the NBD dimerization is driven by ATP binding. In this study, we selected MalKs (NBDs of a maltose transporter) as a representative NBD and calculated the free-energy change upon dimerization using molecular mechanics calculations combined with a statistical thermodynamic theory of liquids, as well as a method to calculate the translational, rotational, and vibrational entropy change. This combined method is applied to a large number of snapshot structures obtained from molecular dynamics simulations containing explicit water molecules. The results suggest that the NBD dimerization proceeds with a large gain of water entropy when ATP molecules bind to the NBDs. The energetic gain arising from direct NBD-NBD interactions is canceled by the dehydration penalty and the configurational-entropy loss. ATP hydrolysis induces a loss of the shape complementarity between the NBDs, which leads to the dissociation of the dimer, due to a decrease in the water-entropy gain and an increase in the configurational-entropy loss. This interpretation of the NBD dimerization mechanism in concert with ATP, especially focused on the water-mediated entropy force, is potentially applicable to a wide variety of the ABC transporters.
Collapse
Affiliation(s)
- Tomohiko Hayashi
- Center for Biological Resources and Informatics, Tokyo Institute of Technology , 4259-B-62, Nagatsuta-cho, Midori-ku, Yokohama 226-8501, Japan
| | | | | | | | | |
Collapse
|
15
|
Suárez D, Díaz N. Direct methods for computing single-molecule entropies from molecular simulations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2014. [DOI: 10.1002/wcms.1195] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Affiliation(s)
- Dimas Suárez
- Departamento de Química Física y Analítica; Universidad de Oviedo; Oviedo Spain
| | - Natalia Díaz
- Departamento de Química Física y Analítica; Universidad de Oviedo; Oviedo Spain
| |
Collapse
|
16
|
Mikulskis P, Genheden S, Ryde U. Effect of explicit water molecules on ligand-binding affinities calculated with the MM/GBSA approach. J Mol Model 2014; 20:2273. [DOI: 10.1007/s00894-014-2273-x] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/12/2013] [Accepted: 04/24/2014] [Indexed: 12/24/2022]
|
17
|
Hayashi T, Oshima H, Mashima T, Nagata T, Katahira M, Kinoshita M. Binding of an RNA aptamer and a partial peptide of a prion protein: crucial importance of water entropy in molecular recognition. Nucleic Acids Res 2014; 42:6861-75. [PMID: 24803670 PMCID: PMC4066790 DOI: 10.1093/nar/gku382] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
It is a central issue to elucidate the new type of molecular recognition accompanied by a global structural change of a molecule upon binding to its targets. Here we investigate the driving force for the binding of R12 (a ribonucleic acid aptamer) and P16 (a partial peptide of a prion protein) during which P16 exhibits the global structural change. We calculate changes in thermodynamic quantities upon the R12–P16 binding using a statistical-mechanical approach combined with molecular models for water which is currently best suited to studies on hydration of biomolecules. The binding is driven by a water-entropy gain originating primarily from an increase in the total volume available to the translational displacement of water molecules in the system. The energy decrease due to the gain of R12–P16 attractive (van der Waals and electrostatic) interactions is almost canceled out by the energy increase related to the loss of R12–water and P16–water attractive interactions. We can explain the general experimental result that stacking of flat moieties, hydrogen bonding and molecular-shape and electrostatic complementarities are frequently observed in the complexes. It is argued that the water-entropy gain is largely influenced by the geometric characteristics (overall shapes, sizes and detailed polyatomic structures) of the biomolecules.
Collapse
Affiliation(s)
- Tomohiko Hayashi
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Hiraku Oshima
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Tsukasa Mashima
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Takashi Nagata
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Masato Katahira
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| | - Masahiro Kinoshita
- Institute of Advanced Energy, Kyoto University, Uji, Kyoto 611-0011, Japan
| |
Collapse
|
18
|
|
19
|
Saito H, Iwayama M, Mizukami T, Kang J, Tateno M, Nagao H. Molecular dynamics study on binding free energy of Azurin–Cytochrome c551 complex. Chem Phys Lett 2013. [DOI: 10.1016/j.cplett.2012.12.016] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/27/2022]
|
20
|
Application of Hydration Thermodynamics to the Evaluation of Protein Structures and Protein-Ligand Binding. ENTROPY 2012. [DOI: 10.3390/e14081443] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|