MD MM
MD MM
MD MM
simulation
Cláudio M. Soares
Protein Modelling Laboratory
ITQB NOVA
Instituto de Tecnologia Química e Biológica António Xavier
Universidade Nova de Lisboa
Oeiras
ITQB NOVA
Instituto de Tecnologia Química e Biológica António Xavier
Universidade Nova de Lisboa
Av. da República,
2780-157 Oeiras
PORTUGAL
Tel: (351-1) 21446 9259 / 21446 9610
e-mail: claudio@itqb.unl.pt
http://www.itqb.unl.pt/pm URL: http://www.itqb.unl.pt/~claudio
Physical Biology
+ and
computational methods
Biotechnology
The aim is to understand the physical world of atoms and molecules at
the microscopic level.
C.M. Soares – Physical Modelling of Biomolecules – 2023
Why Protein Modelling ?
• Structure-Function relationships
Interpretation of experimental results (structural for
example) using physical and empirical models
• Predict (new) molecular properties
Predict structures and (altered) functions of mutant
proteins
• Study properties not accessible to experiments
Protein dynamics in the picosecond time scale
Discrete Continuum
out
in
HF/3-21G
1 ps (10-12 s)
Note: these functions cannot be derived from classical physics, only the
integration of equations of motion is done within this framework
Nb
V( r1 ,.....,rNat ) = 1 2
K ( b − b0 ) +
2 bn n n
n =1
N
1 2
K ( − 0 ) +
n =1 2 n n n
N
1 2
K ( − 0 ) +
n =1 2 n n n
N
1
K 1+ cos( mnn − n +
n =1 2 n
N at
C (i, j ) C (i, j )
qi q j
12 12 − 6 6 + S( rij ) +
i j rij rij 4 0 r rij
Sp ecial terms
C.M. Soares – Physical Modelling of Biomolecules – 2023
Covalent bond potential
Harmonic potential
Potential energy
1 2
K b (bn − b0 ) “Experimental”
bond potential
2 n n
b0n=Equilibrium distance
Interatomic distance
CCSD/cc-pVQZ
CO Carbon monoxide
C.M. Soares – Physical Modelling of Biomolecules – 2023
Bond angle potential
1 2
2
K n ( n − 0 )
n 0n=Equilibrium angle
1 2
Energy
K ( n − 0 ) 0n=Equilibrium angle
2 n n
Planar centers
C
Chiral centers
B D
E
A C A
B
Keeping the angle between E-D in the A-B-E plane Keeping the angle between C-D and the A-B-C plane
C.M. Soares – Physical Modelling of Biomolecules – 2023
Dihedral angle potential
1
𝐾𝜑𝑛 1 + cos( 𝑚𝑛 𝜑𝑛 − 𝛿𝑛 )
Energy
2
Ethane
A D
Three maxima and three minima
C
Triple degenerated system B
12 − 6
equation
•Prevent collision
C6 - Attractive
•Dispersion forces Distance
Induced dipole-Induced dipole
C.M. Soares – Physical Modelling of Biomolecules – 2023
Non-bonded interactions
II.The electrostatic potential - +
qq
i j + +
4 0 r rij
Force field parameters are specific for each atom type (not element) and
for each type of interaction
Many uses
•Monte Carlo
>Metropolis Monte Carlo for instance
•Others
C.M. Soares – Physical Modelling of Biomolecules – 2023
Energy minimization
Starting from a conformation, obtain
conformations with lower energy
for the molecule
•Move atoms
Search conformational space Start
Energy
•Find conformations with lower
energy
Minima in the potential energy
function hipersurface
Energy minimum
“Conformational space”
First conformation
0
-50
Energy (kJ/mol)
-100
-150
-200
After 500 steps -250
-300
-350
-400
0 100 200 300 400 500
Energy minimisation step
Start
•Lots of minima
Energy
Local Energy minimum
“Conformational space”
n =1
dt N
12 K n ( n − 0 ) +
n
2
F = − V( r )
n =1
N
12K n ( n
2
− 0 ) +
n
n =1
N
12Kn 1+cos(mnn − n +
n =1
N at
C12r (12i , j ) − C6r( i6, j ) + 4i j r
qq
S( rij ) +
i j ij ij 0 r ij
Sampling conformational
states
𝐯1 , ⋯ , 𝐯𝑖 , ⋯ , 𝐯𝑛 Velocities F2
F1 v2
𝐚1 , ⋯ , 𝐚𝑖 , ⋯ , 𝐚𝑛 Accelerations
v1
𝑑 2 𝐫𝑖 𝐅𝑖
= 𝐚𝑖 = F4
𝑑𝑡 𝑚𝑖
F3 v4
𝑑𝐯𝑖 𝐅𝑖
= 𝐚𝑖 = v3
𝑑𝑡 𝑚𝑖
𝑑𝐫𝑖
= 𝐯𝑖
𝑑𝑡
C.M. Soares – Physical Modelling of Biomolecules – 2023
Integration of equations of motion
𝐫1 , ⋯ , 𝐫𝑖 , ⋯ , 𝐫𝑛
𝑥𝑖 The Euler integrator
𝐫𝑖 = 𝑦𝑖 𝑑𝐫𝑖
𝑧𝑖 = 𝐯𝑖
𝑑𝑡
𝐯1 , ⋯ , 𝐯𝑖 , ⋯ , 𝐯𝑛 Can be approximated by
𝐚1 , ⋯ , 𝐚𝑖 , ⋯ , 𝐚𝑛 𝐫𝑖 𝑡 + ∆𝑡 − 𝐫𝑖 𝑡
≈ 𝐯𝑖 (𝑡)
∆𝑡
F2
𝑑 2 𝐫𝑖 𝐅𝑖 F1 v2
= 𝐚𝑖 = 𝐫𝑖 𝑡 + ∆𝑡 = 𝐫𝑖 𝑡 + 𝐯𝑖 (𝑡) ∙ ∆𝑡
𝑑𝑡 𝑚𝑖
v1
𝑑𝐯𝑖 𝐅𝑖 𝐅𝑖
= 𝐚𝑖 = 𝐯𝑖 𝑡 + ∆𝑡 = 𝐯𝑖 𝑡 + ∙ ∆𝑡
𝑑𝑡 𝑚𝑖 𝑚𝑖
F4
F3 v4
𝑑𝐫𝑖 • The simplest
= 𝐯𝑖 • But unstable…
𝑑𝑡 v3
C.M. Soares – Physical Modelling of Biomolecules – 2023
Integration of equations of motion
𝐫1 , ⋯ , 𝐫𝑖 , ⋯ , 𝐫𝑛
𝑥𝑖 The Leap-Frog integrator
𝐫𝑖 = 𝑦𝑖 1
𝑧𝑖 𝐫𝑖 𝑡 + ∆𝑡 = 𝐫𝑖 𝑡 + 𝐯𝑖 (𝑡 + ∆𝑡) ∙ ∆𝑡
2
𝐯1 , ⋯ , 𝐯𝑖 , ⋯ , 𝐯𝑛
1 1 𝐅𝑖
𝐯𝑖 𝑡 + ∆𝑡 = 𝐯𝑖 𝑡 − ∆𝑡 + ∙ ∆𝑡
𝐚1 , ⋯ , 𝐚𝑖 , ⋯ , 𝐚𝑛 2 2 𝑚𝑖
F2
• Just a bit more complex
𝑑 2 𝐫𝑖 𝐅𝑖 • Quite stable for the time steps F1 v2
= 𝐚𝑖 =
𝑑𝑡 𝑚𝑖 used in MD
v1
𝑑𝐯𝑖 𝐅𝑖
= 𝐚𝑖 =
𝑑𝑡 𝑚𝑖 F4
F3 v4
𝑑𝐫𝑖
= 𝐯𝑖 Taken from the GROMACS 2016 manual
𝑑𝑡 v3
C.M. Soares – Physical Modelling of Biomolecules – 2023
Including solvent
•Implicitly
Electrostatic screening, Langevin dynamics, hydrophobic terms, etc
Smith, P. E. & Pettitt, B. M. (1994). J.Phys.Chem. 98, 9700-9711; Solmajer, T. & Mehler, E. L. (1992). Int.J.Quan.Chem.
44, 291-299; Still, W. C., Tempczyk, A., Hawley, R. C. & Hendrickson, T. (1990). J.Am.Chem.Soc. 112, 6127-6129;
Wesson, L. & Eisenberg, D. (1992). Prot.Sci. 1, 227-235.
•Explicitly
Inclusion of simulated solvent
molecules
Rc
Use of periodic
R
boundary
caixa
conditions
Ex: Berendsen, H. J. C., Postma, J. P. M., van Gunsteren, W. F., DiNola, A. & Haak, J. R. (1984). Molecular dynamics with coupling to an external bath. J. Chem.
Phys. 81, 3684-3690.
dT (t ) 1
= 2 (T0 − T ) ; =
Heat dt
“Bath” System •Velocity scaling
t T0
vi vi ; = 1+ ( − 1) using
T
has dimensions of time
Temperature coupling constant Similar for Pressure
C.M. Soares – Physical Modelling of Biomolecules – 2023
A recipe for MD
This limits the type of problems that can be solved using “brute force”
•Bond stretching 0.0001 to 0.0001 ns
•Protein sidechain rotation at the surface 0.01 to 0.1 ns
•Hinge bending motion of polypeptide chain segments 0.01 to 100 ns
•Alosteric transitions, folding, etc 10000 to 10000000000 ns
One Tyrosine
C.M. Soares – Physical Modelling of Biomolecules – 2023
Simulating a small system –Tyrosine in water
ABN2 in water
0-20ps
Faster!!!
More time
ABN2 in water
0-1000ps
0-1ns
1ns=10-9 s
ABN2, an endo-1,5-α-L-arabinanase from Bacillus subtilis (443 residues long)
Arabinanase is a glycosyl hydrolase that is able to cleave the glycosidic bonds of α-1,5-L-arabinan, releasing arabino-oligosaccharides and L-arabinose
McVey et al. (2014), JBIC, 19, 505-513.
Even
faster!!!
More time
ABN2 in water
0-10000ps
0-10ns
1ns=10-9 s
ABN2, an endo-1,5-α-L-arabinanase from Bacillus subtilis (443 residues long)
Arabinanase is a glycosyl hydrolase that is able to cleave the glycosidic bonds of α-1,5-L-arabinan, releasing arabino-oligosaccharides and L-arabinose
McVey et al. (2014), JBIC, 19, 505-513.
Exposed
ABN2 in water
0-10000ps
0-10ns
1ns=10-9 s
ABN2, an endo-1,5-α-L-arabinanase from Bacillus subtilis (443 residues long)
Arabinanase is a glycosyl hydrolase that is able to cleave the glycosidic bonds of α-1,5-L-arabinan, releasing arabino-oligosaccharides and L-arabinose
McVey et al. (2014), JBIC, 19, 505-513.
Buried
ABN2 in water
0-10000ps
0-10ns
1ns=10-9 s
ABN2, an endo-1,5-α-L-arabinanase from Bacillus subtilis (443 residues long)
Arabinanase is a glycosyl hydrolase that is able to cleave the glycosidic bonds of α-1,5-L-arabinan, releasing arabino-oligosaccharides and L-arabinose
McVey et al. (2014), JBIC, 19, 505-513.
ABN2 in water
0-10000ps
0-10ns
1ns=10-9 s
ABN2, an endo-1,5-α-L-arabinanase from Bacillus subtilis (443 residues long)
Arabinanase is a glycosyl hydrolase that is able to cleave the glycosidic bonds of α-1,5-L-arabinan, releasing arabino-oligosaccharides and L-arabinose
McVey et al. (2014), JBIC, 19, 505-513.
100 ns
ABN2 in water
10000-110000ps
10-110ns
1ns=10-9 s
ABN2, an endo-1,5-α-L-arabinanase from Bacillus subtilis (443 residues long)
Arabinanase is a glycosyl hydrolase that is able to cleave the glycosidic bonds of α-1,5-L-arabinan, releasing arabino-oligosaccharides and L-arabinose
McVey et al. (2014), JBIC, 19, 505-513.
100 ns
ABN2 in water
10000-110000ps
10-110ns
1ns=10-9 s
ABN2, an endo-1,5-α-L-arabinanase from Bacillus subtilis (443 residues long)
Arabinanase is a glycosyl hydrolase that is able to cleave the glycosidic bonds of α-1,5-L-arabinan, releasing arabino-oligosaccharides and L-arabinose
McVey et al. (2014), JBIC, 19, 505-513.
Ubiquitin
Human
Ubiquitin in water
0-1ns
1ns=10-9s
Ubiquitin in water
0-100ns=0-0.1s
1s=10-6s
•Objective: Understand
the molecular basis of
enzyme stability in
ethanol/water mixtures
using very long
(microsecond) MD
simulations
Diana Lousa
Cytochrome c3
DvH
Sav1866 from S.aureus inserted into a simulated DMPC membrane Oliveira et al. (2011) Proteins, 79, 1977-90
A prototype to
eukaryotic ABC
transporters like CTFR,
where mutations are
responsible for cystic
fibrosis
Oliveira et al. (2010) J Phys Chem B., 114, 5486-96. Damas et al (2011) Prot.Sci., 20, 1220-1230
C.M. Soares – Physical Modelling of Biomolecules – 2023
Conformational changes in bacterial ABC transporters
Simulating membrane proteins
Sav1866
Helix 7
Helix 11
Helix 15
Sav1866
Periplasmic
direction
Cytoplasmic
direction
H+
Electron flow
ADP+Pi ATP Proton flow
Cytoplasm
9hcc
c3 I
c3 II
Hydrogenase
H2
H2
Vitor H. Teixeira
•Reaction field
electrostatics
•Reaction field
electrostatics
•Reaction field
electrostatics
Alternative channel for H2 access to the active site in the Dmb [NiFeSe] Hase
V/
G
0 1
2-phenylpropanol S
+cutinase
(transition state)
R
Nuno Micaelo
<dV/d> kJ/mol
100
-100 In water
-200
-300
0.0 0.2 0.4 0.6 0.8 1.0
Micaelo et al. 2002
•Leach, A. R. (1996). Molecular modelling. Principles and applications. Addison Wesley Longman, Essex.
•Karplus, M. & Petsko, A. (1990). Molecular dynamics simulations in biology. Nature 347, 631-639.
•van Gunsteren, W. F., Weiner, P. K. & Wilkinson, A. J. (1993). Computer simulation of biomolecular systems. Theoretical and experimental
applications. Vol. 2. ESCOM, Leiden.
•van Gunsteren, W. F. & Berendsen, H. J. C. (1990). Computer simulation of molecular dynamics: Methodology, applications, and perspectives in
chemistry. Angew. Chem. Int. 29, 992-1023.
•McCammon, J. A. & Harvey, S. C. (1987). Dynamics of proteins and nucleic acids. Cambridge University Press, Cambridge.
•Brooks, C. L., Karplus, M. & Pettitt, B. M. (1988). Proteins: A theoretical perpective of dynamics, structure, and thermodinamics, 1st edition. John wiley
&Sons, New York.
•Kollman, P. (1993). Free energy calculations: applications to chemical and biochemical phenomena. Chem.Rev. 93, 2395-2417.
•Karplus M., McCammon J. A. (2002) Molecular dynamics simulations of biomolecules. Nat.Struc.Biol. 9, 646-652
•Hansson T., Oostenbrink C., van Gunsteren W. F. (2002) Molecular dynamics simulations. Curr.Opin.Struct.Biol. 12, 190-196
•van Gunsteren, W.F., D. Bakowies, R. Baron, I. Chandrasekhar, M. Christen, X. Daura, P. Gee, D.P. Geerke, A. Glattli, P.H. Hunenberger and others.
2006. Biomolecular modeling: Goals, problems, perspectives. Angew Chem Int Ed Engl 45(25):4064-92.
•Karplus, M., J. Kuriyan. 2005. Molecular dynamics and protein function. Proc Natl Acad Sci U S A 102(19):6679-85.