Przemyslaw Chelminiak

Output-input ratio in thermally fluctuating biomolecular machines.

Biological molecular machines are proteins that operate under isothermal conditions and hence are referred to as free energy transducers. They can be formally considered as enzymes that simultaneously catalyze two chemical reactions: the free energy-donating (input) reaction and the free energy-accepting (output) one. Most if not all biologically active proteins display a slow stochastic dynamics of transitions between a variety of conformational substates composing their native state. This makes the description of the enzymatic reaction kinetics in terms of conventional rate constants insufficient. In the steady state, upon taking advantage of the assumption that each reaction proceeds through a single pair (the gate) of transition conformational substates of the enzyme-substrates complex, the degree of coupling between the output and the input reaction fluxes has been expressed in terms of the mean first-passage times on a conformational transition network between the distinguished substates. The theory is confronted with the results of random-walk simulations on the five-dimensional hypercube. The formal proof is given that, for single input and output gates, the output-input degree of coupling cannot exceed unity. As some experiments suggest such exceeding, looking for the conditions for increasing the degree of coupling value over unity challenges the theory. Performed simulations of random walks on several model networks involving more extended gates indicate that the case of the degree of coupling value higher than 1 is realized in a natural way on critical branching trees extended by long-range shortcuts. Such networks are scale-free and display the property of the small world. For short-range shortcuts, the networks are scale-free and fractal, representing a reasonable model for biomolecular machines displaying tight coupling, i.e., the degree of coupling equal exactly to unity. A hypothesis is stated that the protein conformational transition networks, as just as higher-level biological networks, the protein interaction network, and the metabolic network, have evolved in the process of self-organized criticality.

Stochastic Dynamics of Proteins and the Action of Biological Molecular Machines

It is now well established that most if not all enzymatic proteins display a slow stochastic dynamics of transitions between a variety of conformational substates composing their native state. A hypothesis is stated that the protein conformational transition networks, as just as higher-level biological networks, the protein interaction network, and the metabolic network, have evolved in the process of self-organized criticality. Here, the criticality means that all the three classes of networks are scale-free and, moreover, display a transition from the fractal organization on a small length-scale to the small-world organization on the large length-scale. Good mathematical models of such networks are stochastic critical branching trees extended by long-range shortcuts. Biological molecular machines are proteins that operate under isothermal conditions and hence are referred to as free energy transducers. They can be formally considered as enzymes that simultaneously catalyze two chemical reactions: the free energy-donating (input) reaction and the free energy-accepting (output) one. The far-from-equilibrium degree of coupling between the output and the input reaction fluxes have been studied both theoretically and by means of the Monte Carlo simulations on model networks. For single input and output gates the degree of coupling cannot exceed unity. Study simulations of random walks on model networks involving more extended gates indicate that the case of the degree of coupling value higher than one is realized on the mentioned above critical branching trees extended by long-range shortcuts.

Temperature and Detection-Wavelength Dependence of the Electron Transfer Rates in Initial Stages of Photosynthesis

Unusual temperature behavior, observed in the initial electron transfer stages in the photosynthetic reaction centers of the purple bacteria, and a strong probing pulse wavelength dependence of transfer rates, determined in transient absorption spectroscopy, can easily be explained on assuming that the transfer takes place from dynamically unrelaxed states of protein environment. The transitions from the primary special pair (P) to a single bacteriochlorophyll (B) and next to a bacteriopheophytin (H) are controlled by diffusion down the energy value of underdamped vibrational modes of frequency 200 K, probably determining distances between the succeeding cofactors. The subsequent transition to the quinone A (Q) is controlled by diffusion in the position value of an overdamped conformational mode, probably corresponding to the local polarization. From the fit of available experimental data to simple theoretical formulas, the important physical conclusion arises that the very electronic transitions are fast as compared to the relaxation processes and, in the first approximation, only the latter contribute to the overall times of the initial electron transfer stages in photosynthesis.

RaTrav: a tool for calculating mean first-passage times on biochemical networks

BackgroundThe concept of mean first-passage times (MFPTs) occupies an important place in the theory of stochastic processes, with the methods of their calculation being equally important in theoretical physics, chemistry and biology. We present here a software tool designed to support computational biology studies where Markovian dynamics takes place and MFPTs between initial and single or multiple final states in network-like systems are used. Two methods are made available for which their efficiency is strongly dependent on the topology of the defined network: the combinatorial Hill technique and the Monte Carlo simulation method.ResultsAfter a brief introduction to RaTrav, we highlight the utility of MFPT calculations by providing two examples (accompanied by Additional file 1) where they are deemed to be of importance: analysis of a protein-protein docking funnel and interpretation of the free energy transduction between two coupled enzymatic reactions controlled by the dynamics of transition between enzyme conformational states.ConclusionsRaTrav is a versatile and easy to use software tool for calculating MFPTs across biochemical networks. The user simply prepares a text file with the structure of a given network, along with some additional basic parameters such as transition probabilities, waiting probabilities (if any) and local times (weights of edges), which define explicitly the stochastic dynamics on the network. The RaTrav tool can then be applied in order to compute desired MFPTs. For the provided examples, we were able to find the favourable binding path within a protein-protein docking funnel and to calculate the degree of coupling for two chemical reactions catalysed simultaneously by the same protein enzyme. However, the list of possible applications is much wider.