You are here

Steric clashes in RNA minimization

6 posts / 0 new
Last post
Steric clashes in RNA minimization
#1

I have run into an intermittent problem with RNA_minimizer within the rna_denovo application. I am using identical setups in different directories and shells. I have two calculations currently running without problems. When I try the same calculation in a 3rd directory, I am getting high scores for lk_nonpolar and the resulting pdb structures have steric clashes with overlap side chain atoms. There are no errors reported during the calculation, but the results are not usable.

The calculation is a simple rna hairpin using rna_helix to generate the stem.pdb file:

rna_denovo.macosclangrelease -nstruct 100 -params_file motif.params -fasta motif.fasta -out:file:silent motif.out -include_neighbor_base_stacks -minimize_rna true -s stem.pdb -input_res 6-12 18-24 -cycles 5000 -close_loops -close_loops_after_each_move -filter_lores_base_pairs_early -autofilter -filter_chain_closure_halfway -output_res_num 1-24

Note, I am also having problems with rna_minimize alone. Either the molecules explode or the minimization creates steric clashes between base pairs that are defined in the params file. The rna_minimize application behaves as if it is doing 4D annealing rather than a simple gradient based minimization. Is there anyway to increase the penalty for vdw to prevent steric clashes?

I am using the latest build of Rosetta (2014.35.57232) compiled with clang version 6.0 (600.0.54) in Mac OS X 10.10.0 (currently running 10.10.1). I am preparing to test this in linux next.

I have just tested this in linux (Ubuntu 14.04 kernel 3.13.0-39) with Rosetta compiled under gcc 4.8. I did two runs in different directories with the same setup and both ran fine without the obscenely large lk_nonpolar penalties during the RNA_minimizer step in rna_denovo.

Any suggestions or similar observations are welcome. Thanks. Brian

Category: 
Post Situation: 
Tue, 2014-11-18 18:29
brianmlee

I emailed the person responsible for the rna_denovo code (Rhiju Das). He said that it looks like it's related to a bug he recently fixed. The weekly release 2014_46_57428 or later should include the bugfixes.

Fri, 2014-11-21 09:36
rmoretti

How can I access the newer weekly releases (e.g. 2014.46)? The latest build that appears on the website is 2014.35. Thanks again. Brian

Fri, 2014-11-21 14:35
brianmlee

It might be pre-release and working it's way through. It has to pass certain quality control measures before it makes it to public release stage. (Which that particular release might not.)

When a weekly later than 2014.46 does make it through QA, it will have the bug fixes in question.

Mon, 2014-11-24 15:36
rmoretti

Thanks again. Brian

Mon, 2014-11-24 19:40
brianmlee

Thank you. Brian

Fri, 2014-11-21 14:34
brianmlee