水稻OsGAP1 C2结构域中对钙离子的结合能力强于对钾离子北大核心CSCD

钙离子作为植物细胞的第二信使,广泛参与植物应对不同逆境胁迫的信号调控过程。水稻G蛋白促进蛋白1(Oryza sativa GTPase-activating protein 1,OsGAP1)包含1个C2结构域,而含C2结构域的蛋白质是一类钙离子结合蛋白质,受钙信号的调控。本研究鉴定了水稻OsGAP1的由5个保守性天冬氨酸残基组成的阳离子结合区域。该区域可结合2个钙离子或者钾离子,且其结合钙离子的强度高于其结合钾离子的强度,但是不能结合镁离子。当将其中2个保守的天冬氨酸残基(Asp-23和Asp-28)突变为丙氨酸后,其对钙离子的结合能力减弱。对OsGAP1 C2结构域阳离子结合区域结合金属离子能力的研究,有助于加深对钙信号调控蛋白质的认识,为其在农业生产中的应用提供理论依据。     

PICK1中N端酸性区域对其PDZ结构域脂质结合能力的调节

蛋白激酶Cα相互作用蛋白1(protein interacting with Cα kinase 1, PICK1)是衔接膜上受体和蛋白激酶Cα的重要蛋白.利用荧光光谱结合定点突变技术 、蛋白与脂质覆盖法等方法,分析了PICK1蛋白N末端区域几个酸性氨基酸残基对PDZ 结构域与膜脂结合的影响,以及钙离子结合N末端酸性区域对PDZ脂结合能力的调节. 结果显示, 带有上游酸性区域的PDZ结构域（NPDZ）的脂质结合能力仅相当PDZ结构 域的15%,相比单独的PDZ结构域与脂质的解离常数Kd(PDZ)为1.58×103 μg·L-1, NPDZ与脂质解离常数Kd(NPDZ)为3.3×104μg·L-1,其中在N末端酸性残基中D8与 D12两个天冬氨酸是影响脂质结合能力减弱的关键残基,若将二者分别突变为丙氨酸 后,NPDZ与脂质的解离常数分别为：Kd (D8/A)=4.42×103μg·L-1;Kd (D12/A) =1.73×103μg·L-1接近于PDZ结构域与脂质结合能力;钙离子会增强NPDZ脂结合能力,当钙离子浓度达到30 μmol/L时,NPDZ的脂结合能力提高2.3倍,但只相当于PDZ的50% 的结合能力.     

植物PP2C蛋白磷酸酶ABA信号转导及逆境胁迫调控机制研究进展

蛋白磷酸酶(protein phosphatase,PP)是蛋白质可逆磷酸化调节机制中的关键酶,而PP2C磷酸酶是一类丝氨酸/苏氨酸残基蛋白磷酸酶,是高等植物中最大的蛋白磷酸酶家族,包含76个家族成员,广泛存在于生物体中。迄今为止,在植物体内已经发现了4种PP2C蛋白磷酸酶。蛋白激酶和蛋白磷酸酶协同催化蛋白质可逆磷酸化,在植物体内信号转导和生理代谢中起着重要的调节作用,蛋白质的磷酸化几乎存在于所有的信号转导途径中。大量研究表明,PP2Cs参与多条信号转导途径,包括PP2C参与ABA调控,对干旱、低温、高盐等逆境胁迫的响应,参与植物创伤和种子休眠或萌发等信号途径,其调控机制不同,但酶催化活性都依赖于Mg2+或Mn2+的浓度。植物PP2C蛋白的C端催化结构域高度保守,而N端功能各异。文中还综述了高等植物PP2C的分类、结构、ABA受体与PP2Cs蛋白互作、PP2C基因参与ABA信号途径以及其他逆境信号转导途径的研究进展。     

WRKY转录因子超家族的研究

WRKY转录因子是一类能与W盒特异结合的DNA结合蛋白,最初从植物中分离获得,该家族因子均含有一个或两个保守的WRKY结构域,该结构域约含有60个氨基酸残基,在WRKYGQK残基核心序列之后接有一个C2H2或C2HC类型的锌指基序。WRKY转录因子在高等植物中形成一个庞大的基因家族,基因数量众多。大量的实验证据说明,WRKY蛋白参与植物的抗病反应,并影响植物的衰老、抗胁迫能力以及生长和发育。     

水稻WRKY转录调控因子研究进展

WRKY转录调控因子是植物中最大的转录调控因子家族之一,包含1或2个大约由60个氨基酸残基组成的高度保守的WRKY结构域。WRKY结构域的N端有严格保守的WRKYGQK氨基酸序列,其C端有锌指结构模型Cys(2)-His(2)或Cys(2)-His Cys。WRKY转录调控因子与靶基因启动子区域的DNA序列(T)TGAC(C/T)(W盒)特异性结合,调节目的基因的表达,在调控植物生长发育、物质代谢、抗病耐逆及氧化衰老过程中起重要作用。我们就有关水稻WRKY转录调控因子的结构分类及其生物学功能的研究进展做简要综述。     

抑癌蛋白PTEN及PDZ蛋白对其的调节

pten基因是迄今为止发现的第1个具有双特异性磷酸酶活性的抑癌基因,该基因的编码产物PTEN蛋白,是具有蛋白与脂质磷酸酯酶活性的双特异性磷酸酯酶,作为1种重要的信号分子参与细胞增殖、分化、黏附、迁移、凋亡以及基因转录的调控. 最近,关于PTEN在信号转导中的作用以及细胞内PTEN的调节机制研究较多,尤其是PDZ蛋白对PTEN的调节作用. PTEN蛋白包括1个氨基端（N端）磷酸酯酶区域,1个与脂质结合的C2区域和1个含有PDZ结合序列的羧基端（C端）区域. PDZ结构域通过识别目标蛋白羧基端PDZ结合序列与目标蛋白相互作用,调控多种重要的细胞生理过程和信号传导途径.本文就抑癌基因pten编码产物PTEN蛋白的结构、PTEN的生物学功能和PDZ蛋白对PTEN调节的研究进展进行综述.     

ROPs:植物细胞内多种信号通路的分子开关

植物RHO相关蛋白GTPases(RHO-related GTPases of plants, ROPs)是广泛存在于植物中的一类信号转导G蛋白(又称GTP结合蛋白),其通过结合GDP或GTP在非活性和活性状态间进行切换,进而在细胞极性控制、形态发育、激素水平调控、逆境反应等诸多植物生命活动的信号转导过程中扮演重要的分子开关角色。本文对ROP蛋白的结构域及基于蛋白质结构分类进行了介绍,并对拟南芥、玉米、水稻和大麦中的ROP家族蛋白质进行了系统进化分析。分析结果表明,这些植物中的ROP蛋白根据蛋白质结构域组成可分为Ⅰ类(typeⅠ)和Ⅱ类(typeⅡ)两种类型,而根据蛋白质序列的保守性可将其在植物中的ROP蛋白划分为4个进化枝。本综述不但对ROP蛋白作为分子开关在细胞内调控各种信号通路的机制进行了叙述,还对ROP在花粉管、根毛及植物表皮铺盖细胞极性发育,以及其他抗逆反应中的具体作用和机制及研究进展进行了阐述。本文还对ROP蛋白在ABA、IAA、BR等植物激素信号传导过程中的调控作用及研究进展进行了阐述。本文对植物ROP蛋白研究过程中尚未解决的问题,例如不同的ROP蛋白在同一个信号通路中的作用为何如此不同,以及ROP是如何协调不同的信号通路以共同调控一个植物发育或者生理过程等问题进行了总结,并在此基础上对未来的研究方向进行了展望。     

一种新的平滑肌调控蛋白Calponin

Calponin为一种平滑肌特有的调控蛋白, 分子克隆的证据表明, 它具有两种亚型, α型和β型分别由292和252个氨基酸组成它能与肌动蛋白结合, 抑制肌球蛋白ATP酶活性和平滑肌收缩其与肌动蛋白结合域有38个氨基酸残基(第145～182位), 丝氨酸175在调宁蛋白与肌动蛋白的相互作用中起重要作用, 它还能与钙调蛋白结合, 呈钙依赖性, 其结构域在第52～144位残基调宁蛋白的机能受磷酸化与脱磷酸化的调节.     

高等植物中蛋白磷酸酶2C的结构与功能

蛋白质磷酸化／去磷酸化是生物信号级联传递的重要方式之一,主要通过生化性质互为对立的蛋白激酶和蛋白磷酸酶实现。蛋白磷酸酶2C(PP2C)是蛋白磷酸酶的一个分支,其生化性质、蛋白质组成与结构都和其他磷酸酶显著不同,但都在生物信号传递中扮演重要角色。高等植物中PP2C广泛参与脱落酸(ABA)的各种信号途径,包括ABA诱导的种子萌发／休眠、保卫细胞及离子通道调控和气孔关闭、逆境胁迫等。PP2C也多样地参与植物创伤反应、生长发育以及抗病性等各个途径。作为大多数信号途径的负调控因子,PP2C能直接与激酶结合,与其他调控蛋白结合,以及直接与DNA结合调控相关基因的表达。     

植物RACK1蛋白研究进展

RACK1（蛋白激酶C受体）是一种色氨酸-天门冬氨酸域（WD40结构）重复蛋白。它是一种多功能支架蛋白, 结合来自不同转导通路的信号分子并在多种哺乳动物发育过程中起关键作用。在植物中也存在RACK1同源基因, 如拟南芥基因组有3个编码RACK1蛋白质的基因, 这3个蛋白质与哺乳动物RACK1在氨基酸水平的相似性都超过75%。此外, 植物RACK1蛋白质包含的WD40数量、位置和蛋白激酶C结合位点的结构域在很大程度上是保守的。该文对植物RACK1蛋白的发现、结构及其在信号转导方面的功能进行综述。     

Crystallographic identification of Ca2+ and Sr2+ coordination sites in synaptotagmin I C2B domain

Synaptotagmin I has two tandem Ca(2+)-binding C(2) domains, which are essential for fast synchronous synaptic transmission in the central nervous system. We have solved four crystal structures of the C(2)B domain, one of them in the cation-free form at 1.50 A resolution, two in the Ca(2+)-bound form at 1.04 A (two bound Ca(2+) ions) and 1.65 A (three bound Ca(2+) ions) resolution and one in the Sr(2+)-bound form at 1.18 A (one bound Sr(2+) ion) resolution. The side chains of four highly conserved aspartic acids (D303, D309, D363, and D365) and two main chain oxygens (M302:O and Y364:O), together with water molecules, are in direct contact with two bound Ca(2+) ions (sites 1 and 2). At higher Ca(2+) concentrations, the side chain of N333 rotates and cooperates with D309 to generate a third Ca(2+) coordination site (site 3). Divalent cation binding sites 1 and 2 in the C(2)B domain were previously identified from NMR NOE patterns and titration studies, supplemented by site-directed mutation analysis. One difference between the crystal and NMR studies involves D371, which is not involved in coordination with any of the identified Ca(2+) sites in the crystal structures, while it is coordinated to Ca(2+) in site 2 in the NMR structure. In the presence of Sr(2+), which is also capable of triggering exocytosis, but with lower efficiency, only one cation binding site (site 1) was occupied in the crystallographic structure.     

Insights into how CUB domains can exert specific functions while sharing a common fold: conserved and specific features of the CUB1 domain contribute to the molecular basis of procollagen C-proteinase enhancer-1 activity

Procollagen C-proteinase enhancers (PCPE-1 and -2) are extracellular glycoproteins that can stimulate the C-terminal processing of fibrillar procollagens by tolloid proteinases such as bone morphogenetic protein-1. They consist of two CUB domains (CUB1 and -2) that alone account for PCPE-enhancing activity and one C-terminal NTR domain. CUB domains are found in several extracellular and plasma membrane-associated proteins, many of which are proteases. We have modeled the structure of the CUB1 domain of PCPE-1 based on known three-dimensional structures of CUB-containing proteins. Sequence alignment shows conserved amino acids, notably two acidic residues (Asp-68 and Asp-109) involved in a putative surface-located calcium binding site, as well as a conserved tyrosine residue (Tyr-67). In addition, three residues (Glu-26, Thr-89, and Phe-90) are found only in PCPE CUB1 domains, in putative surface-exposed loops. Among the conserved residues, it was found that mutations of Asp-68 and Asp-109 to alanine almost completely abolished PCPE-1 stimulating activity, whereas mutation of Tyr-67 led to a smaller reduction of activity. Among residues specific to PCPEs, mutation of Glu-26 and Thr-89 had little effect, whereas mutation of Phe-90 dramatically decreased the activity. Changes in activity were paralleled by changes in binding of different PCPE-1 mutants to a mini-procollagen III substrate, as shown by surface plasmon resonance. We conclude that PCPE-stimulating activity requires a calcium binding motif in the CUB1 domain that is highly conserved among CUB-containing proteins but also that PCPEs contain specific sites that could become targets for the development of novel anti-fibrotic therapies.     

Amino acid residues outside of the pore region contribute to N-type calcium channel permeation

It is widely believed that the selectivity of voltage-dependent calcium channels is mainly controlled by amino acid residues contained within four p-loop motifs forming the pore of the channel. An examination of the amino acid sequences of high voltage-activated calcium channels reveals that their domain III S5-H5 regions contain a highly conserved motif with homology to known EF hand calcium binding proteins, hinting that this region may contribute to channel permeation. To test this hypothesis, we used site-directed mutagenesis to replace three conserved negatively charged residues in the N-type calcium channel alpha1B subunit (Glu-1321, Asp-1323, and Glu-1332) with positively charged amino acids (lysine and arginine) and studied their effect on ion selectivity using whole cell and single channel patch clamp recordings. Whereas the wild type channels conducted barium much more effectively than calcium, the mutant displayed nearly equal permeabilities for these two ions. Individual replacement of residue 1332 or a double substitution of residues 1321 and 1323 with lysine and arginine, respectively, were equally effective. Disruption of the putative EF hand motif through replacement of the central glycine residue (1326) with proline resulted in a similar effect, indicating that the responses observed with the triple mutant were not due to changes in the net charge of the channel. Overall, our data indicate that residues outside of the narrow region of the pore have the propensity to contribute to calcium channel permeation. They also raise the possibility that interactions of calcium ions with a putative calcium binding domain at the extracellular side of the channel may underlie the differential permeabilities of the channel for barium and calcium ions.     

Site-directed Mutagenesis Shows the Significance of Interactions with Phospholipids and the G-protein OsYchF1 for the Physiological Functions of the Rice GTPase-activating Protein 1 (OsGAP1)

The C2 domain is one of the most diverse phospholipid-binding domains mediating cellular signaling. One group of C2-domain proteins are plant-specific and are characterized by their small sizes and simple structures. We have previously reported that a member of this group, OsGAP1, is able to alleviate salt stress and stimulate defense responses, and bind to both phospholipids and an unconventional G-protein, OsYchF1. Here we solved the crystal structure of OsGAP1 to a resolution of 1.63 Å. Using site-directed mutagenesis, we successfully differentiated between the clusters of surface residues that are required for binding to phospholipids versus OsYchF1, which, in turn, is critical for its role in stimulating defense responses. On the other hand, the ability to alleviate salt stress by OsGAP1 is dependent only on its ability to bind OsYchF1 and is independent of its phospholipid-binding activity.     

How an epidermal growth factor (EGF)-like domain binds calcium. High resolution NMR structure of the calcium form of the NH2-terminal EGF-like domain in coagulation factor X.

Domains homologous to the epidermal growth factor (EGF) are important building blocks for extracellular proteins. Proteins containing these domains have been shown to function in such diverse biological processes as blood coagulation, complement activation, and the developmental determination of embryonic cell fates. Many of these proteins require calcium for their biological function. In the case of coagulation factors IX and X and anticoagulants proteins C and S, calcium has been found to bind to the EGF-like domains. We have now determined the three-dimensional structure of the calcium-bound form of the NH2-terminal EGF-like domain in coagulation factor X by two-dimensional NMR and simulated folding. Ligands to the calcium ion are the two backbone carbonyls in Gly-47 and Gly-64, as well as the side chains in Gln-49, erythro-beta-hydroxyaspartic acid (Hya) 63, and possibly Asp-46. The conserved Asp-48 is not a ligand in our present structures. The remaining ligands are assumed to be solvent molecules or, in the intact protein, ligands from neighboring domains. Other proteins interacting in a calcium-dependent manner may also contribute ligands. A comparison with the calcium-free form shows that calcium binding induces strictly local structural changes in the domain. Residues corresponding to the side chain ligands in factor X are conserved in many other proteins, such as the integral membrane protein TAN-1 of human lymphocytes and its developmentally important homolog, Notch, in Drosophila. Calcium binding to EGF-like domains may be crucial for numerous protein-protein interactions involving EGF-like domains in coagulation factors, plasma proteins, and membrane proteins. Therefore, there is reason to believe that this novel calcium site plays an important role in the biochemistry of extracellular proteins.     

The role of Asp-49 and other conserved amino acids in phospholipases A2 and their importance for enzymatic activity

The role of aspartic acid-49 (Asp-49) in the active site of porcine pancreatic phospholipase A2 was studied by recombinant DNA techniques: two mutant proteins were constructed containing either glutamic acid (Glu) or lysine (Lys) at position 49. Enzymatic characterization indicated that the presence of Asp-49 is essential for effective hydrolysis of phospholipids. Conversion of Asp-49 to either Glu or Lys strongly reduces the binding of Ca2+ ions, in particular for the lysine mutant, but the affinity for substrate analogues is hardly affected. Extensive purification of naturally occurring Lys-49 phospholipase A2 from the venom of Agkistrodon piscivorus piscivorus yielded a protein that was nearly inactive. Inhibition studies showed that this residual activity was due to a small amount of contaminating enzyme and that the Lys-49 homologue itself has no enzymatic activity. Our results indicate that Asp-49 is essential for the catalytic action of phospholipase A2. The importance of Asp-49 was further evaluated by comparison of the primary sequences of 53 phospholipases A2 and phospholipase homologues showing that substitutions at position 49 are accompanied by structural variations of otherwise conserved residues. The occurrence of several nonconserved substitutions appeared to be a general characteristic of nonactive phospholipase A2 homologues.     

Additional binding sites for anionic phospholipids and calcium ions in the crystal structures of complexes of the C2 domain of protein kinase calpha

The C2 domain of protein kinase Calpha (PKCalpha) corresponds to the regulatory sequence motif, found in a large variety of membrane trafficking and signal transduction proteins, that mediates the recruitment of proteins by phospholipid membranes. In the PKCalpha isoenzyme, the Ca2+-dependent binding to membranes is highly specific to 1,2-sn-phosphatidyl-l-serine. Intrinsic Ca2+ binding tends to be of low affinity and non-cooperative, while phospholipid membranes enhance the overall affinity of Ca2+ and convert it into cooperative binding. The crystal structure of a ternary complex of the PKCalpha-C2 domain showed the binding of two calcium ions and of one 1,2-dicaproyl-sn-phosphatidyl-l-serine (DCPS) molecule that was coordinated directly to one of the calcium ions. The structures of the C2 domain of PKCalpha crystallised in the presence of Ca2+ with either 1,2-diacetyl-sn-phosphatidyl-l-serine (DAPS) or 1,2-dicaproyl-sn-phosphatidic acid (DCPA) have now been determined and refined at 1.9 A and at 2.0 A, respectively. DAPS, a phospholipid with short hydrocarbon chains, was expected to facilitate the accommodation of the phospholipid ligand inside the Ca2+-binding pocket. DCPA, with a phosphatidic acid (PA) head group, was used to investigate the preference for phospholipids with phosphatidyl-l-serine (PS) head groups. The two structures determined show the presence of an additional binding site for anionic phospholipids in the vicinity of the conserved lysine-rich cluster. Site-directed mutagenesis, on the lysine residues from this cluster that interact directly with the phospholipid, revealed a substantial decrease in C2 domain binding to vesicles when concentrations of either PS or PA were increased in the absence of Ca2+. In the complex of the C2 domain with DAPS a third Ca2+, which binds an extra phosphate group, was identified in the calcium-binding regions (CBRs). The interplay between calcium ions and phosphate groups or phospholipid molecules in the C2 domain of PKCalpha is supported by the specificity and spatial organisation of the binding sites in the domain and by the variable occupancies of ligands found in the different crystal structures. Implications for PKCalpha activity of these structural results, in particular at the level of the binding affinity of the C2 domain to membranes, are discussed.     

Identification of a phosphorylation site and functional analysis of conserved aspartic acid residues of OmpR, a transcriptional activator for ompF and ompC in Escherichia coli

In Escherichia coli the OmpR and EnvZ proteins regulate the expression of the outer membrane porin proteins OmpC and OmpF. EnvZ and OmpR belong to a family of sensor/effector protein pairs that control adaptation to a variety of environmental conditions. EnvZ acts as the sensor protein that phosphorylates OmpR, which in turn regulates porin gene expression. The level of phosphorylated OmpR appears to be a determining factor for ompC and ompF regulation. Phosphorylation of OmpR is considered to occur at one or more aspartic acid residues (Asp-11, Asp-12 and/or Asp-55) that are highly conserved among the effector proteins. In this report we biochemically characterized the aspartic acid residue(s) in OmpR that were phosphorylated by EnvZ. Reduction of aspartyl phosphate residues in the amino-terminal domain of OmpR with [³H]-NaBH₄ indicated that Asp-55 was a primary site of modification. We further studied the role of the highly conserved aspartate residues by creating OmpR mutants having aspartate to alanine substitutions at positions 11 (D11A), 12 (D12A) and 55 (D55A). Studies of ompF and ompC expression as well as in vivo and in vitro phosphorylation experiments also demonstrated that while Asp-55 is the primary phosphate acceptor site in OmpR, Asp-11 may also serve as a phosphorylation site, particularly in the absence of Asp-55.     

Evidence for a common metal ion-dependent transition in the 4-carboxyglutamic acid domains of several vitamin K-dependent proteins

The murine monoclonal antibody H-11 binds a conserved epitope found at the amino terminal of the vitamin K-dependent blood proteins prothrombin, factors VII and X, and protein C. The sequence of polypeptide recognized by antibody H-11 contains 2 residues of gamma-carboxyglutamic acid, and binding of the antibody is inhibited by divalent metal ions. By using a solid-phase immunoassay with 125I-labeled antibody and immobilized vitamin K-dependent protein, binding of the antibody to the vitamin K-dependent proteins was inhibited by increasing concentrations of calcium, manganese, and magnesium ion. The transition midpoints for antibody binding were in the millimolar concentration range and were different for each metal ion. In general, the transition midpoints were lowest for manganese ion, intermediate for calcium ion, and highest for magnesium ion. Antibody H-11 bound specifically to a synthetic peptide corresponding to residues 1-12 of human prothrombin that was synthesized as the gamma-carboxyglutamic acid-containing derivative. Binding of the antibody to the peptide was not inhibited by calcium ion. These data suggest that inhibition of antibody H-11 binding by divalent metal ions is not due simply to neutralization of negative charge by Ca2+. This transition which is conserved in vitamin K-dependent proteins containing the H-11 antigenic site is likely due to a structural transition of the amino-terminal polypeptide possibly from a random (accessible) to ordered (inaccessible) structure.     

Molecular dynamics study of Ca(2+) binding loop variants of silver hake parvalbumin with aspartic acid at the "gateway" position

The helix-loop-helix (i.e., EF-hand) Ca(2+) 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 beta-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 Ca(2+) 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.