Microscopically Computing Free-energy Profiles and Transition Path Time of Rare Macromolecular Transitions

We introduce a rigorous method to microscopically compute the observables which characterize the thermodynamics and kinetics of rare macromolecular transitions for which it is possible to identify a priori a slow reaction coordinate. In order to sample the ensemble of statistically significant reaction pathways, we define a biased molecular dynamics (MD) in which barrier-crossing transitions are accelerated without introducing any unphysical external force. In contrast to other biased MD methods, in the present approach the systematic errors which are generated in order to accelerate the transition can be analytically calculated and therefore can be corrected for. This allows for a computationally efficient reconstruction of the free-energy profile as a function of the reaction coordinate and for the calculation of the corresponding diffusion coefficient. The transition path time can then be readily evaluated within the Dominant Reaction Pathways (DRP) approach. We illustrate and test this method by characterizing a thermally activated transition on a two-dimensional energy surface and the folding of a small protein fragment within a coarse-grained model.Comment: Accepted for publication in Physical Review

Transition Path Theory from Biased Simulations

Transition Path Theory (TPT) provides a rigorous framework to investigate the dynamics of rare thermally activated transitions. In this theory, a central role is played by the forward committor function q^+(x), which provides the ideal reaction coordinate. Furthermore, the reactive dynamics and kinetics are fully characterized in terms of two time-independent scalar and vector distributions. In this work, we develop a scheme which enables all these ingredients of TPT to be efficiently computed using the short non-equilibrium trajectories generated by means of a specific combination of enhanced path sampling techniques. In particular, first, we further extend the recently introduced Self-Consistent Path Sampling (SCPS) algorithm in order to compute the committor q^+(x). Next, we show how this result can be exploited in order to define efficient algorithms which enable us to directly sample the transition path ensemble.Comment: Version accepted for publication in J. Chem. Phy

Quantum Interactions Between Non-Perturbative Vacuum Fields

We develop an approach to investigate the non-perturbative dynamics of quantum field theories, in which specific vacuum field fluctuations are treated as the low-energy dynamical degrees of freedom, while all other vacuum field configurations are explicitly integrated out from the path integral. We show how to compute the effective interaction between the vacuum field degrees of freedom both perturbatively (using stochastic perturbation theory) and fully non-perturbatively (using lattice field theory simulations). The present approach holds to all orders in the couplings and does not rely on the semi-classical approximation.Comment: 15 pages, 4 figure

The Role of Non-native Interactions in the Folding of Knotted Proteins

Stochastic simulations of coarse-grained protein models are used to investigate the propensity to form knots in early stages of protein folding. The study is carried out comparatively for two homologous carbamoyltransferases, a natively-knotted N-acetylornithine carbamoyltransferase (AOTCase) and an unknotted ornithine carbamoyltransferase (OTCase). In addition, two different sets of pairwise amino acid interactions are considered: one promoting exclusively native interactions, and the other additionally including non-native quasi-chemical and electrostatic interactions. With the former model neither protein show a propensity to form knots. With the additional non-native interactions, knotting propensity remains negligible for the natively-unknotted OTCase while for AOTCase it is much enhanced. Analysis of the trajectories suggests that the different entanglement of the two transcarbamylases follows from the tendency of the C-terminal to point away from (for OTCase) or approach and eventually thread (for AOTCase) other regions of partly-folded protein. The analysis of the OTCase/AOTCase pair clarifies that natively-knotted proteins can spontaneously knot during early folding stages and that non-native sequence-dependent interactions are important for promoting and disfavoring early knotting events.Comment: Accepted for publication on PLOS Computational Biolog