Despite decades of investigations, the principal mechanisms responsible for the high affinity and specificity of proteins for key physiological cations K+, Na+, and Ca2+ remain a hotly debated topic. At the core of the debate is an apparent need (or lack thereof) for an accurate description of the electrostatic response of the charge distribution in a protein to the binding of an ion. These effects range from partial electronic polarization of the directly ligating atoms to long-range effects related to partial charge transfer and electronic delocalization effects. While accurate modeling of cation recognition by metalloproteins warrants the use of quantum-mechanics (QM) calculations, the most popular approximations used in major biomolecular simulation packages rely on the implicit modeling of electronic polarization effects. That is, high-level QM computations for ion binding to proteins are desirable, but they are often unfeasible, because of the large size of the reactive-site models and the need to sample conformational space exhaustively at finite temperature. Several solutions to this challenge have been proposed in the field, ranging from the recently developed Drude polarizable force-field for simulations of metalloproteins to approximate tight-binding density functional theory (DFTB). To delineate the usefulness of different approximations, we examined the accuracy of three recent and commonly used theoretical models and numerical algorithms, namely, CHARMM C36, the latest developed Drude polarizable force fields, and DFTB3 with the latest 3OB parameters. We performed MD simulations for 30 cation-selective proteins with high-resolution X-ray structures to create ensembles of structures for analysis with different levels of theory, e.g., additive and polarizable force fields, DFTB3, and DFT. The results from DFT computations were used to benchmark CHARMM C36, Drude, and DFTB3 performance. The explicit modeling of quantum effects unveils the key electrostatic properties of the protein sites and the importance of specific ion-protein interactions. One of the most interesting findings is that secondary coordination shells of proteins are noticeably perturbed in a cation-dependent manner, showing significant delocalization and long-range effects of charge transfer and polarization upon binding Ca2+.

The thermodynamics of ion solvation in non-aqueous solvents remains of great significance for understanding cellular transport and ion homeostasis for the design of novel ion-selective materials and applications in molecular pharmacology. Molecular simulations play pivotal roles in connecting experimental measurements to the microscopic structures of liquids. One of the most useful and versatile mimetic systems for understanding biological ion transport is N-methyl-acetamide (NMA). A plethora of theoretical studies for ion solvation in NMA have appeared recently, but further progress is limited by two factors. One is an apparent lack of experimental data on solubility and thermodynamics of solvation for a broad panel of 1 : 1 salts over an appropriate temperature and concentration range. The second concern is more substantial and has to do with the limitations hardwired in the additive (fixed charge) approximations used for most of the existing force-fields. In this submission, we report on the experimental evaluation of LiCl solvation in NMA over a broad range of concentrations and temperatures and compare the results with those of MD simulations with several additive and one polarizable force-field (Drude). By comparing our simulations and experimental results to density functional theory computations, we discuss the limiting factors in existing potential functions. To evaluate the possible implications of explicit and implicit polarizability treatments on ion permeation across biological channels, we performed potential of mean force (PMF) computations for Li(+) transport through a model narrow ion channel with additive and polarizable force-fields.

Electronic properties of carbon nanotubes (CNTs) play an important role in their interactions with nano-structured materials. In this work, interactions of adenosine monophosphate (AMP), a DNA nucleotide, with metallic and semi-conducting CNTs are studied using the density functional tight binding (DFTB) method. The electronic structure of semi-conducting CNTs was found to be changed as they turned to metallic CNTs in a vacuum upon interaction with the nucleotide while metallic CNTs remain metallic. Specifically, the band gap of semi-conducting CNTs was decreased by 0.79 eV on average while nearly no change was found in the metallic tubes. However, our investigations showed that the presence of explicit water molecules prevents the metallicity change and only small changes in the CNT band gap occur. According to our charge analysis, the average negative charge accumulated on CNTs upon interaction with the AMP was determined to be 0.77 e in a vacuum while it was 0.03 e in solution. Therefore, it is essential to include explicit water molecules in simulating complexes formed by DNA nucleotides and CNTs which were ignored in several past studies performed using quantum mechanical approaches.

No abstract