Molecular simulations of charged complex fluids: A review

Theoretical Study of Fe3+ and Ni2+ Ion Interactions in Ethaline as the Deep Eutectic Solvent and Water Solutions Using Molecular Dynamics, Quantum Theory of Atoms in Molecules, and Non-Covalent Interactions

Deep eutectic solvents (DES) have several advantages compared to water and traditional solvents, making them an alternative, especially in applications that require better solvation, greater thermal stability, and a lower environmental impact. This study aimed to analyze the behavior of Fe3+ and Ni2+ ions in two quantities of water (300 and 5580 molecules) and in a eutectic solvent based on choline chloride and ethylene glycol (1ChCl:2EG). The computational methods used involved molecular dynamics, quantum theory of atoms in molecules (QTAIM), and noncovalent interactions simulations. Analysis of the radial distribution function multiplied by the number density [g(r)ρ] and the cumulative number (CN) indicated that the interactions between the metal ions and the water molecules were strongest for the systems with the most water. QTAIM determined the bond critical point, the electron density [ρ(r)], the Laplacian of the electronic density [∇2ρ(r)], and the electron localization function (η, ELF), allowing the interactions to be analyzed. Polarizability of the metal ions in both water and ethaline was compared; the increasing order of polarizability was Fe3+ < Ni2+ < Fe2+, with Fe3+ being the least polarizable due to its high charge and smaller ionic radius. In the mixed systems with Fe2+ or Fe3+ added to Ni2+, the metal species with the same charge competed similarly in water and DES. The intermolecular forces in DES are weaker due to the solvent's lower polarity and dielectric constant than water. In the systems with the highest water content (5580 molecules), Fe3+ ions were surrounded by the largest water molecules, followed by Fe2+ and Ni2+. These results may help to better understand the solvation and behavior of these ions in different media, which has implications for use in electrodeposition, batteries, and corrosion inhibition.

Unraveling the orientation of an enzyme adsorbed onto a metal-organic framework

Bio-conversion of lignocellulosic biomass to bioethanol fuel is a highly desirable yet challenging objective because of the low catalytic activity and high cost of β-glucosidase (BGL). Recently, ZIF-8, an emerging organic porous material, has been proposed as a promising candidate for enzyme immobilization to improve associated activity and stability. However, the underlying interaction mechanism of binding BGL on the ZIF-8 surface is yet to be clarified. Here, the adsorption of BGL onto ZIF-8 is explored for the first time by molecular dynamics simulations. The results show that BGL adsorbs on the ZIF-8 surface with a "back-on" orientation. The adsorption free energy analysis shows that the adsorption process is enthalpy driven. In addition, the electrostatic interaction between negatively charged residues and Zn2+ on the surface of ZIF-8 is found to play a decisive role in surface binding, which accounts for 98% of the total interaction energy. The secondary structure of BGL is not affected despite the strong adsorption, suggesting the good biocompatibility of ZIF-8. This study not only provides a reliable theoretical insight into understanding the interaction mechanism between BGL and ZIF-8, but also helps the rational design of ZIF-8-based materials for bio-related applications.

Finely tuned water structure and transport in functionalized carbon nanotube membranes during desalination

Molecular dynamics simulations were performed to tune the transport of water molecules in nanostructured membrane in a desalination process. Four armchair-type (7,7), (8,8), (9,9) and (10,10) carbon nanotubes (CNTs) with pore diameters around 1 nm were chosen, their interior surfaces were modified with -OH, -CH3 and -F groups. Simulation results show that water transport in nanochannel depends on confined water structures which could be regulated by precisely controlled channel diameter and chemical functionalization. Increasing CNT diameter changes water structures from single-file-like to be square and hexagonal-like, then into a disordered pattern, resulting in a concave-shaped trend of water permeance. The -OH functional groups promote structural ordering of water molecules in (7,7) CNT, but disrupt water structures in (8,8) and (9,9) CNTs, and reduce the order degree of water molecules in (10,10) CNT, moreover, exert an attraction to enhance surface friction inside channel. The -CH3 groups induce more strictly single-file movement of water molecules in (7,7) CNT, turning water structures in (8,8) and (9,9) CNTs into two and triangular column arrangements, improving water transport, however, causing again square-like water structure in (10,10) CNT. Fluorinations of CNT make water structure more disordered in (7,7), (9,9) and (10,10) CNTs, while enhance the square water structure in (8,8) CNT with a lower water permeance. Through changing channel diameter and functionalization, the low tetrahedral order corresponds to a more single-file-like water structure, associated with rapid water diffusion and high permeability; an increase in tetrahedrality results in more ice-like water structures, lower water diffusion coefficients, and permeability. The results of this study demonstrate that water transport could be finely regulated via a functionalized CNT membrane.

Modulating the adsorption orientation of methionine-rich laccase by tailoring the surface chemistry of single-walled carbon nanotubes

Achieving fast electron transfer process between oxidoreductase and electrodes is pivotal for the biocathode of enzymatic biofuel cells (EBFCs). However, in-depth understanding of the interplay mechanism between enzymes and electrode materials remains challenging when designing and constructing EBFCs. Herein, atomic-scale insight into the direct electron transfer (DET) behavior of Thermus thermophilus laccase (TtLac) with a special methionine-rich β-hairpin motif adsorbed on the carboxyl-functionalized carbon nanotube (COOH-CNT) and amino-functionalized carbon nanotube (NH₂-CNT) surfaces were disclosed by multi-scale molecular simulations. Simulation results reveal that electrostatic modification is an effective way to tune the DET behavior for TtLac on the modified-CNTs electrode surface. Surprisingly, the positively charged TtLac can be attracted by both negatively charged COOH-CNT and positively charged NH₂-CNT surfaces, yet only the latter is capable to trigger the DET process due to the 'lying-on' adsorption orientation. Specifically, the T1 copper site is near the methionine-rich β-hairpin motif, which is the key binding site for TtLac binding onto the NH₂-CNT surface via electrostatic interaction, π-π stacking and cation-π interaction. Moreover, TtLac on the NH₂-CNT surface undergoes less conformational changes than those on the COOH-CNT surface, which allows the laccase stability and catalytic efficiency to be well preserved. These findings provide a fundamental guidance for future design and fabrication of methionine-rich laccase-based EBFCs with high power output and long lifespan.