Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink’s group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.

Trivalent actinides and their lanthanide homologues are being scrutinized for their potential health risk when ingested as a result of a range of industrial activities such as mining. Importantly, these ions are known to exhibit high affinity towards calmodulin (CaM). In case of their inadvertent uptake, the holoproteins that are occupied by these cations may block signal transduction pathways or increase the concentration of these ions in intact cells, which could lead to accumulation in human organs. Accordingly, this investigation employed spectroscopy, computational chemistry, calorimetry, and biochemistry to study the results of metal ion substitution on the protein structure, enzymatic activity and chemo- and cytotoxicity of An³⁺/Ln³⁺ ions. As will be demonstrated herein, our data confirm the higher affinity of Cm³⁺ and Eu³⁺ compared to Ca²⁺ to all 4 binding sites of CaM, with one site differing from the remaining three. This higher-affinity site will complex Eu³⁺ in an exothermic fashion; in contrast, ion binding to the three lower-affinity EF-hands was found to be endothermic. The overall endothermic binding process is ascribed to the loss of the hydration shells of the trivalent ions upon protein binding. These findings are supported by extensive quantum chemical calculations of full holo-CaM, which were performed at the MP2 level using the fragment molecular orbital method. The exceptional binding site (EF-hand 3) features fewer negatively charged residues compared to the other EF-hands, thereby allowing Eu³⁺ and Cm³⁺ to carry one or two additional waters compared to Ca²⁺-CaM, while also causing the structure of Cm³⁺/Eu³⁺-CaM to become slightly disordered. Moreover, the enzymatic activity decreases somewhat in comparison to Ca²⁺-CaM. By utilizing a combination of techniques, we were able to generate a comprehensive picture of the CaM-actinide/lanthanide system from the molecular level to its functional impact. Such knowledge could also be applied to other metal-binding proteins.