Elucidating the mechanisms of cooperative calcium-calmodulin interactions: a structural systems biology approach

BACKGROUND: Calmodulin is an important multifunctional molecule that regulates the activities of a large number of proteins in the cell. Calcium binding induces conformational transitions in calmodulin that make it specifically active to particular target proteins. The precise mechanisms underlying calcium binding to calmodulin are still, however, quite poorly understood. RESULTS: In this study, we adopt a structural systems biology approach and develop a mathematical model to investigate various types of cooperative calcium-calmodulin interactions. We compare the predictions of our analysis with physiological dose-response curves taken from the literature, in order to provide a quantitative comparison of the effects of different mechanisms of cooperativity on calcium-calmodulin interactions. The results of our analysis reduce the gap between current understanding of intracellular calmodulin function at the structural level and physiological calcium-dependent calmodulin target activation experiments. CONCLUSION: Our model predicts that the specificity and selectivity of CaM target regulation is likely to be due to the following factors: variations in the target-specific Ca2+ dissociation and cooperatively effected dissociation constants, and variations in the number of Ca2+ ions required to bind CaM for target activation

Protein multi-scale organization through graph partitioning and robustness analysis: Application to the myosin-myosin light chain interaction

Despite the recognized importance of the multi-scale spatio-temporal organization of proteins, most computational tools can only access a limited spectrum of time and spatial scales, thereby ignoring the effects on protein behavior of the intricate coupling between the different scales. Starting from a physico-chemical atomistic network of interactions that encodes the structure of the protein, we introduce a methodology based on multi-scale graph partitioning that can uncover partitions and levels of organization of proteins that span the whole range of scales, revealing biological features occurring at different levels of organization and tracking their effect across scales. Additionally, we introduce a measure of robustness to quantify the relevance of the partitions through the generation of biochemically-motivated surrogate random graph models. We apply the method to four distinct conformations of myosin tail interacting protein, a protein from the molecular motor of the malaria parasite, and study properties that have been experimentally addressed such as the closing mechanism, the presence of conserved clusters, and the identification through computational mutational analysis of key residues for binding.Comment: 13 pages, 7 Postscript figure

Hydrogen-activation mechanism of [Fe] hydrogenase revealed by multi-scale modeling

When investigating the mode of hydrogen activation by [Fe] hydrogenases, not only the chemical reactivity at the active site is of importance but also the large-scale conformational change between the so-called open and closed conformations, which leads to a special spatial arrangement of substrate and iron cofactor. To study H2 activation, a complete model of the solvated and cofactor-bound enzyme in complex with the substrate methenyl-H4MPT+ was constructed. Both the closed and open conformations were simulated with classical molecular dynamics on the 100 ns time scale. Quantum-mechanics/molecular-mechanics calculations on snapshots then revealed the features of the active site that enable the facile H2 cleavage. The hydroxyl group of the pyridinol ligand can easily be deprotonated. With the deprotonated hydroxyl group and the structural arrangement in the closed conformation, H2 coordinated to the Fe center is subject to an ionic and orbital push-pull effect and can be rapidly cleaved with a concerted hydride transfer to methenyl-H4MPT+. An intermediary hydride species is not formed

Structural alphabets derived from attractors in conformational space

Background: The hierarchical and partially redundant nature of protein structures justifies the definition of frequently occurring conformations of short fragments as 'states'. Collections of selected representatives for these states define Structural Alphabets, describing the most typical local conformations within protein structures. These alphabets form a bridge between the string-oriented methods of sequence analysis and the coordinate-oriented methods of protein structure analysis.Results: A Structural Alphabet has been derived by clustering all four-residue fragments of a high-resolution subset of the protein data bank and extracting the high-density states as representative conformational states. Each fragment is uniquely defined by a set of three independent angles corresponding to its degrees of freedom, capturing in simple and intuitive terms the properties of the conformational space. The fragments of the Structural Alphabet are equivalent to the conformational attractors and therefore yield a most informative encoding of proteins. Proteins can be reconstructed within the experimental uncertainty in structure determination and ensembles of structures can be encoded with accuracy and robustness.Conclusions: The density-based Structural Alphabet provides a novel tool to describe local conformations and it is specifically suitable for application in studies of protein dynamics. © 2010 Pandini et al; licensee BioMed Central Ltd

MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations

Molecular dynamics (MD) simulations can model the interactions between macromolecules with high spatiotemporal resolution but at a high computational cost. By combining high-throughput MD with Markov state models (MSMs), it is now possible to obtain long-timescale behavior of small to intermediate biomolecules and complexes. To model the interactions of many molecules at large lengthscales, particle-based reaction-diffusion (RD) simulations are more suitable but lack molecular detail. Thus, coupling MSMs and RD simulations (MSM/RD) would be highly desirable, as they could efficiently produce simulations at large time- and lengthscales, while still conserving the characteristic features of the interactions observed at atomic detail. While such a coupling seems straightforward, fundamental questions are still open: Which definition of MSM states is suitable? Which protocol to merge and split RD particles in an association/dissociation reaction will conserve the correct bimolecular kinetics and thermodynamics? In this paper, we make the first step towards MSM/RD by laying out a general theory of coupling and proposing a first implementation for association/dissociation of a protein with a small ligand (A + B C). Applications on a toy model and CO diffusion into the heme cavity of myoglobin are reported

High-Resolution Optical Studies on C-Phycocyanin via Photochemical Hole Burning

We have shown that both the native C-phycocyanin and its corresponding free biline chromophore undergo reversible, low-temperature photochemistry. We attribute this photochemistry to reversible proton-transfer processes and utilize the observed photoreaction for photochemical hole burning (PHB). Using narrow-band PHB experiments, we have been able to perform high-resolution optical studies and show that the protein-chromophore assembly forms a very rigid structure. The results lead to the conclusion that the light-induced proton transfer occurs most probably in the triplet state