[AMBER] newbie question: simulating miRNA

From: Michiel de Hoon <mdehoon.gsc.riken.jp>
Date: Fri, 31 Aug 2012 23:30:25 +0900

Dear all,

I am trying to simulate the dynamics of a mature microRNA of 22
nucleotides, and I am running into a problem creating the initial structure.

The sequence I want to simulate is "uagcuuaucagacugauguuga".
To create a PDB file with this sequence, I use this .nab file:

molecule m;

m = fd_helix("arna", "auagcuuaucagacugauguuga", "rna");
putpdb( "mirna.pdb", m, "-wwpdb");

I run nab on it, and I run the resulting executable. This gives me the
mirna.pdb file.
Since this PDB file contains a double helix of RNA and I am interested
in only a single strand, I edit the PDB file to remove the second strand.

The sequence I give to fd_helix has one additional nucleotide "a" at the
5' side. This is because I noticed that fd_helix creates an RNA with a
5' hydroxyl group instead of a 5' phosphate group. So I added this dummy
nucleotide, I remove all its atoms except its final oxygen from the PDB
file, and associate the last oxygen with the subsequent "u" nucleotide.
My PDB file then looks as follows:

ATOM 1 OP3 U 1 7.297 -6.145 5.250 1.00
0.00 O
ATOM 2 P U 1 6.882 -5.338 6.560 1.00
0.00 P
ATOM 3 OP1 U 1 7.505 -5.970 7.750 1.00
0.00 O
ATOM 4 OP2 U 1 5.404 -5.201 6.610 1.00
0.00 O
ATOM 5 O5' U 1 7.538 -3.907 6.290 1.00
0.00 O
ATOM 6 C5' U 1 8.968 -3.825 6.160 1.00
0.00 C
ATOM 7 H5' U 1 9.304 -4.443 5.327 1.00
0.00 H
ATOM 8 H5'' U 1 9.430 -4.176 7.083 1.00
0.00 H
ATOM 9 C4' U 1 9.384 -2.375 5.910 1.00
0.00 C
ATOM 10 H4' U 1 10.470 -2.299 5.860 1.00
0.00 H
...
ATOM 698 C2' A 22 7.329 -4.907 64.990 1.00
0.00 C
ATOM 699 H2' A 22 7.159 -4.004 65.576 1.00
0.00 H
ATOM 700 O2' A 22 8.718 -5.197 64.960 1.00
0.00 O
ATOM 701 HO2' A 22 9.154 -4.443 64.556 1.00
0.00 H
ATOM 702 O3' A 22 7.233 -6.221 67.070 1.00
0.00 O
ATOM 703 HO3' A 22 7.285 -5.486 67.686 1.00
0.00 H
TER


Then I start xleap:

xleap -s -f $AMBERHOME/dat/leap/cmd/leaprc.ff99SB

I don't know if ff99SB is the appropriate force field to use to simulate
RNA, but that is what I used.

Inside xleap I then do:

> rna1 = loadpdb "mirna.pdb"
Loading PDB file: ./mirna.pdb
Created a new atom named: OP3 within residue: .R<RU5 1>
Created a new atom named: P within residue: .R<RU5 1>
Created a new atom named: OP1 within residue: .R<RU5 1>
Created a new atom named: OP2 within residue: .R<RU5 1>
  total atoms in file: 703
  Leap added 1 missing atom according to residue templates:
       1 H / lone pairs
  The file contained 4 atoms not in residue templates

Next I add sodium ions for neutralization:
> addions rna1 Na+ 0
21 Na+ ions required to neutralize
Adding 21 counter ions to "rna1" using 1A grid
...

So this is my problem: Only 21 counter ions are added, while the
microRNA contains 22 nucleotides. The output from loadpdb suggests that
the 5' phospate group is not being recognized, which may be the reason
why xleap only adds 21 sodium ions.
What is then the right way to create an RNA sequence with a phosphate on
its 5' side and a hydroxyl group on its 3' side?

Many thanks in advance,
-Michiel


_______________________________________________
AMBER mailing list
AMBER.ambermd.org
http://lists.ambermd.org/mailman/listinfo/amber
Received on Fri Aug 31 2012 - 08:00:03 PDT
Custom Search