1
|
Kurapati R, Natarajan U. Complex role of chemical nature and tacticity in the adsorption free energy of carboxylic acid polymers at the oil-water interface: molecular dynamics simulations. Phys Chem Chem Phys 2023; 25:27783-27797. [PMID: 37814803 DOI: 10.1039/d3cp02754f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/11/2023]
Abstract
Scientific understanding of the molecular structure and adsorption of polymers at oil-water liquid interfaces is very limited. In this study the adsorption free energy at the oil (CCl4)-water interface was estimated using umbrella sampling molecular dynamics simulations for six carboxylate type vinyl polymers differing in hydrophobic nature and tacticity: isotactic and syndiotactic poly(acrylic acid) (i-PAA, s-PAA), isotactic and syndiotactic poly(methacrylic acid) (i-PMA, s-PMA), and atactic and syndiotactic poly(ethylacrylic acid) (a-PEA, s-PEA). ΔGads values are in the order i-PMA < a-PEA < s-PEA < s-PAA < i-PAA < s-PMA. The results show the significant and complex influence of the chemical nature as well as tacticity of the polymer on its adsorption free energy as related to hydrogen bonding and orientation of bonds with respect to oil and water phases. The influence of tacticity is found to be the highest for PMA, which is interpreted to occur due to the balance between interactions among side groups and those occurring between side groups and solvent. Interactions between side-groups are crucial for determining the conformation of PAA (most hydrophilic) and the solvation of the side-group in water is crucial for determining the conformation of PEA (most hydrophobic). The adsorption of PMA represents the transition between these two dominating effects. The molecular contributions to the enthalpy of adsorption indicate that adsorption is favored mainly through two interactions: polymer-CCl4 and water-water.
Collapse
Affiliation(s)
- Raviteja Kurapati
- Macromolecular Modeling and Simulation Laboratory, Department of Chemical Engineering, Indian Institute of Technology (IIT) Madras, Chennai, 600036, India.
| | - Upendra Natarajan
- Macromolecular Modeling and Simulation Laboratory, Department of Chemical Engineering, Indian Institute of Technology (IIT) Madras, Chennai, 600036, India.
| |
Collapse
|
2
|
Dongmo Foumthuim CJ, Giacometti A. Solvent quality and solvent polarity in polypeptides. Phys Chem Chem Phys 2023; 25:4839-4853. [PMID: 36692363 DOI: 10.1039/d2cp05214h] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
Using molecular dynamics and thermodynamic integration, we report on the solvation process of seven polypeptides (GLY, ALA, ILE, ASN, LYS, ARG, GLU) in water and in cyclohexane. The polypeptides are selected to cover the full hydrophobic scale while varying their chain length from tri- to undeca-homopeptides, providing indications on possible non-additivity effects as well as the role of the peptide backbone in the overall stability of the polypeptides. The use of different solvents and different polypeptides allows us to investigate the relation between solvent quality - the capacity of a given solvent to fold/unfold a given biopolymer often described on a scale ranging from "good" to "poor"; and solvent polarity - related to the specific interactions of any solvent with respect to a reference solvent. Undeca-glycine is found to be the only polypeptide to have a stable collapse in water (polar solvent), with the other hydrophobic polypeptides displaying repeated folding and unfolding events in water, with polar polypeptides presenting even more complex behavior. By contrast, all polypeptides are found to keep an extended conformation in cyclohexane, irrespective of their polarity. All considered polypeptides are also found to have favorable solvation free energy independent of the solvent polarity and their intrinsic hydrophobicity, clearly highlighting the prominent stabilizing role of the peptide backbone - with the solvation process largely enthalpically dominated in polar polypeptides and partially entropically driven for hydrophobic polypeptides. Our study thus reveals the complexity of the solvation process of polypeptides defying the common view "like dissolves like", with the solute polarity playing the most prominent role. The absence of mirror symmetry upon the inversion of polarities of both the solvent and the polypeptides is confirmed.
Collapse
Affiliation(s)
- Cedrix J Dongmo Foumthuim
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca' Foscari di Venezia, Campus Scientifico, Edificio Alfa, via Torino 155, 30172 Venezia Mestre, Italy.
| | - Achille Giacometti
- Dipartimento di Scienze Molecolari e Nanosistemi, Università Ca' Foscari di Venezia, Campus Scientifico, Edificio Alfa, via Torino 155, 30172 Venezia Mestre, Italy. .,European Centre for Living Technology (ECLT) Ca Bottacin, Dorsoduro 3911, Calle Crosera 30123 Venice, Italy
| |
Collapse
|
3
|
Yang H, Wang L, Yuan L, Du H, Pan B, Lu K. Antimicrobial Peptides with Rigid Linkers against Gram-Negative Bacteria by Targeting Lipopolysaccharide. JOURNAL OF AGRICULTURAL AND FOOD CHEMISTRY 2022; 70:15903-15916. [PMID: 36511360 DOI: 10.1021/acs.jafc.2c05921] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
A series of hybrid peptides were designed by connecting an antimicrobial peptide Ce(1-8) with a lipopolysaccharide (LPS)-targeting peptide Lf(28-34) via different linkers. Antimicrobial experimental results indicated that linkers play an essential role in the anti-Gram-negative bacterial activity of the hybrid peptides. Among these hybrid peptides, peptide CL5 with dipeptide rigid linker LP exhibited excellent activity and selectivity against Gram-negative bacteria. The minimum inhibitory concentrations of CL5 against the tested Gram-negative bacteria were 4-32 μM, while the toxicity toward HEK-293 cells was relatively low. It was found that the interactions of the peptides with LPS were crucial for peptide activity against Gram-negative bacteria. Antimicrobial mechanistic studies showed that peptide CL5 contributed to the death of Gram-negative bacterial cells by disrupting the integrity of the bacterial membranes. This study revealed the importance of linker selection in the design of hybrid peptides and provides the basis for the further development of antimicrobial peptides.
Collapse
Affiliation(s)
- Hongyan Yang
- College of Food Science and Engineering, Henan University of Technology, Zhengzhou 450001, China
| | - Lan Wang
- School of Chemical Engineering and Food Science, Zhengzhou University of Technology, Zhengzhou 450044, China
| | - Libo Yuan
- School of Chemistry and Chemical Engineering, Henan University of Technology, Zhengzhou 450001, China
| | - Heng Du
- College of Food Science and Engineering, Henan University of Technology, Zhengzhou 450001, China
| | - Boyuan Pan
- School of Chemistry and Chemical Engineering, Henan University of Technology, Zhengzhou 450001, China
| | - Kui Lu
- College of Food Science and Engineering, Henan University of Technology, Zhengzhou 450001, China
- School of Chemical Engineering and Food Science, Zhengzhou University of Technology, Zhengzhou 450044, China
| |
Collapse
|
4
|
Xiong D, Zhao X, Luo S, Zhang JZH, Duan L. Molecular Mechanism of the Non-Covalent Orally Targeted SARS-CoV-2 M pro Inhibitor S-217622 and Computational Assessment of Its Effectiveness against Mainstream Variants. J Phys Chem Lett 2022; 13:8893-8901. [PMID: 36126063 DOI: 10.1021/acs.jpclett.2c02428] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Convenient and efficient therapeutic agents are urgently needed to block the continued spread of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Here, the mechanism for the novel orally targeted SARS-CoV-2 main protease (Mpro) inhibitor S-217622 is revealed through a molecular dynamics simulation. The difference in the movement modes of the S-217622-Mpro complex and apo-Mpro suggested S-217622 could inhibit the motility intensity of Mpro, thus maintaining their stable binding. Subsequent energy calculations showed that the P2 pharmacophore possessed the highest energy contribution among the three pharmacophores of S-217622. Additionally, hot-spot residues H41, M165, C145, E166, and H163 have strong interactions with S-217622. To further investigate the resistance of S-217622 to six mainstream variants, the binding modes of S-217622 with these variants were elucidated. The subtle differences in energy compared to that of the wild type implied that the binding patterns of these systems were similar, and S-217622 still inhibited these variants. We hope this work will provide theoretical insights for optimizing novel targeted Mpro drugs.
Collapse
Affiliation(s)
- Danyang Xiong
- School of Physics and Electronics, Shandong Normal University, Jinan, Shandong 250014, China
| | - Xiaoyu Zhao
- School of Physics and Electronics, Shandong Normal University, Jinan, Shandong 250014, China
| | - Song Luo
- School of Physics and Electronics, Shandong Normal University, Jinan, Shandong 250014, China
| | - John Z H Zhang
- Shenzhen Institute of Synthetic Biology, Faculty of Synthetic Biology, Shenzhen Institute of Advanced Technology, Chinese Academy of Sciences, Shenzhen, Guangdong 518055, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai 200062, China
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai 200062, China
- Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China
| | - Lili Duan
- School of Physics and Electronics, Shandong Normal University, Jinan, Shandong 250014, China
| |
Collapse
|
5
|
Khan MT, Ali S, Zeb MT, Kaushik AC, Malik SI, Wei DQ. Gibbs Free Energy Calculation of Mutation in PncA and RpsA Associated With Pyrazinamide Resistance. Front Mol Biosci 2020; 7:52. [PMID: 32328498 PMCID: PMC7160322 DOI: 10.3389/fmolb.2020.00052] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2019] [Accepted: 03/16/2020] [Indexed: 12/16/2022] Open
Abstract
A central approach for better understanding the forces involved in maintaining protein structures is to investigate the protein folding and thermodynamic properties. The effect of the folding process is often disturbed in mutated states. To explore the dynamic properties behind mutations, molecular dynamic (MD) simulations have been widely performed, especially in unveiling the mechanism of drug failure behind mutation. When comparing wild type (WT) and mutants (MTs), the structural changes along with solvation free energy (SFE), and Gibbs free energy (GFE) are calculated after the MD simulation, to measure the effect of mutations on protein structure. Pyrazinamide (PZA) is one of the first-line drugs, effective against latent Mycobacterium tuberculosis isolates, affecting the global TB control program 2030. Resistance to this drug emerges due to mutations in pncA and rpsA genes, encoding pyrazinamidase (PZase) and ribosomal protein S1 (RpsA) respectively. The question of how the GFE may be a measure of PZase and RpsA stabilities, has been addressed in the current review. The GFE and SFE of MTs have been compared with WT, which were already found to be PZA-resistant. WT structures attained a more stable state in comparison with MTs. The physiological effect of a mutation in PZase and RpsA may be due to the difference in energies. This difference between WT and MTs, depicted through GFE plots, might be useful in predicting the stability and PZA-resistance behind mutation. This study provides useful information for better management of drug resistance, to control the global TB problem.
Collapse
Affiliation(s)
- Muhammad Tahir Khan
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| | - Sajid Ali
- Department of Microbiology, Quaid-i-Azam University Islamabad, Islamabad, Pakistan
| | | | - Aman Chandra Kaushik
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, Shanghai, China
| | - Shaukat Iqbal Malik
- Department of Bioinformatics and Biosciences, Capital University of Science and Technology, Islamabad, Pakistan
| | - Dong-Qing Wei
- State Key Laboratory of Microbial Metabolism, School of Life Sciences and Biotechnology, and Joint Laboratory of International Cooperation in Metabolic and Developmental Sciences, Ministry of Education, Shanghai Jiao Tong University, Shanghai, China
- Peng Cheng Laboratory, Shenzhen, China
| |
Collapse
|
6
|
Ansari N, Laio A, Hassanali A. Spontaneously Forming Dendritic Voids in Liquid Water Can Host Small Polymers. J Phys Chem Lett 2019; 10:5585-5591. [PMID: 31469575 DOI: 10.1021/acs.jpclett.9b02052] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Some liquids are characterized by the presence of large voids with dendritic shapes and for this reason are dubbed transiently porous. By using a battery of data analysis tools, we demonstrate that liquid water and methane are both characterized by transient porosity. We show that the thermodynamics of porosity is distinct from that associated with cavitation á la classical nucleation theory. The shapes of dendritic voids in both liquids with very different chemistries resemble those of small polymers. We further show, using free energy calculations, that the cost of solvating small hydrophobic polymers in water is consistent with the work associated with creating dendritic voids. The entropic and enthalpic contributions associated with hosting these polymers can thus be rationalized by the thermodynamics of fluctuations in bulk water.
Collapse
Affiliation(s)
- Narjes Ansari
- The Abdus Salam International Centre for Theoretical Physics , Strada Costiera 11 , 34151 Trieste , Italy
| | - Alessandro Laio
- The Abdus Salam International Centre for Theoretical Physics , Strada Costiera 11 , 34151 Trieste , Italy
- SISSA , Via Bonomea 265 , I-34136 Trieste , Italy
| | - Ali Hassanali
- The Abdus Salam International Centre for Theoretical Physics , Strada Costiera 11 , 34151 Trieste , Italy
| |
Collapse
|
7
|
Sumi T, Koga K. Theoretical analysis on thermodynamic stability of chignolin. Sci Rep 2019; 9:5186. [PMID: 30914684 PMCID: PMC6435801 DOI: 10.1038/s41598-019-41518-1] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2018] [Accepted: 03/11/2019] [Indexed: 11/28/2022] Open
Abstract
Understanding the dominant factor in thermodynamic stability of proteins remains an open challenge. Kauzmann's hydrophobic interaction hypothesis, which considers hydrophobic interactions between nonpolar groups as the dominant factor, has been widely accepted for about sixty years and attracted many scientists. The hypothesis, however, has not been verified or disproved because it is difficult, both theoretically and experimentally, to quantify the solvent effects on the free energy change in protein folding. Here, we developed a computational method for extracting the dominant factor behind thermodynamic stability of proteins and applied it to a small, designed protein, chignolin. The resulting free energy profile quantitatively agreed with the molecular dynamics simulations. Decomposition of the free energy profile indicated that intramolecular interactions predominantly stabilized collapsed conformations, whereas solvent-induced interactions, including hydrophobic ones, destabilized them. These results obtained for chignolin were consistent with the site-directed mutagenesis and calorimetry experiments for globular proteins with hydrophobic interior cores.
Collapse
Affiliation(s)
- Tomonari Sumi
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-Naka, Kita-ku, Okayama, 700-8530, Japan.
- Department of Chemistry, Faculty of Science, Okayama University, 3-1-1 Tsushima-Naka, Kita-ku, Okayama, 700-8530, Japan.
| | - Kenichiro Koga
- Research Institute for Interdisciplinary Science, Okayama University, 3-1-1 Tsushima-Naka, Kita-ku, Okayama, 700-8530, Japan
- Department of Chemistry, Faculty of Science, Okayama University, 3-1-1 Tsushima-Naka, Kita-ku, Okayama, 700-8530, Japan
| |
Collapse
|
8
|
Jiang Z, Remsing RC, Rego NB, Patel AJ. Characterizing Solvent Density Fluctuations in Dynamical Observation Volumes. J Phys Chem B 2019; 123:1650-1661. [PMID: 30682885 DOI: 10.1021/acs.jpcb.8b11423] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Hydrophobic effects drive diverse aqueous assemblies, such as micelle formation or protein folding, wherein the solvent plays an important role. Consequently, characterizing the free energetics of solvent density fluctuations can lead to important insights into these processes. Although techniques such as the indirect umbrella sampling (INDUS) method can be used to characterize solvent fluctuations in static observation volumes of various sizes and shapes, characterizing how the solvent mediates inherently dynamic processes, such as self-assembly or conformational change, remains a challenge. In this work, we generalize the INDUS method to facilitate the enhanced sampling of solvent fluctuations in dynamical observation volumes, whose positions and shapes can evolve. We illustrate the usefulness of this generalization by characterizing water density fluctuations in dynamical volumes pertaining to the hydration of flexible solutes, the assembly of small hydrophobes, and conformational transitions in a model peptide. We also use the method to probe the dynamics of hard spheres.
Collapse
Affiliation(s)
| | - Richard C Remsing
- Institute for Computational Molecular Science , Temple University , Philadelphia , Pennsylvania 19122 , United States
| | | | | |
Collapse
|
9
|
Seckfort D, Montgomery Pettitt B. Price of disorder in the lac repressor hinge helix. Biopolymers 2019; 110:e23239. [PMID: 30485404 PMCID: PMC6335174 DOI: 10.1002/bip.23239] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2018] [Revised: 09/12/2018] [Accepted: 10/04/2018] [Indexed: 12/26/2022]
Abstract
The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.
Collapse
Affiliation(s)
- Danielle Seckfort
- Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas
| | - B Montgomery Pettitt
- Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology, University of Texas Medical Branch, Galveston, Texas
| |
Collapse
|
10
|
Chaudhari MI, Rempe SB, Pratt LR. Quasi-chemical theory of F -(aq): The "no split occupancies rule" revisited. J Chem Phys 2018; 147:161728. [PMID: 29096480 DOI: 10.1063/1.4986244] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
We use ab initio molecular dynamics (AIMD) calculations and quasi-chemical theory (QCT) to study the inner-shell structure of F-(aq) and to evaluate that single-ion free energy under standard conditions. Following the "no split occupancies" rule, QCT calculations yield a free energy value of -101 kcal/mol under these conditions, in encouraging agreement with tabulated values (-111 kcal/mol). The AIMD calculations served only to guide the definition of an effective inner-shell constraint. QCT naturally includes quantum mechanical effects that can be concerning in more primitive calculations, including electronic polarizability and induction, electron density transfer, electron correlation, molecular/atomic cooperative interactions generally, molecular flexibility, and zero-point motion. No direct assessment of the contribution of dispersion contributions to the internal energies has been attempted here, however. We anticipate that other aqueous halide ions might be treated successfully with QCT, provided that the structure of the underlying statistical mechanical theory is absorbed, i.e., that the "no split occupancies" rule is recognized.
Collapse
Affiliation(s)
- Mangesh I Chaudhari
- Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Susan B Rempe
- Center for Biological and Engineering Sciences, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Lawrence R Pratt
- Department of Chemical and Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, USA
| |
Collapse
|
11
|
Izadi S, Harris RC, Fenley MO, Onufriev AV. Accuracy Comparison of Generalized Born Models in the Calculation of Electrostatic Binding Free Energies. J Chem Theory Comput 2018; 14:1656-1670. [PMID: 29378399 DOI: 10.1021/acs.jctc.7b00886] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
The need for accurate yet efficient representation of the aqueous environment in biomolecular modeling has led to the development of a variety of generalized Born (GB) implicit solvent models. While many studies have focused on the accuracy of available GB models in predicting solvation free energies, a systematic assessment of the quality of these models in binding free energy calculations, crucial for rational drug design, has not been undertaken. Here, we evaluate the accuracies of eight common GB flavors (GB-HCT, GB-OBC, GB-neck2, GBNSR6, GBSW, GBMV1, GBMV2, and GBMV3), available in major molecular dynamics packages, in predicting the electrostatic binding free energies ( ΔΔ Gel) for a diverse set of 60 biomolecular complexes belonging to four main classes: protein-protein, protein-drug, RNA-peptide, and small complexes. The GB flavors are examined in terms of their ability to reproduce the results from the Poisson-Boltzmann (PB) model, commonly used as accuracy reference in this context. We show that the agreement with the PB of ΔΔ Gel estimates varies widely between different GB models and also across different types of biomolecular complexes, with R2 correlations ranging from 0.3772 to 0.9986. A surface-based "R6" GB model recently implemented in AMBER shows the closest overall agreement with reference PB ( R2 = 0.9949, RMSD = 8.75 kcal/mol). The RNA-peptide and protein-drug complex sets appear to be most challenging for all but one model, as indicated by the large deviations from the PB in ΔΔ Gel. Small neutral complexes present the least challenge for most of the GB models tested. The quantitative demonstration of the strengths and weaknesses of the GB models across the diverse complex types provided here can be used as a guide for practical computations and future development efforts.
Collapse
Affiliation(s)
- Saeed Izadi
- Early Stage Pharmaceutical Development , Genentech Inc. , 1 DNA Way , South San Francisco , California 94080 , United States
| | - Robert C Harris
- Department of Pharmaceutical Sciences , University of Maryland School of Pharmacy , Baltimore , Maryland 21201 , United States
| | - Marcia O Fenley
- Institute of Molecular Biophysics , Florida State University , Tallahassee , Florida 32306-3408 , United States
| | | |
Collapse
|
12
|
Asthagiri D, Karandur D, Tomar DS, Pettitt BM. Intramolecular Interactions Overcome Hydration to Drive the Collapse Transition of Gly 15. J Phys Chem B 2017; 121:8078-8084. [PMID: 28774177 DOI: 10.1021/acs.jpcb.7b05469] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Simulations and experiments show oligo-glycines, polypeptides lacking any side chains, can collapse in water. We assess the hydration thermodynamics of this collapse by calculating the hydration free energy at each of the end points of the reaction coordinate, here taken as the end-to-end distance (r) in the chain. To examine the role of the various conformations for a given r, we study the conditional distribution, P(Rg|r), of the radius of gyration for a given value of r. The free energy change versus Rg, -kBT ln P(Rg|r), is found to vary more gently compared to the corresponding variation in the excess hydration free energy. Using this observation within a multistate generalization of the potential distribution theorem, we calculate a tight upper bound for the hydration free energy of the peptide for a given r. On this basis, we find that peptide hydration greatly favors the expanded state of the chain, despite primitive hydrophobic effects favoring chain collapse. The net free energy of collapse is seen to be a delicate balance between opposing intrapeptide and hydration effects, with intrapeptide contributions favoring collapse.
Collapse
Affiliation(s)
- D Asthagiri
- Department of Chemical and Biomolecular Engineering, Rice University , Houston, Texas, United States.,Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch , Galveston, Texas, United States
| | - Deepti Karandur
- Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine , Houston, Texas, United States
| | - Dheeraj S Tomar
- Chemical and Biomolecular Engineering, Johns Hopkins University , Baltimore, Maryland, United States
| | - B Montgomery Pettitt
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch , Galveston, Texas, United States.,Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine , Houston, Texas, United States
| |
Collapse
|
13
|
Drake JA, Harris RC, Pettitt BM. Solvation Thermodynamics of Oligoglycine with Respect to Chain Length and Flexibility. Biophys J 2017; 111:756-767. [PMID: 27558719 DOI: 10.1016/j.bpj.2016.07.013] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2016] [Revised: 06/30/2016] [Accepted: 07/06/2016] [Indexed: 01/24/2023] Open
Abstract
Oligoglycine is a backbone mimic for all proteins and is prevalent in the sequences of intrinsically disordered proteins. We have computed the absolute chemical potential of glycine oligomers at infinite dilution by simulation with the CHARMM36 and Amber ff12SB force fields. We performed a thermodynamic decomposition of the solvation free energy (ΔG(sol)) of Gly2-5 into enthalpic (ΔH(sol)) and entropic (ΔS(sol)) components as well as their van der Waals and electrostatic contributions. Gly2-5 was either constrained to a rigid/extended conformation or allowed to be completely flexible during simulations to assess the effects of flexibility on these thermodynamic quantities. For both rigid and flexible oligoglycine models, the decrease in ΔG(sol) with chain length is enthalpically driven with only weak entropic compensation. However, the apparent rates of decrease of ΔG(sol), ΔH(sol), ΔS(sol), and their elec and vdw components differ for the rigid and flexible models. Thus, we find solvation entropy does not drive aggregation for this system and may not explain the collapse of long oligoglycines. Additionally, both force fields yield very similar thermodynamic scaling relationships with respect to chain length despite both force fields generating different conformational ensembles of various oligoglycine chains.
Collapse
Affiliation(s)
- Justin A Drake
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas
| | - Robert C Harris
- Department of Chemistry, Temple University, Philadelphia, Pennsylvania
| | - B Montgomery Pettitt
- Department of Biochemistry and Molecular Biology, Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas.
| |
Collapse
|
14
|
Matubayasi N. Free-energy analysis of protein solvation with all-atom molecular dynamics simulation combined with a theory of solutions. Curr Opin Struct Biol 2016; 43:45-54. [PMID: 27835803 DOI: 10.1016/j.sbi.2016.10.005] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2016] [Accepted: 10/07/2016] [Indexed: 11/25/2022]
Abstract
The structure of a protein is strongly influenced by solvation. In the present review, the solvation effect is treated within the framework of statistical thermodynamics. The key quantity is the solvation free energy of protein and a fast method for its computation with explicit solvent is introduced. The applications of the method to the structural fluctuation of protein in water, structure determination of protein-protein complex, and urea effect on protein structure are then described, and the roles of solvent water and cosolvent are discussed from the standpoint of energetics.
Collapse
Affiliation(s)
- Nobuyuki Matubayasi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, Toyonaka, Osaka 560-8531, Japan; Elements Strategy Initiative for Catalysts and Batteries, Kyoto University, Katsura, Kyoto 615-8520, Japan.
| |
Collapse
|
15
|
Harris RC, Mackoy T, Fenley MO. Problems of robustness in Poisson-Boltzmann binding free energies. J Chem Theory Comput 2016; 11:705-12. [PMID: 26528091 PMCID: PMC4610304 DOI: 10.1021/ct5005017] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/09/2014] [Indexed: 11/29/2022]
Abstract
Although models based on the Poisson–Boltzmann (PB) equation have been fairly successful at predicting some experimental quantities, such as solvation free energies (ΔG), these models have not been consistently successful at predicting binding free energies (ΔΔG). Here we found that ranking a set of protein–protein complexes by the electrostatic component (ΔΔGel) of ΔΔG was more difficult than ranking the same molecules by the electrostatic component (ΔGel) of ΔG. This finding was unexpected because ΔΔGel can be calculated by combining estimates of ΔGel for the complex and its components with estimates of the ΔΔGel in vacuum. One might therefore expect that if a theory gave reliable estimates of ΔGel, then its estimates of ΔΔGel would be reliable. However, ΔΔGel for these complexes were orders of magnitude smaller than ΔGel, so although estimates of ΔGel obtained with different force fields and surface definitions were highly correlated, similar estimates of ΔΔGel were often not correlated.
Collapse
Affiliation(s)
- Robert C Harris
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, 301 University Boulevard, Galveston, Texas 77555-0304, United States
| | | | | |
Collapse
|
16
|
Harris RC, Pettitt BM. Reconciling the understanding of 'hydrophobicity' with physics-based models of proteins. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2016; 28:083003. [PMID: 26836518 PMCID: PMC5370576 DOI: 10.1088/0953-8984/28/8/083003] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
The idea that a 'hydrophobic energy' drives protein folding, aggregation, and binding by favoring the sequestration of bulky residues from water into the protein interior is widespread. The solvation free energies (ΔGsolv) of small nonpolar solutes increase with surface area (A), and the free energies of creating macroscopic cavities in water increase linearly with A. These observations seem to imply that there is a hydrophobic component (ΔGhyd) of ΔGsolv that increases linearly with A, and this assumption is widely used in implicit solvent models. However, some explicit-solvent molecular dynamics studies appear to contradict these ideas. For example, one definition (ΔG(LJ)) of ΔGhyd is that it is the free energy of turning on the Lennard-Jones (LJ) interactions between the solute and solvent. However, ΔG(LJ) decreases with A for alanine and glycine peptides. Here we argue that these apparent contradictions can be reconciled by defining ΔGhyd to be a near hard core insertion energy (ΔGrep), as in the partitioning proposed by Weeks, Chandler, and Andersen. However, recent results have shown that ΔGrep is not a simple function of geometric properties of the molecule, such as A and the molecular volume, and that the free energy of turning on the attractive part of the LJ potential cannot be computed from first-order perturbation theory for proteins. The theories that have been developed from these assumptions to predict ΔGhyd are therefore inadequate for proteins.
Collapse
Affiliation(s)
- Robert C Harris
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, 301 University Blvd, Galveston, TX 77555-0304, USA
| | | |
Collapse
|
17
|
Karandur D, Pettitt B. The contribution of electrostatic interactions to the collapse of oligoglycine in water. CONDENSED MATTER PHYSICS 2016; 19:23802. [PMID: 28943833 PMCID: PMC5607018 DOI: 10.5488/cmp.19.23802] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Protein solubility and conformational stability are a result of a balance of interactions both within a protein and between protein and solvent. The electrostatic solvation free energy of oligoglycines, models for the peptide backbone, becomes more favorable with an increasing length, yet longer peptides collapse due to the formation of favorable intrapeptide interactions between CO dipoles, in some cases without hydrogen bonds. The strongly repulsive solvent cavity formation is balanced by van der Waals attractions and electrostatic contributions. In order to investigate the competition between solvent exclusion and charge interactions we simulate the collapse of a long oligoglycine comprised of 15 residues while scaling the charges on the peptide from zero to fully charged. We examine the effect this has on the conformational properties of the peptide. We also describe the approximate thermodynamic changes that occur during the scaling both in terms of intrapeptide potentials and peptide-water potentials, and estimate the electrostatic solvation free energy of the system.
Collapse
Affiliation(s)
- D. Karandur
- Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, TX 77030, USA
| | - B.M. Pettitt
- Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, TX 77030, USA
- Sealy Center for Structural Biology and Molecular Biophysics, Department of Biochemistry and Molecular Biology, University of Texas Medical Branch, 301 University Blvd., Galveston, TX 77555-0304, USA
| |
Collapse
|
18
|
Tomar DS, Weber V, Pettitt BM, Asthagiri D. Importance of Hydrophilic Hydration and Intramolecular Interactions in the Thermodynamics of Helix-Coil Transition and Helix-Helix Assembly in a Deca-Alanine Peptide. J Phys Chem B 2015; 120:69-76. [PMID: 26649757 DOI: 10.1021/acs.jpcb.5b09881] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Abstract
For a model deca-alanine peptide the cavity (ideal hydrophobic) contribution to hydration favors the helix state over extended states and the paired helix bundle in the assembly of two helices. The energetic contributions of attractive protein-solvent interactions are separated into quasi-chemical components consisting of a short-range part arising from interactions with solvent in the first hydration shell and the remaining long-range part that is well described by a Gaussian. In the helix-coil transition, short-range attractive protein-solvent interactions outweigh hydrophobic hydration and favor the extended coil states. Analysis of enthalpic effects shows that it is the favorable hydration of the peptide backbone that favors the unfolded state. Protein intramolecular interactions favor the helix state and are decisive in favoring folding. In the pairing of two helices, the cavity contribution outweighs the short-range attractive protein-water interactions. However, long-range, protein-solvent attractive interactions can either enhance or reverse this trend depending on the mutual orientation of the helices. In helix-helix assembly, change in enthalpy arising from change in attractive protein-solvent interactions favors disassembly. In helix pairing as well, favorable protein intramolecular interactions are found to be as important as hydration effects. Overall, hydrophilic protein-solvent interactions and protein intramolecular interactions are found to play a significant role in the thermodynamics of folding and assembly in the system studied.
Collapse
Affiliation(s)
- Dheeraj S Tomar
- Department of Chemical and Biomolecular Engineering, Johns Hopkins University , Baltimore, Maryland 21218, United States
| | - Valéry Weber
- IBM Research, Zurich , CH-8803 Rüschlikon, Switzerland
| | - B Montgomery Pettitt
- Sealy Center for Structural Biology and Molecular Biophysics, Department of Biochemistry and Molecular Biology, University of Texas Medical Branch , Galveston, Texas 77555, United States
| | - D Asthagiri
- Sealy Center for Structural Biology and Molecular Biophysics, Department of Biochemistry and Molecular Biology, University of Texas Medical Branch , Galveston, Texas 77555, United States.,Department of Chemical and Biomolecular Engineering, Rice University , Houston, Texas 77005, United States
| |
Collapse
|
19
|
Cumberworth A, Bui JM, Gsponer J. Free energies of solvation in the context of protein folding: Implications for implicit and explicit solvent models. J Comput Chem 2015; 37:629-40. [DOI: 10.1002/jcc.24235] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2015] [Revised: 09/25/2015] [Accepted: 10/06/2015] [Indexed: 12/13/2022]
Affiliation(s)
| | | | - Jörg Gsponer
- Center for High-Throughput Biology, UBC; Vancouver Canada
| |
Collapse
|
20
|
Harris RC, Pettitt BM. Examining the assumptions underlying continuum-solvent models. J Chem Theory Comput 2015; 11:4593-600. [PMID: 26574250 DOI: 10.1021/acs.jctc.5b00684] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Continuum-solvent models (CSMs) have successfully predicted many quantities, including the solvation-free energies (ΔG) of small molecules, but they have not consistently succeeded at reproducing experimental binding free energies (ΔΔG), especially for protein-protein complexes. Several CSMs break ΔG into the free energy (ΔGvdw) of inserting an uncharged molecule into solution and the free energy (ΔGel) gained from charging. Some further divide ΔGvdw into the free energy (ΔGrep) of inserting a nearly hard cavity into solution and the free energy (ΔGatt) gained from turning on dispersive interactions between the solute and solvent. We show that for 9 protein-protein complexes neither ΔGrep nor ΔGvdw was linear in the solvent-accessible area A, as assumed in many CSMs, and the corresponding components of ΔΔG were not linear in changes in A. We show that linear response theory (LRT) yielded good estimates of ΔGatt and ΔΔGatt, but estimates of ΔΔGatt obtained from either the initial or final configurations of the solvent were not consistent with those from LRT. The LRT estimates of ΔGel differed by more than 100 kcal/mol from the explicit solvent model's (ESM's) predictions, and its estimates of the corresponding component (ΔΔGel) of ΔΔG differed by more than 10 kcal/mol. Finally, the Poisson-Boltzmann equation produced estimates of ΔGel that were correlated with those from the ESM, but its estimates of ΔΔGel were much less so. These findings may help explain why many CSMs have not been consistently successful at predicting ΔΔG for many complexes, including protein-protein complexes.
Collapse
Affiliation(s)
- Robert C Harris
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch , 301 University Blvd, Galveston, Texas 77555-0304, United States
| | - B Montgomery Pettitt
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch , 301 University Blvd, Galveston, Texas 77555-0304, United States
| |
Collapse
|
21
|
Karandur D, Harris RC, Pettitt BM. Protein collapse driven against solvation free energy without H-bonds. Protein Sci 2015; 25:103-10. [PMID: 26174309 DOI: 10.1002/pro.2749] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2015] [Accepted: 07/13/2015] [Indexed: 12/22/2022]
Abstract
Proteins collapse and fold because intramolecular interactions and solvent entropy, which favor collapse, outweigh solute-solvent interactions that favor expansion. Since the protein backbone actively participates in protein folding and some intrinsically disordered proteins are glycine rich, oligoglycines are good models to study the protein backbone as it collapses, both during conformational changes in disordered proteins and during folding. The solvation free energies of short glycine oligomers become increasingly favorable as chain length increases. In contrast, the solubility limits of glycine oligomers decrease with increasing chain length, indicating that peptide-peptide, and potentially solvent-solvent interactions, overcome peptide-solvent interactions to favor aggregation at finite concentrations of glycine oligomers. We have recently shown that hydrogen- (H-) bonds do not contribute significantly to the concentration-based aggregation of pentaglycines but that dipole-dipole (CO) interactions between the amide groups on the backbone do. Here we demonstrate for the collapse of oligoglycines ranging in length from 15 to 25 residues similarly that H-bonds do not contribute significantly to collapse but that CO dipole interactions do. These results illustrate that some intrapeptide interactions that determine the solubility limit of short glycine oligomers are similar to those that drive the collapse of longer glycine peptides.
Collapse
Affiliation(s)
- Deepti Karandur
- Howard Hughes Medical Institute at University of California, Berkeley, California, 94720-3220.,Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas, 77030
| | - Robert C Harris
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas, 77555-0304
| | - B Montgomery Pettitt
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, Galveston, Texas, 77555-0304.,Structural and Computational Biology and Molecular Biophysics, Baylor College of Medicine, Houston, Texas, 77030
| |
Collapse
|
22
|
Harris RC, Drake JA, Pettitt BM. Multibody correlations in the hydrophobic solvation of glycine peptides. J Chem Phys 2015; 141:22D525. [PMID: 25494796 DOI: 10.1063/1.4901886] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Protein collapse during folding is often assumed to be driven by a hydrophobic solvation energy (ΔGvdw) that scales linearly with solvent-accessible surface area (A). In a previous study, we argued that ΔGvdw, as well as its attractive (ΔGatt) and repulsive (ΔGrep) components, was not simply a linear function of A. We found that the surface tensions, γrep, γatt, and γvdw, gotten from ΔGrep, ΔGatt, and ΔGvdw against A for four configurations of deca-alanine differed from those obtained for a set of alkanes. In the present study, we extend our analysis to fifty decaglycine structures and atomic decompositions. We find that different configurations of decaglycine generate different estimates of γrep. Additionally, we considered the reconstruction of the solvation free energy from scaling the free energy of solvation of each atom type, free in solution. The free energy of the isolated atoms, scaled by the inverse surface area the atom would expose in the molecule does not reproduce the γrep for the intact decaglycines. Finally, γatt for the decaglycine conformations is much larger in magnitude than those for deca-alanine or the alkanes, leading to large negative values of γvdw (-74 and -56 cal/mol/Å(2) for CHARMM27 and AMBER ff12sb force fields, respectively). These findings imply that ΔGvdw favors extended rather than compact structures for decaglycine. We find that ΔGrep and ΔGvdw have complicated dependencies on multibody correlations between solute atoms, on the geometry of the molecular surface, and on the chemical identities of the atoms.
Collapse
Affiliation(s)
- Robert C Harris
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, 301 University Blvd, Galveston, Texas 77555-0304, USA
| | - Justin A Drake
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, 301 University Blvd, Galveston, Texas 77555-0304, USA
| | - B Montgomery Pettitt
- Sealy Center for Structural Biology and Molecular Biophysics, University of Texas Medical Branch, 301 University Blvd, Galveston, Texas 77555-0304, USA
| |
Collapse
|
23
|
Hajari T, van der Vegt NFA. Peptide backbone effect on hydration free energies of amino acid side chains. J Phys Chem B 2014; 118:13162-8. [PMID: 25338222 DOI: 10.1021/jp5094146] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
We have studied the hydrophobicity of amino acid side chains by computing conditional solvation free energies that account for effects of the peptide backbone on the side chains' solvent environment. The free energies reported herein correspond to a gas-liquid transfer process, which mimics solvation of the side chain under the condition that the backbone has been solvated already, and have been obtained on the basis of free energy calculations with empirical force field models. We find that the peptide backbone strongly impacts the solvation of nonpolar side chains, while its effect on the polar side chains is less pronounced. The results indicate that, in the presence of the short peptide backbone, nonpolar amino acid side chains are less hydrophobic than what is expected based on small molecule (analogue) solvation data.
Collapse
Affiliation(s)
- Timir Hajari
- Center of Smart Interfaces, Technische Universität Darmstadt , Alarich-Weiss-Straße 10, 64287, Darmstadt, Germany
| | | |
Collapse
|
24
|
Abstract
Inserting an uncharged van der Waals (vdw) cavity into water disrupts the distribution of water and creates attractive dispersion interactions between the solvent and solute. This free-energy change is the hydrophobic solvation energy (ΔG(vdw)). Frequently, it is assumed to be linear in the solvent-accessible surface area, with a positive surface tension (γ) that is independent of the properties of the molecule. However, we found that γ for a set of alkanes differed from that for four configurations of decaalanine, and γ = -5 was negative for the decaalanines. These findings conflict with the notion that ΔG(vdw) favors smaller A. We broke ΔG(vdw) into the free energy required to exclude water from the vdw cavity (ΔG(rep)) and the free energy of forming the attractive interactions between the solute and solvent (ΔG(att)) and found that γ < 0 for the decaalanines because -γ(att) > γ(rep) and γ(att) < 0. Additionally, γ(att) and γ(rep) for the alkanes differed from those for the decaalanines, implying that none of ΔG(att), ΔG(rep), and ΔG(vdw) can be computed with a constant surface tension. We also showed that ΔG(att) could not be computed from either the initial or final water distributions, implying that this quantity is more difficult to compute than is sometimes assumed. Finally, we showed that each atom's contribution to γ(rep) depended on multibody interactions with its surrounding atoms, implying that these contributions are not additive. These findings call into question some hydrophobic models.
Collapse
|