Messages

1526

General discussions (Q&As) /

« on: July 08, 2007, 10:51:31 am »

Good question!

The utility "blocview" was designed to generates a schematic image which combines base block representation with protein ribbon. The image has informative color coding for the nucleic acid part and is set in the "best-view" by default.

For the nucleic acid part, the "primary" object is the planar base rings which are simplified by corresponding rectangular blocks. Base-pairs are only secondary objects, abstracted at a higher level. Moreover, except for normal Watson-Crick base-pairs, there is an issue of how to represent the "faces" of a pair. Also, think of triplets and higher order base associations. Those are the rationales that "blocview" generates only a block representation for each individual base.

I did have the bp block issue in mind when I wrote 3DNA, and it is possible to generate bp blocks for normal Watson-Crick base-pairs with 3DNA v1.5. Using bdl084 as an example, try the following steps:

Code: [Select]

   blocview bdl084.pdb
    cp t.pdb bdl084_bv.pdb
    find_pair bdl084_bv.pdb stdout | analyze

    frame_mol -s ref_frames.dat bdl084_bv.alc

    alc2img -r -n bdl084_bv.alc bdl084_bp.r3d
    cat t.r3d bdl084_bp.r3d > bdl084_all.r3d
    render < bdl084_all.r3d > bdl084.avs

    convert -trim +repage -border 5 -bordercolor white bdl084.avs bdl084.jpg


The generated image is attached below. Knowing what you are doing and taking 3DNA as a toolkit, you should be able to get what you want in a script. A possible simple starting point is by modifying the "blocview" Perl script. You might want to have a try and possibly contribute back once you are happy with your script. Do not forget that there is a Users' Corner section in 3DNA homepage. With enough interest,  I might consider to add a new script in 3DNA v2.0

As a side note, the ".r3d" file "bdl084_all.r3d" can be fed into Pymol directly to generate photo-realistic quality images. Please also refer to another post in this forum for further example images. Of course, you might be aware that "blocview" generated images are used in the NDB, and as a nice surprise (to me), now also  used in the PDB!

HTH,

Xiang-Jun

1527

General discussions (Q&As) /

« on: July 01, 2007, 10:58:56 am »

Thanks for using 3DNA.

Firstly, a detour: you mentioned the "o1p_o2p" utility program which was designed to check for correct labelling of the O1P and O2P atoms of nucleotides in a PDB file. The utility was initially written when I noticed the incorrect O1P/O2P labelling problem in the NDB entry adh008 (9dna) used as a test case in the 3DNA distribution.

To my surprise, even as of today, the problem is still there in the NDB (adh008) and PDB (9dna). To make it concrete, here is the except from corresponding data files in PDB format:

NDB Asymmetric Unit coordinates

[pdb format, Unix compressed(.Z): pdb9dna.ent]:

Code: [Select]

ATOM     24  P     C A   2      35.309  46.286  11.717  1.00 28.71      9DNA  87
ATOM     25  O1P   C A   2      35.789  44.945  11.170  1.00 26.94      9DNA  88
ATOM     26  O2P   C A   2      36.213  47.108  12.492  1.00 25.98      9DNA  89[/list]

NDB Biological Unit coordinates

[pdb format: 9dna.pdb1]:

Code: [Select]

ATOM     20  P     C A   2      35.309  46.286  11.717  1.00 28.71           P  
ATOM     37  OP1   C A   2      36.213  47.108  12.492  1.00 25.98           O  
ATOM     38  OP2   C A   2      35.789  44.945  11.170  1.00 26.94           O  [/list]

PDB coordinates [

9DNA.pdb]:

Code: [Select]

ATOM     24  P     C A   2      35.309  46.286  11.717  1.00 28.71      9DNA  87
ATOM     25  O1P   C A   2      35.789  44.945  11.170  1.00 26.94      9DNA  88
ATOM     26  O2P   C A   2      36.213  47.108  12.492  1.00 25.98      9DNA  89[/list]
Have you noticed the inconsistencies?

I am impressed you cared to run "o1p_o2p" to check for the correct O1P/O2P labelling in your interested DNA structures. From my limited experience, I feel those days people are much more interested in solving problems of "grand importance". Yet, whenever checked seriously, even "simple" issues claimed to have been resolved for so many years still have left so many details to be worked out. From SCHNAaP and SCHNArP, to resolving the discrepancies, to 3DNA and its continuous support/further refinements over the years (in my spare time), I have been trying to solve one specific problem that apparently has made a much bigger impact than I originally expected ten years ago.

Now, to answer your question. As of 3DNA v1.5, the auxiliary.par file does not contain the O1P (and O2P for that matter) xyz coordinates w.r.t. the middle frame. I am getting that info into v2.0 which hopefully could be released in the near future.

In the meantime, please check the file "stacking.pdb" which contains all the info you asked and much more. Of course, you need to have some (script) programming skills to extract the info you need and be careful about other technical details such as residue numbering, and direction of middle frame etc. I would say it's well worth the efforts ...

HTH,

Xiang-Jun

1528

General discussions (Q&As) /

« on: June 15, 2007, 10:00:19 pm »

Firstly, please refer to a related post for background information.

More specifically, your DNA structures are too irregular to be recognized by "find_pair" as you would expect. The output file from "analyze" follows whatever base-pair information produced by "find_pair", and thus the "problems" your noticed. You could run "find_pair" to create an approximate starting input file and modify it accordingly as you see fit. The file format produced by "find_pair" is documented at the 3DNA homepage.

Also note that you can now add attachments to your post in this forum.

HTH,

Xiang-Jun

1529

General discussions (Q&As) /

« on: June 15, 2007, 09:49:09 pm »

As explained in my reply to your original question, what you now propose to do is perfectly fine. It is actually the preferred way for MD trajectories and NMR ensembles, where the DNA structure changes only in conformation not in its base-pairing. You can thus start with a representative structure to derive the correct pairing information with "find_pair" (and make manual changes where necessary), and you then need only to change the first two lines (PDB file and corresponding output file) of the input file to "analyze/cehs". Of course, this is best done automatically with a script.

The program "find_pair" makes analyzing nucleic acid structure so easy that people normally use it together with "analyze/cehs" as an integrated unit. As a matter of fact, "find_pair" has more functionalities (as another application, see SWS from Pascal Auffinger) than just providing input to the analysis routines, and it is best to think them as two separate entities.

Based on my experience, "find_pair" has been working really great for what it was designed for. Over the years, it has been refined and thus become more robust and sophisticated in the coming new release of 3DNA v2.0. At this time, I feel confident to say that whenever "find_pair" (and accordingly "analyze/cehs") produces some weird base-pairing results, it is because the nucleic acid structure is too distorted in the corresponding region(s). Whenever in doubt, it certainly helps to have a look of your structure using a molecular viewing program such as Rasmol/Pymol.

As always, we welcome users' feedback for specific cases where the algorithm could possibly be improved. That's why I asked you to send me some structures you felt that "find_pair" was not working as expected. To make life easier for future communications, I have installed the phpBB2 "Attachment Mod" for this forum. Please make good use of it!

HTH,

Xiang-Jun

1530

General discussions (Q&As) /

« on: June 13, 2007, 09:28:10 pm »

I do not know if I understand exactly what you mean. Why not post a specific (minimal) sample case to illustrate your point?

Also, I find it hard to read a sentence such as "its showng jst 3 residues and in some 4 nd lk wise." -- Maybe that's because English is not my native language? In some sense,   this post might help.

Thanks,

Xiang-Jun

1531

General discussions (Q&As) /

« on: May 16, 2007, 10:34:04 pm »

Thanks for using 3DNA! I am glad that you solved your two previous problems by browsing the forum -- that's exactly why the forum was created in the fist place.

While 3DNA is not a point-and-click package, setting it up should be straightforward and requires minimum effort, i.e., by setting the environment variable X3DNA and put $X3DNA/bin into your command line search path. To make life even easier for inexperienced Linux/Unix users, I have created a Perl script (setup_x3dna) to detect the SHELL and set up these two items automatically -- it will be included in the coming new release.

Now come to your question. The main output file name is derived from your PDB file by changing to (or appending) the ".out" extension. It will be in the same directory as your PDB structure file. Have a look at .inp file -- the 2nd line will tell you where the .out resides. Now I notice a bit of inconsistency here -- all other files are written in the current working directory (CWD), but this .out file is not. Maybe it is more logical to put the .out file in the CWD as well. I will make this change into the new release of 3DNA.

HTH,

Xiang-Jun

1532

General discussions (Q&As) /

« on: May 08, 2007, 10:06:40 pm »

Hi Whitney,

I deliberately skipped the details on calculating bp-reference frame when writing that part of the manual for two reasons: (1) it is identical to the procedure used for the step "middle-frame"; (2) it gives serious users like you an opportunity to figure it out (normal users won't bother tech part anyway). In my experience, this is the best way to have a clear understanding of the mathematics under 3DNA -- which is actually quite simple.

For your verification, the bp frame origin is the direct average of the two bases forming the pair. The bp z-axis is actually the averaged and normalized vector of the two base z-axes. This is not the case for bp x- and y-axes.

Code: [Select]

AG1_o = [15.1632 -0.0362 -4.4678]
AG1_z = [-0.4004 0.4012 0.8238]

BC8_o = [14.9124 0.2803 -4.7498]
BC8_z = [-0.3422 0.5326 0.7741]

bp_o = mean([AG1_o; BC8_o])
         % ====> 15.03780    0.12205   -4.60880

zave = mean([AG1_z; BC8_z])
         % =====> -0.37130   0.46690   0.79895
bp_z = zave / norm(zave)
         % =====> -0.37239   0.46826   0.80128

stagger = dot(AG1_o - BC8_o, bp_z)
         % =====> -0.015638


I am using Matlab (octave) syntax for illustration. In 3DNA, the bp z-axis is calculated along with bp x- and y-axes after symmetrically rotating the two base frames along the buckle-propeller hinge axis. Please have a try and hopefully you would post back step-by-step how you proceed to clarify your own understanding and for the benefit of other 3DNA users.

HTH,

Xiang-Jun

1533

General discussions (Q&As) /

« on: April 16, 2007, 09:50:13 pm »

"find_pair" generates an input file in text format that you can edit to suit your specific need before feeding into "analyze" (or "cehs").

Just imagine the situation without access to "find_pair": how would you prepare (from a nucleic acid containing structure in PDB format) an input file for an analysis program such as Curves etc?

In the 3DNA homepage, I have also maintained a list of citations to 3DNA compiled from Web of Science. It is not  necessarily complete, but could serve as a good starting point to related published work.

So far, I've been the only person in answering questions in the forum. That's fine as a starting point. In the long run, I do hope other colleagues would get more actively involved: not only in asking questions, but also in answering them; moreover, in sharing your experiences, or contributing scripts/programs etc. Your contribution can make a difference!

HTH,

Xiang-Jun

1534

General discussions (Q&As) /

« on: April 13, 2007, 10:02:55 pm »

Hi,

Please check carefully the PDB file format. As I understand it, REMARK 210 contains info about which one is the best representative structure in an NMR ensemble, which is used by "ex_str". If no such info is available, or the PDB is somehow inconsistent for certain entry, the first structure is extracted. That's exactly what my previous message is about. It also contains some example PDB entries that you can verify. These were the samples I tested when adding such functionality in "ex_str".

A better understanding of the PDB format would certainly help. In any event, you can always check the PDB file itself to see what's in there -- it is a text file anyway.

Hope this clarifies the issue a little bit.

Xiang-Jun

1535

General discussions (Q&As) /

« on: April 12, 2007, 09:16:19 pm »

Hi,

Please down the 3DNA v1.5 from its homepage. The Linux version, as I checked it just a moment ago, has the content of file "misc_3dna.par" as follows:

Code: [Select]

4.0 0.0 ON A1 # upper H-bond length limits/atoms, alternative location
15.0          # max. distance between paired base origins (FIND_PAIR)
2.5           # max. vertical distance between paired base origins (FIND_PAIR)
65.0          # max. angle between paired bases [0-90] (FIND_PAIR)
4.5           # MIN. distance between RN9/YN1 atoms (FIND_PAIR)
7.5           # max. distance criterion for helix break (FIND_PAIR)
0.6           # r7 criterion to decide if a helix is curved
3.2           # r8 H-bond distance for water molecules with base N/O atoms
4.5           # r9 maximum O3'--P distance for linkage (REBUILD)


So it seems you are using an earlier version.

HTH,

Xiang-Jun

1536

General discussions (Q&As) /

« on: April 12, 2007, 12:21:05 am »

Hi,

There is a default cut-off threshold (0.6) used in 3DNA to decide if a structure is too curved to calculate the global parameters based on C1'-C1' vector. It is defined in the file "misc_3dna.par" (under directory $X3DNA/BASEPARS) on the 7th line -- increase the value (based on your data) will make 3DNA output the parameters you request.

Needless to say, the file "misc_3dna.par" is a bit cryptic -- the parameters were added while 3DNA was under actively development. In the coming new release of 3DNA v2.0, the file is much easer to understand, as show below.

Code: [Select]

# Section 4: is this double helix curved?
std_curved: 0.6         # criterion to decide if a helix is curved (0.6)

HTH,

Xiang-Jun

1537

General discussions (Q&As) / Re: TA-DNA dataset

« on: April 10, 2007, 08:56:56 pm »

This question has popped up a few times. Now with this forum, I do not need to repeat myself, but refer you to the thread "DNA standards/statistics using 3DNA", which contains a detailed list of the TA-DNA dataset used in the 2003 NAR 3DNA paper.

HTH,

Xiang-Jun

1538

General discussions (Q&As) /

« on: April 09, 2007, 11:38:36 pm »

Hi,

The following is taken directly from the C-source code comment:

Code: [Select]

/* Best representative NMR conformer, use #1 if not assigned in REMARK 210.
 * Typical cases:
 *    BEST REPRESENTATIVE CONFORMER IN THIS ENSEMBLE : NULL [pdb28sp]
 *    BEST REPRESENTATIVE CONFORMER IN THIS ENSEMBLE : 1    [pdb2bj2]
 *    BEST REPRESENTATIVE CONFORMER IN THIS ENSEMBLE : 3    [pdb2lef]
 *    BEST REPRESENTATIVE CONFORMER IN THIS ENSEMBLE : 5    [pdb3php]
 *
 * For structures pdb1qby, pdb1djd & pdb1bjd, the REMARK 210 says the
 *    best representative conformer is #12 or #13, but the data file
 *    actually does not has that model. Reset inum to 1 for such cases.
 */
HTH,

Xiang-Jun

1539

General discussions (Q&As) /

« on: April 05, 2007, 08:52:04 pm »

It is not an error, but for information only. It means that in your nucleic acid structure there is some O3' to P[i+1] linkage that's out of normal range (> 2.5 A is the cut off in 3DNA).  This message is associated with "analyze" instead of "find_pair". For example, it would happen with the "-s" option unless your nucleic acid structure is single strand and continuous (e.g., find_pair -s bdl084.pdb stdout | analyze). Please also refer to the thread on missing P atom in first residues.  

If anything is still unclear, please provide a repeatable example (structure) with a step-by-step description what you did and what you found.

HTH,

Xiang-Jun

1540

General discussions (Q&As) /

« on: March 03, 2007, 12:14:50 am »

Hi,

Welcome to the 3DNA community!

The program "fiber" is your friend to build long regular DNA models. Also check the topic on PDBML in this forum.

HTH,

Xiang-Jun

1541

General discussions (Q&As) /

« on: February 24, 2007, 11:37:56 am »

Then it is clear that you do not have permission set up correctly for $X3DNA/bin executables. I do not know your level of Linux/Unix, and/or access right to the /usr/local directory. But you can try the following:

Code: [Select]

export X3DNA=/usr/local/X3DNA
export PATH=$PATH:$X3DNA/bin
chmod a+rx /usr/local/X3DNA/bin/*
find_pair


Other than that, you need to consult your system administrator or local Linux/Unix expert for help.

Alternatively, you can always install 3DNA under your home directory.

HTH,

Xiang-Jun

1542

General discussions (Q&As) /

« on: February 22, 2007, 10:42:38 pm »

Hi,

Life should not be that hard! Setting up the 3DNA environment has rarely been an issue, indeed yours is the first case I've come cross over the past few years.

I noticed that you set X3DNA as

Code: [Select]

export X3DNA=/usr/local/X3DNA/i.e., with a trailing '/'. In the output error message, you have path

Code: [Select]

/usr/local/X3DNA//binNoticed the two '//' between X3DNA and bin? The error message means you have not put 3DNA bin directory into your command line search path yet.

Try the following in your current terminal, one line at a time, what did you get?

Code: [Select]

echo $SHELL
       # do you get '/bin/bash'? -- verifies you have 'bash' shell
ls /usr/local/X3DNA/bin/find_pair
       # this verifies that you indeed installed 3DNA under /usr/local
export X3DNA=/usr/local/X3DNA
       # without trailing '/'
echo $X3DNA
ls $X3DNA/bin/find_pair
       # If 3DNA is installed under /usr/local, this should print the full path
export PATH=$PATH:$X3DNA/bin
which find_pair

Please report back how it goes.

Xiang-Jun

1543

General discussions (Q&As) /

« on: February 19, 2007, 11:24:57 pm »

Yes, there is a way in 3DNA to calculate all backbone torsion angles of any nucleic acid structure in PDB format. It is as follows:

Code: [Select]

find_pair -s -t PDBFILE stdout | analyze

where option "-s" (or equivalently "-1") is used to treat the whole structure  as a continuous single helix,  needed for getting all backbone torsion angles.

The "-t" option also takes consideration of HETATM records for modified bases. Note that in the coming new release of 3DNA v2.0, "-t" will be set by default, thus it will no longer be needed.

Of course, this method assumes that the nucleotides in the PDB structure are in proper order in the 5'-->3' direction. This should normally be the case. If need arises, I will consider add a new functionality in "find_pair" to make the  nucleotides in order no matter whatever their original positionings in the PDB file are.

HTH,

Xiang-Jun

1544

General discussions (Q&As) /

« on: February 13, 2007, 08:25:04 pm »

Hi,

Regarding applying 3DNA to quadruplex, you might find this post in the forum relevant.

If your are only interested in calculate the backbone torsions and pucker of any nucleic acid structure, the following command will do the trick:

Code: [Select]

find_pair -s -t PDB_FILE stdout | analyze

HTH,

Xiang-Jun

1545

General discussions (Q&As) /

« on: February 13, 2007, 08:13:08 pm »

Hi,

Thanks for using 3DNA.

What happens if you do the following?

Code: [Select]

which find_pair
source ~/.bashrc
which find_pair


In other words, I suspect your setting is not in effect yet before running the source command, or starting another shell.

HTH,

Xiang-Jun

1546

General discussions (Q&As) /

« on: December 27, 2006, 10:36:43 pm »

Hi,

This is weird -- it certainly should not happen as far as 3DNA itself is concerned. I suspect that you have another program called "rebuild" in your command line search path, and it is ahead of $3DNA/bin.

To verify, try:

Code: [Select]

which rebuild from your shell prompt. What do you get?

Alternatively, try to run $3DNA/bin/rebuild instead of just "rebuild".

Have a try and please let me know how it goes.

Xiang-Jun

1547

General discussions (Q&As) /

« on: December 11, 2006, 09:51:08 pm »

Hi,

The problem you have is due to the wrong PDB format. Specifically, with the atom name, as shown below:

Code: [Select]

ATOM     11 N1   CYT     1    -278.189 335.614  -8.159                       N
instead of the correct format as in:
ATOM      9  N1    C A   1      18.770  29.711  23.733  1.00  7.31           N

In the coming new release of 3DNA, I have taken such case into consideration, and the input file following "find_pair" would be:

Code: [Select]

iBdnaVacFINAL_x3dna2.pdb
iBdnaVacFINAL_x3dna2.out
    2         # duplex
   10         # number of base-pairs
    1    1    # explicit bp numbering/hetero atoms
    1   21  0 #    1 | ....>-:...1_:[CYT]C-----G[GUA]:..22_:-<....  0.64  0.09 11.06  9.06 -3.63
    2   20  0 #    2 | ....>-:...2_:[CYT]C-----G[GUA]:..21_:-<....  0.65  0.01  9.87  9.04 -3.84
    3   19  0 #    3 | ....>-:...3_:[ADE]A-----T[THY]:..20_:-<....  0.11  0.09 13.97  9.06 -4.00
    4   18  0 #    4 | ....>-:...4_:[THY]T-----A[ADE]:..19_:-<....  0.14  0.08 11.22  9.06 -4.13
    5   17  0 #    5 | ....>-:...5_:[CYT]C-----G[GUA]:..18_:-<....  0.67  0.08  5.34  9.07 -3.89
    7   16  0 #    6 | ....>-:...7_:[CYT]C-----G[GUA]:..16_:-<....  0.70  0.26 16.15  8.99 -2.98
    8   15  0 #    7 | ....>-:...8_:[THY]T-----A[ADE]:..15_:-<....  0.29  0.29 13.17  9.06 -3.49
    9   14  0 #    8 | ....>-:...9_:[ADE]A-----T[THY]:..14_:-<....  0.21  0.01  1.70  9.21 -4.69
   10   13  0 #    9 | ....>-:..10_:[CYT]C-----G[GUA]:..13_:-<....  0.66  0.05 17.87  9.01 -3.34
   11   12  0 #   10 | ....>-:..11_:[CYT]C-----G[GUA]:..12_:-<....  0.60  0.10 12.40  9.06 -3.58
##### Base-pair criteria used:   4.00 15.00  2.50 65.00  4.50  7.50 ON
##### 0 non-Watson-Crick base-pairs, and 1 helix (0 isolated bps)
##### Helix #1 (10): 1 - 10  ***broken O3'[i] to P[i+1] linkage***


For the time being, you could write a script to change the atom name following the PDB convention.

Xiang-Jun

1548

General discussions (Q&As) /

« on: November 27, 2006, 06:48:43 pm »

Hi Pascal,

The "dv" is the vertical distance between the two base planes -- it is actually stagger in magnitude (i.e., abs(stagger)).

Or it can be thought of as the "vertical distance between base pair origins" projected onto the mean normals of the two bases (z-axis of the bp).

The "base pair origin" is actually the geometric average of the two base origins. In 3DNA, starting from the standard base reference frame (Atomic_A.pdb etc), a least-squares fitting onto corresponding experimental base gives the base reference frame (i.e., its position [origin] and orientation [x-, y- and z-axis]).

A detailed step-by-step description of the ls-fitting procedure is provide at the 3DNA homepage technical details section. Please note that over there, the standard base coordinates used is not in the standard reference frame.

In 3DNA, you can get the output of the each base reference frame by:

Code: [Select]

find_pair -s PDB_FILE stdout | analyze
Then have a look of the file "ref_frames.dat".

HTH,

Xiang-Jun

1549

General discussions (Q&As) /

« on: November 24, 2006, 02:52:11 pm »

Hi Pascal,

Nice to communicate with you here.

The standard base reference frame is described in the Olson et al. (2001) JMB report. Specifically, the origin is defined as the center of what would a perfect Watson-Crick base-pair.

As detailed in the FAQ section of the current 3DNA homepage, the default set of parameters used in 3DNA was based on a survey of the NDB structures and is very generous to account for distorted bps. As an example, how much a H-bond distance cut-off should be choosen? Intuitively, 3.2 A would be a more reasonable value than 4.0 A. However, even with such high cut-off, some users still report "find_pair" missing some bps, which turns out that the shorted possible H-bond distance would be over 5.0 A in such bps.

That's exactly the reason that 3DNA provides the "misc_3dna.par" file for users to tune it for their specific purpose. For example, by setting d_org = 1.0, dv = 0.5, z_ang = 15, one will find only Watson-Crick base-pairs.

As far as how stagger parameter is calculated, it is the projection of the vector linking the two base origins onto the mean base z-axis. A step-by-step description of the procedure can be found in the 1997 JMB SCHNAaP paper. 3DNA uses exactly the same algorithm except for a change in reference frame.

HTH,

Xiang-Jun

1550

General discussions (Q&As) /

« on: November 16, 2006, 12:18:24 am »

Hi,

Thanks for sending me the 3 coordinate files via email -- they are helpful in clarifying the issue.

You might be interested to know that in the future, you no longer need to send us email with attachments: I have set up a file upload forum at URL:
http://3dna.rutgers.edu:8080/upload/ so that users can send related information for us to debug and solve their problems. This will also allow for storage and consistent handling of user-cases.

As to your questions, firstly, you can run find_pair on it after adding the following entries to the file baselist.dat, as documented in the FAQs:

Code: [Select]

DG5     G
DG      G
DC      C
MOL     g
DC3     C
Secondly, given the irregularity of the the G-quartet structure, the normal double strand parameters obviously no longer apply. However, you can run the following to get backbone torsion angle etc parameters (of course, watch out the chi torsion related to isoG residue):

Code: [Select]

find_pair -s L1_avstr_min.pdb stdout | analyzeYou can also

Code: [Select]

find_pair -p L1_avstr_min.pdb  multi.infoThe program identifies the two G-quartets. The file multiplets.pdb contained the structures of the multiplets, which you can extract and generate publication quality images, as referred in my previous reply.

Thirdly, with the utility program blocview, run as follows:

Code: [Select]

blocview -i=L1_view.jpg L1_avstr_min.pdb
blocview -i=L1_view_z40.jpg -z=40 L1_avstr_min.pdb

You will get the following two images:


which is the default,  and
,
with a rotation about z-axis by 40 degress to make it "vertical". The color coding the Gs and Cs make them clearly distinguished, and the black minor groove edges of Gs are also obvious, and of course, the stacking ...

The 3 in 3DNA certainly not only stands for 3-dimension, but also the three integrated parts: analysis, rebuilding and visualizaion. Such block view images have been used in the NDB, and PDB. Somehow, it has not been widely adopted by the 3DNA community at large. I am just taking this opportunity to illustrate some of the not-commonly used features.

HTH,

Xiang-Jun