PROLSQ: Constrained Refinement of Atomic Parameters

archiv specify items to be placed in reflection records

uov overall isotropic thermal parameter

dampsh shift damping factors to control refinement

print print control line

reflim limits on reflections to be included in refinement

rtest set up survey of R values to assess refinement

scale scale factors

scanis anisotropic scale factor

solmod solvent model parameters

stats number and limits of d shells for R value survey

wgtb set thermal restraint weighting

wgtchi set chiral centre restraint weighting

wgtcon set contact distance restraint weighting

wgtdis set distance restraint weighting

wgtpla set plane restraint weighting

wgtref set structure factor restraint weighting

wgtshf set excessive shift restraint weighting

wgtsym set noncrystallographic symmetry restraint weighting

wgttor set torsion angle restraint weighting

noref control atoms included in refinement

PROLSQ Option Code

congugent gradient cycles ncy n 50

Lagrange multiplier scale lag n 10

refinement to fit diffraction data1 xrfapply

nxr skip

refinement to fit idealized geometry 1 geoapply

nge skip

set individ, tf's to overall value tfs set tf's to overall

refine populations if input <1.0 pop refine populations

anomolous dispersion noa do not apply

ano apply dispersion

thermal parameters ova overall isotropic

iso individual isotropic

data set number dsn n 1

constrain origin orx along x no

ory along y no

orz along z no

punch atoms in PDB format pch yes no

1 The default is to use both types of information in the refinement. If both types are set to skip the result will be a simple copying of the bdf from one file to another, i.e. a null refinement.

archiv 1- ID numbers of data items to be added to or deleted from the bdf

Positive ID numbers signal addition, negative deletion.

uov 1 B overall or U overall depending on magnitude U from input bdf

If value is 0 < U < 0.2 the quantity is stored as U. If the quantity is 0.2 < B < 40.0 it is multiplied by 78.9568 and stored as U. This line is ignored unless in overall t.f. mode.

dampsh 1-3 damping factor for x, y , z parameter shifts 0.5

4 damping factor for B parameter shifts 0.5

5 damping factor for population parameter shifts 0.5

6-11 damping factor for B11 to B23 parameter shifts 0.5

These damping factors are to be applied to parameter shifts determined by the least-squares process. These factors will be used to establish the new parameters written to the output bdf and, optionally on unitpch:. Factors between 0 & 4 allowed.

N.B. Unequal damping factors will distort the geometry.

print Option Code

list first m reflections lfir m 10

list reflection data flis m every mth refln 0

list first n input atoms ilis n 99999

list updated atoms alis m 0/1/2 no/all/>shift limit

atom shifts on x, y, z parc d only shifts > d 0.1

temperature factor shifts bcut d 0.25

population param. shifts popc d 0.25

torsion angle printing ntor statistics only

stor all standards and statistics

ctor d single angle contacts 0.1

mtor d multiple 0.1

chiral centres printing nchi statistics only

wchi worst case chiral centres

achi all chiral centres

hydrogen bond printing nhbo statistics only

hbon all possible H bonds

restrained distances1 disc d dists of type code "J" 5

restrained planes plan d all planes with RMS 2

deviation > d*sigp

1 Print only if |DIDEAL-DMODEL| >= d*SIGD(J) If d=0.0, print all.

reflim Option Code

discard reflections fmin x < x 0

smin x sinθ/λ < x 0

smax x sinθ/λ > x 1.5

dmax x d maximum 1000

dmin x d minimum 0.35

sigm x < σ 0

Note: Use either d or sinθ/λ specs but not both.

rtest 1 sampling frequency for included reflections1 10000

2 sampling frequency for excluded reflections2 0

3 sampling offset of included reflection subset 0

4 sampling offset of excluded reflection subset 0

5 number of coordinate, x, y, z, shift steps to be tested 1

6 number of B steps to be tested 1

7-n factors to multiply shifts for testing3 0.5

Set up R testing for reflection subsets. The subsets are defined by the reflection count plus a constant modulo the frequency being zero. After the least-squares refinement structure factors are recalculated for the subsets in order to evaluate the efficacy of the shift damping and progress of the refinement.

1Time consuming if this is set too low and fields 5 & 6 too large.

2 Excluded from the refinement but used in subsequent structure factor R value estimate. This gives a stringent test of convergence, but caution is required or the number of reflections eliminated may be too great.

3 These are the same kind of factors as those in the dampsh line, e.g. x(test) = x(start) +factor*x(shift) The sum of fields 5 & 6 determines the value of n.

scale 1 scale factor input bdf

2 scale group number 1

scanis 1 component of overall anisotropic scaling of 0

2-6 components

solmod 1 linear scale K for partial 1

2 exponential scale U for partial 0

3 linear scale UK for analytic solvent attenuation function 1

4 exp. scale UB for analytic solvent attenuation function 0

5 electron density plateau level for unoccupied model space 0

stats 1-7 resolution limits for reciprocal space shells used for R-factors

Maximum of seven shells; default one. Values supplied must be descending order of d-spacings in Ĺ.

wgtb 1 weight for temperature factor restraints 1

2 σ of type 1 main-chain bond distances 1

3 σ of type 2 main-chain angle distance 1.5

4 σ of type 3 side-chain bond distance 1

5 σ of type 4 side-chain angle distance 1

6 σ of type 5 X-H bond 0.5

7 σ of type 6 X-H angle 1

8 σ of type 7 hydrogen-bond, metal, etc 0

wgtchi 1 weight for chiral centre restraints WC 1

2 σ(χ) 0.15

wgtcon 1 weight for van der Waals contacts WV 0.5

2 σ(vdW) 0.5

3 type 1 distance factor in Ĺ1 -0.3

4 type 2 distance factor in Ĺ1 0

5 type 3 distance factor in Ĺ1 -0.2

6 type 4 distance factor in Ĺ1 0

1 These fields may be used to supply parameters to modify the minimum "theoretical" van der Waals contact distance of four possible types which are defined and calculated in PROTIN. There are four parameters corresponding to:

1) Atoms in which the relative position is determined by only one torsion angle.

2) Two atoms in which >two torsion angles are involved, multiple torsion contact.

3) A possible (X....Y) hydrogen bond pair involving contacts between nitrogen-oxygen pairs, but not N main - N main or O main - O main.

4) A possible (X-H....Y) hydrogen bond.

Note: The distance factors must be negative amounts. They are used to adjust the ideal contact distances supplied from PROTIN in order to establish when a contact should be restrained.

wgtdis 1 weight for distance restraints WD 1

2 σ of type1 1 bonded distances (e.g. Cα - C) 0.03

3 σ of type1 2 angle distances (e.g. N - C) 0.04

4 σ of type1 3 planar 1-4 distances (e.g. O i-1 - Cα1) 0.05

5 σ of type1 4 bond distances involving H (e.g. Cα- H) 0.1

6 σ of type1 5 angle distance involving H (e.g. Hα- C) 0.08

7 σ of type1 6 hydrogen bond, metal 0

8 σ of type1 7 coordination, etc 0

wgtpla 1 weight for plane restraints, WP 1

2 σ of a plane 0.02

The weight applied to the derivative is: (WP/(σ(p))2.

wgtref 1 weighting scheme as 1 through 5 1

2 A term for weights 4

3 B term for weights 0

wgtshf 1 shifts magnitude restraint positional PDEL 0.2

2 shifts magnitude restraint thermal BDEL 3

3 shifts magnitude restraint population QDEL 0.1

4 σ for tight positional restraints 0.05

5 σ for medium positional restraints 0.5

6 σ for loose positional restraints 2.0

7 σ for tight thermal restraints 0.5

8 σ for medium thermal restraints 1

9 σ for loose thermal restraints 0

wgtsym 1 weight scale non-crystallographic symmetry restraints WS 1

wgttor 1 weight scale torsion angle restraints 1

2 σ(T) secondary structure angle φ or ψ 15

3 σ(T) planar angle such as ω 3

4 σ(T) staggered angle such as χ1 15

5 σ(T) orthonormal angle such as χ2 of aromatics 20

noref 1- atom labels for atoms which are not refined

More than one noref line may be entered.