Interactive Computer Program: Packaging DNA into Chromosomes Xiaoli Yang 1, Yifan Cai 1 and Charles Tseng 2 1Department of Electrical and Computer Engineering 2Department of Biological Sciences Purdue University Calumet Hammond, IN, USA Abstract - As part of the interactive program for teaching and serve as a model for STEM (science, technology, engineering, learning genetics, the module on packaging DNA into and mathematics) education via distance learning. chromosomes involves the simultaneous coordination of eyes, mind, and hands for visualization, cognitive feedback, and 2 Model Development manipulation, respectively. Computer modeling of various chromatin structures during packaging is based on OpenTK- Models of various structures were developed based on OpenGL on .Net Platform, which is coupled with an inquiry the following system: OpenTK-OpenGL on .Net Platform. based content design to enhance the efficiency of teaching and The Open Tool Kit (OpenTK) is a free project that allows learning. The prototype has been successfully tested in a developers to use OpenGL, OpenGL|ES, OpenCL, and genetics class at Purdue University Calumet. It should also be OpenAL APIs from a managed language (e.g. VB.NET). applicable to a number of undergraduate biology courses . Features include: • Written in cross-platform C# and usable by all managed Keywords: DNA, Chromosomes, Modeling, Computer languages (F#, Boo, VB.Net, C++/CLI). Program • Consistent, strongly typed bindings, suitable for RAD development. 1 Introduction • Usable standing alone or integrated with Windows.Forms, GTK#, and WPF. From its central role in real-life forensic investigations • Cross-platform binaries that are portable on .Net and to being the basis of major biotechnological applications in Mono without recompilation. medicine, agriculture, and the environment, DNA based • Wide platform support: Windows, Linux, and Mac OS genetics is an essential discipline in the life sciences. As X, with iPhone port in process. fascinating as the subject is, however, teaching and learning genetics has often been fraught with difficulty [1-3]. Confronted with intricate molecular structures, complex 2.1 3D model of double helix DNA packaging schemes, and elaborate mechanisms of action, both DNA is modeled as a double helix. The model is teacher and student are frequently at a loss – the teacher in specified by l, the length of the helix, r, the radius of the helix, how to convey this material in a clear and understandable and w and h, the width and the thickness of one strand of the way, and the student in how to assimilate all the information double helix, respectively (Fig. 1). These parameters generate usefully. To be sure, the abstract and intangible nature of a group of points, which are used to construct the DNA much of the material is the source of the problem. model. The points are linked together to form a sketch of the Traditional methods of teaching genetics, employing double helix. After shading the sketch, a 3D DNA model is classroom lectures, textbook readings, homework created. The double helix model is calculated during runtime assignments, and laboratory exercises, have not proven to be based on the equations below: very effective [4, 5]. Recently, efforts have been made to integrate computer visualization technologies into pedagogy to enhance the learning process [6-8]. Current computer- , 0 (1) based tools, however, do not stress cognitive feedback in their sin designs. The present paper describes an innovative approach cos to teaching and learning genetics, in which students can Where t is the length variable along x-axis, r the radius of visualize a real-time, interactive DNA model, as well as the helix. the angle increment, controlling the smoothness actively control the dynamic process of packaging DNA into of the helix. We chose from the experiment to make the a compact metaphase chromosome. model smooth. determines 6 the initial angle. We chose The objectives of the program are to 1) develop, as part from the experiment to generate double helical of a web-based interactive program, a DNA packaging shapes. From the above equations, module suitable for a wide range of college courses and 2) 0 45 0, ! determine the Cartesian coordinates x, y and ) z in 3D space. , , 1 ( * , 1, 1 By calculating the position of all the points, a helical line & ) can be generated (Fig. 2, left). The quadrupling of the line $ 3 ( * , 3, 3 (2) (Fig. 2, right) is generated by replicating the original line four -./-* ) 0./0* 1./1* , ( , , 2 times. % * * * * $ -3/-4 ) 03/04 13/14 #, * ( * , * , * 2 The result is a complete DNA model (Fig. 5). Fig. 1. Parameters for a helix (perpendicular views) Fig. 2. Left: helical line; Right: doubling of line Fig. 5. Screen snapshot of DNA model from the program After further duplicating the strand with a different value and shading the sketch, two DNA strands of different 2.3 3D model of histone octomer colors are created (Fig. 3). The histone octomer is represented by an elongated ball, which is described by the following equations: sin 5 cos Fig. 3. DNA double helix structure with shading 8 (3) cos 5 7 * , 5 9 0 2.2 3D model of nucleotide bases 6 8 cos 5 * , 5 : 0 We use line segments (cuboid) of different colors to determine sin;sin"the Cartesian coordinate x, y and represent DNA bases. The points that form the DNA strands z in ,a 3 dimensional, space, where r is the radius of the (determined above) are used to calculate the points histone; the angle between the diameter and z-axis, representing the line segments (bases). Assume that p1 (x1, 5 ; the angle between the projection of the y1, z1) and p2 (x2, y2, z2) are corresponding points on = = different strands generated by the same t value, but different diameter5 <7 * ,on*> the plane and the x-axis. and h the height of the elongated ball (Fig. 6). values, p3 (x3, y3, z3) and p4 (x4, y4, z4) are the points , 0, ! "next to p1 and p2, respectively, w is the length of the side, pc1 is the midpoint between p1 and p2, and pc2 is the midpoint between p3 and p4 (Fig. 4). Fig. 6. Sphere coordinate system Fig. 4. Parameters used to calculate the position and shape of These points generate a sketch of the histone octomer. bases. Shading the sketch with a color completes the model (Fig. 7). Then all 8 points needed to describe a cuboid can be calculated as follows. implemented this function by adjusting camera position when we developed the model. The camera was located on the surface of sphere with the target at the center of sphere, so that the distance between camera and target never changes. In Fig. 7. Left: sketch of histone octomer; other words, the size of the target remains the same, so that Right: shaded histone octomer the model size does not change. Camera position (a point on the surface of sphere) is described by φ , θ and r, where r is the radius of sphere. While the value of r never changes, φ 2.4 DNA wrapping-formation of core and θ are variable. Changing φ and θ changes the camera nucleosome position, and they are changed by moving the mouse. Moving the mouse produces component values dx and dy along the x- In this step, DNA is simplified as a line, which can be and y-axis, respectively. Therefore, mapping dx and dy to φ wrapped around the histone octomer. Fig. 8 shows how to and θ is an effective way to adjust camera position. calculate the position of the binding points. 3 The Program Content Design The design of module focuses on inquiry based methods with cognitive feedback and interactive experiences as important components [9]. In every section, a question is followed by observations and measurements, hands-on experiments, and conclusions. In each of the learning steps, dynamic models of DNA molecules, chromatin fibers, and metaphase chromosomes are presented for interaction through visualization, cognition, and operation. Completion of the Fig. 8. Calculation for DNA-protein binding position program requires comprehension of the entire concept and The histone octomer is projected on the xz -plane as a thus ensures the success of the learning experience. The circle. Assume that O is the center of the circle, P is outside computer-based content is summarized below: the circle, r is the radius of the circle, dx and dy are the differences between O and P in x and y components, 3.1 DNA and chromosomes in prokaryotes and respectively, D is the distance between P and O, Pb is the eukaryotes binding point, is the angle between line P0 and the vertical line, and is ? the angle between line P0 and the line PbO . Inside the cell, DNA molecules are packaged, with helped Then and @ can be calculated as follows: of proteins, into thread-like structures called chromosomes. In prokaryotes (such as bacteria), the chromosomal DNA, ? @ when open, is often circular. The total length of a bacterial arcsin C (4) chromosomal DNA (e.g., E. coli DNA) may be a thousand F0 times longer than the cell that contains it. Little is known The value of is D arctan F- about the packaging of bacterial DNA, although a few major (5) DNA regions anchored by proteins at specific sites in the cell Finally, Pb is represented by ( x,y ), where 70H 7 7 D have been noted. In eukaryotes (such as animals and plants), DNA (6) sin molecules are linear. Each eukaryotic species has a fixed When DNA binds to a histone octomer, it starts to wrap cos number of chromosomes. For diploid species (species with 2 around the octomer. After Pb is determined, the points on the sets of chromosomes, one from each parent), chromosomes spiral nearest to it can be calculated (Fig.