ResiCon: a method for the identification of dynamic domains, hinges and interfacial regions in proteins

Granger causality based on vector time series and quaternion algebra with possible applications to molecular dynamics data analysis

Causal analysis plays a significant role in physics, chemistry, and biology. Dynamics of complex (bio)molecular and nanosystems, from the microscopic to the macroscopic scale, are characterized by time-dependent vectors such as positions, forces, momenta, angular momenta, or torques. Identification and analysis of causal relationships between these time-dependent signals is an important problem in the multidimensional time-series analysis and is of great practical importance in describing the properties of such dynamical systems, and to understanding their functionality. For linear stochastic systems characterized by multidimensional scalar signals, Granger proposed a simple procedure to detect causal relationships, called Granger causality. In this study we extended this formalism to vector signals representing physical vector quantities. For this purpose, we used quaternion algebra, where vector signals are treated as time-dependent quaternions. The developed analytical model is based on the autoregressive formalism. This formalism (Q-MVAR) and its numerical implementation were validated using two simple dynamic models: a rigid body model represented by a benzenelike molecular fragment, interacting with a short-range harmonic potential with a wall, as well as a system of three model atomic balls moving inside a soft spherical surface and interacting with long range electrostatic forces. Although the motivation to these studies was the analysis of classical motions in complex (bio)molecular systems, described with a mechanical model and based on molecular dynamics (MD) simulations, in particular coarse-grained ones, it should be noted that the developed extended formalism can be applied to any system composed of many rigid elements that interact with arbitrary potentials and are characterized by complex internal motions. A description of the detailed procedure for calculating causality measures is provided in the Appendices of the Supplemental Material. This formalism and the prototype of its numerical implementation can be further developed and applied in many different fields of physical, natural, and engineering sciences.

ART-RRT: As-Rigid-As-Possible search for protein conformational transition paths

The possible functions of a protein are strongly related to its structural rearrangements in the presence of other molecules or environmental changes. Hence, the evaluation of transition paths of proteins, which encodes conformational changes between stable states, is important since it may reveal the underlying mechanisms of the biochemical processes related to these motions. During the last few decades, different geometry-based methods have been proposed to predict such transition paths. However, in the cases where the solution requires complex motions, these methods, which typically constrain only locally the molecular structures, could produce physically irrelevant solutions involving self-intersection. Recently, we have proposed ART-RRT, an efficient method for finding ligand-unbinding pathways. It relies on the exploration of energy valleys in low-dimensional spaces, taking advantage of some mechanisms inspired from computer graphics to ensure the consistency of molecular structures. This article extends ART-RRT to the problem of finding probable conformational transition between two stable states for proteins. It relies on a bidirectional exploration rooted on the two end states and introduces an original strategy to attempt connections between the explored regions. The resulting method is able to produce at low computational cost biologically realistic paths free from self-intersection. These paths can serve as valuable input to other advanced methods for the study of proteins. A better understanding of conformational changes of proteins is important since it may reveal the underlying mechanisms of the biochemical processes related to such motions. Recently, the ART-RRT method has been introduced for finding ligand-unbinding pathways. This article presents an adaptation of the method for finding probable conformational transition between two stable states of a protein. The method is not only computationally cost-effective but also able to produce biologically realistic paths which are free from self-intersection.

Quantitative Characterization of Domain Motions in Molecular Machines

The work of molecular machines such as the ribosome is accompanied by conformational changes, often characterized by relative motions of their domains. The method we have developed seeks to quantify these motions in a general way, facilitating comparisons of results obtained by different researchers. Typically there are multiple snapshots of a structure in the form of pdb coordinates resulting from flexible fitting of low-resolution density maps, from X-ray crystallography, or from molecular dynamics simulation trajectories. Our objective is to characterize the motion of each domain as a coordinate transformation using moments of inertia tensor, a method we developed earlier. What has been missing until now are ancillary tools that make this task practical, general, and biologically informative. We have provided a comprehensive solution to this task with a set of tools implemented on the VMD platform. These tools address the need for reproducible segmentation of domains, and provide a generalized description of their motions using principal axes of inertia. Although this methodology has been specifically developed for studying ribosome motion, it is applicable to any molecular machine.

Visualizing the nanoscale: protein internal dynamics and neutron spin echo spectroscopy

The most complex molecular machines are proteins found within cells. Protein dynamics, in particular dynamics on nanoscales, presents us with a novel paradigm for cell signaling: the idea that proteins and protein complexes can communicate directly within themselves to effect long-range information transfer, via coupled domains and correlated residue clusters. This idea has been little explored, in large part because of a paucity of experimental techniques that can address the necessary questions. Here we review recent progress in developing a promising new approach, neutron spin echo spectroscopy.

WeBIAS: a web server for publishing bioinformatics applications

Background

One of the requirements for a successful scientific tool is its availability. Developing a functional web service, however, is usually considered a mundane and ungratifying task, and quite often neglected. When publishing bioinformatic applications, such attitude puts additional burden on the reviewers who have to cope with poorly designed interfaces in order to assess quality of presented methods, as well as impairs actual usefulness to the scientific community at large.

Results

In this note we present WeBIAS—a simple, self-contained solution to make command-line programs accessible through web forms. It comprises a web portal capable of serving several applications and backend schedulers which carry out computations. The server handles user registration and authentication, stores queries and results, and provides a convenient administrator interface. WeBIAS is implemented in Python and available under GNU Affero General Public License. It has been developed and tested on GNU/Linux compatible platforms covering a vast majority of operational WWW servers. Since it is written in pure Python, it should be easy to deploy also on all other platforms supporting Python (e.g. Windows, Mac OS X). Documentation and source code, as well as a demonstration site are available at http://bioinfo.imdik.pan.pl/webias.

Conclusions

WeBIAS has been designed specifically with ease of installation and deployment of services in mind. Setting up a simple application requires minimal effort, yet it is possible to create visually appealing, feature-rich interfaces for query submission and presentation of results.