# 1-4 interactions

**URL:** https://gromacs.bioexcel.eu/t/1-4-interactions/1218
**Category:** User discussions
**Created:** [December 8, 2020, 1:39am UTC](https://gromacs.bioexcel.eu/t/1-4-interactions/1218 "2020-12-08T01:39:17Z")
**Posts on this page:** 5
**Page:** 1

<div class="post-metadata">

### Author: ![Qinghua\_Liao](https://avatars.discourse-cdn.com/v4/letter/q/46a35a/32.png) [@Qinghua\_Liao](https://gromacs.bioexcel.eu/u/Qinghua_Liao)
#### Post date: [December 8, 2020, 1:39am UTC](https://gromacs.bioexcel.eu/t/1-4-interactions/1218/1 "2020-12-08T01:39:17Z")

</div>

Hello,

I have a system of a protein and a ligand. When building topology, I used pdb2gmx for the protein and  
other tools for the ligand. I got a itp file for the ligand, then append it in the master topology files.

As Glycam force field (different scaling factor for 1-4 interactions) was used for the ligand, the parameters for 1-4 interactions of the ligand are explicitly included in the itp file, [pairs\_nb].

Then when running simulations, how are the 1-4 interactions within the ligand being calculated, using the explicit parameters in the itp file, or grompp will generate new parameters based on the protein force field (different scaling factor)? Any comments are welcome and appreciated!

All the best,  
Qinghua

---

<div class="post-metadata">

### Author: ![alevilla](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/alevilla/32/439_2.png) [@alevilla](https://gromacs.bioexcel.eu/u/alevilla)
#### Post date: [December 10, 2020, 10:00am UTC](https://gromacs.bioexcel.eu/t/1-4-interactions/1218/2 "2020-12-10T10:00:55Z")

</div>

Hi,  
“Then when running simulations, how are the 1-4 interactions within the ligand being calculated, using the explicit parameters in the itp file, or grompp will generate new parameters based on the protein force field (different scaling factor)? Any comments are welcome and appreciated!”

The last parameters defined in the topology (in your case the one in the directive [pairs\_nb] of “ligand”.itp file) will be used.

Just a comment [pairs\_nb] is intended for free-energy simulations (see [http://manual.gromacs.org/current/reference-manual/topologies/molecule-definition.html#pairinteractions](http://manual.gromacs.org/current/reference-manual/topologies/molecule-definition.html#pairinteractions)). An alternative is [pairs].

Best regards  
Alessandra

---

<div class="post-metadata">

### Author: ![Qinghua\_Liao](https://avatars.discourse-cdn.com/v4/letter/q/46a35a/32.png) [@Qinghua\_Liao](https://gromacs.bioexcel.eu/u/Qinghua_Liao)
#### Post date: [December 10, 2020, 11:49am UTC](https://gromacs.bioexcel.eu/t/1-4-interactions/1218/3 "2020-12-10T11:49:52Z")

</div>

Hello Alessandra,

Thanks a lot for your confirmation on it!  
I guess [pairs\_nb] is equivalent to [pairs], if we know what we intend to do, right?

All the best,  
Qinghua

---

<div class="post-metadata">

### Author: ![alevilla](https://dub1.discourse-cdn.com/flex017/user_avatar/gromacs.bioexcel.eu/alevilla/32/439_2.png) [@alevilla](https://gromacs.bioexcel.eu/u/alevilla)
#### Post date: [December 11, 2020, 8:43am UTC](https://gromacs.bioexcel.eu/t/1-4-interactions/1218/4 "2020-12-11T08:43:50Z")

</div>

Yes, knowing that you have to define 4 parameters  
Best regards  
Alessandra

---

<div class="post-metadata">

### Author: ![Qinghua\_Liao](https://avatars.discourse-cdn.com/v4/letter/q/46a35a/32.png) [@Qinghua\_Liao](https://gromacs.bioexcel.eu/u/Qinghua_Liao)
#### Post date: [December 11, 2020, 9:16am UTC](https://gromacs.bioexcel.eu/t/1-4-interactions/1218/5 "2020-12-11T09:16:23Z")

</div>

Thanks a lot Alessandra!
