1
|
Wu N, Lin Q, Shao F, Chen L, Zhang H, Chen K, Wu J, Wang G, Wang H, Yang Q. Insect cuticle-inspired design of sustainably sourced composite bioplastics with enhanced strength, toughness and stretch-strengthening behavior. Carbohydr Polym 2024; 333:121970. [PMID: 38494224 DOI: 10.1016/j.carbpol.2024.121970] [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/21/2023] [Revised: 02/02/2024] [Accepted: 02/17/2024] [Indexed: 03/19/2024]
Abstract
Insect cuticles that are mainly made of chitin, chitosan and proteins provide insects with rigid, stretchable and robust skins to defend harsh external environment. The insect cuticle therefore provides inspiration for engineering biomaterials with outstanding mechanical properties but also sustainability and biocompatibility. We herein propose a design of high-performance and sustainable bioplastics via introducing CPAP3-A1, a major structural protein in insect cuticles, to specifically bind to chitosan. Simply mixing 10w/w% bioengineered CPAP3-A1 protein with chitosan enables the formation of plastics-like, sustainably sourced chitosan/CPAP3-A1 composites with significantly enhanced strength (∼90 MPa) and toughness (∼20 MJ m -3), outperforming previous chitosan-based composites and most synthetic petroleum-based plastics. Remarkably, these bioplastics exhibit a stretch-strengthening behavior similar to the training living muscles. Mechanistic investigation reveals that the introduction of CPAP3-A1 induce chitosan chains to assemble into a more coarsened fibrous network with increased crystallinity and reinforcement effect, but also enable energy dissipation via reversible chitosan-protein interactions. Further uniaxial stretch facilitates network re-orientation and increases chitosan crystallinity and mechanical anisotropy, thereby resulting in stretch-strengthening behavior. In general, this study provides an insect-cuticle inspired design of high-performance bioplastics that may serve as sustainable and bio-friendly materials for a wide range of engineering and biomedical application potentials.
Collapse
Affiliation(s)
- Nan Wu
- School of Bioengineering, Dalian University of Technology, Dalian, Liaoning 116024, China
| | - Qiaoxia Lin
- School of Bioengineering, Dalian University of Technology, Dalian, Liaoning 116024, China
| | - Fei Shao
- School of Bioengineering, Dalian University of Technology, Dalian, Liaoning 116024, China
| | - Lei Chen
- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China; Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Haoyue Zhang
- School of Bioengineering, Dalian University of Technology, Dalian, Liaoning 116024, China
| | - Kaiwen Chen
- School of Bioengineering, Dalian University of Technology, Dalian, Liaoning 116024, China
| | - Jinrong Wu
- State Key Laboratory of Polymer Materials Engineering, College of Polymer Science and Engineering, Sichuan University, Chengdu 610065, China
| | - Guirong Wang
- State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China; Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China
| | - Huanan Wang
- School of Bioengineering, Dalian University of Technology, Dalian, Liaoning 116024, China.
| | - Qing Yang
- School of Bioengineering, Dalian University of Technology, Dalian, Liaoning 116024, China; State Key Laboratory for Biology of Plant Diseases and Insect Pests, Institute of Plant Protection, Chinese Academy of Agricultural Sciences, Beijing 100193, China; Shenzhen Branch, Guangdong Laboratory of Lingnan Modern Agriculture, Key Laboratory of Synthetic Biology, Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518120, China.
| |
Collapse
|
2
|
Dance I. The binding of reducible N 2 in the reaction domain of nitrogenase. Dalton Trans 2023; 52:2013-2026. [PMID: 36691966 DOI: 10.1039/d2dt03599e] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
The binding of N2 to FeMo-co, the catalytic site of the enzyme nitrogenase, is central to the conversion to NH3, but also has a separate role in promoting the N2-dependent HD reaction (D2 + 2H+ + 2e- → 2HD). The protein surrounding FeMo-co contains a clear channel for ingress of N2, directly towards the exo-coordination position of Fe2, a position which is outside the catalytic reaction domain. This led to the hypothesis [I. Dance, Dalton Trans., 2022, 51, 12717] of 'promotional' N2 bound at exo-Fe2, and a second 'reducible' N2 bound in the reaction domain, specifically the endo-coordination position of Fe2 or Fe6. The range of possibilities for the binding of reducible N2 in the presence of bound promotional N2 is described here, using density functional simulations with a 486 atom model of the active site and surrounding protein. The pathway for ingress of the second N2 through protein, past the first N2 at exo-Fe2, and tumbling into the binding domain between Fe2 and Fe6, is described. The calculations explore 24 structures involving 6 different forms of hydrogenated FeMo-co, including structures with S2BH unhooked from Fe2 but tethered to Fe6. The calculations use the most probable electronic states. End-on (η1) binding of N2 at the endo position of either Fe2 or Fe6 is almost invariably exothermic, with binding potential energies ranging up to -18 kcal mol-1. Many structures have binding energies in the range -6 to -14 kcal mol-1. The relevant entropic penalty for N2 binding from a diffusible position within the protein is estimated to be 4 kcal mol-1, and so the binding free energies for reducible N2 are suitably negative. N2 binding at endo-Fe2 is stronger than at endo-Fe6 in three of the six structure categories. In many cases the reaction domain containing reducible N2 is expanded. These results inform computational simulation of the subsequent steps in which surrounding H atoms transfer to reducible N2.
Collapse
Affiliation(s)
- Ian Dance
- School of Chemistry, UNSW Sydney, Australia.
| |
Collapse
|
3
|
Dance I. The HD Reaction of Nitrogenase: a Detailed Mechanism. Chemistry 2023; 29:e202202502. [PMID: 36274057 PMCID: PMC10099629 DOI: 10.1002/chem.202202502] [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/11/2022] [Indexed: 11/06/2022]
Abstract
Nitrogenase is the enzyme that converts N2 to NH3 under ambient conditions. The chemical mechanism of this catalysis at the active site FeMo-co [Fe7 S9 CMo(homocitrate)] is unknown. An obligatory co-product is H2 , while exogenous H2 is a competitive inhibitor. Isotopic substitution using exogenous D2 revealed the N2 -dependent reaction D2 +2H+ +2e- →2HD (the 'HD reaction'), together with a collection of additional experimental characteristics and requirements. This paper describes a detailed mechanism for the HD reaction, developed and elaborated using density functional simulations with a 486-atom model of the active site and surrounding protein. First D2 binds at one Fe atom (endo-Fe6 coordination position), where it is flanked by H-Fe6 (exo position) and H-Fe2 (endo position). Then there is synchronous transfer of these two H atoms to bound D2 , forming one HD bound to Fe2 and a second HD bound to Fe6. These two HD dissociate sequentially. The final phase is recovery of the two flanking H atoms. These H atoms are generated, sequentially, by translocation of a proton from the protein surface to S3B of FeMo-co and combination with introduced electrons. The first H atom migrates from S3B to exo-Fe6 and the second from S3B to endo-Fe2. Reaction energies and kinetic barriers are reported for all steps. This mechanism accounts for the experimental data: (a) stoichiometry; (b) the N2 -dependence results from promotional N2 bound at exo-Fe2; (c) different N2 binding Km for the HD reaction and the NH3 formation reaction results from involvement of two different sites; (d) inhibition by CO; (e) the non-occurrence of 2HD→H2 +D2 results from the synchronicity of the two transfers of H to D2 ; (f) inhibition of HD production at high pN2 is by competitive binding of N2 at endo-Fe6; (g) the non-leakage of D to solvent follows from the hydrophobic environment and irreversibility of proton introduction.
Collapse
Affiliation(s)
- Ian Dance
- School of ChemistryUNSWSydneyAustralia
| |
Collapse
|
4
|
Kalayan J, Chakravorty A, Warwicker J, Henchman RH. Total free energy analysis of fully hydrated proteins. Proteins 2023; 91:74-90. [PMID: 35964252 PMCID: PMC10087023 DOI: 10.1002/prot.26411] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2022] [Revised: 08/04/2022] [Accepted: 08/09/2022] [Indexed: 12/15/2022]
Abstract
The total free energy of a hydrated biomolecule and its corresponding decomposition of energy and entropy provides detailed information about regions of thermodynamic stability or instability. The free energies of four hydrated globular proteins with different net charges are calculated from a molecular dynamics simulation, with the energy coming from the system Hamiltonian and entropy using multiscale cell correlation. Water is found to be most stable around anionic residues, intermediate around cationic and polar residues, and least stable near hydrophobic residues, especially when more buried, with stability displaying moderate entropy-enthalpy compensation. Conversely, anionic residues in the proteins are energetically destabilized relative to singly solvated amino acids, while trends for other residues are less clear-cut. Almost all residues lose intraresidue entropy when in the protein, enthalpy changes are negative on average but may be positive or negative, and the resulting overall stability is moderate for some proteins and negligible for others. The free energy of water around single amino acids is found to closely match existing hydrophobicity scales. Regarding the effect of secondary structure, water is slightly more stable around loops, of intermediate stability around β strands and turns, and least stable around helices. An interesting asymmetry observed is that cationic residues stabilize a residue when bonded to its N-terminal side but destabilize it when on the C-terminal side, with a weaker reversed trend for anionic residues.
Collapse
Affiliation(s)
- Jas Kalayan
- Division of Pharmacy and Optometry, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, UK
| | - Arghya Chakravorty
- Department of Chemistry and Biophysics, University of Michigan, Ann Arbor, Michigan, USA
| | - Jim Warwicker
- Manchester Institute of Biotechnology and School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, UK
| | - Richard H Henchman
- Sydney Medical School, Faculty of Medicine and Health, University of Sydney, Sydney, Australia
| |
Collapse
|
5
|
Esteve F, Altava B, Luis SV, García-Verdugo E. Basically, nucleophilicity matters little: towards unravelling the supramolecular driving forces in enzyme-like CO 2 conversion. Org Biomol Chem 2022; 20:6637-6645. [PMID: 35929502 DOI: 10.1039/d2ob00948j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The reaction mechanism for the cycloaddition of CO2 to styrene oxide in the presence of macrocyclic pseudopeptides has been studied using DFT methods. Computational calculations indicate that the unprecedented catalytic behaviour previously observed experimentally, in which the most reactive species was not the most nucleophilic but the most basic one, can be associated to the tight cooperativity between several supramolecular interactions promoted by simple peptidomimetics able to display a synzymatic behaviour. This bizarre catalytic performance afforded remarkable conversions of a sluggish substrate like styrene oxide into the desired cyclic carbonate, even under relatively mild reaction conditions, opening the way for the practical use of CO2 as a raw material in the preparation of valuable chemicals. Furthermore, the remote modification of essential structural features of the macrocycle (synzyme engineering) permitted the driving forces of the synzymatic system to be analyzed, stressing the crucial synergic effect between an elegantly preorganized oxyanion hole and additional aromatic interactions.
Collapse
Affiliation(s)
- Ferran Esteve
- Departamento de Química Inorgánica y Orgánica, Universitat Jaume I, Av. Sos Baynat s/n, Castellón, 12071, Spain.
| | - Belén Altava
- Departamento de Química Inorgánica y Orgánica, Universitat Jaume I, Av. Sos Baynat s/n, Castellón, 12071, Spain.
| | - Santiago V Luis
- Departamento de Química Inorgánica y Orgánica, Universitat Jaume I, Av. Sos Baynat s/n, Castellón, 12071, Spain.
| | - Eduardo García-Verdugo
- Departamento de Química Inorgánica y Orgánica, Universitat Jaume I, Av. Sos Baynat s/n, Castellón, 12071, Spain.
| |
Collapse
|
6
|
Dance I. Structures and reaction dynamics of N 2 and H 2 binding at FeMo-co, the active site of nitrogenase. Dalton Trans 2021; 50:18212-18237. [PMID: 34860237 DOI: 10.1039/d1dt03548g] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
The chemical reactions occurring at the Fe7MoS9C(homocitrate) cluster, FeMo-co, the active site of the enzyme nitrogenase (N2 → NH3), are enigmatic. Experimental information collected over a long period reveals aspects of the roles of N2 and H2, each with more than one type of reactivity. This paper reports investigations of the binding of H2 and N2 at intact FeMo-co, using density functional simulations of a large 486 atom relevant portion of the protein, resulting in 27 new structures containing H2 and/or N2 bound at the exo and endo coordination sites of the participating Fe atoms, Fe2 and Fe6. Binding energies and transition states for association/dissociation are determined, and trajectories for the approach, binding and separation of H2/N2 are described, including diffusion of these small molecules through proximal protein. Influences of surrounding amino acids are identified. FeMo-co deforms geometrically when binding H2 or N2, and a procedure for calculating the energy cost involved, the adaptation energy, is introduced here. Adaptation energies, which range from 7 to 36 kcal mol-1 for the reported structures, are influenced by the protonation state of the His195 side chain. Seven N2 structures and three H2 structures have negative binding free energies, which include the estimated entropy penalties for binding of N2, H2 from proximal protein. These favoured structures have N2 bound end-on at exo-Fe2, exo-Fe6 and endo-Fe2 positions of FeMo-co, and H2 bound at the endo-Fe2 position. Various postulated structures with N2 bridging Fe2 and Fe6 revert to end-on-N2 at endo positions. The structures are also assessed via the calculated potential energy barriers for association and dissociation. Barriers to the binding of H2 range from 1 to 20 kcal mol-1 and barriers to dissociation of H2 range from 3 to 18 kcal mol-1. Barriers to the binding of N2, in either side-on or end-on mode, range from 2 to 18 kcal mol-1, while dissociation of bound N2 encounters barriers of 3 to 8 kcal mol-1 for side-on bonding and 7 to 18 kcal mol-1 for end-on bonding. These results allow formulation of mechanisms for the H2/N2 exchange reaction, and three feasible mechanisms for associative exchange and three for dissociative exchange are identified. Consistent electronic structures and potential energy surfaces are maintained throughout. Changes in the spin populations of Fe2 and Fe6 connected with cluster deformation and with metal-ligand bond formation are identified.
Collapse
Affiliation(s)
- Ian Dance
- School of Chemistry, UNSW Sydney, NSW 2051, Australia.
| |
Collapse
|
7
|
Ali HS, Chakravorty A, Kalayan J, de Visser SP, Henchman RH. Energy-entropy method using multiscale cell correlation to calculate binding free energies in the SAMPL8 host-guest challenge. J Comput Aided Mol Des 2021; 35:911-921. [PMID: 34264476 PMCID: PMC8367938 DOI: 10.1007/s10822-021-00406-5] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2021] [Accepted: 06/22/2021] [Indexed: 11/29/2022]
Abstract
Free energy drives a wide range of molecular processes such as solvation, binding, chemical reactions and conformational change. Given the central importance of binding, a wide range of methods exist to calculate it, whether based on scoring functions, machine-learning, classical or electronic structure methods, alchemy, or explicit evaluation of energy and entropy. Here we present a new energy-entropy (EE) method to calculate the host-guest binding free energy directly from molecular dynamics (MD) simulation. Entropy is evaluated using Multiscale Cell Correlation (MCC) which uses force and torque covariance and contacts at two different length scales. The method is tested on a series of seven host-guest complexes in the SAMPL8 (Statistical Assessment of the Modeling of Proteins and Ligands) "Drugs of Abuse" Blind Challenge. The EE-MCC binding free energies are found to agree with experiment with an average error of 0.9 kcal mol-1. MCC makes clear the origin of the entropy changes, showing that the large loss of positional, orientational, and to a lesser extent conformational entropy of each binding guest is compensated for by a gain in orientational entropy of water released to bulk, combined with smaller decreases in vibrational entropy of the host, guest and contacting water.
Collapse
Affiliation(s)
- Hafiz Saqib Ali
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| | - Arghya Chakravorty
- Department of Chemistry, University of Michigan, Ann Arbor, MI, 48109, USA
| | - Jas Kalayan
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| | - Samuel P de Visser
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK
- Department of Chemical Engineering and Analytical Science, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| | - Richard H Henchman
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK.
- Department of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK.
- Sydney Medical School, The University of Sydney, Sydney, NSW, 2006, Australia.
| |
Collapse
|
8
|
Falcioni F, Kalayan J, Henchman RH. Energy-entropy prediction of octanol-water logP of SAMPL7 N-acyl sulfonamide bioisosters. J Comput Aided Mol Des 2021; 35:831-840. [PMID: 34244906 PMCID: PMC8295089 DOI: 10.1007/s10822-021-00401-w] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2021] [Accepted: 06/17/2021] [Indexed: 12/23/2022]
Abstract
Partition coefficients quantify a molecule's distribution between two immiscible liquid phases. While there are many methods to compute them, there is not yet a method based on the free energy of each system in terms of energy and entropy, where entropy depends on the probability distribution of all quantum states of the system. Here we test a method in this class called Energy Entropy Multiscale Cell Correlation (EE-MCC) for the calculation of octanol-water logP values for 22 N-acyl sulfonamides in the SAMPL7 Physical Properties Challenge (Statistical Assessment of the Modelling of Proteins and Ligands). EE-MCC logP values have a mean error of 1.8 logP units versus experiment and a standard error of the mean of 1.0 logP units for three separate calculations. These errors are primarily due to getting sufficiently converged energies to give accurate differences of large numbers, particularly for the large-molecule solvent octanol. However, this is also an issue for entropy, and approximations in the force field and MCC theory also contribute to the error. Unique to MCC is that it explains the entropy contributions over all the degrees of freedom of all molecules in the system. A gain in orientational entropy of water is the main favourable entropic contribution, supported by small gains in solute vibrational and orientational entropy but offset by unfavourable changes in the orientational entropy of octanol, the vibrational entropy of both solvents, and the positional and conformational entropy of the solute.
Collapse
Affiliation(s)
- Fabio Falcioni
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK.
- School of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK.
| | - Jas Kalayan
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK
- School of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
| | - Richard H Henchman
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester, M1 7DN, UK.
- School of Chemistry, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK.
- Sydney Medical School, The University of Sydney, Sydney, NSW, 2006, Australia.
| |
Collapse
|
9
|
Chakravorty A, Higham J, Henchman RH. Entropy of Proteins Using Multiscale Cell Correlation. J Chem Inf Model 2020; 60:5540-5551. [PMID: 32955869 DOI: 10.1021/acs.jcim.0c00611] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A new multiscale method is presented to calculate the entropy of proteins from molecular dynamics simulations. Termed Multiscale Cell Correlation (MCC), the method decomposes the protein into sets of rigid-body units based on their covalent-bond connectivity at three levels of hierarchy: molecule, residue, and united atom. It evaluates the vibrational and topographical entropy from forces, torques, and dihedrals at each level, taking into account correlations between sets of constituent units that together make up a larger unit at the coarser length scale. MCC gives entropies in close agreement with normal-mode analysis and smaller than those using quasiharmonic analysis as well as providing much faster convergence. Moreover, MCC provides an insightful decomposition of entropy at each length scale and for each type of amino acid according to their solvent exposure and whether they are terminal residues. While the residue entropy depends weakly on solvent exposure, there is greater variation in entropy components for larger, more polar amino acids, which have increased conformational entropy but reduced vibrational entropy with greater solvent exposure.
Collapse
Affiliation(s)
- Arghya Chakravorty
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Jonathan Higham
- MRC Human Genetics Unit, Institute of Genetics & Molecular Medicine, The University of Edinburgh, Western General Hospital, Crewe Road South, Edinburgh EH4 2XU, United Kingdom
| | - Richard H Henchman
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom.,Department of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
| |
Collapse
|
10
|
Sung HL, Nesbitt DJ. Sequential Folding of the Nickel/Cobalt Riboswitch Is Facilitated by a Conformational Intermediate: Insights from Single-Molecule Kinetics and Thermodynamics. J Phys Chem B 2020; 124:7348-7360. [DOI: 10.1021/acs.jpcb.0c05625] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Affiliation(s)
- Hsuan-Lei Sung
- JILA, National Institute of Standards and Technology and University of Colorado, Boulder, Colorado 80309, United States,
- Department of Chemistry, University of Colorado, Boulder, Colorado 80309, United States
| | - David J. Nesbitt
- JILA, National Institute of Standards and Technology and University of Colorado, Boulder, Colorado 80309, United States,
- Department of Chemistry, University of Colorado, Boulder, Colorado 80309, United States
- Department of Physics, University of Colorado, Boulder, Colorado 80309, United States
| |
Collapse
|
11
|
Ali HS, Higham J, Henchman RH. Entropy of Simulated Liquids Using Multiscale Cell Correlation. ENTROPY 2019; 21:e21080750. [PMID: 33267464 PMCID: PMC7515279 DOI: 10.3390/e21080750] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/29/2019] [Revised: 07/22/2019] [Accepted: 07/28/2019] [Indexed: 12/16/2022]
Abstract
Accurately calculating the entropy of liquids is an important goal, given that many processes take place in the liquid phase. Of almost equal importance is understanding the values obtained. However, there are few methods that can calculate the entropy of such systems, and fewer still to make sense of the values obtained. We present our multiscale cell correlation (MCC) method to calculate the entropy of liquids from molecular dynamics simulations. The method uses forces and torques at the molecule and united-atom levels and probability distributions of molecular coordinations and conformations. The main differences with previous work are the consistent treatment of the mean-field cell approximation to the approriate degrees of freedom, the separation of the force and torque covariance matrices, and the inclusion of conformation correlation for molecules with multiple dihedrals. MCC is applied to a broader set of 56 important industrial liquids modeled using the Generalized AMBER Force Field (GAFF) and Optimized Potentials for Liquid Simulations (OPLS) force fields with 1.14*CM1A charges. Unsigned errors versus experimental entropies are 8.7 J K - 1 mol - 1 for GAFF and 9.8 J K - 1 mol - 1 for OPLS. This is significantly better than the 2-Phase Thermodynamics method for the subset of molecules in common, which is the only other method that has been applied to such systems. MCC makes clear why the entropy has the value it does by providing a decomposition in terms of translational and rotational vibrational entropy and topographical entropy at the molecular and united-atom levels.
Collapse
Affiliation(s)
- Hafiz Saqib Ali
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester M1 7DN, UK
- School of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Jonathan Higham
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester M1 7DN, UK
- School of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Richard H. Henchman
- Manchester Institute of Biotechnology, The University of Manchester, 131 Princess Street, Manchester M1 7DN, UK
- School of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
- Correspondence: ; Tel.: +44-161-306-5194
| |
Collapse
|
12
|
Higham J, Chou SY, Gräter F, Henchman RH. Entropy of flexible liquids from hierarchical force–torque covariance and coordination. Mol Phys 2018. [DOI: 10.1080/00268976.2018.1459002] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Jonathan Higham
- Manchester Institute of Biotechnology, The University of Manchester , Manchester, United Kingdom
- School of Chemistry, The University of Manchester , Manchester, United Kingdom
| | - Szu-Yu Chou
- Manchester Institute of Biotechnology, The University of Manchester , Manchester, United Kingdom
- School of Chemistry, The University of Manchester , Manchester, United Kingdom
| | - Frauke Gräter
- Heidelberg Institute for Theoretical Studies , Heidelberg, Germany
- Interdisciplinary Center for Scientific Computing (IWR), Mathematikon, Heidelberg University , Heidelberg, Germany
| | - Richard H. Henchman
- Manchester Institute of Biotechnology, The University of Manchester , Manchester, United Kingdom
- School of Chemistry, The University of Manchester , Manchester, United Kingdom
| |
Collapse
|
13
|
Structural characterization of the nitrogenase molybdenum-iron protein with the substrate acetylene trapped near the active site. J Inorg Biochem 2018; 180:129-134. [DOI: 10.1016/j.jinorgbio.2017.12.008] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2017] [Revised: 12/08/2017] [Accepted: 12/10/2017] [Indexed: 11/21/2022]
|
14
|
Chong SH, Ham S. Dynamics of Hydration Water Plays a Key Role in Determining the Binding Thermodynamics of Protein Complexes. Sci Rep 2017; 7:8744. [PMID: 28821854 PMCID: PMC5562873 DOI: 10.1038/s41598-017-09466-w] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2017] [Accepted: 07/25/2017] [Indexed: 01/15/2023] Open
Abstract
Interfacial waters are considered to play a crucial role in protein–protein interactions, but in what sense and why are they important? Here, using molecular dynamics simulations and statistical thermodynamic analyses, we demonstrate distinctive dynamic characteristics of the interfacial water and investigate their implications for the binding thermodynamics. We identify the presence of extraordinarily slow (~1,000 times slower than in bulk water) hydrogen-bond rearrangements in interfacial water. We rationalize the slow rearrangements by introducing the “trapping” free energies, characterizing how strongly individual hydration waters are captured by the biomolecular surface, whose magnitude is then traced back to the number of water–protein hydrogen bonds and the strong electrostatic field produced at the binding interface. We also discuss the impact of the slow interfacial waters on the binding thermodynamics. We find that, as expected from their slow dynamics, the conventional approach to the water-mediated interaction, which assumes rapid equilibration of the waters’ degrees of freedom, is inadequate. We show instead that an explicit treatment of the extremely slow interfacial waters is critical. Our results shed new light on the role of water in protein–protein interactions, highlighting the need to consider its dynamics to improve our understanding of biomolecular bindings.
Collapse
Affiliation(s)
- Song-Ho Chong
- Department of Chemistry, Sookmyung Women's University, Cheongpa-ro 47-gil 100, Yongsan-Ku, Seoul, 04310, Korea
| | - Sihyun Ham
- Department of Chemistry, Sookmyung Women's University, Cheongpa-ro 47-gil 100, Yongsan-Ku, Seoul, 04310, Korea.
| |
Collapse
|
15
|
Chong SH, Ham S. New Computational Approach for External Entropy in Protein-Protein Binding. J Chem Theory Comput 2016; 12:2509-16. [PMID: 27153451 DOI: 10.1021/acs.jctc.6b00174] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Abstract
Molecular recognition through the noncovalent association of biomolecules is of central importance in biology and pharmacology. Developing reliable computational methods for estimating binding thermodynamic parameters is therefore of great practical value. However, considerable uncertainty remains regarding the external entropy that is associated with the reduction in the external (positional and orientational) degrees of freedom upon complex formation. Here, we present a novel statistical mechanical method for computing the external entropy by extending the energetic approach we have developed for unimolecular processes to association processes. We find that, in contrary to what is postulated in most of the previous methods, intrinsic couplings between the internal and external degrees of freedom of bound complex cannot in general be neglected in the determination of the external entropy. Nevertheless, there exists the best choice of the external coordinates with which those couplings are minimized. With such a choice of the external coordinates, the lowest upper bound of the external entropy is obtained from a tractable expression, which serves as an estimate of the external entropy. Our method can be implemented in a straightforward manner with molecular dynamics simulations, and its applicability is demonstrated through the application to the barnase-barstar complex.
Collapse
Affiliation(s)
- Song-Ho Chong
- Department of Chemistry, Sookmyung Women's University , Cheongpa-ro 47-gil 100, Yongsan-Ku, Seoul 04310, Korea
| | - Sihyun Ham
- Department of Chemistry, Sookmyung Women's University , Cheongpa-ro 47-gil 100, Yongsan-Ku, Seoul 04310, Korea
| |
Collapse
|
16
|
A computational UV–Vis spectroscopic study of the chemical speciation of 2-mercaptobenzothiazole corrosion inhibitor in aqueous solution. Theor Chem Acc 2016. [DOI: 10.1007/s00214-016-1839-3] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
17
|
Estimation of the driving force for dioxygen formation in photosynthesis. BIOCHIMICA ET BIOPHYSICA ACTA-BIOENERGETICS 2016; 1857:23-33. [DOI: 10.1016/j.bbabio.2015.09.011] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/28/2015] [Revised: 09/10/2015] [Accepted: 09/30/2015] [Indexed: 11/22/2022]
|
18
|
Abstract
Detailed thermodynamic analysis of fragment binding revealed that unlike drug-like compounds, fragments bind with significant enthalpic preference.
Collapse
Affiliation(s)
- György G. Ferenczy
- Medicinal Chemistry Research Group
- Research Centre for Natural Sciences
- Hungarian Academy of Sciences
- Budapest 1117
- Hungary
| | - György M. Keserű
- Medicinal Chemistry Research Group
- Research Centre for Natural Sciences
- Hungarian Academy of Sciences
- Budapest 1117
- Hungary
| |
Collapse
|
19
|
Hensen U, Gräter F, Henchman RH. Macromolecular Entropy Can Be Accurately Computed from Force. J Chem Theory Comput 2015; 10:4777-81. [PMID: 26584364 DOI: 10.1021/ct500684w] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
A method is presented to evaluate a molecule's entropy from the atomic forces calculated in a molecular dynamics simulation. Specifically, diagonalization of the mass-weighted force covariance matrix produces eigenvalues which in the harmonic approximation can be related to vibrational frequencies. The harmonic oscillator entropies of each vibrational mode may be summed to give the total entropy. The results for a series of hydrocarbons, dialanine and a β hairpin are found to agree much better with values derived from thermodynamic integration than results calculated using quasiharmonic analysis. Forces are found to follow a harmonic distribution more closely than coordinate displacements and better capture the underlying potential energy surface. The method's accuracy, simplicity, and computational similarity to quasiharmonic analysis, requiring as input force trajectories instead of coordinate trajectories, makes it readily applicable to a wide range of problems.
Collapse
Affiliation(s)
- Ulf Hensen
- ETH Zürich , Biosystems Science and Engineering, Mattenstrasse 26, 4058 Basel, Switzerland
| | - Frauke Gräter
- Heidelberg Institute for Theoretical Studies , Schloss-Wolfsbrunnenweg 35, 69118 Heidelberg, Germany
| | - Richard H Henchman
- Manchester Institute of Biotechnology, The University of Manchester , 131 Princess Street, Manchester M1 7DN, United Kingdom.,School of Chemistry, The University of Manchester , Oxford Road, Manchester M13 9PL, United Kingdom
| |
Collapse
|
20
|
Ulucan O, Jaitly T, Helms V. Energetics of Hydrophilic Protein-Protein Association and the Role of Water. J Chem Theory Comput 2015; 10:3512-24. [PMID: 26588315 DOI: 10.1021/ct5001796] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023]
Abstract
Hydrophilic protein-protein interfaces constitute a major part of all protein-protein interfaces and are thus of great importance. However, the quantitative characterization of their association is still an ongoing challenge and the driving force behind their association remains poorly characterized. Here, we have addressed the association of hydrophilic proteins and the role of water by means of extensive molecular dynamics simulations in explicit water using three well studied protein complexes; Barnase-Barstar, cytochrome c-cytochrome c peroxidase, and the N-terminal domain of enzyme I-histidine-containing phosphocarrier. The one-dimensional free energy profiles obtained from umbrella sampling simulations are downhill or, in other words, barrierless. Using these one-dimensional free energy profiles, the computed standard free energies of binding are -12.7 ± 1.1 kcal/mol, -9.4 ± 0.7 kcal/mol, and -8.4 ± 1.9 kcal/mol that are in reasonable to very good agreement with the experimental values of -19.6 kcal/mol, -8.8 kcal/mol, and -7.8 kcal/mol. As expected, analysis of the confined water between the hydrophilic complex partners shows that the density and the orientational order parameter deviate noticeably from the bulk values, especially at close separations of the confining proteins.
Collapse
Affiliation(s)
- Ozlem Ulucan
- Center for Bioinformatics, Saarland University , Saarbruecken, Germany
| | - Tanushree Jaitly
- Center for Bioinformatics, Saarland University , Saarbruecken, Germany
| | - Volkhard Helms
- Center for Bioinformatics, Saarland University , Saarbruecken, Germany
| |
Collapse
|
21
|
Zhao C, Caplan DA, Noskov SY. Evaluations of the Absolute and Relative Free Energies for Antidepressant Binding to the Amino Acid Membrane Transporter LeuT with Free Energy Simulations. J Chem Theory Comput 2015; 6:1900-14. [PMID: 26615849 DOI: 10.1021/ct9006597] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The binding of ligands to protein receptors with high affinity and specificity is central to many cellular processes. The quest for the development of computational models capable of accurately evaluating binding affinity remains one of the main goals of modern computational biophysics. In this work, free energy perturbation/molecular dynamics simulations were used to evaluate absolute and relative binding affinity for three different antidepressants to a sodium-dependent membrane transporter, LeuT, a bacterial homologue of human serotonin and dopamine transporters. Dysfunction of these membrane transporters in mammals has been implicated in multiple diseases of the nervous system, including bipolar disorder and depression. Furthermore, these proteins are key targets for antidepressants including fluoxetine (aka Prozac) and tricyclic antidepressants known to block transport activity. In addition to being clinically relevant, this system, where multiple crystal structures are readily available, represents an ideal testing ground for methods used to study the molecular mechanisms of ligand binding to membrane proteins. We discuss possible pitfalls and different levels of approximation required to evaluate binding affinity, such as the dependence of the computed affinities on the strength of constraints and the sensitivity of the computed affinities to the particular partial charges derived from restrained electrostatic potential fitting of quantum mechanics electrostatic potential. Finally, we compare the effects of different constraint schemes on the absolute and relative binding affinities obtained from free energy simulations.
Collapse
Affiliation(s)
- Chunfeng Zhao
- Institute for Biocomplexity and Informatics and Department of Biological Sciences, University of Calgary, 2500 University Drive, BI558, Calgary, AB, Canada T2N 1N4 and Molecular Structure and Function, Hospital for Sick Children and Department of Biochemistry, University of Toronto, Ontario, Canada
| | - David A Caplan
- Institute for Biocomplexity and Informatics and Department of Biological Sciences, University of Calgary, 2500 University Drive, BI558, Calgary, AB, Canada T2N 1N4 and Molecular Structure and Function, Hospital for Sick Children and Department of Biochemistry, University of Toronto, Ontario, Canada
| | - Sergei Yu Noskov
- Institute for Biocomplexity and Informatics and Department of Biological Sciences, University of Calgary, 2500 University Drive, BI558, Calgary, AB, Canada T2N 1N4 and Molecular Structure and Function, Hospital for Sick Children and Department of Biochemistry, University of Toronto, Ontario, Canada
| |
Collapse
|
22
|
Li J, Farrokhnia M, Rulíšek L, Ryde U. Catalytic Cycle of Multicopper Oxidases Studied by Combined Quantum- and Molecular-Mechanical Free-Energy Perturbation Methods. J Phys Chem B 2015; 119:8268-84. [DOI: 10.1021/acs.jpcb.5b02864] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Jilai Li
- Department
of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
- Institute
of Theoretical Chemistry, Jilin University, Changchun 130023, People’s Republic of China
| | - Maryam Farrokhnia
- Department
of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
- The
Persian Gulf Marine Biotechnology Research Center, The Persian Gulf
Biomedical Sciences Research Institute, Bushehr University of Medical Sciences, Bushehr, Iran
| | - Lubomír Rulíšek
- Institute
of Organic Chemistry and Biochemistry, Gilead Sciences and IOCB Research
Center, Academy of Sciences of the Czech Republic, Flemingovo
náměstí 2, 166
10 Prague 6, Czech Republic
| | - Ulf Ryde
- Department
of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
23
|
Gerogiokas G, Southey MWY, Mazanetz MP, Hefeitz A, Bodkin M, Law RJ, Michel J. Evaluation of water displacement energetics in protein binding sites with grid cell theory. Phys Chem Chem Phys 2015; 17:8416-26. [DOI: 10.1039/c4cp05572a] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
The grid cell theory method was used to elucidate perturbations in water network energetics in a range of protein–ligand complexes.
Collapse
Affiliation(s)
| | | | | | | | | | | | - J. Michel
- EaStCHEM School of Chemistry
- Edinburgh
- UK
| |
Collapse
|
24
|
Sandbeck DJS, Kuntz CM, Luu C, Mondor RA, Ottaviano JG, Rayer AV, Sumon KZ, East ALL. Challenges in Predicting ΔrxnG in Solution: The Mechanism of Ether-Catalyzed Hydroboration of Alkenes. J Phys Chem A 2014; 118:11768-79. [DOI: 10.1021/jp507793v] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Daniel J. S. Sandbeck
- Department of Chemistry and
Biochemistry, University of Regina, 3737 Wascana Parkway, Regina, SK S4S
0A2, Canada
| | - Colin M. Kuntz
- Department of Chemistry and
Biochemistry, University of Regina, 3737 Wascana Parkway, Regina, SK S4S
0A2, Canada
| | - Christine Luu
- Department of Chemistry and
Biochemistry, University of Regina, 3737 Wascana Parkway, Regina, SK S4S
0A2, Canada
| | - Rachelle A. Mondor
- Department of Chemistry and
Biochemistry, University of Regina, 3737 Wascana Parkway, Regina, SK S4S
0A2, Canada
| | - John G. Ottaviano
- Department of Chemistry and
Biochemistry, University of Regina, 3737 Wascana Parkway, Regina, SK S4S
0A2, Canada
| | - Aravind V. Rayer
- Department of Industrial Systems Engineering, University of Regina, 3737 Wascana Parkway, Regina, SK S4S
0A2, Canada
| | - Kazi Z. Sumon
- Department of Industrial Systems Engineering, University of Regina, 3737 Wascana Parkway, Regina, SK S4S
0A2, Canada
| | - Allan L. L. East
- Department of Chemistry and
Biochemistry, University of Regina, 3737 Wascana Parkway, Regina, SK S4S
0A2, Canada
| |
Collapse
|
25
|
Li J, Ryde U. Comparison of the Active-Site Design of Molybdenum Oxo-Transfer Enzymes by Quantum Mechanical Calculations. Inorg Chem 2014; 53:11913-24. [DOI: 10.1021/ic5010837] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Jilai Li
- Department of Theoretical Chemistry, Chemical
Centre, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
- State Key Laboratory of Theoretical and
Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, People’s Republic of China
| | - Ulf Ryde
- Department of Theoretical Chemistry, Chemical
Centre, Lund University, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
26
|
Sasikala WD, Mukherjee A. Single water entropy: hydrophobic crossover and application to drug binding. J Phys Chem B 2014; 118:10553-64. [PMID: 25133988 DOI: 10.1021/jp502852f] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Entropy of water plays an important role in both chemical and biological processes e.g. hydrophobic effect, molecular recognition etc. Here we use a new approach to calculate translational and rotational entropy of the individual water molecules around different hydrophobic and charged solutes. We show that for small hydrophobic solutes, the translational and rotational entropies of each water molecule increase as a function of its distance from the solute reaching finally to a constant bulk value. As the size of the solute increases (0.746 nm), the behavior of the translational entropy is opposite; water molecules closest to the solute have higher entropy that reduces with distance from the solute. This indicates that there is a crossover in translational entropy of water molecules around hydrophobic solutes from negative to positive values as the size of the solute is increased. Rotational entropy of water molecules around hydrophobic solutes for all sizes increases with distance from the solute, indicating the absence of crossover in rotational entropy. This makes the crossover in total entropy (translation + rotation) of water molecule happen at much larger size (>1.5 nm) for hydrophobic solutes. Translational entropy of single water molecule scales logarithmically (Str(QH) = C + kB ln V), with the volume V obtained from the ellipsoid of inertia. We further discuss the origin of higher entropy of water around water and show the possibility of recovering the entropy loss of some hypothetical solutes. The results obtained are helpful to understand water entropy behavior around various hydrophobic and charged environments within biomolecules. Finally, we show how our approach can be used to calculate the entropy of the individual water molecules in a protein cavity that may be replaced during ligand binding.
Collapse
Affiliation(s)
- Wilbee D Sasikala
- Department of Chemistry, Indian Institute of Science Education and Research , Pune, India
| | | |
Collapse
|
27
|
Michel J, Henchman RH, Gerogiokas G, Southey MWY, Mazanetz MP, Law RJ. Evaluation of Host–Guest Binding Thermodynamics of Model Cavities with Grid Cell Theory. J Chem Theory Comput 2014; 10:4055-68. [DOI: 10.1021/ct500368p] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
Affiliation(s)
- Julien Michel
- EaStCHEM
School of Chemistry, Joseph Black Building, King’s Buildings, Edinburgh EH9 3JJ, United Kingdom
| | - Richard H. Henchman
- Manchester
Institute of Biotechnology, The University of Manchester, 131 Princess
Street, Manchester M1 7DN, United Kingdom
- School
of Chemistry, The University of Manchester, Oxford Road, Manchester M13 9PL, United Kingdom
| | - Georgios Gerogiokas
- EaStCHEM
School of Chemistry, Joseph Black Building, King’s Buildings, Edinburgh EH9 3JJ, United Kingdom
| | - Michelle W. Y. Southey
- Evotec (U.K.) Limited, Innovation
Drive 114 Milton Park, Abingdon, Oxfordshire, OX14 4RZ, United Kingdom
| | - Michael P. Mazanetz
- Evotec (U.K.) Limited, Innovation
Drive 114 Milton Park, Abingdon, Oxfordshire, OX14 4RZ, United Kingdom
| | - Richard J. Law
- Evotec (U.K.) Limited, Innovation
Drive 114 Milton Park, Abingdon, Oxfordshire, OX14 4RZ, United Kingdom
| |
Collapse
|
28
|
A quantum-mechanical study of the reaction mechanism of sulfite oxidase. J Biol Inorg Chem 2014; 19:1165-79. [PMID: 24957901 DOI: 10.1007/s00775-014-1172-z] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2013] [Accepted: 06/05/2014] [Indexed: 10/25/2022]
Abstract
The oxidation of sulfite to sulfate by two different models of the active site of sulfite oxidase has been studied. Both protonated and deprotonated substrates were tested. Geometries were optimized with density functional theory (TPSS/def2-SV(P)) and energies were calculated either with hybrid functionals and large basis sets (B3LYP/def2-TZVPD) including corrections for dispersion, solvation, and entropy, or with coupled-cluster theory (LCCSD(T0)) extrapolated toward a complete basis set. Three suggested reaction mechanisms have been compared and the results show that the lowest barriers are obtained for a mechanism where the substrate attacks a Mo-bound oxo ligand, directly forming a Mo-bound sulfate complex, which then dissociates into the products. Such a mechanism is more favorable than mechanisms involving a Mo-sulfite complex with the substrate coordinating either by the S or O atom. The activation energy is dominated by the Coulomb repulsion between the Mo complex and the substrate, which both have a negative charge of -1 or -2.
Collapse
|
29
|
Hansen N, van Gunsteren WF. Practical Aspects of Free-Energy Calculations: A Review. J Chem Theory Comput 2014; 10:2632-47. [PMID: 26586503 DOI: 10.1021/ct500161f] [Citation(s) in RCA: 289] [Impact Index Per Article: 28.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
Free-energy calculations in the framework of classical molecular dynamics simulations are nowadays used in a wide range of research areas including solvation thermodynamics, molecular recognition, and protein folding. The basic components of a free-energy calculation, that is, a suitable model Hamiltonian, a sampling protocol, and an estimator for the free energy, are independent of the specific application. However, the attention that one has to pay to these components depends considerably on the specific application. Here, we review six different areas of application and discuss the relative importance of the three main components to provide the reader with an organigram and to make nonexperts aware of the many pitfalls present in free energy calculations.
Collapse
Affiliation(s)
- Niels Hansen
- Institute of Thermodynamics and Thermal Process Engineering, University of Stuttgart , D-70569 Stuttgart, Germany.,Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH , CH-8093 Zürich, Switzerland
| | - Wilfred F van Gunsteren
- Laboratory of Physical Chemistry, Swiss Federal Institute of Technology, ETH , CH-8093 Zürich, Switzerland
| |
Collapse
|
30
|
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]
|
31
|
Gerogiokas G, Calabro G, Henchman RH, Southey MWY, Law RJ, Michel J. Prediction of Small Molecule Hydration Thermodynamics with Grid Cell Theory. J Chem Theory Comput 2013; 10:35-48. [PMID: 26579889 DOI: 10.1021/ct400783h] [Citation(s) in RCA: 50] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022]
Abstract
An efficient methodology has been developed to quantify water energetics by analysis of explicit solvent molecular simulations of organic and biomolecular systems. The approach, grid cell theory (GCT), relies on a discretization of the cell theory methodology on a three-dimensional grid to spatially resolve the density, enthalpy, and entropy of water molecules in the vicinity of solute(s) of interest. Entropies of hydration are found to converge more efficiently than enthalpies of hydration. GCT predictions of free energies of hydration on a data set of small molecules are strongly correlated with thermodynamic integration predictions. Agreement with the experiment is comparable for both approaches. A key advantage of GCT is its ability to provide from a single simulation insightful graphical analyses of spatially resolved components of the enthalpies and entropies of hydration.
Collapse
Affiliation(s)
- Georgios Gerogiokas
- EaStCHEM School of Chemistry , Joseph Black Building, The King's Buildings, Edinburgh EH9 3JJ, United Kingdom
| | - Gaetano Calabro
- EaStCHEM School of Chemistry , Joseph Black Building, The King's Buildings, Edinburgh EH9 3JJ, United Kingdom
| | - Richard H Henchman
- Manchester Institute of Biotechnology, The University of Manchester , 131 Princess Street, Manchester M1 7DN, United Kingdom and School of Chemistry, The University of Manchester , Oxford Road, Manchester M13 9PL, United Kingdom
| | - Michelle W Y Southey
- Evotec (U.K.) Ltd. , 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | - Richard J Law
- Evotec (U.K.) Ltd. , 114 Innovation Drive, Milton Park, Abingdon, Oxfordshire OX14 4SA, United Kingdom
| | - Julien Michel
- EaStCHEM School of Chemistry , Joseph Black Building, The King's Buildings, Edinburgh EH9 3JJ, United Kingdom
| |
Collapse
|
32
|
Ashworth IW, Hillier IH, Nelson DJ, Percy JM, Vincent MA. Olefin Metathesis by Grubbs–Hoveyda Complexes: Computational and Experimental Studies of the Mechanism and Substrate-Dependent Kinetics. ACS Catal 2013. [DOI: 10.1021/cs400164w] [Citation(s) in RCA: 51] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Affiliation(s)
- Ian W. Ashworth
- Global Research and Development, AstraZeneca, Silk Road Business Park, Charter Way,
Macclesfield SK10 2NA, U.K
| | - Ian H. Hillier
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13
9PL, U.K
| | - David J. Nelson
- WestCHEM, Department
of Pure and
Applied Chemistry, University of Strathclyde, 295 Cathedral Street, Glasgow G1 1XL, U.K
| | - Jonathan M. Percy
- WestCHEM, Department
of Pure and
Applied Chemistry, University of Strathclyde, 295 Cathedral Street, Glasgow G1 1XL, U.K
| | - Mark A. Vincent
- School of Chemistry, University of Manchester, Oxford Road, Manchester M13
9PL, U.K
| |
Collapse
|
33
|
Huber RG, Fuchs JE, von Grafenstein S, Laner M, Wallnoefer HG, Abdelkader N, Kroemer RT, Liedl KR. Entropy from state probabilities: hydration entropy of cations. J Phys Chem B 2013; 117:6466-72. [PMID: 23651109 PMCID: PMC3668459 DOI: 10.1021/jp311418q] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Entropy is an important energetic quantity determining the progression of chemical processes. We propose a new approach to obtain hydration entropy directly from probability density functions in state space. We demonstrate the validity of our approach for a series of cations in aqueous solution. Extensive validation of simulation results was performed. Our approach does not make prior assumptions about the shape of the potential energy landscape and is capable of calculating accurate hydration entropy values. Sampling times in the low nanosecond range are sufficient for the investigated ionic systems. Although the presented strategy is at the moment limited to systems for which a scalar order parameter can be derived, this is not a principal limitation of the method. The strategy presented is applicable to any chemical system where sufficient sampling of conformational space is accessible, for example, by computer simulations.
Collapse
Affiliation(s)
- Roland G Huber
- Department of Theoretical Chemistry, Faculty for Chemistry and Pharmacy, Center for Molecular Biosciences Innsbruck, Leopold-Franzens University Innsbruck, Innrain 80/82, A-6020 Innsbruck, Austria
| | | | | | | | | | | | | | | |
Collapse
|
34
|
Rosado LA, Vasconcelos IB, Palma MS, Frappier V, Najmanovich RJ, Santos DS, Basso LA. The mode of action of recombinant Mycobacterium tuberculosis shikimate kinase: kinetics and thermodynamics analyses. PLoS One 2013; 8:e61918. [PMID: 23671579 PMCID: PMC3646032 DOI: 10.1371/journal.pone.0061918] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2012] [Accepted: 03/14/2013] [Indexed: 12/03/2022] Open
Abstract
Tuberculosis remains as one of the main cause of mortality worldwide due to a single infectious agent, Mycobacterium tuberculosis. The aroK-encoded M. tuberculosis Shikimate Kinase (MtSK), shown to be essential for survival of bacilli, catalyzes the phosphoryl transfer from ATP to the carbon-3 hydroxyl group of shikimate (SKH), yielding shikimate-3-phosphate and ADP. Here we present purification to homogeneity, and oligomeric state determination of recombinant MtSK. Biochemical and biophysical data suggest that the chemical reaction catalyzed by monomeric MtSK follows a rapid-equilibrium random order of substrate binding, and ordered product release. Isothermal titration calorimetry (ITC) for binding of ligands to MtSK provided thermodynamic signatures of non-covalent interactions to each process. A comparison of steady-state kinetics parameters and equilibrium dissociation constant value determined by ITC showed that ATP binding does not increase the affinity of MtSK for SKH. We suggest that MtSK would more appropriately be described as an aroL-encoded type II shikimate kinase. Our manuscript also gives thermodynamic description of SKH binding to MtSK and data for the number of protons exchanged during this bimolecular interaction. The negative value for the change in constant pressure heat capacity (ΔCp) and molecular homology model building suggest a pronounced contribution of desolvation of non-polar groups upon binary complex formation. Thermodynamic parameters were deconvoluted into hydrophobic and vibrational contributions upon MtSK:SKH binary complex formation. Data for the number of protons exchanged during this bimolecular interaction are interpreted in light of a structural model to try to propose the likely amino acid side chains that are the proton donors to bulk solvent following MtSK:SKH complex formation.
Collapse
Affiliation(s)
- Leonardo Astolfi Rosado
- Centro de Pesquisas em Biologia Molecular e Funcional (CPBMF), Instituto Nacional de Ciência e Tecnologia em Tuberculose (INCT-TB), Pontifícia Universidade Católica do Rio Grande do Sul (PUCRS), Porto Alegre, RS, Brazil
- Programa de Pós-Graduação em Biologia Celular e Molecular, PUCRS, Porto Alegre, RS, Brazil
| | - Igor Bordin Vasconcelos
- Centro de Pesquisas em Biologia Molecular e Funcional (CPBMF), Instituto Nacional de Ciência e Tecnologia em Tuberculose (INCT-TB), Pontifícia Universidade Católica do Rio Grande do Sul (PUCRS), Porto Alegre, RS, Brazil
- Programa de Pós-Graduação em Biologia Celular e Molecular, PUCRS, Porto Alegre, RS, Brazil
| | - Mário Sérgio Palma
- Laboratório de Biologia Estrutural e Zooquímica, Centro de Estudos de Insetos Sociais, Departamento de Biologia, Instituto de Biociências de Rio Claro, Universidade Estadual Paulista (UNESP), Rio Claro, SP, Brazil
| | - Vincent Frappier
- Department of Biochemistry, Faculty of Medicine, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - Rafael Josef Najmanovich
- Department of Biochemistry, Faculty of Medicine, Université de Sherbrooke, Sherbrooke, Québec, Canada
| | - Diógenes Santiago Santos
- Centro de Pesquisas em Biologia Molecular e Funcional (CPBMF), Instituto Nacional de Ciência e Tecnologia em Tuberculose (INCT-TB), Pontifícia Universidade Católica do Rio Grande do Sul (PUCRS), Porto Alegre, RS, Brazil
- Programa de Pós-Graduação em Medicina e Ciências da Saúde, PUCRS, Porto Alegre, RS, Brazil
- Programa de Pós-Graduação em Biologia Celular e Molecular, PUCRS, Porto Alegre, RS, Brazil
| | - Luiz Augusto Basso
- Centro de Pesquisas em Biologia Molecular e Funcional (CPBMF), Instituto Nacional de Ciência e Tecnologia em Tuberculose (INCT-TB), Pontifícia Universidade Católica do Rio Grande do Sul (PUCRS), Porto Alegre, RS, Brazil
- Programa de Pós-Graduação em Medicina e Ciências da Saúde, PUCRS, Porto Alegre, RS, Brazil
- Programa de Pós-Graduação em Biologia Celular e Molecular, PUCRS, Porto Alegre, RS, Brazil
| |
Collapse
|
35
|
Affiliation(s)
- Riccardo Baron
- Department of Medicinal Chemistry, College of Pharmacy, and The Henry Eyring Center for Theoretical Chemistry, The University of Utah, Salt Lake City, Utah 84112-5820;
| | - J. Andrew McCammon
- Howard Hughes Medical Institute, Department of Chemistry and Biochemistry, Department of Pharmacology, and Center for Theoretical Biological Physics, University of California, San Diego, La Jolla, California 92093-0365;
| |
Collapse
|
36
|
Li JL, Mata RA, Ryde U. Large Density-Functional and Basis-Set Effects for the DMSO Reductase Catalyzed Oxo-Transfer Reaction. J Chem Theory Comput 2013; 9:1799-807. [DOI: 10.1021/ct301094r] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Affiliation(s)
- Ji-Lai Li
- Department of Theoretical Chemistry,
Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
- State Key Laboratory of Theoretical
and Computational Chemistry, Institute of Theoretical Chemistry, Jilin
University, Changchun 130023, People’s Republic of China
| | - Ricardo A. Mata
- Institut für Physikalische
Chemie, Georg-August-Universität Göttingen, Tammannstrasse
6, D-37077, Göttingen, Germany
| | - Ulf Ryde
- Department of Theoretical Chemistry,
Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| |
Collapse
|
37
|
Chen S, Ma J. Substituent effects on packing entropy and film morphologies in the nucleation of functionalized pentacenes on SiO2 substrate: Molecular dynamics simulations. J Chem Phys 2012; 137:074708. [DOI: 10.1063/1.4745899] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
38
|
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]
|
39
|
Purisima EO, Hogues H. Protein-ligand binding free energies from exhaustive docking. J Phys Chem B 2012; 116:6872-9. [PMID: 22432509 DOI: 10.1021/jp212646s] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We explore the use of exhaustive docking as an alternative to Monte Carlo and molecular dynamics sampling for the direct integration of the partition function for protein-ligand binding. We enumerate feasible poses for the ligand and calculate the Boltzmann factor contribution of each pose to the partition function. From the partition function, the free energy, enthalpy, and entropy can be derived. All our calculations are done with a continuum solvation model that includes solving the Poisson equation. In contrast to Monte Carlo and molecular dynamics simulations, exhaustive docking avoids (within the limitations of a discrete sampling) the question of "Have we run long enough?" due to its deterministic complete enumeration of states. We tested the method on the T4 lysozyme L99A mutant, which has a nonpolar cavity that can accommodate a number of small molecules. We tested two electrostatic models. Model 1 used a solute dielectric of 2.25 for the complex apoprotein and free ligand and 78.5 for the solvent. Model 2 used a solute dielectric of 2.25 for the complex and apoprotein but 1.0 for the free ligand. For our test set of eight molecules, we obtain a reasonable correlation with a Pearson r(2) = 0.66 using model 1. The trend in binding affinity ranking is generally preserved with a Kendall τ = 0.64 and Spearman ρ = 0.83. With model 2, the correlation is improved with a Pearson r(2) = 0.83, Kendall τ = 0.93, and Spearman ρ = 0.98. This suggests that the energy function and sampling method adequately captured most of the thermodynamics of binding of the nonpolar ligands to T4 lysozyme L99A.
Collapse
Affiliation(s)
- Enrico O Purisima
- Biotechnology Research Institute, National Research Council of Canada, 6100 Royalmount Avenue, Montreal, Quebec H4P 2R2, Canada.
| | | |
Collapse
|
40
|
Polyansky AA, Zubac R, Zagrovic B. Estimation of conformational entropy in protein-ligand interactions: a computational perspective. Methods Mol Biol 2012; 819:327-53. [PMID: 22183546 DOI: 10.1007/978-1-61779-465-0_21] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Conformational entropy is an important component of the change in free energy upon binding of a ligand to its target protein. As a consequence, development of computational techniques for reliable estimation of conformational entropies is currently receiving an increased level of attention in the context of computational drug design. Here, we review the most commonly used techniques for conformational entropy estimation from classical molecular dynamics simulations. Although by-and-large still not directly used in practical drug design, these techniques provide a golden standard for developing other, computationally less-demanding methods for such applications, in addition to furthering our understanding of protein-ligand interactions in general. In particular, we focus on the quasi-harmonic approximation and discuss different approaches that can be used to go beyond it, most notably, when it comes to treating anharmonic and/or correlated motions. In addition to reviewing basic theoretical formalisms, we provide a concrete set of steps required to successfully calculate conformational entropy from molecular dynamics simulations, as well as discuss a number of practical issues that may arise in such calculations.
Collapse
Affiliation(s)
- Anton A Polyansky
- Laboratory of Computational Biophysics, Department of Structural and Computational Biology, Max Perutz Laboratories, Vienna, Austria
| | | | | |
Collapse
|
41
|
Echeverria I, Amzel LM. Disaccharide binding to galectin-1: free energy calculations and molecular recognition mechanism. Biophys J 2011; 100:2283-92. [PMID: 21539798 DOI: 10.1016/j.bpj.2011.03.032] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2010] [Revised: 03/21/2011] [Accepted: 03/22/2011] [Indexed: 01/19/2023] Open
Abstract
Galectin-1, a member of the conserved family of carbohydrate-binding proteins with affinity for β-galactosides, is a key modulator of diverse cell functions such as immune response and regulation. The binding affinity and specificity of galectin-1 for eight different β-galactosyl terminal disaccharides was studied using molecular-dynamics simulations in which the ligand was pulled away from the binding site using a mechanical force. We present what we believe to be a novel procedure, based on combinations of multistep trajectories, that was used to estimate the binding free energy (ΔG) of each disaccharide. The computed binding free energy differences show excellent correlation with experimental values determined previously. The small differences in affinity among the disaccharides are the result of an exquisite balance between the strengths of the galectin-sugar H-bonds and the H-bonds the protein and the disaccharides make with the solvent. Analysis of the free energies along the reaction coordinate shows that disaccharide unbinding/binding presents no energetic barrier and, therefore, is diffusion-limited. In addition, the calculations revealed that as the ligand is undocked from the binding site, breaking of protein-disaccharide H-bonds takes place in stages with intermediate states in which the interactions are bridged by water molecules.
Collapse
Affiliation(s)
- Ignacia Echeverria
- Department of Biophysics and Biophysical Chemistry, The Johns Hopkins University, School of Medicine, Baltimore, Maryland, USA
| | | |
Collapse
|
42
|
Li X, Bachmanov AA, Maehashi K, Li W, Lim R, Brand JG, Beauchamp GK, Reed DR, Thai C, Floriano WB. Sweet taste receptor gene variation and aspartame taste in primates and other species. Chem Senses 2011; 36:453-75. [PMID: 21414996 DOI: 10.1093/chemse/bjq145] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023] Open
Abstract
Aspartame is a sweetener added to foods and beverages as a low-calorie sugar replacement. Unlike sugars, which are apparently perceived as sweet and desirable by a range of mammals, the ability to taste aspartame varies, with humans, apes, and Old World monkeys perceiving aspartame as sweet but not other primate species. To investigate whether the ability to perceive the sweetness of aspartame correlates with variations in the DNA sequence of the genes encoding sweet taste receptor proteins, T1R2 and T1R3, we sequenced these genes in 9 aspartame taster and nontaster primate species. We then compared these sequences with sequences of their orthologs in 4 other nontasters species. We identified 9 variant sites in the gene encoding T1R2 and 32 variant sites in the gene encoding T1R3 that distinguish aspartame tasters and nontasters. Molecular docking of aspartame to computer-generated models of the T1R2 + T1R3 receptor dimer suggests that species variation at a secondary, allosteric binding site in the T1R2 protein is the most likely origin of differences in perception of the sweetness of aspartame. These results identified a previously unknown site of aspartame interaction with the sweet receptor and suggest that the ability to taste aspartame might have developed during evolution to exploit a specialized food niche.
Collapse
Affiliation(s)
- Xia Li
- Monell Chemical Senses Center, Philadelphia, PA 19104, USA
| | | | | | | | | | | | | | | | | | | |
Collapse
|
43
|
Irudayam SJ, Henchman RH. Prediction and interpretation of the hydration entropies of monovalent cations and anions. Mol Phys 2011. [DOI: 10.1080/00268976.2010.532162] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
44
|
Ryde U, Mata RA, Grimme S. Does DFT-D estimate accurate energies for the binding of ligands to metal complexes? Dalton Trans 2011; 40:11176-83. [DOI: 10.1039/c1dt10867k] [Citation(s) in RCA: 73] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
45
|
Hayes JM, Skamnaki VT, Archontis G, Lamprakis C, Sarrou J, Bischler N, Skaltsounis AL, Zographos SE, Oikonomakos NG. Kinetics, in silico docking, molecular dynamics, and MM-GBSA binding studies on prototype indirubins, KT5720, and staurosporine as phosphorylase kinase ATP-binding site inhibitors: the role of water molecules examined. Proteins 2010; 79:703-19. [PMID: 21287607 DOI: 10.1002/prot.22890] [Citation(s) in RCA: 70] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2010] [Revised: 09/15/2010] [Accepted: 09/16/2010] [Indexed: 11/10/2022]
Abstract
With an aim toward glycogenolysis control in Type 2 diabetes, we have investigated via kinetic experiments and computation the potential of indirubin (IC₅₀ > 50 μM), indirubin-3'-oxime (IC₅₀ = 144 nM), KT5720 (K(i) = 18.4 nM) and staurosporine (K(i) = 0.37 nM) as phosphorylase kinase (PhKγtrnc) ATP-binding site inhibitors, with the latter two revealed as potent inhibitors in the low nM range. Because of lack of structural information, we have exploited information from homologous kinase complexes to direct in silico calculations (docking, molecular dynamics, and MMGBSA) to predict the binding characteristics of the four ligands. All inhibitors are predicted to bind in the same active site area as the ATP adenine ring, with binding dominated by hinge region hydrogen bonds to Asp104:O and Met106:O (all four ligands) and also Met106:NH (for the indirubins). The PhKγtrnc-staurosporine complex has the greatest number of receptor-ligand hydrogen bonds, while for the indirubin-3'-oxime and KT5720 complexes there is an important network of interchanging water molecules bridging inhibitor-enzyme contacts. The MM-GBSA results revealed the source of staurosporine's low nM potency to be favorable electrostatic interactions, while KT5720 has strong van der Waals contributions. KT5720 interacts with the greatest number of protein residues either by direct or 1-water bridged hydrogen bond interactions, and the potential for more selective PhK inhibition based on a KT5720 analogue has been established. Including receptor flexibility in Schrödinger induced-fit docking calculations in most cases correctly predicted the binding modes as compared with the molecular dynamics structures; the algorithm was less effective when there were key structural waters bridging receptor-ligand contacts.
Collapse
Affiliation(s)
- Joseph M Hayes
- Institute of Organic and Pharmaceutical Chemistry, National Hellenic Research Foundation, 48 Vassileos Constantinou Avenue, Athens 11635, Greece.
| | | | | | | | | | | | | | | | | |
Collapse
|
46
|
Henchman RH, Irudayam SJ. Topological hydrogen-bond definition to characterize the structure and dynamics of liquid water. J Phys Chem B 2010; 114:16792-810. [PMID: 21114302 DOI: 10.1021/jp105381s] [Citation(s) in RCA: 67] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
A definition that equates a hydrogen bond topologically with a local energy well in the potential energy surface is used to study the structure and dynamics of liquid water. We demonstrate the robustness of this hydrogen-bond definition versus the many other definitions which use fixed, arbitrary parameters, do not account for variable molecular environments, and cannot effectively resolve transition states. Our topology definition unambiguously shows that most water molecules are double acceptors but sizable proportions are single or triple acceptors. Almost all hydrogens are found to take part in hydrogen bonds. Broken hydrogen bonds only form when two molecules try to form two hydrogen bonds between them. The double acceptors have tetrahedral geometry, lower potential energy, entropy, and density, and slower dynamics. The single and triple acceptors have trigonal and trigonal bipyramidal geometry and when considered together have higher density, potential energy, and entropy, faster dynamics, and a tendency to cluster. These calculations use an extended theory for the entropy of liquid water that takes into account the variable number of hydrogen bonds. Hydrogen-bond switching is shown to depend explicitly on the variable number of hydrogen bonds accepted and the presence of interstitial water molecules. Transition state theory indicates that the switching of hydrogen bonds is a mildly activated process, requiring only a moderate distortion of hydrogen bonds. Three main types of switching events are observed depending on whether the donor and acceptor are already sharing a hydrogen bond. The switch may proceed with no intermediate or via a bifurcated-oxygen or cyclic dimer, both of which have a broken hydrogen bond and symmetric and asymmetric forms. Switching is found to be strongly coupled to whole-molecule vibration, particularly for the more mobile single and triple acceptors. Our analysis suggests that even though water is heterogeneous in terms of the number of hydrogen bonds, the coupling between neighbors on various length and time scales brings about greater continuity in its properties.
Collapse
Affiliation(s)
- Richard H Henchman
- Manchester Interdisciplinary Biocentre, The University of Manchester, 131 Princess Street, Manchester M1 7DN, United Kingdom.
| | | |
Collapse
|
47
|
Sawada T, Fedorov DG, Kitaura K. Role of the Key Mutation in the Selective Binding of Avian and Human Influenza Hemagglutinin to Sialosides Revealed by Quantum-Mechanical Calculations. J Am Chem Soc 2010; 132:16862-72. [DOI: 10.1021/ja105051e] [Citation(s) in RCA: 71] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Toshihiko Sawada
- Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, Core Research for Evolutional Science and Technology (CREST), Japan Science and Technology Agency (JST), 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan, and Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| | - Dmitri G. Fedorov
- Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, Core Research for Evolutional Science and Technology (CREST), Japan Science and Technology Agency (JST), 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan, and Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| | - Kazuo Kitaura
- Nanosystem Research Institute, National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba, Ibaraki 305-8568, Japan, Core Research for Evolutional Science and Technology (CREST), Japan Science and Technology Agency (JST), 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan, and Graduate School of Pharmaceutical Sciences, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan
| |
Collapse
|
48
|
Valiente PA, Gil A, Batista PR, Caffarena ER, Pons T, Pascutti PG. New parameterization approaches of the LIE method to improve free energy calculations of PlmII-Inhibitors complexes. J Comput Chem 2010; 31:2723-34. [PMID: 20839299 DOI: 10.1002/jcc.21566] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Pedro A Valiente
- Laboratorio de Biología Computacional y Diseño de Proteínas, Centro de Estudios de Proteínas, Facultad de Biología, Universidad de La Habana, Cuba.
| | | | | | | | | | | |
Collapse
|
49
|
Abstract
We use explicit solvent molecular dynamics simulations to estimate free energy, enthalpy, and entropy changes along the cavity-ligand association coordinate for a set of seven model systems with varying physicochemical properties. Owing to the simplicity of the considered systems we can directly investigate the role of water thermodynamics in molecular recognition. A broad range of thermodynamic signatures is found in which water (rather than cavity or ligand) enthalpic or entropic contributions appear to drive cavity-ligand binding or rejection. The unprecedented, nanoscale picture of hydration thermodynamics can help the interpretation and design of protein-ligand binding experiments. Our study opens appealing perspectives to tackle the challenge of solvent entropy estimation in complex systems and for improving molecular simulation models.
Collapse
|
50
|
Laurence C, Graton J, Berthelot M, Besseau F, Le Questel JY, Luçon M, Ouvrard C, Planchat A, Renault E. An Enthalpic Scale of Hydrogen-Bond Basicity. 4. Carbon π Bases, Oxygen Bases, and Miscellaneous Second-Row, Third-Row, and Fourth-Row Bases and a Survey of the 4-Fluorophenol Affinity Scale. J Org Chem 2010; 75:4105-23. [DOI: 10.1021/jo100461z] [Citation(s) in RCA: 75] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Christian Laurence
- Laboratoire CEISAM, UMR CNRS 6230, Faculté des Sciences et des Techniques, Université de Nantes, 2, rue de la Houssinière, BP 92208, 44322 Nantes, Cedex 3, France
| | - Jérôme Graton
- Laboratoire CEISAM, UMR CNRS 6230, Faculté des Sciences et des Techniques, Université de Nantes, 2, rue de la Houssinière, BP 92208, 44322 Nantes, Cedex 3, France
| | - Michel Berthelot
- Laboratoire CEISAM, UMR CNRS 6230, Faculté des Sciences et des Techniques, Université de Nantes, 2, rue de la Houssinière, BP 92208, 44322 Nantes, Cedex 3, France
| | - François Besseau
- Laboratoire CEISAM, UMR CNRS 6230, Faculté des Sciences et des Techniques, Université de Nantes, 2, rue de la Houssinière, BP 92208, 44322 Nantes, Cedex 3, France
| | - Jean-Yves Le Questel
- Laboratoire CEISAM, UMR CNRS 6230, Faculté des Sciences et des Techniques, Université de Nantes, 2, rue de la Houssinière, BP 92208, 44322 Nantes, Cedex 3, France
| | - Maryvonne Luçon
- Laboratoire CEISAM, UMR CNRS 6230, Faculté des Sciences et des Techniques, Université de Nantes, 2, rue de la Houssinière, BP 92208, 44322 Nantes, Cedex 3, France
| | - Carole Ouvrard
- Laboratoire CEISAM, UMR CNRS 6230, Faculté des Sciences et des Techniques, Université de Nantes, 2, rue de la Houssinière, BP 92208, 44322 Nantes, Cedex 3, France
| | - Aurélien Planchat
- Laboratoire CEISAM, UMR CNRS 6230, Faculté des Sciences et des Techniques, Université de Nantes, 2, rue de la Houssinière, BP 92208, 44322 Nantes, Cedex 3, France
| | - Eric Renault
- Laboratoire CEISAM, UMR CNRS 6230, Faculté des Sciences et des Techniques, Université de Nantes, 2, rue de la Houssinière, BP 92208, 44322 Nantes, Cedex 3, France
| |
Collapse
|