Towards predicting Ca2+-binding sites with different coordination numbers in proteins with atomic resolution

Identification of the modulatory Ca2+-binding sites of acid-sensing ion channel 1a

Acid-sensing ion channels (ASICs) are neuronal Na+-permeable ion channels activated by extracellular acidification. ASICs are involved in learning, fear sensing, pain sensation and neurodegeneration. Increasing the extracellular Ca2+ concentration decreases the H+ sensitivity of ASIC1a, suggesting a competition for binding sites between H+ and Ca2+ ions. Here, we predicted candidate residues for Ca2+ binding on ASIC1a, based on available structural information and our molecular dynamics simulations. With functional measurements, we identified several residues in cavities previously associated with pH-dependent gating, whose mutation reduced the modulation by extracellular Ca2+ of the ASIC1a pH dependence of activation and desensitization. This occurred likely owing to a disruption of Ca2+ binding. Our results link one of the two predicted Ca2+-binding sites in each ASIC1a acidic pocket to the modulation of channel activation. Mg2+ regulates ASICs in a similar way as does Ca2+. We show that Mg2+ shares some of the binding sites with Ca2+. Finally, we provide evidence that some of the ASIC1a Ca2+-binding sites are functionally conserved in the splice variant ASIC1b. Our identification of divalent cation-binding sites in ASIC1a shows how Ca2+ affects ASIC1a gating, elucidating a regulatory mechanism present in many ion channels.

The structural basis of tRNA recognition by arginyl-tRNA-protein transferase

Arginyl-tRNA-protein transferase 1 (ATE1) is a master regulator of protein homeostasis, stress response, cytoskeleton maintenance, and cell migration. The diverse functions of ATE1 arise from its unique enzymatic activity to covalently attach an arginine onto its protein substrates in a tRNA-dependent manner. However, how ATE1 (and other aminoacyl-tRNA transferases) hijacks tRNA from the highly efficient ribosomal protein synthesis pathways and catalyzes the arginylation reaction remains a mystery. Here, we describe the three-dimensional structures of Saccharomyces cerevisiae ATE1 with and without its tRNA cofactor. Importantly, the putative substrate binding domain of ATE1 adopts a previously uncharacterized fold that contains an atypical zinc-binding site critical for ATE1 stability and function. The unique recognition of tRNA^Arg by ATE1 is coordinated through interactions with the major groove of the acceptor arm of tRNA. Binding of tRNA induces conformational changes in ATE1 that helps explain the mechanism of substrate arginylation.

Calcium regulates acid-sensing ion channel 3 activation by competing with protons in the channel pore and at an allosteric binding site

The extracellular Ca²⁺ concentration changes locally under certain physiological and pathological conditions. Such variations affect the function of ion channels of the nervous system and consequently also neuronal signalling. We investigated here the mechanisms by which Ca²⁺ controls the activity of acid-sensing ion channel (ASIC) 3. ASICs are neuronal, H⁺-gated Na⁺ channels involved in several physiological and pathological processes, including the expression of fear, learning, pain sensation and neurodegeneration after ischaemic stroke. It was previously shown that Ca²⁺ negatively modulates the ASIC pH dependence. While protons are default activators of ASIC3, this channel can also be activated at pH7.4 by the removal of the extracellular Ca²⁺. Two previous studies concluded that low pH opens ASIC3 by displacing Ca²⁺ ions that block the channel pore at physiological pH. We show here that an acidic residue, distant from the pore, together with pore residues, controls the modulation of ASIC3 by Ca²⁺. Our study identifies a new regulatory site in ASIC3 and demonstrates that ASIC3 activation involves an allosteric mechanism together with Ca²⁺ unbinding from the channel pore. We provide a molecular analysis of a regulatory mechanism found in many ion channels.

Insight into Calcium-Binding Motifs of Intrinsically Disordered Proteins

Motifs within proteins help us categorize their functions. Intrinsically disordered proteins (IDPs) are rich in short linear motifs, conferring them many different roles. IDPs are also frequently highly charged and, therefore, likely to interact with ions. Canonical calcium-binding motifs, such as the EF-hand, often rely on the formation of stabilizing flanking helices, which are a key characteristic of folded proteins, but are absent in IDPs. In this study, we probe the existence of a calcium-binding motif relevant to IDPs. Upon screening several carefully selected IDPs using NMR spectroscopy supplemented with affinity quantification by colorimetric assays, we found calcium-binding motifs in IDPs which could be categorized into at least two groups-an Excalibur-like motif, sequentially similar to the EF-hand loop, and a condensed-charge motif carrying repetitive negative charges. The motifs show an affinity for calcium typically in the ~100 μM range relevant to regulatory functions and, while calcium binding to the condensed-charge motif had little effect on the overall compaction of the IDP chain, calcium binding to Excalibur-like motifs resulted in changes in compaction. Thus, calcium binding to IDPs may serve various structural and functional roles that have previously been underreported.

Structural Changes beyond the EF-Hand Contribute to Apparent Calcium Binding Affinities: Insights from Parvalbumins

Members of the parvalbumin (PV) family of calcium (Ca²⁺) binding proteins (CBPs) share a relatively high level of sequence similarity. However, their Ca²⁺ affinities and selectivities against competing ions like Mg²⁺ can widely vary. We conducted molecular dynamics simulations of several α-parvalbumin (αPV) constructs with micromolar to nanomolar Ca²⁺ affinities to identify structural and dynamic features that contribute to their binding of ions. Specifically, we examined a D94S/G98E construct with a lower Ca²⁺ affinity (≈-18 kcal/mol) relative to the wild type (WT) (≈-22 kcal/mol) and an S55D/E59D variant with enhanced affinity (≈-24 kcal/mol). Additionally, we also examined the binding of Mg²⁺ to these isoforms, which is much weaker than Ca²⁺. We used mean spherical approximation (MSA) theory to evaluate ion binding thermodynamics within the proteins' EF-hand domains to account for the impact of ions' finite sizes and the surrounding electrolyte composition. While the MSA scores differentiated Mg²⁺ from Ca²⁺, they did not indicate that Ca²⁺ binding affinities at the binding loop differed between the PV isoforms. Instead, molecular mechanics generalized Born surface area (MM/GBSA) approximation energies, which we used to quantify the thermodynamic cost of structural rearrangement of the proteins upon binding ions, indicated that S55D/E59D αPV favored Ca²⁺ binding by -20 kcal/mol relative to WT versus 30 kcal/mol for D94S/G98E αPV. Meanwhile, Mg²⁺ binding was favored for the S55D/E59D αPV and D94S/G98E αPV variants by -18.32 and -1.65 kcal/mol, respectively. These energies implicate significant contributions to ion binding beyond oxygen coordination at the binding loop, which stemmed from changes in α-helicity, β-sheet character, and hydrogen bonding. Hence, Ca²⁺ affinity and selectivity against Mg²⁺ are emergent properties stemming from both local effects within the proteins' ion binding sites as well as non-local contributions elsewhere. Our findings broaden our understanding of the molecular bases governing αPV ion binding that are likely shared by members of the broad family of CBPs.

Structural Aspects and Prediction of Calmodulin-Binding Proteins

Calmodulin (CaM) is an important intracellular protein that binds Ca²⁺ and functions as a critical second messenger involved in numerous biological activities through extensive interactions with proteins and peptides. CaM’s ability to adapt to binding targets with different structures is related to the flexible central helix separating the N- and C-terminal lobes, which allows for conformational changes between extended and collapsed forms of the protein. CaM-binding targets are most often identified using prediction algorithms that utilize sequence and structural data to predict regions of peptides and proteins that can interact with CaM. In this review, we provide an overview of different CaM-binding proteins, the motifs through which they interact with CaM, and shared properties that make them good binding partners for CaM. Additionally, we discuss the historical and current methods for predicting CaM binding, and the similarities and differences between these methods and their relative success at prediction. As new CaM-binding proteins are identified and classified, we will gain a broader understanding of the biological processes regulated through changes in Ca²⁺ concentration through interactions with CaM.

Structural Mechanism of Cooperative Regulation of Calcium-Sensing Receptor-Mediated Cellular Signaling

Calcaium sensing receptors (CaSRs) play a central role in regulating extracellular calcium (Ca²⁺) homeostasis and many (patho)physiological processes. This regulation is primarily orchestrated in response to extracellular stimuli via the extracellular domain (ECD). This paper first reviews the modeled structure of the CaSR ECD and the prediction and investigation of the Ca²⁺ and amino acid binding sites. Several recently solved X-ray structures are then compared to support a proposed CaSR activation model involving functional cooperativity. The review also discusses recent implications for drug development. These studies provide new insights into the molecular basis of diseases and the design of therapeutic agents that target CaSR and other family C G protein-coupled receptors (cGPCRs).

Irritant-evoked activation and calcium modulation of the TRPA1 receptor

The TRPA1 ion channel is expressed by primary afferent nerve fibers, where it functions as a low threshold sensor for structurally diverse electrophilic irritants ranging from small volatile environmental toxicants to endogenous algogenic lipids¹. TRPA1 is also a ‘receptor-operated’ channel whose activation downstream of metabotropic receptors elicits inflammatory pain or itch, making it an attractive target for novel analgesic therapies². However, we lack mechanistic insight into how TRPA1 recognizes and responds to electrophiles or cytoplasmic second messengers. Here, we show that electrophiles act through a two-step process in which modification of a highly reactive cysteine (C621) promotes reorientation of a cytoplasmic loop to enhance nucleophilicity and modification of a nearby cysteine (C665), thereby stabilizing the loop in an activating configuration. These actions modulate two restrictions controlling ion permeation, including widening of the selectivity filter to enhance calcium permeability and opening of a canonical gate at the cytoplasmic end of the pore. We propose a model to explain functional coupling between electrophile action and these control points. We also characterize a calcium binding pocket that is remarkably conserved across TRP channel subtypes and accounts for all aspects of calcium-dependent TRPA1 regulation, including potentiation, desensitization, and activation by metabotropic receptors. These findings provide a structural framework for understanding how a broad-spectrum irritant receptor is controlled by endogenous and exogenous agents that elicit or exacerbate pain and itch.

International Union of Basic and Clinical Pharmacology. CVIII. Calcium-Sensing Receptor Nomenclature, Pharmacology, and Function

The calcium-sensing receptor (CaSR) is a class C G protein-coupled receptor that responds to multiple endogenous agonists and allosteric modulators, including divalent and trivalent cations, L-amino acids, γ-glutamyl peptides, polyamines, polycationic peptides, and protons. The CaSR plays a critical role in extracellular calcium (Ca2+ o) homeostasis, as demonstrated by the many naturally occurring mutations in the CaSR or its signaling partners that cause Ca2+ o homeostasis disorders. However, CaSR tissue expression in mammals is broad and includes tissues unrelated to Ca2+ o homeostasis, in which it, for example, regulates the secretion of digestive hormones, airway constriction, cardiovascular effects, cellular differentiation, and proliferation. Thus, although the CaSR is targeted clinically by the positive allosteric modulators (PAMs) cinacalcet, evocalcet, and etelcalcetide in hyperparathyroidism, it is also a putative therapeutic target in diabetes, asthma, cardiovascular disease, and cancer. The CaSR is somewhat unique in possessing multiple ligand binding sites, including at least five putative sites for the "orthosteric" agonist Ca2+ o, an allosteric site for endogenous L-amino acids, two further allosteric sites for small molecules and the peptide PAM, etelcalcetide, and additional sites for other cations and anions. The CaSR is promiscuous in its G protein-coupling preferences, and signals via Gq/11, Gi/o, potentially G12/13, and even Gs in some cell types. Not surprisingly, the CaSR is subject to biased agonism, in which distinct ligands preferentially stimulate a subset of the CaSR's possible signaling responses, to the exclusion of others. The CaSR thus serves as a model receptor to study natural bias and allostery. SIGNIFICANCE STATEMENT: The calcium-sensing receptor (CaSR) is a complex G protein-coupled receptor that possesses multiple orthosteric and allosteric binding sites, is subject to biased signaling via several different G proteins, and has numerous (patho)physiological roles. Understanding the complexities of CaSR structure, function, and biology will aid future drug discovery efforts seeking to target this receptor for a diversity of diseases. This review summarizes what is known to date regarding key structural, pharmacological, and physiological features of the CaSR.

Investigation of bound and unbound phosphoserine phosphatase conformations through elastic network models and molecular dynamics simulations

The human phosphoserine phosphatase (hPSP) catalyses the last step in the biosynthesis of L-serine. It involves conformational changes of the enzyme lid once the substrate, phosphoserine (PSer), is bound in the active site. Here, Elastic Network Model (ENM) is applied to the crystal structure of hPSP to probe the transition between open and closed conformations of hPSP. Molecular Dynamics (MD) simulations are carried out on several PSer-hPSP systems to characterise the intermolecular interactions and their effect on the dynamics of the enzyme lid. Systems involving either Ca⁺⁺ or Mg⁺⁺ are considered. The first ENM normal mode shows that an open-closed transition can be explained from a simple description of the enzyme in terms of harmonic potentials. Principal Component Analyses applied to the MD trajectories also highlight a trend for a closing/opening motion. Different PSer orientations inside the enzyme cavity are identified, i.e. either the carboxylate, the phosphate group of PSer, or both, are oriented towards the cation. The interaction patterns are analysed in terms of hydrogen bonds, electrostatics, and bond critical points of the electron density distributions. The latter approach yields a global description of the bonding intermolecular interactions. The PSer orientation determines the content of the cation coordination shell and the mobility of the substrate, while Lys158 and Thr182, involved in the reaction mechanism, are always in interaction with the substrate. Closed enzyme conformations involve Met52-Gln204, Arg49-Glu29, and Arg50-Glu29 interactions. Met52, as well as Arg49 and Arg50, also stabilize PSer inside the cavity. Communicated by Ramaswamy H. Sarma.

Design of Calcium-Binding Proteins to Sense Calcium

Calcium controls numerous biological processes by interacting with different classes of calcium binding proteins (CaBP’s), with different affinities, metal selectivities, kinetics, and calcium dependent conformational changes. Due to the diverse coordination chemistry of calcium, and complexity associated with protein folding and binding cooperativity, the rational design of CaBP’s was anticipated to present multiple challenges. In this paper we will first discuss applications of statistical analysis of calcium binding sites in proteins and subsequent development of algorithms to predict and identify calcium binding proteins. Next, we report efforts to identify key determinants for calcium binding affinity, cooperativity and calcium dependent conformational changes using grafting and protein design. Finally, we report recent advances in designing protein calcium sensors to capture calcium dynamics in various cellular environments.

Structure functional insights into calcium binding during the activation of coagulation factor XIII A

The dimeric FXIII-A₂, a pro-transglutaminase is the catalytic part of the heterotetrameric coagulation FXIII-A₂B₂ complex that upon activation by calcium binding/thrombin cleavage covalently cross-links preformed fibrin clots protecting them from premature fibrinolysis. Our study characterizes the recently disclosed three calcium binding sites of FXIII-A concerning evolution, mutual crosstalk, thermodynamic activation profile, substrate binding, and interaction with other similarly charged ions. We demonstrate unique structural aspects within FXIII-A calcium binding sites that give rise to functional differences making FXIII unique from other transglutaminases. The first calcium binding site showed an antagonistic relationship towards the other two. The thermodynamic profile of calcium/thrombin-induced FXIII-A activation explains the role of bulk solvent in transitioning its zymogenic dimeric form to an activated monomeric form. We also explain the indirect effect of solvent ion concentration on FXIII-A activation. Our study suggests FXIII-A calcium binding sites could be putative pharmacologically targetable regions.

Designing Calcium-Binding Proteins for Molecular MR Imaging

Early diagnosis, noninvasive detection, and staging of various diseases, remain one of the major clinical barriers to effective medical treatment and prevention of disease progression toward major clinical consequences. Molecular imaging technologies play an indispensable role in the clinical field in overcoming these major barriers. The increasing application of imaging techniques and agents in early detection of different diseases such as cancer has resulted in improved treatment response and clinical patient management. In this chapter we will first introduce criteria for the design and engineering of calcium-binding protein (CaBP) parvalbumin as a protein Gd-MRI contrast agent (ProCA) with unprecedented metal selectivity for Gd³⁺ over physiological metal ions. We will then discuss the further development of targeted MRI contrast agent for molecular imaging of PSMA biomarker for early detection of prostate cancer.

The Different Facets of Extracellular Calcium Sensors: Old and New Concepts in Calcium-Sensing Receptor Signalling and Pharmacology

The current interest of the scientific community for research in the field of calcium sensing in general and on the calcium-sensing Receptor (CaR) in particular is demonstrated by the still increasing number of papers published on this topic. The extracellular calcium-sensing receptor is the best-known G-protein-coupled receptor (GPCR) able to sense external Ca2+ changes. Widely recognized as a fundamental player in systemic Ca2+ homeostasis, the CaR is ubiquitously expressed in the human body where it activates multiple signalling pathways. In this review, old and new notions regarding the mechanisms by which extracellular Ca2+ microdomains are created and the tools available to measure them are analyzed. After a survey of the main signalling pathways triggered by the CaR, a special attention is reserved for the emerging concepts regarding CaR function in the heart, CaR trafficking and pharmacology. Finally, an overview on other Ca2+ sensors is provided.

Molecular determinants of Ca2+ sensitivity at the intersubunit interface of the BK channel gating ring

The large-conductance calcium-activated K⁺ (BK) channel contains two intracellular tandem Ca²⁺-sensing RCK domains (RCK1 and RCK2), which tetramerize into a Ca²⁺ gating ring that regulates channel opening by conformational expansion in response to Ca²⁺ binding. Interestingly, the gating ring’s intersubunit assembly interface harbors the RCK2 Ca²⁺-binding site, known as the Ca²⁺ bowl. The gating ring’s assembly interface is made in part by intersubunit coordination of a Ca²⁺ ion between the Ca²⁺ bowl and an RCK1 Asn residue, N449, and by apparent intersubunit electrostatic interactions between E955 in RCK2 and R786 and R790 in the RCK2 of the adjacent subunit. To understand the role of the intersubunit assembly interface in Ca²⁺ gating, we performed mutational analyses of these putative interacting residues in human BK channels. We found that N449, despite its role in Ca²⁺ coordination, does not set the channel’s Ca²⁺ sensitivity, whereas E955 is a determinant of Ca²⁺ sensitivity, likely through intersubunit electrostatic interactions. Our findings provide evidence that the intersubunit assembly interface contains molecular determinants of Ca²⁺-sensitivity in BK channels.

Defining potential roles of Pb(2+) in neurotoxicity from a calciomics approach

Metal ions play crucial roles in numerous biological processes, facilitating biochemical reactions by binding to various proteins. An increasing body of evidence suggests that neurotoxicity associated with exposure to nonessential metals (e.g., Pb(2+)) involves disruption of synaptic activity, and these observed effects are associated with the ability of Pb(2+) to interfere with Zn(2+) and Ca(2+)-dependent functions. However, the molecular mechanism behind Pb(2+) toxicity remains a topic of debate. In this review, we first discuss potential neuronal Ca(2+) binding protein (CaBP) targets for Pb(2+) such as calmodulin (CaM), synaptotagmin, neuronal calcium sensor-1 (NCS-1), N-methyl-d-aspartate receptor (NMDAR) and family C of G-protein coupled receptors (cGPCRs), and their involvement in Ca(2+)-signalling pathways. We then compare metal binding properties between Ca(2+) and Pb(2+) to understand the structural implications of Pb(2+) binding to CaBPs. Statistical and biophysical studies (e.g., NMR and fluorescence spectroscopy) of Pb(2+) binding are discussed to investigate the molecular mechanism behind Pb(2+) toxicity. These studies identify an opportunistic, allosteric binding of Pb(2+) to CaM, which is distinct from ionic displacement. Together, these data suggest three potential modes of Pb(2+) activity related to molecular and/or neural toxicity: (i) Pb(2+) can occupy Ca(2+)-binding sites, inhibiting the activity of the protein by structural modulation, (ii) Pb(2+) can mimic Ca(2+) in the binding sites, falsely activating the protein and perturbing downstream activities, or (iii) Pb(2+) can bind outside of the Ca(2+)-binding sites, resulting in the allosteric modulation of the protein activity. Moreover, the data further suggest that even low concentrations of Pb(2+) can interfere at multiple points within the neuronal Ca(2+) signalling pathways to cause neurotoxicity.

Modulation of the FMRFamide-gated Na+ channel by external Ca2

FMRFamide-gated Na⁺ channel (FaNaC) is a member of the DEG/ENaC family. Amino acid sequence of the second transmembrane region (TM2) of FaNaC is quite similar to that of the acid-sensing ion channels (ASIC) of the same family. In the upper part of TM2, there are two aspartate residues (D552 and D556 in Aplysia FaNaC, AkFaNaC) which construct two negative rings in the external vestibule. In the present study, we examined the function of D552/D556 mutants of AkFaNaC in Xenopus oocytes with special interest in Ca²⁺ sensitivity of FaNaC. The FMRFamide-evoked current through AkFaNaC was depressed by submillimolar Ca²⁺ such that the current in Ca²⁺-free condition was 2-3-fold larger than that in the control solution which contained 1.8 mM CaCl ₂. Both D552 and D556 were found to be indispensable for the sensitivity of FaNaC to submillimolar Ca²⁺. Unexpectedly, however, both acidic residues were not essential for the inhibition by millimolar Ca²⁺ concentrations. The Ca²⁺-sensitive gating of FaNaC was recapitulated by an allosteric model in which Ca²⁺-bound channels are more difficult to open. The desensitization of FaNaC was also inhibited by Ca²⁺, which was abolished in some D552/D556 mutants. Structural models of FaNaC made by homology modeling showed that the distance between oxygen atoms of D552 and D556 on the adjacent subunits is close enough to coordinate Ca²⁺ in the nonconducting desensitized channel but not in the open channel. The results suggest that Ca²⁺ coordination between oxygen atoms of D552 and D556 disturbs the opening transition as well as the desensitization of FaNaC.

Molecular Basis for Modulation of Metabotropic Glutamate Receptors and Their Drug Actions by Extracellular Ca2

Metabotropic glutamate receptors (mGluRs) associated with the slow phase of the glutamatergic signaling pathway in neurons of the central nervous system have gained importance as drug targets for chronic neurodegenerative diseases. While extracellular Ca²⁺ was reported to exhibit direct activation and modulation via an allosteric site, the identification of those binding sites was challenged by weak binding. Herein, we review the discovery of extracellular Ca²⁺ in regulation of mGluRs, summarize the recent developments in probing Ca²⁺ binding and its co-regulation of the receptor based on structural and biochemical analysis, and discuss the molecular basis for Ca²⁺ to regulate various classes of drug action as well as its importance as an allosteric modulator in mGluRs.

Molecular Basis of the Extracellular Ligands Mediated Signaling by the Calcium Sensing Receptor

Ca²⁺-sensing receptors (CaSRs) play a central role in regulating extracellular calcium concentration ([Ca²⁺]_o) homeostasis and many (patho)physiological processes in multiple organs. This regulation is orchestrated by a cooperative response to extracellular stimuli such as small changes in Ca²⁺, Mg²⁺, amino acids, and other ligands. In addition, CaSR is a pleiotropic receptor regulating several intracellular signaling pathways, including calcium mobilization and intracellular calcium oscillation. Nearly 200 mutations and polymorphisms have been found in CaSR in relation to a variety of human disorders associated with abnormal Ca²⁺ homeostasis. In this review, we summarize efforts directed at identifying binding sites for calcium and amino acids. Both homotropic cooperativity among multiple calcium binding sites and heterotropic cooperativity between calcium and amino acid were revealed using computational modeling, predictions, and site-directed mutagenesis coupled with functional assays. The hinge region of the bilobed Venus flytrap (VFT) domain of CaSR plays a pivotal role in coordinating multiple extracellular stimuli, leading to cooperative responses from the receptor. We further highlight the extensive number of disease-associated mutations that have also been shown to affect CaSR's cooperative action via several types of mechanisms. These results provide insights into the molecular bases of the structure and functional cooperativity of this receptor and other members of family C of the G protein-coupled receptors (cGPCRs) in health and disease states, and may assist in the prospective development of novel receptor-based therapeutics.

Surface-Directed Assembly of Sequence-Defined Synthetic Polymers into Networks of Hexagonally Patterned Nanoribbons with Controlled Functionalities

The exquisite self-assembly of proteins and peptides in nature into highly ordered functional materials has inspired innovative approaches to the design and synthesis of biomimetic materials. While sequence-defined polymers hold great promise to mimic proteins and peptides for functions, controlled assembly of them on surfaces still remains underdeveloped. Here, we report the assembly of 12-mer peptoids containing alternating acidic and aromatic monomers into networks of hexagonally patterned nanoribbons on mica surfaces. Ca(2+)-carboxylate coordination creates peptoid-peptoid and peptoid-mica interactions that control self-assembly. In situ atomic force microscopy (AFM) shows that peptoids first assemble into discrete nanoparticles; these particles then transform into hexagonally patterned nanoribbons on mica surfaces. AFM-based dynamic force spectroscopy studies show that peptoid-mica interactions are much stronger than peptoid-peptoid interactions, illuminating the driving forces for mica-directed peptoid assembly. We further demonstrate the display of functional domains at the N-terminus of assembling peptoids to produce extended networks with similar hierarchical structures. This research demonstrates that surface-directed peptoid assembly can be used as a robust platform to develop biomimetic coating materials for applications.

Vasohibins: new transglutaminase-like cysteine proteases possessing a non-canonical Cys-His-Ser catalytic triad

Vasohibin-1 and Vasohibin-2 regulate angiogenesis, tumour growth and metastasis. Their molecular functions, however, were previously unknown, in large part owing to their perceived lack of homology to proteins of known structure and function. To identify their functional amino acids and domains, their molecular activity and their evolutionary history, we undertook an in-depth analysis of Vasohibin sequences. We find that Vasohibin proteins are previously undetected members of the transglutaminase-like cysteine protease superfamily, and all possess a non-canonical Cys-His-Ser catalytic triad. We further propose a calcium-dependent activation mechanism for Vasohibin proteins. These findings can now be used to design constructs for protein structure determination and to develop enzyme inhibitors as angiogenic regulators to treat metastasis and tumour growth.

CONTACT

luis.sanchezpulido@dpag.ox.ac.uk

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Fast kinetics of calcium signaling and sensor design

Fast calcium signaling is regulated by numerous calcium channels exhibiting high spatiotemporal profiles which are currently measured by fluorescent calcium sensors. There is still a strong need to improve the kinetics of genetically encoded calcium indicators (sensors) to capture calcium dynamics in the millisecond time frame. In this review, we summarize several major fast calcium signaling pathways and discuss the recent developments and application of genetically encoded calcium indicators to detect these pathways. A new class of genetically encoded calcium indicators designed with site-directed mutagenesis on the surface of beta-barrel fluorescent proteins to form a pentagonal bipyramidal-like calcium binding domain dramatically accelerates calcium binding kinetics. Furthermore, novel genetically encoded calcium indicators with significantly increased fluorescent lifetime change are advantageous in deep-field imaging with high light-scattering and notable morphology change.

Dimerization of peptides by calcium ions: investigation of a calcium-binding motif

We investigated calcium-binding motifs of peptides and their recognition of active functionalities for coordination. This investigation generates the fundamentals to design carrier material for calcium-bound peptide-peptide interactions. Interactions of different peptides with active calcium domains were investigated. Evaluation of selectivity was performed by electrospray ionization mass spectrometry by infusing solutions containing two different peptides (P₁ and P₂) in the presence of calcium ions. In addition to signals for monomer species, intense dimer signals are observed for the heterodimer ions (P₁ ⋯ Ca²⁺ ⋯ P₂) (⋯ represents the noncovalent binding of calcium with the peptide) in the positive ion mode and for ions ([P₁-2H]²⁻ ⋯ Ca²⁺ ⋯ [P₂-2H]²⁻) in the negative ion mode. Monitoring of the dissociation from these mass selected dimer ions via the kinetic method provides information on the calcium affinity order of different peptide sequences.

Prediction and analysis of canonical EF hand loop and qualitative estimation of Ca²⁺ binding affinity

The diversity of functions carried out by EF hand-containing calcium-binding proteins is due to various interactions made by these proteins as well as the range of affinity levels for Ca²⁺ displayed by them. However, accurate methods are not available for prediction of binding affinities. Here, amino acid patterns of canonical EF hand sequences obtained from available crystal structures were used to develop a classifier that distinguishes Ca²⁺-binding loops and non Ca²⁺-binding regions with 100% accuracy. To investigate further, we performed a proteome-wide prediction for E. histolytica, and classified known EF-hand proteins. We compared our results with published methods on the E. histolytica proteome scan, and demonstrated our method to be more specific and accurate for predicting potential canonical Ca²⁺-binding loops. Furthermore, we annotated canonical EF-hand motifs and classified them based on their Ca²⁺-binding affinities using support vector machines. Using a novel method generated from position-specific scoring metrics and then tested against three different experimentally derived EF-hand-motif datasets, predictions of Ca²⁺-binding affinities were between 87 and 90% accurate. Our results show that the tool described here is capable of predicting Ca²⁺-binding affinity constants of EF-hand proteins. The web server is freely available at http://202.41.10.46/calb/index.html.

Design of ProCAs (protein-based Gd(3+) MRI contrast agents) with high dose efficiency and capability for molecular imaging of cancer biomarkers

Magnetic resonance imaging (MRI) is the leading imaging technique for disease diagnostics, providing high resolution, three-dimensional images noninvasively. MRI contrast agents are designed to improve the contrast and sensitivity of MRI. However, current clinically used MRI contrast agents have relaxivities far below the theoretical upper limit, which largely prevent advancing molecular imaging of biomarkers with desired sensitivity and specificity. This review describes current progress in the development of a new class of protein-based MRI contrast agents (ProCAs) with high relaxivity using protein design to optimize the parameters that govern relaxivity. Further, engineering with targeting moiety allows these contrast agents to be applicable for molecular imaging of prostate cancer biomarkers by MRI. The developed protein-based contrast agents also exhibit additional in vitro and in vivo advantages for molecular imaging of disease biomarkers, such as high metal-binding stability and selectivity, reduced toxicity, proper blood circulation time, and higher permeability in tumor tissue in addition to improved relaxivities.

Identification of an L-phenylalanine binding site enhancing the cooperative responses of the calcium-sensing receptor to calcium

Functional positive cooperative activation of the extracellular calcium ([Ca(2+)]o)-sensing receptor (CaSR), a member of the family C G protein-coupled receptors, by [Ca(2+)]o or amino acids elicits intracellular Ca(2+) ([Ca(2+)]i) oscillations. Here, we report the central role of predicted Ca(2+)-binding site 1 within the hinge region of the extracellular domain (ECD) of CaSR and its interaction with other Ca(2+)-binding sites within the ECD in tuning functional positive homotropic cooperativity caused by changes in [Ca(2+)]o. Next, we identify an adjacent L-Phe-binding pocket that is responsible for positive heterotropic cooperativity between [Ca(2+)]o and L-Phe in eliciting CaSR-mediated [Ca(2+)]i oscillations. The heterocommunication between Ca(2+) and an amino acid globally enhances functional positive homotropic cooperative activation of CaSR in response to [Ca(2+)]o signaling by positively impacting multiple [Ca(2+)]o-binding sites within the ECD. Elucidation of the underlying mechanism provides important insights into the longstanding question of how the receptor transduces signals initiated by [Ca(2+)]o and amino acids into intracellular signaling events.

Molecular imaging of EGFR/HER2 cancer biomarkers by protein MRI contrast agents

Epidermal growth factor receptor (EGFR) and HER2 are major prognosis biomarkers and drug targets overexpressed in various types of cancer cells. There is a pressing need to develop MRI contrast agents capable of enhancing the contrast between normal tissues and tumors with high relaxivity, capable of targeting tumors, and with high intratumoral distribution and minimal toxicity. In this review, we first discuss EGFR signaling and its role in tumor progression as a major drug target. We then report our progress in the development of protein contrast agents with significant improvement of both r1 and r2 relaxivities, pharmacokinetics, in vivo retention time, and in vivo dose efficiency. Finally, we report our effort in the development of EGFR-targeted protein contrast agents with the capability to cross the endothelial boundary and with good tissue distribution across the entire tumor mass. The noninvasive capability of MRI to visualize spatially and temporally the intratumoral distribution as well as quantify the levels of EGFR and HER2 would greatly improve our ability to track changes of the biomarkers during tumor progression, monitor treatment efficacy, aid in patient selection, and further develop novel targeted therapies for clinical application.

Extracellular calcium modulates actions of orthosteric and allosteric ligands on metabotropic glutamate receptor 1α

Metabotropic glutamate receptor 1α (mGluR1α), a member of the family C G protein-coupled receptors, is emerging as a potential drug target for various disorders, including chronic neuronal degenerative diseases. In addition to being activated by glutamate, mGluR1α is also modulated by extracellular Ca(2+). However, the underlying mechanism is unknown. Moreover, it has long been challenging to develop receptor-specific agonists due to homologies within the mGluR family, and the Ca(2+)-binding site(s) on mGluR1α may provide an opportunity for receptor-selective targeting by therapeutics. In the present study, we show that our previously predicted Ca(2+)-binding site in the hinge region of mGluR1α is adjacent to the site where orthosteric agonists and antagonists bind on the extracellular domain of the receptor. Moreover, we found that extracellular Ca(2+) enhanced mGluR1α-mediated intracellular Ca(2+) responses evoked by the orthosteric agonist l-quisqualate. Conversely, extracellular Ca(2+) diminished the inhibitory effect of the mGluR1α orthosteric antagonist (S)-α-methyl-4-carboxyphenylglycine. In addition, selective positive (Ro 67-4853) and negative (7-(hydroxyimino)cyclopropa[b]chromen-1a-carboxylate ethyl ester) allosteric modulators of mGluR1α potentiated and inhibited responses to extracellular Ca(2+), respectively, in a manner similar to their effects on the response of mGluR1α to glutamate. Mutations at residues predicted to be involved in Ca(2+) binding, including E325I, had significant effects on the modulation of responses to the orthosteric agonist l-quisqualate and the allosteric modulator Ro 67-4853 by extracellular Ca(2+). These studies reveal that binding of extracellular Ca(2+) to the predicted Ca(2+)-binding site in the extracellular domain of mGluR1α modulates not only glutamate-evoked signaling but also the actions of both orthosteric ligands and allosteric modulators on mGluR1α.

Calciomics: integrative studies of Ca2+-binding proteins and their interactomes in biological systems

Calcium ion (Ca(2+)), the fifth most common chemical element in the earth's crust, represents the most abundant mineral in the human body. By binding to a myriad of proteins distributed in different cellular organelles, Ca(2+) impacts nearly every aspect of cellular life. In prokaryotes, Ca(2+) plays an important role in bacterial movement, chemotaxis, survival reactions and sporulation. In eukaryotes, Ca(2+) has been chosen through evolution to function as a universal and versatile intracellular signal. Viruses, as obligate intracellular parasites, also develop smart strategies to manipulate the host Ca(2+) signaling machinery to benefit their own life cycles. This review focuses on recent advances in applying both bioinformatic and experimental approaches to predict and validate Ca(2+)-binding proteins and their interactomes in biological systems on a genome-wide scale (termed "calciomics"). Calmodulin is used as an example of Ca(2+)-binding protein (CaBP) to demonstrate the role of CaBPs on the regulation of biological functions. This review is anticipated to rekindle interest in investigating Ca(2+)-binding proteins and Ca(2+)-modulated functions at the systems level in the post-genomic era.

Tailoring encodable lanthanide-binding tags as MRI contrast agents

Lanthanide-binding tags (LBTs), peptide-based coexpression tags with high affinity for lanthanide ions, have previously been applied as luminescent probes to provide phasing for structure determination in X-ray crystallography and to provide restraints for structural refinement and distance information in NMR. The native affinity of LBTs for Gd(3+) indicates their potential as the basis for engineering of peptide-based MRI agents. However, the lanthanide coordination state that enhances luminescence and affords tightest binding would not be ideal for applications of LBTs as contrast agents, due to the exclusion of water from the inner coordination sphere. Herein, we use structurally defined LBTs as the starting point for re-engineering the first coordination shell of the lanthanide ion to provide for high contrast through direct coordination of water to Gd(3+) (resulting in the single LBT peptide, m-sLBT). The effectiveness of LBTs as MRI contrast agents was examined in vitro through measurement of binding affinity and proton relaxivity. For imaging applications that require targeted observation, fusion to specific protein partners is desirable. However, a fusion protein comprising a concatenated double LBT (dLBT) as an N-terminal tag for the model protein ubiquitin had reduced relaxivity compared with the free dLBT peptide. This limitation was overcome by the use of a construct based on the m-sLBT sequence (q-dLBT-ubiquitin). The structural basis for the enhanced contrast was examined by comparison of the X-ray crystal structure of xq-dLBT-ubiquitin (wherein two tryptophan residues are replaced with serine), to that of dLBT-ubiquitin. The structure shows that the backbone conformational dynamics of the MRI variant may allow enhanced water exchange. This engineered LBT represents a first step in expanding the current base of specificity-targeted agents available.

Predicting Ca2+ -binding sites using refined carbon clusters

Identifying Ca(2+) -binding sites in proteins is the first step toward understanding the molecular basis of diseases related to Ca(2+) -binding proteins. Currently, these sites are identified in structures either through X-ray crystallography or NMR analysis. However, Ca(2+) -binding sites are not always visible in X-ray structures due to flexibility in the binding region or low occupancy in a Ca(2+) -binding site. Similarly, both Ca(2+) and its ligand oxygens are not directly observed in NMR structures. To improve our ability to predict Ca(2+) -binding sites in both X-ray and NMR structures, we report a new graph theory algorithm (MUG(C) ) to predict Ca(2+) -binding sites. Using carbon atoms covalently bonded to the chelating oxygen atoms, and without explicit reference to side-chain oxygen ligand co-ordinates, MUG(C) is able to achieve 94% sensitivity with 76% selectivity on a dataset of X-ray structures composed of 43 Ca(2+) -binding proteins. Additionally, prediction of Ca(2+) -binding sites in NMR structures was obtained by MUG(C) using a different set of parameters, which were determined by the analysis of both Ca(2+) -constrained and unconstrained Ca(2+) -loaded structures derived from NMR data. MUG(C) identified 20 of 21 Ca(2+) -binding sites in NMR structures inferred without the use of Ca(2+) constraints. MUG(C) predictions are also highly selective for Ca(2+) -binding sites as analyses of binding sites for Mg(2+) , Zn(2+) , and Pb(2+) were not identified as Ca(2+) -binding sites. These results indicate that the geometric arrangement of the second-shell carbon cluster is sufficient not only for accurate identification of Ca(2+) -binding sites in NMR and X-ray structures but also for selective differentiation between Ca(2+) and other relevant divalent cations.

Structural characterization of a conserved, calcium-dependent periplasmic protease from Legionella pneumophila

The bacterial dinucleotide second messenger c-di-GMP has emerged as a central molecule in regulating bacterial behavior, including motility and biofilm formation. Proteins for the synthesis and degradation of c-di-GMP and effectors for its signal transmission are widely used in the bacterial domain. Previous work established the GGDEF-EAL domain-containing receptor LapD as a central switch in Pseudomonas fluorescens cell adhesion. LapD senses c-di-GMP inside the cytosol and relays this signal to the outside by the differential recruitment of the periplasmic protease LapG. Here we identify the core components of an orthologous system in Legionella pneumophila. Despite only moderate sequence conservation at the protein level, key features concerning the regulation of LapG are retained. The output domain of the LapD-like receptor from L. pneumophila, CdgS9, binds the LapG ortholog involving a strictly conserved surface tryptophan residue. While the endogenous substrate for L. pneumophila LapG is unknown, the enzyme processed the corresponding P. fluorescens substrate, indicating a common catalytic mechanism and substrate recognition. Crystal structures of L. pneumophila LapG provide the first atomic models of bacterial proteases of the DUF920 family and reveal a conserved calcium-binding site important for LapG function.

Predicting the disruption by UO2(2+) of a protein-ligand interaction

The uranyl cation (UO(2) (2+)) can be suspected to interfere with the binding of essential metal cations to proteins, underlying some mechanisms of toxicity. A dedicated computational screen was used to identify UO(2) (2+) binding sites within a set of nonredundant protein structures. The list of potential targets was compared to data from a small molecules interaction database to pinpoint specific examples where UO(2) (2+) should be able to bind in the vicinity of an essential cation, and would be likely to affect the function of the corresponding protein. The C-reactive protein appeared as an interesting hit since its structure involves critical calcium ions in the binding of phosphorylcholine. Biochemical experiments confirmed the predicted binding site for UO(2) (2+) and it was demonstrated by surface plasmon resonance assays that UO(2) (2+) binding to CRP prevents the calcium-mediated binding of phosphorylcholine. Strikingly, the apparent affinity of UO(2) (2+) for native CRP was almost 100-fold higher than that of Ca(2+). This result exemplifies in the case of CRP the capability of our computational tool to predict effective binding sites for UO(2) (2+) in proteins and is a first evidence of calcium substitution by the uranyl cation in a native protein.

A protein sequence meta-functional signature for calcium binding residue prediction

The diversity of characterized protein functions found amongst experimentally interrogated proteins suggests that a vast array of unknown functions remains undiscovered. These protein functions are imparted by specific geometric distributions of amino acid residue chemical moieties, each contributing a functional interaction. We hypothesize that individual residue function contributions are predictable through sequence analytic knowledge based algorithms, and that they can be recombined to understand composite protein function by predicting spatial relation in tertiary structure. We assess the former by training a meta-functional signature algorithm to specifically predict calcium ion binding residues from protein sequence. We estimate the latter by testing for match between predictive contribution of positions in predicted secondary structures and patterns of side chain proximity forced by secondary structure moieties. Specific training for calcium binding results in 83% area under the receiver operator characteristic curve added value over random (AUCoR) and p<10(-300) significance as measured by Kendall's τ in ten fold cross validation for parallel sets of 811 residues in 336 proteins and 696 residues in 299 proteins. Training for generalized function results in 63% AUCoR and p≅10(-221) for the same tests. Including inference of side chain proximity improves predictive ability by 2% AUCoR consistently. The results demonstrate that protein meta-functional signatures can be trained to predict specific protein functions by considering amino acid identity and structural features accessible from sequence, laying the groundwork for composite sequence based function site prediction.

Analysis and prediction of calcium-binding pockets from apo-protein structures exhibiting calcium-induced localized conformational changes

Calcium binding in proteins exhibits a wide range of polygonal geometries that relate directly to an equally diverse set of biological functions. The binding process stabilizes protein structures and typically results in local conformational change and/or global restructuring of the backbone. Previously, we established the MUG program, which utilized multiple geometries in the Ca(2+)-binding pockets of holoproteins to identify such pockets, ignoring possible Ca(2+)-induced conformational change. In this article, we first report our progress in the analysis of Ca(2+)-induced conformational changes followed by improved prediction of Ca(2+)-binding sites in the large group of Ca(2+)-binding proteins that exhibit only localized conformational changes. The MUG(SR) algorithm was devised to incorporate side chain torsional rotation as a predictor. The output from MUG(SR) presents groups of residues where each group, typically containing two to five residues, is a potential binding pocket. MUG(SR) was applied to both X-ray apo structures and NMR holo structures, which did not use calcium distance constraints in structure calculations. Predicted pockets were validated by comparison with homologous holo structures. Defining a "correct hit" as a group of residues containing at least two true ligand residues, the sensitivity was at least 90%; whereas for a "correct hit" defined as a group of residues containing at least three true ligand residues, the sensitivity was at least 78%. These data suggest that Ca(2+)-binding pockets are at least partially prepositioned to chelate the ion in the apo form of the protein.

Elucidation of a novel extracellular calcium-binding site on metabotropic glutamate receptor 1{alpha} (mGluR1{alpha}) that controls receptor activation

Metabotropic glutamate receptor 1α (mGluR1α) exerts important effects on numerous neurological processes. Although mGluR1α is known to respond to extracellular Ca(2+) ([Ca(2+)](o)) and the crystal structures of the extracellular domains (ECDs) of several mGluRs have been determined, the calcium-binding site(s) and structural determinants of Ca(2+)-modulated signaling in the Glu receptor family remain elusive. Here, we identify a novel Ca(2+)-binding site in the mGluR1α ECD using a recently developed computational algorithm. This predicted site (comprising Asp-318, Glu-325, and Asp-322 and the carboxylate side chain of the receptor agonist, Glu) is situated in the hinge region in the ECD of mGluR1α adjacent to the reported Glu-binding site, with Asp-318 involved in both Glu and calcium binding. Mutagenesis studies indicated that binding of Glu and Ca(2+) to their distinct but partially overlapping binding sites synergistically modulated mGluR1α activation of intracellular Ca(2+) ([Ca(2+)](i)) signaling. Mutating the Glu-binding site completely abolished Glu signaling while leaving its Ca(2+)-sensing capability largely intact. Mutating the predicted Ca(2+)-binding residues abolished or significantly reduced the sensitivity of mGluR1α not only to [Ca(2+)](o) and [Gd(3+)](o) but also, in some cases, to Glu. The dual activation of mGluR1α by [Ca(2+)](o) and Glu has important implications for the activation of other mGluR subtypes and related receptors. It also opens up new avenues for developing allosteric modulators of mGluR function that target specific human diseases.

Integration of Diverse Research Methods to Analyze and Engineer Ca-Binding Proteins: From Prediction to Production

In recent years, increasingly sophisticated computational and bioinformatics tools have evolved for the analyses of protein structure, function, ligand interactions, modeling and energetics. This includes the development of algorithms to recursively evaluate side-chain rotamer permutations, identify regions in a 3D structure that meet some set of search parameters, calculate and minimize energy values, and provide high-resolution visual tools for theoretical modeling. Here we discuss the interdependency between different areas of bioinformatics, the evolution of different algorithm design approaches, and finally the transition from theoretical models to real-world design and application as they relate to Ca(2+)-binding proteins. Within this context, it has become evident that significant pre-experimental design and calculations can be modeled through computational methods, thus eliminating potentially unproductive research and increasing our confidence in the correlation between real and theoretical models. Moving from prediction to production, it is anticipated that bioinformatics tools will play an increasingly significant role in research and development, improving our ability to both understand the physiological roles of Ca(2+) and other metals and to extend that knowledge to the design of function-specific synthetic proteins capable of fulfilling different roles in medical diagnostics and therapeutics.