Supplemental Material for Global Langevin model of multidimensional biomolecular dynamics
|
|
- Robert Poole
- 5 years ago
- Views:
Transcription
1 Supplemental Material for Global Langevin model of multidimensional biomolecular dynamics Norbert Schaudinnus, Benjamin Lickert, Mithun Biswas and Gerhard Stock Biomolecular Dynamics, Institute of Physics, Albert Ludwigs University, 74 Freiburg, Germany (Dated: September, 6) The first two authors contributed equally.
2 I. ALANINE DIPEPTIDE (ALAD) A. MD simulation details We conducted MD simulations of alanine dipeptide (Ac-Ala-NHCH 3, AlaD) at 3K, using the AMBER ff99sb force field in the GROMACS program suite. The structures were immersed in a dodecahedral box of TIP3P 3 water molecules such that there was at least Å of separation between the solute and the edge of the box. Periodic boundary conditions were used and the long range electrostatics was treated with the particle mesh Ewald method 4 with a Fourier spacing of. Å and sixth order interpolation to compute potential and forces between the grid points. For both the van der Waals and Coulomb interactions the cut-off was set to Å. All bonds were constrained using the LINCS algorithm. 5 The integration time step was fs. and coordinates were saved every steps during the simulations. A high frequency for saving the frames was necessary to provide sufficient data for Langevin modeling. First the water molecules were energy minimized using the steepest decent method. To equilibrate the system at the target temperature and to find an optimum box size, the system was simulated for ns in a NPT ensemble. Pressure was kept constant with a Berendsen barostat 6 having coupling constant τ P = ps and a compressibility of β = bar. For temperature control, the Bussi-Donadio-Parinello thermostat 7 with a coupling time τ T = ps was used. The last 5 ns of the NPT run was used to calculate the mean size of the water box. Finally, the production run for the MD simulations was performed in the NVT ensemble. B. Definition of reaction coordinate The backbone dihedral angles φ, ψ are often used to describe conformational dynamics of peptides. The Ramachandran plot of AlaD in Fig. Sa reveals four distinct metastable states: Left- and right-handed helical conformation α L and α R, the β sheet and the polyproline II state P II. While the most frequent conformational transition occurs between β sheet and P II, the α L state is very sparsely visited. In the course of our 5 ns MD simulation we observe only 3 transitions. Hence, in order to establish a simple one-dimensional model we restrained our system to the framed area in Fig. Sa, by only considering parts of the trajectory that stay in that region. Since in this area the dynamics of the coordinates φ and ψ are in good approximation uncorrelated, we focused on the essential coordinate ψ. In order to achieve a path-connected coordinate space, we shifted our system cyclically along the φ direction by φ = 5. Along the ψ direction, we excluded the very few transitions across the highest free energy barrier at ψ = to construct a simple linear representation of the cyclic coordinate. For convenience, we introduced the shifted coordinate x = ψ 7.
3 An excerpt of the resulting MD time series x(t) is shown in Fig. Sb. ψ - β PII α R α L a 8 4 x - b φ t [ns] FIG. S: (a) Ramachandran plot of alanine dipeptide. We focus on the marked area and impose cyclic shifts φ = 5 and ψ = 7 to achieve a path-connected collective coordinate. (b) The resulting time series resolves the transition α R β/p II along the shifted coordinate x = ψ 7. C. Calculation of first passage time and transition path time The free energy of AlaD projected on x = ψ 7 shows two minima centered at x 83 and x 85, respectively, which correspond to the states α R and β/p II, see Fig. Sa. To define the extension of the states, we chose a radius of r = 5, i.e. x i ± r define the borders of the states. Based on the definition of the states (centers and radii), we follow the trajectory after exiting one state till it enters the other state which directly yields the transition path time τ tp. To calculate the first passage time τ FP, we first assign all frames along the trajectory without a state definition to the last state visited. Subsequently, for each frame in the initial state, the first passage time τ FP denotes the number of steps it takes to reach the other state. 3
4 D. Properties of the noise model Using various time steps δt for the propagation of the dle, Fig. S4 shows (a) the autocorrelation function and (b) mean and standard deviation of the noise (for calculation, seeref. 8 ). Whilealltimestepsproducezeromeanandunitstandarddeviation,weneedδt fs to have the autocorrelation function decay to (almost) zero in the first time step (which reflects delta-correlated white noise). Note that the behavior of the noise autocorrelation function for longer times does not affect the dle simulations. noise autocorrelation a 4fs fs fs 4fs noise mean/sd 3 - b 4fs fs fs 4fs timestep ψ FIG. S: Properties of the noise model of AlaD, obtained for various time steps δt of the dle simulations. Shown are (a) autocorrelation function and (b) mean and standard deviation of the noise. 4
5 II. REPRESENTATION OF MULTIDIMENSIONAL FIELDS As explained in the main text, we obtain from the dle various fields, e.g., the friction ˆΓ(x)andthenoiseamplitude ˆK(x). Duetofinitestatistics, thesefieldestimatesaretypically noisy, in particular for multivariate observables. Since the construction of a global model requires sufficiently smooth and well sampled fields, in the following we propose two simple algorithms that smoothen and extrapolate a given field. As an illustrative example, we consider a scalar field, b(x,x ), depending on two variables. By using some suitable grid (i,j) the field is histogrammed into bins, which gives for each bin n(i,j) raw data points with an average value b org (i,j). Smoothing of this field is achieved by employing a weighted average over all neighboring bins that contain at least a single entry, b(i,j) = n tot i+ j+ k=i l=j n(k,l)b org (k,l), () where n tot = i+ j+ k=i l=j n(k,l) denotes the total number of data points. The procedure may be repeated several times in order to achieve sufficiently smooth fields. We applied the proposed algorithm to two-dimensional projections of the dle field estimates for the three-dimensional model of heptaalanine (Ala 7 ). As a representative example, Fig. S3 shows the friction fields Γ (x,x ) and Γ (x,x ) before and after smoothing. Note that the shape and the value range are preserved by the algorithm. By comparing the value ranges, the figure also reveals that the off-diagonal part of the friction matrix Γ is much smaller that the diagonal part. In a second step, one may improve the representation of sparsely sampled regions by replacing the field values of empty bins by the average trend of the surrounding bins. That is, given some empty bin with n(i,j) =, we scan its neighborhood for nonempty bins b(k,l) which satisfy k i and l j. In order to estimate the value of b(i,j), we make first-order Taylor expansions for the neighboring points. E.g., starting at point (i, j) and extrapolating in the direction of the first index, this approximation reads b(i,j)+(b(i,j) b(i,j)). Averaging over all possible directions, the estimate for b(i,j) can be written as b(i,j) = A i+ j+ k=i l=j b(k,l)+b ijkl, b ijkl = b(i+k i, j+l j) δ ijkl, () where δ ijkl δ(n(k,l))δ(n(i+k i,j+l j)) assures that the bins n(k,l) and n(i+k i,j+l j) are not empty and A = i+ j+ k=i l=j δ ijkl denotes the number of contributing terms. Subsequent to the extrapolation, another smoothening step is applied as described above. 5
6 a x b x x c x d x. x x - x FIG. S 3: Effect of the smoothening algorithm demonstrated for friction field Γ of the threedimensional model of Ala7. The diagonal element Γ (x, x ) is shown (a) from a binning of the dle estimates and (b) after applying the smoothening algorithm to it. Similarly, the off-diagonal element Γ (x, x ) is shown (c) before and (d) after smoothening. 6
7 III. HEPTAALANINE (ALA 7 ) A. Definition of reaction coordinates As for AlaD, the dihedral angles ψ i reveal the essential dynamical information of Ala 7. So, we again chose shifted coordinates x i = ψ i 7 (i =,,3) of the three inner residues as reaction coordinates. For purposes of presentation (see Fig. 4a in the main text), we rotated our data by α = 45 clockwise around the axis v = (,,) T. This rotation is described by x = R(α) x with 3 sin ( α)+cos(α) 3 sin ( α)+ 3 sin(α) 3 sin ( α) 3 sin(α) R(α) = 3 sin ( α) 3 sin(α) 3 sin ( α)+cos(α) 3 sin ( α)+ 3 sin(α). (3) 3 sin ( α)+ 3 sin(α) 3 sin ( α) 3 sin(α) 3 sin ( α)+cos(α) B. EVB model parameters The EVB model of the free energy of Ala 7 [Eq. (6) in the main text] is described by the energy matrix V(x,x,x 3 ) with diagonal elements (in units of k B T) V (x,x,x 3 ) =.43(x +43.) +.445(x +3.6) +.44(x ) V (x,x,x 3 ) =.48(x 4.59) +.94(x 38.8) +.99(x ) V 33 (x,x,x 3 ) =.5(x +9.6) +.49(x 64.65) +.97(x ) V 44 (x,x,x 3 ) =.6(x +5.3) +.55(x 6.9) +.4(x ) V 55 (x,x,x 3 ) =.73(x +68.9) +.659(x 8.4) +.646(x ) V 66 (x,x,x 3 ) =.6(x.3) +.6(x 97.) +.6(x ) +7.4 V 77 (x,x,x 3 ) =.95(x 33.4) +.5(x ) +.46(x ) V 88 (x,x,x 3 ) =.56(x 79.3) +.9(x +84.5) +.54(x ) and the off-diagonal part V coup =
8 C. Properties of the noise model Using various time steps δt for the propagation of the dle, Fig. S?? shows (a) the autocorrelation function and (b) mean and standard deviation of the noise. The minimum of δt = ps is determined by the write-out time of the MD simulation, as a maximum we choose δt = ps. In all cases, the noise is found to be delta-correlated and shows zero mean and unit standard deviation noise autocorrelation ps ps 4 ps ps a timestep 4 noise mean/stdev ps ps 4 ps ps b x FIG. S4: Properties of the noise model of Ala 7, obtained for various time steps δt of the dle simulations. Shown are (a) autocorrelation function and (b) mean and standard deviation of the noise. 8
9 D. Calculation of transition rates To describe the conformational dynamics of the three-dimensional model of Ala 7, we calculate the transition rates k i j between the eight states of the system. As written in the main text, the eight states can be labeled as abc with a,b,c = α,β in line with the structures of the three inner residues of Ala 7. In the rotated space (see above), the eight states are centered at x (ααα) = (,, ) T, x (βαα) = (45,, 4) T, x (αβα) = ( 4,45, ) T, x (ββα) = (,3, 6) T, x (αββ) = ( 6,,3) T, x (βββ) = (8,8,8) T, x (ααβ) = (, 4,45) T and x (βαβ) = (3, 6, ) T. Approximating the shape of the states by a sphere, we furthermore choose a radius r = 5 to define the region of each state. By knowing centers and widths of these states, a state trajectory n(t) was calculated that assigns a state n to every MD frame. Based on this definition, we counted the number of transitions N i j between state i (start) and state j (ending) for all states which occur within a certain lag time τ lag. We chose τ lag = ps, which is also the time step of the MD and dle (and hence the most natural choice). The transition probability i j for the transition p i j is then p i j = N i j / k N i k, (5) with the summation over k being performed over all eight ending states (so, i = k is included too). Then, the transition rates k i j can be calculated by k i j = p i j /τ lag. (6) Errors are estimated by cutting the trajectory into pieces and calculating the standard deviation of the rates for the different pieces. 9
10 E. Modifications of the model As explained in the main text, we constructed a modified model of Ala 7, in which the four barriers are increased over which transitions of the central angle occur, i.e., aαb aβb where a,b {α,β}. Fig. S5a shows the resulting free energy landscape of the modified EVB model. As a further illustration, Fig. S5b shows the outcome of a mle simulation starting in the lower part of the free energy landscape. Due to the increase of the barriers, the dynamics is restricted to the lower part. x x a x x b FIG. S5: (a) Free energy landscape of a modified EVB model of Ala 7. (b) Due to the high barriers between the upper and lower part of the energy landscape, a mle mle simulation starting in the lower part remains restricted to the lower part. F. Calculation of mean first passage times To calculate mean first passage times τ MFP of Ala 7 listed in Table I of the main text, we need to generalize the simple method for two states (see Sec. IC) to the case of several states. As a simple example, we consider three states a,b,c. To calculate the mean first passage time from a c, we consider an example trajectory with total time T taken to reach from a to c and is written as: T = t a, +t b, +t a, ; where t m,i is the lifetime in state m during the i-th visit. Since a trajectory can start in state a in any one of the multiple visits, we obtain a distribution of passage times from a c. The mean first passage time τ MFP is the probability of starting a trajectory during the i-th visit times the averaged first passage time from this visit to c. In the example above, we find t a, τ MFP = ( t a, t a, +t a, +t t a, b, +t a, )+ ( t a, t a, +t a, ) (7)
11 V. Hornak, R. Abel, A. Okur, B. Strockbine, A. Roitberg, and C. Simmerling, Comparison of multiple Amber force fields and development of improved protein backbone parameters, Proteins 65, 7 (6). B. Hess, C. Kutzner, D. van der Spoel, and E. Lindahl, Gromacs 4: Algorithms for highly efficient, load-balanced, and scalable molecular simulation, Journal of Chemical Theory and Computation 4, 435 (8). 3 W. L. Jorgensen, J. Chandrasekhar, J. D. Madura, R. W. Impey, and M. Klein, Comparison of simple potential functions for simulating liquid water, J. Chem. Phys. 79, 96 (983). 4 T. Darden, D. York, and L. Petersen, Particle mesh Ewald: An N log(n) method for Ewald sums in large systems, J. Chem. Phys. 98, 89 (993). 5 B. Hess, H. Bekker, H. J. C. Berendsen, and J. G. E. M. Fraaije, LINCS: A linear constraint solver for molecular simulations, J. Comp. Chem. 8, 463 (997). 6 H. J. C. Berendsen, J. P. M. Postma, W. F. van Gunsteren, A. Dinola, and J. R. Haak, Molecular dynamics with coupling to an external bath, J. Chem. Phys. 8, 3684 (984). 7 G. Bussi, D. Donadio, and M. Parrinello, Canonical sampling through velocity rescaling, J. Chem. Phys. 6, 4 (7). 8 N. Schaudinnus, A. J. Rzepiela, R. Hegger, and G. Stock, Data driven Langevin modeling of biomolecular dynamics, J. Chem. Phys. 38, 46 (3).
Molecular modeling. A fragment sequence of 24 residues encompassing the region of interest of WT-
SUPPLEMENTARY DATA Molecular dynamics Molecular modeling. A fragment sequence of 24 residues encompassing the region of interest of WT- KISS1R, i.e. the last intracellular domain (Figure S1a), has been
More informationRoutine access to millisecond timescale events with accelerated molecular dynamics
Routine access to millisecond timescale events with accelerated molecular dynamics Levi C.T. Pierce, Romelia Salomon-Ferrer, Cesar Augusto F. de Oliveira #, J. Andrew McCammon #, Ross C. Walker * SUPPORTING
More informationSupporting Information
Supporting Information Structure and Dynamics of Uranyl(VI) and Plutonyl(VI) Cations in Ionic Liquid/Water Mixtures via Molecular Dynamics Simulations Katie A. Maerzke, George S. Goff, Wolfgang H. Runde,
More informationMedical Research, Medicinal Chemistry, University of Leuven, Leuven, Belgium.
Supporting Information Towards peptide vaccines against Zika virus: Immunoinformatics combined with molecular dynamics simulations to predict antigenic epitopes of Zika viral proteins Muhammad Usman Mirza
More informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
More informationChoosing weights for simulated tempering
PHYSICAL REVIEW E 7, 173 7 Choosing weights for simulated tempering Sanghyun Park 1 and Vijay S. Pande 1, 1 Department of Chemistry, Stanford University, Stanford, California 935, USA Department of Structural
More informationSupporting Information for. Hydrogen Bonding Structure at Zwitterionic. Lipid/Water Interface
Supporting Information for Hydrogen Bonding Structure at Zwitterionic Lipid/Water Interface Tatsuya Ishiyama,, Daichi Terada, and Akihiro Morita,, Department of Applied Chemistry, Graduate School of Science
More informationMARTINI simulation details
S1 Appendix MARTINI simulation details MARTINI simulation initialization and equilibration In this section, we describe the initialization of simulations from Main Text section Residue-based coarsegrained
More informationAccurate Structure Prediction and Conformational Analysis of Cyclic Peptides with Residue-Specific Force Fields
Accurate Structure Prediction and Conformational Analysis of Cyclic Peptides with Residue-Specific Force Fields Hao Geng, Fan Jiang, *, and Yun-Dong Wu *,, School of Chemical Biology and Biotechnology,
More informationSupporting Information
Supporting Information ph-responsive self-assembly of polysaccharide through a rugged energy landscape Brian H. Morrow, Gregory F. Payne, and Jana Shen Department of Pharmaceutical Sciences, School of
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationSupporting Information: Improved Parametrization. of Lithium, Sodium, Potassium, and Magnesium ions. for All-Atom Molecular Dynamics Simulations of
Supporting Information: Improved Parametrization of Lithium, Sodium, Potassium, and Magnesium ions for All-Atom Molecular Dynamics Simulations of Nucleic Acid Systems. Jejoong Yoo and Aleksei Aksimentiev,,
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 informationSupporting information to: Time-resolved observation of protein allosteric communication. Sebastian Buchenberg, Florian Sittel and Gerhard Stock 1
Supporting information to: Time-resolved observation of protein allosteric communication Sebastian Buchenberg, Florian Sittel and Gerhard Stock Biomolecular Dynamics, Institute of Physics, Albert Ludwigs
More informationMolecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains
Molecular Dynamics Investigation of the ω-current in the Kv1.2 Voltage Sensor Domains Fatemeh Khalili-Araghi, Emad Tajkhorshid, Benoît Roux, and Klaus Schulten Department of Physics, Department of Biochemistry,
More informationImproved Resolution of Tertiary Structure Elasticity in Muscle Protein
Improved Resolution of Tertiary Structure Elasticity in Muscle Protein Jen Hsin and Klaus Schulten* Department of Physics and Beckman Institute, University of Illinois at Urbana-Champaign, Urbana, Illinois
More informationMolecular dynamics simulation of Aquaporin-1. 4 nm
Molecular dynamics simulation of Aquaporin-1 4 nm Molecular Dynamics Simulations Schrödinger equation i~@ t (r, R) =H (r, R) Born-Oppenheimer approximation H e e(r; R) =E e (R) e(r; R) Nucleic motion described
More informationJournal of Pharmacology and Experimental Therapy-JPET#172536
A NEW NON-PEPTIDIC INHIBITOR OF THE 14-3-3 DOCKING SITE INDUCES APOPTOTIC CELL DEATH IN CHRONIC MYELOID LEUKEMIA SENSITIVE OR RESISTANT TO IMATINIB Manuela Mancini, Valentina Corradi, Sara Petta, Enza
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationExercise 2: Solvating the Structure Before you continue, follow these steps: Setting up Periodic Boundary Conditions
Exercise 2: Solvating the Structure HyperChem lets you place a molecular system in a periodic box of water molecules to simulate behavior in aqueous solution, as in a biological system. In this exercise,
More informationScrutinizing Molecular Mechanics Force Fields on the Microsecond Timescale With NMR Data. Oliver Lange, David van der Spoel, and Bert L de Groot
Biophysical Journal, Volume 99 Supporting Material Scrutinizing Molecular Mechanics Force Fields on the Microsecond Timescale With NMR Data Oliver Lange, David van der Spoel, and Bert L de Groot Supplement:
More informationSupporting Information
pk a values of titrable amino acids at the water/membrane interface Vitor H. Teixeira, Diogo Vila-Viçosa, Pedro B. P. S. Reis, and Miguel Machuqueiro Centro de Química e Bioquímica, Departamento de Química
More informationarxiv: v1 [cond-mat.stat-mech] 28 Mar 2008
Canonical sampling through velocity-rescaling Giovanni Bussi, Davide Donadio, and Michele Parrinello Computational Science, Department of Chemistry and Applied Biosciences, ETH Zürich, USI Campus, Via
More informationThe Dominant Interaction Between Peptide and Urea is Electrostatic in Nature: A Molecular Dynamics Simulation Study
Dror Tobi 1 Ron Elber 1,2 Devarajan Thirumalai 3 1 Department of Biological Chemistry, The Hebrew University, Jerusalem 91904, Israel 2 Department of Computer Science, Cornell University, Ithaca, NY 14853
More informationWhat is Classical Molecular Dynamics?
What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential functions Newton s equations of motion are integrated
More informationA Nobel Prize for Molecular Dynamics and QM/MM What is Classical Molecular Dynamics? Simulation of explicit particles (atoms, ions,... ) Particles interact via relatively simple analytical potential
More informationUniversal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information
Universal Repulsive Contribution to the Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information B. Shadrack Jabes, Dusan Bratko, and Alenka Luzar Department of Chemistry,
More informationPeptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia
University of Groningen Peptide folding in non-aqueous environments investigated with molecular dynamics simulations Soto Becerra, Patricia IMPORTANT NOTE: You are advised to consult the publisher's version
More informationarxiv: v1 [physics.chem-ph] 30 Jul 2012
Is ion channel selectivity mediated by confined water? Diego Prada-Gracia and Francesco Rao School of Soft Matter, Freiburg Institute for Advanced Studies, Freiburg, Germany arxiv:27.6953v [physics.chem-ph]
More informationProbing α/β Balances in Modified Amber Force Fields from a Molecular Dynamics Study on a ββα Model Protein (1FSD)
Benchmark Study of Modified Amber Force Fields on ββα Protein Bull. Korean Chem. Soc. 214, Vol. 35, No. 6 1713 http://dx.doi.org/1.512/bkcs.214.35.6.1713 Probing α/β Balances in Modified Amber Force Fields
More informationHydrophobic Aided Replica Exchange: an Efficient Algorithm for Protein Folding in Explicit Solvent
19018 J. Phys. Chem. B 2006, 110, 19018-19022 Hydrophobic Aided Replica Exchange: an Efficient Algorithm for Protein Folding in Explicit Solvent Pu Liu, Xuhui Huang, Ruhong Zhou,, and B. J. Berne*,, Department
More informationarxiv: v2 [physics.chem-ph] 23 Aug 2013
The quest for self-consistency in hydrogen bond definitions Diego Prada-Gracia, Roman Shevchuk, and Francesco Rao Freiburg Institute for Advanced Studies, School of Soft Matter Research, Albertstrasse
More informationCrystal nucleation and growth both from melts and
pubs.acs.org/jpcl How Crystals Nucleate and Grow in Aqueous NaCl Solution Debashree Chakraborty and G. N. Patey* Department of Chemistry, University of British Columbia, Vancouver, British Columbia, Canada
More informationMolecular Dynamics. A very brief introduction
Molecular Dynamics A very brief introduction Sander Pronk Dept. of Theoretical Physics KTH Royal Institute of Technology & Science For Life Laboratory Stockholm, Sweden Why computer simulations? Two primary
More informationIntroduction to Classical Molecular Dynamics. Giovanni Chillemi HPC department, CINECA
Introduction to Classical Molecular Dynamics Giovanni Chillemi g.chillemi@cineca.it HPC department, CINECA MD ingredients Coordinates Velocities Force field Topology MD Trajectories Input parameters Analysis
More informationSupporting Information for: Physics Behind the Water Transport through. Nanoporous Graphene and Boron Nitride
Supporting Information for: Physics Behind the Water Transport through Nanoporous Graphene and Boron Nitride Ludovic Garnier, Anthony Szymczyk, Patrice Malfreyt, and Aziz Ghoufi, Institut de Physique de
More informationCarlo Camilloni, Andrea Cavalli,, and Michele Vendruscolo INTRODUCTION
pubs.acs.org/jpcb Assessment of the Use of NMR Chemical Shifts as Replica-Averaged Structural Restraints in Molecular Dynamics Simulations to Characterize the Dynamics of Proteins Carlo Camilloni, Andrea
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 informationSupporting Information. for. Influence of Cononsolvency on the. Aggregation of Tertiary Butyl Alcohol in. Methanol-Water Mixtures
Supporting Information for Influence of Cononsolvency on the Aggregation of Tertiary Butyl Alcohol in Methanol-Water Mixtures Kenji Mochizuki,, Shannon R. Pattenaude, and Dor Ben-Amotz Research Institute
More informationAnalysis of the simulation
Analysis of the simulation Marcus Elstner and Tomáš Kubař January 7, 2014 Thermodynamic properties time averages of thermodynamic quantites correspond to ensemble averages (ergodic theorem) some quantities
More informationSupporting Information. for An atomistic modeling of F-actin mechanical responses and determination of. mechanical properties
Supporting Information for An atomistic modeling of F-actin mechanical responses and determination of mechanical properties Authors: Si Li 1, Jin Zhang 2, Chengyuan Wang 1, Perumal Nithiarasu 1 1. Zienkiewicz
More informationQuantum Interference and Selectivity through Biological Ion Channels
1 SUPPLEMENTARY INFORMATION Quantum Interference and Selectivity through Biological Ion Channels Vahid Salari 1,3,*, Hamidreza Naeij 2, Afshin Shafiee 2,3 1 Department of Physics, Isfahan University of
More informationSupporting information
Electronic Supplementary Material ESI for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2018 Supporting information Understanding three-body contributions to coarse-grained force
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationBiomolecular modeling II
2015, December 16 System boundary and the solvent Biomolecule in solution typical MD simulations molecular system in aqueous solution task make the system as small as possible (reduce cost) straightforward
More informationProtecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations
pubs.acs.org/jctc Protecting High Energy Barriers: A New Equation to Regulate Boost Energy in Accelerated Molecular Dynamics Simulations William Sinko,*,, Cesar Augusto F. de Oliveira,*,,, Levi C. T. Pierce,
More informationSupplemntary Infomation: The nanostructure of. a lithium glyme solvate ionic liquid at electrified. interfaces
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 207 Supplemntary Infomation: The nanostructure of a lithium glyme solvate ionic liquid
More informationAtomistic nature of NaCl nucleation at the solid-liquid interface
THE JOURNAL OF CHEMICAL PHYSICS 126, 044708 2007 Atomistic nature of NaCl nucleation at the solid-liquid interface Yong Yang Institute of Physics, Chinese Academy of Sciences, Box 603, Beijing 100080,
More informationSupporting Information for the Effects of Electrostatic Interaction and Chirality on the
Supporting Information for the Effects of Electrostatic Interaction and Chirality on the Friction Coefficient of Water Flow Inside Single-Walled Carbon Nanotubes and Boron Nitride Nanotubes Xingfei Wei
More informationSupporting information:
Supporting information: A simple method to measure protein side-chain mobility using NMR chemical shifts. Mark V. Berjanskii, David S. Wishart* Department of Computing Science, University of Alberta, Edmonton,
More informationarxiv: v1 [q-bio.bm] 9 Jul 2007
Exploring the Protein G Helix Free Energy Surface by Solute Tempering Metadynamics Carlo Camilloni, 1, 2 Davide Provasi, 1, 2 Guido Tiana, 1, 2, and Ricardo A. Broglia 1,2,3 1 Department of Physics, University
More informationFree energy calculations and the potential of mean force
Free energy calculations and the potential of mean force IMA Workshop on Classical and Quantum Approaches in Molecular Modeling Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
More informationSimulation of molecular systems by molecular dynamics
Simulation of molecular systems by molecular dynamics Yohann Moreau yohann.moreau@ujf-grenoble.fr November 26, 2015 Yohann Moreau (UJF) Molecular Dynamics, Label RFCT 2015 November 26, 2015 1 / 35 Introduction
More informationPotential Energy (hyper)surface
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = " d dx U(x) Conformation
More informationGeneration of topology files of a protein chain and simulations of a dipeptide
Chapter-36 Generation of topology files of a protein chain and simulations of a dipeptide After studying this chapter, you will be able to 1) Generate a topology file of a dipeptide molecule. 2) Generate
More informationBiomolecular modeling I
2015, December 15 Biomolecular simulation Elementary body atom Each atom x, y, z coordinates A protein is a set of coordinates. (Gromacs, A. P. Heiner) Usually one molecule/complex of interest (e.g. protein,
More informationSupporting Information
Supporting Information Comparative Study of Materials-Binding Peptide Interactions with Gold and Silver Surfaces and Nanostructures: A Thermodynamic Basis for Biological Selectivity of Inorganic Materials
More informationProceedings of the ASME 2009 International Mechanical Engineering Congress & Exposition
Proceedings of the ASME 9 International Mechanical Engineering Congress & Exposition IMECE9 November 3-9, Lake Buena Vista, Florida, USA Proceedings of the ASME International Mechanical Engineering Congress
More informationSignificantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field
Significantly Improved Protein Folding Thermodynamics Using a Dispersion-Corrected Water Model and a New Residue-Specific Force Field Hao-Nan Wu, Fan Jiang, *, and Yun-Dong Wu, Laboratory of Computational
More informationChemical Physics Letters
Chemical Physics Letters 542 (2012) 37 41 Contents lists available at SciVerse ScienceDirect Chemical Physics Letters journal homepage: www.elsevier.com/locate/cplett Thermal conductivity, shear viscosity
More informationarxiv: v1 [physics.chem-ph] 24 Apr 2018
arxiv:1804.08913v1 [physics.chem-ph] 24 Apr 2018 Fast-Forward Langevin Dynamics with Momentum Flips Mahdi Hijazi, 1 David M. Wilkins, 1, a) and Michele Ceriotti 1 Laboratory of Computational Science and
More informationForce fields, thermo- and barostats. Berk Hess
Force fields, thermo- and barostats Berk Hess What is a force field? A force field usually consists of three parts: a set of functional forms parameters for the functional forms that, usually, depend on
More informationOn the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study
On the calculation of solvation free energy from Kirkwood- Buff integrals: A large scale molecular dynamics study Wynand Dednam and André E. Botha Department of Physics, University of South Africa, P.O.
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 informationAvailable online at
Available online at www.sciencedirect.com Chemical Physics Letters 452 (2008) 78 83 www.elsevier.com/locate/cplett Geometry and excitation energy fluctuations of NMA in aqueous solution with CHARMM, AMBER,
More informationHow NMR Works Section II
How NMR Works ----- Section II by Shaoxiong Wu 02-20-2010 In the section I we have evaluated one spin system by using an one spin Hamiltonian act on a EigenFunction and then get the EigenValue, i.e. the
More informationMultiple time step Monte Carlo simulations: Application to charged systems with Ewald summation
JOURNAL OF CHEMICAL PHYSICS VOLUME 11, NUMBER 1 1 JULY 004 Multiple time step Monte Carlo simulations: Application to charged systems with Ewald summation Katarzyna Bernacki a) Department of Chemistry
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for Chemical Science. This journal is The Royal Society of Chemistry 2016 Electronic Supplementary Information Active Site Structure and Absorption Spectrum of Channelrhodopsin-2
More informationReplica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2)
pubs.acs.org/jpcb Replica Exchange with Solute Scaling: A More Efficient Version of Replica Exchange with Solute Tempering (REST2) Lingle Wang, Richard A. Friesner, and B. J. Berne* Department of Chemistry,
More informationThe Molecular Dynamics Method
H-bond energy (kcal/mol) - 4.0 The Molecular Dynamics Method Fibronectin III_1, a mechanical protein that glues cells together in wound healing and in preventing tumor metastasis 0 ATPase, a molecular
More informationExploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials
Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials Dr. V. Subramanian Chemical Laboratory, IPC Division CSIR-Central Leather Research
More informationWhy Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar.
Why Proteins Fold? (Parts of this presentation are based on work of Ashok Kolaskar) CS490B: Introduction to Bioinformatics Mar. 25, 2002 Molecular Dynamics: Introduction At physiological conditions, the
More informationMolecular Dynamics Simulation of Methanol-Water Mixture
Molecular Dynamics Simulation of Methanol-Water Mixture Palazzo Mancini, Mara Cantoni University of Urbino Carlo Bo Abstract In this study some properties of the methanol-water mixture such as diffusivity,
More informationAb initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
More informationMOLECULAR DYNAMICS SIMULATION OF THE STRUCTURE OF C6 ALKANES INTRODUCTION. A. V. Anikeenko, A. V. Kim, and N. N. Medvedev UDC 544.2: 544.
Journal of Structural Chemistry. Vol. 51, No. 6, pp. 1090-1096, 2010 Original Russian Text Copyright 2010 by A. V. Anikeenko, A. V. Kim, and N. N. Medvedev MOLECULAR DYNAMICS SIMULATION OF THE STRUCTURE
More informationForce Field for Water Based on Neural Network
Force Field for Water Based on Neural Network Hao Wang Department of Chemistry, Duke University, Durham, NC 27708, USA Weitao Yang* Department of Chemistry, Duke University, Durham, NC 27708, USA Department
More informationExploring the Free Energy Surface of Short Peptides by Using Metadynamics
John von Neumann Institute for Computing Exploring the Free Energy Surface of Short Peptides by Using Metadynamics C. Camilloni, A. De Simone published in From Computational Biophysics to Systems Biology
More informationMolecular Dynamics Studies of Human β-glucuronidase
American Journal of Applied Sciences 7 (6): 83-88, 010 ISSN 1546-939 010Science Publications Molecular Dynamics Studies of Human β-lucuronidase Ibrahim Ali Noorbatcha, Ayesha Masrur Khan and Hamzah Mohd
More informationSupplementary Information. Surface Microstructure Engenders Unusual Hydrophobicity in. Phyllosilicates
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Supplementary Information Surface Microstructure Engenders Unusual Hydrophobicity in Phyllosilicates
More informationGromacs Workshop Spring CSC
Gromacs Workshop Spring 2007 @ CSC Erik Lindahl Center for Biomembrane Research Stockholm University, Sweden David van der Spoel Dept. Cell & Molecular Biology Uppsala University, Sweden Berk Hess Max-Planck-Institut
More information6 Similarities between principal components of protein dynamics and random diffusion
6 Similarities between principal components of protein dynamics and random diffusion B. Hess Phys. Rev. E 6:8438-8448 () Principal component analysis, also called essential dynamics, is a powerful tool
More informationWater Dynamics at Interfaces and Solutes: Disentangling Free Energy and Diffusivity Contributions
J Stat Phys (2011) 145:240 252 DOI 10.1007/s10955-011-0338-0 Water Dynamics at Interfaces and Solutes: Disentangling Free Energy and Diffusivity Contributions Felix Sedlmeier Yann von Hansen Liang Mengyu
More informationPrinciples and Applications of Molecular Dynamics Simulations with NAMD
Principles and Applications of Molecular Dynamics Simulations with NAMD Nov. 14, 2016 Computational Microscope NCSA supercomputer JC Gumbart Assistant Professor of Physics Georgia Institute of Technology
More informationDiffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study
Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Julia Deitz, Yeneneh Yimer, and Mesfin Tsige Department of Polymer Science University of Akron
More informationSupporting information for the article. Nuclear Spin Relaxation due to Chemical Shift Anisotropy of Gas-Phase 129 Xe
Supporting information for the article Nuclear Spin Relaxation due to Chemical Shift Anisotropy of Gas-Phase 9 Xe Matti Hanni, Perttu Lantto, and Juha Vaara NMR Research Group Department of Physics P.O.
More informationMolecular Dynamics in Principal Component Space
pubs.acs.org/jpcb Molecular Dynamics in Principal Component Space Servaas Michielssens,*, Titus S. van Erp, Carsten Kutzner, Arnout Ceulemans, and Bert L. de Groot Department of Chemistry and Centre for
More informationOn the accurate calculation of the dielectric constant and the diffusion coefficient from molecular dynamics simulations: the case of SPC/E water
On the accurate calculation of the dielectric constant and the diffusion coefficient from molecular dynamics simulations: the case of SPC/E water Orsolya Gereben and László Pusztai Research Institute for
More informationk θ (θ θ 0 ) 2 angles r i j r i j
1 Force fields 1.1 Introduction The term force field is slightly misleading, since it refers to the parameters of the potential used to calculate the forces (via gradient) in molecular dynamics simulations.
More informationSUPPLEMENTARY INFORMATION. Very Low Temperature Membrane Free Desalination. by Directional Solvent Extraction
SUPPLEMENTARY INFORMATION Very Low Temperature Membrane Free Desalination by Directional Solvent Extraction Anurag Bajpayee, Tengfei Luo, Andrew J. Muto, and Gang Chen Department of Mechanical Engineering,
More informationDihedral Angles. Homayoun Valafar. Department of Computer Science and Engineering, USC 02/03/10 CSCE 769
Dihedral Angles Homayoun Valafar Department of Computer Science and Engineering, USC The precise definition of a dihedral or torsion angle can be found in spatial geometry Angle between to planes Dihedral
More informationA MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES
A MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES Anastassia N. Rissanou b,c*, Vagelis Harmandaris a,b,c* a Department of Applied Mathematics, University of Crete, GR-79, Heraklion, Crete,
More informationDetermination of Kamlet-Taft parameters for selected solvate ionic liquids.
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Determination of Kamlet-Taft parameters for selected solvate ionic liquids. Daniel
More informationAb initio molecular dynamics and nuclear quantum effects
Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and
More informationSUPPLEMENTARY INFORMATION
DOI: 10.1038/NCHEM.1821 Cloud-based simulations on Google Exacycle reveal ligand-modulation of GPCR activation pathways Kai J. Kohlhoff 1,4 *, Diwakar Shukla 1,2, Morgan Lawrenz 2, Gregory R. Bowman 2,
More informationLecture 11: Potential Energy Functions
Lecture 11: Potential Energy Functions Dr. Ronald M. Levy ronlevy@temple.edu Originally contributed by Lauren Wickstrom (2011) Microscopic/Macroscopic Connection The connection between microscopic interactions
More informationOrganization of NAMD Tutorial Files
Organization of NAMD Tutorial Files .1.1. RMSD for individual residues Objective: Find the average RMSD over time of each residue in the protein using VMD. Display the protein with the residues colored
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 informationTailoring the Properties of Quadruplex Nucleobases for Biological and Nanomaterial Applications
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2014 Supporting Information for: Tailoring the Properties of Quadruplex Nucleobases
More informationEffects of bundling on the properties of the SPC water model
Theor Chem Acc (21) 125:335 344 DOI 1.17/s214-9-59-4 REGULAR ARTICLE Effects of bundling on the properties of the SPC water model Marc Fuhrmans Æ Barbara P. Sanders Æ Siewert-Jan Marrink Æ Alex H. de Vries
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More information