Elvira Tarasova

All-Atom Molecular Dynamics Simulations of Entire Virus Capsid Reveal the Role of Ion Distribution in Capsid’s Stability

Present experimental methods do not have sufficient resolution to investigate all processes in virus particles at atomistic details. We report the results of molecular dynamics simulations and analyze the connection between the number of ions inside an empty capsid of PCV2 virus and its stability. We compare the crystallographic structures of the capsids with unresolved N-termini and without them in realistic conditions (room temperature and aqueous solution) and show that the structure is preserved. We find that the chloride ions play a key role in the stability of the capsid. A low number of chloride ions results in loss of the native icosahedral symmetry, while an optimal number of chloride ions create a neutralizing layer next to the positively charged inner surface of the capsid. Understanding the dependence of the capsid stability on the distribution of the ions will help clarify the details of the viral life cycle that is ultimately connected to the role of packaged viral genome inside the capsid.

Two-phase flow analogy as an effective boundary condition for modelling liquids at atomistic resolution

A hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework based on the analogy with two-phase hydrodynamics has been extended to dynamically tracking the feature of interest at all-atom resolution. In the model, the hydrodynamics description is used as an effective boundary condition to close the molecular dynamics solution without resorting to standard periodic boundary conditions. The approach is implemented in a popular Molecular Dynamics package GROMACS and results for two biomolecular systems are reported. A small peptide dialanine and a complete capsid of a virus porcine circovirus 2 in water are considered and shown to reproduce the structural and dynamic properties compared to those obtained in theory, purely atomistic simulations, and experiment.