minAmber.pl - Revision history

Edit: /* Options */

2009-07-31T03:44:34Z

Options

Edit: /* Description */

2009-07-31T03:41:53Z

Description

Edit: /* Description */

2009-07-31T03:41:32Z

Description

Edit: /* Description */

2009-07-31T03:40:45Z

Description

Edit: /* Description */

2009-07-31T03:40:11Z

Description

Feig at 15:55, 8 August 2008

2008-08-08T15:55:40Z

66.232.105.72: /* Usage */

2008-05-27T14:08:24Z

Usage

58.65.234.1: /* Options */

2007-03-08T10:19:04Z

Options

192.26.250.230: /* Description */

2006-08-12T17:53:21Z

Description

Feig at 14:17, 12 August 2006

2006-08-12T14:17:38Z

New page

__NOTOC__
== Usage ==

<mmtsbToolUsage cmd="minAmber.pl"></mmtsbToolUsage>

== Description ==

This script is used to minimize a protein structure from a PDB file with
Amber. Without any options a 50-step steepest descent minimization
is followed by the a conjugate gradient minimization algorithm over 500
steps or until the energy decrease between steps becomes less than
1.0E-5 kcal/mol. With the default protocol the minimization is carried in
vacuum (constant dielectric of 1.0) with a 14.0 A cutoff for non-bonded interactions. 
Command-line options are available to change the default protocol, introduce
restraints and specify output files. 
Minimization parameters can be set with -par followed by key=value pairs.
A complete list of parameters is available
<A HREF="http://mmtsb.scripps.edu/cgi-bin/amberpar" TARGET=apar>here</A>.
The most relevant options are explained in the following:
The number of steps in the final conjugate gradient minimization run is set
with minsteps, a different energy tolerance is selected with
minetol. The number of steps during the initial steepest descent
minimization is changed with sdsteps.
The type of dielectric environment can be set as constant with dielec=cdie or
distance dependent with dielec=rdie. The value of epsilon can be changed
using epsilon. The non-bonded cutoff radius is set with cutoff or
switched off completely with nocut.
The hsd and hse parameters are available to specify explicitly
the protonation site for histidine residues. By default HSD is assumed
(protonated at ND1). 
By default, the PARAM94 force field is used. Other Amber force fields may be
selected with -param. More information on differences between these force fields
is available from the <A HREF="http://www.amber.ucsf.edu">Amber web site</A>. 
More realistic than constant or distance-dependent dielectrics but less expensive than
explicit solvent is an implicit solvent approximation based on a generalized Born (GB)
formalism. This function is available with the gb switch of the -par
option. It needs to be set to one of the methods implemented in Amber. Possible options
are: tc: Tsui & Case, tc6 Tsui & Case (original radii),
obc: Onufriev, Bashford, Case, jsb: Jayaram, Sprous, Beveridge,
mgb: modified Jayaram, Sprous, Beveridge. Again, more information is available
on the Amber web site. 
Different types of restraints may be specified to keep parts of the structure
near their initial positions or near coordinates from a reference file. This is
useful for modeling of fragments or loops in otherwise known structure templates.
Atom-based restraints are specified with -cons or -l.
-cons requires a list of residue ranges where the restraints should be
applied and corresponding force constants. By default all heavy atoms will be
restrained to their initial positions. The keywords ca,
cb, or cab can be specified to select only C-alpha, C-beta,
or C-alpha/C-beta restraints, respectively. If a file name is given, the reference
atom positions from this file are used instead of the initial coordinates from the
input file. 
In loop modeling it is more convenient to specify the residues in the
loop rather than the residues outside the loop for which restraints would be requested.
The option -l is available to specify a range of residues that are kept
flexible, while restraints are generated implicitly for all other residues
outside the given residue range(s). The syntax is otherwise similar to
-cons. However, only a single force constant can be
supplied as an extra argument. The restraints are generated with continuously
increasing force constants from the beginning and ends of the loop regions up to
the specified force constant in a similar fashion as the force constants generated by
<docmark>genrestr.pl</docmark> for lattice simulations. 
The final structure after minimization is written to standard output if
input is read from standard input or to a file name derived from the input
file name by adding min. An alternate output
file name can be specified with -out. 
The options -log and -elog are available to save the complete
output from Amber or only the energy values at each printed minimization
step, respectively 
The location of the Amber binary and data directory are expected
from the environment variables AMBERHOME, SANDEREXEC and LEAPEXEC. At least
AMBERHOME should be set accordingly. Otherwise an
attempt is made to locate the Amber binary from the search path and deduce the data
directory from the binary location. The current version of this utility expects
Amber 7 for full functionality. 
A similar script <docmark>minCHARMM.pl</docmark> is available for minimizations
using CHARMM.

== Options ==

; -help : usage information

== Examples ==

<mmtsbToolExample cmd="minAmber.pl" set="/apps/mmtsb/bench/minAmber.pl-test/test1"></mmtsbToolExample>

<mmtsbToolExample cmd="minAmber.pl" set="/apps/mmtsb/bench/minAmber.pl-test/test2"></mmtsbToolExample>

<mmtsbToolExample cmd="minAmber.pl" set="/apps/mmtsb/bench/minAmber.pl-test/test3"></mmtsbToolExample>

<mmtsbToolExample cmd="minAmber.pl" set="/apps/mmtsb/bench/minAmber.pl-test/test4"></mmtsbToolExample>

@@ Line 34: / Line 34: @@
 By default, the PARAM94 force field is used. Other Amber force fields may be
 selected with <B>-param</B>. More information on differences between these force fields
-is available from the [http://www.amber.ucsf.edu Amber web site].<BR>
+is available from the [http://ambermd.org Amber web site].<BR>
 More realistic than constant or distance-dependent dielectrics but less expensive than

@@ Line 12: / Line 12: @@
 .0E-5 kcal/mol. With the default protocol the minimization is carried in
 vacuum (constant dielectric of 1.0) with a 14.0 A cutoff for non-bonded interactions.<BR>
 Command-line options are available to change the default protocol, introduce
-restraints and specify output files.<BR><BR>
+restraints and specify output files.<BR>
 Minimization parameters can be set with <B>-par</B> followed by key=value pairs.
 A complete list of parameters is available
-<A HREF="http://mmtsb.scripps.edu/cgi-bin/amberpar" TARGET=apar>here</A>.
+[http://feig.bch.msu.edu/mmtsb/Amber_Parameters here].
 The most relevant options are explained in the following:
 The number of steps in the final conjugate gradient minimization run is set
@@ Line 28: / Line 30: @@
 The <B>hsd</B> and <B>hse</B> parameters are available to specify explicitly
 the protonation site for histidine residues. By default HSD is assumed
-(protonated at ND1).<BR><BR>
+(protonated at ND1).<BR>
 By default, the PARAM94 force field is used. Other Amber force fields may be
 selected with <B>-param</B>. More information on differences between these force fields
-is available from the <A HREF="http://www.amber.ucsf.edu">Amber web site</A>.<BR><BR>
+is available from the [http://www.amber.ucsf.edu Amber web site].<BR>
 More realistic than constant or distance-dependent dielectrics but less expensive than
 explicit solvent is an implicit solvent approximation based on a generalized Born (GB)
@@ Line 51: / Line 55: @@
 atom positions from this file are used instead of the initial coordinates from the
 input file.<BR>
 In loop modeling it is more convenient to specify the residues in the
 loop rather than the residues outside the loop for which restraints would be requested.
@@ Line 60: / Line 65: @@
 increasing force constants from the beginning and ends of the loop regions up to
 the specified force constant in a similar fashion as the force constants generated by
-<docmark>genrestr.pl</docmark> for lattice simulations.<BR><BR>
+<docmark>genrestr.pl</docmark> for lattice simulations.<BR>
 The final structure after minimization is written to standard output if
 input is read from standard input or to a file name derived from the input
 file name by adding <font color=#508060>min</font>. An alternate output
 file name can be specified with <B>-out</B>.<BR>
 The options <B>-log</B> and <B>-elog</B> are available to save the complete
 output from Amber or only the energy values at each printed minimization
-step, respectively<BR><BR>
+step, respectively<BR>
 The location of the Amber binary and data directory are expected
 from the environment variables AMBERHOME, SANDEREXEC and LEAPEXEC. At least
@@ Line 73: / Line 81: @@
 attempt is made to locate the Amber binary from the search path and deduce the data
 directory from the binary location. The current version of this utility expects
-Amber 7 for full functionality.<BR><BR>
+Amber 7 for full functionality.<BR>
-A similar script <docmark>minCHARMM.pl</docmark> is available for minimizations
 using CHARMM.

←Older revision		Revision as of 10:19, 8 March 2007
Line 80:		Line 80:

	; -help : usage information		; -help : usage information
−

	== Examples ==		== Examples ==

@@ Line 89: / Line 89: @@
 ; -help : usage information
+; -par key=value[,...] : Amber simulation paramters
 == Examples ==

@@ Line 1: / Line 1: @@
 __NOTOC__
-== Usage ==
+t3Yjg  <a href="http://nkkhneuzofua.com/">nkkhneuzofua</a>, [url=http://apacvxhegerq.com/]apacvxhegerq[/url], [link=http://lewsjbipnerj.com/]lewsjbipnerj[/link], http://rlaaaawiotil.com/
 == Description ==