# SMD Pulling of DNA Toward Asymmetric Bilayer

**URL:** https://gromacs.bioexcel.eu/t/smd-pulling-of-dna-toward-asymmetric-bilayer/12038
**Category:** User discussions
**Tags:** simulation-setup
**Created:** [May 6, 2025, 8:14am UTC](https://gromacs.bioexcel.eu/t/smd-pulling-of-dna-toward-asymmetric-bilayer/12038 "2025-05-06T08:14:06Z")
**Posts on this page:** 1
**Showing post:** 9

<div class="post-metadata">

### Author: ![mostafagh98](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/mostafagh98/32/4214_2.png) [@mostafagh98](https://gromacs.bioexcel.eu/u/mostafagh98)
#### Post date: [May 8, 2025, 9:14am UTC](https://gromacs.bioexcel.eu/t/smd-pulling-of-dna-toward-asymmetric-bilayer/12038/9 "2025-05-08T09:14:32Z")

</div>

Dear @obZehn, thank you very much for your thoughtful and detailed feedback.

> [@obZehn](#):
>
> Do you expect water molecules to be there?

I definitely do not expect water molecules to enter the bilayer core. I tried removing them manually after solvation, but they reappear during equilibration. I’ve read about possible strategies like:  
Flat-bottomed position restraints on water molecules,  
But I’m unsure whether this method would introduce artifacts or negatively affect the membrane dynamics. Do you have experience with such techniques or any recommendations for keeping water out of the bilayer?  
Do you think I may have contributed to this by setting compressibility = 0.0 in the z-direction? I found it based on a suggestion [here](https://gromacs.bioexcel.eu/t/issue-with-direction-periodic-and-periodicity-during-umbrella-sampling/11552/4).

> [@obZehn](#):
>
> I do not know the nature of the membrane.

My membrane is not a nanodisc. It’s composed of:

Upper leaflet: 128 POPC, 27 Cholesterol, 15 GM3  
Lower leaflet: 110 POPC, 27 Cholesterol, 33 POPS  
(170 lipids per leaflet, total 340)

I made sure the area per lipid (APL) is balanced between the two leaflets before and after equilibration and that the bilayer plane spans at least 1.5 nm more than the DNA helical axis on each side. Do you think I should still increase the membrane size?

> [@obZehn](#):
>
> can a DNA strand physiologically partition through a bilayer?

> [@obZehn](#):
>
> is this expected to happen? Is this a physiological process/an experimental set up of some kind?

You’re absolutely right to ask — no, DNA does not physiologically partition through a bilayer on its own. In experiments, we use naked DNA as a negative control.

My actual goal is to later simulate DNA + peptide complexes, where different numbers of peptides are bound to the DNA. The current naked DNA simulation serves as a baseline for PMF comparison, to show that the DNA alone has a higher free energy barrier compared to the complexes.

I’ve seen similar PMF studies done using symmetric bilayers. But in those cases, the neutral charge of the membrane tends to favor DNA over complexes. That’s why I’m using an asymmetric bilayer to reflect more realistic surface charge conditions.

> [@obZehn](#):
>
> What is the expected behavior that you want to simulate?

I’m intending to pull the DNA along its long (helical) axis.

Regarding your point about applying force only to the center of mass (COM) of the DNA, I understand that this provides little control over the molecule’s orientation. I’ve seen COM-based pulling used in several studies, but I’m open to better alternatives.  
Would you recommend a more localized pulling group or a more refined collective variable?

The membrane rotation (from z-normal to y-normal) is concerning.

I’m considering adding the following lines in the `.mdp` to suppress membrane rotation:

```auto
nstcomm = 100
comm_mode = angular
comm_grps = MEMB

```

Would this be appropriate for a bilayer system, or would it introduce unintended side effects?

Thanks again for all your valuable insights.

---

_[View the full topic](https://gromacs.bioexcel.eu/t/smd-pulling-of-dna-toward-asymmetric-bilayer/12038)._