Rendering volumetric data in molecular systems

A raster-based computer graphics method of imaging volumetric data (data sampled on a three-dimensional grid) has been added to an existing molecular rendering program. Any molecular property expressed as a 3D grid of scalar data values can be displayed around and inside of the molecular van der Waals surface as a collection of colored clouds and transparent surfaces. Different antical properties, such as color and opacity, are assigned to ranges of the molecular property and are rendered by a ray-tracing technique simulating the shading and shadowing of real objects, making the final images highly eadable. We have found useful applications in the visualization of crystallographic electron density, electrostatic potentials and protein active sites.     

The molecular surface package

The Molecular Surface Package is a reimplementation, in C, of a set of earlier FORTRAN programs for computing analytical molecular surfaces, areas, volumes, polyhedral molecular surfaces, and surface curvatures. The software does not do interactive molecular graphics, but it will produce pixel maps of smooth molecular surfaces. The polyhedral molecular surfaces are suited to display on graphics systems with real-time rendering of polyhedra.     

Realistic Real-Time Outdoor Rendering in Augmented Reality

Realistic rendering techniques of outdoor Augmented Reality (AR) has been an attractive topic since the last two decades considering the sizeable amount of publications in computer graphics. Realistic virtual objects in outdoor rendering AR systems require sophisticated effects such as: shadows, daylight and interactions between sky colours and virtual as well as real objects. A few realistic rendering techniques have been designed to overcome this obstacle, most of which are related to non real-time rendering. However, the problem still remains, especially in outdoor rendering. This paper proposed a much newer, unique technique to achieve realistic real-time outdoor rendering, while taking into account the interaction between sky colours and objects in AR systems with respect to shadows in any specific location, date and time. This approach involves three main phases, which cover different outdoor AR rendering requirements. Firstly, sky colour was generated with respect to the position of the sun. Second step involves the shadow generation algorithm, Z-Partitioning: Gaussian and Fog Shadow Maps (Z-GaF Shadow Maps). Lastly, a technique to integrate sky colours and shadows through its effects on virtual objects in the AR system, is introduced. The experimental results reveal that the proposed technique has significantly improved the realism of real-time outdoor AR rendering, thus solving the problem of realistic AR systems.     

Photorealistic image generation of molecular structure on PC screen using the ray-tracing technique

A PC version of three-dimensional molecular graphics package has been developed to run under MS-DOS environment on IBM PC-compatible computers equipped with a VGA graphics board. The program consists of two parts: a menu-driven interactive system module in EGA mode, and a ray-tracing module in VGA mode. In the 256-color VGA mode, ray-tracing images are represented with a 4-color map, with 64 levels for each color: 32 levels of illuminance and 32 levels of saturation. Molecular structure can be analyzed along various directions with various light sources. Ray-tracing images are also represented in a 16-color EGA mode with the half-toning method, which can display 76 gray levels for each color. To obtain good photo-realistic images in an efficient way, we have used two light sources, with an intensity ratio of 7:3, which are located in front of the top right and bottom left corners of the screen.     

Electron microscope tomography of Balbiani Ring hnRNP substructure

Three-dimensional (3-D) reconstructions, by electron microscope tomography, of selectively stained, contrast enhanced Balbiani Ring (BR) hnRNP granules reveal a complex spatial arrangement of RNA-rich domains. This particulate substructure was examined by volume rendering computer graphics. Modeling the arrangement of RNA-rich domains is made difficult by apparent structural flexibility and/or heterogeneity of composition. Formulation of a consensus 3-D arrangement of RNA-rich domains will require an expanded data base of reconstructed BR granules and the development of new image manipulation and analysis techniques. This study demonstrates the potential for ultra-structural cell biology of combining several new techniques: selective nucleic acid staining, electron spectroscopic imaging to enhance contrast, electron microscope tomography and volume rendering computer graphics.Abbreviations BR Balbiani Ring - EMT electron microscope tomography - ESI electron spectroscopic imaging - hnRNP heterogeneous nuclear ribonucleoprotein - OA-B osmium ammine-B - kb kilobases by P.B. Moens     

Texture mapping: A new tool for molecular graphics

The real-time texture mapping capabilities of modern graphics workstations are explored with respect to their applications in a variety of relevant scenarios in interactive molecular modeling techniques. The common usage of texture mapping to reduce geometric complexity while enhancing realism is extended, opening new ways to visualize large amounts of molecular data in a comprehensive fashion. Thus, texture mapping may be employed to (1) display and filter multichannel information of structural properties on molecular surfaces, (2) improve the quality and accuracy of highly complex isodensity contours, (3) increase the rendering speed of space-filling atomic representations by two orders of magnitude and (4) apply volume-rendering techniques to large, three-dimensional density distributions in real time. Implementation of these novel techniques requires only moderate modifications or extensions to existing molecular modeling applications.     

ProteinShader: illustrative rendering of macromolecules

Background

Cartoon-style illustrative renderings of proteins can help clarify structural features that are obscured by space filling or balls and sticks style models, and recent advances in programmable graphics cards offer many new opportunities for improving illustrative renderings.

Results

The ProteinShader program, a new tool for macromolecular visualization, uses information from Protein Data Bank files to produce illustrative renderings of proteins that approximate what an artist might create by hand using pen and ink. A combination of Hermite and spherical linear interpolation is used to draw smooth, gradually rotating three-dimensional tubes and ribbons with a repeating pattern of texture coordinates, which allows the application of texture mapping, real-time halftoning, and smooth edge lines. This free platform-independent open-source program is written primarily in Java, but also makes extensive use of the OpenGL Shading Language to modify the graphics pipeline.

Conclusion

By programming to the graphics processor unit, ProteinShader is able to produce high quality images and illustrative rendering effects in real-time. The main feature that distinguishes ProteinShader from other free molecular visualization tools is its use of texture mapping techniques that allow two-dimensional images to be mapped onto the curved three-dimensional surfaces of ribbons and tubes with minimum distortion of the images.     

Game On,Science - How Video Game Technology May Help Biologists Tackle Visualization Challenges

The video games industry develops ever more advanced technologies to improve rendering, image quality, ergonomics and user experience of their creations providing very simple to use tools to design new games. In the molecular sciences, only a small number of experts with specialized know-how are able to design interactive visualization applications, typically static computer programs that cannot easily be modified. Are there lessons to be learned from video games? Could their technology help us explore new molecular graphics ideas and render graphics developments accessible to non-specialists? This approach points to an extension of open computer programs, not only providing access to the source code, but also delivering an easily modifiable and extensible scientific research tool. In this work, we will explore these questions using the Unity3D game engine to develop and prototype a biological network and molecular visualization application for subsequent use in research or education. We have compared several routines to represent spheres and links between them, using either built-in Unity3D features or our own implementation. These developments resulted in a stand-alone viewer capable of displaying molecular structures, surfaces, animated electrostatic field lines and biological networks with powerful, artistic and illustrative rendering methods. We consider this work as a proof of principle demonstrating that the functionalities of classical viewers and more advanced novel features could be implemented in substantially less time and with less development effort. Our prototype is easily modifiable and extensible and may serve others as starting point and platform for their developments. A webserver example, standalone versions for MacOS X, Linux and Windows, source code, screen shots, videos and documentation are available at the address: http://unitymol.sourceforge.net/.     

Dynamic modeling of chemical reactions: the diels-alder cycloaddition

In order to visualize the stereochemical aspects of the Diels-Alder model cycloaddition of ethylene on butadiene, we developed a computer graphics animated model. The structural data base was deduced from MINDO/3 calculations, and the application program makes it possible to display in detail the different steps of the reaction mechanism. The scope of the application has been enlarged by a similar representation of the Diels-Alder cycloaddition of bis (methylene)—2,3 bicyclo [2.2.1]heptane to ethylene. Both of these examples suggest that molecular graphics is an ideal tool for visualizing and understanding the stereochemistry of complex chemical reactions.     

Role of network branching in eliciting differential short-term signaling responses in the hypersensitive epidermal growth factor receptor mutants implicated in lung cancer

We study the effects of EGFR inhibition in wild-type and mutant cell lines upon tyrosine kinase inhibitor TKI treatment through a systems level deterministic and spatially homogeneous model to help characterize the hypersensitive response of the cancer cell lines harboring constitutively active mutant kinases to inhibitor treatment. By introducing a molecularly resolved branched network systems model (the molecular resolution is introduced for EGFR reactions and interactions in order to distinguish differences in activation between wild-type and mutants), we are able to quantify differences in (1) short-term signaling in downstream ERK and Akt activation, (2) the changes in the cellular inhibition EC50 associated with receptor phosphorylation (i.e., 50% inhibition of receptor phosphorylation in the cellular context), and (3) EC50 for the inhibition of activated downstream markers ERK-(p) and Akt-(p), where (p) denotes phosphorylated, upon treatment with the inhibitors in cell lines carrying both wild-type and mutant forms of the receptor. Using the branched signaling model, we illustrate a possible mechanism for preferential Akt activation in the cell lines harboring the oncogenic mutants of EGFR implicated in non-small-cell lung cancer and the enhanced efficacy of the inhibitor erlotinib especially in ablating the cellular Akt-(p) response. Using a simple phenomenological model to describe the effect of Akt activation on cellular decisions, we discuss how this preferential Akt activation is conducive to cellular oncogene addiction and how its disruption can lead to dramatic apoptotic response and hence remarkable inhibitor efficacies. We also identify key network nodes of our branched signaling model through sensitivity analysis as those rendering the network hypersensitive to enhanced ERK-(p) and Akt-(p); intriguingly, the identified nodes have a strong correlation with species implicated in oncogenic transformations in human cancers as well as in drug resistance mechanisms identified for the inhibitors in non-small-cell lung cancer therapy.     

Multiplexed-Replica Exchange Molecular Dynamics Method for Protein Folding Simulation

Simulating protein folding thermodynamics starting purely from a protein sequence is a grand challenge of computational biology. Here, we present an algorithm to calculate a canonical distribution from molecular dynamics simulation of protein folding. This algorithm is based on the replica exchange method where the kinetic trapping problem is overcome by exchanging noninteracting replicas simulated at different temperatures. Our algorithm uses multiplexed-replicas with a number of independent molecular dynamics runs at each temperature. Exchanges of configurations between these multiplexed-replicas are also tried, rendering the algorithm applicable to large-scale distributed computing (i.e., highly heterogeneous parallel computers with processors having different computational power). We demonstrate the enhanced sampling of this algorithm by simulating the folding thermodynamics of a 23 amino acid miniprotein. We show that better convergence is achieved compared to constant temperature molecular dynamics simulation, with an efficient scaling to large number of computer processors. Indeed, this enhanced sampling results in (to our knowledge) the first example of a replica exchange algorithm that samples a folded structure starting from a completely unfolded state.     

Application of laser-capture microdissection to analysis of gene expression in the testis

The isolation and molecular analysis of highly purified cell populations from complex, heterogeneous tissues has been a challenge for many years. Spermatogenesis in the testis is a particularly difficult process to study given the unique multiple cellular associations within the seminiferous epithelium, making the isolation of specific cell types difficult. Laser-capture microdissection (LCM) is a recently developed technique that enables the isolation of individual cell populations from complex tissues. This technology has enhanced our ability to directly examine gene expression in enriched testicular cell populations by routine methods of gene expression analysis, such as real-time RT-PCR, differential display, and gene microarrays. The application of LCM has however introduced methodological hurdles that have not been encountered with more conventional molecular analyses of whole tissue. In particular, tissue handling (i.e. fixation, storage, and staining), consumables (e.g. slide choice), staining reagents (conventional H&E vs. fluorescence), extraction methods, and downstream applications have all required re-optimisation to facilitate differential gene expression analysis using the small amounts of material obtained using LCM. This review will discuss three critical issues that are essential for successful procurement of cells from testicular tissue sections; tissue morphology, capture success, and maintenance of molecular integrity. The importance of these issues will be discussed with specific reference to the two most commonly used LCM systems; the Arcturus PixCell IIe and PALM systems. The rat testis will be used as a model, and emphasis will be placed on issues of tissue handling, processing, and staining methods, including the application of fluorescence techniques to assist in the identification of cells of interest for the purposes of mRNA expression analysis.     

Variations of Surface Areas and Volumes in Distinct Molecular Surfaces of Biomolecules

The different surfaces usually defined in molecular modeling are explored focusing the attention on comparing the two common surfaces which incorporate the presence of the solvent: the Accessibility Molecular Surface (AMS) and the Richards Molecular Surface (RMS). The information contained in these surfaces is analyzed by studying the changes associated to the systematic variation of the radius of the solvent, represented by a probe sphere. An approximation to the effective fractal dimension associated with the RMS is also obtained. Four biomolecules, hyaluronic acid (HYA) and three proteins, crambin (CRN), lysozyme (LYZ), and thaumatin-I (THI), with sizes ranging from around 100 to 1500 atoms, have been chosen for this study. Graphical ray-traced renderings of these biomolecules, obtained by adapting public domain ray-tracing software are presented.     

Advancing simulations of biological materials: applications of coarse-grained models on graphics processing unit hardware

The timescales of biological processes, primarily those inherent to the molecular mechanisms of disease, are long (>μs) and involve complex interactions of systems consisting of many atoms (>10⁶). Simulating these systems requires an advanced computational approach, and as such, coarse-grained (CG) models have been developed and highly optimised for accelerator hardware, primarily graphics processing units (GPUs). In this review, I discuss the implementation of CG models for biologically relevant systems, and show how such models can be optimised and perform well on GPU-accelerated hardware. Several examples of GPU implementations of CG models for both molecular dynamics and Monte Carlo simulations on purely GPU and hybrid CPU/GPU architectures are presented. Both the hardware and algorithmic limitations of various models, which depend greatly on the application of interest, are discussed.     

Personal computer-based visualization of three-dimensional scalar and vector fields: An application to molecular graphics

A desktop PC-based graphics package, UNIVIS, for visualization of three-dimensional numerical data is described. Apart from routine molecular model visualization, the package provides for a host of other features such as extraction of various surfaces, planar cross-sections of the three-dimensional data, and property texturing. Fast rendering and transparency are the strengths of the present package. These features are comprehensively discussed. The salient features of UNIVIS are presented in the form of visualization of a variety of molecular properties, which are of immense importance in understanding molecular structure and reactivity patterns.     

FPV: fast protein visualization using Java 3D

MOTIVATION: Many tools have been developed to visualize protein structures. Tools that have been based on Java 3D((TM)) are compatible among different systems and they can be run remotely through web browsers. However, using Java 3D for visualization has some performance issues with it. The primary concerns about molecular visualization tools based on Java 3D are in their being slow in terms of interaction speed and in their inability to load large molecules. This behavior is especially apparent when the number of atoms to be displayed is huge, or when several proteins are to be displayed simultaneously for comparison. RESULTS: In this paper we present techniques for organizing a Java 3D scene graph to tackle these problems. We have developed a protein visualization system based on Java 3D and these techniques. We demonstrate the effectiveness of the proposed method by comparing the visualization component of our system with two other Java 3D based molecular visualization tools. In particular, for van der Waals display mode, with the efficient organization of the scene graph, we could achieve up to eight times improvement in rendering speed and could load molecules three times as large as the previous systems could. AVAILABILITY: EPV is freely available with source code at the following URL: http://www.cs.ucsb.edu/~tcan/fpv/     

Three-dimensional surface reconstruction software system for IBM personal computers

Surface and volumetric three-dimensional imaging methods have found application in fields as diverse as diagnostic medical imaging and paleontological research. The acquisition, modeling, classification, and computer graphics rendering of discrete image volumes will be introduced. Applications in diagnosis (craniofacial, orthopedic, cardiovascular, and others) as well as reconstruction methods for generic serial sections will be described. C language software for three-dimensional reconstruction which operates on an IBM PC/AT clone is described.