Messages

1

General discussions (Q&As) / Re: basepair centers

« on: April 29, 2020, 10:59:52 am »

Hi Dr. Xiang-Jun,

I'm rather used to modeling cows as discs in a 2d space  , so I can't say I'm uncomfortable with this simplification of the DNA, but I see your point.

Thank you for taking the time to help me. Nice tool by the way, it's been incredibly helpful for analyzing the results of my simulations.

Best regards,

Daniel

2

General discussions (Q&As) / Re: basepair centers

« on: April 28, 2020, 11:22:38 pm »

Hi Dr. Xian-Jun,

thank you for your response.

I think I've finally understood how to obtain the information I'm looking for, which is the origin of each basepair. I am attaching one snapshot of the simulations I've performed in pdb (and the result of find_pair). In this case, what I'm interpreting as the end-to-end distance corresponds to the distance between the origin of the A-T base pair at (166.994, 190.835, 4.669) and the G-C basepair at (146.156, 129.39, 144.990), does this make sense? Could this be interpreted as the center of a basepair?

Best regards,

Daniel

3

General discussions (Q&As) / basepair centers

« on: April 28, 2020, 06:10:00 pm »

Hi everyone,


I hope you are all doing ok on these crazy days.

I have performed some molecular dynamics simulations for a DNA-protein system. I am trying to understand the results by analyzing the DNA's structure before and after the protein binds. In particular I'm trying to calculate a parameter called the crookedness (beta) of the DNA which is defined as cosine(beta) = x / y , where x is the end-to-end distance of the molecule and y is the sum of the individual base-pair-center distance (PHYSICAL REVIEW LETTERS 122, 048102 (2019)). Reading the x3dna documentation I've found that the distance between two consecutive base-pair centers is equal to the square root of (shift^2 + slide² + rise^2), nevertheless I still need to calculate the end-to-end distance of the entire DNA (or the centers for the corresponding basepairs). Can anyone please help me to obtain this information with x3dna? I have already run x3dna_ensemble analyze for a list of MD snapshots and can run x3dna_ensemble extract stuff but do not know where to find the information for the coordinates of the basepairs.

Best regards.

Daniel