SasCalc Periodic Boundary Conditions Ian Hunt-Isaak, Joseph E. - - PowerPoint PPT Presentation

sascalc periodic boundary conditions
SMART_READER_LITE
LIVE PREVIEW

SasCalc Periodic Boundary Conditions Ian Hunt-Isaak, Joseph E. - - PowerPoint PPT Presentation

SasCalc Periodic Boundary Conditions Ian Hunt-Isaak, Joseph E. Curtis, Steven C. Howell Road Map 1. Question to answer 2. Issue with current approach 3. Alternative approaches 4. Results 5. Next Steps Why are concentrated proteins so


slide-1
SLIDE 1

SasCalc Periodic Boundary Conditions

Ian Hunt-Isaak, Joseph E. Curtis, Steven C. Howell

slide-2
SLIDE 2

Road Map

1. Question to answer 2. Issue with current approach 3. Alternative approaches 4. Results 5. Next Steps

slide-3
SLIDE 3

Why are concentrated proteins so viscous?

https://openclipart.org/image/2400px/svg_to_png/2058 86/remigho-syringe.png
slide-4
SLIDE 4

Small-Angle Scattering (SAS)

slide-5
SLIDE 5

Small-Angle Scattering (SAS)

Form Factor, P(Q) I(Q) ∝ P(Q) S(Q) Structure Factor, S(Q)

slide-6
SLIDE 6

SAS of dilute proteins used to determine shape

slide-7
SLIDE 7

Concentration Affects the Scattering

slide-8
SLIDE 8

Concentration Affects the Scattering

Form Factor Form Factor

slide-9
SLIDE 9

Co-solutes modulate protein interactions

Change in scattering is not dominated by P(q) or S(q) alone.

slide-10
SLIDE 10

Molecular Simulation

Adapted from http://nmschneider.com/blog/2014/5/9/nearest-image-for-particle-simulations-in-matlab
slide-11
SLIDE 11

Molecular Simulation

Adapted from http://nmschneider.com/blog/2014/5/9/nearest-image-for-particle-simulations-in-matlab
slide-12
SLIDE 12

Fourier Transform: g(r) → S(Q)

  • Most general
  • Slow
  • Limited to spherical systems
slide-13
SLIDE 13

Finite size effects

slide-14
SLIDE 14

Finite size effects

slide-15
SLIDE 15

Desired Features of Calculator

  • Simultaneously calculate both P(Q) and S(Q)
  • Avoid finite size effects
  • Fast
slide-16
SLIDE 16

Debye Formula

Not universally applicable: vs

slide-17
SLIDE 17

Explicit Fourier Transform

  • Most general

○ Contains P(Q) and S(Q)

slide-18
SLIDE 18

Removing finite size effects

  • Fancy Stat Mech - g(r)

○ Limited to spherical systems

  • Bigger Box

○ Num Atoms ~ (Box Length)3

  • Calculate Box scattering and remove
slide-19
SLIDE 19

First Test System - Lennard Jones Particles

Radial Distribution Function for Liquid Argon at 85 K - Yarnell, et al.

Verify Radial Distribution of Simulation

slide-20
SLIDE 20

Removing Box Effects

slide-21
SLIDE 21

Removing Box Effects

slide-22
SLIDE 22

Removing Box Effects

slide-23
SLIDE 23

Integrating g(r) vs Box Subtraction

slide-24
SLIDE 24

Next Steps

  • Apply this to a periodic box of proteins (lysozyme, mAb)
  • Automate algorithm to subtract box effects
  • Parallelize using GPUs
slide-25
SLIDE 25

Thank you!

slide-26
SLIDE 26

Backup Slides

slide-27
SLIDE 27

Existing Calculators

Experiment Simulation Watch out S(Q) ≠ I(Q).

slide-28
SLIDE 28

S(Q) via Fourier Transform

  • g(r) simple to

calculate

  • Can extend h(r)

via use of Ornstein-Zernike Equation

slide-29
SLIDE 29

SasCalc (Golden Vector)

  • Multiple Molecules?
  • Separation of S and P?

Two questions:

  • Can Extract S?
  • Periodic Boundary Conditions effects?
slide-30
SLIDE 30

Scattering of test system

slide-31
SLIDE 31

Current Limitations

10 mg/ml 100 mg/ml Of Medical interest? YES YES Can we calculate Scattering? YES NO (somewhat)

Simulation Boxes:

slide-32
SLIDE 32

Integrating g(r) vs Box Subtraction

slide-33
SLIDE 33

Integrating g(r) vs Box Subtraction