Symmetrical stabilization of bound Ca2+ ions in a cooperative pair of EF-hands through hydrogen bonding of coordinating water molecules in calbindin D(9k)

Charge dependent retardation of amyloid β aggregation by hydrophilic proteins

The aggregation of amyloid β peptides (Aβ) into amyloid fibrils is implicated in the pathology of Alzheimer's disease. In light of the increasing number of proteins reported to retard Aβ fibril formation, we investigated the influence of small hydrophilic model proteins of different charge on Aβ aggregation kinetics and their interaction with Aβ. We followed the amyloid fibril formation of Aβ40 and Aβ42 using thioflavin T fluorescence in the presence of six charge variants of calbindin D9k and single-chain monellin. The formation of fibrils was verified with transmission electron microscopy. We observe retardation of the aggregation process from proteins with net charge +8, +2, -2, and -4, whereas no effect is observed for proteins with net charge of -6 and -8. The single-chain monellin mutant with the highest net charge, scMN+8, has the largest retarding effect on the amyloid fibril formation process, which is noticeably delayed at as low as a 0.01:1 scMN+8 to Aβ40 molar ratio. scMN+8 is also the mutant with the fastest association to Aβ40 as detected by surface plasmon resonance, although all retarding variants of calbindin D9k and single-chain monellin bind to Aβ40.

Calculation of Free Energy Landscape in Multi-Dimensions with Hamiltonian-Exchange Umbrella Sampling on Petascale Supercomputer

An extremely scalable computational strategy is described for calculations of the potential of mean force (PMF) in multidimensions on massively distributed supercomputers. The approach involves coupling thousands of umbrella sampling (US) simulation windows distributed to cover the space of order parameters with a Hamiltonian molecular dynamics replica-exchange (H-REMD) algorithm to enhance the sampling of each simulation. In the present application, US/H-REMD is carried out in a two-dimensional (2D) space and exchanges are attempted alternatively along the two axes corresponding to the two order parameters. The US/H-REMD strategy is implemented on the basis of parallel/parallel multiple copy protocol at the MPI level, and therefore can fully exploit computing power of large-scale supercomputers. Here the novel technique is illustrated using the leadership supercomputer IBM Blue Gene/P with an application to a typical biomolecular calculation of general interest, namely the binding of calcium ions to the small protein Calbindin D9k. The free energy landscape associated with two order parameters, the distance between the ion and its binding pocket and the root-mean-square deviation (rmsd) of the binding pocket relative the crystal structure, was calculated using the US/H-REMD method. The results are then used to estimate the absolute binding free energy of calcium ion to Calbindin D9k. The tests demonstrate that the 2D US/H-REMD scheme greatly accelerates the configurational sampling of the binding pocket, thereby improving the convergence of the potential of mean force calculation.

Coupling of ligand binding and dimerization of helix-loop-helix peptides: spectroscopic and sedimentation analyses of calbindin D9k EF-hands

Isolated Ca2+-binding EF-hand peptides have a tendency to dimerize. This study is an attempt to account for the coupled equilibria of Ca2+-binding and peptide association for two EF-hands with strikingly different loop sequence and net charge. We have studied each of the two separate EF-hand fragments from calbindin D9k. A series of Ca2+-titrations at different peptide concentrations were monitored by CD and fluorescence spectroscopy. All data were fitted simultaneously to both a complete model of all possible equilibrium intermediates and a reduced model not including dimerization in the absence of Ca2+. Analytical ultracentrifugation shows that the peptides may occur as monomers or dimers depending on the solution conditions. Our results show strikingly different behavior for the two EF-hands. The fragment containing the N-terminal EF-hand shows a strong tendency to dimerize in the Ca2+-bound state. The average Ca2+-affinity is 3.5 orders of magnitude lower than for the intact protein. We observe a large apparent cooperativity of Ca2+ binding for the overall process from Ca2+-free monomer to fully loaded dimer, showing that a Ca2+-free EF-hand folds upon dimerization to a Ca2+-bound EF-hand, thereby presenting a preformed binding site to the second Ca2+-ion. The C-terminal EF-hand shows a much smaller tendency to dimerize, which may be related to its larger net negative charge. In spite of the differences in dimerization behavior, the Ca2+ affinities of both EF-hand fragments are similar and in the range lgK = 4.6-5.3.