Fork me on GitHub

Real Space Rotamer Optimization

Discrete optimization using a many-body expansion of a real space target function and elimination expressions.


Usage: ffxc realspace.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] [--wA=1.0] [--window=7] [-a=0] [-L=2] [-s=-1] [-x=-1] [-X=<data> <data>]... files...
Discrete optimization using a many-body expansion and elimination expressions.
    files...                            PDB and Real Space input files.
    --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.
    --wA, --dataWeight=1.0              The weight of the real space data (wA).
    --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.
-X, --data=<data> <data>                Specify input data filename, and its weight (wA) (e.g. -X filename 1.0).
-z, --revert                            Revert unfavorable changes.