Molecular graphics

Molecular graphics

Molecular graphics (MG) is the discipline and philosophy of studying molecules and their properties through graphical representation.[1] IUPAC limits the definition to representations on a "graphical display device".[2] Ever since Dalton's atoms and Kekulé's benzene, there has been a rich history of hand-drawn atoms and molecules, and these representations have had an important influence on modern molecular graphics. This article concentrates on the use of computers to create molecular graphics. Note, however, that many molecular graphics programs and systems have close coupling between the graphics and editing commands or calculations such as in molecular modelling.

Contents

Relation to molecular models

Fig. 1. Key: Hydrogen = white, carbon = grey, nitrogen = blue, oxygen = red, and phosphorus = orange.

There has been a long tradition of creating molecular models from physical materials. Perhaps the best known is Crick and Watson's model of DNA built from rods and planar sheets, but the most widely used approach is to represent all atoms and bonds explicitly using the "ball and stick" approach. This can demonstrate a wide range of properties, such as shape, relative size, and flexibility. Many chemistry courses expect that students will have access to ball and stick models. One goal of mainstream molecular graphics has been to represent the "ball and stick" model as realistically as possible and to couple this with calculations of molecular properties.

Figure 1 shows a small molecule (NH3CH2CH2C(OH)(PO3H)(PO3H)-), as drawn by the Jmol program. It is important to realize that the colors and shapes are purely a convention, as individual atoms are not colored, nor do they have hard surfaces. Bonds between atoms are also not rod-shaped.

Comparison of physical models with molecular graphics

Physical models and computer models have partially complementary strengths and weaknesses. Physical models can be used by those without access to a computer and now can be made cheaply out of plastic materials. Their tactile and visual aspects cannot be easily reproduced by computers (although haptic devices have occasionally been built). On a computer screen, the flexibility of molecules is also difficult to appreciate; illustrating the pseudorotation of cyclohexane is a good example of the value of mechanical models.

However, it is difficult to build large physical molecules, and all-atom physical models of even simple proteins could take weeks or months to build. Moreover, physical models are not robust and they decay over time. Molecular graphics is particularly valuable for representing global and local properties of molecules, such as electrostatic potential. Graphics can also be animated to represent molecular processes and chemical reactions, a feat that is not easy to reproduce physically.

History

Initially the rendering was on early Cathode ray tube screens or through plotters drawing on paper. Molecular structures have always been an attractive choice for developing new computer graphics tools, since the input data are easy to create and the results are usually highly appealing. The first example of MG was a display of a protein molecule (Project MAC, 1966) by Cyrus Levinthal and Robert Langridge. Among the milestones in high-performance MG was the work of Nelson Max in "realistic" rendering of macromolecules using reflecting spheres.

By about 1980 many laboratories both in academia and industry had recognized the power of the computer to analyse and predict the properties of molecules, especially in materials science and the pharmaceutical industry. The discipline was often called "molecular graphics" and in 1982 a group of academics and industrialists in the UK set up the Molecular Graphics Society (MGS). Initially much of the technology concentrated either on high-performance 3D graphics, including interactive rotation or 3D rendering of atoms as spheres (sometimes with radiosity). During the 1980s a number of programs for calculating molecular properties (such as molecular dynamics and quantum mechanics) became available and the term "molecular graphics" often included these. As a result the MGS has now changed its name to the Molecular Graphics and Modelling Society (MGMS).

The requirements of macromolecular crystallography also drove MG because the traditional techniques of physical model-building could not scale. Alwyn Jones' FRODO program (and later "O") were developed to overlay the molecular electron density determined from X-ray crystallography and the hypothetical molecular structure.

Art, science and technology in molecular graphics

Fig. 2. Image of hemagglutinin with alpha helices depicted as cylinders and the rest of the chain as silver coils. The individual protein atoms (several thousand) have been hidden. All of the non-hydrogen atoms in the two ligands (presumably sialic acid) have been shown near the top of the diagram. Key: Carbon = grey, oxygen = red, nitrogen = blue.

Both computer technology and graphic arts have contributed to molecular graphics. The development of structural biology in the 1950s led to a requirement to display molecules with thousands of atoms. The existing computer technology was limited in power, and in any case a naive depiction of all atoms left viewers overwhelmed. Most systems therefore used conventions where information was implicit or stylistic. Two vectors meeting at a point implied an atom or (in macromolecules) a complete residue (10-20 atoms).

The macromolecular approach was popularized by Dickerson and Geis' presentation of proteins and the graphic work of Jane Richardson through high-quality hand-drawn diagrams such as the "ribbon" representation. In this they strove to capture the intrinsic 'meaning' of the molecule. This search for the "messages in the molecule" has always accompanied the increasing power of computer graphics processing. Typically the depiction would concentrate on specific areas of the molecule (such as the active site) and this might have different colors or more detail in the number of explicit atoms or the type of depiction (e.g., spheres for atoms).

In some cases the limitations of technology have led to serendipitous methods for rendering. Most early graphics devices used vector graphics, which meant that rendering spheres and surfaces was impossible. Michael Connolly's program "MS" calculated points on the surface-accessible surface of a molecule, and the points were rendered as dots with good visibility using the new vector graphics technology, such as the Evans and Sutherland PS300 series. Thin sections ("slabs") through the structural display showed very clearly the complementarity of the surfaces for molecules binding to active sites, and the "Connolly surface" became a universal metaphor.

The relationship between the art and science of molecular graphics is shown in the exhibitions sponsored by the Molecular Graphics Society. Some exhibits are created with molecular graphics programs alone, while others are collages, or involve physical materials. An example from Mike Hann (1994), inspired by Magritte's painting Ceci n'est pas une pipe, uses an image of a salmeterol molecule.

"Ceci n'est pas une molecule," writes Mike Hann, "serves to remind us that all of the graphics images presented here are not molecules, not even pictures of molecules, but pictures of icons which we believe represent some aspects of the molecule's properties."

Space-filling models

Fig. 4. Space-filling model of formic acid. Key: Hydrogen = white, carbon = black, oxygen = red.

Fig. 4 is a "space-filling" representation of formic acid, where atoms are drawn to suggest the amount of space they occupy. This is necessarily an icon: in the quantum mechanical representation of molecules, there are only (positively charged) nuclei and a "cloud" of negative electrons. The electron cloud defines an approximate size for the molecule, though there can be no single precise definition of size. For many years the size of atoms has been approximated by mechanical models (CPK), where the atoms have been represented by plastic spheres whose radius (van der Waals radius) describes a sphere within which "most" of the electron density can be found. These spheres could be clicked together to show the steric aspects of the molecule rather than the positions of the nuclei. Fig. 4 shows the intricacy required to make sure that all spheres intersect correctly, and also demonstrates a reflective model.

Fig. 5. A molecule (zirconocene) where part (left) is rendered as ball-and-stick and part (right) as an isosurface.

Since the atomic radii (e.g. in Fig. 4) are only slightly less than the distance between bonded atoms, the iconic spheres intersect, and in the CPK models, this was achieved by planar truncations along the bonding directions, the section being circular. When raster graphics became affordable, one of the common approaches was to replicate CPK models in silico. It is relatively straightforward to calculate the circles of intersection, but more complex to represent a model with hidden surface removal. A useful side product is that a conventional value for the molecular volume can be calculated.

The use of spheres is often for convenience, being limited both by graphics libraries and the additional effort required to compute complete electronic density or other space-filling quantities. It is now relatively common to see images of surfaces that have been colored to show quantities such as electrostatic potential. Common surfaces in molecular visualization include solvent-accessible ("Lee-Richards") surfaces, solvent-excluded ("Connolly") surfaces, and isosurfaces. The isosurface in Fig. 5 appears to show the electrostatic potential, with blue colors being negative and red/yellow (near the metal) positive (there is no absolute convention of coloring, and red/positive, blue/negative are often reversed). Opaque isosurfaces do not allow the atoms to be seen and identified and it is not easy to deduce them. Because of this, isosurfaces are often drawn with a degree of transparency.

Technology

Early interactive molecular computer graphics systems were vector graphics machines, which used stroke-writing vector monitors, sometimes even oscilloscopes. The electron beam does not sweep left-and-right as in a raster display. The display hardware followed a sequential list of digital drawing instructions (the display list), directly drawing at an angle one stroke for each molecular bond. When the list was complete, drawing would begin again from the top of the list, so if the list was long (a large number of molecular bonds), the display would flicker heavily. Later vector displays could rotate complex structures with smooth motion, since the orientation of all of the coordinates in the display list could be changed by loading just a few numbers into rotation registers in the display unit, and the display unit would multiply all coordinates in the display list by the contents of these registers as the picture was drawn.

The early black-and white vector displays could somewhat distinguish for example a molecular structure from its surrounding electron density map for crystallographic structure solution work by drawing the molecule brighter than the map. Color display makes them easier to tell apart. During the 1970s two-color stroke-writing Penetron tubes were available, but not used in molecular computer graphics systems. In about 1980 Evans & Sutherland made the first practical full-color vector displays for molecular graphics, typically attached to an E&S PS-300 display. This early color tube was expensive, because its was originally engineered to withstand the shaking of a flight-simulator motion base.

Color raster graphics display of molecular models began around 1978 as seen in this paper by Porter [3] on spherical shading of atomic models. Early raster molecular graphics systems displayed static images that could take around a minute to generate. Dynamically rotating color raster molecular display phased in during 1982-1985 with the introduction of the Ikonas programmable raster display.

Molecular graphics has always pushed the limits of display technology, and has seen a number of cycles of integration and separation of compute-host and display. Early systems like Project MAC were bespoke and unique, but in the 1970s the MMS-X and similar systems used (relatively) low-cost terminals, such as the Tektronix 4014 series, often over dial-up lines to multi-user hosts. The devices could only display static pictures but were able to evangelize MG. In the late 1970s, it was possible for departments (such as crystallography) to afford their own hosts (e.g., PDP-11) and to attach a display (such as Evans & Sutherland's MPS) directly to the bus. The display list was kept on the host, and interactivity was good since updates were rapidly reflected in the display—at the cost of reducing most machines to a single-user system.

In the early 1980s, Evans & Sutherland (E&S) decoupled their PS300 display, which contained its own display information transformable through a dataflow architecture. Complex graphical objects could be downloaded over a serial line (e.g. 9600 baud) and then manipulated without impact on the host. The architecture was excellent for high performance display but very inconvenient for domain-specific calculations, such as electron-density fitting and energy calculations. Many crystallographers and modellers spent arduous months trying to fit such activities into this architecture.

The benefits for MG were considerable, but by the later 1980s, UNIX workstations such as Sun-3 with raster graphics (initially at a resolution of 256 by 256) had started to appear. Computer-assisted drug design in particular required raster graphics for the display of computed properties such as atomic charge and electrostatic potential. Although E&S had a high-end range of raster graphics (primarily aimed at the aerospace industry) they failed to respond to the low-end market challenge where single users, rather than engineering departments, bought workstations. As a result the market for MG displays passed to Silicon Graphics, coupled with the development of minisupercomputers (e.g., CONVEX and Alliant) which were affordable for well-supported MG laboratories. Silicon Graphics provided a graphics language, IrisGL, which was easier to use and more productive than the PS300 architecture. Commercial companies (e.g., Biosym, Polygen/MSI) ported their code to Silicon Graphics, and by the early 1990s, this was the "industry standard". Dial boxes were often used as control devices.

Stereoscopic displays were developed based on liquid crystal polarized spectacles, and while this had been very expensive on the PS300, it now became a commodity item. A common alternative was to add a polarizable screen to the front of the display and to provide viewers with extremely cheap spectacles with orthogonal polarization for separate eyes. With projectors such as Barco, it was possible to project stereoscopic display onto special silvered screens and supply an audience of hundreds with spectacles. In this way molecular graphics became universally known within large sectors of chemical and biochemical science, especially in the pharmaceutical industry. Because the backgrounds of many displays were black by default, it was common for modelling sessions and lectures to be held with almost all lighting turned off.

In the last decade almost all of this technology has become commoditized. IrisGL evolved to OpenGL so that molecular graphics can be run on any machine. In 1992, Roger Sayle released his RasMol program into the public domain. RasMol contained a very high-performance molecular renderer that ran on Unix/X Window, and Sayle later ported this to the Windows and Macintosh platforms. The Richardsons developed kinemages and the Mage software, which was also multi-platform. By specifying the chemical MIME type, molecular models could be served over the Internet, so that for the first time MG could be distributed at zero cost regardless of platform. In 1995, Birkbeck College's crystallography department used this to run "Principles of Protein Structure", the first multimedia course on the Internet, which reached 100 to 200 scientists.

Porin.qutemol.dl.png
Porin.qutemol.ao.png
Fig. 6. A molecule of Porin (protein) shown without ambient occlusion (left) and with (right). Advanced rendering effects can improve the comprehension of the 3D shape of a molecule.

MG continues to see innovation that balances technology and art, and currently zero-cost or open source programs such as PyMOL and Jmol have very wide use and acceptance.

Recently the wide spread diffusion of advanced graphics hardware, has improved the rendering capabilities of the visualization tools. The capabilities of current shading languages allow the inclusion of advanced graphic effects (like ambient occlusion, cast shadows and non-photorealistic rendering techniques) in the interactive visualization of molecules. These graphic effects, beside being eye candy, can improve the comprehension of the three dimensional shapes of the molecules. An example of the effects that can be achieved exploiting recent graphics hardware can be seen in the simple open source visualization system QuteMol.

Algorithms

Reference frames

Drawing molecules requires a transformation between molecular coordinates (usually, but not always, in Angstrom units) and the screen. Because many molecules are chiral it is essential that the handedness of the system (almost always right-handed) is preserved. In molecular graphics the origin (0, 0) is usually at the lower left, while in many computer systems the origin is at top left. If the z-coordinate is out of the screen (towards the viewer) the molecule will be referred to right-handed axes, while the screen display will be left-handed.

Molecular transformations normally require:

  • scaling of the display (but not the molecule).
  • translations of the molecule and objects on the screen.
  • rotations about points and lines.

Conformational changes (e.g. rotations about bonds) require rotation of one part of the molecule relative to another. The programmer must decide whether a transformation on the screen reflects a change of view or a change in the molecule or its reference frame.

Simple

Fig. 7. Stick model of caffeine drawn in Jmol.


In early displays only vectors could be drawn e.g. (Fig. 7) which are easy to draw because no rendering or hidden surface removal is required.

On vector machines the lines would be smooth but on raster devices Bresenham's algorithm is used (note the "jaggies" on some of the bonds, which can be largely removed with antialiasing software.)

Atoms can be drawn as circles, but these should be sorted so that those with the largest z-coordinates (nearest the screen) are drawn last. Although imperfect, this often gives a reasonably attractive display. Other simple tricks which do not include hidden surface algorithms are:

  • coloring each end of a bond with the same color as the atom to which it is attached (Fig. 7).
  • drawing less than the whole length of the bond (e.g. 10%-90%) to simulate the bond sticking out of a circle.
  • adding a small offset white circle within the circle for an atom to simulate reflection.

Typical pseudocode for creating Fig. 7 (to fit the molecule exactly to the screen):

// assume:
// atoms with x, y, z coordinates (Angstrom) and elementSymbol
// bonds with pointers/references to atoms at ends
// table of colors for elementTypes
// find limits of molecule in molecule coordinates as xMin, yMin, xMax, yMax
scale = min(xScreenMax/(xMax-xMin), yScreenMax/(yMax-yMin))
xOffset = -xMin * scale; yOffset = -yMin * scale
for (bond in $bonds) {
  atom0 = bond.getAtom(0)
  atom1 = bond.getAtom(1)
  x0 = xOffset+atom0.getX()*scale; y0 = yOffset+atom0.getY()*scale // (1)
  x1 = xOffset+atom1.getX()*scale; y1 = yOffset+atom1.getY()*scale // (2)
  x1 = atom1.getX();  y1 = atom1.getY()
  xMid = (x0 + x1) /2;  yMid = (y0 + y1) /2;
  color0 = ColorTable.getColor(atom0.getSymbol())
  drawLine (color0, x0, y0, xMid, yMid)
  color1 = ColorTable.getColor(atom1.getSymbol())
  drawLine (color1, x1, y1, xMid, yMid)
}

Note that this assumes the origin is in the bottom left corner of the screen, with Y up the screen. Many graphics systems have the origin at the top left, with Y down the screen. In this case the lines (1) and (2) should have the y coordinate generation as:

 y0 = yScreenMax -(yOffset+atom0.getY()*scale) // (1)
 y1 = yScreenMax -(yOffset+atom1.getY()*scale) // (2)

Changes of this sort change the handedness of the axes so it is easy to reverse the chirality of the displayed molecule unless care is taken.

Advanced

For greater realism and better comprehension of the 3D structure of a molecule many computer graphics algorithms can be used. For many years molecular graphics has stressed the capabilities of graphics hardware and has required hardware-specific approaches. With the increasing power of machines on the desktop, portability is more important and programs such as Jmol have advanced algorithms that do not rely on hardware. On the other hand recent graphics hardware is able to interactively render very complex molecule shapes with a quality that would not be possible with standard software techniques.

Chronology

Timeline

Developer(s) Approximate date Technology Comments
Crystallographers < 1960 Hand-drawn Crystal structures, with hidden atom and bond removal. Often clinographic projections.
Johnson, Motherwell ca 1970 Pen plotter ORTEP, PLUTO. Very widely deployed for publishing crystal structures.
Cyrus Levinthal, Bob Langridge, Ward, Stots [4] 1966 Project MAC display system, two-degree of freedom, spring-return velocity joystick for rotating the image. First protein display on screen. System for interactively building protein structures.
Barry [5] 1969 Link 300 computer with a dual trace oscilloscope display. Interactive molecular structure viewing system. Early examples of dynamic rotation, intensity depth·cueing, and side-by-side stereo. Early use of the small angle approximations (a = sin a, 1 = cos a) to speed up graphical rotation calculations.
Ortony [6] 1971 Designed a stereo viewer (British patent appl. 13844/70) for molecular computer graphics. Horizontal two-way (half-silvered) mirror combines images drawn on the upper and lower halves of a CRT. Crossed polarizers isolate the images to each eye.
Ortony [7] 1971 Light pen, knob. Interactive molecular structure viewing system. Select bond by turning another knob until desired bond lights up in sequence, a technique later used on the MMS-4 system below, or by picking with the light pen. Points in space are specified with a 3-D ”bug” under dynamic control.
Barry, Graesser, Marshall [8] 1971 CHEMAST: LINK 300 computer driving an oscilloscope. Two-axis joystick, similar to one used later by GRIP-75 (below). Interactive molecular structure viewing system. Structures dynamically rotated using the joystick.
Tountas and Katz [9] 1971 Adage AGT/50 display Interactive molecular structure viewing system. Mathematics of nested rotation and for laboratory-space rotation.
Perkins, Piper, Tattam, White [10] 1971 Honeywell DDP 516 computer, EAL TR48 analog computer, Lanelec oscilloscope, 7 linear potentiometers. Stereo. Interactive molecular structure viewing system.
Wright [11][12][13] 1972 GRIP-71 at UNC-CH: IBM System/360 Model 40 time-shared computer, IBM 2250 display, buttons, light pen, keyboard. Discrete manipulation and energy relaxation of protein structures. Program code became the foundation of the GRIP-75 system below.
Barry and North [14] 1972 Oxford Univ.: Ferranti Argus 500 computer, Ferranti model 30 display, keyboard, track ball, one knob. Stereo. Prototype large-molecule crystallographic structure solution system. Track ball rotates a bond, knob brightens the molecule vs. electron density map.
North, Ford, Watson Early 1970s Leeds Univ.: DEC PDP·11/40 computer, Hewlett-Packard display. 16 knobs, keyboard, spring-return joystick. Stereo. Prototype large-molecule crystallographic structure solution system. Six knobs rotate and translate a small molecule.
Barry, Bosshard, Ellis, Marshall, Fritch, Jacobi 1974 MMS-4:[15][16] Washington Univ. at St. Louis, Link 300 computer and an LDS-1 / Link 300 display, custom display modules. Rotation joystick, knobs. Stereo. Prototype large-molecule crystallographic structure solution system. Select bond to rotate by turning another knob until desired bond lights up in sequence.
Cohen and Feldmann [17] 1974 DEC PDP-10 computer, Adage display, push buttons, keyboard, knobs Prototype large-molecule crystallographic structure solution system.
Stellman [18] 1975 Princeton: PDP-10 computer, LDS-1 display, knobs Prototype large-molecule crystallographic structure solution system. Electron density map not shown; instead an "H Factor" figure of merit is updated as the molecular structure is manipulated.
Collins, Cotton, Hazen, Meyer, Morimoto 1975 CRYSNET ,[19] Texas A&M Univ. DEC PDP-11/40 computer, Vector General Series 3 display, knobs, keyboard. Stereo. Prototype large-molecule crystallographic structure solution system. Variety of viewing modes: rocking, spinning, and several stereo display modes.
Cornelius and Kraut 1976 (approx.) Univ, of Calif. at San Diego: DEC PDP-11/40 emulator (CalData 135), Evans and Sutherland Picture System display, keyboard, 6 knobs. Stereo. Prototype large-molecule crystallographic structure solution system.
(Yale Univ.) 1976 (approx.) PIGS: DEC PDP-11/70 computer, Evans and Sutherland Picture System 2 display, data tablet, knobs. Prototype large-molecule crystallographic structure solution system. The tablet was used for most interactions.
Feldmann and Porter 1976 NIH: DEC PDP—11/70 computer. Evans and Sutherland Picture System 2 display, knobs. Stereo. Interactive molecular structure viewing system. Intended to display interactively molecular data from the AMSOM – Atlas of Macromolecular Structure on Microfiche.[20]
Rosenberger et al. 1976 MMS-X:[21] Washington Univ. at St. Louis, TI 980B computer, Hewlett-Packard 1321A display, Beehive video terminal, custom display modules, pair of 3-D spring-return joysticks, knobs. Prototype (and later successful) large-molecule crystallographic structure solution system. Successor to the MMS-4 system above. The 3-D spring-return joysticks either translate and rotate the molecular structure for viewing or a molecular substructure for fitting, mode controlled by a toggle switch.
Britton, Lipscomb, Pique, Wright, Brooks 1977 GRIP-75 [13][22][23][24][25] at UNC-CH: Time-shared IBM System/360 Model 75 computer, DEC PDP 11/45 computer, Vector General Series 3 display, 3-D movement box from A.M. Noll and 3-D spring return joystick for substructure manipulation, Measurement Systems nested joystick, knobs, sliders, buttons, keyboard, light pen. First large-molecule crystallographic structure solution.[26]
Jones 1978 FRODO and RING [27][28] Max Planck Inst., Germany, RING: DEC PDP-11/40 and Siemens 4004 computers, Vector General 3404 display, 6 knobs. Large-molecule crystallographic structure solution. FRODO may have run on a DEC VAX-780 as a follow-on to RING.
Diamond 1978 Bilder [29] Cambridge, England, DEC PDP-11/50 computer, Evans and Sutherland Picture System display, tablet. Large-molecule crystallographic structure solution. All input is by data tablet. Molecular structures built on-line with ideal geometry. Later passes stretch bonds with idealization.
Langridge, White, Marshall Late 1970s Departmental systems (PDP-11, Tektronix displays or DEC-VT11, e.g. MMS-X) Mixture of commodity computing with early displays.
Davies, Hubbard Mid-1980s CHEM-X, HYDRA Laboratory systems with multicolor, raster and vector devices (Sigmex, PS300).
Biosym, Tripos, Polygen Mid-1980s PS300 and lower cost dumb terminals (VT200, SIGMEX) Commercial integrated modelling and display packages.
Silicon Graphics, Sun Late 1980s IRIS GL (UNIX) workstations Commodity-priced single-user workstations with stereoscopic display.
EMBL - WHAT IF 1989, 2000 Machine independent Nearly free, multifunctional, still fully supported, many free servers based on it
Sayle, Richardson 1992, 1993 RasMol, Kinemage Platform-independent MG.
MDL (van Vliet, Maffett, Adler, Holt) 1995–1998 Chime proprietary C++ ; free browser plugin for Mac (OS9) and PCs
ChemAxon 1998- MarvinSketch & MarvinView. MarvinSpace (2005) proprietary Java applet or stand-alone application.
Community efforts 2000- Jmol, PyMol, Avogadro, PDB Open-source Java applet or stand-alone application.
NOCH 2002- NOC Powerful and open source code molecular structure explorer
LION Bioscience / EMBL 2004- SRS 3D Free, open-source system based on Java3D. Integrates 3D structures with sequence and feature data (domains, SNPs, etc.).
San Diego Supercomputer Center 2006- Sirius Free for academic/non-profit institutions
Weizmann Institute of Science - Community efforts 2008- Proteopedia Collaborative, 3D wiki encyclopedia of proteins & other molecules

Electronic Richards Box Systems

Before computer graphics could do the job, mechanical methods were used to fit large molecules to their electron density maps. Using techniques of X-ray crystallography crystal of a substance were bombarded with x-rays, and the diffracted beams that came off were assembled by computer using a Fourier transform into a usually blurry 3-D image of the molecule, made visible by drawing contour circles around high electron density to produce a contoured electron density map.

In the earliest days contoured electron density maps were hand drawn on large plastic sheets. Sometimes bingo chips were placed on the plastic sheets where atoms were interpreted to be.

This was superseded by the Richards Box (see picture here) [1] in which an adjustable brass Kendrew molecular model was placed front of a 2-way mirror, behind which were plastic sheets of the electron density map. This optically superimposed the molecular model and the electron density map. The model was moved to within the contour lines of the superimposed map. Then, atomic coordinates were recorded using a plumb bob and a meter stick.

Computer graphics held out the hope of vastly speeding up this process, as well as giving a clearer view in many ways.

A noteworthy attempt to overcome the low speed of graphics displays of the time took place at Washington University in St. Louis. Dave Barry's group attempted to leapfrog the state of the art in graphics displays by making custom display hardware to display images complex enough for large-molecule crystallographic structure solution, fitting molecules to their electron-density maps. The MMS-4 (table above) display modules were slow and expensive, so a second generation of modules was produced for the MMS-X (table above) system.

The first large molecule whose atomic structure was partly determined on a molecular computer graphics system was Transfer RNA by Sung-Hou Kim's team in 1976.[30][31] after initial fitting on a mechanical Richards Box. The first large molecule whose atomic structure was entirely determined on a molecular computer graphics system is said to be neurotoxin A from venom of the Philippines sea snake, by Tsernoglou, Petsko, and Tu,[32] with a statement of being first [26] in 1977. The Richardson group published partial atomic structure results [33] of the protein superoxide dismutase the same year, in 1977. All of these were done using the GRIP-75 system.

Other structure fitting systems, FRODO, RING, Builder, MMS-X, etc. (table above) succeeded as well within three years and became dominant.

The reason that most of these systems succeeded in just those years, not earlier or later, and within a short timespan had to do with the arrival of commercial hardware that was up to the job. Two things were needed and arrived at about the same time. First, electron density maps are large and require either a computer with at least a 24-bit address space or a combination of a computer with a lesser 16-bit address space plus several years to overcome the difficulties of an address space that is smaller than the data. The second arrival was that of interactive computer graphics displays that were fast enough to display electron-density maps, whose contour circles require the display of numerous short vectors. The first such displays were the Vector General Series 3 and the Evans and Sutherland Picture System 2, MultiPicture System, and PS-300.

For further reading about early computer graphics systems for large-molecule crystallographic structure solutions see Duncan McRee.

Nowadays fitting of the molecular structure to the electron density map is largely automated by algorithms with computer graphics a guide to the process. An example is the XtalView XFit program.

See also

References

  1. ^ Dickerson, R.E.; Geis, I. (1969). The structure and action of proteins. Menlo Park, CA: W.A. Benjamin. 
  2. ^ IUPAC, Compendium of Chemical Terminology, 2nd ed. (the "Gold Book") (1997). Online corrected version:  (1997) "molecular graphics".
  3. ^ Porter TK, Spherical shading, ACM SIGGRAPH Computer Graphics, 12, 3 (August 1978), 282-285.
  4. ^ Levinthal, C. 1966. Molecular Model-building by Computer. Scientific American, 214, #6 (June), 42-52.
  5. ^ Barry, C. D., Ellis, R. A., Graesser, S. M., and Marshall, G. R. 1969. Display and Manipulation in Three Dimensions. Pertinent Comcepts in Computer Graphics, Univ. of Ill. Press, 104-153.
  6. ^ Ortony, A. 1971a. A System for Stereo Viewing. Computing Journal, 14, #2 (May), 140-144.. Also appears in: Conference on Displays, Institution of Electrical Engineers Conf. Pub. No. 80 (7–10 September 1971), C. Baldwin Ltd., 225-232.
  7. ^ Ortony, A. 1971b. Interactive Stereographics Conference on Displays, Institution of Electrical Engineers Conf. Pub. No. 80 (7–10 September), C. Baldwin Ltd., 185-193.
  8. ^ Barry, C. D., Ellis, R. A., Graesser, S. M., and Marshall, G. R. 1971. CHEMAST: A Computer Program for Modeling Molecular Structures. Proc. 1971 IFIP, 1552-1558.
  9. ^ Tountas, C. and Katz, L. 1971. Interactive Graphics in Molecular Biology. Real·time Three-dimensional Rotations of Images and Image Fragments. Proc. Summer Computer Simulation Conf., 1, 241-247.
  10. ^ Perkins, W. J., Piper, E. A., Tattam, F. G., and White, J. C. 1971. Interactive Stereoscopic Computer Displays for Biomedical Research. Computers and Biomedical Research, 4, #3, 249-261.
  11. ^ Wright, W. V. 1972a. An Interactive Computer Graphic System for Molecular Studies. Ph. D. Dissertation, University of North Carolina, Chapel Hill, North Carolina.
  12. ^ Wright, W. V. 1972b. The Two-Dimensional Interface of an Interactive System for Molecular Studies. ACM SIGPLAN notices, 1, #10 (October), 76-85.
  13. ^ a b Brooks FP Jr. The Computer "Scientist" as Toolsmith: Studies in Interactive Computer Graphics. Proc. IFIP, 625-634 (1977).
  14. ^ Barry, C .D. and North, A. C. T. 1972. The Use of a Computer-Controlled Display System in the Study of Molecular Conformations. Cold Springe Harbor Symp. Quant. Biol., 3p, 577.
  15. ^ Barry, C. D., Bosshard, R. A., Ellis, R. A., and Marshall, G. R. 1974. Evolving Macromodular Modeling System. Federation of American Societiee for Experimental Biology, Federation Proc., 33, #12 (December), 2368-2374.
  16. ^ Fritch, J. M., Ellis, R. A., Jacobi T. H., and Marshall, G. R. 1975. A Macromolecular Graphics System for Protein Structure Research. Computers and Graphics, 1, #2/3 (September), 271-278.
  17. ^ Cohen, G. H. and Feldmann, R. J. 1974. MAP - An Interactive Graphics Computer Program for the Manipulation and Fitting of Protein Molecules to Electron Density Maps. Am. Crystallography. Assoc. Spring 23, (Abstr.).
  18. ^ Stellman, S. D. 1975. Three-dimensional Interactive Graphics in X-ray Crystalographic Analysis. Computers and Graphics, #2/3 (September), 279-288.
  19. ^ Collins, D. M., Cotton, F. A., Hazen, E. E., Jr., Meyer, E. F., Jr., and Morimoto, C. N. 1975. Protein Crystal Structures: Quicker, Cheaper Approaches. Science, 190, #4219 (12 December), 1047-1053.
  20. ^ Feldmann, R. J. 1976. AMSOM – Atlas of Macromolecular Structure on Microfiche.. Maryland: Tracor Jitco Inc.
  21. ^ Rosenberger, F. U., et al. 1976. Extracts from 1976 NIH Annual Report. Technical Memorandum No. 230, Computer Science Laboratory, Washington University, St. Louis, Missouri.
  22. ^ Lipscomb, JS. Three-dimensional cues for a molecular computer graphics system. Ph.D. Dissertation, University of North Carolina at Chapel Hill, North Carolina. (1981)
  23. ^ Britton E, Lipscomb JS, Pique ME, Wright, WV, Brooks FP Jr, Pique ME. The GRIP-75 Man-machine Interface. ACM SIGGRAPH Video Review, (4), (Aug. 1981).
  24. ^ Britton, E. G. 1977. A Methodology for the Ergonomic Design of Interactive Computer Graphics Systems, and its Application to Crystallography. Ph. D. Dissertation, University of North Carolina, Chapel Hill, North Carolina.
  25. ^ Pique, M. E. 1980. Nested Dynamic Rotations for Computer Graphics. M. S. Thesis, University of North Carolina, Chapel Hill, North Carolina..
  26. ^ a b Tsernoglou D, Petsko GA, Tu AT. Protein sequencing by computer graphics. Biochem. and Biophys. ACTA, 491(2):605–608 (Apr 25, 1977).Link to abstract
  27. ^ Jones, T. A. 1978a. A Graphics Model Building and Refinement System for Macromolecules. Journal of Applied Crystallography, 11, 268-272.
  28. ^ Jones, T. A. 1978b. The RING [user manual]. Max-Planck-Institut fur Biochemie, 8033 Martinsried bei Muchen, Germany.
  29. ^ Diamond, R. 1978. Bilder. A computer graphics program for bipolymers and its application to interpretation of structure of tobacco mosaic virus protein disks at 2-A resolution. Proc. International Union of Pure and Applied Biochemistry: International Symposium on Structure, Conformation, Function, and Evolution. Madras, India, (4 January), Pergamon Press.
  30. ^ Sussman, JL. and Kim, SH. Idealized atomic coordinates of yeast phenylalanine Transfer RNA. Biochem. and Biophys. Res. Commun., 68, 1, (1976) 89-96.
  31. ^ Sussman, JL. and Kim, SH. Three-dimensional structure of a Transfer RNA common in two crystal forms. Science, 192, 4242 (20 May 1976), 835-838.
  32. ^ Tsernoglou, D. and Petsko, GA., Three-dimensional structure of neurotoxin A from venom of the Philippines sea snake. Proc. National Academy of Sciences USA, 74, 3 (March 1977), 971-974.
  33. ^ Richardson D.C , 1977. Three-dimensional structure of cu, zn superoxide dismutase. Superoxide and superoxide disumases. Ed. By A.M. Michelson, J.M. McCord, and I. Fridoivich. London, NY: Academic Press.

External links


Wikimedia Foundation. 2010.

Игры ⚽ Нужен реферат?

Look at other dictionaries:

  • Journal of Molecular Graphics & Modelling —   Titre abrégé J. Mol. Graph. Discipline Chimie …   Wikipédia en Français

  • Molecular design software — is software for molecular modeling, that provides special support for developing molecular models de novo. In contrast to the usual molecular modeling programs such as the molecular dynamics and quantum chemistry programs, such software directly… …   Wikipedia

  • Molecular dynamics — (MD) is a computer simulation of physical movements of atoms and molecules. The atoms and molecules are allowed to interact for a period of time, giving a view of the motion of the atoms. In the most common version, the trajectories of molecules… …   Wikipedia

  • Molecular Biology and Evolution —   Titre abrégé Mol. Biol. Evol. Discipline Biologie Langue …   Wikipédia en Français

  • Molecular Cell —   Titre abrégé Mol. Cell. Discipline Biochimie, biologie moléculaire, biologie cellulaire …   Wikipédia en Français

  • Molecular and Cellular Biology —   Titre abrégé Mol. Cell. Biol. Discipline Biochimie, biologie moléculaire, biol …   Wikipédia en Français

  • Molecular BioSystems —   Titre abrégé Mol. Biosyst. Discipline Biochimie Langue …   Wikipédia en Français

  • Molecular model — A molecular model, in this article, is a physical model that represents molecules and their processes. The creation of mathematical models of molecular properties and behaviour is molecular modelling, and their graphical depiction is molecular… …   Wikipedia

  • Molecular modelling — The backbone dihedral angles are included in the molecular model of a protein. Modelling of ionic li …   Wikipedia

  • Molecular mechanics — A force field is used to minimize the bond stretching energy of this ethane molecule. Molecular mechanics uses Newtonian mechanics to model molecular systems. The potential energy of all systems in molecular mechanics is calculated using force… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”