Interactive fitting augmented by force-feedback and virtual reality

Interactive Flexible-Receptor Molecular Docking in Virtual Reality Using DockIT

Interactive docking enables the user to guide and control the docking of two biomolecules into a binding pose. It is of particular use when the binding site is known and is thought to be applicable to structure-based drug design (SBDD) and educating students about biomolecular interactions. For SBDD, it enables expertise and intuition to be brought to bear in the drug design process. In education, it can teach students about the most basic level of biomolecular function. Here, we introduce DockIT for virtual reality (VR) that uses a VR headset and hand-held controllers. Using the method of linear response on explicit solvent molecular dynamics simulations, DockIT can model both global and local conformational changes within the receptor due to forces of interaction with the ligand. It has real-time flexible molecular surface rendering and can show the real-time formation and breaking of hydrogen bonds, both between the ligand and receptor and within the receptor itself as it smoothly changes conformation.

Beyond the Black Box: Interactive Global Docking of Protein Complexes

A key to understanding the function of biological systems is the visualization of their natural state, ideally in a natural environment. At a molecular level, this is challenging. Traditional experimental techniques, like X-ray crystallography, can provide the atomic structure of proteins, but only by removing them from their native surroundings and forcing them into crystals. Over the past decade, microscopy techniques have emerged as alternatives to these traditional structure determination methods, with the advantage of visualizing molecules in a near-native state. Given the current focus of structural biology on interactions between proteins and better understanding of large protein complexes, cryo-electron microscopy (cryo-EM) has become a valuable tool.

Multi-scale Visualization of Molecular Architecture Using Real-Time Ambient Occlusion in Sculptor

The modeling of large biomolecular assemblies relies on an efficient rendering of their hierarchical architecture across a wide range of spatial level of detail. We describe a paradigm shift currently under way in computer graphics towards the use of more realistic global illumination models, and we apply the so-called ambient occlusion approach to our open-source multi-scale modeling program, Sculptor. While there are many other higher quality global illumination approaches going all the way up to full GPU-accelerated ray tracing, they do not provide size-specificity of the features they shade. Ambient occlusion is an aspect of global lighting that offers great visual benefits and powerful user customization. By estimating how other molecular shape features affect the reception of light at some surface point, it effectively simulates indirect shadowing. This effect occurs between molecular surfaces that are close to each other, or in pockets such as protein or ligand binding sites. By adding ambient occlusion, large macromolecular systems look much more natural, and the perception of characteristic surface features is strongly enhanced. In this work, we present a real-time implementation of screen space ambient occlusion that delivers realistic cues about tunable spatial scale characteristics of macromolecular architecture. Heretofore, the visualization of large biomolecular systems, comprising e.g. hundreds of thousands of atoms or Mega-Dalton size electron microscopy maps, did not take into account the length scales of interest or the spatial resolution of the data. Our approach has been uniquely customized with shading that is tuned for pockets and cavities of a user-defined size, making it useful for visualizing molecular features at multiple scales of interest. This is a feature that none of the conventional ambient occlusion approaches provide. Actual Sculptor screen shots illustrate how our implementation supports the size-dependent rendering of molecular surface features.

In this work, we present an implementation of screen space ambient occlusion (SSAO), which supports the visualization and modeling of multi-scale biophysical data, such as atomic structures and 3D density maps, at multiple scales of interest. The ever-growing size of macromolecular assemblies presents a formidable challenge to molecular modeling programs. Ambient occlusion (AO) has recently received a lot of attention in high-quality rendering engines as well as in video games. The technique provides an acceptable real-time approximation of global illumination, by adding realistic cues about the indirect shading of the surface topology. SSAO leverages the computational power and flexibility of graphical processors, which are present on almost all modern video cards. By using SSAO, we put an emphasis on compatibility and speed, while delivering the visual benefits of AO at a user-selected level of detail. Our work is particularly timely because many molecular graphics packages are currently adopting global illumination schemes, but our work is unique in providing size-specificity of the shaded features.

Automated tracing of filaments in 3D electron tomography reconstructions using Sculptor and Situs

The molecular graphics program Sculptor and the command-line suite Situs are software packages for the integration of biophysical data across spatial resolution scales. Herein, we provide an overview of recently developed tools relevant to cryo-electron tomography (cryo-ET), with an emphasis on functionality supported by Situs 2.7.1 and Sculptor 2.1.1. We describe a work flow for automatically segmenting filaments in cryo-ET maps including denoising, local normalization, feature detection, and tracing. Tomograms of cellular actin networks exhibit both cross-linked and bundled filament densities. Such filamentous regions in cryo-ET data sets can then be segmented using a stochastic template-based search, VolTrac. The approach combines a genetic algorithm and a bidirectional expansion with a tabu search strategy to localize and characterize filamentous regions. The automated filament segmentation by VolTrac compares well to a manual one performed by expert users, and it allows an efficient and reproducible analysis of large data sets. The software is free, open source, and can be used on Linux, Macintosh or Windows computers.

Using Sculptor and Situs for simultaneous assembly of atomic components into low-resolution shapes

We describe an integrated software system called Sculptor that combines visualization capabilities with molecular modeling algorithms for the analysis of multi-scale data sets. Sculptor features extensive special purpose visualization techniques that are based on modern GPU programming and are capable of representing complex molecular assemblies in real-time. The integration of graphics and modeling offers several advantages. The user interface not only eases the usually steep learning curve of pure algorithmic techniques, but it also permits instant analysis and post-processing of results, as well as the integration of results from external software. Here, we implemented an interactive peak-selection strategy that enables the user to explore a preliminary score landscape generated by the colors tool of Situs. The interactive placement of components, one at a time, is advantageous for low-resolution or ambiguously shaped maps, which are sometimes difficult to interpret by the fully automatic peak selection of colors. For the subsequent refinement of the preliminary models resulting from both interactive and automatic peak selection, we have implemented a novel simultaneous multi-body docking in Sculptor and Situs that softly enforces shape complementarities between components using the normalization of the cross-correlation coefficient. The proposed techniques are freely available in Situs version 2.6 and Sculptor version 2.0.

Evolutionary tabu search strategies for the simultaneous registration of multiple atomic structures in cryo-EM reconstructions

A structural characterization of multi-component cellular assemblies is essential to explain the mechanisms governing biological function. Macromolecular architectures may be revealed by integrating information collected from various biophysical sources - for instance, by interpreting low-resolution electron cryomicroscopy reconstructions in relation to the crystal structures of the constituent fragments. A simultaneous registration of multiple components is beneficial when building atomic models as it introduces additional spatial constraints to facilitate the native placement inside the map. The high-dimensional nature of such a search problem prevents the exhaustive exploration of all possible solutions. Here we introduce a novel method based on genetic algorithms, for the efficient exploration of the multi-body registration search space. The classic scheme of a genetic algorithm was enhanced with new genetic operations, tabu search and parallel computing strategies and validated on a benchmark of synthetic and experimental cryo-EM datasets. Even at a low level of detail, for example 35-40 A, the technique successfully registered multiple component biomolecules, measuring accuracies within one order of magnitude of the nominal resolutions of the maps. The algorithm was implemented using the Sculptor molecular modeling framework, which also provides a user-friendly graphical interface and enables an instantaneous, visual exploration of intermediate solutions.

Problems in fitting high resolution structures into electron microscopic reconstructions

Great advances have been made in electron microscopy (EM) over the past decade, with the result that a number of protein complexes have been solved at near-atomic resolution using EM imaging. However, only a limited number of such complexes are expected to have the high degree of internal order needed to achieve this type of resolution. Many other complexes and polymers will be visualized and reconstructed by EM at an intermediate level of resolution, where the polypeptide chain cannot be directly traced. Crystal and nuclear magnetic resonance structures for components or subunits of these higher-order assemblies are frequently available. One of the greatest strengths of EM continues to be the ability to dock high-resolution structures of components into low or intermediate resolution reconstructions of assemblies to build pseudoatomic models for quaternary structure. This review discusses the strengths and limitations of this approach, with particular emphasis on protein polymers. I discuss how limitations in resolution can lead to ambiguities in building models, and these cannot be always be resolved with available data. The use of homology models for quaternary structure are particularly problematic, given accumulating evidence for the divergence of quaternary structures at the same time that tertiary structure can be conserved.

Visualization software for molecular assemblies

Software for viewing three-dimensional models and maps of viruses, ribosomes, filaments, and other molecular assemblies is advancing on many fronts. New developments include molecular representations that offer better control over level of detail, lighting that improves the perception of depth, and two-dimensional projections that simplify data interpretation. Programmable graphics processors offer quality, speed, and visual effects not previously possible, while 3D printers, haptic interaction devices, and auto-stereo displays show promise in more naturally engaging our senses. Visualization methods are developed by diverse groups of researchers with differing goals: experimental biologists, database developers, computer scientists, and package developers. We survey recent developments and problems faced by the developer community in bringing innovative visualization methods into widespread use.

Multi-resolution anchor-point registration of biomolecular assemblies and their components

An atomic scale interpretation facilitates the assignment of functional properties to 3D reconstructions of macromolecular assemblies in electron microscopy (EM). Such a high-resolution interpretation is typically achieved by docking the known atomic structures of components into the volumetric EM maps. Docking locations are often determined by maximizing the cross-correlation coefficient of the two objects in a slow, exhaustive search. If time is of essence, such as in related visualization and image processing fields, the matching of data is accelerated by incorporating feature points that form a compact description of 3D objects. The complexity reduction afforded by the feature point representation enables a near-instantaneous matching. We show that such reduced matching can also deliver robust and accurate results in the presence of noise or artifacts. We therefore propose a novel multi-resolution registration technique employing feature-based shape descriptions of the volumetric and structural data. The pattern-matching algorithm carries out a hierarchical alignment of the point sets generated by vector quantization. The search-space complexity is reduced by an integrated tree-pruning technique, which permits the detection of subunits in large macromolecular assemblies in real-time. The efficiency and accuracy of the novel algorithm are validated on a standard test system of homo-oligomeric assemblies.