Messages

1

MD simulations / Re: Stacking of two neighboring non-base-paired bases

« on: February 19, 2013, 06:30:38 pm »

Dear Dr. Liu,

Thank you very much for your response. I'm trying to put find_pair in a for loop to write out separate outputs for each of my pdb files. find_pair seems to be working fine, but my outputs aren't being written. I was wondering if there was a bug that you knew of that wouldn't allow me to use find_pair in such a fashion. I realize that my questions may be outside the scope of this forum; I just wanted to bring it up in case this was an unknown bug with find_pair. Please do let me know if this is the case, and I very much appreciate all your help.

Thanks again, and have a great day.

2

MD simulations / Re: Stacking of two neighboring non-base-paired bases

« on: February 17, 2013, 10:36:06 am »

Dear Dr. Liu,

I have a follow-up question on this discussion: I've been able to use the protocol that you suggested; however, I'm wondering if you have a suggestion as to how to optimize it for analyzing each frame of my trajectory. I have about 25,004 frames in my MD trajectory, and I would have to run find_pair -s for each of those frames, modify by hand each of those outputs and then run analyze on each of those outputs. Is this possible?

Thank you so much, and have a great day.

3

MD simulations / Re: Stacking of two neighboring non-base-paired bases

« on: January 28, 2013, 04:32:06 pm »

Thank you so much, Dr. Lu! I really appreciate it, as well as all your help!

Have a great day!

4

MD simulations / Re: Stacking of two neighboring non-base-paired bases

« on: January 28, 2013, 07:27:53 am »

Dear Dr. Lu,

Since I am conducting this analysis in preparation for a manuscript, would you mind removing the image of the whole structure that you posted, in case a competitor might recognize it? I may be being overly cautious, but I would prefer that to the alternative.

Thank you for your consideration, and have a great day.

5

MD simulations / Re: Stacking of two neighboring non-base-paired bases

« on: January 28, 2013, 07:21:15 am »

Dear Dr. Lu,

Thank you so much for your response. I did think that perhaps my request would be outside x3DNA's normal usage, and that some parameters, like the "Rise" parameter are meaningless. It does seem to me, however, that the slide parameter might still have meaning, as you suggest.

Thank you very much again, and have a great day.

6

MD simulations / Re: Stacking of two neighboring non-base-paired bases

« on: January 27, 2013, 12:25:43 pm »

Dear Dr. Lu,

Thank you very much for the details on x3dna_ensemble; I'm just wondering how one would specify the pdb files of the ensemble? Or if one runs x3dna_ensemble analyze and x3dna_ensemble extract in the same directory as the pdb files to be analyzed, will x3dna_ensemble just run through all the pdb files available? I'm sorry for these questions; I've looked in the x3dna-v2.1/examples/ensemble/md directory, but haven't found a README file, and so am a bit confused about the syntax. Please don't hesitate to direct me to other threads or other sources if they exist and I haven't found them yet.

I'm sorry for not attaching the PDB file before; I only just now understood what you meant. I've attached it to this post.

I was wondering if you could provide more information on how one could manually modify the output from find_pair? As I said, the G1 and U24 bases are not in base pairs themselves, so I'm not sure if adding them to the find_pair output would help.

Thank you very much again, and sorry for more questions.

7

MD simulations / Re: Stacking of two neighboring non-base-paired bases

« on: January 27, 2013, 06:55:10 am »

Dear Dr. Lu,

Thank you for your response.

I want to analyze two bases in a ribozyme, G1 and U24, neither of which are in any base pairs. I want to track the sliding interaction of these two bases across the total length of my trajectory, so that I can correlate the sliding to various other distances and angles that I measure across the length of the trajectory as well. I can see how one can extract the slide parameter for a base pair step, but I don't understand if and how one can extract the slide parameter for two unpaired bases.

First I run find_pair as my first step, and you can see that neither G1 or U24 is listed in the output, which I have attached.

Next, if I use the 'extract' option for x3dna_ensemble, then if my understanding is correct, I would first need to run 'analyze' whose input I believe is from find_pair. I don't think this will work for me because my two bases of interest aren't listed inthe find_pair output.

I also see that there is an 'analyze' option for x3dna_ensemble. I am not sure if this will work for me, but I haven't found the correct syntax to enter my range of PDB files (I have generated one PDB file per frame of the trajectory that I wan't to analyze.). I've tried each of the following 3 commands:

> x3dna_ensemble analyze snaps/10000/10/*
> x3dna_ensemble -v analyze snaps/10000/10/*
> x3dna_ensemble -h analyze snaps/10000/10/*

where "snaps/10000/10/" is the directory that contains the first 10 pdb files from the first 10 frames of my trajectory.

And this is the message I get every time:

> Specify only one of: ensemble|models|pattern|list|one

Please let me know if you need further detail, Dr. Lu, and thank you so much again!

8

MD simulations / Stacking of two neighboring non-base-paired bases

« on: January 26, 2013, 05:54:53 pm »

Dear Dr. Lu and X3DNA users,

I am relatively new to X3DNA, and am trying to use the x3dna_md.rb and extract_par.rb to analyze stacking and sliding of two neighboring bases in my RNA MD simulations. I've successfully used x3dna_md.rb on PDB files of each frame of one of my MD simulations, but I'm thinking this might not be the best way to analyze my data.  I believe that since the two bases I'm interesting in monitoring are not in base pairs themselves, they would be missed by find_pair. I was thinking that the analyze script might be a way to go, but if I'm not mistaken, analyze also takes the output of find_pair for its input.

Any help that can be provided would be much appreciated! Please also do let me know if there are threads that have already addressed my concern; I have looked for these but may have missed them.