Computational Modeling of Ion Transport in Bulk and through a Nanopore Using the Drude Polarizable Force Field
Jigneshkumar Dahyabhai Prajapati,Crystal Mele,Mehmet Alphan Aksoyoglu,Mathias Winterhalter,Ulrich Kleinekathöfer,Ulrich Kleinekathöfer
DOI: https://doi.org/10.1021/acs.jcim.0c00389
IF: 6.162
2020-06-01
Journal of Chemical Information and Modeling
Abstract:In the past two decades, molecular dynamics simulations have become the method of choice for elucidating the transport mechanisms of ions through various membrane channels. Often, these simulations heavily rely on classical nonpolarizable force fields (FFs), which lack electronic polarizability in the treatment of the electrostatics. The recent advancements in the Drude polarizable FF lead to a complete set of parameters for water, ions, protein, and lipids, allowing for a more realistic modeling of membrane proteins. However, the quality of these Drude FFs remains untested for such systems. Here, we examine the quality of this FF set in two ways, i.e., (i) in simple ionic aqueous solution simulations and (ii) in more complex membrane channel simulations. First, the aqueous solutions of KCl, NaCl, MgCl<sub>2</sub>, and CaCl<sub>2</sub> salts are simulated using the polarizable Drude and the nonpolarizable CHARMM36 FFs. The bulk conductivity has been estimated for both FF sets using applied-field simulations for several concentrations and temperatures in the case of all investigated salts and compared to experimental findings. An excellent improvement in the ability of the Drude FF to reproduce the experimental bulk conductivities for KCl, NaCl, and MgCl<sub>2</sub> solutions can be observed but not in the case of CaCl<sub>2</sub>. Moreover, the outer membrane channel OmpC from the bacterium <i>Escherichia coli</i> has been employed to examine the ability of the polarizable and nonpolarizable FFs to reproduce ion transport-related quantities known from experiment. Unbiased and applied-field simulations have been performed in the presence of KCl using both FF sets. Unlike for the bulk systems of aqueous salt solutions, it has been found that the Drude FF is not accurate in modeling KCl transport properties across the OmpC porin.The Supporting Information is available free of charge at <a class="ext-link" href="/doi/10.1021/acs.jcim.0c00389?goto=supporting-info">https://pubs.acs.org/doi/10.1021/acs.jcim.0c00389</a>.Figures of concentration- and temperature-dependent bulk conductivity for cations and anions and dynamics of loop L3 observed in presence of external voltages (<a class="ext-link" href="/doi/suppl/10.1021/acs.jcim.0c00389/suppl_file/ci0c00389_si_001.pdf">PDF</a>)This article has not yet been cited by other publications.
chemistry, multidisciplinary, medicinal,computer science, interdisciplinary applications, information systems