Fork me on GitHub

Many-Body Rotamer Optimization

Run Many-Body rotamer optimization on a system.

Usage: ffxc ManyBody [-EhOTz] [--dee] [--out] [--bB=0.0] [--bC=1] [--bL=20.0] [--ch=-1] [--eR=none] [--fi=-1] [--fR=-1,-1] [--increment=3] [--Ln=Richardson] [--lR=none] [--mC=-1] [--nB=3,3,3] [--pr=1] [--radius=2.0] [--tC=3.0] [--thC=3.0] [--window=7] [-a=0] [-L=2] [-s=-1] [-x=-1] files...
Run ManyBody algorithm on a system.
    files...                            PDB input file.
    --bB, --boxBorderSize=0.0           Extent of overlap between optimization boxes in Angstroms.
    --bC, --boxInclusionCriterion=1     Criterion to use for adding a residue to a box: (1) uses C alpha only (N1/9 for nucleic acids), (2) uses any atom, and (3) uses any rotamer
    --bL, --approxBoxLength=20.0        Approximate side lengths of boxes to be constructed (over-rides numXYZBoxes).
    --ch, --chain=-1                    Single character chain ID of the residues to optimize.
    --dee, --deadEnd                    Use dead-end elimination criteria instead of Goldstein criteria.
    --eR, --energyRestart=none          Load energy restart file from a previous run (requires that all parameters are the same).
    --fi, --final=-1                    Final residue to perform the optimization on (-1 exits). For box optimization, final box to optimize.
    --fR, --forceResidues=-1,-1         Force residues in this range to be considered for sliding window radii, regardless of whether they lack rotamers.
    --increment=3                       Sliding window increment.
    --Ln, --libraryNucleic=Richardson   Nucleic acid library to select: [Richardson]
    --lR, --listResidues=none           Choose a list of individual residues to optimize (eg. A11, A24,B40).
    --mC, --monteCarlo=-1               Follow elimination criteria with (n) Monte Carlo steps, or enumerate all remaining conformations, whichever is smaller.
    --nB, --numBoxes=3,3,3              Specify number of boxes along X, Y, and Z (default: 3,3,3)
    --out, --output                     Save eliminated singles and eliminated pairs to a text file.
    --pr, --prune=1                     Prune no clashes (0), only single clashes (1), or all clashes (2)
    --radius=2.0                        The sliding window cutoff radius (Angstroms).
    --tC, --twoBodyCutoff=3.0           Cutoff distance for two body interactions.
    --thC, --threeBodyCutoff=3.0        Cutoff distance for three-body interactions.
    --window=7                          Size of the sliding window with respect to adjacent residues.
-a, --algorithm=0                       Algorithm: default automatic settings (0), independent residues (1), all with rotamer elimination (2), all brute force (3), sliding window (4), or box optimization (5)
-E, --decompose                         Print energy decomposition for the input structure (no optimization!).
-h, --help                              Print this help message.
-L, --library=2                         Ponder and Richards (1) or Richardson (2) rotamer library.
-O, --noOriginal                        Do not include starting coordinates as their own rotamer.
-s, --start=-1                          Starting residue to perform the optimization on (-1 exits). For box optimization, first box to optimize.
-T, --threeBody                         Include 3-Body interactions in the elimination criteria.
-x, --all=-1                            Optimize all residues beginning from the passed value (overrides other options); for box optimization, optimizes all boxes beginning from the passed index.
-z, --revert                            Revert unfavorable changes.