170 likes | 380 Views
Molecular modeling group Bioengineering department Biology Faculty M.V.Lomonosov Moscow State University. Calculation of interaction energy between voltage-gated potassium channel Kv1.2 and blocker agitoxin Valery N. Novoseletsky Maria A. Bolshakova
E N D
Molecular modeling group Bioengineering department Biology Faculty M.V.Lomonosov Moscow State University Calculation of interaction energy between voltage-gated potassium channel Kv1.2 and blocker agitoxin Valery N. Novoseletsky Maria A. Bolshakova Konstantin V. Shaitan Moscow 2013
Voltage-Dependent Potassium Channels Pore-forming membrane proteins that open or close in response to changes in the transmembrane voltage. Potassium channels are the most widely distributed type of ion channel and are found in virtually all living organisms. Action potential is a short-lasting event in which the electrical membrane potential of a cell rapidly rises and falls, following a consistent trajectory. During action potentials, voltage-dependent potassium channelsplay a crucial role in returning the depolarized cell to a resting state.
Crystal Structure of a Mammalian Voltage-Dependent K+ Channel (Long, Campbell, MacKinnon, 2005)
Crystal Structure of a Mammalian Voltage-Dependent K+ Channel Transmembrane domain scheme: helices S1-S4 form the voltage-sensing domain, helicesS5 and S6 take place in the ions conduction.
Crystal Structure of a Mammalian Voltage-Dependent K+ Channel Side view: voltage-sensing domain, pore domain and selective filter.
Scorpion venom toxin – agitoxin Agitoxin structure was determined by NMR(Kresel, Kasibhatha, 1995). Structure has highly conserved motif consisting of - 3 antiparallelβ-sheets - α-helix - 3 disulfide bridges Agitoxin specifically blocks Kv1.2 channel with high affinity (Kd < 1 nmol/L)
Kv1.2 channel complex with agitoxin Extracellular and membrane view. Highly conserved residues involved in binding are marked (pdb-codes1AGT for agitoxin, 2A79 for channel).
Potential of mean force with umbrella sampling 1. Pulling the toxin away fromthe channel using steered molecular dynamics Toxin-channel distance increases with time from 2.5 nm at 0 ps to 7 nm at 500 ps.
Potential of mean force with umbrella sampling 2. Frame selection 3.Configuration space sampling
Potential of mean force with umbrella sampling 4. PMF extraction and ΔG calculation Dependence of PMF on toxin-channel distance in 0.1 MNaCl solution.
Potential of mean force with umbrella sampling 5. Error estimation 22 ± 2kcal/mol Dependence of PMF on toxin-channel distance in 0.1 MNaCl solution. Binding energy is equal to the max PMF value and consists 22 ± 2kcal/mol.
Effect of ionic strength on the free energy of binding 0, 025M 25 20 15 10 5 0 0, 05M 0,1M 0,2M 0,4M E (kcal/mol) z (nm) Binding energy drops with the increase of ionic strength.
Comparison of the results * Zachariae et al, 2008 ** Our data *** Khabiri et al, 2011
Conclusion • Structural model of agitoxin-channel complex was constructed • Free energy binding was calculated using potential of mean force and umbrella sampling • Obtained results differ from experimental data. Probable reason is in the mistake of system preparation