Scalable and fast heterogeneous molecular simulation with predictive parallelization schemes

GENESIS CGDYN: large-scale coarse-grained MD simulation with dynamic load balancing for heterogeneous biomolecular systems

Residue-level coarse-grained (CG) molecular dynamics (MD) simulation is widely used to investigate slow biological processes that involve multiple proteins, nucleic acids, and their complexes. Biomolecules in a large simulation system are distributed non-uniformly, limiting computational efficiency with conventional methods. Here, we develop a hierarchical domain decomposition scheme with dynamic load balancing for heterogeneous biomolecular systems to keep computational efficiency even after drastic changes in particle distribution. These schemes are applied to the dynamics of intrinsically disordered protein (IDP) droplets. During the fusion of two droplets, we find that the changes in droplet shape correlate with the mixing of IDP chains. Additionally, we simulate large systems with multiple IDP droplets, achieving simulation sizes comparable to those observed in microscopy. In our MD simulations, we directly observe Ostwald ripening, a phenomenon where small droplets dissolve and their molecules redeposit into larger droplets. These methods have been implemented in CGDYN of the GENESIS software, offering a tool for investigating mesoscopic biological processes using the residue-level CG models.

Orchestrating Extracellular Vesicle With Dual Reporters for Imaging and Capturing in Mammalian Cell Culture

Background: Recent technological advancements have enabled live-cell imaging of intracellular organelles to monitor their biogenesis in mammalian cells. However, applying this method to gain insight into extracellular organelles, such as extracellular vesicles (EVs), presents unique challenges that require special considerations in design and engineering.

Results: We have developed a dual-reporter system that combines genetic fusion, fluorescence microcopy and magnetic beads capture of EVs to study the biogenesis of EVs in mammalian cell cultures. First, we genetically produced a series of reporters by fusing a green fluorescent protein (GFP) and an affinity peptide (6xHis), with either the endogenous transmembrane protein, CD63, or EVs targeting vesicular stomatitis viral glycoprotein (VSVG). Transfection of these reporters into human 293T cells resulted in expression and integration of these reporters into pre-exosome compartments, which were subsequently released into the culture medium. Confocal imaging and nano-particle tracking analysis demonstrated that EVs were appropriately labeled and exhibited a single dominant peak in the 80–110 nm size range, indicating that isolated EVs were comprised of micro-vesicles and/or exosome subpopulations. Incubation of isolated EVs with nickel-coated magnetic beads resulted in successful capture of GFP-positive EVs. Finally, addition of EVs into culture medium was able to reveal the cellular uptake of GFP-labeled EVs by recipient cells. Taken together, our dual-reporter system provides a powerful method for both monitoring and capturing of EVs in mammalian cell culture systems.

Conclusion: A dual-reporter system provides a robust tool to study the life cycle of EVs in mammalian cells from biogenesis and excretion to cellular uptake.

SWINGER: a clustering algorithm for concurrent coupling of atomistic and supramolecular liquids

In this contribution, we review recent developments and applications of a dynamic clustering algorithm SWINGER tailored for the multiscale molecular simulations of biomolecular systems. The algorithm on-the-fly redistributes solvent molecules among supramolecular clusters. In particular, we focus on its applications in combination with the adaptive resolution scheme, which concurrently couples atomistic and coarse-grained molecular representations. We showcase the versatility of our multiscale approach on a few applications to biomolecular systems coupling atomistic and supramolecular water models such as the well-established MARTINI and dissipative particle dynamics models and provide an outlook for future work.

Open-Boundary Molecular Dynamics of a DNA Molecule in a Hybrid Explicit/Implicit Salt Solution

The composition and electrolyte concentration of the aqueous bathing environment have important consequences for many biological processes and can profoundly affect the behavior of biomolecules. Nevertheless, because of computational limitations, many molecular simulations of biophysical systems can be performed only at specific ionic conditions: either at nominally zero salt concentration, i.e., including only counterions enforcing the system's electroneutrality, or at excessive salt concentrations. Here, we introduce an efficient molecular dynamics simulation approach for an atomistic DNA molecule at realistic physiological ionic conditions. The simulations are performed by employing the open-boundary molecular dynamics method that allows for simulation of open systems that can exchange mass and linear momentum with the environment. In our open-boundary molecular dynamics approach, the computational burden is drastically alleviated by embedding the DNA molecule in a mixed explicit/implicit salt-bathing solution. In the explicit domain, the water molecules and ions are both overtly present in the system, whereas in the implicit water domain, only the ions are explicitly present and the water is described as a continuous dielectric medium. Water molecules are inserted and deleted into/from the system in the intermediate buffer domain that acts as a water reservoir to the explicit domain, with both water molecules and ions free to enter or leave the explicit domain. Our approach is general and allows for efficient molecular simulations of biomolecules solvated in bathing salt solutions at any ionic strength condition.