Molecular dynamics study of Ca2+ binding loop variants of silver hake parvalbumin with aspartic acid at the "gateway" position

Kamau Fahie, Rebecca Pitts, Kelly M. Elkins, Donald J. Nelson

Research output: Contribution to journalArticle

Abstract

The helix-loop-helix (i.e., EF-hand) Ca2+ ion binding motif is characteristic of a large family of high-affinity calcium ion binding proteins, including the parvalbumins, oncomodulins and calmodulins. In this work we describe a set of molecular dynamics computations on the major parvalbumin from the silver hake (SHPV-B) and on functional fragments of this protein, consisting of the first four helical regions (the ABCD fragment), and the internal helix-loop-helix region (the CD fragment). In both whole protein and protein fragments (i.e., ABCD and CD fragments), the 9th loop residue in the calcium ion binding site in the CD helix-loop-helix region (the so-called "gateway" position) has been mutated from glutamic acid to aspartic acid. Aspartic acid is one of the most common residues found at the gateway position in other (non-parvalbumin) EF-hand proteins, but has never been found at the gateway position of any parvalbumin. (Interestingly, aspartic acid does occur at the gateway position in the closely related rat and human oncomodulins.) Consistent with experimental observations, the results of our molecular dynamics simulations show that incorporation of aspartic acid at the gateway position is very disruptive to the structural integrity of the calcium ion coordination site in the whole protein. The aspartic acid mutation is somewhat less disruptive to the calcium ion coordination sites in the two parvalbumin fragments (i.e., the ABCD and CD fragments), presumably due to the higher degree of motional freedom allowable in these protein fragments. One problem associated with the E59D whole protein variant is a prohibitively close approach of the aspartate carboxyl group to the CD calcium ion observed in the energy-minimized (pre-molecular dynamics) structure. This steric situation does not emerge during energy-minimization of the wild-type protein. The damage to the structural integrity of the calcium ion coordination site in the whole protein E59D variant is not relieved during the molecular dynamics simulation. In fact, during the course of the 300 picosecond simulation, all of the calcium ion ligands leave the primary coordination sphere. In addition, the conserved hydrogen-bonds (in the short β-sheet structure) that links the CD site to the symmetry-related EF site (in the non-mutated whole protein) is also somewhat disrupted in the E59D whole protein variant. These results suggest that the Ca2+ ion binding deficiencies in the CD loop are related, at least in part, to the unique interaction that exists between the paired CD and EF hands in the whole protein. Our theoretical results correlate well with previous studies on engineered EF-hand proteins and with all of our experimental evidence on whole silver hake parvalbumin and enzymatically-generated parvalbumin fragments.

Original languageEnglish (US)
Pages (from-to)821-837
Number of pages17
JournalJournal of Biomolecular Structure and Dynamics
Volume19
Issue number5
DOIs
StatePublished - Jan 1 2002
Externally publishedYes

Fingerprint

Parvalbumins
Molecular Dynamics Simulation
Silver
Aspartic Acid
Ions
EF Hand Motifs
Proteins
Calcium
Molecular Computers
Calcium-Binding Proteins
Calmodulin
Molecular Structure
Glutamic Acid
Hydrogen

ASJC Scopus subject areas

  • Structural Biology
  • Molecular Biology

Cite this

Molecular dynamics study of Ca2+ binding loop variants of silver hake parvalbumin with aspartic acid at the "gateway" position. / Fahie, Kamau; Pitts, Rebecca; Elkins, Kelly M.; Nelson, Donald J.

In: Journal of Biomolecular Structure and Dynamics, Vol. 19, No. 5, 01.01.2002, p. 821-837.

Research output: Contribution to journalArticle

@article{fcbe6bd9a6d641cb9f8cfee3a99254a4,
title = "Molecular dynamics study of Ca2+ binding loop variants of silver hake parvalbumin with aspartic acid at the {"}gateway{"} position",
abstract = "The helix-loop-helix (i.e., EF-hand) Ca2+ ion binding motif is characteristic of a large family of high-affinity calcium ion binding proteins, including the parvalbumins, oncomodulins and calmodulins. In this work we describe a set of molecular dynamics computations on the major parvalbumin from the silver hake (SHPV-B) and on functional fragments of this protein, consisting of the first four helical regions (the ABCD fragment), and the internal helix-loop-helix region (the CD fragment). In both whole protein and protein fragments (i.e., ABCD and CD fragments), the 9th loop residue in the calcium ion binding site in the CD helix-loop-helix region (the so-called {"}gateway{"} position) has been mutated from glutamic acid to aspartic acid. Aspartic acid is one of the most common residues found at the gateway position in other (non-parvalbumin) EF-hand proteins, but has never been found at the gateway position of any parvalbumin. (Interestingly, aspartic acid does occur at the gateway position in the closely related rat and human oncomodulins.) Consistent with experimental observations, the results of our molecular dynamics simulations show that incorporation of aspartic acid at the gateway position is very disruptive to the structural integrity of the calcium ion coordination site in the whole protein. The aspartic acid mutation is somewhat less disruptive to the calcium ion coordination sites in the two parvalbumin fragments (i.e., the ABCD and CD fragments), presumably due to the higher degree of motional freedom allowable in these protein fragments. One problem associated with the E59D whole protein variant is a prohibitively close approach of the aspartate carboxyl group to the CD calcium ion observed in the energy-minimized (pre-molecular dynamics) structure. This steric situation does not emerge during energy-minimization of the wild-type protein. The damage to the structural integrity of the calcium ion coordination site in the whole protein E59D variant is not relieved during the molecular dynamics simulation. In fact, during the course of the 300 picosecond simulation, all of the calcium ion ligands leave the primary coordination sphere. In addition, the conserved hydrogen-bonds (in the short β-sheet structure) that links the CD site to the symmetry-related EF site (in the non-mutated whole protein) is also somewhat disrupted in the E59D whole protein variant. These results suggest that the Ca2+ ion binding deficiencies in the CD loop are related, at least in part, to the unique interaction that exists between the paired CD and EF hands in the whole protein. Our theoretical results correlate well with previous studies on engineered EF-hand proteins and with all of our experimental evidence on whole silver hake parvalbumin and enzymatically-generated parvalbumin fragments.",
author = "Kamau Fahie and Rebecca Pitts and Elkins, {Kelly M.} and Nelson, {Donald J.}",
year = "2002",
month = "1",
day = "1",
doi = "10.1080/07391102.2002.10506787",
language = "English (US)",
volume = "19",
pages = "821--837",
journal = "Journal of Biomolecular Structure and Dynamics",
issn = "0739-1102",
publisher = "Adenine Press",
number = "5",

}

TY - JOUR

T1 - Molecular dynamics study of Ca2+ binding loop variants of silver hake parvalbumin with aspartic acid at the "gateway" position

AU - Fahie, Kamau

AU - Pitts, Rebecca

AU - Elkins, Kelly M.

AU - Nelson, Donald J.

PY - 2002/1/1

Y1 - 2002/1/1

N2 - The helix-loop-helix (i.e., EF-hand) Ca2+ ion binding motif is characteristic of a large family of high-affinity calcium ion binding proteins, including the parvalbumins, oncomodulins and calmodulins. In this work we describe a set of molecular dynamics computations on the major parvalbumin from the silver hake (SHPV-B) and on functional fragments of this protein, consisting of the first four helical regions (the ABCD fragment), and the internal helix-loop-helix region (the CD fragment). In both whole protein and protein fragments (i.e., ABCD and CD fragments), the 9th loop residue in the calcium ion binding site in the CD helix-loop-helix region (the so-called "gateway" position) has been mutated from glutamic acid to aspartic acid. Aspartic acid is one of the most common residues found at the gateway position in other (non-parvalbumin) EF-hand proteins, but has never been found at the gateway position of any parvalbumin. (Interestingly, aspartic acid does occur at the gateway position in the closely related rat and human oncomodulins.) Consistent with experimental observations, the results of our molecular dynamics simulations show that incorporation of aspartic acid at the gateway position is very disruptive to the structural integrity of the calcium ion coordination site in the whole protein. The aspartic acid mutation is somewhat less disruptive to the calcium ion coordination sites in the two parvalbumin fragments (i.e., the ABCD and CD fragments), presumably due to the higher degree of motional freedom allowable in these protein fragments. One problem associated with the E59D whole protein variant is a prohibitively close approach of the aspartate carboxyl group to the CD calcium ion observed in the energy-minimized (pre-molecular dynamics) structure. This steric situation does not emerge during energy-minimization of the wild-type protein. The damage to the structural integrity of the calcium ion coordination site in the whole protein E59D variant is not relieved during the molecular dynamics simulation. In fact, during the course of the 300 picosecond simulation, all of the calcium ion ligands leave the primary coordination sphere. In addition, the conserved hydrogen-bonds (in the short β-sheet structure) that links the CD site to the symmetry-related EF site (in the non-mutated whole protein) is also somewhat disrupted in the E59D whole protein variant. These results suggest that the Ca2+ ion binding deficiencies in the CD loop are related, at least in part, to the unique interaction that exists between the paired CD and EF hands in the whole protein. Our theoretical results correlate well with previous studies on engineered EF-hand proteins and with all of our experimental evidence on whole silver hake parvalbumin and enzymatically-generated parvalbumin fragments.

AB - The helix-loop-helix (i.e., EF-hand) Ca2+ ion binding motif is characteristic of a large family of high-affinity calcium ion binding proteins, including the parvalbumins, oncomodulins and calmodulins. In this work we describe a set of molecular dynamics computations on the major parvalbumin from the silver hake (SHPV-B) and on functional fragments of this protein, consisting of the first four helical regions (the ABCD fragment), and the internal helix-loop-helix region (the CD fragment). In both whole protein and protein fragments (i.e., ABCD and CD fragments), the 9th loop residue in the calcium ion binding site in the CD helix-loop-helix region (the so-called "gateway" position) has been mutated from glutamic acid to aspartic acid. Aspartic acid is one of the most common residues found at the gateway position in other (non-parvalbumin) EF-hand proteins, but has never been found at the gateway position of any parvalbumin. (Interestingly, aspartic acid does occur at the gateway position in the closely related rat and human oncomodulins.) Consistent with experimental observations, the results of our molecular dynamics simulations show that incorporation of aspartic acid at the gateway position is very disruptive to the structural integrity of the calcium ion coordination site in the whole protein. The aspartic acid mutation is somewhat less disruptive to the calcium ion coordination sites in the two parvalbumin fragments (i.e., the ABCD and CD fragments), presumably due to the higher degree of motional freedom allowable in these protein fragments. One problem associated with the E59D whole protein variant is a prohibitively close approach of the aspartate carboxyl group to the CD calcium ion observed in the energy-minimized (pre-molecular dynamics) structure. This steric situation does not emerge during energy-minimization of the wild-type protein. The damage to the structural integrity of the calcium ion coordination site in the whole protein E59D variant is not relieved during the molecular dynamics simulation. In fact, during the course of the 300 picosecond simulation, all of the calcium ion ligands leave the primary coordination sphere. In addition, the conserved hydrogen-bonds (in the short β-sheet structure) that links the CD site to the symmetry-related EF site (in the non-mutated whole protein) is also somewhat disrupted in the E59D whole protein variant. These results suggest that the Ca2+ ion binding deficiencies in the CD loop are related, at least in part, to the unique interaction that exists between the paired CD and EF hands in the whole protein. Our theoretical results correlate well with previous studies on engineered EF-hand proteins and with all of our experimental evidence on whole silver hake parvalbumin and enzymatically-generated parvalbumin fragments.

UR - http://www.scopus.com/inward/record.url?scp=0036250233&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=0036250233&partnerID=8YFLogxK

U2 - 10.1080/07391102.2002.10506787

DO - 10.1080/07391102.2002.10506787

M3 - Article

C2 - 11922838

AN - SCOPUS:0036250233

VL - 19

SP - 821

EP - 837

JO - Journal of Biomolecular Structure and Dynamics

JF - Journal of Biomolecular Structure and Dynamics

SN - 0739-1102

IS - 5

ER -