Accurate prediction of macromolecular structure and function is a key challenge in many fields, from designing new nanomaterials to drug discovery. Nanoparticles, nanochannels, proteins, and other biomolecules are typically charged, and this charge dictates their interactions and function. It is intuitively expected that two identical particles will contain the same charge. However, we find that when multiple particles aggregate, a spontaneous and global redistribution of charge can cause initially identical particles to take on different roles.
Molecules and particles commonly attain charge via dissociation or association of ions – the same process that governs acid and base reactions. However, this process of charge transfer between the molecules and the solution is perturbed when other charged entities are present in the vicinity, for example, a DNA molecule can alter the charge distribution of a nearby protein. We have wondered to what extent these perturbations could add up, leading to possible new modes of interactions and new structures.
To find out, we implemented a simulation method that dynamically resolves these “charge regulation” effects. We expected that the charge will redistribute within individual particles, which has been proposed using theoretical arguments already in the 1950s. Surprisingly, however, the simulations predict that initially, identical nanoparticles will spontaneously exchange charge among themselves, resulting in disparate charging and the corresponding formation of asymmetric aggregates. Following this result, we conclude that charge regulation must be accounted for to accurately predict the structure formation of molecular and nano-scale systems. Our simulation method should significantly improve the accuracy of calculations used for drug discovery and the design of nanomaterials.
The computational method that dynamically simulates charge-regulation effects has been made freely accessible as a module for the LAMMPS molecular dynamics simulation package at "https://docs.lammps.org/fix_charge_regulation.html". Moreover, the simulation method is currently being implemented within Ludwig Lattice-Boltzmann software to enable accurate predictions of charge transport in nano-channels.