BioNMR
NMR aggregator & online community since 2003
BioNMR    
Learn or help to learn NMR - get free NMR books!
 

Go Back   BioNMR > Software > Molecular modeling software > NMR structure calculation
Advanced Search



Jobs Groups Conferences Literature Pulse sequences Software forums Programs Sample preps Web resources BioNMR issues


Webservers
NMR processing:
MDD
NMR assignment:
Backbone:
Autoassign
MARS
UNIO Match
PINE
Side-chains:
UNIO ATNOS-Ascan
NOEs:
UNIO ATNOS-Candid
UNIO Candid
ASDP
Structure from NMR restraints:
Ab initio:
GeNMR
Cyana
XPLOR-NIH
ASDP
UNIO ATNOS-Candid
UNIO Candid
Fragment-based:
BMRB CS-Rosetta
Rosetta-NMR (Robetta)
Template-based:
GeNMR
I-TASSER
Refinement:
Amber
Structure from chemical shifts:
Fragment-based:
WeNMR CS-Rosetta
BMRB CS-Rosetta
Homology-based:
CS23D
Simshift
Torsion angles from chemical shifts:
Preditor
TALOS
Promega- Proline
Secondary structure from chemical shifts:
CSI (via RCI server)
TALOS
MICS caps, β-turns
d2D
PECAN
Flexibility from chemical shifts:
RCI
Interactions from chemical shifts:
HADDOCK
Chemical shifts re-referencing:
Shiftcor
UNIO Shiftinspector
LACS
CheckShift
RefDB
NMR model quality:
NOEs, other restraints:
PROSESS
PSVS
RPF scores
iCing
Chemical shifts:
PROSESS
CheShift2
Vasco
iCing
RDCs:
DC
Anisofit
Pseudocontact shifts:
Anisofit
Protein geomtery:
Resolution-by-Proxy
PROSESS
What-If
iCing
PSVS
MolProbity
SAVES2 or SAVES4
Vadar
Prosa
ProQ
MetaMQAPII
PSQS
Eval123D
STAN
Ramachandran Plot
Rampage
ERRAT
Verify_3D
Harmony
Quality Control Check
NMR spectrum prediction:
FANDAS
MestReS
V-NMR
Flexibility from structure:
Backbone S2
Methyl S2
B-factor
Molecular dynamics:
Gromacs
Amber
Antechamber
Chemical shifts prediction:
From structure:
Shiftx2
Sparta+
Camshift
CH3shift- Methyl
ArShift- Aromatic
ShiftS
Proshift
PPM
CheShift-2- Cα
From sequence:
Shifty
Camcoil
Poulsen_rc_CS
Disordered proteins:
MAXOCC
Format conversion & validation:
CCPN
From NMR-STAR 3.1
Validate NMR-STAR 3.1
NMR sample preparation:
Protein disorder:
DisMeta
Protein solubility:
camLILA
ccSOL
Camfold
camGroEL
Zyggregator
Isotope labeling:
UPLABEL
Solid-state NMR:
sedNMR


Reply
Thread Tools Search this Thread Rate Thread Display Modes
  #1  
Unread 08-15-2005, 01:51 AM
nmrlearner's Avatar
Senior Member
 
Join Date: Jan 2005
Posts: 23,751
Points: 193,617, Level: 100
Points: 193,617, Level: 100 Points: 193,617, Level: 100 Points: 193,617, Level: 100
Level up: 0%, 0 Points needed
Level up: 0% Level up: 0% Level up: 0%
Activity: 50.7%
Activity: 50.7% Activity: 50.7% Activity: 50.7%
Last Achievements
Award-Showcase
NMR Credits: 0
NMR Points: 193,617
Downloads: 0
Uploads: 0
Default

The following script for CNS simulated annealing was submitted by trent. He recommends to use this script as a part of RECOORD structure determination and water refinement protocol to prevent wrong geometry of aromatic rings (as reported here and here). You may find this script useful if you have a ligand in your model that makes CNS go crazy and produce protein structures with collapsed aromatic rings.

Thanks trent!

Mark


************************************************** *****************

remarks file sa_cns.inp
remarks Simulated annealing protocol for NMR structure determination with cns.
remarks The starting structure for this protocol can be any structure with
remarks a reasonable geometry, such as randomly assigned torsion angles or
remarks extended strands.
remarks based on sa_new.inp from xplor
remarks high_steps and cool_steps increased 4 fold to improve convergence
remarks uses parallhdg5.3.pro and topallhdg5.3.pro from recoord without error
remarks works with additional ligand files
remarks Author: Trent Bjorndahl

{====>}
evaluate ($init_t = 1000 )
{====>}
evaluate ($high_steps= 24000 )
{====>}
evaluate ($cool_steps = 12000 )

parameter
{====>}
@parallhdg5.3.pro
!@yourparameterfilehere.par
end

{====>}
structure @yourstructurefile.mtf end {*Read the structure file.*}

{====>}
coordinates @yourcoordinatefile.pdb {*Read the coordinates.*}


noe
{====>}
nres=10000 {*Estimate greater than the actual number of NOEs.*}
class all
{====>}
@noes.tbl {*Read NOE distance ranges.*}
end

{====>}
restraints dihedral reset
nass = 1000
@dihedrals.tbl
end

flags exclude * include bonds angle impr vdw noe cdih end

{*Friction coefficient for MD heatbath, in 1/ps. *}
do (fbeta=10) (all)
{*Uniform heavy masses to speed molecular dynamics.*}
do (mass=100) (all)

noe {*Parameters for NOE effective energy term.*}
ceiling=1000
averaging * cent
potential * soft
scale * 50.
sqoffset * 0.0
sqconstant * 1.0
sqexponent * 2
soexponent * 1
asymptote * 0.1 {*Initial value--modified later.*}
rswitch * 0.5
end

parameter {*Parameters for the repulsive energy term.*}
nbonds
repel=1. {*Initial value for repel--modified later.*}
rexp=2 irexp=2 rcon=1.
nbxmod=3
wmin=0.01
cutnb=4.5 ctonnb=2.99 ctofnb=3.
tolerance=0.5
end
end

restraints dihedral
scale=5.
end

{====>}
evaluate ($end_count=10) {*Loop through a family of 10 structures.*}

coor copy end

evaluate ($count = 0)
while ($count < $end_count ) loop main

evaluate ($count=$count+1)

coor swap end
coor copy end

{* Initial minimization.*}
restraints dihedral scale=5. end
noe asymptote * 0.1 end
parameter nbonds repel=1. end end
igroup interaction
(all) (all) weights * 1 vdw 0.002 end end
minimize powell nstep=50 drop=10. nprint=25 end


{* High-temperature dynamics.*}
igroup interaction (all) (all)
weights * 1 angl 0.4 impr 0.1 vdw 0.002 end end

evaluate ($nstep1=int($high_steps * 2. / 3. ) )
evaluate ($nstep2=int($high_steps * 1. / 3. ) )

dynamics cartesian
nstep=$nstep1
timestep=0.003
temperature=$init_t
tcoup=true
nprint=50
ntrfrq=9999
cmremove=true
cmperiodic=0
end


{* Tilt the asymptote and increase weights on geometry.*}
noe asymptote * 1.0 end

igroup interaction
(all) (all) weights * 1 vdw 0.002 end end

{* Bring scaling factor for S-S bonds back *}
parameter
bonds ( name SG ) ( name SG ) 1000. TOKEN
angle ( name CB ) ( name SG ) ( name SG ) 500. TOKEN
end

dynamics cartesian
nstep=$nstep2
timestep=0.003
tcoup=true
temperature=$init_t
nprint=50
ntrfrq=9999
cmremove=true
cmperiodic=0
end

{* Cool the system.*}

restraints dihedral scale=200. end

evaluate ($final_t = 100) { K }
evaluate ($tempstep = 50) { K }

evaluate ($ncycle = ($init_t-$final_t)/$tempstep)
evaluate ($nstep = int($cool_steps/$ncycle))

evaluate ($ini_rad = 0.9) evaluate ($fin_rad = 0.75)
evaluate ($ini_con= 0.003) evaluate ($fin_con= 4.0)

evaluate ($bath = $init_t)
evaluate ($k_vdw = $ini_con)
evaluate ($k_vdwfact = ($fin_con/$ini_con)^(1/$ncycle))
evaluate ($radius= $ini_rad)
evaluate ($radfact = ($fin_rad/$ini_rad)^(1/$ncycle))

evaluate ($i_cool = 0)
while ($i_cool < $ncycle) loop cool
evaluate ($i_cool=$i_cool+1)

evaluate ($bath = $bath - $tempstep)
evaluate ($k_vdw=min($fin_con,$k_vdw*$k_vdwfact))
evaluate ($radius=max($fin_rad,$radius*$radfact))

parameter nbonds repel=$radius end end
igroup interaction (all) (all)
weights * 1. vdw $k_vdw end end

dynamics cartesian
nstep=$nstep
time=0.005
tcoup=true
temperature=$bath
nprint=$nstep
ntrfrq=9999
cmremove=true
cmperiodic=0
end

{====>} {*Abort condition.*}
evaluate ($critical=$temp/$bath)
if ($critical > 10. ) then
display ****&&&& rerun job with smaller timestep (i.e., 0.003)
stop
end if

end loop cool

{* ================================================= Final minimization.*}

igroup interaction (all) (all) weights * 1. vdw 1. end end
parameter
nbonds
repel=0.80
rexp=2 irexp=2 rcon=1.
nbxmod=3
wmin=0.01
cutnb=6.0 ctonnb=2.99 ctofnb=3.
tolerance=1.5
end
end

minimize powell nstep=1000 drop=10.0 nprint=25 end

{* =================================== Write out the final structure(s).*}
print threshold=0.5 noe
evaluate ($rms_noe=$result)
evaluate ($violations_noe=$violations)
print threshold=5. cdih
evaluate ($rms_cdih=$result)
evaluate ($violations_cdih=$violations)
print thres=0.05 bonds
evaluate ($rms_bonds=$result)
print thres=5. angles
evaluate ($rms_angles=$result)
print thres=5. impropers
evaluate ($rms_impropers=$result)
remarks ================================================== =============
remarks overall,bonds,angles,improper,vdw,noe,cdih
remarks energies: $ener, $bond, $angl, $impr, $vdw, $noe, $cdih
remarks ================================================== =============
remarks bonds,angles,impropers,noe,cdih
remarks rms-d: $rms_bonds,$rms_angles,$rms_impropers,$rms_noe,$rm s_cdih
remarks ================================================== =============
remarks noe, cdih
remarks violations.: $violations_noe, $violations_cdih
remarks ================================================== =============

{====>} {*Name(s) of the family of final structures.*}
evaluate ($filename="sa_cns_"+encode($count)+".pdb")

write coordinates output =$filename end

end loop main

stop
Reply With Quote


Did you find this post helpful? Yes | No

Reply
Similar Threads
Thread Thread Starter Forum Replies Last Post
[NMRpipe Yahoo group] Structures Refine by DYNAMO
Structures Refine by DYNAMO Hello! Anyone*could help me with*some information*about*this*step in*DYNAMO? Any information will be helpful! 1. Additional structure calculations may be More...
NMRpipe Yahoo group news News from other NMR forums 0 02-06-2012 03:54 PM
[CNS Yahoo group] Re: Refine.inp
Re: Refine.inp ... Most likely your CNS_SOLVE environmental variable is not set properly. Check your cns_solve_env. To verify, you can also go back to the command line and More...
nmrlearner News from other NMR forums 0 02-06-2012 03:54 PM
[CNS Yahoo group] Refine.inp
Refine.inp Hi i am trying to run the CNS with the refine.inp but it is giving below error. Please help if anybody knows how to break it. EVALUATE: symbol $LOG_LEVEL set More...
nmrlearner News from other NMR forums 0 02-06-2012 03:54 PM
[NMRpipe Yahoo group] Structures Refine by DYNAMO
Structures Refine by DYNAMO Hello! Anyone*could help me with*some information*about*this*step in*DYNAMO? Any information will be helpful! 1. Additional structure calculations may be More...
NMRpipe Yahoo group news News from other NMR forums 0 01-31-2012 08:29 AM
[CNS Yahoo group] Re: Refine.inp
Re: Refine.inp ... Most likely your CNS_SOLVE environmental variable is not set properly. Check your cns_solve_env. To verify, you can also go back to the command line and More...
nmrlearner News from other NMR forums 0 12-28-2011 05:16 AM
[CNS Yahoo group] Refine.inp
Refine.inp Hi i am trying to run the CNS with the refine.inp but it is giving below error. Please help if anybody knows how to break it. EVALUATE: symbol $LOG_LEVEL set More...
nmrlearner News from other NMR forums 0 12-25-2011 10:35 PM
Hydrogen-bonding potential to refine NMR structure
An Empirical Backbone-Backbone Hydrogen-Bonding Potential in Proteins and Its Applications to NMR Structure Refinement and Validation Alexander Grishaev and Ad Bax J. Am. Chem. Soc.; 2004; 126(23) pp 7281 - 7292 http://pubs.acs.org./isubscribe/journals/jacsat/126/i23/figures/ja0319994n00001.gif Abstract: A new multidimensional potential is described that encodes for the relative spatial arrangement of the peptidyl backbone units as observed within a large database of high-resolution X-ray structures. The detailed description afforded by such an analysis provides an opportunity to study...
nmrlearner Journal club 0 06-29-2005 04:16 AM
CNS water refinement scripts from RECOORD website
CNS scripts for water refinement of NMR-derived protein structures are available on the RECOORD website. Reference: RECOORD: A recalculated coordinate database of 500+ proteins from the PDB using restraints from the BioMagResBank Aart J. Nederveen 1, Jurgen F. Doreleijers 2, Wim Vranken 3, Zachary Miller 4, Chris A.E.M. Spronk 5, Sander B. Nabuurs 5, Peter Güntert 6, Miron Livny 4, John L. Markley 2, Michael Nilges 7, Eldon L. Ulrich 2, Robert Kaptein 1, Alexandre M.J.J. Bonvin 1 * 1Bijvoet Center for Biomolecular Research, Utrecht University, Utrecht, The Netherlands 2Center...
nmrlearner NMR structure calculation 0 05-21-2005 07:45 PM


Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes Rate This Thread
Rate This Thread:

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is On
Trackbacks are Off
Pingbacks are Off
Refbacks are Off



BioNMR advertisements to pay for website hosting and domain registration. Nobody does it for us.



Powered by vBulletin® Version 3.7.3
Copyright ©2000 - 2024, Jelsoft Enterprises Ltd.
Copyright, BioNMR.com, 2003-2013
Search Engine Friendly URLs by vBSEO 3.6.0

All times are GMT. The time now is 10:09 AM.


Map