Disulfide.world
Go to the documentation of this file.
1 
3 [CONSTANTS]
4 
5 % Link 1 is a N-C group fixing the loop
6 l1d1:=sqrt(39.8786)
7 l1d2:=sqrt(31.7352)
8 l1d3:=sqrt( 2.3892)
9 l1d4:=sqrt( 2.1050)
10 l1d5:=sqrt(57.3503)
11 l1d6:=sqrt(45.0087)
12 
13 l1x1:=l1d4
14 l1y1:=0
15 l1z1:=0
16 
17 l1x2:=(l1d5^2-l1d6^2+l1x1^2)/(2*l1x1)
18 l1y2:=sqrt(l1d5^2-l1x2^2)
19 l1z2:=0
20 
21 l1x3:=(l1d1^2-l1d2^2+l1x1^2)/(2*l1x1)
22 l1y3:=(l1d2^2-l1d3^2-(l1x3-l1x1)^2+(l1x3-l1x2)^2+l1y2^2)/(2*l1y2)
23 l1z3:=sqrt(l1d1^2-l1x3^2-l1y3^2)
24 
25 % Link 2 is a C
26 l2d4:=sqrt(2.4106)
27 l2d5:=sqrt(2.1050)
28 l2d6:=sqrt(6.1279)
29 
30 l2x1:=l2d4
31 l2y1:=0
32 l2z1:=0
33 
34 l2x2:=(l2d5^2-l2d6^2+l2x1^2)/(2*l2x1)
35 l2y2:=sqrt(l2d5^2-l2x2^2)
36 l2z2:=0
37 
38 % Link s3 is a C
39 l3d4:=sqrt(3.2178)
40 l3d5:=sqrt(2.4106)
41 l3d6:=sqrt(7.9290)
42 
43 l3x1:=l3d4
44 l3y1:=0
45 l3z1:=0
46 
47 l3x2:=(l3d5^2-l3d6^2+l3x1^2)/(2*l3x1)
48 l3y2:=sqrt(l3d5^2-l3x2^2)
49 l3z2:=0
50 
51 % Link 4 is a S-S
52 l4d1:=sqrt(10.0384)
53 l4d2:=sqrt(16.6496)
54 l4d3:=sqrt( 3.2178)
55 l4d4:=sqrt( 3.2133)
56 l4d5:=sqrt( 4.1544)
57 l4d6:=sqrt(10.0348)
58 
59 l4x1:=l4d4
60 l4y1:=0
61 l4z1:=0
62 
63 l4x2:=(l4d5^2-l4d6^2+l4x1^2)/(2*l4x1)
64 l4y2:=sqrt(l4d5^2-l4x2^2)
65 l4z2:=0
66 
67 l4x3:=(l4d1^2-l4d2^2+l4x1^2)/(2*l4x1)
68 l4y3:=(l4d2^2-l4d3^2-(l4x3-l4x1)^2+(l4x3-l4x2)^2+l4y2^2)/(2*l4y2)
69 l4z3:=sqrt(l4d1^2-l4x3^2-l4y3^2)
70 
71 % Link 5 is a C
72 l5d4:=sqrt(2.4105)
73 l5d5:=sqrt(3.2133)
74 l5d6:=sqrt(7.9244)
75 
76 l5x1:=l5d4
77 l5y1:=0
78 l5z1:=0
79 
80 l5x2:=(l5d5^2-l5d6^2+l5x1^2)/(2*l5x1)
81 l5y2:=sqrt(l5d5^2-l5x2^2)
82 l5z2:=0
83 
84 % Link 6 is a C
85 l6d4:=sqrt(2.3892)
86 l6d5:=sqrt(2.4105)
87 l6d6:=sqrt(6.6307)
88 
89 l6x1:=l6d4
90 l6y1:=0
91 l6z1:=0
92 
93 l6x2:=(l6d5^2-l6d6^2+l6x1^2)/(2*l6x1)
94 l6y2:=sqrt(l6d5^2-l6x2^2)
95 l6z2:=0
96 
97 
98 % Radious of the spheres representing atoms (0.25% Van der Waals radii)
99 r:=0.25
100 rC:=1.7*r
101 rN:=1.55*r
102 rS:=1.8*r
103 
104 % Radii of the cylinders representing bonds
105 rCC:=0.15
106 rCN:=0.15
107 rCS:=0.15
108 rSS:=0.15
109 
110 % Colors for the atoms (Carbon/Nitrogen/Sulfur)
111 Cr:=0
112 Cg:=0
113 Cb:=0
114 
115 Nr:=0
116 Ng:=0
117 Nb:=1
118 
119 Sr:=1
120 Sg:=1
121 Sb:=0
122 
123 
124 [LINKS]
125 
126 C1 : sphere rN (0,0,0) color (Nr,Ng,Nb)
127  sphere rC (l1x2,l1y2,l1z2) color (Cr,Cg,Cb)
128  cylinder rCN (0,0,0) (l1x1/2,l1y1/2,l1z1/2) color (Nr,Ng,Nb) decoration
129  cylinder rCC (l1x2,l1y2,l1z2) (l1x2+(l1x3-l1x2)/2,l1y2+(l1y3-l1y2)/2,l1z2+(l1z3-l1z2)/2) color (Cr,Cg,Cb) decoration
130 
131 C2 : sphere rC (0,0,0) color (Cr,Cg,Cb)
132  cylinder rCC (0,0,0) (l2x1/2,l2y1/2,l2z1/2) color (Cr,Cg,Cb) decoration
133  cylinder rCC (0,0,0) (l2x2/2,l2y2/2,l2z2/2) color (Cr,Cg,Cb) decoration
134 
135 C3 : sphere rC (0,0,0) color (Cr,Cg,Cb)
136  cylinder rCC (0,0,0) (l3x1/2,l3y1/2,l3z1/2) color (Cr,Cg,Cb) decoration
137  cylinder rCC (0,0,0) (l3x2/2,l3y2/2,l3z2/2) color (Cr,Cg,Cb) decoration
138 
139 C4 : sphere rS (0,0,0) color (Sr,Sg,Sb)
140  sphere rS (l4x2,l4y2,l4z2) color (Sr,Sg,Sb)
141  cylinder rCS (0,0,0) (l4x1/2,l4y1/2,l4z1/2) color (Sr,Sg,Sb) decoration
142  cylinder rSS (0,0,0) (l4x2,l4y2,l4z2) color (Sr,Sg,Sb) decoration
143  cylinder rCS (l4x2,l4y2,l4z2) (l4x2+(l4x3-l4x2)/2,l4y2+(l4y3-l4y2)/2,l4z2+(l4z3-l4z2)/2) color (Sr,Sg,Sb) decoration
144 
145 C5 : sphere rC (0,0,0) color (Cr,Cg,Cb)
146  cylinder rCC (0,0,0) (l5x1/2,l5y1/2,l5z1/2) color (Cr,Cg,Cb) decoration
147  cylinder rCC (0,0,0) (l5x2/2,l5y2/2,l5z2/2) color (Cr,Cg,Cb) decoration
148 
149 C6 : sphere rC (0,0,0) color (Cr,Cg,Cb)
150  cylinder rCC (0,0,0) (l6x1/2,l6y1/2,l6z1/2) color (Cr,Cg,Cb) decoration
151  cylinder rCC (0,0,0) (l6x2/2,l6y2/2,l6z2/2) color (Cr,Cg,Cb) decoration
152 
153 [JOINTS]
154 
155 revolute: C1 (l1x1,l1y1,l1z1) (0,0,0)
156  C2 (0,0,0) (l2x2,l2y2,l2z2)
157 
158 revolute: C2 (l2x1,l2y1,l2z1) (0,0,0)
159  C3 (0,0,0) (l3x2,l3y2,l3z2)
160 
161 revolute: C3 (l3x1,l3y1,l3z1) (0,0,0)
162  C4 (l4x2,l4y2,l4z2) (l4x3,l4y3,l4z3)
163 
164 revolute: C4 (l4x1,l4y1,l4z1) (0,0,0)
165  C5 (0,0,0) (l5x2,l5y2,l5z2)
166 
167 revolute: C5 (l5x1,l5y1,l5z1) (0,0,0)
168  C6 (0,0,0) (l6x2,l6y2,l6z2)
169 
170 revolute: C6 (l6x1,l6y1,l6z1) (0,0,0)
171  C1 (l1x2,l1y2,l1z2) (l1x3,l1y3,l1z3)
172