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/     

Volume-rendering representations of ab initio electron densities

Many molecular properties are three-dimensional scalar fields, making the representation of their values throughout the molecular volume a four-dimensional problem. Volume rendering offers a means of graphically displaying the value of a property throughout a volume. In this paper we present the result of an interface of the ab initio electronic structure program, GAUSSIAN 85, and the interactive volume-rendering program, VOXEL-VIEW PLUS, and its initial application: the representation of electron densities for acetylene and difluoroacetylene, and transition structures for the corresponding vinylidene isomerizations.     

PyContact: Rapid,Customizable, and Visual Analysis of Noncovalent Interactions in MD Simulations

Molecular dynamics (MD) simulations have become ubiquitous in all areas of life sciences. The size and model complexity of MD simulations are rapidly growing along with increasing computing power and improved algorithms. This growth has led to the production of a large amount of simulation data that need to be filtered for relevant information to address specific biomedical and biochemical questions. One of the most relevant molecular properties that can be investigated by all-atom MD simulations is the time-dependent evolution of the complex noncovalent interaction networks governing such fundamental aspects as molecular recognition, binding strength, and mechanical and structural stability. Extracting, evaluating, and visualizing noncovalent interactions is a key task in the daily work of structural biologists. We have developed PyContact, an easy-to-use, highly flexible, and intuitive graphical user interface-based application, designed to provide a toolkit to investigate biomolecular interactions in MD trajectories. PyContact is designed to facilitate this task by enabling identification of relevant noncovalent interactions in a comprehensible manner. The implementation of PyContact as a standalone application enables rapid analysis and data visualization without any additional programming requirements, and also preserves full in-program customization and extension capabilities for advanced users. The statistical analysis representation is interactively combined with full mapping of the results on the molecular system through the synergistic connection between PyContact and VMD. We showcase the capabilities and scientific significance of PyContact by analyzing and visualizing in great detail the noncovalent interactions underlying the ion permeation pathway of the human P2X₃ receptor. As a second application, we examine the protein-protein interaction network of the mechanically ultrastable cohesin-dockering complex.     

ESSA: an integrated and interactive computer tool for analysing RNA secondary structure.

With ESSA, we propose an approach of RNA secondary structure analysis based on extensive viewing within a friendly graphical interface. This computer program is organized around the display of folding models produced by two complementary methods suitable to draw long RNA molecules. Any feature of interest can be managed directly on the display and highlighted by a rich combination of colours and symbols with emphasis given to structural probe accessibilities. ESSA also includes a word searching procedure allowing easy visual identification of structural features even complex and degenerated. Analysis functions make it possible to calculate the thermodynamic stability of any part of a folding using several models and compare homologous aligned RNA both in primary and secondary structure. The predictive capacities of ESSA which brings together the experimental, thermodynamic and comparative methods, are increased by coupling it with a program dedicated to RNA folding prediction based on constraints management and propagation. The potentialities of ESSA are illustrated by the identification of a possible tertiary motif in the LSU rRNA and the visualization of a pseudoknot in S15 mRNA.     

XmMol: An X11 and motif program for macromolecular visualization and modeling

XmMol is a desktop tool designed to provide both interactive molecular graphics on X11 displays and easy interface with external applications. A kernel provides an interactive wire-frame display of macromolecules. It supports depth cueing, 3D clipping, and stereo. Various representations, coloring, and labeling modes are proposed. Docking and interactive back-bone deformation tools are also supported. Communication protocols allow the user to develop new external features or to use XmMol as a visualization tool for external numerical programs.     

MD display: An interactive graphics program for visualization of molecular dynamics trajectories

MD Display was developed as a means of visualizing molecular dynamic trajectories generated by Amber.¹ The program runs on Silicon Graphics workstations, and features a simple user interface, and convenient display and analysis options. The program has now been extended to accept input from several other molecular dynamics programs.     

Approximation and visualization of large-scale motion of protein surfaces

We present a method for the approximation and real-time visualization of large-scale motion of protein surfaces. A molecular surface is represented by an expansion of spherical harmonic functions, and the motion of protein atoms around their equilibrium positions is computed by normal mode analysis. The motion of the surface is approximated by projecting the normal mode vectors of the solvent-accessible atoms to the spherical harmonic representation of the molecular surface. These surface motion vectors are represented by a separated spherical harmonic expansion. Representing the surface geometry and the surface motion vectors by spherical harmonic expansions allows variable-resolution analysis and real-time display of the large-scale surface motion. This technique has been applied to interactive visualization, interactive surface manipulation, and animation.     

Alignment editing and identification of consensus secondary structures for nucleic acid sequences: interactive use of dot matrix representations.

We present a computer-aided approach for identifying and aligning consensus secondary structure within a set of functionally related oligonucleotide sequences aligned by sequence. The method relies on visualization of secondary structure using a generalization of the dot matrix representation appropriate for consensus sequence data sets. An interactive computer program implementing such a visualization of consensus structure has been developed. The program allows for alignment editing, data and display filtering and various modes of base pair representation, including co-variation. The utility of this approach is demonstrated with four sample data sets derived from in vitro selection experiments and one data set comprising tRNA sequences.     

Interactive program for visualization and modelling of proteins,nucleic acids and small molecules

The graphics package Insight for the DEC VAX and Evans and Sutherland PS300, created as part of a joint university-industry research project, provides a broad set of capabilities which allow the user to display molecular models in stick figure and surface representation. The Insight program allows the user to model and manipulate proteins, nucleic acids and small molecules. The software accepts coordinate input from several possible sources and provides both a command and menu interface for manipulation of the graphics objects. The command language and program structure make it easy for the biochemist or molecular biologist to use.     

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 new program OPAL for molecular dynamics simulations and energy refinements of biological macromolecules

Summary A new program for molecular dynamics (MD) simulation and energy refinement of biological macromolecules, OPAL, is introduced. Combined with the supporting program TRAJEC for the analysis of MD trajectories, OPAL affords high efficiency and flexibility for work with diferent force fields, and offers a user-friendly interface and extensive trajectory analysis capabilities. Salient features are computational speeds of up to 1.5 GFlops on vector supercomputers such as the NEC SX-3, ellipsoidal boundaries to reduce the system size for studies in explicit solvents, and natural treatment of the hydrostatic pressure. Practical applications of OPAL are illustrated with MD simulations of pure water, energy minimization of the NMR structure of the mixed disulfide of a mutant E. coli glutaredoxin with glutathione in different solvent models, and MD simulations of a small protein, pheromone Er-2, using either instantaneous or time-averaged NMR restraints, or no restraints.Abbreviations D diffusion constant in cm²/s - Er-2 pheromone 2 from Euplotes raikovi - GFlop one billion floating point operations per second - Grx(C14S)-SG mixed disulfide between a mutant E. coli glutaredoxin, with Cys¹⁴ replaced by Ser, and glutathione - MD molecular dynamics - NOE nuclear Overhauser enhancement - rmsd root-mean-square deviation - density in g/cm³     

Recent advances in macromolecular hydrodynamic modeling

The modern implementation of the boundary element method [23] has ushered unprecedented accuracy and precision for the solution of the Stokes equations of hydrodynamics with stick boundary conditions. This article begins by reviewing computations with the program BEST of smooth surface objects such as ellipsoids, the dumbbell, and cylinders that demonstrate that the numerical solution of the integral equation formulation of hydrodynamics yields very high precision and accuracy. When BEST is used for macromolecular computations, the limiting factor becomes the definition of the molecular hydrodynamic surface and the implied effective solvation of the molecular surface. Studies on 49 different proteins, ranging in molecular weight from 9 to over 400kDa, have shown that a model using a 1.1? thick hydration layer describes all protein transport properties very well for the overwhelming majority of them. In addition, this data implies that the crystal structure is an excellent representation of the average solution structure for most of them. In order to investigate the origin of a handful of significant discrepancies in some multimeric proteins (about -20% observed in the intrinsic viscosity), the technique of Molecular Dynamics simulation (MD) has been incorporated into the research program. A preliminary study of dimeric α-chymotrypsin using approximate implicit water MD is presented. In addition I describe the successful validation of modern protein force fields, ff03 and ff99SB, for the accurate computation of solution structure in explicit water simulation by comparison of trajectory ensemble average computed transport properties with experimental measurements. This work includes small proteins such as lysozyme, ribonuclease and ubiquitin using trajectories around 10ns duration. We have also studied a 150kDa flexible monoclonal IgG antibody, Trastuzumab, with multiple independent trajectories encompassing over 320ns of simulation. The close agreement within experimental error of the computed and measured properties allows us to conclude that MD does produce structures typical of those in solution, and that flexible molecules can be properly described using the method of ensemble averaging over a trajectory. We review similar work on the study of a transfer RNA molecule and DNA oligomers that demonstrate that within 3% a simple uniform hydration model 1.1? thick provides agreement with experiment for these nucleic acids. In the case of linear oligomers, the precision can be improved close to 1% by a non-uniform hydration model that hydrates mainly in the DNA grooves, in agreement with high resolution X-ray diffraction. We conclude with a vista on planned improvements for the BEST program to decrease its memory requirements and increase its speed without sacrificing accuracy.     

Nucleic acid visualization with UCSF Chimera

With the increase in the number of large, 3D, high-resolution nucleic acid structures, particularly of the 30S and 50S ribosomal subunits and the intact bacterial ribosome, advancements in the visualization of nucleic acid structural features are essential. Large molecular structures are complicated and detailed, and one goal of visualization software is to allow the user to simplify the display of some features and accent others. We describe an extension to the UCSF Chimera molecular visualization system for the purpose of displaying and highlighting nucleic acid characteristics, including a new representation of sugar pucker, several options for abstraction of base geometries that emphasize stacking and base pairing, and an adaptation of the ribbon backbone to accommodate the nucleic acid backbone. Molecules are displayed and manipulated interactively, allowing the user to change the representations as desired for small molecules, proteins and nucleic acids. This software is available as part of the UCSF Chimera molecular visualization system and thus is integrated with a suite of existing tools for molecular graphics.     

MDMovie: A molecular dynamics viewing tool

The graphics program MDMovie (Molecular Dynamics Movie), written in C using IRIS GL graphics library calls, is designed to facilitate the visualization and interpretation of empirical force field data. MDMovie was created and initially adapted in accord with the needs of physical chemists and thereafter became an expandable analysis tool. Capabilities include the display of chemical structure, animation of molecular dynamics and Monte Carlo trajectories, and the visual representation of various vector and scalar dynamical properties. In addition to being a research tool, MDMovie has features for creating presentation videos and hardcopy output. A library is also available for linking to Fortran simulation codes running on a remote machine and connecting to MDMovie via a socket connection. MDMovie continues to be an ongoing research project and new features are actively being added in collaboration with various research groups. Future plans include porting to OpenGL and the design of an X11-based user interface.     

MOLMOL: A program for display and analysis of macromolecular structures

MOLMOL is a molecular graphics program for display, analysis, and manipulation of three-dimensional structures of biological macromolecules, with special emphasis on nuclear magnetic resonance (NMR) solution structures of proteins and nucleic acids. MOLMOL has a graphical user interface with menus, dialog boxes, and on-line help. The display possibilities include conventional presentation, as well as novel schematic drawings, with the option of combining different presentations in one view of a molecule. Covalent molecular structures can be modified by addition or removal of individual atoms and bonds, and three-dimensional structures can be manipulated by interactive rotation about individual bonds. Special efforts were made to allow for appropriate display and analysis of the sets of typically 20–40 conformers that are conventionally used to represent the result of an NMR structure determination, using functions for superimposing sets of conformers, calculation of root mean square distance (RMSD) values, identification of hydrogen bonds, checking and displaying violations of NMR constraints, and identification and listing of short distances between pairs of hydrogen atoms.     

Key structural features of the actin filament Arp2/3 complex branch junction revealed by molecular simulation

We investigated the structure, properties and dynamics of the actin filament branch junction formed by actin-related protein (Arp) 2/3 complex using all-atom molecular dynamics (MD) simulations based on a model fit to a reconstruction from electron tomograms. Simulations of the entire structure consisting of 31 protein subunits together with solvent molecules containing ～3 million atoms were performed for an aggregate time of 175 ns. One 75-ns simulation of the original reconstruction was compared to two 50-ns simulations of alternate structures, showing that the hypothesized branch junction structure is very stable. Our simulations revealed that the interface between Arp2/3 complex and the mother actin filament features a large number of salt bridges and hydrophobic contacts, many of which are dynamic and formed/broken on the timescale of the simulation. The simulations suggest that the DNase binding loops in Arp3, and possibly Arp2, form stabilizing contacts with the mother filament. Unbiased comparison of models sampled from the MD simulation trajectory with the primary experimental electron tomography data identified regions were snapshots from the simulation provide atomic details of the model structures and also pinpoints regions where the initial modeling based on the electron tomogram reconstruction may be suboptimal.