1 / 37

Introduction to Molecular Simulation

Introduction to Molecular Simulation. Chih-Hao Lu Associate Professor China Medical University. Movies. http://www.youtube.com/watch?v=Vj8ri57GE_M http://www.youtube.com/watch?v=lm-dAvbl330&feature=bf_prev&list=PL738AECAE0F368423&lf=results_main

rossman
Download Presentation

Introduction to Molecular Simulation

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Introduction to Molecular Simulation Chih-Hao Lu Associate Professor China Medical University

  2. Movies • http://www.youtube.com/watch?v=Vj8ri57GE_M • http://www.youtube.com/watch?v=lm-dAvbl330&feature=bf_prev&list=PL738AECAE0F368423&lf=results_main • http://www.youtube.com/watch?v=o5-a39tbT9w&feature=bf_prev&list=PL738AECAE0F368423&lf=results_main • http://www.youtube.com/watch?v=meNEUTn9Atg&feature=related • http://www.youtube.com/watch?v=kW2oajY-B6E&feature=related

  3. Molecular dynamics (MD) is a form of computer simulation in which atoms and molecules are allowed to interact for a period of time by approximations of known physics, giving a view of the motion of the particles.

  4. This kind of simulation is frequently used in the study of proteins and biomolecules. It is tempting, though not entirely accurate, to describe the technique as a "virtual microscope" with high temporal and spatial resolution.

  5. Richard Phillips Feynman (1918–1988) Richard Feynman once said that "If we were to name the most powerful assumption of all, which leads one on and on in an attempt to understand life, it is that all things are made of atoms, and that everything that living things do can be understood in terms of the jigglings and wigglings of atoms."

  6. Molecular dynamics lets scientists peer into the motion of individual atoms in a way which is not possible in laboratory experiments.

  7. John Desmond Bernal (1901-1971) in 1962: "... I took a number of rubber balls and stuck them together with rods of a selection of different lengths ranging from 2.75 to 4 inches. I tried to do this in the first place as casually as possible, working in my own office, being interrupted every five minutes or so and not remembering what I had done before the interruption."

  8. Because molecular systems generally consist of a vast number of particles, it is in general impossible to find the properties of such complex systems analytically.

  9. MD simulation circumvents the analytical intractability by using numerical methods.

  10. Babylonian clay tablet YBC 7289 (1800–1600 BC) ? Tip: Sexagesimal

  11. Babylonian clay tablet YBC 7289 (1800–1600 BC) 42,25,35 = 42/60 + 25/602 + 35/603 = 0.707106

  12. long MD simulations are mathematically ill-conditioned, generating cumulative errors in numerical integration that can be minimized with proper selection of algorithms and parameters, but not eliminated entirely.

  13. Current potential functions are, in many cases, not sufficiently accurate to reproduce the dynamics of molecular systems, so the much more computationally demanding ab Initio Molecular Dynamics method must be used.

  14. Design of a molecular dynamics simulation should account for the available computational power. Simulation size (n=number of particles), timestep and total time duration must be selected so that the calculation can finish within a reasonable time period.

  15. However, the simulations should be long enough to be relevant to the time scales of the natural processes being studied.

  16. Most scientific publications about the dynamics of proteins and DNA use data from simulations spanning nanoseconds (1E-9 s) to microseconds (1E-6 s).

  17. A molecular dynamics simulation requires the definition of a potential function, or a description of the terms by which the particles in the simulation will interact. In chemistry and biology this is usually referred to as a force field.

  18. Most force fields consist of a summation of bonded forces associated with chemical bonds, bond angles, and bond dihedrals, and non-bonded forces associated with van der Waals forces and electrostatic charge.

  19. These potentials contain free parameters such as atomic charge, van der Waals parameters reflecting estimates of atomic radius, and equilibrium bond length, angle, and dihedral

  20. In addition to the functional form of the potentials, a force field defines a set of parameters for each type of atom.

  21. For example, a force field would include distinct parameters for an oxygen atom in a carbonyl functional group and in a hydroxyl group.

  22. How many types of bonds? C-C 2 C-H 8 10 bonds ?? angle terms ?? torsional terms ?? non-bonded interactions

  23. How many types of angles? C-C-C 1 C-C-H 10 H-C-H 7 10 bonds 18 angle terms ?? torsional terms ?? non-bonded interactions

  24. How many types of torsions? H-C-C-H 12 H-C-C-C 6 10 bonds 18 angle terms 18 torsional terms ?? non-bonded interactions

  25. How many types of nonbonded ineractions? H-H 28 C-H 16C-C 1 10 bonds 18 angle terms 18 torsional terms 45 non-bonded interactions

  26. ? Bonds ? C-H ? C-C ? Angle terms ? C-C-C ? H-C-H ? H-C-C ? Torsional terms ? H-C-C-H ? H-C-C-C ? Non-bonded interactions ? C-H ? C-C ? H-H

  27. 13 Bonds 10 C-H 3 C-C 24 Angle terms 3 C-C-C 9 H-C-H 12 H-C-C 27 Torsional terms 9 H-C-C-H 18 H-C-C-C 78 Non-bonded interactions 30 C-H 3 C-C 45 H-H

  28. Thank you! Chih-Hao Lu chlu@mail.cmu.edu.tw

More Related