Analysis and Visualization of Molecular Structures
The analysis of complex molecular structures plays an important role in biochemistry, biophysics, molecular biology and pharmaceutics. Therefore, many visualization techniques have been developed for depicting and interactively analyzing molecules. In this project we address two related problems in molecular analysis:
(i) characterization and visualization of steric hindrance or spatial accessibility of molecules and
(ii) computation and depiction of long molecular paths.
Steric effects in molecules arise from the fact that each atom occupies a certain amount of space. Bringing atoms too close together, requires a high amount of energy due to the Pauli repulsion between overlapping electron clouds. The spatial accessibility of molecules (approximated by sets of spherical atoms) is usually visualized by molecular surfaces. We develop improved and accelerated algorithms for computing such molecular surfaces. Furthermore, we propose a new type of molecular surface that better characterizes the mutual spatial accessibility of molecules.
In molecular dynamics (MD) computations it is difficult to bridge the gap between the fs time scale of the oscillations of covalent bonds and the µs timescale of the interesting biological processes. For instance, often it is computationally too expensive to compute tracks of substrates to binding sites of proteins. We address this problem by computing sets of geometrically correct molecular paths in the presence of steric hindrances resulting from the vanderWaals radii of the atoms. These sets of geometrically possible paths constrain the subsets of dynamically possible physical paths. With specialized path visualization techniques, we achieve fast and perceptually optimized interactive representations.
Furthermore, we develop visual analytics tools for uncovering possible molecular paths in complex, dynamic protein configurations. This allows us to analyze the paths of small molecules through transport proteins in cell membranes that become possible due to conformational changes.
Introduction
The analysis and visualization of molecular structures is a large research area. In our project, we focus on methods based on the vanderWaals spheres of the atoms. We developed fast computation and rendering methods for molecular surfaces. These surfaces are most often used to study molecular behaviour and interactions. They show regions accessible for other molecules like substrates. Furthermore, we work on the extraction of geometric molecular paths. We developed a fast approach to compute all of these paths inside a molecule. Our novel visualization techniques support the analysis and offer easy views into cavities and tunnels.
Surface Definitions
Two very important molecular surfaces are the solvent excluded surface (SES) and the molecular skin surface (MSS). Both surfaces are based on the vanderWaals spheres of the atoms. The SES is defined as the surface enclosing all points, that are not accessible to a given probe sphere. This probe typically approximates solvents, most often water is used. The SES can be seen as the track of the probe which rolls over the atom spheres. In contrast, the MSS has no chemical background and is defined by a shrink factor which decomposes the surface into quadric patches. If the shrink factor is close to 0, the surface becomes nearly the convex hull of the atom spheres. For a shrink factor equal to 1, the surface is identical to the surface described by the atom spheres, called vanderWaals surface. The MSS is complete tangent continuous for a shrink factor smaller than 1. In contrast, the SES can have singularities represented by sharp edges and consists of spherical patches as well as toroidal patches of degree four.
Surface Computation
The analytical description of the SES as well as the MSS can be computed very fast by weighted Voronoi diagrams given by the position of the atoms and their vanderWaals radii. It is not necessary to compute the correct Voronoi diagram, but an approximation of the Voronoi cells is sufficient, where each cell depends only on the neighboring atoms. Varshney et al. (1994) presented a definition for the neighborhood of atoms for the SES and a fast parallel algorithm. We developed the neighborhood definition for the MSS to use the same concept.
Surface Rendering
Most visualization tools triangulate the molecular surfaces. This takes several seconds for molecules with a few thousand atoms. For large molecules, the surfaces can consist of several million triangles, so it is not possible to visualize them interactively. We developed a GPUbased raycasting technique for both surfaces, similar to Chavent et al. (2009) and Krone et al (2010). This technique allows the user to visualize surfaces with about 100.000 atoms at interactive frame rates. Combined with our fast computation we are able to render dynamic surfaces for molecules with several thousand atoms. We implemented also classical visualization techniques like coloring, silhouettes, depth darkening and surface blending.
Path Computation
Currently, we compute molecular paths based on the geometry of the vanderWaals spheres of the atoms. The edges of the Voronoi diagram of the atom spheres consists of all possible paths with maximal circular cross section. We developed data structures for a very fast computation of the Voronoi diagram of spheres with the edge tracing algorithm by Kim et al. (2005). A lot of paths of the Voronoi diagram are redundant, because there are more significant paths. Furthermore, there are many paths outside of the domain of the molecule and paths with too small or negative cross sections. We developed a novel filtering technique to keep only the most significant paths of the molecule. The overall computation takes less than one minute, even for molecules with approximately 100.000 atoms.
Path Visualization
Path visualization in dense geometry is a challenging problem in computer graphics. We developed several techniques to enable a fast analysis of the molecular structure. In a first step, we place a lot of small point lights along the extracted paths, which highlight the molecular surface in cavities and tunnels. To keep interactive frame rates, we use deferred shading with high dynamic range rendering to avoid artifacts. We also implemented screen space ambient occlusion for a better depth impression. Furthermore, we provide surface clipping for selected paths, which allows one to analyze the dimension and the cavities of the selected paths.
Dynamic Paths
In order to analyze cavities and channels in dynamic molecular data from simulations, we compute all molecular paths for all time steps of the trajectory. The user is allowed to select one or more path components of the overall path graph for an arbitrary time step. Note that each path component corresponds to exactly one cavity or channel. The selected path components will be traced over time by analyzing the intersections of the corresponding cavities with all cavities in the consecutive time step. A valid tracing is found if an intersection is large enough so that a given probe sphere is able to move from one cavity into a cavity of the consecutive time step. During this analysis, cavity splits and merges are detected and visualized in a 2dimensional graph. Each path in this graph describes a possible geometrical dynamic behavior of the cavity. Thus, dynamic molecular channels can be detected. The dynamic cavities and channels can be visualized either by the skin surface or volume rendering. We tested our tool for the membrane protein bacteriorhodopsin, which works as a light driven proton pump.
Publications
2019 

Norbert Lindow, Daniel Baum, Morgan Leborgne, HansChristian Hege  Interactive Visualization of RNA and DNA Structures  IEEE Transactions on Visualization and Computer Graphics, 25(1), pp. 967976, 2019 (preprint available as ZIBReport 1833) 
PDF (ZIBReport)
BibTeX DOI 
2018 

Péter Koltai, Christof Schütte  A multiscale perturbation expansion approach for Markov state modeling of nonstationary molecular dynamics  SIAM J. Multiscale Model. Simul., 16(4), pp. 14551485, 2018 (preprint available as ) 
BibTeX
DOI 
Norbert Lindow, Daniel Baum, HansChristian Hege  Atomic Accessibility Radii for Molecular Dynamics Analysis  Workshop on Molecular Graphics and Visual Analysis of Molecular Data, 2018, ISBN: 9783038680611 (preprint available as ZIBReport 1818) 
PDF (ZIBReport)
BibTeX DOI 
Manuel Dibak, Mauricio J. del Razo, David De Sancho, Christof Schütte, Frank Noé  MSM/RD: Coupling Markov state models of molecular kinetics with reactiondiffusion simulations  Journal of Chemical Physics, 148(214107), 2018 
BibTeX
DOI 
2016 

Norbert Lindow  Visual Analysis of Atomic Structures Based on the HardSphere Model  Doctoral thesis, Freie Universität Berlin, Christof Schütte (Advisor), 2016 
PDF
BibTeX URN 
Michael Krone, Barbora Kozlíková, Norbert Lindow, Marc Baaden, Daniel Baum, Julius Parulek, HansChristian Hege, Ivan Viola  Visual Analysis of Biomolecular Cavities: State of the Art  Computer Graphics Forum, 35(3), pp. 527551, 2016 (preprint available as ZIBReport 1642) 
PDF (ZIBReport)
BibTeX DOI 
Barbora Kozlíková, Michael Krone, Martin Falk, Norbert Lindow, Marc Baaden, Daniel Baum, Ivan Viola, Julius Parulek, HansChristian Hege  Visualization of Biomolecular Structures: State of the Art Revisited  Computer Graphics Forum, 36(8), pp. 178204, 2016 
BibTeX
DOI 
2015 

Zoe Cournia, Toby W. Allen, Ioan Andricioaei, Bruno Antonny, Daniel Baum, Grace Brannigan, NicolaeViorel Buchete, Jason T. Deckman, Lucie Delemotte, Coral del Val, Ran Friedman, Paraskevi Gkeka, HansChristian Hege, Jérôme Hénin, Marina A. Kasimova, Antonios Kolocouris, Michael L. Klein, Syma Khalid, Joanne Lemieux, Norbert Lindow, Mahua Roy, Jana Selent, Mounir Tarek, Florentina Tofoleanu, Stefano Vanni, Sinisa Urban, David J. Wales, Jeremy C. Smith, AnaNicoleta Bondar  Membrane Protein Structure, Function and Dynamics: A Perspective from Experiments and Theory  Journal of Membrane Biology, 248(4), pp. 611640, 2015 
BibTeX
DOI 
Barbora Kozlikova, Michael Krone, Norbert Lindow, Martin Falk, Marc Baaden, Daniel Baum, Ivan Viola, Julius Parulek, HansChristian Hege  Visualization of Biomolecular Structures: State of the Art  EuroVis 2015 STARS Proceedings, pp. 6181, 2015 (preprint available as ZIBReport 1563) 
PDF (ZIBReport)
BibTeX DOI 
2014 

Norbert Lindow, Daniel Baum, HansChristian Hege  Ligand Excluded Surface: A New Type of Molecular Surface  IEEE Transactions on Visualization and Computer Graphics, 20(12), pp. 24862495, 2014 (preprint available as ZIBReport 1427) 
PDF (ZIBReport)
BibTeX DOI 
HansChristian Hege  Visual analysis of molecular dynamics data using geometric and topological methods  Forum "Math for Industry" 2014, pp. 4546, Vol.57, MI Lecture Notes, Kyushu University, 2014 
BibTeX

2013 

Norbert Lindow, Daniel Baum, AnaNicoleta Bondar, HansChristian Hege  Exploring cavity dynamics in biomolecular systems  BMC Bioinformatics, Vol.14, 2013 
BibTeX
DOI 
2012 

Norbert Lindow, Daniel Baum, AnaNicoleta Bondar, HansChristian Hege  Dynamic Channels in Biomolecular Systems: Path Analysis and Visualization  Proceedings of IEEE Symposium on Biological Data Visualization (biovis’12), pp. 99106, 2012 
BibTeX
DOI 
2011 

Norbert Lindow, Daniel Baum, HansChristian Hege  VoronoiBased Extraction and Visualization of Molecular Paths  IEEE Transactions on Visualization and Computer Graphics, 17(12), pp. 20252034, 2011 
BibTeX
DOI 
2010 

Norbert Lindow, Daniel Baum, Steffen Prohaska, HansChristian Hege  Accelerated Visualization of Dynamic Molecular Surfaces  Comput. Graph. Forum, Vol.29, pp. 943952, 2010 
BibTeX
DOI 