From: Mehdi Bagherpour (mehdi.bpour_at_gmail.com)
Date: Tue Jun 13 2017 - 12:30:29 CDT
Dear NAMD users,
I am using the ABF simulation for protein-DNA system in NAMD. Actually I
have dimer-protein that are binding to DNA and I want to pull first monomer
of protein vertically in order to get some idea about its binding affinity
to DNA. My selected colvar is the distance between COM of first monomer and
COM of its related interacted DNA. I have use *distanceZ* option as below:
The total mass of protein monomer = 8105.1
The total mass of related DNA = 4309.73
width 0.1 # 0.1 Angstrom of estimated fluctuation width
# difference of two distances
axis (-13.28, -6.63, -6.3)
# atoms of first monomer of protein
# atoms of first related DNA
In simulation at first protein monomer is departing a bit from DNA (around
2 Angstrom) and then instead of protein monomer, DNA is staring to bend in
order to be apart of protein.
In fact this is not the case and I want to see the just effect of protein
leaving DNA without any force on DNA.
My question is: *Do you have any idea what can I do in ABF without putting
any external constrains on DNA in order to prevent DNA bending?*
Looking forward to hearing your ideas.
This archive was generated by hypermail 2.1.6 : Mon Dec 31 2018 - 23:20:22 CST