if (_applet) {
  script https://proteopedia.org/wiki/images/6/65/Molecular_playground_helpers.spt
  load https://proteopedia.org/wiki/images/e/ea/Cellulose.pngj
} else {
  script Molecular_playground_helpers.spt
  load UvrB_1d9z.pngj
}

spin y 40
spin on
delay2play(2, "1D fiber")

moveto 1.0 { -560 -530 -637 129.92} 465.24 0.0 0.0 {12.774068181818182 24.421159090909086 16.725386363636364} 29.39059369311233 {0 0 0} 0 0 0 3.0 0.0 0.0;
moveto 1.0 { -526 -593 -610 123.68} 140.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
display chain=M						# 1D
select not (4-5 and chain=M);wireframe -0.2;
select all; hbond on; hbond 0.15; color HBONDS white
delay2play(2, "2D fiber")

spin off;
moveto 1.0 { -40 704 709 177.29} 140.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
moveto 1.0 { -526 -593 -610 123.68} 140.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
display chain=B or chain=G or chain=M or chain=R 	# 2D
delay2play(2, "3D fiber")

color chain; select (4-5 and chain=M);set bondmode OR; wireframe -0.3
moveto 1.0 { -40 704 709 177.29} 140.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
display all						# 3D
moveto 2.0 { 8 14 -1000 95.17} 140.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
delay2play(2, "space filling")


select all;spacefill only; delay 0.5;
moveto 2.0 { -674 624 -396 141.7} 140.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
moveto 2.0 { 685 -617 -388 139.45} 140.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
moveto 2.0 { 8 14 -1000 95.17} 140.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
delay2play(2, "nonpolar contacts")

display (chain=M and 4-5) or (3-5 and (:F, :H, :L, :N))
moveto 1.0 { 8 14 -1000 95.17} 300.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0; # zoomto 1.0 300
delay 0.5
select (4-5 and chain=M); wireframe -0.2;
contact ID "contact1" (4-5 and chain=M) (not(:G, :M, :R)) surface;
contact ID "contact1" fill noMesh noDots notFrontOnly frontlit;
color $"contact1""roygb" range -0.5 1;
delay 0.5;
select not (4-5 and chain=M);color cpk; wireframe -0.2
moveto 1.0 { -586 548 -596 123.98} 300.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
moveto 1.0 { -48 998 -41 91.43} 300.0 0.0 0.0 {12.7961052631579 24.454105263157892 16.85585380116959} 29.27020811620052 {0 0 0} 0 0 0 3.0 0.0 0.0;
delay2play(2, "repeating unit")

display chain=M and 4-5 				# none
moveto 1.0 { -49 639 767 175.38} 465.24 0.0 0.0 {12.774068181818182 24.421159090909086 16.725386363636364} 29.39059369311233 {0 0 0} 0 0 0 3.0 0.0 0.0;
spin on; delay 0.5
select 4-5 and chain=M;set bondmode OR; wireframe -0.2; spacefill 25%; color cpk
contact off


if (not _applet) {
	banner("For next molecule, shrink this one")
}
message MP_done