Coarse Graining in Conformational Space. Feliks Nüske and Cecilia Clementi Rice University Sep 14, 2017
|
|
- Alyson Ellis
- 5 years ago
- Views:
Transcription
1 Coarse Graining in Conformational Space and Cecilia Clementi Rice University Sep 14, 2017
2 Coarse-Graining coarse-graining in conformation space coarse-graining in structural space
3 Biophysical systems are often characterized by different timescales, separated by large gaps Example: helix-coil transition coil helix typical timescales: helix-coil transition t slow µs diffusion within a state t fast ns
4 Reaction Coordinates should capture the (rare) barrier-crossing events. Courtesy of Michele Ceriotti
5 I. Optimal Reaction Coordinates
6 Basic Assumptions Dynamical Model is a Markov process X t,t 0 The process is ergodic if spatial and temporal averages coincide: E µ (f(x 0,X,...,X L )) = = lim K!1 K 1 L K L 1 X k=0 f(x k ),...,X (k+l) ) The process is reversible if there is no preferred direction: P µ (X 0 2 A, X t 2 B) =P µ (X 0 2 B,X t 2 A).
7 Basic Assumptions Dynamical Model is a Markov process X t,t 0 A distribution is stationary if it is invariant under the dynamics: P µ (X t 2 A) =µ(a) For an ergodic process, there is a unique stationary distribution.
8 Transfer Operators Transfer Operators describe the evolution of probability distributions: Z Z A [T t f](x)µ(dx) = S p t (x, A)f(x)µ(dx). Here, A is a set, µ is a measure, f is a density w.r.t. µ, and p t (x, A) is the stochastic transition function. Lasota, Mackey: Chaos, Fractals and Noise, Springer (1993)
9 Properties If µ is stationary, transfer op. are defined on all spaces L p µ. Semigroup property: T s+t = T t T s,s,t 0. For reversible processes, the transfer operators are self-adjoint on. L 2 µ Lasota, Mackey: Chaos, Fractals and Noise, Springer (1993)
10 Spectral Decomposition Using the infinitesimal Generator: Lf =lim t!0 T t f f t and its spectral components:, L j = apple j j we have the spectral mapping theorem: e (L) (T ) e (L) [ {0}, Pazy: Semigroups of linear operators and applications to partial differential equations, Springer (1983) and the spectral decomposition: T f = 1X j=1 e apple j hf, j i µ j
11 Slow Processes In many biophysical applications, we have: such that, for 0=apple 0 < apple 1 <...<apple M apple M+1 1 apple M+1 T f MX j=1 e apple j hf, j i µ j The presence of separated timescales is directly linked to the structure of the spectrum.
12 Slow Processes Using the spectral decomposition, we find for the auto-correlation function: E t (f(x t )f(x t+ )) = E µ (f(x 0 )f(x )) = ht f,fi µ 1X = e applej hf, j i 2 µ j=1
13 Variational Theorem Variational Theorem for self-adjoint operators: MX e apple j =sup j=1 j=1 MX ht f j,f j i µ, =sup MX j=1 E t (f j (X t )f j (X t+ )), hf j,f j 0i µ = j,j 0. Dominant eigenfunctions are slowest to de-correlate. They are optimal descriptors of slow dynamics. Noé and Nüske, SIAM Multiscale Model. Simul. (2013)
14 Galerkin Projection Proposition 3.2. Let D be a space of N linearly independent ansatz functions f i, i = 1,..., N. The set of M Nmutuallyorthonormalfunctionsf am, m = 1,..., M which maximize the Rayleigh trace restricted to D, is given bythefirstmeigenvectors of the generalized eigenvalue problem where the matrices C τ, C 0 are given by C τ a m = ˆλ m C 0 a m, (3.1.7) C τ (i, j) = T τ f i, f j µ (3.1.8) C 0 (i, j) = f i, f j µ. (3.1.9) Noé and Nüske, SIAM Multiscale Model. Simul. (2013)
15 Galerkin Projection By ergodicity of the process: Theorem 3.4. Let τ = L t beanintegermultipleofthediscretetimestep. Ifthe Markov process is initially distributed according to the unique invariant measure µ, then for P t -a.s. all trajectories (X k t ) k=0,wehave: C τ (i, j) = lim K 1 K L C 0 (i, j) = lim K 1 K L K L 1 k=0 K L 1 k=0 f i (X k t ) f j (X (k+l) t ), (3.2.1) f i (X k t ) f j (X k t ). (3.2.2) Noé and Nüske, SIAM Multiscale Model. Simul. (2013)
16 MSMs Markov State Models (MSM): Apply previous formalism with a basis of step functions: Schütte et. al. J. Comput. Phys. (1999) Prinz et. al. J. Chem. Phys. (2011)! c ij = 1 T i(x) = X T k=1 i(x k ) j (X k+ ) Normalize to stochastic matrix: T( ) = Schütte and Sarich, Metastability and Markov State Models in Molecular Dynamics: Modeling, Analysis, Algorithmic Approaches, Courant Lecture Notes (2013) Bowman et. al. (Eds) An Introduction to Markov State Models and Their Application to Long Timescale Molecular Simulation, Springer (2014) ( 1, x2 S i 0, else cij c i ij
17 MSMs Identification of metastable states: and the reaction field method for computing nonbonded forces. The simulation setup is described in detail in the SI Appendix. The simulated structures were aligned onto the native structure and then clustered finely into 1,734 kinetically connected and wellpopulated clusters. A transition matrix T(τ) was constructed by counting transitions between these clusters at a lag time of τ = 2 ns (see Theory). It was verified that T(τ) is a good model for the long-time kinetics (details on the Markov model construction and validation are given in the SI Appendix). All properties computed from the Markov model are associated with statistical uncertainty resulting from the fact that only a finite amount of simulation data has been used to construct the model. These uncertainties are computed by using a Bayesian inference method Transition Path Theory 6 h-bonds in the native state, has 3 h-bonds in the native sta the transition matrix T(τ) betw tor probabilities and the foldin in Theory. In order to obtain a view unbiased by defining reaction must be considered individua decomposed into individual p of them the times when hairp were computed. Formation average number of hydrogen b state, but variations of this th in high-dimensional spaces where binning of all relevant coordithe space spanned by the US coordinate and the second IC was discretized into 100 Voronoi cells with the k-means algorithm (Fig. nates is not an option. We first analyze pure US simulations with 150 umbrella windows 5A). This number of states is far smaller than the number of bins used to sample the position of benzamidine between the bound pose that would be required with a binned estimator such as WHAM or discrete TRAM. A count matrix ckij was estimated for every umand a prebinding site (Fig. 5A, structures i iv; details in Appendix). brella at a lag time of 11 ns, and the largest strongly To detect rare events in the unbiased coordinates, time-lagged P connected component S of the summed count matrix cij = k ckij was deindependent component analysis (TICA) (29, 30) was used with termined. The initial set was strongly disconnected, and we therethe Cartesian coordinates of residues around the binding site. fore adaptively started new umbrella simulations in nine rounds, to The first independent component (IC) is strongly correlated with improve the connectivity (Appendix). In the complete dataset, some the US coordinate. From the remaining ICs, two had timescales clusters are still disconnected (red clusters in Fig. 5A). implied by the TICA eigenvalues larger than the trajectory In particular, these disconnected states include structures in length, indicating undesirable metastable transitions orthogonal which the binding site is occluded by a tryptophan side chain, while to the umbrella coordinate. The second IC corresponds to benzamidine is still inside, and structures in which the binding site closing of the binding pocket by the Trp 215 side chain (Fig. 5D, attempts to close during the exit pathway. TRAM is applied on the structures i and iii). The third IC corresponds to an isomerization connected subset of states (white clusters in Fig. 5A). The TRAM of the disulfide bond between Cys 191 and Cys 220. An analysis results show that the Trp-occluded conformation is a local miniusing MBAR or WHAM is thus unfeasible or inefficient, as the mum inofthe free-energy landscape This is confirmed by ustrative Markov state model analysis of a pentapeptide. (a) Projection a trajectory onto the(fig. slow5c). collective coordinates (independent global equilibrium assumption is strongly violated. refs. 13 and 61 where the Trp-occluded conformation is shown to Markov model, ICs) shows rare transitions between different metastable states. (b) Implied relaxation time scales as a function of the One strategy to deal with this very common problem is to re- be a metastable conformation of the protein. In contrast, this local hows that strain a lag coordinates time τ = 1orthogonal ns is suitable. Shown coordinate, as shadedtoregions 95% confidence intervals. (c) several A Chapman Kolmogorov test shows to the umbrella avoid are minimum is not found by MBAR, and disconnected minundesirable degrees of freedom switchingon (45). Although this scales. 1 ns model accurately predicts thefrom behavior longer time (d)spuriously The free energy(boxes landscape is computed from the MSM as a ima are estimated in Fig. 5B). approach is useful for computing energy differences between end To analyze full high-dimensional binding mechanism and the two slowest ICs. The metastable states can be visually distinguished as freethe energy minima. (e) Probability distributions are given for states, it may change or restrain the transition mechanism and estimate unbiased kinetics, we must go beyond US simulations. gest livingartificially metastable states. The resulting assignment between energy minima and state numbers is shown in panel d. (f) The rate model increase free-energy barriers along the pathway. With We therefore used TRAM to combine the US data with up to from a hidden model based of the Rates areunbiased given MD in nanoseconds; 10 structures have been sampled for TICA Markov and TRAM, we now have coarse-graining the possibility to allow thesemsm. or μs of data (details in Appendix). The unbiased started the unbound state, such that many dynamics to happen,shown and to treat these events able state thogonal from the distributions in panel c. Allexplicitly. subfigurestrajectories except for the in molecular structure images havebinding been generated with Deuflhard and Weber, Linear Algebra and its Applications (2005) Multi-Ensemble MSMs: Molecular structures were rendered with VMD.134 D connected set disconnected (iii) or generalized ensemble simulations (e.g., replica-exchange MD), and combination of such data with direct MD via MSMs. Will contain estimators such as WHAM104 and TRAM.97,98,100 mtms (planned). Markov transition models.93 Continuous state kinetic models with a Markovian kernel, i.e., a means to compute the transition density p(x,y) or rate density k(x,y) between20000 two continuous60points of state space. IC2 (i) time scales 4as a function of some model parameter, (iii) representation of the kinetic model as a network. 3 (planned). Discrete state models with multiple ther(i) c states. This2 package is currently developed separately thub.com/markovmodel/thermotools) and will be 1 in PyEMMA in the near future. It will enable the 0 data from biased simulations (e.g., umbrella sampling), 10 B US coordinate / Å C 18 C IC (iv) (ii) free energy / kt US coordinate / Å US coordinate / Å 0 (ii) (iv) / cgi / doi / / pnas Noé et al. PNAS (2009) E TRAM MSM % of unbiased data used Fig. 5. Thermodynamics and kinetics of all-atom protein ligand binding model for trypsin benzamidine. (A C) US simulations. (D and E) MEMM using both unbiased and US simulations. (A) Trajectories projected on the space of the umbrella sampling (US) coordinate and the second independent component (IC). The US coordinate describes a transition from benzamidine bound to Asp-189 to benzamidine located outside the binding pocket on the surface of trypsin. The second IC corresponds to concerted opening of loop (Trp-215-Gln-221) and flipping of Trp-215. The Voronoi centers of the Markov states are shows as disks. Markov states that are irreversibly connected to the data set are shown as red disks and are excluded from the MEMM. (B) Potential of mean force (PMF) in the same coordinate space computed with unbounddoi: /acs.jctc.5b00743 J. Chem. Theory Comput. XXXX, XXX, XXX XXX 1000 probab. of finding koff A5 Wu et al. PNAS (2016)
18 Diffusion Maps Courtesy of Ralf Banisch
19 Diffusion Maps Pick a symmetric, positive kernel k. Define a Markov chain by p(x,y) = k(x,y) d(x). d(x) = X k(x,y)dµ(y) Introduce the diffusion distance: { } D t (x, y) 2 p t (x, ) p t (y, ) 2 L 2 (X,dµ/π) = By spectral expansion: ( ( D t (x, y) = ( that as { l 1 λ 2t l X ( pt (x, u) p t (y, u) ) 2 dµ(u) π(u). ( ψl (x) ψ l (y) ) ) In the following space, Euclidean distance equals diffusion dist. Ψ t (x) is constant, we have omitted the λ t 1 ψ 1(x) λ t 2 ψ 2(x).. λ t s(δ,t) ψ s(δ,t)(x). Coffman et. al. Appl. Comput. Harmon. Anal. (2006)
20 Diffusion Maps (1) Fix α R and a rotation-invariant kernel k ε (x, y) = h ( x y 2 ε (2) Let q ε (x) = k ε (x, y)q(y) dy ). X and form the new kernel k (α) ε (x, y) = k ε(x, y) q α ε (x)qα ε (y). (3) Apply the weighted graph Laplacian normalization to this kernel by setting d ε (α) (x) = k ε (α) (x, y)q(y) dy X and by defining the anisotropic transition kernel p ε, α (x, y) = k(α) ε (x, y) d (α) ε (x). Coffman et. al. Appl. Comput. Harmon. Anal. (2006)
21 Diffusion Maps Theorem 2. Let L ε,α = I P ε,α ε be the infinitesimal generator of the Markov chain. Then for a fixed K>0,wehaveonE K lim L ε,αf = (f q1 α ) ε 0 q 1 α (q1 α ) q 1 α f. For =0.5, this operator becomes the generator of a diffusion process Coffman et. al. Appl. Comput. Harmon. Anal. (2006)
22 II. Effective Dynamics on Reaction Coordinates
23 How can we define dynamics on the space of reaction coordinates? Assume the process is given by an SDE dx t = b(x t )+ 2β 1 σdb t. Then the generator is a differential operator: Lf(x) = n i=1 b i (x) f x i (x)+ 1 β n i=1 σ 2 ii 2 f x 2 i (x). Reaction coordinates are represented by smooth mapping: : R n 7! R m
24 Projection Formalism Projection operator: Pf(x) = 1 ν(z) ˆ Σ z f(x)µ(x)j 1/2 (x)dσ z (x). Projection is an orthogonal projection onto the space: H 0 = { } f L 2 µ : f = f(z) =f(ξ(x)) Zhang et. al., Faraday Discuss. (2016)
25 Projection Formalism Projected generator is again the generator of a diffusion: L ξ = PLP, m L ξ f(z) = P(Lξ l )(z) f l=1 + 1 β m l,r=1 P z l ( n i=1 σ 2 ii ξ l x i ξ r x i ) 2 f z l z r. Effective drift and diffusion can be estimated from data: e bl (z) = lim E l (x(s)) z l x(0) µ z, 1 apple l apple m, s!0+ s ea lk (z) = lim E ( l (x(s)) z l )( k (x(s)) z k ) x(0) µ z 2 s!0+ s Legoll and Lelièvre, Nonlinearity (2010) Zhang et. al., Faraday Discuss. (2016)
26 Approximation Quality Do these effective dynamics preserve the dominant timescales of the original process? If we use the optimal reaction coordinates, yes! More precisely, if the eigenfunctions only depend on z: j = j (z) the dominant timescales are exactly preserved: PLP j (z) = apple j j (z) Zhang et. al., Faraday Discuss. (2016)
27 Approximation Quality What if we only have an approximation? Proposition 1 If, for each dominant eigenfunction ψ i,i=2,...,m,thereisafunctiong i H 0 that approximates ψ i in H 1 -norm, i.e. g i ψ i 2 H 1 = g i ψ i 2 L 2 + N j=1 g i x j ψ i x j 2 L 2 ϵ2, (31) then the maximal relative error between the dominant eigenvalues of L ξ and L is bounded by max i=2,...,m ω i κ i ω i (M 1)κ 1/2 2 σ β 1 ϵ. (32) Nüske and Clementi (in preparation)
28 Numerical Example - Two-d landscape - Dynamics in x-direction is slow. - Study several rotations of the system. - Always use x as reaction coordinate. Nüske and Clementi (in preparation)
29 Numerical Example Nüske and Clementi (in preparation) Figure 36. Implied timescales of the one-dimensional effective dynamics obtained by projecting the dynamics Eq. (C9) onto the original coordinate x 0,fordifferent values of θ shown on the left. The columns correspond
30 Numerical Example II Use the polar angle as reaction coordinate. Nüske and Clementi (in preparation)
31 Numerical Example II Figure 26. Implied timescale tests for the effective dynamicssimulationsusingthedriftanddiffusion estimated at different values of the offset s (shown on top of each plot). The dashed lines correspond to the sixleading Nüske and Clementi (in preparation)
32 ACKNOWLEDGMENTS Prof. Cecilia Clementi Dr. Giovanni Pinamonti Lorenzo Boninsegna Alexander Kluber Justin Chen Eugen Hruska Our collaborators: Dr. Ralf Banisch (Freie Universität Berlin) Dr. Péter Koltai (Freie Universität Berlin) Prof. Frank Noé (Freie Universität Berlin)
arxiv: v2 [stat.ml] 22 Jan 2017
Variational Koopman models: slow collective variables and molecular kinetics from short off-equilibrium simulations Hao Wu, Feliks Nüske, Fabian Paul, Stefan Klus, Péter Koltai, and Frank Noé Department
More informationMulti-Ensemble Markov Models and TRAM. Fabian Paul 21-Feb-2018
Multi-Ensemble Markov Models and TRAM Fabian Paul 21-Feb-2018 Outline Free energies Simulation types Boltzmann reweighting Umbrella sampling multi-temperature simulation accelerated MD Analysis methods
More informationFinding slow modes and accessing very long timescales in molecular dynamics. Frank Noé (FU Berlin)
Finding slow modes and accessing very long timescales in molecular dynamics Frank Noé (FU Berlin) frank.noe@fu-berlin.de Cell Proteins McGufee and Elcock, PloS Comput Biol 2010 Protein-Protein binding
More informationRELIABLE APPROXIMATION OF LONG RELAXATION TIMESCALES IN MOLECULAR DYNAMICS
RELIABLE APPROXIMATION OF LONG RELAXATION TIMESCALES IN MOLECULAR DYNAMICS WEI ZHANG AND CHRISTOF SCHÜTTE,2 Abstract. Many interesting rare events in molecular systems like ligand association, protein
More informationSupplementary Figures:
Supplementary Figures: Supplementary Figure 1: The two strings converge to two qualitatively different pathways. A) Models of active (red) and inactive (blue) states used as end points for the string calculations
More informationEffective dynamics for the (overdamped) Langevin equation
Effective dynamics for the (overdamped) Langevin equation Frédéric Legoll ENPC and INRIA joint work with T. Lelièvre (ENPC and INRIA) Enumath conference, MS Numerical methods for molecular dynamics EnuMath
More informationMarkov State Models for Rare Events in Molecular Dynamics
Entropy 214, 16, 258-286; doi:1.339/e161258 OPEN ACCESS entropy ISSN 199-43 www.mdpi.com/journal/entropy Article Markov State Models for Rare Events in Molecular Dynamics Marco Sarich 1, *, Ralf Banisch
More informationVariational Tensor Approach for Approximating the Rare-Event Kinetics of Macromolecular Systems
TECHNISCHE UNIVERSITÄT BERLIN Variational Tensor Approach for Approximating the Rare-Event Kinetics of Macromolecular Systems Feliks Nüske Reinhold Schneider Francesca Vitalini Frank Noé Preprint 2015/30
More informationThe Monte Carlo Computation Error of Transition Probabilities
Zuse Institute Berlin Takustr. 7 495 Berlin Germany ADAM NILSN The Monte Carlo Computation rror of Transition Probabilities ZIB Report 6-37 (July 206) Zuse Institute Berlin Takustr. 7 D-495 Berlin Telefon:
More informationMarch 13, Paper: R.R. Coifman, S. Lafon, Diffusion maps ([Coifman06]) Seminar: Learning with Graphs, Prof. Hein, Saarland University
Kernels March 13, 2008 Paper: R.R. Coifman, S. Lafon, maps ([Coifman06]) Seminar: Learning with Graphs, Prof. Hein, Saarland University Kernels Figure: Example Application from [LafonWWW] meaningful geometric
More informationSupplementary information for cloud computing approaches for prediction of ligand binding poses and pathways
Supplementary information for cloud computing approaches for prediction of ligand binding poses and pathways Morgan Lawrenz 1, Diwakar Shukla 1,2 & Vijay S. Pande 1,2 1 Department of Chemistry, Stanford
More informationDIFFUSION MAPS, REDUCTION COORDINATES AND LOW DIMENSIONAL REPRESENTATION OF STOCHASTIC SYSTEMS
DIFFUSION MAPS, REDUCTION COORDINATES AND LOW DIMENSIONAL REPRESENTATION OF STOCHASTIC SYSTEMS R.R. COIFMAN, I.G. KEVREKIDIS, S. LAFON, M. MAGGIONI, AND B. NADLER Abstract. The concise representation of
More informationAPPROXIMATING SELECTED NON-DOMINANT TIMESCALES
APPROXIMATING SELECTED NON-DOMINANT TIMESCALES BY MARKOV STATE MODELS Marco Sarich Christof Schütte Institut für Mathematik, Freie Universität Berlin Arnimallee 6, 95 Berlin, Germany Abstract We consider
More informationCompeting sources of variance reduction in parallel replica Monte Carlo, and optimization in the low temperature limit
Competing sources of variance reduction in parallel replica Monte Carlo, and optimization in the low temperature limit Paul Dupuis Division of Applied Mathematics Brown University IPAM (J. Doll, M. Snarski,
More informationLecture 12: Detailed balance and Eigenfunction methods
Miranda Holmes-Cerfon Applied Stochastic Analysis, Spring 2015 Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility),
More informationApplication of the Markov State Model to Molecular Dynamics of Biological Molecules. Speaker: Xun Sang-Ni Supervisor: Prof. Wu Dr.
Application of the Markov State Model to Molecular Dynamics of Biological Molecules Speaker: Xun Sang-Ni Supervisor: Prof. Wu Dr. Jiang Introduction Conformational changes of proteins : essential part
More informationDownloaded 11/16/15 to Redistribution subject to SIAM license or copyright; see
MULTISCALE MODEL. SIMUL. Vol. 8, No. 4, pp. 1154 1177 c 21 Society for Industrial and Applied Mathematics Downloaded 11/16/15 to 162.15.68.31. Redistribution subject to SIAM license or copyright; see http://www.siam.org/journals/ojsa.php
More informationUsing Spectral Clustering to Sample Molecular States and Pathways
Using Spectral Clustering to Sample Molecular States and Pathways Surl-Hee Ahn 1, a) 2, b) and Johannes Birgmeier 1) Chemistry Department, Stanford University, Stanford, California 94305, USA 2) Computer
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationAsymptotic stability of an evolutionary nonlinear Boltzmann-type equation
Acta Polytechnica Hungarica Vol. 14, No. 5, 217 Asymptotic stability of an evolutionary nonlinear Boltzmann-type equation Roksana Brodnicka, Henryk Gacki Institute of Mathematics, University of Silesia
More informationLearning gradients: prescriptive models
Department of Statistical Science Institute for Genome Sciences & Policy Department of Computer Science Duke University May 11, 2007 Relevant papers Learning Coordinate Covariances via Gradients. Sayan
More informationMetastable Conformations via successive Perron-Cluster Cluster Analysis of dihedrals
Konrad-Zuse-Zentrum für Informationstechnik Berlin Takustraße 7 D-14195 Berlin-Dahlem Germany F. CORDES, M. WEBER, J. SCHMIDT-EHRENBERG Metastable Conformations via successive Perron-Cluster Cluster Analysis
More informationarxiv: v2 [stat.ml] 11 Dec 2017
Variational approach for learning Markov processes from time series data Hao Wu 1, and Frank Noé 1, arxiv:1707.04659v2 [stat.ml] 11 Dec 2017 1: Freie Universität Berlin, Dept. of Mathematics and Computer
More informationThermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning
Thermostatic Controls for Noisy Gradient Systems and Applications to Machine Learning Ben Leimkuhler University of Edinburgh Joint work with C. Matthews (Chicago), G. Stoltz (ENPC-Paris), M. Tretyakov
More informationLecture 12: Detailed balance and Eigenfunction methods
Lecture 12: Detailed balance and Eigenfunction methods Readings Recommended: Pavliotis [2014] 4.5-4.7 (eigenfunction methods and reversibility), 4.2-4.4 (explicit examples of eigenfunction methods) Gardiner
More informationModeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase. Andy Long Kridsanaphong Limtragool
Modeling the Free Energy Landscape for Janus Particle Self-Assembly in the Gas Phase Andy Long Kridsanaphong Limtragool Motivation We want to study the spontaneous formation of micelles and vesicles Applications
More informationLecture 14: Advanced Conformational Sampling
Lecture 14: Advanced Conformational Sampling Dr. Ronald M. Levy ronlevy@temple.edu Multidimensional Rough Energy Landscapes MD ~ ns, conformational motion in macromolecules ~µs to sec Interconversions
More informationOptimal Data-Driven Estimation of Generalized Markov State Models for Non-Equilibrium Dynamics
computation Article Optimal Data-Driven Estimation of Generalized Markov State Models for Non-Equilibrium Dynamics Péter Koltai, ID, Hao Wu, Frank Noé and Christof Schütte, Department of Mathematics and
More information2012 NCTS Workshop on Dynamical Systems
Barbara Gentz gentz@math.uni-bielefeld.de http://www.math.uni-bielefeld.de/ gentz 2012 NCTS Workshop on Dynamical Systems National Center for Theoretical Sciences, National Tsing-Hua University Hsinchu,
More informationConstrained Ulam Dynamic Mode Decomposition: Approximation of the Perron-Frobenius Operator for Deterministic and Stochastic Systems
Constrained Ulam Dynamic Mode Decomposition: Approximation of the Perron-Frobenius Operator for Deterministic and Stochastic Systems Debdipta Goswami, Emma Thackray 2, and Derek A. Paley 3 Abstract Dynamical
More informationSTA 4273H: Statistical Machine Learning
STA 4273H: Statistical Machine Learning Russ Salakhutdinov Department of Computer Science! Department of Statistical Sciences! rsalakhu@cs.toronto.edu! h0p://www.cs.utoronto.ca/~rsalakhu/ Lecture 7 Approximate
More informationFree energy simulations
Free energy simulations Marcus Elstner and Tomáš Kubař January 14, 2013 Motivation a physical quantity that is of most interest in chemistry? free energies Helmholtz F or Gibbs G holy grail of computational
More informationFIRST PASSAGE TIMES AND RELAXATION TIMES OF UNFOLDED PROTEINS AND THE FUNNEL MODEL OF PROTEIN FOLDING
FIRST PASSAGE TIMES AND RELAXATION TIMES OF UNFOLDED PROTEINS AND THE FUNNEL MODEL OF PROTEIN FOLDING By WEI DAI A dissertation submitted to the Graduate School New Brunswick Rutgers, The State University
More informationDefinition and basic properties of heat kernels I, An introduction
Definition and basic properties of heat kernels I, An introduction Zhiqin Lu, Department of Mathematics, UC Irvine, Irvine CA 92697 April 23, 2010 In this lecture, we will answer the following questions:
More informationBeyond Scalar Affinities for Network Analysis or Vector Diffusion Maps and the Connection Laplacian
Beyond Scalar Affinities for Network Analysis or Vector Diffusion Maps and the Connection Laplacian Amit Singer Princeton University Department of Mathematics and Program in Applied and Computational Mathematics
More informationSupporting Information. RNA hairpin
Supporting Information RNA hairpin 1. Initial Configurations We started our ST simulations from two different initial configurations as shown in Figure S1: a near-native state and a random coil. The near-native
More informationIdentification and Analysis of Transition and Metastable Markov States
Identification and Analysis of Transition and Metastable Markov States Linda Martini, 1 Adam Kells, 1 Gerhard Hummer, 2 Nicolae-Viorel Buchete, 3 and Edina Rosta 1 1 King s College London, Department of
More informationTheory and applications of random Poincaré maps
Nils Berglund nils.berglund@univ-orleans.fr http://www.univ-orleans.fr/mapmo/membres/berglund/ Random Dynamical Systems Theory and applications of random Poincaré maps Nils Berglund MAPMO, Université d'orléans,
More informationBayesian reconstruction of free energy profiles from umbrella samples
Bayesian reconstruction of free energy profiles from umbrella samples Thomas Stecher with Gábor Csányi and Noam Bernstein 21st November 212 Motivation Motivation and Goals Many conventional methods implicitly
More informationStatistical mechanics of random billiard systems
Statistical mechanics of random billiard systems Renato Feres Washington University, St. Louis Banff, August 2014 1 / 39 Acknowledgements Collaborators: Timothy Chumley, U. of Iowa Scott Cook, Swarthmore
More informationStabilization and Acceleration of Algebraic Multigrid Method
Stabilization and Acceleration of Algebraic Multigrid Method Recursive Projection Algorithm A. Jemcov J.P. Maruszewski Fluent Inc. October 24, 2006 Outline 1 Need for Algorithm Stabilization and Acceleration
More informationReaction Pathways of Metastable Markov Chains - LJ clusters Reorganization
Reaction Pathways of Metastable Markov Chains - LJ clusters Reorganization Eric Vanden-Eijnden Courant Institute Spectral approach to metastability - sets, currents, pathways; Transition path theory or
More informationModeling Biological Systems Opportunities for Computer Scientists
Modeling Biological Systems Opportunities for Computer Scientists Filip Jagodzinski RBO Tutorial Series 25 June 2007 Computer Science Robotics & Biology Laboratory Protein: πρώτα, "prota, of Primary Importance
More information3.320 Lecture 23 (5/3/05)
3.320 Lecture 23 (5/3/05) Faster, faster,faster Bigger, Bigger, Bigger Accelerated Molecular Dynamics Kinetic Monte Carlo Inhomogeneous Spatial Coarse Graining 5/3/05 3.320 Atomistic Modeling of Materials
More informationLECTURE NOTE #11 PROF. ALAN YUILLE
LECTURE NOTE #11 PROF. ALAN YUILLE 1. NonLinear Dimension Reduction Spectral Methods. The basic idea is to assume that the data lies on a manifold/surface in D-dimensional space, see figure (1) Perform
More informationControlled Diffusions and Hamilton-Jacobi Bellman Equations
Controlled Diffusions and Hamilton-Jacobi Bellman Equations Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 2014 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationConformational Studies of UDP-GlcNAc in Environments of Increasing Complexity
John von Neumann Institute for Computing Conformational Studies of UDP-GlcNAc in Environments of Increasing Complexity M. Held, E. Meerbach, St. Hinderlich, W. Reutter, Ch. Schütte published in From Computational
More informationITERATED FUNCTION SYSTEMS WITH CONTINUOUS PLACE DEPENDENT PROBABILITIES
UNIVERSITATIS IAGELLONICAE ACTA MATHEMATICA, FASCICULUS XL 2002 ITERATED FUNCTION SYSTEMS WITH CONTINUOUS PLACE DEPENDENT PROBABILITIES by Joanna Jaroszewska Abstract. We study the asymptotic behaviour
More informationProtein Dynamics, Allostery and Function
Protein Dynamics, Allostery and Function Lecture 3. Protein Dynamics Xiaolin Cheng UT/ORNL Center for Molecular Biophysics SJTU Summer School 2017 1 Obtaining Dynamic Information Experimental Approaches
More informationThe parallel replica method for Markov chains
The parallel replica method for Markov chains David Aristoff (joint work with T Lelièvre and G Simpson) Colorado State University March 2015 D Aristoff (Colorado State University) March 2015 1 / 29 Introduction
More informationMetastability of Markovian systems
Metastability of Markovian systems A transfer operator based approach in application to molecular dynamics Vom Fachbereich Mathematik und Informatik der Freien Universität Berlin zur Erlangung des akademischen
More informationThe Kramers problem and first passage times.
Chapter 8 The Kramers problem and first passage times. The Kramers problem is to find the rate at which a Brownian particle escapes from a potential well over a potential barrier. One method of attack
More informationMarkovian Description of Irreversible Processes and the Time Randomization (*).
Markovian Description of Irreversible Processes and the Time Randomization (*). A. TRZĘSOWSKI and S. PIEKARSKI Institute of Fundamental Technological Research, Polish Academy of Sciences ul. Świętokrzyska
More informationBrownian Motion. An Undergraduate Introduction to Financial Mathematics. J. Robert Buchanan. J. Robert Buchanan Brownian Motion
Brownian Motion An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Background We have already seen that the limiting behavior of a discrete random walk yields a derivation of
More informationA mathematical framework for Exact Milestoning
A mathematical framework for Exact Milestoning David Aristoff (joint work with Juan M. Bello-Rivas and Ron Elber) Colorado State University July 2015 D. Aristoff (Colorado State University) July 2015 1
More informationFinite Frames and Graph Theoretical Uncertainty Principles
Finite Frames and Graph Theoretical Uncertainty Principles (pkoprows@math.umd.edu) University of Maryland - College Park April 13, 2015 Outline 1 Motivation 2 Definitions 3 Results Outline 1 Motivation
More informationSupporting Information
Supporting Information Boehr et al. 10.1073/pnas.0914163107 SI Text Materials and Methods. R 2 relaxation dispersion experiments. 15 NR 2 relaxation dispersion data measured at 1 H Larmor frequencies of
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationLECTURE 15 Markov chain Monte Carlo
LECTURE 15 Markov chain Monte Carlo There are many settings when posterior computation is a challenge in that one does not have a closed form expression for the posterior distribution. Markov chain Monte
More informationAdaptive simulation of hybrid stochastic and deterministic models for biochemical systems
publications Listing in reversed chronological order. Scientific Journals and Refereed Proceedings: A. Alfonsi, E. Cances, G. Turinici, B. Di Ventura and W. Huisinga Adaptive simulation of hybrid stochastic
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
More informationMD Thermodynamics. Lecture 12 3/26/18. Harvard SEAS AP 275 Atomistic Modeling of Materials Boris Kozinsky
MD Thermodynamics Lecture 1 3/6/18 1 Molecular dynamics The force depends on positions only (not velocities) Total energy is conserved (micro canonical evolution) Newton s equations of motion (second order
More informationIntroduction to asymptotic techniques for stochastic systems with multiple time-scales
Introduction to asymptotic techniques for stochastic systems with multiple time-scales Eric Vanden-Eijnden Courant Institute Motivating examples Consider the ODE {Ẋ = Y 3 + sin(πt) + cos( 2πt) X() = x
More informationComputing Reaction Rates Using Macro-states
Computing Reaction Rates Using Macro-states Badi Abdul-Wahid 1 Ronan Costaouec 2 Eric Darve 2 Michelle Feng 1 Jesús Izaguirre 1 1 Notre-Dame University 2 Stanford University Tuesday July 3rd 2012 14:30
More informationA Geometric Interpretation of the Metropolis Hastings Algorithm
Statistical Science 2, Vol. 6, No., 5 9 A Geometric Interpretation of the Metropolis Hastings Algorithm Louis J. Billera and Persi Diaconis Abstract. The Metropolis Hastings algorithm transforms a given
More informationFunctional central limit theorem for super α-stable processes
874 Science in China Ser. A Mathematics 24 Vol. 47 No. 6 874 881 Functional central it theorem for super α-stable processes HONG Wenming Department of Mathematics, Beijing Normal University, Beijing 1875,
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More informationTransition Paths of Met-Enkephalin from Markov State Modeling of a Molecular Dynamics Trajectory
pubs.acs.org/jpcb Transition Paths of Met-Enkephalin from Markov State Modeling of a Molecular Dynamics Trajectory Rahul Banerjee and Robert I. Cukier* Department of Chemistry, Michigan State University,
More informationContribution from: Springer Verlag Berlin Heidelberg 2005 ISBN
Contribution from: Mathematical Physics Studies Vol. 7 Perspectives in Analysis Essays in Honor of Lennart Carleson s 75th Birthday Michael Benedicks, Peter W. Jones, Stanislav Smirnov (Eds.) Springer
More informationUnfolding CspB by means of biased molecular dynamics
Chapter 4 Unfolding CspB by means of biased molecular dynamics 4.1 Introduction Understanding the mechanism of protein folding has been a major challenge for the last twenty years, as pointed out in the
More informationThe Continuity of SDE With Respect to Initial Value in the Total Variation
Ξ44fflΞ5» ο ffi fi $ Vol.44, No.5 2015 9" ADVANCES IN MATHEMATICS(CHINA) Sep., 2015 doi: 10.11845/sxjz.2014024b The Continuity of SDE With Respect to Initial Value in the Total Variation PENG Xuhui (1.
More informationDNA Bubble Dynamics. Hans Fogedby Aarhus University and Niels Bohr Institute Denmark
DNA Bubble Dynamics Hans Fogedby Aarhus University and Niels Bohr Institute Denmark Principles of life Biology takes place in wet and warm environments Open driven systems - entropy decreases Self organization
More informationSpectra of Large Random Stochastic Matrices & Relaxation in Complex Systems
Spectra of Large Random Stochastic Matrices & Relaxation in Complex Systems Reimer Kühn Disordered Systems Group Department of Mathematics, King s College London Random Graphs and Random Processes, KCL
More informationDeterminantal point processes and random matrix theory in a nutshell
Determinantal point processes and random matrix theory in a nutshell part I Manuela Girotti based on M. Girotti s PhD thesis and A. Kuijlaars notes from Les Houches Winter School 202 Contents Point Processes
More information25.1 Ergodicity and Metric Transitivity
Chapter 25 Ergodicity This lecture explains what it means for a process to be ergodic or metrically transitive, gives a few characterizes of these properties (especially for AMS processes), and deduces
More informationIdentification of slow molecular order parameters for Markov model construction
Identification of slow molecular order parameters for Markov model construction Guillermo Pérez-Hernández, Fabian Paul, Toni Giorgino, Gianni De Fabritiis, and Frank Noé Citation: The Journal of Chemical
More informationOverview & Applications. T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 04 June, 2015
Overview & Applications T. Lezon Hands-on Workshop in Computational Biophysics Pittsburgh Supercomputing Center 4 June, 215 Simulations still take time Bakan et al. Bioinformatics 211. Coarse-grained Elastic
More information215 Problem 1. (a) Define the total variation distance µ ν tv for probability distributions µ, ν on a finite set S. Show that
15 Problem 1. (a) Define the total variation distance µ ν tv for probability distributions µ, ν on a finite set S. Show that µ ν tv = (1/) x S µ(x) ν(x) = x S(µ(x) ν(x)) + where a + = max(a, 0). Show that
More informationMATH 590: Meshfree Methods
MATH 590: Meshfree Methods Chapter 2: Radial Basis Function Interpolation in MATLAB Greg Fasshauer Department of Applied Mathematics Illinois Institute of Technology Fall 2010 fasshauer@iit.edu MATH 590
More informationNumerical methods in molecular dynamics and multiscale problems
Numerical methods in molecular dynamics and multiscale problems Two examples T. Lelièvre CERMICS - Ecole des Ponts ParisTech & MicMac project-team - INRIA Horizon Maths December 2012 Introduction The aim
More informationEnergy Barriers and Rates - Transition State Theory for Physicists
Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle
More information1. Introduction Systems evolving on widely separated time-scales represent a challenge for numerical simulations. A generic example is
COMM. MATH. SCI. Vol. 1, No. 2, pp. 385 391 c 23 International Press FAST COMMUNICATIONS NUMERICAL TECHNIQUES FOR MULTI-SCALE DYNAMICAL SYSTEMS WITH STOCHASTIC EFFECTS ERIC VANDEN-EIJNDEN Abstract. Numerical
More informationWhy study protein dynamics?
Why study protein dynamics? Protein flexibility is crucial for function. One average structure is not enough. Proteins constantly sample configurational space. Transport - binding and moving molecules
More informationGeometry on Probability Spaces
Geometry on Probability Spaces Steve Smale Toyota Technological Institute at Chicago 427 East 60th Street, Chicago, IL 60637, USA E-mail: smale@math.berkeley.edu Ding-Xuan Zhou Department of Mathematics,
More informationp 1 ( Y p dp) 1/p ( X p dp) 1 1 p
Doob s inequality Let X(t) be a right continuous submartingale with respect to F(t), t 1 P(sup s t X(s) λ) 1 λ {sup s t X(s) λ} X + (t)dp 2 For 1 < p
More informationWeb Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D.
Web Appendix for Hierarchical Adaptive Regression Kernels for Regression with Functional Predictors by D. B. Woodard, C. Crainiceanu, and D. Ruppert A. EMPIRICAL ESTIMATE OF THE KERNEL MIXTURE Here we
More informationAPPLICATIONS OF THE KANTOROVICH-RUBINSTEIN MAXIMUM PRINCIPLE IN THE THEORY OF MARKOV OPERATORS
12 th International Workshop for Young Mathematicians Probability Theory and Statistics Kraków, 20-26 September 2009 pp. 43-51 APPLICATIONS OF THE KANTOROVICH-RUBINSTEIN MAIMUM PRINCIPLE IN THE THEORY
More informationNonparametric Drift Estimation for Stochastic Differential Equations
Nonparametric Drift Estimation for Stochastic Differential Equations Gareth Roberts 1 Department of Statistics University of Warwick Brazilian Bayesian meeting, March 2010 Joint work with O. Papaspiliopoulos,
More informationChapter 7. Markov chain background. 7.1 Finite state space
Chapter 7 Markov chain background A stochastic process is a family of random variables {X t } indexed by a varaible t which we will think of as time. Time can be discrete or continuous. We will only consider
More informationHierarchical Metastable States and Kinetic Transition Networks: Trajectory Mapping and Clustering
Hierarchical Metastable States and Kinetic Transition Networks: Trajectory Mapping and Clustering Xin Zhou xzhou@gucas.ac.cn Graduate University of Chinese Academy of Sciences, Beijing 2012.6.5 KITP, Santa
More informationDiffusion Maps, Spectral Clustering and Eigenfunctions of Fokker-Planck Operators
Diffusion Maps, Spectral Clustering and Eigenfunctions of Fokker-Planck Operators Boaz Nadler Stéphane Lafon Ronald R. Coifman Department of Mathematics, Yale University, New Haven, CT 652. {boaz.nadler,stephane.lafon,ronald.coifman}@yale.edu
More informationStochastic Particle Methods for Rarefied Gases
CCES Seminar WS 2/3 Stochastic Particle Methods for Rarefied Gases Julian Köllermeier RWTH Aachen University Supervisor: Prof. Dr. Manuel Torrilhon Center for Computational Engineering Science Mathematics
More informationMarkov State Models. Gregory R. Bowman Miller Fellow University of California, Berkeley
Markov State Models Gregory R. Bowman Miller Fellow University of California, Berkeley Tuesday, August 16, 2011 Defining The Protein Folding Problem Defining The Protein Folding Problem Why Mechanism?
More informationLinear algebra and applications to graphs Part 1
Linear algebra and applications to graphs Part 1 Written up by Mikhail Belkin and Moon Duchin Instructor: Laszlo Babai June 17, 2001 1 Basic Linear Algebra Exercise 1.1 Let V and W be linear subspaces
More informationA tutorial on numerical transfer operator methods. numerical analysis of dynamical systems
A tutorial on transfer operator methods for numerical analysis of dynamical systems Gary Froyland School of Mathematics and Statistics University of New South Wales, Sydney BIRS Workshop on Uncovering
More informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationStructural and mechanistic insight into the substrate. binding from the conformational dynamics in apo. and substrate-bound DapE enzyme
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 215 Structural and mechanistic insight into the substrate binding from the conformational
More informationLimitations of temperature replica exchange (T-REMD) for protein folding simulations
Limitations of temperature replica exchange (T-REMD) for protein folding simulations Jed W. Pitera, William C. Swope IBM Research pitera@us.ibm.com Anomalies in protein folding kinetic thermodynamic 322K
More information