F. Piazza Center for Molecular Biophysics and University of Orléans, France. Selected topic in Physical Biology. Lectures 2-3
|
|
- Tyrone Sparks
- 5 years ago
- Views:
Transcription
1 F. Piazza Center for Molecular Biophysics and University of Orléans, France Selected topic in Physical Biology Lectures 2-3 Elastic network models of proteins. Theory, applications and much more than this
2 The force field The first thing: a force field
3 Structure-energy relation
4 Energy landscapes: local versus global minima 2D representation of a 3nD-dimensional landscape Energy minimization
5 Molecular dynamics simulations: stuck in a local minimum Computational costs at present only allow to explore local minima in general. For small system the ever-growing computational power allows to observe a few barrier-crossing events
6 Molecular dynamics simulations: integrate Newton s equation Nothing more than F = ma
7 Many force fields and full simulation packages are available You can install one of these packages on your local work-station and run your molecular dynamics simulation right away
8 A recent innovation: ANTON A microchip designed and optimized only to run molecular dynamics simulations
9 Coarse-grained models W. G. Noid, Perspective: Coarse-grained models for biomolecular systems The Journal of Chemical Physics 139, (2013) Given the system coordinates, a coordinate mapping, M, determines the configuration, R, of the CG model as a function of the configuration, r, of an underlying atomistic model. The Cartesian coordinates, R I, of site I are typically determined as a linear combination of atomic Cartesian coordinates, r i, with constant, positive coefficients that often correspond to, e.g., the center of mass or geometry for the associated atomic group.
10 A long story, starting in 1975 and recently culminated in a Nobel prize! Abstract A new and very simple representation of protein conformations has been used together with energy minimisation and thermalisation to simulate protein folding. Under certain conditions, the method succeeds in renaturing bovine pancreatic trypsin inhibitor from an open-chain conformation into a folded conformation close to that of the native molecule.
11 The Go model: the birth of native-centric modeling strategies Independently, in that same year Nobuhiro Go and his collaborators proposed a model an (even simpler) model where the chain of beads is mounted on on a lattice (initially they took a 2D lattice). Each bead would correspond to a residue or even to a secondary structure element of a protein (e.g. an -helix). The protein is restricted to fluctuate about the imposed structure: native-centric
12 Native-centric and off-lattice: the class of Elastic Network Models (ENMs) The story begins in 1996 Monique Tirion shows that the low frequency normal modes (NM) of a protein are not significantly altered when Interatomic interactions are replaced by identical Hookean springs (one-parameter model)!
13 Normal Mode Analysis (NMA) Typical all-atom potential energy V For small enough displacements about the equilibrium position:
14 Normal Mode Analysis (NMA) Displacements The equations of motion reads! (1)! where we have introduced the Hessian matrix! The system of coupled ODEs (1) can be transformed into a set of uncoupled ODEs by the Normal Mode coordinate transformation! Mass-weighted Hessian! Diagonalization! Coordinate transformation: the Normal Modes!
15 Normal Mode Analysis (NMA) Perform the change of variable and use the orthonormality and completeness relations! We get 3N uncoupled ODEs. It is the harmonic oscillator! The solutions are readily determined! So that the bead displacements can be written as!
16 Normal Mode Analysis (NMA) We note that the NM change of variables is such that the total potential energy and kinetic energy quadratic forms be diagonalized simultaneously. As a consequence, one has! 1 2 m i u 2 iα iα iα,jβ where we have introduced the normal mode energies! k = 1 2 Q 2 k + ω 2 kq 2 k u iα K αβ ij u jβ = k k At thermodynamic equilibrium! k = k B T H = 1 2 Ckω 2 k 2 = C k = k 2kB T ω k The amplitude of the k-th mode goes! as the inverse of its frequency. In general,! for proteins modes with frequencies! below cm -1 are responsible! for % of displacements.!
17 Elastic Network Models (ENMs)
18 Elastic Network Models (ENMs): the potential energy V A = 1 2 i>j k ij rij r 0 ij 2 Anisotropic Network Model (ANM) Nearly pure central forces V G = 1 2 = 1 2 i>j k ij r ij r 0 ij 2 k ij u j u i 2 Gaussian Network Model (ANM) Angular forces of the same order as central ones i>j i and j can be atoms (Tirion) or groups of atoms e.g. aggregated particles, e.g. amino acids r ij = r i r j r 0 ij = r 0 i r 0 j r ij r 0 ij = u j u i
19 Elastic Network Models (ENMs): the cutoff issue k ij = kf( r 0 ij ) is the force constant between particles i and j Popular choices for the force constants include: f( r 0 ij ) =c ij 1 for r 0 ij <R c 0 otherwise f( r 0 ij ) exp[ ( r 0 ij /σ) 2 ] f( r 0 ij ) r 0 ij α, α > 0 The sharp cutoff model R c vary between 8 and 16 Ang! Two-parameter models, in the sense that there is 1. One physical force scale gauging homogeneously interparticle force constants 2. One reference length specifying the range of interparticle interactions
20 The Gaussian network model is equivalent to a scalar model The GNM is intrinsically a scalar model!
21 Elastic Network Models (ENMs): the forces. The ANM scheme
22 Elastic Network Models (ENMs): the forces. The GNM ANM scheme
23 Elastic Network Models (ENMs): angular versus central forces Case (a) is meant to illustrate the magnitude of angle-bending forces, case (b) illustrates bond stretching. We want to compare the two kind of forces in the ANM and in the GNM schemes!
24 Elastic Network Models (ENMs): angular versus central forces ANM! GNM! Coupling along the direction of displacement! ANM! GNM! f (a) iy f (b) ix f (a) iy f (b) ix = O[(u/R) 2 ] = O[1]
25
26 Building an ENM An energy minimization is no longer required. This means that the equilibrium structure is assumed to coincide with the experimentally resolved structure (X-ray Crystallography, NMR). This can be risky! 1 First step: go to a structure repository and download a file containing the atomic coordinates of the macromolecules h"p://pdb.org/pdb/home/home.do The extension of such files is.pdb (for Protein Data Bank) 2 Use a computer program to read the coordinates from the PDB file, which has its own specific format. This will be your r 0 i vectors
27 PDB files contain a wealth of information on the protein x y z Atom number Atom type B factor Amino acid (three- leeer code) Chain number Chain
28 Calculating the mass-weighted Hessian is simple for ENMs! ANM H αβ ij = k ij k ij ˆRα mi m ij ˆRβ ij δ ij j m ˆR α im ˆR β im ˆR α ij r0 i,α r0 j,α r 0 ij Cartesian components ( =x,y,z) of the equilibrium unit vectors of inter-particle bonds GNM V G = 1 2 u T αku α α GNM H αβ ij = 1 2 ut (I 3 K)u = K ijδ αβ mi m j K = u is the 3N-dimensional vector of particle displacement k ij j=i K ij for i = j for i = j
29 Inter-particle correlation: the covariance matrix. The covariance matrix can be computed analytically!! Z u is the partition function! Let us introduce the matrix of eigenvectors of the mass-weighted Hessian matrix! Then we have!
30 Inter-particle correlation: the covariance matrix Let us perform the following change of variables! where!
31 Inter-particle correlation: the covariance matrix where!
32 Inter-particle correlation: the covariance matrix
33 The crystallographic B-factors: definition The crystallographic Debye-Waller factors (so-called isotropic B-factors) are related to atomic fluctuations. In fact, this is not as simple as that, as they are indeed refining parameters used in fitting the X-ray diffraction spectra, measuring line-widths. As such, they also contain (large!) contributions from roto-translations of the protein as a whole in the crystal and static disorder (the atoms of the protein in different crystal cells are not exactly in the very same position) In the GNM the problem is effectively N-dimensional and not 3N-dimensional B i = 8π 3 α = 8πk BT 3m i u 2 iα α (H 1 ) αα ii (H 1 ) αβ ij = m i m j K 1 ij δ αβ = B i =8πk B T K 1 ii
34 Atomic fluctuations in the GNM PnB- Esterase 13 It appears that the protein is almost rigid. With the excepion that the C- terminal region exhibit the largest correlated movement within the whole protein. This could be a linker or entrance mechanism to the acive site. InteresIngly, the residues around the proposed acive site (HIS 399, GLU 308 and SER 187) measured also high slow modes (in contrast of the whole protein)
35 Atomic fluctuations provide insight into biological function β- lactamase (1BLC) (a) and penicillopepsin (1BXO) (b), illustra7ng the mobility of residues in the first (lowest frequency) GNM mode. The color code is blue- red- yellow- green in the order of increasing mobility. Both enzymes contain an inhibitor (shown in space filling, gray) bound near the most constrained (lowest mobility) region. (c) and (d) Corresponding square fluctua7on profiles and posi7ons of cataly7c and inhibitor- binding residues. Residues directly involved in cataly1c func1on at ac1ve sites are shown by the green open circles, inhibitor- binding residues are shown by the red squares and residues serving both cataly1c and inhibitor- binding func1ons are marked by the orange diamond. Cataly1c residues tend to lie in the s1ffest por1ons of the protein structures ElasIc network models for understanding biomolecular machinery: from enzymes to supramolecular assemblies Chakra Chennubhotla, A J Rader, Lee- Wei Yang and Ivet Bahar Phys. Biol. 2 (2005) S173- S180
36 Catalytic Residues Coincide or communicate with Global Hinge Regions Fluctuation profiles in the global mode (k = 1) and position of catalytic and inhibitor binding residues illustrated for six enzymes. Residues involved in catalytic function are marked with an open circle, inhibitors binding sites are marked with a closed square, and residues serving both catalytic and inhibitor binding functions are marked with a closed circle. Arrows indicate the hinge sites. Lee-Wei Yang and Ivet Bahar, Coupling between Catalytic Site and Collective Dynamics: A Requirement for Mechanochemical Activity of Enzymes, Structure, 13, , June, 2005
37 With care: B factors do not only contain internal motions!
38 Conformational changes in proteins Many proteins exist in open (apo) and closed (holo, liganded) form. The biological function is closely related to the conformational change brought about by the apo-holo transition Open (top and middle) and closed (bottom) forms of lysine-arginine-ornithine (LAO) binding protein as shown usually (top) and modeled as an ENM coarse-grained at the C level with a cutoff R c = 8 Ang. F. Tama and Y.-H. Sanejouand, Conformational changes in proteins arising from normal mode calculations, Prot. Eng., 14, 1-6 (2001). One can ask the following question: How many NMs (a complete orthonormal basis in E 3N ) will be required to describe (reconstruct) the conformational change (a given percentage of it)? Will I require a number order N or order 1?
39 Conformational changes in proteins reconstructed through NMs: the overlap coefficients Let us consider our NMs as normalized If I use the whole basis of NMs, the (normalized) conformational change is reconstructed exactly ˆR iα Rholo iα Rapo 3N 6 iα R holo R apo = k=1 I k a k iα I k iα ˆR iα a k iα Overlap coefficients
40 The surprising answer There exists a single mode m for which I m is of the order !!
41 Normal modes of the open conformations perform better than those of the closed conformations In the open conformations domains are better separated, hence better defined. Therefore a coarse-grained description of the large-scale dynamics works better.
42 This is usually true for collective motions (low-frequency NMs) but There exist conformational changes that are captured by more localized modes I m = 0.3
43 I k iα ˆR iα a k iα Overlap coefficients Collectivity index The more collective the conformational change, the better the one-mode overlap c k = 1 3N iα (ak iα ak )( R iα R) σ(a k )σ( R) Correlation coefficients η = 1 3N exp jβ ˆR jβ 2 log ˆR jβ 2
44 Case study: predicting active sites in enzymes We have seen catalytic sites seem to have a tendency to lie at hinge-like regions in enzyme structures. This indication comes from the fact that these sites surprisingly often coincide with nodes of low-frequency normal modes.! Is it possible to devise specific indicators that help identify and predict active sites?! Simple tools from network theory: the connectivity graph! C i = j c ij Connectivity. How many neighbors at each node! CC i = j ij 1 Closeness centrality Inverse of the sum of the shortest paths from a given node to all other nodes.! χ i = k S hf ξ k i 2 Spectral stiffness Contribution of a reduced subset of! high-frequency NMs to the local fluctuation of a given node!
45 Predicting active sites in enzymes: high-pass filter Filtering procedure! We compute the indicator patterns but apply a high-pass filter, so as to only retain a reduced number of peaks! Arginin Glycineaminotransferase. PDB code 1JDW!
46 Predicting active sites in enzymes: the cutoff lensing idea Study the patterns of our structural indicators as functions of the cutoff used to build the connectivity graph. For the sake of the argument, we can also push it to values that may be thought of unphysical (excessively connected structures)! R c = 10 Å R c = 20 Å Cutoff lensing effect Some irrelevant peaks disappear and an additional peak appears flagging an active site.! Arginin Kinase. PDB code 1BG0!
47 Predicting active sites in enzymes: varying the cutoff The three indicators behave differently when the cutoff is let increase. In particular the connectivity pattern becomes less interesting at high cutoff as the structures become more and more connected.! We need to study the number of peaks as a function of the cutoff! Arginin Glycineaminotransferase. PDB code 1JDW!
48 Analysis of the enzyme database: the catalytic site atlas (CSA) n Fraction of catalytic residues within sites from! the nearest peak versus cutoff, as computed over the ensemble of enzymes from the CSA! Average peak fraction (number of peaks! divided by number of residues) computed over the CSA!
49 Reliability of the stiffness indicator The reliability is defined as the fraction of predicted catalytic sites (within namino acids along the sequence) divided by the fraction of stiffness peaks (number of peaks per amino acid).! Average number of peaks in the reduced! stiffness patterns per catalytic site!! The optimal cutoff corresponds to nearly! one peak per catalytic site.!! Extreme predictive precision!
50 Predicting active sites in enzymes: size matters! n fraction of catalytic sites within sites from the nearest peak of the reduced stiffness! patterns computed over three different size classes in the CSA database.!
51 The best predictions can be obtained by combining the 3 indicators at optimal cutoff in a sequential way The connectivity profiles should be examined first. These are the ones with the largest number! of peaks, often coalescing to highlight extended regions. The search should be subsequently! narrowed down with the corresponding closeness profile, typically featuring more localized peaks,! albeit many of them likely to be orphan ones. The prediction should then be refined through the! reduced stiffness patterns, the ones with the least number of peaks.!!"##$%&'()*+ R c = 20!"!!"##$%&'()*+ R c = 20!"! #$!%! ++!%! '(!%! )*!%! #'!%! '#!%! ))!%! )-!%! #&!%! )+!%! ''!%! &'!%! ''!%! ',!%!!/".$#$..+ R c = 28!"!!n!= 1!,&-#$..+ R c = 22!"!!/".$#$..+ R c = 28!"!!n!= 2!,&-#$..+ R c = 22!"!
52 Scientific Reports 5, Article number: (2015) doi: srep14874!
53 Normal modes strictly refer to very low temperature The thermal overlap coefficients
54 The thermal overlap coefficients
55 Aggregated spectral weight over a few modes is a good temperature-insensitive indicator There is redistribution of spectral weight at the working temperature. CAUTION IN USING T=0 normal modes!
56 The overall shape matters for spectral reconstructions of conformational fluctuations at non-zero temperature The less globular, the less cooperative, the worse 1G2F!
57 Proteins live immersed in a solvent. There is friction! (1) (2)
58 (3) (4) (5) (3)
59 Langevin modes is the block covariance matrix (6) (7) (8)
60 The block covariance matrix can be shown to obey the following solution (this is a straightforward consequence of imposing a Gaussian ansatz for the Fokker-Planck equation) (9) Eq. (9) has the following solution (see e.g. book by Risken, The Fokker-Planck equation ) (10) (11) (12)
61 (12) (12) (13) (10) (14) (11) (15) (15)
62 Intramolecular energy flux in proteins or within other complex three-dimensional molecular structures. General question!! 2-electrode setup: inject energy! at some site and monitor energy! outflux at a different site.!! 1. What are the energy transduction pathways?! 2. Are specific site pairs characterized by low impedance?!!!!
63 A typical non-equilibrium setup! Hot thermostat T 1! Generalized impedance! i j i j Cold thermostat T 2!
64 Energy flux We need first to introduce a measure of local energy flow in protein structures.! This comes naturally if we take the time derivative of the local energies!! Taking the time derivative!
65 Let us introduce the current from i to j, from i to j), then! (positive if energy flows from which leads to! Note that at equilibrium the following relations hold! Hence the total incoming and outgoing energy current is zero at each site!
66 In the harmonic approximation the expression for the energy current can be simplified further. One has! Hot thermostat T 1! i j Cold thermostat T 2! This expression can be used to measure the energy current between two given sites in a nonequilibrium setting such as the one we are interested into!
67 J i j ss = 1 2 αβ K αβ ij [ u iαu jβ ss u jβ u iα ss ] where (ss = steady-state) is an average computed with respect to the non-equilibrium steady-state measure! β 1 i Γ ij The imposed temperature field!
68 Practically, the sequence of operations is the following: 1. β 1 i Γ ij The imposed temperature field! 2. Compute the non-equilibrium covariance matrix! 3. Isolate the 3x3 blocks in the off-diagonal block that correspond to the two electrode beads and compute the current! J i j ss = 1 2 αβ K αβ ij [ u iαu jβ ss u jβ u iα ss ]
69 Questions!!?!
Protein Simulations in Confined Environments
Critical Review Lecture Protein Simulations in Confined Environments Murat Cetinkaya 1, Jorge Sofo 2, Melik C. Demirel 1 1. 2. College of Engineering, Pennsylvania State University, University Park, 16802,
More informationNORMAL MODE ANALYSIS AND GAUSSIAN NETWORK MODEL. Department of Mathematics Iowa State University, Ames, Iowa, 50010, USA
NORMAL MODE ANALYSIS AND GAUSSIAN NETWORK MODEL JUN KOO PARK Department of Mathematics Iowa State University, Ames, Iowa, 50010, USA Email:jkpark@iastate.edu Abstract. The functions of large biological
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 informationarxiv: v1 [q-bio.bm] 31 Jul 2017
Multiscale virtual particle based elastic network model (MVP-ENM) for biomolecular normal mode analysis arxiv:1708.07039v1 [q-bio.bm] 31 Jul 2017 Kelin Xia 1,2 1 Division of Mathematical Sciences, School
More informationUnderstanding protein motions by computational modeling and statistical approaches
Retrospective Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 28 Understanding protein motions by computational modeling and statistical approaches Lei Yang Iowa State
More informationIntroduction to Computational Structural Biology
Introduction to Computational Structural Biology Part I 1. Introduction The disciplinary character of Computational Structural Biology The mathematical background required and the topics covered Bibliography
More informationThe TEN Package. User Guide
================================================================= The TEN Package Tools for Elastic Networks User Guide Lars Ackermann, Silke A. Wieninger and G. Matthias Ullmann =================================================================
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 informationElastic network models for understanding biomolecular machinery: from enzymes to supramolecular assemblies
INSIUE OFPHYSICS PUBLISHING Phys. Biol. 2 (05) S173 S180 PHYSICAL BIOLOGY doi:10.1088/1478-3975/2/4/s12 Elastic network models for understanding biomolecular machinery: from enzymes to supramolecular assemblies
More informationSolid-state NMR and proteins : basic concepts (a pictorial introduction) Barth van Rossum,
Solid-state NMR and proteins : basic concepts (a pictorial introduction) Barth van Rossum, 16.02.2009 Solid-state and solution NMR spectroscopy have many things in common Several concepts have been/will
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 informationProtein dynamics. Folding/unfolding dynamics. Fluctuations near the folded state
Protein dynamics Folding/unfolding dynamics Passage over one or more energy barriers Transitions between infinitely many conformations B. Ozkan, K.A. Dill & I. Bahar, Protein Sci. 11, 1958-1970, 2002 Fluctuations
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 informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationProtein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall How do we go from an unfolded polypeptide chain to a
Lecture 11: Protein Folding & Stability Margaret A. Daugherty Fall 2004 How do we go from an unfolded polypeptide chain to a compact folded protein? (Folding of thioredoxin, F. Richards) Structure - Function
More informationNMR, X-ray Diffraction, Protein Structure, and RasMol
NMR, X-ray Diffraction, Protein Structure, and RasMol Introduction So far we have been mostly concerned with the proteins themselves. The techniques (NMR or X-ray diffraction) used to determine a structure
More informationPredicting the order in which contacts are broken during single molecule protein stretching experiments
Biochemistry, Biophysics and Molecular Biology Publications Biochemistry, Biophysics and Molecular Biology 4-2008 Predicting the order in which contacts are broken during single molecule protein stretching
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 informationF. Piazza Center for Molecular Biophysics and University of Orléans, France. Selected topic in Physical Biology. Lecture 1
Zhou Pei-Yuan Centre for Applied Mathematics, Tsinghua University November 2013 F. Piazza Center for Molecular Biophysics and University of Orléans, France Selected topic in Physical Biology Lecture 1
More informationComparisons of Experimental and Computed Protein Anisotropic Temperature Factors
Biochemistry, Biophysics and Molecular Biology Publications Biochemistry, Biophysics and Molecular Biology 7-2009 Comparisons of Experimental and Computed Protein Anisotropic Temperature Factors Lei Yang
More informationNormal Mode Analysis. Jun Wan and Willy Wriggers School of Health Information Sciences & Institute of Molecular Medicine University of Texas Houston
Normal Mode Analysis Jun Wan and Willy Wriggers School of Health Information Sciences & Institute of Molecular Medicine University of Teas Houston Contents Introduction --- protein dynamics (lowfrequency
More informationIsothermal experiments characterize time-dependent aggregation and unfolding
1 Energy Isothermal experiments characterize time-dependent aggregation and unfolding Technical ote Introduction Kinetic measurements have, for decades, given protein scientists insight into the mechanisms
More informationIntroduction The gramicidin A (ga) channel forms by head-to-head association of two monomers at their amino termini, one from each bilayer leaflet. Th
Abstract When conductive, gramicidin monomers are linked by six hydrogen bonds. To understand the details of dissociation and how the channel transits from a state with 6H bonds to ones with 4H bonds or
More informationDispersion relation for transverse waves in a linear chain of particles
Dispersion relation for transverse waves in a linear chain of particles V. I. Repchenkov* It is difficult to overestimate the importance that have for the development of science the simplest physical and
More informationLecture 11: Protein Folding & Stability
Structure - Function Protein Folding: What we know Lecture 11: Protein Folding & Stability 1). Amino acid sequence dictates structure. 2). The native structure represents the lowest energy state for a
More informationProtein Folding & Stability. Lecture 11: Margaret A. Daugherty. Fall Protein Folding: What we know. Protein Folding
Lecture 11: Protein Folding & Stability Margaret A. Daugherty Fall 2003 Structure - Function Protein Folding: What we know 1). Amino acid sequence dictates structure. 2). The native structure represents
More informationHands-on Course in Computational Structural Biology and Molecular Simulation BIOP590C/MCB590C. Course Details
Hands-on Course in Computational Structural Biology and Molecular Simulation BIOP590C/MCB590C Emad Tajkhorshid Center for Computational Biology and Biophysics Email: emad@life.uiuc.edu or tajkhors@uiuc.edu
More informationX-ray Crystallography
2009/11/25 [ 1 ] X-ray Crystallography Andrew Torda, wintersemester 2009 / 2010 X-ray numerically most important more than 4/5 structures Goal a set of x, y, z coordinates different properties to NMR History
More informationVisualization of Macromolecular Structures
Visualization of Macromolecular Structures Present by: Qihang Li orig. author: O Donoghue, et al. Structural biology is rapidly accumulating a wealth of detailed information. Over 60,000 high-resolution
More informationMolecular dynamics simulations of anti-aggregation effect of ibuprofen. Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov
Biophysical Journal, Volume 98 Supporting Material Molecular dynamics simulations of anti-aggregation effect of ibuprofen Wenling E. Chang, Takako Takeda, E. Prabhu Raman, and Dmitri Klimov Supplemental
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 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 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 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 informationVIII.B Equilibrium Dynamics of a Field
VIII.B Equilibrium Dynamics of a Field The next step is to generalize the Langevin formalism to a collection of degrees of freedom, most conveniently described by a continuous field. Let us consider the
More informationQuantitative Stability/Flexibility Relationships; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 12
Quantitative Stability/Flexibility Relationships; Donald J. Jacobs, University of North Carolina at Charlotte Page 1 of 12 The figure shows that the DCM when applied to the helix-coil transition, and solved
More informationMolecular Modeling lecture 2
Molecular Modeling 2018 -- lecture 2 Topics 1. Secondary structure 3. Sequence similarity and homology 2. Secondary structure prediction 4. Where do protein structures come from? X-ray crystallography
More informationProtein crystallography. Garry Taylor
Protein crystallography Garry Taylor X-ray Crystallography - the Basics Grow crystals Collect X-ray data Determine phases Calculate ρ-map Interpret map Refine coordinates Do the biology. Nitrogen at -180
More informationComputational Structural Biology and Molecular Simulation. Introduction to VMD Molecular Visualization and Analysis
Computational Structural Biology and Molecular Simulation Introduction to VMD Molecular Visualization and Analysis Emad Tajkhorshid Department of Biochemistry, Beckman Institute, Center for Computational
More informationAdvantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations
The Hilltop Review Volume 7 Issue 1 Winter 2014 Article 10 December 2014 Advantages of a Finite Extensible Nonlinear Elastic Potential in Lattice Boltzmann Simulations Tai-Hsien Wu Western Michigan University
More informationPrediction and refinement of NMR structures from sparse experimental data
Prediction and refinement of NMR structures from sparse experimental data Jeff Skolnick Director Center for the Study of Systems Biology School of Biology Georgia Institute of Technology Overview of talk
More informationNormal mode analysis is a powerful tool for describing the
A minimalist network model for coarse-grained normal mode analysis and its application to biomolecular x-ray crystallography Mingyang Lu* and Jianpeng Ma* *Department of Biochemistry and Molecular Biology,
More informationAlpha-helical Topology and Tertiary Structure Prediction of Globular Proteins Scott R. McAllister Christodoulos A. Floudas Princeton University
Alpha-helical Topology and Tertiary Structure Prediction of Globular Proteins Scott R. McAllister Christodoulos A. Floudas Princeton University Department of Chemical Engineering Program of Applied and
More informationPHASE TRANSITIONS IN SOFT MATTER SYSTEMS
OUTLINE: Topic D. PHASE TRANSITIONS IN SOFT MATTER SYSTEMS Definition of a phase Classification of phase transitions Thermodynamics of mixing (gases, polymers, etc.) Mean-field approaches in the spirit
More information3.091 Introduction to Solid State Chemistry. Lecture Notes No. 5a ELASTIC BEHAVIOR OF SOLIDS
3.091 Introduction to Solid State Chemistry Lecture Notes No. 5a ELASTIC BEHAVIOR OF SOLIDS 1. INTRODUCTION Crystals are held together by interatomic or intermolecular bonds. The bonds can be covalent,
More informationIntroduction to" Protein Structure
Introduction to" Protein Structure Function, evolution & experimental methods Thomas Blicher, Center for Biological Sequence Analysis Learning Objectives Outline the basic levels of protein structure.
More informationStructuralBiology. October 18, 2018
StructuralBiology October 18, 2018 1 Lecture 18: Computational Structural Biology CBIO (CSCI) 4835/6835: Introduction to Computational Biology 1.1 Overview and Objectives Even though we re officially moving
More informationMolecular Biophysics III - DYNAMICS
References Molecular Biophysics III - DYNAMICS Molecular Biophysics Structures in Motion, Michel Daune, Oxford U Press, 1999 (reprinted 2003, 2004). Computational Molecular Dynamics: Challenges, Methods,
More informationA Unification of the Elastic Network Model and the Gaussian Network Model for Optimal Description of Protein Conformational Motions and Fluctuations
Biophysical Journal Volume 94 May 2008 3853 3857 3853 A Unification of the Elastic Network Model and the Gaussian Network Model for Optimal Description of Protein Conformational Motions and Fluctuations
More informationTHE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION
THE TANGO ALGORITHM: SECONDARY STRUCTURE PROPENSITIES, STATISTICAL MECHANICS APPROXIMATION AND CALIBRATION Calculation of turn and beta intrinsic propensities. A statistical analysis of a protein structure
More informationEnergy functions and their relationship to molecular conformation. CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror
Energy functions and their relationship to molecular conformation CS/CME/BioE/Biophys/BMI 279 Oct. 3 and 5, 2017 Ron Dror Yesterday s Nobel Prize: single-particle cryoelectron microscopy 2 Outline Energy
More informationDynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics
for resubmission Dynamic force matching: Construction of dynamic coarse-grained models with realistic short time dynamics and accurate long time dynamics Aram Davtyan, 1 Gregory A. Voth, 1 2, a) and Hans
More informationResidue Network Construction and Predictions of Elastic Network Models
Residue Network Construction and Predictions of Elastic Network Models Canan Atilgan, Ibrahim Inanc, and Ali Rana Atilgan Faculty of Engineering and Natural Sciences, Sabanci University, 34956 Istanbul,
More informationInvestigating the effects of different force fields on spring-based normal mode analysis
Graduate Theses and Dissertations Iowa State University Capstones, Theses and Dissertations 2016 Investigating the effects of different force fields on spring-based normal mode analysis Jaekyun Song Iowa
More informationMolecular Mechanics. I. Quantum mechanical treatment of molecular systems
Molecular Mechanics I. Quantum mechanical treatment of molecular systems The first principle approach for describing the properties of molecules, including proteins, involves quantum mechanics. For example,
More informationExamples of Protein Modeling. Protein Modeling. Primary Structure. Protein Structure Description. Protein Sequence Sources. Importing Sequences to MOE
Examples of Protein Modeling Protein Modeling Visualization Examination of an experimental structure to gain insight about a research question Dynamics To examine the dynamics of protein structures To
More informationGeometrical Concept-reduction in conformational space.and his Φ-ψ Map. G. N. Ramachandran
Geometrical Concept-reduction in conformational space.and his Φ-ψ Map G. N. Ramachandran Communication paths in trna-synthetase: Insights from protein structure networks and MD simulations Saraswathi Vishveshwara
More informationOn the Conservation of the Slow Conformational Dynamics within the Amino Acid Kinase Family: NAGK the Paradigm
On the Conservation of the Slow Conformational Dynamics within the Amino Acid Kinase Family: NAGK the Paradigm Enrique Marcos 1, Ramon Crehuet 1 *, Ivet Bahar 2 * 1 Department of Biological Chemistry and
More informationProtein Dynamics, Allostery and Function
Protein Dynamics, Allostery and Function Lecture 2. Protein Dynamics Xiaolin Cheng UT/ORNL Center for Molecular Biophysics SJTU Summer School 2017 1 Functional Protein Dynamics Proteins are dynamic and
More informationOverview - Macromolecular Crystallography
Overview - Macromolecular Crystallography 1. Overexpression and crystallization 2. Crystal characterization and data collection 3. The diffraction experiment 4. Phase problem 1. MIR (Multiple Isomorphous
More informationAA 242B / ME 242B: Mechanical Vibrations (Spring 2016)
AA 242B / ME 242B: Mechanical Vibrations (Spring 206) Solution of Homework #3 Control Tab Figure : Schematic for the control tab. Inadequacy of a static-test A static-test for measuring θ would ideally
More informationProtein Structure Determination from Pseudocontact Shifts Using ROSETTA
Supporting Information Protein Structure Determination from Pseudocontact Shifts Using ROSETTA Christophe Schmitz, Robert Vernon, Gottfried Otting, David Baker and Thomas Huber Table S0. Biological Magnetic
More informationAssignment 2 Atomic-Level Molecular Modeling
Assignment 2 Atomic-Level Molecular Modeling CS/BIOE/CME/BIOPHYS/BIOMEDIN 279 Due: November 3, 2016 at 3:00 PM The goal of this assignment is to understand the biological and computational aspects of macromolecular
More informationATP Binding Enables Broad Antibiotic Selectivity of Aminoglycoside Phosphotransferase(3 )-IIIa: An Elastic Network Analysis
doi:10.1016/j.jmb.2011.03.061 J. Mol. Biol. (2011) 409, 450 465 Contents lists available at www.sciencedirect.com Journal of Molecular Biology journal homepage: http://ees.elsevier.com.jmb ATP Binding
More informationStructural characterization of NiV N 0 P in solution and in crystal.
Supplementary Figure 1 Structural characterization of NiV N 0 P in solution and in crystal. (a) SAXS analysis of the N 32-383 0 -P 50 complex. The Guinier plot for complex concentrations of 0.55, 1.1,
More informationSupplementary Figures
1 Supplementary Figures Supplementary Figure 1 Type I FGFR1 inhibitors (a) Chemical structures of a pyrazolylaminopyrimidine inhibitor (henceforth referred to as PAPI; PDB-code of the FGFR1-PAPI complex:
More informationPROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS
TASKQUARTERLYvol.20,No4,2016,pp.353 360 PROTEIN-PROTEIN DOCKING REFINEMENT USING RESTRAINT MOLECULAR DYNAMICS SIMULATIONS MARTIN ZACHARIAS Physics Department T38, Technical University of Munich James-Franck-Str.
More information1) NMR is a method of chemical analysis. (Who uses NMR in this way?) 2) NMR is used as a method for medical imaging. (called MRI )
Uses of NMR: 1) NMR is a method of chemical analysis. (Who uses NMR in this way?) 2) NMR is used as a method for medical imaging. (called MRI ) 3) NMR is used as a method for determining of protein, DNA,
More informationMany proteins spontaneously refold into native form in vitro with high fidelity and high speed.
Macromolecular Processes 20. Protein Folding Composed of 50 500 amino acids linked in 1D sequence by the polypeptide backbone The amino acid physical and chemical properties of the 20 amino acids dictate
More informationTime-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics
Time-Dependent Statistical Mechanics 5. The classical atomic fluid, classical mechanics, and classical equilibrium statistical mechanics c Hans C. Andersen October 1, 2009 While we know that in principle
More informationCONFORMATIONAL DYNAMICS OF PROTEINS: INSIGHTS FROM STRUCTURAL AND COMPUTATIONAL STUDIES
CONFORMATIONAL DYNAMICS OF PROTEINS: INSIGHTS FROM STRUCTURAL AND COMPUTATIONAL STUDIES By Lin Liu BS in Biological Science, University of Science and Technology of China, 2004 Submitted to the Graduate
More informationProteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability
Proteins are not rigid structures: Protein dynamics, conformational variability, and thermodynamic stability Dr. Andrew Lee UNC School of Pharmacy (Div. Chemical Biology and Medicinal Chemistry) UNC Med
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 informationNitrogenase MoFe protein from Clostridium pasteurianum at 1.08 Å resolution: comparison with the Azotobacter vinelandii MoFe protein
Acta Cryst. (2015). D71, 274-282, doi:10.1107/s1399004714025243 Supporting information Volume 71 (2015) Supporting information for article: Nitrogenase MoFe protein from Clostridium pasteurianum at 1.08
More informationFondamenti di Chimica Farmaceutica. Computer Chemistry in Drug Research: Introduction
Fondamenti di Chimica Farmaceutica Computer Chemistry in Drug Research: Introduction Introduction Introduction Introduction Computer Chemistry in Drug Design Drug Discovery: Target identification Lead
More informationBetter Bond Angles in the Protein Data Bank
Better Bond Angles in the Protein Data Bank C.J. Robinson and D.B. Skillicorn School of Computing Queen s University {robinson,skill}@cs.queensu.ca Abstract The Protein Data Bank (PDB) contains, at least
More informationFigure 1. Molecules geometries of 5021 and Each neutral group in CHARMM topology was grouped in dash circle.
Project I Chemistry 8021, Spring 2005/2/23 This document was turned in by a student as a homework paper. 1. Methods First, the cartesian coordinates of 5021 and 8021 molecules (Fig. 1) are generated, in
More informationApplication of Finite Element Method in Protein Normal Mode Analysis
University of Massachusetts Amherst ScholarWorks@UMass Amherst Masters Theses 1911 - February 2014 2013 Application of Finite Element Method in Protein Normal Mode Analysis Chiung-fang Hsu University of
More informationing equilibrium i Dynamics? simulations on AchE and Implications for Edwin Kamau Protein Science (2008). 17: /29/08
Induced-fit d or Pre-existin ing equilibrium i Dynamics? Lessons from Protein Crystallography yand MD simulations on AchE and Implications for Structure-based Drug De esign Xu Y. et al. Protein Science
More informationI690/B680 Structural Bioinformatics Spring Protein Structure Determination by NMR Spectroscopy
I690/B680 Structural Bioinformatics Spring 2006 Protein Structure Determination by NMR Spectroscopy Suggested Reading (1) Van Holde, Johnson, Ho. Principles of Physical Biochemistry, 2 nd Ed., Prentice
More informationSupplementary Material
1 2 3 Topological defects in confined populations of spindle-shaped cells by G. Duclos et al. Supplementary Material 4 5 6 7 8 9 10 11 12 13 Supplementary Note 1: Characteristic time associated with the
More informationMolecular Biology Course 2006 Protein Crystallography Part II
Molecular Biology Course 2006 Protein Crystallography Part II Tim Grüne University of Göttingen Dept. of Structural Chemistry December 2006 http://shelx.uni-ac.gwdg.de tg@shelx.uni-ac.gwdg.de Overview
More informationBIOMOLECULAR DYNAMICS REVEALED BY ELASTIC NETWORK MODELS AND THE STUDY OF MECHANICAL KEY SITES FOR LIGAND BINDING. Lee-Wei Yang
BIOMOLECULAR DYNAMICS REVEALED BY ELASTIC NETWORK MODELS AND THE STUDY OF MECHANICAL KEY SITES FOR LIGAND BINDING by Lee-Wei Yang BS, National Taiwan University, 1997 MS, National Tsing Hua University,
More informationBioengineering 215. An Introduction to Molecular Dynamics for Biomolecules
Bioengineering 215 An Introduction to Molecular Dynamics for Biomolecules David Parker May 18, 2007 ntroduction A principal tool to study biological molecules is molecular dynamics simulations (MD). MD
More informationPersistent homology analysis of protein structure, flexibility and folding
Persistent homology analysis of protein structure, flexibility and folding arxiv:1412.2779v1 [q-bio.bm] 8 Dec 2014 Kelin Xia 1,2 Guo-Wei Wei 1,2,3,4 1 Department of Mathematics Michigan State University,
More informationSecondary structure stability, beta-sheet formation & stability
Protein Physics 2016 Lecture 6, February 5 Secondary structure stability, beta-sheet formation & stability Magnus Andersson magnus.andersson@scilifelab.se Theoretical & Computational Biophysics Recap of
More informationInterpreting and evaluating biological NMR in the literature. Worksheet 1
Interpreting and evaluating biological NMR in the literature Worksheet 1 1D NMR spectra Application of RF pulses of specified lengths and frequencies can make certain nuclei detectable We can selectively
More informationChapter 2. Dielectric Theories
Chapter Dielectric Theories . Dielectric Theories 1.1. Introduction Measurements of dielectric properties of materials is very important because it provide vital information regarding the material characteristics,
More informationarxiv: v1 [q-bio.bm] 26 Oct 2015
Multiscale Gaussian network model (mgnm) and multiscale anisotropic network model (manm) arxiv:1510.07699v1 [q-bio.bm] 26 Oct 2015 Kelin Xia 1, Kristopher Opron 2 and Guo-Wei Wei 3 1 Department of Mathematics
More informationThe Gaussian Network Model: Theory and Applications
3 The Gaussian Network Model: Theory and Applications A.J. Rader, Chakra Chennubhotla, Lee-Wei Yang, and Ivet Bahar CONTENTS 3.1 Introduction... 41 3.1.1 Conformational Dynamics: A Bridge Between Structure
More informationDesign of a Novel Globular Protein Fold with Atomic-Level Accuracy
Design of a Novel Globular Protein Fold with Atomic-Level Accuracy Brian Kuhlman, Gautam Dantas, Gregory C. Ireton, Gabriele Varani, Barry L. Stoddard, David Baker Presented by Kate Stafford 4 May 05 Protein
More informationAn introduction to Molecular Dynamics. EMBO, June 2016
An introduction to Molecular Dynamics EMBO, June 2016 What is MD? everything that living things do can be understood in terms of the jiggling and wiggling of atoms. The Feynman Lectures in Physics vol.
More informationProtein Dynamics. The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron.
Protein Dynamics The space-filling structures of myoglobin and hemoglobin show that there are no pathways for O 2 to reach the heme iron. Below is myoglobin hydrated with 350 water molecules. Only a small
More informationA. Reaction Mechanisms and Catalysis (1) proximity effect (2) acid-base catalysts (3) electrostatic (4) functional groups (5) structural flexibility
(P&S Ch 5; Fer Ch 2, 9; Palm Ch 10,11; Zub Ch 9) A. Reaction Mechanisms and Catalysis (1) proximity effect (2) acid-base catalysts (3) electrostatic (4) functional groups (5) structural flexibility B.
More informationAb-initio protein structure prediction
Ab-initio protein structure prediction Jaroslaw Pillardy Computational Biology Service Unit Cornell Theory Center, Cornell University Ithaca, NY USA Methods for predicting protein structure 1. Homology
More informationBridging Structure and Function, Experiments and Computations. Ivet Bahar
Day 2 Bridging Structure and Function, Experiments and Computations Ivet Bahar Department of 1 Computational and Systems Biology School of Medicine, University of Pittsburgh, Pittsburgh, PA 15260 Summary
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 informationLecture 27. Transition States and Enzyme Catalysis
Lecture 27 Transition States and Enzyme Catalysis Reading for Today: Chapter 15 sections B and C Chapter 16 next two lectures 4/8/16 1 Pop Question 9 Binding data for your thesis protein (YTP), binding
More informationHelpful resources for all X ray lectures Crystallization http://www.hamptonresearch.com under tech support: crystal growth 101 literature Spacegroup tables http://img.chem.ucl.ac.uk/sgp/mainmenu.htm Crystallography
More informationCoarse-Grained Models Reveal Functional Dynamics - I. Elastic Network Models Theories, Comparisons and Perspectives
REVIEW Coarse-Grained Models Reveal Functional Dynamics - I. Elastic Network Models Theories, Comparisons and Perspectives Lee-Wei Yang and Choon-Peng Chng Institute of Molecular and Cellular Biosciences,
More information