Retention of conformational entropy upon calmodulin binding to target peptides is driven by transient salt bridges

Calcium mediated static and dynamic allostery in S100A12: Implications for target recognition by S100 proteins

Structure and functions of S100 proteins are regulated by two distinct calcium binding EF hand motifs. In this work, we used solution-state NMR spectroscopy to investigate the cooperativity between the two calcium binding sites and map the allosteric changes at the target binding site. To parse the contribution of the individual calcium binding events, variants of S100A12 were designed to selectively bind calcium to either the EF-I (N63A) or EF-II (E31A) loop, respectively. Detailed analysis of the backbone chemical shifts for wildtype protein and its mutants indicates that calcium binding to the canonical EF-II loop is the principal trigger for the conformational switch between 'closed' apo to the 'open' Ca2+ -bound conformation of the protein. Elimination of binding in S100-specific EF-I loop has limited impact on the calcium binding affinity of the EF-II loop and the concomitant structural rearrangement. In contrast, deletion of binding in the EF-II loop significantly attenuates calcium affinity in the EF-I loop and the structure adopts a 'closed' apo-like conformation. Analysis of experimental amide nitrogen (15 N) relaxation rates (R1 , R2 , and 15 N-{1 H} NOE) and molecular dynamics (MD) simulations demonstrate that the calcium bound state is relatively floppy with pico-nanosecond motions induced in functionally relevant domains responsible for target recognition such as the hinge domain and the C-terminal residues. Experimental relaxation studies combined with MD simulations show that while calcium binding in the EF-I loop alone does not induce significant motions in the polypeptide chain, EF-I regulates fluctuations in the polypeptide in the presence of bound calcium in the EF-II loop. These results offer novel insights into the dynamic regulation of target recognition by calcium binding and unravels the role of cooperativity between the two calcium binding events in S100A12.

Markov state modelling reveals heterogeneous drug-inhibition mechanism of Calmodulin

Calmodulin (CaM) is a calcium sensor which binds and regulates a wide range of target-proteins. This implicitly enables the concentration of calcium to influence many downstream physiological responses, including muscle contraction, learning and depression. The antipsychotic drug trifluoperazine (TFP) is a known CaM inhibitor. By binding to various sites, TFP prevents CaM from associating to target-proteins. However, the molecular and state-dependent mechanisms behind CaM inhibition by drugs such as TFP are largely unknown. Here, we build a Markov state model (MSM) from adaptively sampled molecular dynamics simulations and reveal the structural and dynamical features behind the inhibitory mechanism of TFP-binding to the C-terminal domain of CaM. We specifically identify three major TFP binding-modes from the MSM macrostates, and distinguish their effect on CaM conformation by using a systematic analysis protocol based on biophysical descriptors and tools from machine learning. The results show that depending on the binding orientation, TFP effectively stabilizes features of the calcium-unbound CaM, either affecting the CaM hydrophobic binding pocket, the calcium binding sites or the secondary structure content in the bound domain. The conclusions drawn from this work may in the future serve to formulate a complete model of pharmacological modulation of CaM, which furthers our understanding of how these drugs affect signaling pathways as well as associated diseases.

Calmodulin (CaM) is a calcium-sensing protein which makes other proteins dependent on the surrounding calcium concentration by binding to these proteins. Such protein-protein interactions with CaM are vital for calcium to control many physiological pathways within the cell. The antipsychotic drug trifluoperazine (TFP) inhibits CaM’s ability to bind and regulate other proteins. Here, we use molecular dynamics simulations together with Markov state modeling and machine learning to understand the structural and dynamical features by which TFP bound to the one domain of CaM prevents association to other proteins. We find that TFP encourages CaM to adopt a conformation that is like the one stabilized in absence of calcium: depending on the binding orientation of TFP, the drug indeed either affects the CaM hydrophobic binding pocket, the calcium binding sites or the secondary structure content in the domain. Understanding TFP binding is a first step towards designing better drugs targeting CaM.

Assessing the Role of Calmodulin's Linker Flexibility in Target Binding

Calmodulin (CaM) is a highly-expressed Ca2+ binding protein known to bind hundreds of protein targets. Its binding selectivity to many of these targets is partially attributed to the protein’s flexible alpha helical linker that connects its N- and C-domains. It is not well established how its linker mediates CaM’s binding to regulatory targets yet. Insights into this would be invaluable to understanding its regulation of diverse cellular signaling pathways. Therefore, we utilized Martini coarse-grained (CG) molecular dynamics simulations to probe CaM/target assembly for a model system: CaM binding to the calcineurin (CaN) regulatory domain. The simulations were conducted assuming a ‘wild-type’ calmodulin with normal flexibility of its linker, as well as a labile, highly-flexible linker variant to emulate structural changes that could be induced, for instance, by post-translational modifications. For the wild-type model, 98% of the 600 simulations across three ionic strengths adopted a bound complex within 2 μs of simulation time; of these, 1.7% sampled the fully-bound state observed in the experimentally-determined crystallographic structure. By calculating the mean-first-passage-time for these simulations, we estimated the association rate to be ka= 8.7 × 108 M−1 s−1, which is similar to the diffusion-limited, experimentally-determined rate of 2.2 × 108 M−1 s−1. Furthermore, our simulations recapitulated its well-known inverse relationship between the association rate and the solution ionic strength. In contrast, although over 97% of the labile linker simulations formed tightly-bound complexes, only 0.3% achieved the fully-bound configuration. This effect appears to stem from a difference in the ensembles of extended and collapsed states which are controlled by the linker flexibility. Therefore, our simulations suggest that variations in the CaM linker’s propensity for alpha helical secondary structure can modulate the kinetics of target binding.

Effect of Ca2+ on the promiscuous target-protein binding of calmodulin

Calmodulin (CaM) is a calcium sensing protein that regulates the function of a large number of proteins, thus playing a crucial part in many cell signaling pathways. CaM has the ability to bind more than 300 different target peptides in a Ca2+-dependent manner, mainly through the exposure of hydrophobic residues. How CaM can bind a large number of targets while retaining some selectivity is a fascinating open question. Here, we explore the mechanism of CaM selective promiscuity for selected target proteins. Analyzing enhanced sampling molecular dynamics simulations of Ca2+-bound and Ca2+-free CaM via spectral clustering has allowed us to identify distinct conformational states, characterized by interhelical angles, secondary structure determinants and the solvent exposure of specific residues. We searched for indicators of conformational selection by mapping solvent exposure of residues in these conformational states to contacts in structures of CaM/target peptide complexes. We thereby identified CaM states involved in various binding classes arranged along a depth binding gradient. Binding Ca2+ modifies the accessible hydrophobic surface of the two lobes and allows for deeper binding. Apo CaM indeed shows shallow binding involving predominantly polar and charged residues. Furthermore, binding to the C-terminal lobe of CaM appears selective and involves specific conformational states that can facilitate deep binding to target proteins, while binding to the N-terminal lobe appears to happen through a more flexible mechanism. Thus the long-ranged electrostatic interactions of the charged residues of the N-terminal lobe of CaM may initiate binding, while the short-ranged interactions of hydrophobic residues in the C-terminal lobe of CaM may account for selectivity. This work furthers our understanding of the mechanism of CaM binding and selectivity to different target proteins and paves the way towards a comprehensive model of CaM selectivity.

Computer Simulations of Intrinsically Disordered Proteins

The investigation of intrinsically disordered proteins (IDPs) is a new frontier in structural and molecular biology that requires a new paradigm to connect structural disorder to function. Molecular dynamics simulations and statistical thermodynamics potentially offer ideal tools for atomic-level characterizations and thermodynamic descriptions of this fascinating class of proteins that will complement experimental studies. However, IDPs display sensitivity to inaccuracies in the underlying molecular mechanics force fields. Thus, achieving an accurate structural characterization of IDPs via simulations is a challenge. It is also daunting to perform a configuration-space integration over heterogeneous structural ensembles sampled by IDPs to extract, in particular, protein configurational entropy. In this review, we summarize recent efforts devoted to the development of force fields and the critical evaluations of their performance when applied to IDPs. We also survey recent advances in computational methods for protein configurational entropy that aim to provide a thermodynamic link between structural disorder and protein activity.

Multilevel Changes in Protein Dynamics upon Complex Formation of the Calcium-Loaded S100A4 with a Nonmuscle Myosin IIA Tail Fragment

Dysregulation of Ca²⁺ -binding S100 proteins plays important role in various diseases. The asymmetric complex of Ca²⁺ -bound S100A4 with nonmuscle myosin IIA has high stability and highly increased Ca²⁺ affinity. Here we investigated the possible causes of this allosteric effect by NMR spectroscopy. Chemical shift-based secondary-structure analysis did not show substantial changes for the complex. Backbone dynamics revealed slow-timescale local motions in the H1 helices of homodimeric S100A4; these were less pronounced in the complex form and might be accompanied by an increase in dimer stability. Different mobilities in the Ca²⁺ -coordinating EF-hand sites indicate that they communicate by an allosteric mechanism operating through changes in protein dynamics; this must be responsible for the elevated Ca²⁺ affinity. These multilevel changes in protein dynamics as conformational adaptation allow S100A4 fine-tuning of its protein-protein interactions inside the cell during Ca²⁺ signaling.

Molecular dynamics simulations provide insights into the substrate specificity of FAOX family members

Enzymatic assays based on Fructosyl Amino Acid Oxidases (FAOX) represent a potential, rapid and economical strategy to measure glycated hemoglobin (HbA1c), which is in turn a reliable method to monitor the insurgence and the development of diabetes mellitus.

Effect of graphene oxide on the conformational transitions of amyloid beta peptide: A molecular dynamics simulation study

The interactions between nanomaterials (NMs) and amyloid proteins are central to the nanotechnology-based diagnostics and therapy in neurodegenerative disorders such as Alzheimer's and Parkinson's. Graphene oxide (GO) and its derivatives have shown to modulate the aggregation pattern of disease causing amyloid beta (Aβ) peptide. However, the mechanism is still not well understood. Using molecular dynamics simulations, the effect of graphene oxide (GO) and reduced graphene oxide (rGO) having carbon:oxygen ratio of 4:1 and 10:1, respectively, on the conformational transitions (alpha-helix to beta-sheet) and the dynamics of the peptide was investigated. GO and rGO decreased the beta-strand propensity of amino acid residues in Aβ. The peptide displayed different modes of adsorption on GO and rGO. The adsorption on GO was dominated by electrostatic interactions, whereas on rGO, both van der Waals and electrostatic interactions contributed in the adsorption of the peptide. Our study revealed that the slight increase in the hydrophobic patches on rGO made it more effective inhibitor of conformational transitions in the peptide. Alpha helix-beta sheet transition in Aβ peptide could be one of the plausible mechanism by which graphene oxide may inhibit amyloid fibrillation.

Conformational frustration in calmodulin-target recognition

Calmodulin (CaM) is a primary calcium (Ca(2+) )-signaling protein that specifically recognizes and activates highly diverse target proteins. We explored the molecular basis of target recognition of CaM with peptides representing the CaM-binding domains from two Ca(2+) -CaM-dependent kinases, CaMKI and CaMKII, by employing experimentally constrained molecular simulations. Detailed binding route analysis revealed that the two CaM target peptides, although similar in length and net charge, follow distinct routes that lead to a higher binding frustration in the CaM-CaMKII complex than in the CaM-CaMKI complex. We discovered that the molecular origin of the binding frustration is caused by intermolecular contacts formed with the C-domain of CaM that need to be broken before the formation of intermolecular contacts with the N-domain of CaM. We argue that the binding frustration is important for determining the kinetics of the recognition process of proteins involving large structural fluctuations.

Distinguishing unfolding and functional conformational transitions of calmodulin using ultraviolet resonance Raman spectroscopy

Calmodulin (CaM) is a ubiquitous moderator protein for calcium signaling in all eukaryotic cells. This small calcium-binding protein exhibits a broad range of structural transitions, including domain opening and folding-unfolding, that allow it to recognize a wide variety of binding partners in vivo. While the static structures of CaM associated with its various binding activities are fairly well-known, it has been challenging to examine the dynamics of transition between these structures in real-time, due to a lack of suitable spectroscopic probes of CaM structure. In this article, we examine the potential of ultraviolet resonance Raman (UVRR) spectroscopy for clarifying the nature of structural transitions in CaM. We find that the UVRR spectral change (with 229 nm excitation) due to thermal unfolding of CaM is qualitatively different from that associated with opening of the C-terminal domain in response to Ca(2+) binding. This spectral difference is entirely due to differences in tertiary contacts at the interdomain tyrosine residue Tyr138, toward which other spectroscopic methods are not sensitive. We conclude that UVRR is ideally suited to identifying the different types of structural transitions in CaM and other proteins with conformation-sensitive tyrosine residues, opening a path to time-resolved studies of CaM dynamics using Raman spectroscopy.

Modeling the interactions between MC2R and ACTH models from human

Melanocortin system is composed of four peptide hormones namely α-, β-, -γ, and adrenocorticotropic hormone (ACTH), derived from post-translational cleavage of a polypeptide precursor 'proopiomelanocortin (POMC).' Among these hormones, ACTH, a 38 amino acid residue peptide fragment is an important hormone as it is involved in steroid secretion. In addition to this, to cite a few, this hormone is also known to induce variety of other effects, such as alterations in motor/sexual behavior, improvement in memory, and anti-inflammatory effects. To date, five melanocortin receptors (MC1R-MC5R) have been characterized with tissue-specific expression patterns and different binding affinities for each of the melanocortin hormones to regulate various biological functions. In the present work, three-dimensional (3D) models of MC2R and ACTH from human have been predicted, followed by docking and molecular dynamics simulation. While the 3D model of MC2R receptor has been predicted through threading approach, structure of ACTH was built based on ab initio technique. The MC2R model was later successfully docked onto the ACTH structure. Molecular dynamics (MD) simulation for 20 ns was used to compute the binding free energy of MC2R with ACTH model under implicit solvent conditions.

Effect of subdomain interactions on methyl group dynamics in the hydrophobic core of villin headpiece protein

Thermostable villin headpiece protein (HP67) consists of the N-terminal subdomain (residues 10-41) and the autonomously folding C-terminal subdomain (residues 42-76) which pack against each other to form a structure with a unified hydrophobic core. The X-ray structures of the isolated C-terminal subdomain (HP36) and its counterpart in HP67 are very similar for the hydrophobic core residues. However, fine rearrangements of the free energy landscape are expected to occur because of the interactions between the two subdomains. We detect and characterize these changes by comparing the µs-ms time scale dynamics of the methyl-bearing side chains in isolated HP36 and in HP67. Specifically, we probe three hydrophobic side chains at the interface of the two subdomains (L42, V50, and L75) as well as at two residues far from the interface (L61 and L69). Solid-state deuteron NMR techniques are combined with computational modeling for the detailed characterization of motional modes in terms of their kinetic and thermodynamic parameters. The effect of interdomain interactions on side chain dynamics is seen for all residues but L75. Thus, changes in dynamics because of subdomain interactions are not confined to the site of perturbation. One of the main results is a two- to threefold increase in the value of the activation energies for the rotameric mode of motions in HP67 compared with HP36. Detailed analysis of configurational entropies and heat capacities complement the kinetic view of the degree of the disorder in the folded state.