molecular dynamics simulations of cro proteins mutation l.
Skip this Video
Download Presentation
Molecular Dynamics Simulations of Cro Proteins: Mutation!

Loading in 2 Seconds...

play fullscreen
1 / 44

Molecular Dynamics Simulations of Cro Proteins: Mutation! - PowerPoint PPT Presentation

  • Uploaded on

Molecular Dynamics Simulations of Cro Proteins: Mutation! . Max Shokhirev Miyashita-Tama Group 5-14-08. Background Image from 1rzs1.pdb courtesy of PDB. Overview. Background Evolution of Cro Proteins and what they are Ideas behind Molecular Dynamics (MD) Alanine Scanning Simulations

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
Download Presentation

PowerPoint Slideshow about 'Molecular Dynamics Simulations of Cro Proteins: Mutation!' - slone

Download Now 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.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
molecular dynamics simulations of cro proteins mutation
Molecular Dynamics Simulations of Cro Proteins:Mutation!

Max Shokhirev

Miyashita-Tama Group


Background Image from 1rzs1.pdb courtesy of PDB

  • Background
    • Evolution of Cro Proteins and what they are
    • Ideas behind Molecular Dynamics (MD)
  • Alanine Scanning Simulations
  • Conclusions
evolution of protein structure
Evolution of Protein Structure

Neutral Sequence Networks1

1= ancestor

2= same fold descendant

3= different fold via unstable mutations (relaxed)

4= frameshift descendant

5= different fold via stable mutations

cro proteins
Cro Proteins?
  • DNA-binding proteins
    • Initiate lytic pathway in bacteria3
    • Ancestral forms have 5 α-helices, with the 2nd and 3rd forming a helix-turn-helix DNA-binding motif (P22 Cro is an example)
    • Bacteriophage λ Cro consists of 3 α-helices and the 4th and 5th helices are replaced by a β-hairpin.
p22 vs cro
P22 vs λ Cro

P22 Croλ Cro

two approaches
Two approaches…
  • The Cro protein family has been studied with Alanine-Scanning Mutagenesis and Hybrid-Scanning Mutagenesis1
  • Computational approach
    • Molecular Dynamics
    • Data-mining 4
    • Etc.
molecular dynamics md
Molecular Dynamics (MD)
  • Deterministic
    • Given initial conditions and parameters it is possible to calculate the conditions at any other point in time.
  • Iterative (Discrete)
    • Repeat force calculations at each time step and move particles accordingly.
    • Need to pick Δt such that the particles move continuously
velocity verlet integrator
Velocity-Verlet Integrator
  • Scheme for calculating new position, velocity, and acceleration at each time step:
  • Compute New Position
  • Compute Half Velocity
  • Compute Force
  • Compute Velocity

Time step

-1 -.5 0 .5 1




initial conditions
Initial Conditions…
  • Initial Positions
    • Extracted from PDB file
  • Bonding Interactions
    • Bonding information from PDB
      • Direct bonds, allowed angles, allowed dihedrals
  • Velocity?
    • Generated using genVel based on equipartition theory at a specified temperature.
  • Other parameters
    • Masses, LJ types, Specific LJs, general simulation parameters
initial temperature
Initial Temperature…
  • The temperature is proportional to the average speed of particles in a system. We can assign temperatures based on the Maxwell-Boltzman velocity distribution function:
  • Vi = (Normalized Gaussian Random number) * sqrt((Kb*Na*T)/Mi)
temperature control
Temperature Control…
  • System is coupled to a virtual heat bath:
    • Vnew=Vold*sqrt(1-(ts/tau)*(1- Ttarget/Tcurrent))
      • ts = time step length
      • tau = coupling coefficient
force field
Force Field
  • Force on each particle calculated from components
    • Direct bond
    • Angle
    • Dihedral
    • Specific LJ
    • Non-specific LJ
bond interactions
Bond Interactions
  • V = ½k(Xi-X0)2
  • Fi = k*(Xi-X0)/Xi
lennard jones interactions
Lennard-Jones Interactions


  • Non-specific LJ
    • By atom type (6-12)
  • Specific(native) LJ
    • 6-12
    • 10-12
thus far
Thus far…
  • Phase I
    • Create a program for flexible MD simulations using a Go-like potential
    • Simulator seems to be working for bond, angle, dihedral, LJ (10-12 and 6-12). Cro proteins are folding/unfolding!
  • Phase II
    • Results from honors thesis
  • Phase III
    • Mutational studies of Cro proteins
phase ii honors thesis
Phase II – Honors Thesis
  • Cro folding and unfolding
  • Melting temperature simulations
  • Comparison of 6-12 and 10-12 LJ interactions
  • Alanine Scanning for P22 and Lambda Cro
cro folding and unfolding
Cro Folding and Unfolding

Temp = 350 Temp = 800

P22 Cro

λ Cro

cro folding and unfolding21
Cro Folding and Unfolding

T = 300

T = 300

T = 1000

calculating melting temp
Calculating Melting Temp
  • Run simulation(s) at different temps
  • Calculate Q values for each temp
    • At Tm Q values fluctuate around 0.5
    • Can plot histogram of Q values
    • Free energy profile for each temp
      • E = -Kb*T*log(P(q))
  • Calculate Specific Heat
    • Derivative of total energy plot at each temp.
  • Values are not scaled to real-world values
cro melting temperatures
λ Cro Melting Temperatures

Purple = 10-12 LJ

Orange = 6-12 LJ

melting temperature from specific heat
Melting Temperature from Specific Heat
  • We can obtain the melting temperature by plotting the specific heat as a function of simulation temperature
  • The specific heat is the derivative of the total energy function with respect to temperature
specific heats
Specific Heats


P22 Cro ~ T=750

λ Cro ~ T= 685

real melting temperatures
Real Melting Temperatures
  • λ Cro
    • 334 K1
      • Oligomer with Tm <= 313 K1
      • λ Cro A33W/F58D pure monomer
  • P22 Cro
    • 327 K1
melting temperature conc
Melting Temperature Conc.
  • P22 Cro ~ 745/750
  • λ Cro ~ 690/685
  • P22 Cro is a 2-state folder, λ Cro is not!

P22 λ Cro

test effect of lj10 12 pot
Test Effect of LJ10-12 pot.
  • Simulations performed on P22 Cro and λ Cro under nearly identical conditions
    • Change the Lennard-Jones potential from a 6-12 pot to a 10-12 potential.
    • This should theoretically increase “cooperativity” of folding2
lj10 12 results
LJ10-12 Results

6-12 LJ Potential 10-12 LJ Potential



lj observations
LJ Observations…
  • The melting temperatures decreased when using a 10-12 LJ potential.
  • The 10-12 LJ Potential shows a higher degree of cooperativity (esp for P22)
alanine scanning
Alanine Scanning
  • Mutate the structurally divergent residues to alanine.
  • Remove the native contacts for each residue.
  • Simulations at the folding temperature of each Cro protein.
  • Average Q values for each residue
alanine scanning results
Alanine Scanning Results
  • Alanine Scanning simulations match melting temperature data
  • Alanine Scanning simulations show regions that decrease stability, which does not match the real data.
phase iii cro mutation studies
Phase III – Cro Mutation Studies
  • What drives structural stability?
    • Native interactions
    • Native interactions (between divergent and not divergent domains)
    • Dihedral Interactions
    • Angle Interactions (the future)
removing native dihedrals
Removing native + dihedrals

1rzs: mykkdvidhf gtqravakal gisdaavsqw kevipekday rleivtagal kyqenayrqa a

5cro: meqritlkdyamrf gqtktakdlg vyqsainka- --ihagrkif ltinadgsvy aeevkpfpsn kktta

removing native mixing dihedrals
Removing Native/Mixing Dihedrals

1rzs: mykkdvidhf gtqravakal gisdaavsqw kevipekday rleivtagal kyqenayrqa a

5cro: meqritlkdyamrf gqtktakdlg vyqsainka- --ihagrkif ltinadgsvy aeevkpfpsn kktta

removed inter domain native cont
Removed Inter-domain native cont.

1rzs: mykkdvidhf gtqravakal gisdaavsqw kevipekday rleivtagal kyqenayrqa a

5cro: meqritlkdyamrf gqtktakdlg vyqsainka- --ihagrkif ltinadgsvy aeevkpfpsn kktta

Purple Lambda 6-12 LJ

Gray Lambda 10-12 LJ

Red P22 10-12 LJ

Black P22 6-12 LJ

removing dihedral angles only
Removing Dihedral Angles Only

1rzs: mykkdvidhf gtqravakal gisdaavsqw kevipekday rleivtagal kyqenayrqa a

5cro: meqritlkdyamrf gqtktakdlg vyqsainka- --ihagrkif ltinadgsvy aeevkpfpsn kktta

  • An MD Simulation program was written to study Cro proteins
  • P22 has been shown to unfold and refold as a function of temperature.
  • Folding temperatures observed from free energy profile and specific heat data.
  • λ Cro has only one free energy minimum at its folding temperature, while 2 minima are observed for P22 Cro.
  • The 10-12 LJ interaction allows for higher cooperativity.
  • Alanine scanning simulations qualitatively match real data.
  • Dihedral angle interactions are essential to stability of mutants

Dr. Osamu Miyashita

Dr. Florence Tama

M-T Group

  • "Relationship between sequence determinants of stability for two natural homologous proteins with different folds", L.O. Van Dorn, T. Newlove, S. Chang, W.M. Ingram, and M.H.J. Cordes. Biochemistry.45, 10542–10553 (2006).
  • “Scrutinizing the squeezed exponential kinetics observed in the folding simulation of an off-lattice Go-like protein model”, H. K. Nakamura, M.Sasai, M Takano. Chemical Physics.307 259–267 (2004).
  • “Mechanism of action of the cro protein of bacteriophage lambda.” A Johnson, B J Meyer, and M Ptashne. Proc Natl Acad Sci U S A. 75(4): 1783–1787 (1978).
  • "High polar content of long buried blocks of sequence in protein domains suggests selection against amyloidogenic nonpolar sequences", A.U. Patki, A.C. Hausrath, and M.H.J. Cordes. Journal of Molecular Biology. 362, 800–809 (2006).

Images Used: