# DNA-DNA hybridization simulating

**URL:** <https://gromacs.bioexcel.eu/t/dna-dna-hybridization-simulating/8162>\
**Category:** User discussions\
**Created:** [January 28, 2024, 10:35am UTC](https://gromacs.bioexcel.eu/t/dna-dna-hybridization-simulating/8162 "2024-01-28T10:35:18Z")\
**Posts on this page:** 4\
**Page:** 1

<div class="post-metadata">

**Author:** ![Youzy](https://avatars.discourse-cdn.com/v4/letter/y/f08c70/32.png) [@Youzy](https://gromacs.bioexcel.eu/u/Youzy)\
**Post date:** [January 28, 2024, 10:35am UTC](https://gromacs.bioexcel.eu/t/dna-dna-hybridization-simulating/8162/1 "2024-01-28T10:35:18Z")

</div>

GROMACS version: 2020.6  
GROMACS modification: No

Hello guys. I’m a beginner who has recently entered this field and has limited knowledge about GROMACS. I just finished the beginner’s tutorial on lysozyme in water.

I’m interested in simulating DNA-DNA hybridization. In simple terms, I have a 100 nt single-stranded DNA (ssDNA) with a secondary structure (hairpin), as well as several 20-30 nt ssDNAs with sequences complementary to different regions of the 100 nt hairpin DNA. The term “hairpin” refers to the structure of the ssDNA resembling a hairpin, which is significant in some nucleic acid molecular pathways. My major is Nucleic Acid Molecular Circuit Engineering.

 ![hairpin](https://europe1.discourse-cdn.com/flex017/uploads/bioexcel1/original/2X/7/737e7f619fa1a63429ce2c5b7743f2643f6105e4.png)

I aim to use GROMACS MMPBSA to calculate the key base interactions when these 20-30 nt ssDNAs bind to the 100 nt hairpin DNA with a secondary structure.

I currently have two main questions:

1. My ssDNA has a hairpin structure, and in a sense, its ideal secondary structure is predetermined (designed based on the subsequent reactions I want to trigger). I’m wondering if I can directly create a pdb file with the secondary structure. However, I’m unsure which software can accomplish this part. Any suggestions would be appreciated.
2. If option 1 is not feasible, do I need to perform energy minimization first and then proceed with simulating the interaction with other structures? Or are there alternative approaches? Your guidance would be greatly appreciated. (Note: The temperature here is not room temperature.)

I’m certain that I’m looking for a specific DNA-DNA tutorial, but I can’t find relevant literature for my specific question. Most people seem to use GROMACS for simulations related to proteins. Perhaps someone involved in DNA or nucleic acid research can point me in the right direction and suggest resources that would aid my learning?

If you could provide valuable guidance, I would be extremely grateful.

Thank you in advance!!!

---

<div class="post-metadata">

**Author:** ![milosz.wieczor](https://avatars.discourse-cdn.com/v4/letter/m/b19c9b/32.png) [@milosz.wieczor](https://gromacs.bioexcel.eu/u/milosz.wieczor)\
**Post date:** [January 28, 2024, 9:33pm UTC](https://gromacs.bioexcel.eu/t/dna-dna-hybridization-simulating/8162/2 "2024-01-28T21:33:09Z")

</div>

My experience is by no means complete here, so perhaps I’m missing some promising approaches, but I’d say a combination of folding and hybridization looks like a use case for a general coarse-grained model such as 3SPN.2C (see [the FF parametrization paper](https://pubs.aip.org/aip/jcp/article-abstract/139/14/144903/352010/An-experimentally-informed-coarse-grained-3-site?redirectedFrom=fulltext)). There is also a DNA version of [Martini](https://pubs.acs.org/doi/full/10.1021/acs.jctc.5b00286) but I’m not aware of its widespread use for simulating DNA annealing/hybridization.

While Martini has a [Gromacs port and tutorials](https://www.compchems.com/gromacs-tutorial-coarse-grained-simulations-with-martini-3-force-field/), 3SPN.2C is available in other packages, most notably [OpenMM](https://github.com/cabb99/open3spn2), [Genesis](https://www.r-ccs.riken.jp/labs/cbrt/tutorials2019/tutorial-18-2/) and [CafeMol](https://www.cafemol.org/develop/).

I myself had always used [X3DNA](https://x3dna.org/) to create custom nucleic acid models (and continue to do so), but the software was recently incorporated into DSSR and commercialized.

---

<div class="post-metadata">

**Author:** ![Youzy](https://avatars.discourse-cdn.com/v4/letter/y/f08c70/32.png) [@Youzy](https://gromacs.bioexcel.eu/u/Youzy)\
**Post date:** [January 30, 2024, 1:59pm UTC](https://gromacs.bioexcel.eu/t/dna-dna-hybridization-simulating/8162/3 "2024-01-30T13:59:32Z")

</div>

> [@milosz.wieczor](#):
>
> I myself had always used [X3DNA](https://x3dna.org/) to create custom nucleic acid models (and continue to do so), but the software was recently incorporated into DSSR and commercialized.

Thank U!!! I heard X3DNA is uesful but our group… I find that this website [http://biophy.hust.edu.cn/new/3dRNA/create](http://biophy.hust.edu.cn/new/3dRNA/create) ^\_^ By the way, I’d like to use the force field AMBER99-ILDN

---

<div class="post-metadata">

**Author:** ![keziah](https://avatars.discourse-cdn.com/v4/letter/k/c0e974/32.png) [@keziah](https://gromacs.bioexcel.eu/u/keziah)\
**Post date:** [April 12, 2024, 12:57pm UTC](https://gromacs.bioexcel.eu/t/dna-dna-hybridization-simulating/8162/4 "2024-04-12T12:57:17Z")

</div>

hey…were u able to find a solution for this? As i work in a similar research area of nucleic acid networks and was struggling to find a simulation tool
