Searching low-energy conformations of two elastin sequences

The three-dimensional structures of two common repeat motifs Val$^1$-Pro$^2$-Gly$^3$-Val$^4$-Gly$^5$ and Val$^1$-Gly$^2$-Val$^3$-Pro$^4$-Gly$^5$-Val$^6$-Gly$^7$-Val$^8$-Pro$^9$ of tropoelastin are investigated by using the multicanonical simulation procedure. By minimizing the energy structures along the trajectory the thermodynamically most stable low-energy microstates of the molecule are determined. The structural predictions are in good agreement with X-ray diffraction experiments.Comment: 16 pages, 5 figure

Simple flexible polymers in a spherical cage

We report the results of Monte Carlo simulations investigating the effect of a spherical confinement within a simple model for a flexible homopolymer. We use the parallel tempering method combined with multi-histogram reweighting analysis and multicanonical simulations to investigate thermodynamical observables over a broad range of temperatures, which enables us to describe the behavior of the polymer and to locate the freezing and collapse transitions. We find a strong effect of the spherical confinement on the location of the collapse transition, whereas the freezing transition is hardly effected.Comment: 7 pages, 4 figure

Structure of the Energy Landscape of Short Peptides

We have simulated, as a showcase, the pentapeptide Met-enkephalin (Tyr-Gly-Gly-Phe-Met) to visualize the energy landscape and investigate the conformational coverage by the multicanonical method. We have obtained a three-dimensional topographic picture of the whole energy landscape by plotting the histogram with respect to energy(temperature) and the order parameter, which gives the degree of resemblance of any created conformation with the global energy minimum (GEM).Comment: 17 pages, 4 figure

Accurate implementation of leaping in space: The spatial partitioned-leaping algorithm

There is a great need for accurate and efficient computational approaches that can account for both the discrete and stochastic nature of chemical interactions as well as spatial inhomogeneities and diffusion. This is particularly true in biology and nanoscale materials science, where the common assumptions of deterministic dynamics and well-mixed reaction volumes often break down. In this article, we present a spatial version of the partitioned-leaping algorithm (PLA), a multiscale accelerated-stochastic simulation approach built upon the tau-leaping framework of Gillespie. We pay special attention to the details of the implementation, particularly as it pertains to the time step calculation procedure. We point out conceptual errors that have been made in this regard in prior implementations of spatial tau-leaping and illustrate the manifestation of these errors through practical examples. Finally, we discuss the fundamental difficulties associated with incorporating efficient exact-stochastic techniques, such as the next-subvolume method, into a spatial-leaping framework and suggest possible solutions.Comment: 15 pages, 9 figures, 2 table

A "partitioned leaping" approach for multiscale modeling of chemical reaction dynamics

We present a novel multiscale simulation approach for modeling stochasticity in chemical reaction networks. The approach seamlessly integrates exact-stochastic and "leaping" methodologies into a single "partitioned leaping" algorithmic framework. The technique correctly accounts for stochastic noise at significantly reduced computational cost, requires the definition of only three model-independent parameters and is particularly well-suited for simulating systems containing widely disparate species populations. We present the theoretical foundations of partitioned leaping, discuss various options for its practical implementation and demonstrate the utility of the method via illustrative examples.Comment: v4: 12 pages, 5 figures, final accepted version. Error found and fixed in Appendi

Maximizing Maximal Angles for Plane Straight-Line Graphs

Let $G=(S, E)$ be a plane straight-line graph on a finite point set $S\subset\R^2$ in general position. The incident angles of a vertex $p \in S$ of $G$ are the angles between any two edges of $G$ that appear consecutively in the circular order of the edges incident to $p$. A plane straight-line graph is called $\phi$-open if each vertex has an incident angle of size at least $\phi$. In this paper we study the following type of question: What is the maximum angle $\phi$ such that for any finite set $S\subset\R^2$ of points in general position we can find a graph from a certain class of graphs on $S$ that is $\phi$-open? In particular, we consider the classes of triangulations, spanning trees, and paths on $S$ and give tight bounds in most cases.Comment: 15 pages, 14 figures. Apart of minor corrections, some proofs that were omitted in the previous version are now include

Epigenetic Chromatin Silencing: Bistability and Front Propagation

The role of post-translational modification of histones in eukaryotic gene regulation is well recognized. Epigenetic silencing of genes via heritable chromatin modifications plays a major role in cell fate specification in higher organisms. We formulate a coarse-grained model of chromatin silencing in yeast and study the conditions under which the system becomes bistable, allowing for different epigenetic states. We also study the dynamics of the boundary between the two locally stable states of chromatin: silenced and unsilenced. The model could be of use in guiding the discussion on chromatin silencing in general. In the context of silencing in budding yeast, it helps us understand the phenotype of various mutants, some of which may be non-trivial to see without the help of a mathematical model. One such example is a mutation that reduces the rate of background acetylation of particular histone side-chains that competes with the deacetylation by Sir2p. The resulting negative feedback due to a Sir protein depletion effect gives rise to interesting counter-intuitive consequences. Our mathematical analysis brings forth the different dynamical behaviors possible within the same molecular model and guides the formulation of more refined hypotheses that could be addressed experimentally.Comment: 19 pages, 5 figure

Vibrational Enhancement of the Effective Donor - Acceptor Coupling

The paper deals with a simple three sites model for charge transfer phenomena in an one-dimensional donor (D) - bridge (B) - acceptor (A) system coupled with vibrational dynamics of the B site. It is found that in a certain range of parameters the vibrational coupling leads to an enhancement of the effective donor - acceptor electronic coupling as a result of the formation of the polaron on the B site. This enhancement of the charge transfer efficiency is maximum at the resonance, where the effective energy of the fluctuating B site coincides with the donor (acceptor) energy.Comment: 5 pages, 3 figure