Messages

1

MD simulations / Re: Unnatural base pair helical parameters

« on: October 03, 2020, 07:50:35 am »

Hi

Finally succeed in calculating helical base parameters for unnatural base pairs.
I used following procedure
1) converted gromacs trajectory to amber trajectory format using Amber CPPTRAJ
2)  Used modified cpptraj for unnatural base pairs
https://github.com/Amber-MD/cpptraj/pull/806

Thanks
Regards
Shaikh

2

MD simulations / Re: Unnatural base pair helical parameters

« on: September 10, 2020, 11:57:23 am »

Hi

I used x3dna-dssr to get json file using following command. Pdb conatins 100 frames.

x3dna-dssr -i=bdna1-100.pdb -o=bdna.json --more --nmr --json

Once I got json file, I tried to use Markus's R script to extract the data. I installed dplyr and tidyjson library required for the script.
http://forum.x3dna.org/rna-structures/dssr-analyzing-nmr-structures-overwritten-output-files/

However, I am not getting any output.
R -s -f dssr_json_example-v3.R

Following is the output
Attaching package: ‘dplyr’

The following objects are masked from ‘package:stats’:

    filter, lag

The following objects are masked from ‘package:base’:

    intersect, setdiff, setequal, union


Attaching package: ‘tidyjson’

The following object is masked from ‘package:stats’:

    filter

Warning message:
In gather_array(.) :
  array.index column name already exists, changing to array.index.2
# A tbl_json: 0 x 17 tibble with a "JSON" attribute
# … with 17 variables: ..JSON <chr>, document.id <int>, array.index <int>,
#   array.index.2 <int>, chain <chr>, id <chr>, resi <chr>, resn <chr>,
#   alpha <dbl>, beta <dbl>, gamma <dbl>, delta <dbl>, epsilon <dbl>,
#   zeta <dbl>, chi <dbl>, baseSugar_conf <chr>, sugar_class <chr>

I don't know what is wrong. I am new to R language.
Any suggestion from your side.
Thank you
A Shaikh

3

MD simulations / Re: Unnatural base pair helical parameters

« on: September 10, 2020, 05:23:05 am »

Dear Dr. Xiang-Jun
Thank you for reply.
Can we apply this method to gromacs trajectory or amber trajectory. Or can we use multiple pdb frames to calculate helical parameters.
Thanks for DSSR, I have installed the program, going through the manual.

Thanks
Regards
Shaikh

4

MD simulations / Unnatural base pair helical parameters

« on: September 09, 2020, 05:46:09 pm »

Dear Dr. Xiang,

I did a molecular dynamics simulation for DNA with modified base pairs using gromacs. I wanted to analyze base pair parameters for each base pair.

I tried to use 3DNA also, but it also gave errors.
find_pair -p onesb_frame1.pdb allbp_list

handling file <onesb_frame1.pdb>
Match 'IUG' to 'a' for residue IUG    6  on chain   [#6]
    check it & consider to add line 'IUG     a' to file <baselist.dat>
Match 'JSP' to 'u' for residue JSP   19  on chain   [#19]
    check it & consider to add line 'JSP     u' to file <baselist.dat>
missing ' N1 ' atom : residue name 'JSP', chain  , number [  19 ]
missing ' C5 ' atom : residue name 'JSP', chain  , number [  19 ]
Cannot identify RN9/YN1 in [JSP] -- use atom [ C1 ] instead
missing ' N1 ' atom in base ring atoms
missing ' C5 ' atom in base ring atoms
missing ' N1 ' atom : residue name 'JSP', chain  , number [  19 ]


Do let me know how to calculate base pair parameters including major and minor groove parameters. There was one example in your homepage, but its not working in my case.

Thank you very much
Regards
A Shaikh