#CreatePredAlphaCost pred_alpha2k alpha11 T0331.t2k.alpha.rdb 2 alpha alpha_prev sum
#CreatePredAlphaCost pred_alpha04 alpha11 T0331.t04.alpha.rdb 2 alpha alpha_prev sum
#CreatePredAlphaCost pred_alpha06 alpha11 T0331.t06.alpha.rdb 2 alpha alpha_prev sum

SetCost wet6.5 15 near_backbone 5 way_back 5 dry5 15 dry6.5 20 dry8 15 dry12 5 \
	phobic_fit 2 \
	sidechain 5 \
	n_ca_c	5 bad_peptide 10 \
	bystroff 5 \
	soft_clashes 20  backbone_clashes 2 \
	break 50  \
	constraints 10 \
	hbond_geom 5 \
	hbond_geom_backbone 10 \
	hbond_geom_beta 50 \
	hbond_geom_beta_pair 100 \
	missing_atoms 1 \
	maybe_metal 0.5 \
	maybe_ssbond 1

// remove maybe_ssbond weight if protein known to be in reducing environment.
// remove maybe_metal weight if protein known to have disulfides or
//	known not to bind metal ions.

// Add the following for extra-cellular proteins:
#	known_ssbond 0.5	# and put ssbond command in constraints file

// For comparative modeling, increase weight of break and hbond_geom...

// put any constraints here, such as helix constraints, strand constraints,
//	sheet constraints, ssbonds, hbonds, or just arbitrary pairs of atoms.

KnownBreak M150

#include T0331.dssp-ehl2.constraints
#include T0331.undertaker-align.sheets
#include rr.constraints

SheetConstraint (T0331)K14 (T0331)D22	(T0331)I34 (T0331)N26	hbond (T0331)G16	1
SheetConstraint (T0331)G16 (T0331)T20	(T0331)A66 (T0331)V62	hbond (T0331)V17	1
SheetConstraint (T0331)A36 (T0331)N38	(T0331)F43 (T0331)G41	hbond (T0331)N38	1
SheetConstraint (T0331)G41 (T0331)S47	(T0331)L149 (T0331)H143	hbond (T0331)I42	1
SheetConstraint (T0331)Q60 (T0331)A66	(T0331)A83 (T0331)V77	hbond (T0331)V62	1
SheetConstraint (T0331)T82 (T0331)V86	(T0331)L149 (T0331)E145	hbond (T0331)R84	1

HelixConstraint (T0331)M1	(T0331)V15	1
HelixConstraint (T0331)T50	(T0331)D59	1
HelixConstraint (T0331)E87	(T0331)N98	1
HelixConstraint (T0331)P99	(T0331)K106	1
HelixConstraint (T0331)I118	(T0331)Q130	1
HelixConstraint (T0331)L128	(T0331)K133	1