Multiplicative random walk Metropolis-Hastings on the real line

In this article we propose multiplication based random walk Metropolis Hastings (MH) algorithm on the real line. We call it the random dive MH (RDMH) algorithm. This algorithm, even if simple to apply, was not studied earlier in Markov chain Monte Carlo literature. The associated kernel is shown to have standard properties like irreducibility, aperiodicity and Harris recurrence under some mild assumptions. These ensure basic convergence (ergodicity) of the kernel. Further the kernel is shown to be geometric ergodic for a large class of target densities on $\mathbb{R}$. This class even contains realistic target densities for which random walk or Langevin MH are not geometrically ergodic. Three simulation studies are given to demonstrate the mixing property and superiority of RDMH to standard MH algorithms on real line. A share-price return data is also analyzed and the results are compared with those available in the literature

The Goldbeter-Koshland switch in the first-order region and its response to dynamic disorder

In their classical work (Proc. Natl. Acad. Sci. USA, 1981, 78:6840-6844), Goldbeter and Koshland mathematically analyzed a reversible covalent modification system which is highly sensitive to the concentration of effectors. Its signal-response curve appears sigmoidal, constituting a biochemical switch. However, the switch behavior only emerges in the "zero-order region", i.e. when the signal molecule concentration is much lower than that of the substrate it modifies. In this work we showed that the switching behavior can also occur under comparable concentrations of signals and substrates, provided that the signal molecules catalyze the modification reaction in cooperation. We also studied the effect of dynamic disorders on the proposed biochemical switch, in which the enzymatic reaction rates, instead of constant, appear as stochastic functions of time. We showed that the system is robust to dynamic disorder at bulk concentration. But if the dynamic disorder is quasi-static, large fluctuations of the switch response behavior may be observed at low concentrations. Such fluctuation is relevant to many biological functions. It can be reduced by either increasing the conformation interconversion rate of the protein, or correlating the enzymatic reaction rates in the network.Comment: 23 pages, 4 figures, accepted by PLOS ON

Dynamic Disorder in Quasi-Equilibrium Enzymatic Systems

Conformations and catalytic rates of enzymes fluctuate over a wide range of timescales. Despite these fluctuations, there exist some limiting cases in which the enzymatic catalytic rate follows the macroscopic rate equation such as the Michaelis-Menten law. In this paper we investigate the applicability of macroscopic rate laws for fluctuating enzyme systems in which catalytic transitions are slower than ligand binding-dissociation reactions. In this quasi-equilibrium limit, for an arbitrary reaction scheme we show that the catalytic rate has the same dependence on ligand concentrations as obtained from mass-action kinetics even in the presence of slow conformational fluctuations. These results indicate that the timescale of conformational dynamics – no matter how slow – will not affect the enzymatic rate in quasi-equilibrium limit. Our numerical results for two enzyme-catalyzed reaction schemes involving multiple substrates and inhibitors further support our general theory

Selective Targeting of Tumorigenic Cancer Cell Lines by Microtubule Inhibitors

For anticancer drug therapy, it is critical to kill those cells with highest tumorigenic potential, even when they comprise a relatively small fraction of the overall tumor cell population. We have used the established NCI/DTP 60 cell line growth inhibition assay as a platform for exploring the relationship between chemical structure and growth inhibition in both tumorigenic and non-tumorigenic cancer cell lines. Using experimental measurements of “take rate” in ectopic implants as a proxy for tumorigenic potential, we identified eight chemical agents that appear to strongly and selectively inhibit the growth of the most tumorigenic cell lines. Biochemical assay data and structure-activity relationships indicate that these compounds act by inhibiting tubulin polymerization. Yet, their activity against tumorigenic cell lines is more selective than that of the other microtubule inhibitors in clinical use. Biochemical differences in the tubulin subunits that make up microtubules, or differences in the function of microtubules in mitotic spindle assembly or cell division may be associated with the selectivity of these compounds

Reconstruction of bacterial transcription-coupled repair at single-molecule resolution

International audienc

Investigating the robustness of the classical enzyme kinetic equations in small intracellular compartments

<p>Abstract</p> <p>Background</p> <p>Classical descriptions of enzyme kinetics ignore the physical nature of the intracellular environment. Main implicit assumptions behind such approaches are that reactions occur in compartment volumes which are large enough so that molecular discreteness can be ignored and that molecular transport occurs via diffusion. Though these conditions are frequently met in laboratory conditions, they are not characteristic of the intracellular environment, which is compartmentalized at the micron and submicron scales and in which active means of transport play a significant role.</p> <p>Results</p> <p>Starting from a master equation description of enzyme reaction kinetics and assuming metabolic steady-state conditions, we derive novel mesoscopic rate equations which take into account (i) the intrinsic molecular noise due to the low copy number of molecules in intracellular compartments (ii) the physical nature of the substrate transport process, i.e. diffusion or vesicle-mediated transport. These equations replace the conventional macroscopic and deterministic equations in the context of intracellular kinetics. The latter are recovered in the limit of infinite compartment volumes. We find that deviations from the predictions of classical kinetics are pronounced (hundreds of percent in the estimate for the reaction velocity) for enzyme reactions occurring in compartments which are smaller than approximately 200 nm, for the case of substrate transport to the compartment being mediated principally by vesicle or granule transport and in the presence of competitive enzyme inhibitors.</p> <p>Conclusion</p> <p>The derived mesoscopic rate equations describe subcellular enzyme reaction kinetics, taking into account, for the first time, the simultaneous influence of both intrinsic noise and the mode of transport. They clearly show the range of applicability of the conventional deterministic equation models, namely intracellular conditions compatible with diffusive transport and simple enzyme mechanisms in several hundred nanometre-sized compartments. An active transport mechanism coupled with large intrinsic noise in enzyme concentrations is shown to lead to huge deviations from the predictions of deterministic models. This has implications for the common approach of modeling large intracellular reaction networks using ordinary differential equations and also for the calculation of the effective dosage of competitive inhibitor drugs.</p

How to find simple and accurate rules for viral protease cleavage specificities

<p>Abstract</p> <p>Background</p> <p>Proteases of human pathogens are becoming increasingly important drug targets, hence it is necessary to understand their substrate specificity and to interpret this knowledge in practically useful ways. New methods are being developed that produce large amounts of cleavage information for individual proteases and some have been applied to extract cleavage rules from data. However, the hitherto proposed methods for extracting rules have been neither easy to understand nor very accurate. To be practically useful, cleavage rules should be accurate, compact, and expressed in an easily understandable way.</p> <p>Results</p> <p>A new method is presented for producing cleavage rules for viral proteases with seemingly complex cleavage profiles. The method is based on orthogonal search-based rule extraction (OSRE) combined with spectral clustering. It is demonstrated on substrate data sets for human immunodeficiency virus type 1 (HIV-1) protease and hepatitis C (HCV) NS3/4A protease, showing excellent prediction performance for both HIV-1 cleavage and HCV NS3/4A cleavage, agreeing with observed HCV genotype differences. New cleavage rules (consensus sequences) are suggested for HIV-1 and HCV NS3/4A cleavages. The practical usability of the method is also demonstrated by using it to predict the location of an internal cleavage site in the HCV NS3 protease and to correct the location of a previously reported internal cleavage site in the HCV NS3 protease. The method is fast to converge and yields accurate rules, on par with previous results for HIV-1 protease and better than previous state-of-the-art for HCV NS3/4A protease. Moreover, the rules are fewer and simpler than previously obtained with rule extraction methods.</p> <p>Conclusion</p> <p>A rule extraction methodology by searching for multivariate low-order predicates yields results that significantly outperform existing rule bases on out-of-sample data, but are more transparent to expert users. The approach yields rules that are easy to use and useful for interpreting experimental data.</p

Visualizing chemical states and defects induced magnetism of graphene oxide by spatially-resolved-X-ray microscopy and spectroscopy

[[abstract]]This investigation studies the various magnetic behaviors of graphene oxide (GO) and reduced graphene oxides (rGOs) and elucidates the relationship between the chemical states that involve defects therein and their magnetic behaviors in GO sheets. Magnetic hysteresis loop reveals that the GO is ferromagnetic whereas photo-thermal moderately reduced graphene oxide (M-rGO) and heavily reduced graphene oxide (H-rGO) gradually become paramagnetic behavior at room temperature. Scanning transmission X-ray microscopy and corresponding X-ray absorption near-edge structure spectroscopy were utilized to investigate thoroughly the variation of the C 2p(π*) states that are bound with oxygen-containing and hydroxyl groups, as well as the C 2p(σ*)-derived states in flat and wrinkle regions to clarify the relationship between the spatially-resolved chemical states and the magnetism of GO, M-rGO and H-rGO. The results of X-ray magnetic circular dichroism further support the finding that C 2p(σ*)-derived states are the main origin of the magnetism of GO. Based on experimental results and first-principles calculations, the variation in magnetic behavior from GO to M-rGO and to H-rGO is interpreted, and the origin of ferromagnetism is identified as the C 2p(σ*)- derived states that involve defects/vacancies rather than the C 2p(π*) states that are bound with oxygen-containing and hydroxyl groups on GO sheets.[[notice]]補正完