All-atom molecular dynamics study of hepatitis B virus containing pregenome RNA in solution

Understanding Virus Structure and Dynamics through Molecular Simulations

Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.

An exa-scale high-performance molecular dynamics simulation program: MODYLAS

A new version of the highly parallelized general-purpose molecular dynamics (MD) simulation program MODYLAS with high performance on the Fugaku computer was developed. A benchmark test using Fugaku indicated highly efficient communication, single instruction, multiple data (SIMD) processing, and on-cache arithmetic operations. The system's performance deteriorated only slightly, even under high parallelization. In particular, a newly developed minimum transferred data method, requiring a significantly lower amount of data transfer compared to conventional communications, showed significantly high performance. The coordinates and forces of 101 810 176 atoms and the multipole coefficients of the subcells could be distributed to the 32 768 nodes (1 572 864 cores) in 2.3 ms during one MD step calculation. The SIMD effective instruction rates for floating-point arithmetic operations in direct force and fast multipole method (FMM) calculations measured on Fugaku were 78.7% and 31.5%, respectively. The development of a data reuse algorithm enhanced the on-cache processing; the cache miss rate for direct force and FMM calculations was only 2.74% and 1.43%, respectively, on the L1 cache and 0.08% and 0.60%, respectively, on the L2 cache. The modified MODYLAS could complete one MD single time-step calculation within 8.5 ms for the aforementioned large system. Additionally, the program contains numerous functions for material research that enable free energy calculations, along with the generation of various ensembles and molecular constraints.

Dynamic Monte Carlo calculation generating particle trajectories that satisfy the diffusion equation for heterogeneous systems with a position-dependent diffusion coefficient and free energy

A series of new Monte Carlo (MC) transition probabilities was investigated that could produce molecular trajectories statistically satisfying the diffusion equation with a position-dependent diffusion coefficient and potential energy. The MC trajectories were compared with the numerical solution of the diffusion equation by calculating the time evolution of the probability distribution and the mean first passage time, which exhibited excellent agreement. The method is powerful when investigating, for example, the long-distance and long-time global transportation of a molecule in heterogeneous systems by coarse-graining them into one-particle diffusive molecular motion with a position-dependent diffusion coefficient and free energy. The method can also be applied to many-particle dynamics.