1
|
Na H, Song G. Coarse-Graining Waters: Unveiling The Effective Hydrophilicity/Hydrophobicity of Individual Protein Atoms and The Roles of Waters' Hydrogens. J Chem Theory Comput 2023; 19:7307-7323. [PMID: 37782694 PMCID: PMC10601925 DOI: 10.1021/acs.jctc.3c00700] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2023] [Indexed: 10/04/2023]
Abstract
There have been many coarse-graining methods developed that aim to reduce the sizes of simulated systems and their computational costs. In this work, we develop a new coarse-graining method, called coarse-graining-delta (or δ-CG in short), that reduces the degrees of freedom of the potential energy surface by coarse-graining relative locations of atoms from their unit centers. Our method extends and generalizes the methods used in the coarse-grained normal mode analysis and enables us to study the roles of the individual removed atoms in a system, which have been difficult to study in molecular dynamics simulations. By applying δ-CG to coarse-grain three-point water molecules into single-point solvent particles, we successfully identify the effective hydrophilicity and hydrophobicity of all the individual protein atom types, which collectively correlate well with the known hydrophilic, hydrophobic, and amphipathic characteristics of amino acids. Moreover, our investigation shows that water's hydrogens have two roles in interacting with protein atoms. First, water molecules adjust their poses around different amino acids and their atoms, and the statistical preferences of the hydrogen poses near the atoms determine the effective hydrophilicity and hydrophobicity of the atoms, which have not been successfully addressed before. Second, the collective dynamics of the hydrogens assist the water molecules in escaping from the potential energy wells of the hydrophilic atoms. Our method also shows that coarse-graining a system mathematically leads to breaking antisymmetry of the nonbonded interactions; as a result, two interacting coarse-grained units exert different forces on each other. Our study indicates that the accuracy of coarse-grained force fields, such as the MARTINI force field and the UNRES force field, can be improved in two ways: (i) refining their potential energy functions and coefficients by analyzing the coarse-grained potential energy surface using δ-CG, and (ii) introducing non-antisymmetric interactions.
Collapse
Affiliation(s)
- Hyuntae Na
- Department
of Computer Science, Penn State Harrisburg, Middletown, Pennsylvania 17057, United States
| | - Guang Song
- Department
of Mathematics and Computer Science, Westmont
College, Santa
Barbara, California 93108, United States
| |
Collapse
|
3
|
Ibarra-Tandi B, Moreno-Razo J, Munguía-Valadez J, López-Lemus J, Chávez-Rojo M. Effects of the repulsive and attractive forces on phase equilibrium and critical properties of two-dimensional non-conformal simple fluids. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2020.115234] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
4
|
de Souza RM, Ratochinski RH, Karttunen M, Dias LG. Self-Assembly of Phosphocholine Derivatives Using the ELBA Coarse-Grained Model: Micelles, Bicelles, and Reverse Micelles. J Chem Inf Model 2020; 60:522-536. [PMID: 31714768 DOI: 10.1021/acs.jcim.9b00790] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
The ELBA coarse-grained force field was originally developed for lipids, and its water model is described as a single-site Lennard-Jones particle with electrostatics modeled by an embedded point-dipole, while other molecules in this force field have a three (or four)-to-one mapping scheme. Here, ELBA was applied to investigate the self-assembly processes of dodecyl-phosphocholine (DPC) micelle, 1,2-dipalmitoyl-sn-glycero-3-phosphocholine/1,2-dihexaoyl-sn-glycero-3-phosphocholine (DPPC/DHPC) bicelles, and DPPC/cyclohexane/water reverse micelles through coarse-grained molecular dynamics (MD) simulations. New parameters were obtained using a simplex algorithm-based calibration procedure to determine the Lennard-Jones parameters for cyclohexane, dodecane, and cyclohexane-dodecane cross-interactions. Density, self-diffusion coefficient, surface tension, and mixture excess volume were found to be in fair agreement with experimental data. These new parameters were used in the simulations, and the obtained structures were analyzed for shape, size, volume, and surface area. Except for the shape of DPC micelles, all other properties match well with available experimental data and all-atom simulations. Remarkably, in agreement with experiments the rodlike shape of the DPPC reverse micelle is well described by ELBA, while all-atom data in the literature predicts a disclike shape. To further check the consistency of the force field in reproducing the correct shapes of reverse micelles, additional simulations were performed doubling the system size. Two distinct reverse micelles were obtained both presenting the rodlike shape and correct aggregation number.
Collapse
Affiliation(s)
- R M de Souza
- Department of Chemistry , The University of Western Ontario , London , Ontario , Canada N6A 3K7.,Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil.,The Center for Advanced Materials and Biomaterials Research , The University of Western Ontario , London , Ontario , Canada N6K 3K7
| | - R H Ratochinski
- Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil
| | - Mikko Karttunen
- Department of Chemistry , The University of Western Ontario , London , Ontario , Canada N6A 3K7.,The Center for Advanced Materials and Biomaterials Research , The University of Western Ontario , London , Ontario , Canada N6K 3K7.,Department of Applied Mathematics , The University of Western Ontario , London , Ontario , Canada N6A 5B7
| | - L G Dias
- Departamento de Química, FFCLRP , Universidade de São Paulo , Avenida Bandeirantes 3900 , 14040-901 Ribeirão Preto , SP , Brazil
| |
Collapse
|
5
|
Khalak Y, Baumeier B, Karttunen M. Improved general-purpose five-point model for water: TIP5P/2018. J Chem Phys 2018; 149:224507. [PMID: 30553255 DOI: 10.1063/1.5070137] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
A new five point potential for liquid water, TIP5P/2018, is presented along with the techniques used to derive its charges from ab initio per-molecule electrostatic potentials in the liquid phase using the split charge equilibration of Nistor et al. [J. Chem. Phys. 125, 094108 (2006)]. By taking the density and diffusion dependence on temperature as target properties, significant improvements to the behavior of isothermal compressibility were achieved along with improvements to other thermodynamic and rotational properties. While exhibiting a dipole moment close to ab initio values, TIP5P/2018 suffers from a too small quadrupole moment due to the charge assignment procedure and results in an overestimation of the dielectric constant.
Collapse
Affiliation(s)
- Yuriy Khalak
- Department of Mathematics and Computer Science and Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands
| | - Björn Baumeier
- Department of Mathematics and Computer Science and Institute for Complex Molecular Systems, Eindhoven University of Technology, P.O. Box 513, 5600MB Eindhoven, The Netherlands
| | - Mikko Karttunen
- Department of Chemistry, The University of Western Ontario, 1151 Richmond Street, London, Ontario N6A 3K7, Canada
| |
Collapse
|