These centroid score terms appear in the score0-score5 score functions used in the Abinitio protocol.
The env
score term is a context-dependent one-body energy term that describes the solvation of a particular residue (based on the hydrophobic effect). It is based on the probability of a residue having the specified type given its number of neighboring residues. (Rohl et al. 2004)
The pair
score term is a two-body energy term for residue pair interactions (electrostatics and disulfide bonds). For each pair of residues, it is based on the probability that both of these two residues will have their specified types given their sequence separation and the physical distance between them, normalized by the product of the probabilities that each residue will have its specified type given the same information. (Rohl et al. 2004)
The cbeta term is another solvation term intended to correct for the excluded volume effect introduced by the simulation and favor compact structures. It is based on the ratio of probabilities of a residue having a given number of neighbors in a compact structure vs. random coil and summed over all residues. (Rohl et al. 2004)
The vdw term represents only steric repulsion and not attractive van der Waals forces (those are modeled in terms rewarding compact structures, such as the rg term; local interactions are implicitly included from fragments). It is calculated over pairs of atoms only in cases where: 1. the interatomic distance is less than the sum of the atoms' van der Waals radii, and 2. the interatomic distance does not depend on the torsion angles of a single residue. (Rohl et al. 2004)
The rg term favors compact structures and is calculated as the root mean square distance between residue centroids. (Rohl et al. 2004)
The hs_pair term describes packing between strands and helices. It is based on the probability that two pairs of residues (1 pair in the sheet and 1 pair in the helix) will have their current dihedral angles given the separation (in sequence and physical distance) between the helix and the strand. (Rohl et al. 2004)
The ss_pair term describes hydrogen bonding between beta strands.
It is often used with the STRAND_STRAND_WEIGHTS statement (e.g. STRAND_STRAND_WEIGHTS 1 6
):
The sheet term favors the arrangement of individual beta strands into sheets. It is derived from the probability that a structure with a given number of beta strands will have the current number of beta sheets and lone beta strands. (Rohl et al. 2004)
Scores strand pairs based on the distance between them and the register of the two strands. (Shmygelska and Levitt 2008)
Contact order (Shmygelska and Levitt 2008)
As in full-atom scoring, the rama term describes the probability of a residue having a given set of torsion angles given its identity. (Shmygelska and Levitt 2008)
Short-range backbone-backbone hydrogen bonding energy term (Shmygelska and Levitt 2008)
Long-range backbone-backbone hydrogen bonding energy term (Shmygelska and Levitt 2008)
Centroid score terms that have been optimized with smoothed data and a larger dataset than the original.
These terms are commonly used in centroid docking protocols.
These terms are used in conjunction with centroid rotamer modeling (see the CenRotModel page). NOTE: To use these terms, you must load the centroid rotamer library using the -corrections:score:cenrot
flag.
See env
See pair
???
See cbeta
Dunbrack term for centroid rotamer modeling. Scores internal energy of centroid rotamers.
These score terms are specialized for use with membrane proteins and appear in membrane-specific score functions. Descriptions from Alford et al. 2015.
Knowledge-based pairwise interaction potential between two residues some distance apart at a given depth in the membrane. Originally published in Yarov-Yarovoy et al. 2006 (see Epair description).
Knowledge-based potential describing propensity for a single residue to be at a given depth in the membrane and burial by residues. Originally published in Yarov-Yarovoy et al. 2006 (see Eenv description).
Knowledge-based residue density potential based on number of neighbor residues and conditional upon number of transmembrane helices. Originally published in Yarov-Yarovoy et al. 2006 (see Edensity description). See also Rohl et al. 2004 for more information on residue density calculation.
Scores agreement between predicted lipophilicity (from LIPS server) and the model. See Adamian et al. 2005 for additional information on lipophilicity predictions.
Penalty for non-helical secondary structure in the membrane. Originally published in Yarov-Yarovoy et al. 2006 (see "Search for Embeddings" section).
Penalty for residues outside of the hydrophobic layer of the membrane. Originally published in Yarov-Yarovoy et al. 2006 (see "Search for Embeddings" section).
Penalty for transmembrane helices that project outside of the membrane. Originally published in Yarov-Yarovoy et al. 2006 (see "Search for Embeddings" section).
These score terms are used in coarse-grained modeling of RNA.
See the rg term above.
See the vdw term above.
Centroid score function that uses only the env, pair, cbeta, and vdw score terms with equal weights.
Score function used in the first stage of the ClassicAbInitio protocol. It is intended for use with fragment insertion; the only active term is a low vdw (repulsive) weight to penalize clashes, and the fragments themselves are intended to provide the majority of the information.
Score function used in the second stage of the ClassicAbInitio protocol. Uses the env, pair, vdw, hs_pair, and sheet terms with equal weights and the ss_pair term with a lower weight.
Score function used in the third stage of the ClassicAbInitio protocol. Relative to score1, it turns on the cbeta and cenpack terms at low weights and upweights the ss_pair term. It also decreases the sequence distance cutoff for STRAND_STRAND_WEIGHTS from 11 to 6.
Score function used in the fourth stage of the ClassicAbInitio protocol. Upweights the cbeta and cenpack terms and turns on the rg and rsigma terms. Sets the value of ss_lowstrand (from STRAND_STRAND_WEIGHTS) to 0 (no distance score for strand pairs).
While there is no base score4 score function, it is available as a patch and in smooth and cenrot versions. It introduces hydrogen bonding terms (hbond_lr_bb and hbond_sr_bb) and a rama term (ramachandran probability-based scores for residues based on phi/psi angles).
Score function used in the third stage of the ClassicAbInitio protocol. Nearly identical to score2, but uses different STRAND_STRAND_WEIGHTS cutoff (see above)
Formerly score6. Introduces a rama term, turns off scoring of rg, rsigma, sheet, and hs_pair, and reweights several terms (increases cenpack and cbeta; decreases pair, ss_pair, and vdw).
These score functions are intended for use with a particular protocol or type of sampling.
Commonly used in centroid docking; basically cen_std, but only evaluated between two chains and not within a single chain.
Very coarse score function used in SlidingWindowLoopClosure. Includes only env, pair, and vdw centroid terms plus the linear_chainbreak and overlap_chainbreak terms.
Standard weight set for the low-resolution membrane score function.
Used for centroid docking of membrane proteins (MP_Dock protocol)
Used for centroid symmetric docking of membrane proteins (MP_SymDock protocol)
Formerly used by the MPDockingMover (replaced by mpframework_docking_cen_2015).
These may be seen as suffixes appended to the name of a score function. Some other suffixes (relax, design, etc.) may also appear; these indicate the protocol for which the score function is intended.
These score functions use the optimized "smooth" centroid weights described above.
These score functions are intended for use with centroid rotamer modeling (see the cen_rot terms section above).
These score functions/patches have increased repulsive (vdw) weights.
These score functions/patches are intended for use with cartesian minimization and have the cart_bonded weight set.
Doubles the repulsive (vdw) weight of the score function.
Sets all of the interchain weights to the values that they have in the interchain_cen score function. Does not affect other score terms.
Turns on scoring for short-range and long-range backbone hydrogen bonding, a ramachandran term (phi/psi given residue type), and a chainbreak term; sets rg weight to 2. L is for loop.
Sets all of the standard centroid score term weights to their values in the listed score function. (Interestingly, the vdw weight in the score0 patch is 1.0, while in the main score0 score function, it is 0.1. Is this a bug??)
The best centroid score function for your protocol will (of course) be protocol-dependent.