hw6 - on the x-axis to keep things simple (l12,0,0). For...

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
Biological Computation 20.181 Homework 6 In this homework, you will begin the basic problem of building 3D molecules from an internal reference frame: the rotational angles of each rotatable chemical bond (aka, "torsion space"). However, we will approach these problems from basic geometry. Building a polymer. We need to determine what information we will need from the chemistry of a molecule to determine its 3D structure. For a 'polymer' of a single atom, there are no real constraints -- we can put the atom anywhere we choose. For simplicity, let's put it at the origin of our reference frame (0,0,0). If we have two atoms, we can still put the first one at the origin, but the position of the second atom is constrained by the bond length between atoms 1 and 2 (l12). We can choose to place the second atom
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: on the x-axis to keep things simple (l12,0,0). For the third atom, we need to consider the bond length to atom 2 (l23) and the bond angle formed by atoms 1-3 (a123). We still have some freedom in that we can choose to put these atoms in the x-y plane. If we do so, what will the position of the third atom be? Write a Python program which returns coordinates for three atoms in the following form: (0,0,0) (?,0,0) (?,?,?) given two bond lengths (l12 and l23) and one bond angle (a123), which it will read from a file in this format: l12 1.5 l23 1.45 a123 119.5 Before the next lecture start to think about what information you will need to place a fourth atom. .....
View Full Document

This note was uploaded on 11/11/2011 for the course BIO 20.010j taught by Professor Lindagriffith during the Spring '06 term at MIT.

Ask a homework question - tutors are online