Jan Byška

MoleCollar and Tunnel Heat Map Visualizations for Conveying Spatio-Temporo-Chemical Properties Across and Along Protein Voids

Studying the characteristics of proteins and their inner void space, including their geometry, physico‐chemical properties and dynamics are instrumental for evaluating the reactivity of the protein with other small molecules. The analysis of long simulations of molecular dynamics produces a large number of voids which have to be further explored and evaluated. In this paper we propose three new methods: two of them convey important properties along the long axis of a selected void during molecular dynamics and one provides a comprehensive picture across the void. The first two proposed methods use a specific heat map to present two types of information: an overview of all detected tunnels in the dynamics and their bottleneck width and stability over time, and an overview of a specific tunnel in the dynamics showing the bottleneck position and changes of the tunnel length over time. These methods help to select a small subset of tunnels, which are explored individually and in detail. For this stage we propose the third method, which shows in one static image the temporal evolvement of the shape of the most critical tunnel part, i.e., its bottleneck. This view is enriched with abstract depictions of different physicochemical properties of the amino acids surrounding the bottleneck. The usefulness of our newly proposed methods is demonstrated on a case study and the feedback from the domain experts is included. The biochemists confirmed that our novel methods help to convey the information about the appearance and properties of tunnels in a very intuitive and comprehensible manner.

AnimoAminoMiner: Exploration of Protein Tunnels and their Properties in Molecular Dynamics

In this paper we propose a novel method for the interactive exploration of protein tunnels. The basic principle of our approach is that we entirely abstract from the 3D/4D space the simulated phenomenon is embedded in. A complex 3D structure and its curvature information is represented only by a straightened tunnel centerline and its width profile. This representation focuses on a key aspect of the studied geometry and frees up graphical estate to key chemical and physical properties represented by surrounding amino acids. The method shows the detailed tunnel profile and its temporal aggregation. The profile is interactively linked with a visual overview of all amino acids which are lining the tunnel over time. In this overview, each amino acid is represented by a set of colored lines depicting the spatial and temporal impact of the amino acid on the corresponding tunnel. This representation clearly shows the importance of amino acids with respect to selected criteria. It helps the biochemists to select the candidate amino acids for mutation which changes the protein function in a desired way. The AnimoAminoMiner was designed in close cooperation with domain experts. Its usefulness is documented by their feedback and a case study, which are included.

Comparative visualization of protein secondary structures

BackgroundProtein function is determined by many factors, namely by its constitution, spatial arrangement, and dynamic behavior. Studying these factors helps the biochemists and biologists to better understand the protein behavior and to design proteins with modified properties. One of the most common approaches to these studies is to compare the protein structure with other molecules and to reveal similarities and differences in their polypeptide chains.ResultsWe support the comparison process by proposing a new visualization technique that bridges the gap between traditionally used 1D and 3D representations. By introducing the information about mutual positions of protein chains into the 1D sequential representation the users are able to observe the spatial differences between the proteins without any occlusion commonly present in 3D view. Our representation is designed to serve namely for comparison of multiple proteins or a set of time steps of molecular dynamics simulation.ConclusionsThe novel representation is demonstrated on two usage scenarios. The first scenario aims to compare a set of proteins from the family of cytochromes P450 where the position of the secondary structures has a significant impact on the substrate channeling. The second scenario focuses on the protein flexibility when by comparing a set of time steps our representation helps to reveal the most dynamically changing parts of the protein chain.

Geometry-based Algorithm for Detection of Asymmetric Tunnels in Protein Molecules

We present a novel geometry-based method for computing asymmetric tunnels and voids in proteins, approximating their real shape with selected precision. Our method combines ideas from Voronoi and grid based approaches for protein analysis. We represent tunnels in protein using voxel data grid which allows us to store their shape more accurately. Our algorithm employs a tunnel skeleton computed using Voronoi diagram. The skeleton allows us to perform grid computation in a bounded space, with lower time and memory demands, and easily identify and measure individual tunnels.

Watergate: Visual Exploration of Water Trajectories in Protein Dynamics

The function of proteins is tightly related to their interactions with other molecules. The study of such interactions often requires to track the molecules that enter or exit specific regions of the proteins. This is investigated with molecular dynamics simulations, producing the trajectories of thousands of water molecules during hundreds of thousands of time steps. To ease the exploration of such rich spatio-temporal data, we propose a novel workflow for the analysis and visualization of large sets of water-molecule trajectories. Our solution consists of a set of visualization techniques, which help biochemists to classify, cluster, and filter the trajectories and to explore the properties and behavior of selected subsets in detail. Initially, we use an interactive histogram and a time-line visualization to give an overview of all water trajectories and select the interesting ones for further investigation. Further, we depict clusters of trajectories in a novel 2D representation illustrating the flows of water molecules. These views are interactively linked with a 3D representation where we show individual paths, including their simplification, as well as extracted statistical information displayed by isosurfaces. The proposed solution has been designed in tight collaboration with experts to support specific tasks in their scientific workflows. They also conducted several case studies to evaluate the usability and effectiveness of our new solution with respect to their research scenarios. These confirmed that our proposed solution helps in analyzing water trajectories and in extracting the essential information out of the large amount of input data.

Interactive Exploration of Ligand Transportation through Protein Tunnels

BackgroundProtein structures and their interaction with ligands have been in the focus of biochemistry and structural biology research for decades. The transportation of ligand into the protein active site is often complex process, driven by geometric and physico-chemical properties, which renders the ligand path full of jitter and impasses. This prevents understanding of the ligand transportation and reasoning behind its behavior along the path.ResultsTo address the needs of the domain experts we design an explorative visualization solution based on a multi-scale simplification model. It helps to navigate the user to the most interesting parts of the ligand trajectory by exploring different attributes of the ligand and its movement, such as its distance to the active site, changes of amino acids lining the ligand, or ligand “stuckness”. The process is supported by three linked views – 3D representation of the simplified trajectory, scatterplot matrix, and bar charts with line representation of ligand-lining amino acids.ConclusionsThe usage of our tool is demonstrated on molecular dynamics simulations provided by the domain experts. The tool was tested by the domain experts from protein engineering and the results confirm that it helps to navigate the user to the most interesting parts of the ligand trajectory and to understand the ligand behavior.

Visibility-based approach to surface detection of tunnels in proteins

Structural properties of proteins substantially influence their reactivity. Thus, the presence of pathways serving for transportation of a small molecule to the protein active site is crucial. These pathways, called tunnels, are defined by their surroundings --- tunnel lining amino acids (or residues). In consequence, studying these amino acids and their properties is tightly connected with protein reactivity. The set of tunnel lining amino acids detected for a given tunnel can differ with respect to selected algorithm for their computation. The criteria for evaluating their biochemical relevance are different as well, as they can depend on various physico-chemical properties. In this paper we firstly present a novel approach to the detection of tunnel lining amino acids. This approach is more robust in comparison with the existing methods. It is based on the visibility of atoms of the amino acids from the tunnel surface, often derived from the Voronoi diagram. Moreover, the detected set of amino acids is further utilized for an automatic detection of the asymmetric tunnel surface. The results are compared with already existing approaches and the benefits are discussed.

Visualizing movements of protein tunnels in molecular dynamics simulations

Analysis and visualization of molecules and their structural features help biochemists and biologists to better understand protein behavior. Studying these structures in molecular dynamics simulations enhances this understanding. In this paper we introduce three approaches for animating specific inner pathways composed of an empty space between atoms, called tunnels. These tunnels facilitate the transport of small molecules, water solvent and ions in many proteins. They help researchers understand the structure-function relationships of proteins and the knowledge of tunnel properties improves the design of new inhibitors. Our methods are derived from selected tunnel representations when each stresses some of the important tunnel properties---width, shape, mapping of physico-chemical properties, etc. Our methods provide smooth animation of the movement of tunnels as they change their length and shape throughout the simulation.

A de novo Ser111Thr variant in aquaporin-4 in a patient with intellectual disability, transient signs of brain ischemia, transient cardiac hypertrophy, and progressive gait disturbance

Aquaporin-4, encoded by AQP4, is the major water channel in the central nervous system and plays an important role in the brains water balance, including edema formation and clearance. Using genomic copy-number analysis and trio-exome sequencing, we investigated a male patient with intellectual disability, hearing loss, and progressive gait dysfunction and found a de novo missense change Ser111Thr in AQP4 as the only suspicious finding. Perinatally, signs of brain ischemia were detected in relation to acute collapse 2 h after birth that resolved a few days later. At the age of 3 mo, cardiac hypertrophy was detected that persisted through childhood but was completely resolved by age 16. In theory, this neurodevelopmental disorder with transient cardiomyopathy could be caused by a disturbance of cellular water balance. Ser111 is an extremely conserved residue in the short cytoplasmic loop between AQP4 transmembrane helix 2 and 3, present across all AQP isoforms from plants to mammals, and it does not appear to be a phosphorylation site. We found that the Ser111Thr change does not affect water permeability or protein stability, suggesting another and possibly regulatory role. Although causality remains unproven, this case study draws attention to AQP4 as a candidate gene for a unique developmental disorder and to a specific serine as a residue of possibly great functional importance in many AQPs.

CAVER Analyst 2.0: analysis and visualization of channels and tunnels in protein structures and molecular dynamics trajectories

Motivation: Studying the transport paths of ligands, solvents, or ions in transmembrane proteins and proteins with buried binding sites is fundamental to the understanding of their biological function. A detailed analysis of the structural features influencing the transport paths is also important for engineering proteins for biomedical and biotechnological applications. Results: CAVER Analyst 2.0 is a software tool for quantitative analysis and real‐time visualization of tunnels and channels in static and dynamic structures. This version provides the users with many new functions, including advanced techniques for intuitive visual inspection of the spatiotemporal behavior of tunnels and channels. Novel integrated algorithms allow an efficient analysis and data reduction in large protein structures and molecular dynamic simulations. Availability and implementation: CAVER Analyst 2.0 is a multi‐platform standalone Java‐based application. Binaries and documentation are freely available at www.caver.cz. Supplementary information: Supplementary data are available at Bioinformatics online.