1
|
Hossain M, Huda N, Bhuyan AK. A three-state mechanism for trifluoroethanol denaturation of an intrinsically disordered protein (IDP). J Biochem 2023; 174:519-531. [PMID: 37709541 DOI: 10.1093/jb/mvad067] [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: 04/27/2023] [Revised: 09/05/2023] [Accepted: 09/13/2023] [Indexed: 09/16/2023] Open
Abstract
Relating the amino acid composition and sequence to chain folding and binding preferences of intrinsically disordered proteins (IDPs) has emerged as a huge challenge. While globular proteins have respective 3D structures that are unique to their individual functions, IDPs violate this structure-function paradigm because rather than having a well-defined structure an ensemble of rapidly interconverting disordered structures characterize an IDP. This work measures 2,2,2-trifluoroethanol (TFE)-induced equilibrium transitions of an IDP called AtPP16-1 (Arabidopsis thaliana phloem protein type 16-1) by using fluorescence, circular dichroism, infrared and nuclear magnetic resonance (NMR) methods at pH 4, 298 K. Low TFE reversibly removes the tertiary structure to produce an ensemble of obligate intermediate ($\mathrm{I}$) retaining the native-state ($\mathrm{N}$) secondary structure. The intermediate $\mathrm{I}$ is preceded by a non-obligate tryptophan-specific intermediate ${\mathrm{I}}_{\mathrm{w}}$ whose population is detectable for AtPP16-1 specifically. Accumulation of such non-obligate intermediates is discriminated according to the sequence composition of the protein. In all cases, however, a tertiary structure-unfolded general obligate intermediate $\mathrm{I}$ is indispensable. The $\mathrm{I}$ ensemble has higher helical propensity conducive to the acquisition of an exceedingly large level of α-helices by a reversible denaturation transition of $\mathrm{I}$ to the denatured state $\mathrm{D}$ as the TFE level is increased. Strikingly, it is the same $\mathrm{N}\rightleftharpoons \mathrm{I}\rightleftharpoons \mathrm{D}$ scheme typifying the TFE transitions of globular proteins. The high-energy state $\mathrm{I}$ characterized by increased helical propensity is called a universal intermediate encountered in both genera of globular and disordered proteins. Neither $\mathrm{I}$ nor $\mathrm{D}$ strictly show molten globule (MG)-like properties, dismissing the belief that TFE promotes MGs.
Collapse
Affiliation(s)
- Mujahid Hossain
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Noorul Huda
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Abani K Bhuyan
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| |
Collapse
|
2
|
Pérez-Trejo I, Dominguez L. GaMD simulations as an alternative in the TFE-water mixture description. J Mol Model 2023; 29:352. [PMID: 37906368 PMCID: PMC10618327 DOI: 10.1007/s00894-023-05749-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2023] [Accepted: 10/04/2023] [Indexed: 11/02/2023]
Abstract
CONTEXT 2,2,2-Trifluoroethanol has been widely used to study the structure and dynamic properties of intrinsically disordered proteins. Experimentally, it is known that TFE-water mixtures stabilize secondary structures of IDPs, and therefore, it allows the studying of conformational ensembles of these proteins. In the last decades, molecular dynamic simulations have helped study the IDPs' conformational ensemble. Unfortunately, conventional MD requires very long simulation times to describe the properties of IDPs. Therefore, a variety of accelerated sampling techniques have been developed and employed. The TFE-water mixture arrangement description through MD has faced substantial difficulties since emulating the TFE nanocrowding at certain TFE:H[Formula: see text]O ratios (around 15-40% of TFE). In this work, we determine the most suitable conditions that reproduce experimentally reported properties of TFE-water mixtures. We compared the employment of conventional MD and GaMD simulations and various water parameters. Our results show that the combination of parameters that better reproduce the experimental information is the combination of the TIP4PD water model and GaMD simulations. Therefore, these conditions help accurately describe the structural ensemble of IDPs in TFE-water mixtures. METHODS Conventional MD and GaMD simulations were performed under AMBER 18 software. The TFE and water molecules were described using GAFF2 and a variety of water models, such as TIP3P, TIP4P2005, TIP4PD, and TIP5P, respectively. The systems were simulated a 100 ns at 298 K.
Collapse
Affiliation(s)
- Itzel Pérez-Trejo
- Departamento de Fisicoquímica, Facultad de Química, Universidad Nacional Autónoma de México, Mexico City, 04510, Mexico
| | - Laura Dominguez
- Departamento de Fisicoquímica, Facultad de Química, Universidad Nacional Autónoma de México, Mexico City, 04510, Mexico.
| |
Collapse
|
3
|
Biswas S, Kaur S, Myers CA, Chen AA, Welch JT. Aggregation in Aqueous Solutions of 2-(Tetrafluoro(trifluoromethyl)-λ 6-sulfanyl-ethan-1-ol (CF 3SF 4-ethanol)): A Comparison with Aqueous Trifluoroethanol and Hexafluoroisopropanol Using Molecular Dynamics Simulations and Dynamic Light Scattering Experiments. ACS OMEGA 2023; 8:30037-30047. [PMID: 37636933 PMCID: PMC10448670 DOI: 10.1021/acsomega.3c02019] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 03/25/2023] [Accepted: 06/26/2023] [Indexed: 08/29/2023]
Abstract
2-Tetrafluoro(trifluoromethyl)-λ6-sulfanylethan-1-ol (CF3SF4-ethanol) combines the polar hydrophobicity of tetrafluoro(trifluoromethyl)-λ6-sulfanyl (CF3SF4) group with the polarity of simple alcohols. The properties of aqueous solutions of the well-known fluorinated alcohols 2,2,2-trifluoroethanol (TFE) and 1,1,1,3,3,3-hexafluoroisopropanol (HFIP) were compared with those of aqueous solutions of the novel CF3SF4-ethanol. Those properties were computed using all atom molecular dynamics simulations with OPLS-compatible parameters. DFT ab initio calculations were used to accurately describe the nonsymmetrical, hypervalent sulfur in CF3SF4-ethanol. Although the molecular and conformational characteristics of CF3SF4-ethanol are like those of both TFE and HFIP, the greater hydrophobicity and lower polarity of CF3SF4-ethanol resulted in solution phase aggregation at a much lower concentration. The properties computed for TFE and HFIP in this work were consistent with published computational and experimental studies. CF3SF4-ethanol is predicted to be environmentally benign and hence an excellent green solvent candidate while possessing many of the same properties as TFE or HFIP.
Collapse
Affiliation(s)
- Samadrita Biswas
- Department
of Chemistry, University at Albany, State
University of New York, 1400 Washington Ave, Albany, New York 12222, United States
| | - Simi Kaur
- Department
of Chemistry, University at Albany, State
University of New York, 1400 Washington Ave, Albany, New York 12222, United States
- RNA
Institute, University at Albany, State University
of New York, 1400 Washington Ave, Albany, New York 12222, United States
| | - Christopher A. Myers
- RNA
Institute, University at Albany, State University
of New York, 1400 Washington Ave, Albany, New York 12222, United States
- Department
of Physics, University at Albany, State
University of New York, 1400 Washington Ave, Albany, New York 12222, United States
| | - Alan A. Chen
- Department
of Chemistry, University at Albany, State
University of New York, 1400 Washington Ave, Albany, New York 12222, United States
- RNA
Institute, University at Albany, State University
of New York, 1400 Washington Ave, Albany, New York 12222, United States
| | - John T. Welch
- Department
of Chemistry, University at Albany, State
University of New York, 1400 Washington Ave, Albany, New York 12222, United States
| |
Collapse
|
4
|
Hossain M, Huda N, Bhuyan AK. A surprisingly simple three-state generic process for reversible protein denaturation by trifluoroethanol. Biophys Chem 2022; 291:106895. [PMID: 36182744 DOI: 10.1016/j.bpc.2022.106895] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2022] [Revised: 09/12/2022] [Accepted: 09/14/2022] [Indexed: 12/30/2022]
Abstract
Despite the rich knowledge of the influence of 2,2,2-trifluoroethanol (TFE) on the structure and conformation of peptides and proteins, the mode(s) of TFE-protein interactions and the mechanism by which TFE reversibly denatures a globular protein remain elusive. This study systematically examines TFE-induced equilibrium transition curves for six paradigmatic globular proteins by using basic fluorescence and circular dichroism measurements under neutral pH conditions. The results are remarkably simple. Low TFE invariably unfolds the tertiary structure of all proteins to produce the obligate intermediate (I) which retains nearly all of native-state secondary structure, but enables the formation of extra α-helices as the level of TFE is raised higher. Inspection of the transitions at once reveals that the tertiary structure unfolding is always a distinct process, necessitating the inclusion of at least one obligate intermediate in the TFE-induced protein denaturation. It appears that the intermediate in the minimal unfolding mechanism N⇌I⇌D somehow acquires higher α-helical propensity to generate α-helices in excess of that in the native state to produce the denatured state (D), also called the TFE state. The low TFE-populated intermediate I may be called a universal intermediate by virtue of its α-helical propensity. Contrary to many earlier suggestions, this study dismisses molten globule (MG)-like attribute of I or D.
Collapse
Affiliation(s)
- Mujahid Hossain
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Noorul Huda
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India
| | - Abani K Bhuyan
- School of Chemistry, University of Hyderabad, Hyderabad 500046, India.
| |
Collapse
|
5
|
Engelbrecht LDV, Ji X, Carbonaro CM, Laaksonen A, Mocci F. MD simulations explain the excess molar enthalpies in pseudo-binary mixtures of a choline chloride-based deep eutectic solvent with water or methanol. Front Chem 2022; 10:983281. [DOI: 10.3389/fchem.2022.983281] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2022] [Accepted: 10/28/2022] [Indexed: 11/16/2022] Open
Abstract
The addition of molecular liquid cosolvents to choline chloride (ChCl)-based deep eutectic solvents (DESs) is increasingly investigated for reducing the inherently high bulk viscosities of the latter, which represent a major obstacle for potential industrial applications. The molar enthalpy of mixing, often referred to as excess molar enthalpy HE—a property reflecting changes in intermolecular interactions upon mixing—of the well-known ChCl/ethylene glycol (1:2 molar ratio) DES mixed with either water or methanol was recently found to be of opposite sign at 308.15 K: Mixing of the DES with water is strongly exothermic, while methanol mixtures are endothermic over the entire mixture composition range. Knowledge of molecular-level liquid structural changes in the DES following cosolvent addition is expected to be important when selecting such “pseudo-binary” mixtures for specific applications, e.g., solvents. With the aim of understanding the reason for the different behavior of selected DES/water or methanol mixtures, we performed classical MD computer simulations to study the changes in intermolecular interactions thought to be responsible for the observed HE sign difference. Excess molar enthalpies computed from our simulations reproduce, for the first time, the experimental sign difference and composition dependence of the property. We performed a structural analysis of simulation configurations, revealing an intriguing difference in the interaction modes of the two cosolvents with the DES chloride anion: water molecules insert between neighboring chloride anions, forming ionic hydrogen-bonded bridges that draw the anions closer, whereas dilution of the DES with methanol results in increased interionic separation. Moreover, the simulated DES/water mixtures were found to contain extended hydrogen-bonded structures containing water-bridged chloride pair arrangements, the presence of which may have important implications for solvent applications.
Collapse
|
6
|
Kolocouris A, Arkin I, Glykos NM. A proof-of-concept study of the secondary structure of influenza A, B M2 and MERS- and SARS-CoV E transmembrane peptides using folding molecular dynamics simulations in a membrane mimetic solvent. Phys Chem Chem Phys 2022; 24:25391-25402. [PMID: 36239696 DOI: 10.1039/d2cp02881f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Here, we have carried out a proof-of-concept molecular dynamics (MD) simulation with adaptive tempering in a membrane mimetic environment to study the folding of single-pass membrane peptides. We tested the influenza A M2 viroporin, influenza B M2 viroporin, and protein E from coronaviruses MERS-Cov-2 and SARS-CoV-2 peptides with known experimental secondary structures in membrane bilayers. The two influenza-derived peptides are significantly different in the peptide sequence and secondary structure and more polar than the two coronavirus-derived peptides. Through a total of more than 50 μs of simulation time that could be accomplished in trifluoroethanol (TFE), as a membrane model, we characterized comparatively the folding behavior, helical stability, and helical propensity of these transmembrane peptides that match perfectly their experimental secondary structures, and we identified common motifs that reflect their quaternary organization and known (or not) biochemical function. We showed that BM2 is organized into two structurally distinct parts: a significantly more stable N-terminal half, and a fast-converting C-terminal half that continuously folds and unfolds between α-helical structures and non-canonical structures, which are mostly turns. In AM2, both the N-terminal half and C-terminal half are very flexible. In contrast, the two coronavirus-derived transmembrane peptides are much more stable and fast helix-formers when compared with the influenza ones. In particular, the SARS-derived peptide E appears to be the fastest and most stable helix-former of all the four viral peptides studied, with a helical structure that persists almost without disruption for the whole of its 10 μs simulation. By comparing the results with experimental observations, we benchmarked TFE in studying the conformation of membrane and hydrophobic peptides. This work provided accurate results suggesting a methodology to run long MD simulations and predict structural properties of biologically important membrane peptides.
Collapse
Affiliation(s)
- Antonios Kolocouris
- Laboratory of Medicinal Chemistry, Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Panepistimiopolis, 15771, Greece.
| | - Isaiah Arkin
- Department of Biological Chemistry, The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, Edmond J. Safra Campus Givat-Ram, Jerusalem, 91904, Israel
| | - Nicholas M Glykos
- Department of Molecular Biology and Genetics, Democritus University of Thrace, University Campus, Alexandroupolis, 68100, Greece.
| |
Collapse
|
7
|
Wu X, Yang S, Xu S, Zhang X, Ren Y. Measurement and correlation of the solubility of sodium acetate in eight pure and binary solvents. Chin J Chem Eng 2022. [DOI: 10.1016/j.cjche.2021.06.029] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
8
|
Kolář MH, Nagy G, Kunkel J, Vaiana SM, Bock LV, Grubmüller H. Folding of VemP into translation-arresting secondary structure is driven by the ribosome exit tunnel. Nucleic Acids Res 2022; 50:2258-2269. [PMID: 35150281 PMCID: PMC8887479 DOI: 10.1093/nar/gkac038] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 11/30/2021] [Accepted: 01/25/2022] [Indexed: 01/09/2023] Open
Abstract
The ribosome is a fundamental biomolecular complex that synthesizes proteins in cells. Nascent proteins emerge from the ribosome through a tunnel, where they may interact with the tunnel walls or small molecules such as antibiotics. These interactions can cause translational arrest with notable physiological consequences. Here, we studied the arrest caused by the regulatory peptide VemP, which is known to form α-helices inside the ribosome tunnel near the peptidyl transferase center under specific conditions. We used all-atom molecular dynamics simulations of the entire ribosome and circular dichroism spectroscopy to study the driving forces of helix formation and how VemP causes the translational arrest. To that aim, we compared VemP dynamics in the ribosome tunnel with its dynamics in solution. We show that the VemP peptide has a low helical propensity in water and that the propensity is higher in mixtures of water and trifluorethanol. We propose that helix formation within the ribosome is driven by the interactions of VemP with the tunnel and that a part of VemP acts as an anchor. This anchor might slow down VemP progression through the tunnel enabling α-helix formation, which causes the elongation arrest.
Collapse
Affiliation(s)
- Michal H Kolář
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, 370 77 Göttingen, Germany
- Department of Physical Chemistry, University of Chemistry and Technology in Prague, Technická 5, 166 28 Prague, Czech Republic
| | - Gabor Nagy
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, 370 77 Göttingen, Germany
| | - John Kunkel
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, AZ 85287, USA
| | - Sara M Vaiana
- Department of Physics and Center for Biological Physics, Arizona State University, Tempe, AZ 85287, USA
| | - Lars V Bock
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, 370 77 Göttingen, Germany
| | - Helmut Grubmüller
- Department of Theoretical and Computational Biophysics, Max Planck Institute for Multidisciplinary Sciences, Am Fassberg 11, 370 77 Göttingen, Germany
| |
Collapse
|
9
|
Lee M, Perry SL, Hayward RC. Complex Coacervation of Polymerized Ionic Liquids in Non-aqueous Solvents. ACS POLYMERS AU 2021; 1:100-106. [PMID: 36855425 PMCID: PMC9954202 DOI: 10.1021/acspolymersau.1c00017] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Oppositely charged polymerized ionic liquids (PILs) were used to form complex coacervates in two different organic solvents, 2,2,2-trifluoroethanol (TFE) and hexafluoro-2-propanol (HFIP), and the corresponding phase diagrams were constructed using UV-vis, NMR, and turbidity experiments. While previous studies on complex coacervates have focused almost exclusively on aqueous environments, the use of PILs in the current work enabled studies in solvents with substantially lower dielectric constants (27.0 for TFE, 16.7 for HFIP). The critical salt concentration required to induce complete miscibility was roughly 2-fold larger in HFIP compared with TFE, and two different PIL complexes, solidlike precipitates and liquidlike coacervates, were found in both systems. This study provides insight into the effects of low-dielectric-constant solvents on complex coacervation, which has not been widely studied because of the limited solubility of conventional polyelectrolytes in these media.
Collapse
Affiliation(s)
- Minjung Lee
- Department
of Polymer Science and Engineering, University
of Massachusetts, Amherst, 120 Governors Drive, Amherst, Massachusetts 01003-9263, United States
| | - Sarah L. Perry
- Department
of Chemical Engineering, University of Massachusetts,
Amherst, 686 North Pleasant
Street, Amherst, Massachusetts 01003-9303, United States
| | - Ryan C. Hayward
- Department
of Polymer Science and Engineering, University
of Massachusetts, Amherst, 120 Governors Drive, Amherst, Massachusetts 01003-9263, United States,Department
of Chemical and Biological Engineering, University of Colorado Boulder, 596
UCB, Boulder, Colorado 80309, United States,
| |
Collapse
|
10
|
Ploetz EA, Karunaweera S, Bentenitis N, Chen F, Dai S, Gee MB, Jiao Y, Kang M, Kariyawasam NL, Naleem N, Weerasinghe S, Smith PE. Kirkwood-Buff-Derived Force Field for Peptides and Proteins: Philosophy and Development of KBFF20. J Chem Theory Comput 2021; 17:2964-2990. [PMID: 33878263 DOI: 10.1021/acs.jctc.1c00075] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
A new classical nonpolarizable force field, KBFF20, for the simulation of peptides and proteins is presented. The force field relies heavily on the use of Kirkwood-Buff theory to provide a comparison of simulated and experimental Kirkwood-Buff integrals for solutes containing the functional groups common in proteins, thus ensuring intermolecular interactions that provide a good balance between the peptide-peptide, peptide-solvent, and solvent-solvent distributions observed in solution mixtures. In this way, it differs significantly from other biomolecular force fields. Further development and testing of the intermolecular potentials are presented here. Subsequently, rotational potentials for the ϕ/ψ and χ dihedral degrees of freedom are obtained by analysis of the Protein Data Bank, followed by small modifications to provide a reasonable balance between simulated and observed α and β percentages for small peptides. This, the first of two articles, describes in detail the philosophy and development behind KBFF20.
Collapse
Affiliation(s)
- Elizabeth A Ploetz
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Sadish Karunaweera
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nikolaos Bentenitis
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Feng Chen
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Shu Dai
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Moon B Gee
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Yuanfang Jiao
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Myungshim Kang
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nilusha L Kariyawasam
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | - Nawavi Naleem
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| | | | - Paul E Smith
- Department of Chemistry, Kansas State University, 213 CBC Building, 1212 Mid-Campus Drive North, Manhattan, Kansas 66506, United States
| |
Collapse
|
11
|
Dey A, Banik R, Ghosh S. Temperature Comparative Studies on
Self‐Assembly
of Sodium Dodecyl Sulphate and Didodecyl Dimethyl Ammonium Bromide in Aqueous, Brine, and Trifluoroethanol Media. J SURFACTANTS DETERG 2021. [DOI: 10.1002/jsde.12385] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Apensu Dey
- Centre for Surface Science, Physical Chemistry Section, Department of Chemistry Jadavpur University Kolkata 700032 India
| | - Rajesh Banik
- Centre for Surface Science, Physical Chemistry Section, Department of Chemistry Jadavpur University Kolkata 700032 India
| | - Soumen Ghosh
- Centre for Surface Science, Physical Chemistry Section, Department of Chemistry Jadavpur University Kolkata 700032 India
| |
Collapse
|
12
|
Czarnecka K, Wojasiński M, Ciach T, Sajkiewicz P. Solution Blow Spinning of Polycaprolactone-Rheological Determination of Spinnability and the Effect of Processing Conditions on Fiber Diameter and Alignment. MATERIALS (BASEL, SWITZERLAND) 2021; 14:1463. [PMID: 33802725 PMCID: PMC8002481 DOI: 10.3390/ma14061463] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 02/12/2021] [Revised: 03/03/2021] [Accepted: 03/04/2021] [Indexed: 11/17/2022]
Abstract
The growing popularity of solution blow spinning as a method for the production of fibrous tissue engineering scaffolds and the vast range of polymer-solvent systems available for the method raises the need to study the effect of processing conditions on fiber morphology and develop a method for its qualitative assessment. Rheological approaches to determine polymer solution spinnability and image analysis approaches to describe fiber diameter and alignment have been previously proposed, although in a separate manner and mostly for the widely known, well-researched electrospinning method. In this study, a series of methods is presented to determine the processing conditions for the development of submicron fibrous scaffolds. Rheological methods are completed with extensive image analysis to determine the spinnability window for a polymer-solvent system and qualitatively establish the influence of polymer solution concentration and collector rotational speed on fiber morphology, diameter, and alignment. Process parameter selection for a tissue engineering scaffold target application is discussed, considering the varying structural properties of the native extracellular matrix of the tissue of interest.
Collapse
Affiliation(s)
- Katarzyna Czarnecka
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Pawinskiego 5b, 02-106 Warsaw, Poland;
| | - Michał Wojasiński
- Faculty of Chemical and Process Engineering, Warsaw University of Technology, Warynskiego 1, 00-645 Warsaw, Poland; (M.W.); (T.C.)
| | - Tomasz Ciach
- Faculty of Chemical and Process Engineering, Warsaw University of Technology, Warynskiego 1, 00-645 Warsaw, Poland; (M.W.); (T.C.)
- Warsaw University of Technology, CEZAMAT, Poleczki 19, 02-822 Warsaw, Poland
| | - Pawel Sajkiewicz
- Institute of Fundamental Technological Research, Polish Academy of Sciences, Pawinskiego 5b, 02-106 Warsaw, Poland;
| |
Collapse
|
13
|
Stylianakis I, Shalev A, Scheiner S, Sigalas MP, Arkin IT, Glykos N, Kolocouris A. The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide. J Comput Chem 2020; 41:2177-2188. [PMID: 32735736 DOI: 10.1002/jcc.26381] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2020] [Revised: 06/17/2020] [Accepted: 06/18/2020] [Indexed: 11/07/2022]
Abstract
The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala25 ). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala25 . DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala25 . This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala25 . Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.
Collapse
Affiliation(s)
- Ioannis Stylianakis
- Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Athens, Greece
| | - Ariella Shalev
- Department of Biological Chemistry, The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, Edmond J. Safra Campus Givat-Ram, Jerusalem, Israel
| | - Steve Scheiner
- Department of Chemistry and Biochemistry, Utah State University, Logan, Utah, USA
| | - Michael P Sigalas
- Department of Chemistry, Laboratory of Applied Quantum Chemistry, Aristotle University of Thessaloniki, Thessaloniki, Greece
| | - Isaiah T Arkin
- Department of Biological Chemistry, The Alexander Silberman Institute of Life Sciences, The Hebrew University of Jerusalem, Edmond J. Safra Campus Givat-Ram, Jerusalem, Israel
| | - Nikolas Glykos
- Department of Molecular Biology and Genetics, Democritus University of Thrace, Alexandroupolis, Greece
| | - Antonios Kolocouris
- Section of Pharmaceutical Chemistry, Department of Pharmacy, National and Kapodistrian University of Athens, Athens, Greece
| |
Collapse
|
14
|
Nayar D. Small crowder interactions can drive hydrophobic polymer collapse as well as unfolding. Phys Chem Chem Phys 2020; 22:18091-18101. [PMID: 32760995 DOI: 10.1039/d0cp02402c] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Biomolecules evolve and function in the intracellular crowded environment that is densely packed with macromolecules. Yet, a microscopic understanding of the effects of such an environment on the conformational preferences of biomolecules remains elusive. While prior investigations have attributed crowding effects mainly to the excluded volume (size) effects of the crowders, very little is known about the effects exerted due to their chemical interactions. In this study, crowding effects of tri-alanine peptides on the collapse equilibria of generic hydrophobic polymer are investigated using molecular dynamics simulations. The role of weak, non-specific, attractive polymer-crowder interactions in modulating the polymer collapse equilibria is examined. The results highlight that crowding effects can lead to polymer compaction as well as unfolding depending on the strength of polymer-crowder interaction energy. Strongly interacting crowders weaken hydrophobic collapse (or unfold the polymer) at high volume fractions and induce polymer collapse only under dilute conditions. Weakly interacting crowders induce polymer collapse at all crowder concentrations. Interestingly, the thermodynamic driving forces for polymer collapse are remarkably different in the two cases. Strongly and weakly interacting crowders induce collapse by preferential adsorption and preferential depletion respectively. The findings provide new insights into the possible effects of interplay of intermolecular interactions in a crowded environment. The results have implications in understanding the impact of crowding in altering free energy landscapes of proteins.
Collapse
Affiliation(s)
- Divya Nayar
- Centre for Computational and Data Sciences, Indian Institute of Technology Kharagpur, West Bengal 721302, India.
| |
Collapse
|
15
|
Nair AS, Kumar S, Acharya S, Bagchi B. Rotation of small diatomics in water–ethanol mixture: Multiple breakdowns of hydrodynamic predictions. J Chem Phys 2020; 153:014504. [DOI: 10.1063/5.0005160] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Anjali S. Nair
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| | - Shubham Kumar
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| | - Subhajit Acharya
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| | - Biman Bagchi
- Solid State and Structural Chemistry Unit, Indian Institute of Science, Bangalore 560012, India
| |
Collapse
|
16
|
Pillai SA, Pal H, Ray D, Aswal VK, Bahadur P. Ethanol and 2, 2, 2- trifluoroethanol modified supramolecular assemblies of star block copolymers: A SANS and fluorescence study. J Mol Liq 2019. [DOI: 10.1016/j.molliq.2019.111117] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
17
|
Lupachev EV, Polkovnichenko AV, Kvashnin SY, Lotkhov VA, Kulov NN. Investigation of the Phase Equilibrium in Ternary Systems Formed by 2,2,2-Trifluoroethanol, Water, and Organic Alcohols. THEORETICAL FOUNDATIONS OF CHEMICAL ENGINEERING 2019. [DOI: 10.1134/s0040579519030060] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|
18
|
Blinov N, Wishart DS, Kovalenko A. Solvent Composition Effects on the Structural Properties of the Aβ42 Monomer from the 3D-RISM-KH Molecular Theory of Solvation. J Phys Chem B 2019; 123:2491-2506. [PMID: 30811210 DOI: 10.1021/acs.jpcb.9b00480] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Structural characterization of amyloid (A)β peptides implicated in Alzheimer's disease is a challenging problem due to their intrinsically disordered nature and their high propensity for aggregation. Only limited information is currently available from experiments on conformational properties and aggregation pathways of the peptides in cellular environments. In silico modeling complements experimental information, providing atomistic insight into structure and dynamics of different Aβ species. All-atom explicit solvent molecular dynamics (MD) simulations with a properly selected force field can deliver reliable structural and dynamic information. In the case of intrinsically disordered Aβ peptides, enhanced sampling simulations beyond the nanosecond time scale are required to obtain statistically meaningful results even for simple solvent conditions. To overcome the challenges of conformational sampling in crowded cellular environments, alternative approaches have to be used, including postprocessing of MD data. In this study, we employ the statistical-mechanical, three-dimensional reference interaction site model with the Kovalenko-Hirata closure integral equation molecular theory of solvation to describe solvent composition effects on the conformational equilibrium in a structural ensemble of the Aβ42 (covering residues 1-42) monomer based on a statistical reweighting technique. The methodology enables a computationally efficient prediction on how different factors in the cellular environment, such as solvent composition, nonpolar solvation, and macromolecular crowding, affect the structural properties of the monomer. Similarities have been identified between changes in the structural ensemble caused by nonpolar solvation and crowded environments modeled by ionic solution with large negative ions. In particular, both solvent conditions reduce the random coil content and enhance the helical structure content of the monomer. In contrast to the previous studies, which reported increased α-helical content of peptides in crowded environments, this work attributes these structural features to the difference in solvent exposure of hydrophilic residues of the monomer for different secondary structure elements, rather than to (entropic) excluded volume effects.
Collapse
Affiliation(s)
- Nikolay Blinov
- Department of Mechanical Engineering , Edmonton , Alberta T6G 1H9 , Canada.,Nanotechnology Research Centre , Edmonton , Alberta T6G 2M9 , Canada
| | - David S Wishart
- Departments of Computing Science and Biological Sciences , University of Alberta , Edmonton , Alberta T6G 2E8 , Canada
| | - Andriy Kovalenko
- Department of Mechanical Engineering , Edmonton , Alberta T6G 1H9 , Canada.,Nanotechnology Research Centre , Edmonton , Alberta T6G 2M9 , Canada
| |
Collapse
|
19
|
Hunt JR, Tseng C, Dawlaty JM. Donor-acceptor preassociation, excited state solvation threshold, and optical energy cost as challenges in chemical applications of photobases. Faraday Discuss 2019; 216:252-268. [PMID: 31025987 DOI: 10.1039/c8fd00215k] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Photobases are molecules with increased pKa in the excited state that can serve to transduce light energy into proton removal capability. They can be used to control chemical reactions using light, such as removing protons from a catalytic site in reactions that are rate-limited by proton transfer. We identify and explore several major challenges toward their practical applications. Two important challenges are the need for pre-association (or ground state hydrogen bonding) between the proton donor and the photobase, and the need for excited state solvation of the photogenerated products. We investigate these two challenges with the photobase 5-methoxyquinoline as the proton acceptor and a low-pKa alcohol, 2,2,2-trifluoroethanol, as the proton donor. We vary the concentration of the donor in a background non-hydrogen-bonding solvent. Using absorption spectroscopy, we have identified that the donor-acceptor concentration ratio must exceed 100 : 1 to achieve appreciable ground state hydrogen bonding. Interestingly, emission spectroscopy reveals that the onset of ground state hydrogen bonding does not guarantee successful excited state proton transfer. It takes an additional order of magnitude increase in donor-acceptor ratio to achieve that goal, revealing that it is necessary to have excess donor molecules to reach the solvation threshold for the photogenerated products. The next challenge is reducing the large ground-excited state energy gap, which often requires UV photons to drive proton transfer. We show experimental and computational data comparing the photobasicity and optical energy gap for a few N-aromatic heterocyclic photobases. In general, we find that reducing the energy gap by increasing the conjugation size necessarily reduces photobasicity, while adding substituents of varying electron-withdrawing strength allows some fine-tuning of this effect. The combination of these two factors provide a preliminary design space for creating new photobasic molecules.
Collapse
Affiliation(s)
| | - Cindy Tseng
- University of Southern California, 920 Bloom Walk, Los Angeles, USA.
| | - Jahan M Dawlaty
- University of Southern California, 920 Bloom Walk, Los Angeles, USA.
| |
Collapse
|
20
|
Zeindlhofer V, Berger M, Steinhauser O, Schröder C. A shell-resolved analysis of preferential solvation of coffee ingredients in aqueous mixtures of the ionic liquid 1-ethyl-3-methylimidazolium acetate. J Chem Phys 2018; 148:193819. [PMID: 30307218 DOI: 10.1063/1.5009802] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Abstract
Ionic liquids increase the solubility of various coffee ingredients in aqueous solution but little is known about the underlying mechanism. Kirkwood-Buff integrals as well as the potential of mean force indicate that the imidazolium cations are accumulated at the surface of the solutes, removing water molecules from the solute surface. Although hydrogen bonding of the anions to hydroxy groups of the solutes can be detected, their concentration at the surface is less enhanced compared to the cations. The decomposition into solvation shells by Voronoi tessellation reveals that structural features are only observed in the first solvation shell. Nevertheless, the depletion of water and the excess concentration of the ions and, in particular, of the cations are visible in the next solvation shells as well. Therefore, classical arguments of hydrotropic theory fail to explain this behavior.
Collapse
Affiliation(s)
- Veronika Zeindlhofer
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Vienna A-1090, Austria
| | - Magdalena Berger
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Vienna A-1090, Austria
| | - Othmar Steinhauser
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Vienna A-1090, Austria
| | - Christian Schröder
- Faculty of Chemistry, Department of Computational Biological Chemistry, University of Vienna, Vienna A-1090, Austria
| |
Collapse
|
21
|
Abstract
Methanol is the simplest alcohol and possible energy carrier because it is easier to store than hydrogen and burns cleaner than fossil fuels. It is a colorless liquid, completely miscible with water and organic solvents and is very hygroscopic. Here, simple two-dimensional models of methanol, based on Mercedes-Benz (MB) model of water, are examined by Monte Carlo simulations. Methanol particles are modeled as dimers formed by an apolar Lennard-Jones disk, mimicking the methyl group, and a sphere with two hydrogen bonding arms for the hydroxyl group. The used models are the one proposed by Hribar-Lee and Dill (Acta Chimica Slovenica, 53:257, 2006.) with the overlapping discs and a new model with tangentially fused dimers. The comparison was done between the models, in connection to the MB water, as well as with experimental results and with new simulations done for 3D models of methanol. Both 2D models show similar trends in structuring and thermodynamics. The difference is the most pronounced at lower temperatures, where the smaller model exhibits spontaneous crystallization, while the larger model shows metastable states. The 2D structural organization represents well the clustering tendency observed in 3D models, as well as in experiments. The models qualitatively agree with the bulk methanol thermodynamic properties like density and isothermal compressibility, however, heat capacity at the constant pressure shows trend more similar to the water behavior. This work on the smallest amphiphilic organic solute provides a simple testing ground to study the competition between polar and non-polar effects within the molecule and physical properties.
Collapse
Affiliation(s)
- Tomislav Primorac
- Faculty of Science, University of Split, Rudjera Boškovića 33, 21000 Split, Croatia
- Fakultät für Maschinenbau, Universität Paderborn, Warburger Str. 100, 33098 Paderborn, Germany
| | - Martina Požar
- Faculty of Science, University of Split, Rudjera Boškovića 33, 21000 Split, Croatia
- Laboratoire de Physique Théorique de la Matière Condensée (UMR CNRS 7600), Université Pierre et Marie Curie, 4 Place Jussieu, F75252, Paris cedex 05, France
| | - Franjo Sokolić
- Faculty of Science, University of Split, Rudjera Boškovića 33, 21000 Split, Croatia
| | - Larisa Zoranić
- Faculty of Science, University of Split, Rudjera Boškovića 33, 21000 Split, Croatia
| | - Tomaz Urbic
- Faculty of Chemistry and Chemical Technology, University of Ljubljana, Večna pot 113, 1000 Ljubljana, Slovenia
| |
Collapse
|
22
|
Mondal S, Biswas B, Nandy T, Singh PC. Understanding the Role of Hydrophobic Terminal in the Hydrogen Bond Network of the Aqueous Mixture of 2,2,2-Trifluoroethanol: IR, Molecular Dynamics, Quantum Chemical as Well as Atoms in Molecules Studies. J Phys Chem B 2018; 122:6616-6626. [DOI: 10.1021/acs.jpcb.8b04365] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Saptarsi Mondal
- Department of Spectroscopy, Indian Association for the Cultivation of Science, Kolkata 700032, India
| | - Biswajit Biswas
- Department of Spectroscopy, Indian Association for the Cultivation of Science, Kolkata 700032, India
| | - Tonima Nandy
- Department of Spectroscopy, Indian Association for the Cultivation of Science, Kolkata 700032, India
| | - Prashant Chandra Singh
- Department of Spectroscopy, Indian Association for the Cultivation of Science, Kolkata 700032, India
| |
Collapse
|
23
|
TFE-induced local unfolding and fibrillation of SOD1: bridging the experiment and simulation studies. Biochem J 2018; 475:1701-1719. [DOI: 10.1042/bcj20180085] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Revised: 04/05/2018] [Accepted: 04/23/2018] [Indexed: 01/03/2023]
Abstract
Misfolding and aggregation of Cu, Zn Superoxide dismutase (SOD1) is involved in the neurodegenerative disease, amyotrophic lateral sclerosis. Many studies have shown that metal-depleted, monomeric form of SOD1 displays substantial local unfolding dynamics and is the precursor for aggregation. Here, we have studied the structure and dynamics of different apo monomeric SOD1 variants associated with unfolding and aggregation in aqueous trifluoroethanol (TFE) through experiments and simulation. TFE induces partially unfolded β-sheet-rich extended conformations in these SOD1 variants, which subsequently develops aggregates with fibril-like characteristics. Fibrillation was achieved more easily in disulfide-reduced monomeric SOD1 when compared with wild-type and mutant monomeric SOD1. At higher concentrations of TFE, a native-like structure with the increase in α-helical content was observed. The molecular dynamics simulation results illustrate distinct structural dynamics for different regions of SOD1 variants and show uniform local unfolding of β-strands. The strands protected by the zinc-binding and electrostatic loops were found to unfold first in 20% (v/v) TFE, leading to a partial unfolding of β-strands 4, 5, and 6 which are prone to aggregation. Our results thus shed light on the role of local unfolding and conformational dynamics in SOD1 misfolding and aggregation.
Collapse
|
24
|
Kumar V, Prakash A, Lynn AM. Alterations in local stability and dynamics of A4V SOD1 in the presence of trifluoroethanol. Biopolymers 2018; 109:e23102. [DOI: 10.1002/bip.23102] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2017] [Revised: 12/26/2017] [Accepted: 01/03/2018] [Indexed: 12/13/2022]
Affiliation(s)
- Vijay Kumar
- Centre for Interdisciplinary Research in Basic Sciences, Jamia Millia Islamia, Jamia Nagar; New Delhi 110025 India
| | - Amresh Prakash
- School of Computational and Integrative Sciences; Jawaharlal Nehru University; New Delhi 110067 India
| | - Andrew M. Lynn
- School of Computational and Integrative Sciences; Jawaharlal Nehru University; New Delhi 110067 India
| |
Collapse
|
25
|
Adamidou T, Arvaniti KO, Glykos NM. Folding Simulations of a Nuclear Receptor Box-Containing Peptide Demonstrate the Structural Persistence of the LxxLL Motif Even in the Absence of Its Cognate Receptor. J Phys Chem B 2017; 122:106-116. [DOI: 10.1021/acs.jpcb.7b10292] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Triantafyllia Adamidou
- Department of Molecular Biology
and Genetics, Democritus University of Thrace, University campus, 68100 Alexandroupolis, Greece
| | - Konstantina-Olympia Arvaniti
- Department of Molecular Biology
and Genetics, Democritus University of Thrace, University campus, 68100 Alexandroupolis, Greece
| | - Nicholas M. Glykos
- Department of Molecular Biology
and Genetics, Democritus University of Thrace, University campus, 68100 Alexandroupolis, Greece
| |
Collapse
|
26
|
Požar M, Perera A. Evolution of the micro-structure of aqueous alcohol mixtures with cooling: A computer simulation study. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.10.039] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
27
|
Floris FM. Excess Volumes from the Pressure Derivative of the Excess Chemical Potential: Testing Simple Models for Cavity Formation in Water. ACS OMEGA 2017; 2:6424-6436. [PMID: 31457245 PMCID: PMC6644935 DOI: 10.1021/acsomega.7b01157] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/09/2017] [Accepted: 09/22/2017] [Indexed: 06/10/2023]
Abstract
Excess volumes and excess compressibilities for hard spheres in water were computed by pressure derivatives of the excess chemical potential, which is equivalent to the work of cavity formation. This is relevant to the application of continuum solvation methods at various pressures. The excess chemical potential was modeled within phenomenological expressions for curved surfaces plus a pressure-volume term, for which two approaches were adopted, differing for the radius of the spherical volume. This implies a different dependence on pressure of parameters. In all cases, in the surface term, for the pressure derivative of parameters of the curvature function, use was made of the previously proposed expressions for the first two moments obtained from the density and radial distribution of oxygens in liquid water. Only for the parameter which has the dimension of surface tension (γ̃) was explicit dependence on pressure considered and results are affected by the specific polynomial used. In agreement with what inferred from simulation results obtained for cavities in TIP4P water, negative and positive adsorptions at the contact radius were extrapolated for a very large cavity at 1 and 8000 atm, respectively. The expressions here employed for the excess chemical potential predict the zero value of asymptotic adsorption to be at a pressure between 500 and 800 atm, which can be compared to results from the revised scaled particle theory. In the same range, for a nanometer-sized cavity, a change of behavior occurs regarding the ratio between the excess Helmholtz free energy and the product between pressure and excess volume.
Collapse
|
28
|
Takekiyo T, Yoshida K, Funahashi Y, Nagata S, Abe H, Yamaguchi T, Yoshimura Y. Helix-forming ability of proteins in alkylammonium nitrate. J Mol Liq 2017. [DOI: 10.1016/j.molliq.2017.08.066] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
29
|
Prakash A, Kumar V, Pandey P, Bharti DR, Vishwakarma P, Singh R, Hassan MI, Lynn AM. Solvent sensitivity of protein aggregation in Cu, Zn superoxide dismutase: a molecular dynamics simulation study. J Biomol Struct Dyn 2017; 36:2605-2617. [PMID: 28782426 DOI: 10.1080/07391102.2017.1364670] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
Abstract
Misfolding and aggregation of Cu, Zn Superoxide Dismutase (SOD1) is often found in amyotrophic lateral sclerosis (ALS) patients. The central apo SOD1 barrel was involved in protein maturation and pathological aggregation in ALS. In this work, we employed atomistic molecular dynamics (MD) simulations to study the conformational dynamics of SOD1barrel monomer in different concentrations of trifluoroethanol (TFE). We find concentration dependence unusual structural and dynamical features, characterized by the local unfolding of SOD1barrel. This partially unfolded structure is characterized by the exposure of hydrophobic core, is highly dynamic in nature, and is the precursor of aggregation seen in SOD1barrel. Our computational studies supports the hypothesis of the formation of aggregation 'building blocks' by means of local unfolding of apo monomer as the mechanism of SOD1 fibrillar aggregation. The non-monotonic TFE concentration dependence of protein conformational changes was explored through simulation studies. Our results suggest that altered protein conformation and dynamics within its structure may underlie the aggregation of SOD1 in ALS.
Collapse
Affiliation(s)
- Amresh Prakash
- a School of Computational & Integrative Sciences , Jawaharlal Nehru University , New Delhi , 110067 , India
| | - Vijay Kumar
- b Centre for Interdisciplinary Research in Basic Sciences , Jamia Millia Islamia , Jamia Nagar, New Delhi , 110025 , India
| | - Preeti Pandey
- a School of Computational & Integrative Sciences , Jawaharlal Nehru University , New Delhi , 110067 , India
| | - Deepak R Bharti
- a School of Computational & Integrative Sciences , Jawaharlal Nehru University , New Delhi , 110067 , India
| | - Poonam Vishwakarma
- a School of Computational & Integrative Sciences , Jawaharlal Nehru University , New Delhi , 110067 , India
| | - Ruhar Singh
- a School of Computational & Integrative Sciences , Jawaharlal Nehru University , New Delhi , 110067 , India
| | - Md Imtaiyaz Hassan
- b Centre for Interdisciplinary Research in Basic Sciences , Jamia Millia Islamia , Jamia Nagar, New Delhi , 110025 , India
| | - Andrew M Lynn
- a School of Computational & Integrative Sciences , Jawaharlal Nehru University , New Delhi , 110067 , India
| |
Collapse
|
30
|
Dammers C, Reiss K, Gremer L, Lecher J, Ziehm T, Stoldt M, Schwarten M, Willbold D. Pyroglutamate-Modified Amyloid-β(3-42) Shows α-Helical Intermediates before Amyloid Formation. Biophys J 2017; 112:1621-1633. [PMID: 28445753 PMCID: PMC5406372 DOI: 10.1016/j.bpj.2017.03.007] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2016] [Revised: 03/02/2017] [Accepted: 03/08/2017] [Indexed: 01/15/2023] Open
Abstract
Pyroglutamate-modified amyloid-β (pEAβ) has been described as a relevant Aβ species in Alzheimer's-disease-affected brains, with pEAβ (3-42) as a dominant isoform. Aβ (1-40) and Aβ (1-42) have been well characterized under various solution conditions, including aqueous solutions containing trifluoroethanol (TFE). To characterize structural properties of pEAβ (3-42) possibly underlying its drastically increased aggregation propensity compared to Aβ (1-42), we started our studies in various TFE-water mixtures and found striking differences between the two Aβ species. Soluble pEAβ (3-42) has an increased tendency to form β-sheet-rich structures compared to Aβ (1-42), as indicated by circular dichroism spectroscopy data. Kinetic assays monitored by thioflavin-T show drastically accelerated aggregation leading to large fibrils visualized by electron microscopy of pEAβ (3-42) in contrast to Aβ (1-42). NMR spectroscopy was performed for backbone and side-chain chemical-shift assignments of monomeric pEAβ (3-42) in 40% TFE solution. Although the difference between pEAβ (3-42) and Aβ (1-42) is purely N-terminal, it has a significant impact on the chemical environment of >20% of the total amino acid residues, as revealed by their NMR chemical-shift differences. Freshly dissolved pEAβ (3-42) contains two α-helical regions connected by a flexible linker, whereas the N-terminus remains unstructured. We found that these α-helices act as a transient intermediate to β-sheet and fibril formation of pEAβ (3-42).
Collapse
Affiliation(s)
- Christina Dammers
- Institute of Complex Systems (ICS-6) Structural Biochemistry, Forschungszentrum Jülich, Jülich, Germany
| | - Kerstin Reiss
- Institute of Complex Systems (ICS-6) Structural Biochemistry, Forschungszentrum Jülich, Jülich, Germany
| | - Lothar Gremer
- Institute of Complex Systems (ICS-6) Structural Biochemistry, Forschungszentrum Jülich, Jülich, Germany; Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Justin Lecher
- Institute of Complex Systems (ICS-6) Structural Biochemistry, Forschungszentrum Jülich, Jülich, Germany; Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Tamar Ziehm
- Institute of Complex Systems (ICS-6) Structural Biochemistry, Forschungszentrum Jülich, Jülich, Germany
| | - Matthias Stoldt
- Institute of Complex Systems (ICS-6) Structural Biochemistry, Forschungszentrum Jülich, Jülich, Germany; Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany
| | - Melanie Schwarten
- Institute of Complex Systems (ICS-6) Structural Biochemistry, Forschungszentrum Jülich, Jülich, Germany
| | - Dieter Willbold
- Institute of Complex Systems (ICS-6) Structural Biochemistry, Forschungszentrum Jülich, Jülich, Germany; Institut für Physikalische Biologie, Heinrich-Heine-Universität Düsseldorf, Düsseldorf, Germany.
| |
Collapse
|
31
|
Gimenez D, Dose A, Robson NL, Sandford G, Cobb SL, Coxon CR. 2,2,2-Trifluoroethanol as a solvent to control nucleophilic peptide arylation. Org Biomol Chem 2017; 15:4081-4085. [DOI: 10.1039/c7ob00295e] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Trifluoroethanol enables selective perfluoro-heteroaromatic peptide arylation of cysteine residues.
Collapse
Affiliation(s)
| | - Anica Dose
- Durham University
- Department of Chemistry
- Durham
- UK
| | | | | | | | - Christopher R. Coxon
- School of Pharmacy and Biomolecular Sciences
- Byrom Street Campus
- Liverpool John Moores University
- Liverpool
- UK
| |
Collapse
|
32
|
Chand A, Chowdhuri S. Behaviour of aqueous N-methylacetamide solution in presence of ethanol and 2,2,2 tri-fluoroethanol: Hydrogen bonding structure and dynamics. J Mol Liq 2016. [DOI: 10.1016/j.molliq.2016.10.129] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
33
|
Ohgi H, Imamura H, Yonenaga K, Morita T, Nishikawa K, Westh P, Koga Y. The effect of 2,2,2-trifluoroethanol on water studied by using third derivatives of Gibbs energy, G. J Mol Liq 2016. [DOI: 10.1016/j.molliq.2016.09.101] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
34
|
Biswas B, Mondal S, Singh PC. Role of Hydrogen Bond in the Solvation Behavior of the Binary Mixtures Containing Fluorocarbon Alcohol Molecules and Chloroform: An Experimental and Theoretical Study. ChemistrySelect 2016. [DOI: 10.1002/slct.201601118] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022]
Affiliation(s)
- Biswajit Biswas
- Department of Spectroscopy; Indian Association for the Cultivation of Science; Kolkata 700032 India
| | - Saptarsi Mondal
- Department of Spectroscopy; Indian Association for the Cultivation of Science; Kolkata 700032 India
| | - Prashant C. Singh
- Department of Spectroscopy; Indian Association for the Cultivation of Science; Kolkata 700032 India
| |
Collapse
|
35
|
Amyloid-like aggregates formation by bovine apo-carbonic anhydrase in various alcohols: A comparative study. Int J Biol Macromol 2016; 92:573-580. [DOI: 10.1016/j.ijbiomac.2016.07.083] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2016] [Revised: 07/11/2016] [Accepted: 07/24/2016] [Indexed: 12/31/2022]
|
36
|
|
37
|
Excess properties of non-ideal binary mixtures containing water, methanol and ethanol by molecular simulation. J Mol Liq 2015. [DOI: 10.1016/j.molliq.2015.08.058] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
38
|
|
39
|
Hydrophobicity enhancement in micelles of Triton X-165 by the presence of the cosolvent 2,2,2 trifluoroethanol (TFE). J Mol Liq 2014. [DOI: 10.1016/j.molliq.2014.08.017] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
40
|
Wang Z, Han X, He N, Chen Z, Brooks CL. Environmental Effect on Surface Immobilized Biological Molecules. J Phys Chem B 2014; 118:12176-85. [DOI: 10.1021/jp508550d] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Zunliang Wang
- State
Key Laboratory of Bioelectronics, School of Biological Science and
Medical Engineering, Southeast University, Si Pai Lou 2, Nanjing 210096, China
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Xiaofeng Han
- State
Key Laboratory of Bioelectronics, School of Biological Science and
Medical Engineering, Southeast University, Si Pai Lou 2, Nanjing 210096, China
| | - Nongyue He
- State
Key Laboratory of Bioelectronics, School of Biological Science and
Medical Engineering, Southeast University, Si Pai Lou 2, Nanjing 210096, China
| | - Zhan Chen
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Charles L. Brooks
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
41
|
Di Pierro M, Mugnai ML, Elber R. Optimizing potentials for a liquid mixture: a new force field for a tert-butanol and water solution. J Phys Chem B 2014; 119:836-49. [PMID: 25066823 PMCID: PMC4306501 DOI: 10.1021/jp505401m] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
![]()
A technology for optimization of
potential parameters from condensed-phase simulations (POP) is discussed
and illustrated. It is based on direct calculations of the derivatives
of macroscopic observables with respect to the potential parameters.
The derivatives are used in a local minimization scheme, comparing
simulated and experimental data. In particular, we show that the Newton
trust region protocol allows for more accurate and robust optimization.
We apply the newly developed technology to study the liquid mixture
of tert-butanol and water. We are able to obtain,
after four iterations, the correct phase behavior and accurately predict
the value of the Kirkwood Buff (KB) integrals. We further illustrate
that a potential that is determined solely by KB information, or the
pair correlation function, is not necessarily unique.
Collapse
Affiliation(s)
- Michele Di Pierro
- Institute for Computational Engineering and Sciences and ‡Department of Chemistry, University of Texas at Austin , Austin, Texas 78712, United States
| | | | | |
Collapse
|
42
|
Mijaković M, Polok KD, Kežić B, Sokolić F, Perera A, Zoranić L. A comparison of force fields for ethanol–water mixtures. MOLECULAR SIMULATION 2014. [DOI: 10.1080/08927022.2014.923567] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
43
|
Giambaşu GM, Luchko T, Herschlag D, York DM, Case DA. Ion counting from explicit-solvent simulations and 3D-RISM. Biophys J 2014; 106:883-94. [PMID: 24559991 PMCID: PMC3944826 DOI: 10.1016/j.bpj.2014.01.021] [Citation(s) in RCA: 83] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2013] [Revised: 12/19/2013] [Accepted: 01/10/2014] [Indexed: 12/30/2022] Open
Abstract
The ionic atmosphere around nucleic acids remains only partially understood at atomic-level detail. Ion counting (IC) experiments provide a quantitative measure of the ionic atmosphere around nucleic acids and, as such, are a natural route for testing quantitative theoretical approaches. In this article, we replicate IC experiments involving duplex DNA in NaCl(aq) using molecular dynamics (MD) simulation, the three-dimensional reference interaction site model (3D-RISM), and nonlinear Poisson-Boltzmann (NLPB) calculations and test against recent buffer-equilibration atomic emission spectroscopy measurements. Further, we outline the statistical mechanical basis for interpreting IC experiments and clarify the use of specific concentration scales. Near physiological concentrations, MD simulation and 3D-RISM estimates are close to experimental results, but at higher concentrations (>0.7 M), both methods underestimate the number of condensed cations and overestimate the number of excluded anions. The effect of DNA charge on ion and water atmosphere extends 20-25 Å from its surface, yielding layered density profiles. Overall, ion distributions from 3D-RISMs are relatively close to those from corresponding MD simulations, but with less Na(+) binding in grooves and tighter binding to phosphates. NLPB calculations, on the other hand, systematically underestimate the number of condensed cations at almost all concentrations and yield nearly structureless ion distributions that are qualitatively distinct from those generated by both MD simulation and 3D-RISM. These results suggest that MD simulation and 3D-RISM may be further developed to provide quantitative insight into the characterization of the ion atmosphere around nucleic acids and their effect on structure and stability.
Collapse
Affiliation(s)
- George M Giambaşu
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey
| | - Tyler Luchko
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey
| | - Daniel Herschlag
- Department of Biochemistry, Stanford University, Stanford, California
| | - Darrin M York
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey.
| | - David A Case
- Department of Chemistry and Chemical Biology and BioMaPS Institute, Rutgers University, Piscataway, New Jersey.
| |
Collapse
|
44
|
Gerig JT. Toward a Molecular Dynamics Force Field for Simulations of 40% Trifluoroethanol–Water. J Phys Chem B 2014; 118:1471-80. [DOI: 10.1021/jp408879g] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- J. T. Gerig
- Department of Chemistry and
Biochemistry, University of California, Santa Barbara, Santa Barbara, California 93106, United States
| |
Collapse
|
45
|
Lin B, Lopes PEM, Roux B, MacKerell AD. Kirkwood-Buff analysis of aqueous N-methylacetamide and acetamide solutions modeled by the CHARMM additive and Drude polarizable force fields. J Chem Phys 2013; 139:084509. [PMID: 24007020 PMCID: PMC3772949 DOI: 10.1063/1.4818731] [Citation(s) in RCA: 30] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2013] [Accepted: 08/04/2013] [Indexed: 11/14/2022] Open
Abstract
Kirkwood-Buff analysis was performed on aqueous solutions of N-methylacetamide and acetamide using the Chemistry at HARvard Molecular Mechanics additive and Drude polarizable all-atom force fields. Comparison of a range of properties with experimental results, including Kirkwood-Buff integrals, excess coordination numbers, solution densities, partial molar values, molar enthalpy of mixing, showed both models to be well behaved at higher solute concentrations with the Drude model showing systematic improvement at lower solution concentrations. However, both models showed difficulties reproducing experimental activity derivatives and the excess Gibbs energy, with the Drude model performing slightly better. At the molecular level, the improved agreement of the Drude model at low solute concentrations is due to increased structure in the solute-solute and solute-solvent interactions. The present results indicate that the explicit inclusion of electronic polarization leads to improved modeling of dilute solutions even when those properties are not included as target data during force field optimization.
Collapse
Affiliation(s)
- Bin Lin
- Department of Pharmaceutical Sciences, University of Maryland, School of Pharmacy, 20 Penn Street HSFII, Baltimore, Maryland 21201, USA
| | | | | | | |
Collapse
|
46
|
Equilibrium and folding simulations of NS4B H2 in pure water and water/2,2,2-trifluoroethanol mixed solvent: examination of solvation models. J Mol Model 2013; 19:3931-9. [DOI: 10.1007/s00894-013-1933-6] [Citation(s) in RCA: 10] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2013] [Accepted: 06/23/2013] [Indexed: 10/26/2022]
|
47
|
Kežić B, Perera A. Revisiting aqueous-acetone mixtures through the concept of molecular emulsions. J Chem Phys 2012; 137:134502. [DOI: 10.1063/1.4755816] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
48
|
Perera A, Mazighi R, Kežić B. Fluctuations and micro-heterogeneity in aqueous mixtures. J Chem Phys 2012; 136:174516. [DOI: 10.1063/1.4707745] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
49
|
Fyta M, Netz RR. Ionic force field optimization based on single-ion and ion-pair solvation properties: Going beyond standard mixing rules. J Chem Phys 2012; 136:124103. [DOI: 10.1063/1.3693330] [Citation(s) in RCA: 121] [Impact Index Per Article: 10.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
50
|
Fyta M. Structural and technical details of the Kirkwood-Buff integrals from the optimization of ionic force fields: focus on fluorides. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2012; 35:1-12. [PMID: 22434150 DOI: 10.1140/epje/i2012-12021-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/13/2011] [Revised: 03/01/2012] [Accepted: 03/05/2012] [Indexed: 05/31/2023]
Abstract
Results on the structural details of Kirkwood-Buff integrals obtained from the optimization of ionic force fields are presented. We have proposed and make use of an optimization scheme for ionic force fields, which is based on the modification of the cation-anion mixing rules, the calculation of the thermodynamics properties of various monovalent salt solutions according to the Kirkwood-Buff theory of solutions and the comparison to relevant experimental findings. Here, we complete and extend our calculations and analysis as we focus on the technical details of this optimization procedure and the case of fluorides, which have been proven difficult to handle. Important insight is given on the dependence of the radial distribution functions, the short-ranged potentials of mean force, and the Kirkwood-Buff integrals of the salt solutions on the different scaling factors in the mixing rules. Specifically, the way the structural details and inherent characteristics of the above properties are affected by the quantitative and qualitative differences in the mixing rules for a variety of common biologically relevant monovalent salts is mainly addressed. We conclude on the efficiency of this scheme, again with a focus on the fluorides. In the end, we provide a variation of the ion-pair mixing rules scaling factors with salt concentration to identify regimes for which different mixing rules prefactors lead to well-optimized force fields. All results are obtained through Molecular Dynamics simulations using previously optimized force fields for the monovalent ions.
Collapse
Affiliation(s)
- M Fyta
- Physics Department (T37), Technical University of Munich, Garching, Germany.
| |
Collapse
|