Models for twistable elastic polymers in Brownian dynamics, and their implementation for LAMMPS

C. A. Brackley,A. N. Morozov,D. Marenduzzo

DOI: https://doi.org/10.1063/1.4870088

2014-04-07

Abstract:An elastic rod model for semi-flexible polymers is presented. Theory for a continuum rod is reviewed, and it is shown that a popular discretised model used in numerical simulations gives the correct continuum limit. Correlation functions relating to both bending and twisting of the rod are derived for both continuous and discrete cases, and results are compared with numerical simulations. Finally, two possible implementations of the discretised model in the multi-purpose molecular dynamics software package LAMMPS are described.

Biological Physics,Soft Condensed Matter,Biomolecules

What problem does this paper attempt to address?

### What problems does this paper attempt to solve?

This paper aims to solve the challenge of accurately describing the bending and torsional properties of semi - flexible polymers (such as DNA) in coarse - grained simulations. Specifically, the paper focuses on the following issues:

1. **Limitations of existing models**:
- Many modern and widely - used molecular dynamics codes (such as LAMMPS) do not yet include force fields capable of handling the supercoiled structure of DNA.
- There are multiple methods for studying torsionally elastic chains in the literature, but there is a lack of systematic analysis of the relationships between these methods.

2. **Combination of theory and implementation**:
- Provide detailed continuous elastic rod theory and link it to the discretized energy function to ensure that the discrete model is correct in the continuous limit.
- Derive relevant functions for describing the bending and torsion of elastic rods and compare their theoretical results with numerical simulation results.

3. **Implementation in practical applications**:
- Describe two possible ways of implementing discrete models in LAMMPS to ensure that these models can be used in existing molecular dynamics software.
- Verify the performance of these models in simulations, especially whether they can correctly capture the bending and torsional properties of DNA.

4. **Biological significance**:
- Emphasize the importance of the supercoiled structure of DNA in bacteria and eukaryotic cells, including processes such as gene regulation, transcription, and chromatin structure.
- Point out the formation and regulation mechanisms of the supercoiled structure and the roles of related enzymes (such as topoisomerases and DNA gyrases).

### Specific problems

- **How to establish an elastic rod model that can describe both bending and torsion?**
- **How to ensure the correctness of the discretized model in the continuous limit?**
- **How to implement these models in existing molecular dynamics software (such as LAMMPS)?**
- **How to verify the effectiveness of these models through theoretical calculations and numerical simulations?**

### Solutions

The author has solved the above problems through the following steps:

1. **Theoretical derivation**:
- Review the theory of continuous elastic rods and prove that commonly used discretized models are correct in the continuous limit.
- Derive relevant functions for elastic rod bending and torsion in continuous and discrete cases and compare the results with numerical simulations.

2. **Model implementation**:
- Describe two possible ways of implementing discretized models in LAMMPS: one uses dihedral angle interactions and "patch" beads (Model 1), and the other uses spherical atoms (Model 2).
- Provide detailed formulas and algorithms to ensure that these models can run correctly in LAMMPS.

3. **Verification and comparison**:
- Verify the performance of these models through numerical simulations, especially in terms of average bending angle, bending correlation length, average torsional angle, and torsional correlation length.
- Compare theoretical predictions with simulation results to verify the validity and accuracy of the models.

Through these efforts, the author has not only filled the gaps in the literature but also provided reliable tools and methods for future research.