Negative Charge Cluster Research Articles

Charge cluster occurrence in land plants’ mitochondrial proteomes with functional and structural insights

The Charge Clusters (CCs) are involved in key functions and are distributed according to the organism, the protein’s type, and the charge of amino acids. In the present study, we have explored the occurrence, position, and annotation as a first large-scale study of the CCs in land plants mitochondrial proteomes. A new python script was used for data curation. The Finding Clusters Charge in Protein Sequences Program was performed after adjusting the reading window size. A 44316 protein sequences belonging to 52 species of land plants were analysed. The occurrence of Negative Charge Clusters (NCCs) (1.2%) is two times more frequent than the Positive Charge Clusters (PCCs) (0.64%). Moreover, 39 and 30 NCCs were conserved in 88 and 41 proteins in intra and in inter proteomes respectively, while 14 and 21 PCCs were conserved in 53 and 85 protein sequences in intra and inter proteomes consecutively. Sequences carrying mixed CCs are rare (0.12%). Despite this low abundance, CCs play a crucial role in protein function. The CCs tend to be located mainly in the terminal regions of proteins which guarantees specific protein targeting and import into the mitochondria. In addition, the functional annotation of CCs according to Gene Ontology shows that CCs are involved in binding functions of either proteins or macromolecules which are deployed in different metabolic and cellular processes such as RNA editing and transcription. This study may provide valuable information while considering the CCs in understanding the environmental adaptation of plants. Communicated by Ramaswamy H. Sarma

Comprehensive Insights into the Catalytic Mechanism of Middle East Respiratory Syndrome 3C-Like Protease and Severe Acute Respiratory Syndrome 3C-Like Protease.

Coronavirus 3C-like protease (3CLPro) is a highly conserved cysteine protease employing a catalytic dyad for its functions. 3CLPro is essential to the viral life cycle and, therefore, is an attractive target for developing antiviral agents. However, the detailed catalytic mechanism of coronavirus 3CLPro remains largely unknown. We took an integrated approach of employing X-ray crystallography, mutational studies, enzyme kinetics study, and inhibitors to gain insights into the mechanism. Such experimental work is supplemented by computational studies, including the prereaction state analysis, the ab initio calculation of the critical catalytic step, and the molecular dynamic simulation of the wild-type and mutant enzymes. Taken together, such studies allowed us to identify a residue pair (Glu-His) and a conserved His as critical for binding; a conserved GSCGS motif as important for the start of catalysis, a partial negative charge cluster (PNCC) formed by Arg-Tyr-Asp as essential for catalysis, and a conserved water molecule mediating the remote interaction between PNCC and catalytic dyad. The data collected and our insights into the detailed mechanism have allowed us to achieve a good understanding of the difference in catalytic efficiency between 3CLPro from SARS and MERS, conduct mutational studies to improve the catalytic activity by 8-fold, optimize existing inhibitors to improve the potency by 4-fold, and identify a potential allosteric site for inhibitor design. All such results reinforce each other to support the overall catalytic mechanism proposed herein.

Charge clusters signatures in prokaryotic proteomes: Temperature-dependence and distribution

The identification and the screening of Charged Clusters (CCs) residues in proteins is a key analysis to assess any quantitative structure-function correlation in proteins. Here, we present a proteome wide scan for the occurrence of (CCs) in 99292 proteins using a new tool. Finding Clusters Charge in Protein Sequences Program (FCCP). The FCCP has been employed to search CCs in 35 prokaryotic proteomes (7 Psychrophiles, 10 Mesophiles, 9 thermophiles and for 9 hyperthermophiles). A new repository of 855 CC has been created. Results showed that the mapped proteins containing positive and negative charge clusters are mostly transmembrane proteins while the conserved CCs within the same proteome are transposases or involved in DNA binding and integration. Interestingly, the negative charged cluster was associated to bacteria growth's temperature (p=0.002) acting as proteins' core signature. Taken together the various results provide a consistent picture of these screened CCs in terms of its potentials functional roles.

Screening for clusters of charge in human virus proteomes.

BackgroundThe identification of charge clusters (runs of charged residues) in proteins and their mapping within the protein structure sequence is an important step toward a comprehensive analysis of how these particular motifs mediate, via electrostatic interactions, various molecular processes such as protein sorting, translocation, docking, orientation and binding to DNA and to other proteins. Few algorithms that specifically identify these charge clusters have been designed and described in the literature. In this study, 197 distinctive human viral proteomes were screened for the occurrence of charge clusters (CC) using a new computational approach.ResultsThree hundred and seventy three CC have been identified within the 2549 viral protein sequences screened. The number of protein sequences that are CC-free is 2176 (85.3 %) while 150 and 180 proteins contained positive charge (PCC) and negative charge clusters (NCC), respectively. The NCCs (211 detected) were more prevalent than PCC (162). PCC-containing proteins are significantly longer than those having NCCs (p = 2.10-16). The most prevalent virus families having PCC and NCC were Herpesviridae followed by Papillomaviridae. However, the single-strand RNA group has in average three times more NCC than PCC. According to the functional domain classification, a significant difference in distribution was observed between PCC and NCC (p = 2. 10−8) with the occurrence of NCCs being more frequent in C-terminal region while PCC more often fall within functional domains. Only 29 proteins sequences contained both NCC and PCC. Moreover, 101 NCC were conserved in 84 proteins while only 62 PCC were conserved in 60 protein sequences. To understand the mechanism by which the membrane translocation functionalities are embedded in viral proteins, we screened our PCC for sequences corresponding to cell-penetrating peptides (CPPs) using two online databases: CellPPd and CPPpred. We found that all our PCCs, having length varying from 7 to 30 amino-acids were predicted as CPPs. Experimental validation is required to improve our understanding of the role of these PCCs in viral infection process.ConclusionsScreening distinctive cluster charges in viral proteomes suggested a functional role of these protein regions and might provide potential clues to improve the current understanding of viral diseases in order to tailor better preventive and therapeutic approaches.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3086-3) contains supplementary material, which is available to authorized users.

Negatively charged residues of the segment linking the enzyme and cytolysin moieties restrict the membrane-permeabilizing capacity of adenylate cyclase toxin

The whooping cough agent, Bordetella pertussis, secretes an adenylate cyclase toxin-hemolysin (CyaA) that plays a crucial role in host respiratory tract colonization. CyaA targets CR3-expressing cells and disrupts their bactericidal functions by delivering into their cytosol an adenylate cyclase enzyme that converts intracellular ATP to cAMP. In parallel, the hydrophobic domain of CyaA forms cation-selective pores that permeabilize cell membrane. The invasive AC and pore-forming domains of CyaA are linked by a segment that is unique in the RTX cytolysin family. We used mass spectrometry and circular dichroism to show that the linker segment forms α-helical structures that penetrate into lipid bilayer. Replacement of the positively charged arginine residues, proposed to be involved in target membrane destabilization by the linker segment, reduced the capacity of the toxin to translocate the AC domain across cell membrane. Substitutions of negatively charged residues then revealed that two clusters of negative charges within the linker segment control the size and the propensity of CyaA pore formation, thereby restricting the cell-permeabilizing capacity of CyaA. The ‘AC to Hly-linking segment’ thus appears to account for the smaller size and modest cell-permeabilizing capacity of CyaA pores, as compared to typical RTX hemolysins.

Exploring proteome-wide occurrence of clusters of charged residues in eukaryotes.

Clusters of charged residues are one of the key features of protein primary structure since they have been associated to important functions of proteins. Here, we present a proteome wide scan for the occurrence of Charge Clusters in Protein sequences using a new search tool (FCCP) based on a score-based methodology. The FCCP was run to search charge clusters in seven eukaryotic proteomes: Arabidopsis thaliana, Caenorhabditis elegans, Danio rerio, Drosophila melanogaster, Homo sapiens, Mus musculus, and Saccharomyces cerevisiae. We found that negative charge clusters (NCCs) are three to four times more frequent than positive charge clusters (PCCs). The Drosophila proteome is on average the most charged, whereas the human proteome is the least charged. Only 3 to 8% of the studied protein sequences have negative charge clusters, while 1.6 to 3% having PCCs and only 0.07 to 0.6% have both types of clusters. NCCs are localized predominantly in the N-terminal and C-terminal domains, while PCCs tend to be localized within the functional domains of the protein sequences. Furthermore, the gene ontology classification revealed that the protein sequences with negative and PCCs are mainly binding proteins.

Asymmetric functional contributions of acidic and aromatic side chains in sodium channel voltage-sensor domains.

Voltage-gated sodium (NaV) channels mediate electrical excitability in animals. Despite strong sequence conservation among the voltage-sensor domains (VSDs) of closely related voltage-gated potassium (KV) and NaV channels, the functional contributions of individual side chains in Nav VSDs remain largely enigmatic. To this end, natural and unnatural side chain substitutions were made in the S2 hydrophobic core (HC), the extracellular negative charge cluster (ENC), and the intracellular negative charge cluster (INC) of the four VSDs of the skeletal muscle sodium channel isoform (NaV1.4). The results show that the highly conserved aromatic side chain constituting the S2 HC makes distinct functional contributions in each of the four NaV domains. No obvious cation-pi interaction exists with nearby S4 charges in any domain, and natural and unnatural mutations at these aromatic sites produce functional phenotypes that are different from those observed previously in Kv VSDs. In contrast, and similar to results obtained with Kv channels, individually neutralizing acidic side chains with synthetic derivatives and with natural amino acid substitutions in the INC had little or no effect on the voltage dependence of activation in any of the four domains. Interestingly, countercharge was found to play an important functional role in the ENC of DI and DII, but not DIII and DIV. These results suggest that electrostatic interactions with S4 gating charges are unlikely in the INC and only relevant in the ENC of DI and DII. Collectively, our data highlight domain-specific functional contributions of highly conserved side chains in NaV VSDs.

S1–S3 counter charges in the voltage sensor module of a mammalian sodium channel regulate fast inactivation

The movement of positively charged S4 segments through the electric field drives the voltage-dependent gating of ion channels. Studies of prokaryotic sodium channels provide a mechanistic view of activation facilitated by electrostatic interactions of negatively charged residues in S1 and S2 segments, with positive counterparts in the S4 segment. In mammalian sodium channels, S4 segments promote domain-specific functions that include activation and several forms of inactivation. We tested the idea that S1–S3 countercharges regulate eukaryotic sodium channel functions, including fast inactivation. Using structural data provided by bacterial channels, we constructed homology models of the S1–S4 voltage sensor module (VSM) for each domain of the mammalian skeletal muscle sodium channel hNaV1.4. These show that side chains of putative countercharges in hNaV1.4 are oriented toward the positive charge complement of S4. We used mutagenesis to define the roles of conserved residues in the extracellular negative charge cluster (ENC), hydrophobic charge region (HCR), and intracellular negative charge cluster (INC). Activation was inhibited with charge-reversing VSM mutations in domains I–III. Charge reversal of ENC residues in domains III (E1051R, D1069K) and IV (E1373K, N1389K) destabilized fast inactivation by decreasing its probability, slowing entry, and accelerating recovery. Several INC mutations increased inactivation from closed states and slowed recovery. Our results extend the functional characterization of VSM countercharges to fast inactivation, and support the premise that these residues play a critical role in domain-specific gating transitions for a mammalian sodium channel.

Revisiting the contribution of negative charges on the chaperonin cage wall to the acceleration of protein folding

Chaperonin GroEL mediates the folding of protein encapsulated in a GroES-sealed cavity (cage). Recently, a critical role of negative charge clusters on the cage wall in folding acceleration was proposed based on experiments using GroEL single-ring (SR) mutants SR1 and SRKKK2 [Tang YC, et al. (2006) Cell 125:903-914; Chakraborty K, et al. (2010) Cell 142:112-122]. Here, we revisited these experiments and discovered several inconsistencies. (i) SR1 was assumed to bind to GroES stably and to mediate single-round folding in the cage. However, we show that SR1 repeats multiple turnovers of GroES release/binding coupled with ATP hydrolysis. (ii) Although the slow folding observed for a double-mutant of maltose binding protein (DMMBP) by SRKKK2 was attributed to mutations that neutralize negative charges on the cage wall, we found that the majority of DMMBP escape from SRKKK2 and undergo spontaneous folding in the bulk medium. (iii) An osmolyte, trimethylamine N-oxide, was reported to accelerate SRKKK2-mediated folding of DMMBP by mimicking the effect of cage-wall negative charges of WT GroEL and ordering the water structure to promote protein compaction. However, we demonstrate that in-cage folding by SRKKK2 is unaffected by trimethylamine N-oxide. (iv) Although it was reported that SRKKK2 lost the ability to assist the folding of ribulose-1,5-bisphosphate carboxylase/oxygenase, we found that SRKKK2 retains this ability. Our results argue against the role of the negative charges on the cage wall of GroEL in protein folding. Thus, in chaperonin studies, folding kinetics need to be determined from the fraction of the real in-cage folding.

Chaperonin-Catalyzed Rescue of Kinetically Trapped States in Protein Folding

GroEL and GroES form a chaperonin nano-cage for single protein molecules to fold in isolation. The folding properties that render a protein chaperonin dependent are not yet understood. Here, we address this question using a double mutant of the maltose-binding protein DM-MBP as a substrate. Upon spontaneous refolding, DM-MBP populates a kinetically trapped intermediate that is collapsed but structurally disordered. Introducing two long-range disulfide bonds into DM-MBP reduces the entropic folding barrier of this intermediate and strongly accelerates native state formation. Strikingly, steric confinement of the protein in the chaperonin cage mimics the kinetic effect of constraining disulfides on folding, in a manner mediated by negative charge clusters in the cage wall. These findings suggest that chaperonin dependence correlates with the tendency of proteins to populate entropically stabilized folding intermediates. The capacity to rescue proteins from such folding traps may explain the uniquely essential role of chaperonin cages within the cellular chaperone network.

Comparative analysis of the Epstein-Barr virus encoded nuclear proteins of EBNA-3 family

It is known that the EBNA-3 family proteins (EBNA-3, -4 and -6, alternative nomenclature EBNA-3A, B and C correspondingly) show a limited sequence similarity. We have analyzed EBNA-3 proteins both at the primary sequence and secondary structure levels. EBNA-3 and EBNA-4 were structurally more similar compared to other combinations with EBNA-6. We found “Stonin Homology Domain” profile in EBNA-4 and “Proline Rich Domain” in all EBNA-3 family of proteins. We have also found positive and negative charge clusters in all three proteins and mixed charge clusters in EBNA-3. Charged clusters are believed to play an important role in interactions with DNA or signaling proteins. Additionally, unique primary sequence repeats were found in all three proteins.

Mechanical Force and Cytoplasmic Ca2+ Activate Yeast TRPY1 in Parallel

The ability to sense mechanical and osmotic stimuli is vital to all organisms from mammals to bacteria. Members of the transient receptor potential (TRP) ion-channel family have attracted intense attention for their involvement in mechanosensation. The yeast homologue TRPY1 can clearly be activated by hypertonic shock in vivo and by stretch force under patch clamp. Like its animal counterparts, TRPY1 is polymodal, being gated by membrane stretch force and by cytoplasmic Ca(2+). Here, we investigated how these two gating principles interact. We found that stretch force can induce some channel activation without cytoplasmic Ca(2+). Tens of micromolar Ca(2+) greatly enhance the observed force-induced activities, with open probabilities following well the Boltzmann distribution, in which the two gating energies are summed as exponents. To map this formalism to structures, we found Ca(2+)-binding proteins such as calmodulin or calcineurin to be unnecessary. However, removing a dense cluster of negative charges in the C-terminal cytoplasmic domain of TRPY1 greatly diminishes the Ca(2+) activation as well as its influence on force activation. We also found a strategic point upstream of this charge cluster, at which insertion of amino acids weakens Ca(2+) activation considerably but leaves the mechanosensitivity nearly intact. These results led to a structure-function model in which Ca(2+) binding to the cytoplasmic domain and stretching of the membrane-embedded domain both generate gating force, reaching the gate in parallel.

Gangliosides Potentially Inhibit Extracellular Nucleotide Metabolism

Gangliosides are glycolipids that contain sialic acid and they are mainly located on the outer leaflet of the cellular plasma membrane of most vertebrate and some invertebrate cells. Because they have structurally diverse, bulky and negatively charged oligosaccharide moieties, gangliosides endow cell membranes with unique molecular characteristics. Although they are abundant in the central nervous system (CNS), the complete loss of gangliosides in mice does not result in gross morphological abnormalities of the CNS. However, mutant mice develop neurodegenerative diseases and die soon after birth, suggesting that gangliosides are required for the maintenance and development of a stable CNS and are crucial to sustain life. At the cellular level, gangliosides influence cell growth and death, probably because they are involved in the lipid-mediated assembly of signaling molecules such as growth factor receptors or integrins on the membranes. This article addresses the structural similarity between the tandem sialic acid residues of gangliosides and nicotinamide adenine dinucleotide (NAD(+)) determined from biochemical data showing that gangliosides inhibit NAD(+) glycohydrolase activity and theoretical considerations. An essential feature of the structural similarity resides in a negative charge cluster formed by the two carboxyl groups in the tandem sialic acid residues and the diphosphate moiety of NAD(+). The potential physiological role(s) of gangliosides on the regulation of extracellular nucleotide metabolism are discussed.

Structural basis for unique mechanisms of folding and hemoglobin binding by a malarial protease.

Falcipain-2 (FP2), the major cysteine protease of the human malaria parasite Plasmodium falciparum, is a hemoglobinase and promising drug target. Here we report the crystal structure of FP2 in complex with a protease inhibitor, cystatin. The FP2 structure reveals two previously undescribed cysteine protease structural motifs, designated FP2(nose) and FP2(arm), in addition to details of the active site that will help focus inhibitor design. Unlike most cysteine proteases, FP2 does not require a prodomain but only the short FP2(nose) motif to correctly fold and gain catalytic activity. Our structure and mutagenesis data suggest a molecular basis for this unique mechanism by highlighting the functional role of two Tyr within FP2(nose) and a conserved Glu outside this motif. The FP2(arm) motif is required for hemoglobinase activity. The structure reveals topographic features and a negative charge cluster surrounding FP2(arm) that suggest it may serve as an exo-site for hemoglobin binding. Motifs similar to FP2(nose) and FP2(arm) are found only in related plasmodial proteases, suggesting that they confer malaria-specific functions.

Structures and Electronic States of Gadolinium Oxide Clusters

We have carried out ab initio quantum chemical calculations for gadolinium oxide clusters and clarified the most suitable combination of a theoretical method, basis sets, and effective core potentials for large systems of gadolinium oxides. We confirmed the reliability of our calculation by computing the electronic state of GdO molecule whose bond length and binding energy are well known. In order to examine how the electronic states depend on the size of clusters, we calculated several cluster series: GdOn clusters (n=2–4) with negative net charge, (GdO)n colinear and zigzag chain clusters (n=1–4), (GdO)n square pillar clusters (n=4,6,8), and large clusters based on the Gd2O3 cubic crystal (Gd4O1824-·8Q3+ and Gd4O1926-·8Q3+), where Q denotes the point charge. New images of electronic stresses in the clusters can be visualized in terms of the electronic stress tensor density based on the regional density functional theory.

Neuroblastoma cell-adapted yellow fever virus: mutagenesis of the E protein locus involved in persistent infection and its effects on virus penetration and spread

Persistent infection of mouse neuroblastoma NB41A3 cells with yellow fever 17D virus generates viral variants which exhibit defective cell penetration, poor cell-to-cell spread, small plaque size and reduced growth efficiency, caused by substitution of glycine for aspartic acid or glutamic acid at positions 360 and 362 in the envelope protein. These positions occur within a charge cluster, Asp360-Asp361-Glu362, located in domain III, near its interface with domain I. To characterize further the molecular basis for the variant phenotype, a series of mutant viruses containing substitutions at position 360, 361 and 362, were studied for effects on the cell culture properties typical of the neuroblastoma-adapted variant. Most substitutions at position 360 gave rise to viruses that were very defective in cell penetration, growth efficiency and cell-to-cell spread, whereas substitution with glutamic acid yielded a virus indistinguishable from parental yellow fever 17D. Substitution with lysine was not tolerated and substitution with asparagine resulted in frequent wild-type revertants. A glycine residue was not tolerated at position 361, but substitution at 362 yielded a small plaque virus, similar to the effect of substitution at position 360. These data indicate that the yellow fever virus E protein contains a locus within domain III where a negative-charge cluster is important for optimal function of this domain in virus-cell interactions beyond the stage of virus attachment. Modelling predictions suggest that the mutations alter the local properties of the loop within domain III, and may compromise interactions of this domain with an adjacent region of domain I during conformational changes that occur in the E protein in association with virus entry.

Critical Regions for Assembly of Vertebrate Nonmuscle Myosin II

Myosin II molecules assemble and form filaments through their C-terminal rod region, and the dynamic filament assembly-disassembly process of nonmuscle myosin II molecules is important for cellular activities. To estimate the critical region for filament formation of vertebrate nonmuscle myosin II, we assessed the solubility of a series of truncated recombinant rod fragments of nonmuscle myosin IIB at various concentrations of NaCl. A C-terminal 248-residue rod fragment (Asp 1729-Glu 1976) was shown by its solubility behavior to retain native assembly features, and two regions within it were found to be necessary for assembly: 35 amino acid residues from Asp 1729 to Thr 1763 and 39 amino acid residues from Ala 1875 to Ala 1913, the latter containing a sequence similar to the assembly competence domain (ACD) of skeletal muscle myosin. Fragments lacking either of the two regions were soluble at any NaCl concentration. We referred to these two regions as nonmuscle myosin ACD1 (nACD1) and nACD2, respectively. In addition, we constructed an alpha-helical coiled-coil model of the rod fragment, and found that a remarkable negative charge cluster (termed N1) and a positive charge cluster (termed P2) were present within nACD1 and nACD2, respectively, besides another positive charge cluster (termed P1) in the amino-terminal vicinity of nACD2. From these results, we propose two major electrostatic interactions that are essential for filament formation of nonmuscle myosin II: the antiparallel interaction between P2 and N1 which is essential for the nucleation step and the parallel interaction between P1 and N1 which is important for the elongation step.

Acidic-rich region of amyloid precursor protein induces glial cell apoptosis

Amyloid precursor protein (APP) has several caspase cleavage sites in its C-terminal cytoplasmic domain and N-terminal extracellular domain. Caspase cleavages of APP at its cytosolic tail may result in releasing the domain and inducing cell death. During apoptosis, the N-terminal domain may also be processed at amino acids 197 and 219 by caspases leading to unmasking of an acidic-rich region (AR). In this study, AR-exposing APP was shown to inhibit cell growth after transfection into RBA-1 astrocytes and BV-2 microglial cells. The recombinant AR from residue 220 to 288 of APP (APP220-288) was produced and its biological activities were analyzed. APP220-288 induced morphological changes, cell death, and DNA fragmentation in BV-2 and RBA-1 cells. However, AR was determined to have no apparent effects in suspension cells, erythroleukemia K562 cells, and Jurkat T cells. The cytotoxicity was depending on negative charge cluster and the apoptotic activity of AR was attributed to the inhibition of cell adhesion. In BV-2 microglial cells, AR significantly stimulated Fas expression, although expressions of the pro-inflammatory cytokine genes were not detected. APP220-288 also induced nitric oxide synthase (iNOS) expression and nitric oxide (NO) production. These findings indicate that the acidic-rich domain of APP may have apoptotic activity due to inhibition of cell adhesion and induction of iNOS and Fas expressions. Moreover, unmasking the apoptosis-induced AR may activate and exacerbate glial cells which in turn lead to further progression of the death program.

Solution Structure of the DFF-C Domain of DFF45/ICAD. A Structural Basis for the Regulation of Apoptotic DNA Fragmentation

DFF45/ICAD has dual functions in the final stage of apoptosis, by acting as both a folding chaperone and a DNase inhibitor of DFF40/CAD. Here, we present the solution structure of the C-terminal domain of DFF45, which is essential for its chaperone-like activity. The structure of this domain (DFF-C) consists of four α helices, which are folded in a novel helix-packing arrangement. The 3D structure reveals a large cluster of negatively charged residues on the molecular surface of DFF-C. This observation suggests that charge complementation plays an important role in the interaction of DFF-C with the positively charged catalytic domain of DFF40, and thus for the chaperone activity of DFF45. The structure of DFF-C also provides a rationale for the loss of the chaperone activity in DFF35, a short isoform of DFF45. Indeed, in DFF35, the amino acid sequence is truncated in the middle of the second α helix constituting the structure of DFF-C, and thus both the hydrophobic core and the cluster of negative charges are disrupted.

A paramagnetic probe to localize residues next to carboxylates on protein surfaces.

It is shown that the paramagnetic properties of lanthanides can be exploited to obtain information on specific parts of a protein surface. Owing to the high affinity of coordinatively unsaturated lanthanide complexes for oxygen donors, carboxylate groups can be used as preferential targets for the interaction. The DO3A ligand is particularly useful in these studies, as it coordinates lanthanides in a heptadentate fashion, leaving two sites available for exogenous donors. A solution of a (15)N-labeled sample protein, calbindin D(9k) (75 residues), was titrated with up to 200% of Gd(III)-DO3A complex, and an inversion recovery (15)N-(1)H HSQC experiment was used to measure the paramagnetic contributions to the longitudinal relaxation rates of the amide protons. Relaxation data were used as distance constraints to estimate the number of interacting complexes and the occupancies of their binding sites. Four preferential interaction sites on the protein surface are found. Inspection of the various carboxylate side chains on the surface of the protein indicates that Gd(III)-DO3A interacts preferentially with carboxylate-rich regions, rather than with isolated carboxylates, suggesting the possibility of chelation of one Gd(III)-DO3A molecule by two carboxylate groups. Gd(III)-DO3A is thus a valuable semi-selective probe for clusters of negative charges on the protein surface.