Pair Of Acidic Residues Research Articles

Characterizing conformational states in GPCR structures using machine learning

G protein-coupled receptors (GPCRs) play a pivotal role in signal transduction and represent attractive targets for drug development. Recent advances in structural biology have provided insights into GPCR conformational states, which are critical for understanding their signaling pathways and facilitating structure-based drug discovery. In this study, we introduce a machine learning approach for conformational state annotation of GPCRs. We represent GPCR conformations as high-dimensional feature vectors, incorporating information about amino acid residue pairs involved in the activation pathway. Using a dataset of GPCR conformations in inactive and active states obtained through molecular dynamics simulations, we trained machine learning models to distinguish between inactive-like and active-like conformations. The developed model provides interpretable predictions and can be used for the large-scale analysis of molecular dynamics trajectories of GPCRs.

Molecular basis for non-invasive diagnostics of cardiac amyloids using bone tracers.

Amyloid diseases including Alzheimer's, Parkinson's and over 30 others are incurable life-threatening disorders caused by abnormal protein deposition as fibrils in various organs. Cardiac amyloidosis is particularly challenging to diagnose and treat. Identification of the fibril-forming protein, which in the heart is usually amyloid transthyretin (ATTR) or amyloid immunoglobulin light chain (AL), is paramount to treatment. A transformative non-invasive diagnostic modality is imaging using technetium-labeled pyrophosphate or diphosphonate bone tracers, 99mTc-PYP/DPD/HMDP. For unknown reasons, these tracers show preferential uptake by ATTR deposits. The tracer-binding moiety is unknown and potentially involves amyloid fibrils and/or amyloid-associated calcific deposits. We propose that, like in the bone, the tracers chelate to surface-bound Ca2+ in amyloid. In high-affinity protein sites, Ca2+ is coordinated by pairs of acidic residues. To identify such residues on amyloids, we harnessed atomic structures of patient-derived cardiac amyloids determined using cryogenic electron microscopy since 2019. These structures help explain why most but not all ATTR deposits uptake 99mTc-PYP/DPD/HMDP radiotracers, while in AL the opposite is true. Moreover, fibril structures help explain greater microcalcification observed in ATTR vs. AL deposits. These findings may aid the diagnostics and therapeutic targeting of cardiac amyloidosis and are relevant to other amyloids.

Resolving Protein Conformational Plasticity and Substrate Binding via Machine Learning.

A long-standing target in elucidating the biomolecular recognition process is the identification of binding-competent conformations of the receptor protein. However, protein conformational plasticity and the stochastic nature of the recognition processes often preclude the assignment of a specific protein conformation to an individual ligand-bound pose. Here, we demonstrate that a computational framework coined as RF-TICA-MD, which integrates an ensemble decision-tree-based Random Forest (RF) machine learning (ML) technique with an unsupervised dimension reduction approach time-structured independent component analysis (TICA), provides an efficient and unambiguous solution toward resolving protein conformational plasticity and the substrate binding process. In particular, we consider multimicrosecond-long molecular dynamics (MD) simulation trajectories of a ligand recognition process in solvent-inaccessible cavities of archetypal proteins T4 lysozyme and cytochrome P450cam. We show that in a scenario in which clear correspondence between protein conformation and binding-competent macrostates could not be obtained via an unsupervised dimension reduction approach, an a priori decision-tree-based supervised classification of the simulated recognition trajectories via RF would help characterize key amino acid residue pairs of the protein that are deemed sensitive for ligand binding. A subsequent unsupervised dimensional reduction of the selected residue pairs via TICA would then delineate a conformational landscape of protein which is able to demarcate ligand-bound poses from unbound ones. The proposed RF-TICA-MD approach is shown to be data agnostic and found to be robust when using other ML-based classification methods such as XGBoost. As a promising spinoff of the protocol, the framework is found to be capable of identifying distal protein locations which would be allosterically important for ligand binding and would characterize their roles in recognition pathways. A Python implementation of a proposed ML workflow is available in GitHub https://github.com/navjeet0211/rf-tica-md.

Enhancing protein contact map prediction accuracy via ensembles of inter-residue distance predictors

Protein contact maps capture coevolutionary interactions between amino acid residue pairs that are spatially within certain proximity threshold. Predicted contact maps are used in many protein related problems that include drug design, protein design, protein function prediction, and protein structure prediction. Contact map prediction has achieved significant progress lately but still further challenges remain with prediction of contacts between residues that are separated in the amino acid residue sequence by large numbers of other residues. In this paper, with experimental results on 5 standard benchmark datasets that include membrane proteins, we show that contact map prediction could be significantly enhanced by using ensembles of various state-of-the-art short distance predictors and then by converting predicted distances into contact probabilities. Our program along with its data is available from https://gitlab.com/mahnewton/ecp.

Randomizing of Oligopeptide Conformations by Nearest Neighbor Interactions between Amino Acid Residues.

Flory’s random coil model assumes that conformational fluctuations of amino acid residues in unfolded poly(oligo)peptides and proteins are uncorrelated (isolated pair hypothesis, IPH). This implies that conformational energies, entropies and solvation free energies are all additive. Nearly 25 years ago, analyses of coil libraries cast some doubt on this notion, in that they revealed that aromatic, but also β-branched side chains, could change the 3J(HNHCα) coupling of their neighbors. Since then, multiple bioinformatical, computational and experimental studies have revealed that conformational propensities of amino acids in unfolded peptides and proteins depend on their nearest neighbors. We used recently reported and newly obtained Ramachandran plots of tetra- and pentapeptides with non-terminal homo- and heterosequences of amino acid residues to quantitatively determine nearest neighbor coupling between them with a Ising type model. Results reveal that, depending on the choice of amino acid residue pairs, nearest neighbor interactions either stabilize or destabilize pairs of polyproline II and β-strand conformations. This leads to a redistribution of population between these conformations and a reduction in conformational entropy. Interactions between residues in polyproline II and turn(helix)-forming conformations seem to be cooperative in most cases, but the respective interaction parameters are subject to large statistical errors.

Coevolutionary analysis reveals a distal amino acid residue pair affecting the catalytic activity of GH5 processive endoglucanase from Bacillus subtilis BS-5.

EG5C-1, processive endoglucanase from Bacillus subtilis, is a typical bifunctional cellulase with endoglucanase and exoglucanase activities. The engineering of processive endoglucanase focuses on the catalytic pocket or carbohydrate-binding module tailoring based on sequence/structure information. Herein, a computational strategy was applied to identify the desired mutants in the enzyme molecule by evolutionary-coupling analysis; subsequently, four residue pairs were selected as evolutionary mutational hotspots. Based on iterative-saturation mutagenesis and subsequent enzymatic activity analysis, a superior mutant K51T/L93T has been identified away from the active center. This variant had increased specific activity from 4170 U/µmol of wild-type (WT) to 5678 U/µmol towards carboxymethyl cellulose-Na and an increase towards the substrate Avicel from 320 U/µmol in WT to 521 U/µmol. In addition, kinetic measurements suggested that superior mutant K51T/L93T had a high substrate affinity (Km ) and a remarkable improvement in catalytic efficiency (kcat /Km ). Furthermore, molecular dynamics simulations revealed that the K51T/L93T mutation altered the spatial conformation at the active site cleft, enhancing the interaction frequency between active site residues and substrate, and improving catalytic efficiency and substrate affinity. The current studies provided some perspectives on the effects of distal residue substitution, which might assist in the engineering of processive endoglucanase or other glycoside hydrolases.

In Silico Analysis of Glucose Oxidase from Aspergillus niger: Potential Cysteine Mutation Sites for Enhancing Protein Stability.

Glucose oxidase (GOx) holds considerable advantages for various applications. Nevertheless, the thermal instability of the enzyme remains a grand challenge, impeding the success in applications outside the well-controlled laboratories, particularly in practical bioelectronics. Many strategies to modify GOx to achieve better thermal stability have been proposed. However, modification of this enzyme by adding extra disulfide bonds is yet to be explored. This work describes the in silico bioengineering of GOx from Aspergillus niger by judiciously analyzing characteristics of disulfide bonds found in the Top8000 protein database, then scanning for amino acid residue pairs that are suitable to be replaced with cysteines in order to establish disulfide bonds. Next, we predicted and assessed the mutant GOx models in terms of disulfide bond quality (bond length and α angles), functional impact by means of residue conservation, and structural impact as indicated by Gibbs free energy. We found eight putative residue pairs that can be engineered to form disulfide bonds. Five of these are located in less conserved regions and, therefore, are unlikely to have a deleterious impact on functionality. Finally, two mutations, Pro149Cys and His158Cys, showed potential for stabilizing the protein structure as confirmed by a structure-based stability analysis tool. The findings in this study highlight the opportunity of using disulfide bond modification as a new alternative technique to enhance the thermal stability of GOx.

Elucidating important structural features for the binding affinity of spike - SARS-CoV-2 neutralizing antibody complexes.

The coronavirus disease 2019 (COVID‐19) has affected the lives of millions of people around the world. In an effort to develop therapeutic interventions and control the pandemic, scientists have isolated several neutralizing antibodies against SARS‐CoV‐2 from the vaccinated and convalescent individuals. These antibodies can be explored further to understand SARS‐CoV‐2 specific antigen–antibody interactions and biophysical parameters related to binding affinity, which can be utilized to engineer more potent antibodies for current and emerging SARS‐CoV‐2 variants. In the present study, we have analyzed the interface between spike protein of SARS‐CoV‐2 and neutralizing antibodies in terms of amino acid residue propensity, pair preference, and atomic interaction energy. We observed that Tyr residues containing contacts are highly preferred and energetically favorable at the interface of spike protein–antibody complexes. We have also developed a regression model to relate the experimental binding affinity for antibodies using structural features, which showed a correlation of 0.93. Moreover, several mutations at the spike protein–antibody interface were identified, which may lead to immune escape (epitope residues) and improved affinity (paratope residues) in current/emerging variants. Overall, the work provides insights into spike protein–antibody interactions, structural parameters related to binding affinity and mutational effects on binding affinity change, which can be helpful to develop better therapeutics against COVID‐19.

Topological analyses of the L-lysine exporter LysO reveal a critical role for a conserved pair of intramembrane solvent-exposed acidic residues

LysO, a prototypical member of the LysO family, mediates export of L-lysine (Lys) and resistance to the toxic Lys antimetabolite, L-thialysine (Thl) in Escherichia coli. Here, we have addressed unknown aspects of LysO function pertaining to its membrane topology and the mechanism by which it mediates Lys/Thl export. Using substituted cysteine (Cys) accessibility, here we delineated the membrane topology of LysO. Our studies support a model in which both the N- and C-termini of LysO are present at the periplasmic face of the membrane with a transmembrane (TM) domain comprising eight TM segments (TMSs) between them. In addition, a feature of intramembrane solvent exposure in LysO is inferred with the identification of membrane-located solvent-exposed Cys residues. Isosteric substitutions of a pair of conserved acidic residues, one E233, located in the solvent-exposed TMS7 and the other D261, in a solvent-exposed intramembrane segment located between TMS7 and TMS8, abolished LysO function in vivo. Thl, but not Lys, elicited proton release in inside-out membrane vesicles, a process requiring the presence of both E233 and D261. We postulate that Thl may be exported in antiport with H+ and that Lys may be a low-affinity export substrate. Our findings are compatible with a physiological scenario wherein in vivo LysO exports the naturally occurring antimetabolite Thl with higher affinity over the essential cellular metabolite Lys, thus affording protection from Thl toxicity and limiting wasteful export of Lys.

Structure of the E. coli agmatinase, SPEB.

Agmatine amidinohydrolase, or agmatinase, catalyzes the conversion of agmatine to putrescine and urea. This enzyme is found broadly across kingdoms of life and plays a critical role in polyamine biosynthesis and the regulation of agmatine concentrations. Here we describe the high-resolution X-ray crystal structure of the E. coli agmatinase, SPEB. The data showed a relatively high degree of pseudomerohedral twinning, was ultimately indexed in the P31 space group and led to a final model with eighteen chains, corresponding to three full hexamers in the asymmetric unit. There was a solvent content of 38.5% and refined R/Rfree values of 0.166/0.216. The protein has the conserved fold characteristic of the agmatine ureohydrolase family and displayed a high degree of structural similarity among individual protomers. Two distinct peaks of electron density were observed in the active site of most of the eighteen chains of SPEB. As the activity of this protein is known to be dependent upon manganese and the fold is similar to other dinuclear metallohydrolases, these peaks were modeled as manganese ions. The orientation of the conserved active site residues, in particular those amino acids that participate in binding the metal ions and a pair of acidic residues (D153 and E274 in SPEB) that play a role in catalysis, are similar to other agmatinase and arginase enzymes and is consistent with a hydrolytic mechanism that proceeds via a metal-activated hydroxide ion.

Elucidating the Mechanism of Membrane Destabilization by the Preferred Modes of Insertion of the Sars-Cov2 Fusion Peptide

Cell penetration of the SARS-Cov2 virus after recognition by the ACE2 receptor requires the fusion of the viral envelope membrane with cellular membranes. The Spike protein (S) of the virus harbors a region identified as the “fusion peptide” (FP) which is liberated at its N-terminal site by a specific cleavage occurring in concert with the interaction of the receptor binding domain of the Spike. Studies have shown that the SARS-Cov2-FP binds Ca2+ ions and perturbs membranes in a calcium-dependent fashion. But the mechanisms of membrane insertion and destabilization remained unclear. We have identified the preferred modes of SARS-Cov2-FP insertion and the role of Ca2+ ions in mediating peptide-membrane interactions from extensive atomistic molecular dynamics (MD) simulations and trajectory analyses. Multiple Ca2+ binding modes in which the ion engaged with different pairs of acidic residues on the FP were observed. A systematic sampling of the interactions of these Ca2+-bound peptide models with lipid membranes showed that SARS-Cov2-FP penetrated the bilayer only in two modes involving different structural domains. In the first mode the hydrophobic residues F833/I834 from the middle region of the peptide are inserted. The second, and more prevalent mode of penetration involves residues L822/F823 in the freed N-terminus of the peptide. The latter insertion was related to the specific mode of Ca2+ association with the peptide in which one Ca2+ is bound to the D830/D839 pair and another to E819/D820. These findings explain the mechanistic role of the necessary cleavage (termed S2’) that frees the N-terminus of the peptide to effectuate the cell-entry of SARS-Cov2. Computational resources provided by the COVID-19 HPC Consortium are gratefully acknowledged.

Transmembrane features governing Fc receptor CD16A assembly with CD16A signaling adaptor molecules

Many activating immunoreceptors associate with signaling adaptor molecules like FcεR1γ or CD247. FcεR1γ and CD247 share high sequence homology and form disulphide-linked homodimers that contain a pair of acidic aspartic acid residues in their transmembrane (TM) domains that mediate assembly, via interaction with an arginine residue at a similar register to these aspartic acids, with the activating immunoreceptors. However, this model cannot hold true for receptors like CD16A, whose TM domains do not contain basic residues. We have carried out an extensive site-directed mutagenesis analysis of the CD16A receptor complex and now report that the association of receptor with the signaling adaptor depends on a network of polar and aromatic residues along the length of the TM domain. Molecular modeling indicates that CD16A TM residues F202, D205, and T206 form the core of the membrane-embedded trimeric interface by establishing highly favorable contacts to the signaling modules through rearrangement of a hydrogen bond network previously identified in the CD247 TM dimer solution NMR structure. Strikingly, the amino acid D205 also regulates the turnover and surface expression of CD16A in the absence of FcεR1γ or CD247. Modeling studies indicate that similar features underlie the association of other activating immune receptors, including CD64 and FcεR1α, with signaling adaptor molecules, and we confirm experimentally that equivalent F, D, and T residues in the TM domain of FcεR1α markedly influence the biology of this receptor and its association with FcεR1γ.

EHD3 Protein Is Required for Tubular Recycling Endosome Stabilization, and an Asparagine-Glutamic Acid Residue Pair within Its Eps15 Homology (EH) Domain Dictates Its Selective Binding to NPF Peptides

An elaborate network of dynamic lipid membranes, termed tubular recycling endosomes (TRE), coordinates the process of endocytic recycling in mammalian cells. The C-terminal Eps15 homology domain (EHD)-containing proteins have been implicated in the bending and fission of TRE, thus regulating endocytic recycling. EHD proteins have an EH domain that interacts with proteins containing an NPF motif. We found that NPF-containing EHD1 interaction partners such as molecules interacting with CasL-like1 (MICAL-L1) and Syndapin2 are essential for TRE biogenesis. Also crucial for TRE biogenesis is the generation of phosphatidic acid, an essential lipid component of TRE that serves as a docking point for MICAL-L1 and Syndapin2. EHD1 and EHD3 have 86% amino acid identity; they homo- and heterodimerize and partially co-localize to TRE. Despite their remarkable identity, they have distinct mechanistic functions. EHD1 induces membrane vesiculation, whereas EHD3 supports TRE biogenesis and/or stabilization by an unknown mechanism. While using phospholipase D inhibitors (which block the conversion of glycerophospholipids to phosphatidic acid) to deplete cellular TRE, we observed that, upon inhibitor washout, there was a rapid and dramatic regeneration of MICAL-L1-marked TRE. Using this "synchronized" TRE biogenesis system, we determined that EHD3 is involved in the stabilization of TRE rather than in their biogenesis. Moreover, we identify the residues Ala-519/Asp-520 of EHD1 and Asn-519/Glu-520 of EHD3 as defining the selectivity of these two paralogs for NPF-containing binding partners, and we present a model to explain the atomic mechanism and provide new insight for their differential roles in vesiculation and tubulation, respectively.

Specificity of broad protein interaction surfaces for proteins with multiple binding partners

Analysis of protein-protein interaction networks has revealed the presence of proteins with multiple interaction ligand proteins, such as hub proteins. For such proteins, multiple ligands would be predicted as interacting partners when predicting all-to-all protein-protein interactions (PPIs). In this work, to obtain a better understanding of PPI mechanisms, we focused on protein interaction surfaces, which differ between protein pairs. We then performed rigid-body docking to obtain information of interfaces of a set of decoy structures, which include many possible interaction surfaces between a certain protein pair. Then, we investigated the specificity of sets of decoy interactions between true binding partners in each case of alpha-chymotrypsin, actin, and cyclin-dependent kinase 2 as test proteins having multiple true binding partners. To observe differences in interaction surfaces of docking decoys, we introduced broad interaction profiles (BIPs), generated by assembling interaction profiles of decoys for each protein pair. After cluster analysis, the specificity of BIPs of true binding partners was observed for each receptor. We used two types of BIPs: those involved in amino acid sequences (BIP-seqs) and those involved in the compositions of interacting amino acid residue pairs (BIP-AAs). The specificity of a BIP was defined as the number of group members including all true binding partners. We found that BIP-AA cases were more specific than BIP-seq cases. These results indicated that the composition of interacting amino acid residue pairs was sufficient for determining the properties of protein interaction surfaces.

Modeling of interaction between cytochrome c and the WD domains of Apaf-1: bifurcated salt bridges underlying apoptosome assembly.

BackgroundBinding of cytochrome c, released from the damaged mitochondria, to the apoptotic protease activating factor 1 (Apaf-1) is a key event in the apoptotic signaling cascade. The binding triggers a major domain rearrangement in Apaf-1, which leads to oligomerization of Apaf-1/cytochrome c complexes into an apoptosome. Despite the availability of crystal structures of cytochrome c and Apaf-1 and cryo-electron microscopy models of the entire apoptosome, the binding mode of cytochrome c to Apaf-1, as well as the nature of the amino acid residues of Apaf-1 involved remain obscure.ResultsWe investigated the interaction between cytochrome c and Apaf-1 by combining several modeling approaches. We have applied protein-protein docking and energy minimization, evaluated the resulting models of the Apaf-1/cytochrome c complex, and carried out a further analysis by means of molecular dynamics simulations. We ended up with a single model structure where all the lysine residues of cytochrome c that are known as functionally-relevant were involved in forming salt bridges with acidic residues of Apaf-1. This model has revealed three distinctive bifurcated salt bridges, each involving a single lysine residue of cytochrome c and two neighboring acidic resides of Apaf-1. Salt bridge-forming amino acids of Apaf-1 showed a clear evolutionary pattern within Metazoa, with pairs of acidic residues of Apaf-1, involved in bifurcated salt bridges, reaching their highest numbers in the sequences of vertebrates, in which the cytochrome c-mediated mechanism of apoptosome formation seems to be typical.ConclusionsThe reported model of an Apaf-1/cytochrome c complex provides insights in the nature of protein-protein interactions which are hard to observe in crystallographic or electron microscopy studies. Bifurcated salt bridges can be expected to be stronger than simple salt bridges, and their formation might promote the conformational change of Apaf-1, leading to the formation of an apoptosome. Combination of structural and sequence analyses provides hints on the evolution of the cytochrome c-mediated apoptosis.ReviewersThis article was reviewed by Andrei L. Osterman, Narayanaswamy Srinivasan, Igor N. Berezovsky, and Gerrit Vriend (nominated by Martijn Huynen).Electronic supplementary materialThe online version of this article (doi:10.1186/s13062-015-0059-4) contains supplementary material, which is available to authorized users.

External pH modulates EAG superfamily K+ channels through EAG-specific acidic residues in the voltage sensor

The Ether-a-go-go (EAG) superfamily of voltage-gated K+ channels consists of three functionally distinct gene families (Eag, Elk, and Erg) encoding a diverse set of low-threshold K+ currents that regulate excitability in neurons and muscle. Previous studies indicate that external acidification inhibits activation of three EAG superfamily K+ channels, Kv10.1 (Eag1), Kv11.1 (Erg1), and Kv12.1 (Elk1). We show here that Kv10.2, Kv12.2, and Kv12.3 are similarly inhibited by external protons, suggesting that high sensitivity to physiological pH changes is a general property of EAG superfamily channels. External acidification depolarizes the conductance–voltage (GV) curves of these channels, reducing low threshold activation. We explored the mechanism of this high pH sensitivity in Kv12.1, Kv10.2, and Kv11.1. We first examined the role of acidic voltage sensor residues that mediate divalent cation block of voltage activation in EAG superfamily channels because protons reduce the sensitivity of Kv12.1 to Zn2+. Low pH similarly reduces Mg2+ sensitivity of Kv10.1, and we found that the pH sensitivity of Kv11.1 was greatly attenuated at 1 mM Ca2+. Individual neutralizations of a pair of EAG-specific acidic residues that have previously been implicated in divalent block of diverse EAG superfamily channels greatly reduced the pH response in Kv12.1, Kv10.2, and Kv11.1. Our results therefore suggest a common mechanism for pH-sensitive voltage activation in EAG superfamily channels. The EAG-specific acidic residues may form the proton-binding site or alternatively are required to hold the voltage sensor in a pH-sensitive conformation. The high pH sensitivity of EAG superfamily channels suggests that they could contribute to pH-sensitive K+ currents observed in vivo.

Prediction of conformational epitopes with the use of a knowledge-based energy function and geometrically related neighboring residue characteristics

BackgroundA conformational epitope (CE) in an antigentic protein is composed of amino acid residues that are spatially near each other on the antigen's surface but are separated in sequence; CEs bind their complementary paratopes in B-cell receptors and/or antibodies. CE predication is used during vaccine design and in immuno-biological experiments. Here, we develop a novel system, CE-KEG, which predicts CEs based on knowledge-based energy and geometrical neighboring residue contents. The workflow applied grid-based mathematical morphological algorithms to efficiently detect the surface atoms of the antigens. After extracting surface residues, we ranked CE candidate residues first according to their local average energy distributions. Then, the frequencies at which geometrically related neighboring residue combinations in the potential CEs occurred were incorporated into our workflow, and the weighted combinations of the average energies and neighboring residue frequencies were used to assess the sensitivity, accuracy, and efficiency of our prediction workflow.ResultsWe prepared a database containing 247 antigen structures and a second database containing the 163 non-redundant antigen structures in the first database to test our workflow. Our predictive workflow performed better than did algorithms found in the literature in terms of accuracy and efficiency. For the non-redundant dataset tested, our workflow achieved an average of 47.8% sensitivity, 84.3% specificity, and 80.7% accuracy according to a 10-fold cross-validation mechanism, and the performance was evaluated under providing top three predicted CE candidates for each antigen.ConclusionsOur method combines an energy profile for surface residues with the frequency that each geometrically related amino acid residue pair occurs to identify possible CEs in antigens. This combination of these features facilitates improved identification for immuno-biological studies and synthetic vaccine design. CE-KEG is available at http://cekeg.cs.ntou.edu.tw.

Bulky Substitutions in the Thumb Domain of ASIC1a Restrict Proton-Gating

Acid-sensing ion channels (ASICs) are trimeric proton-gated cation selective ion channels expressed in the nervous system. Each monomer consists of two transmembrane helices with short cytoplasmic tails and a large extracellular region organized in distinct domains named the thumb, finger, knuckle, β-ball, and palm. In this report we investigated the contribution to mouse ASIC1a activation and desensitization of three pairs of acidic residues, Glu238-Asp345, Asp237-Asp349 and Glu219-Asp407. Residues in each pair are located in close proximity and together form an acidic cluster at the interface between the thumb and β-ball domains. We examined the reactivity toward methanethiosulfonate (MTS) reagents of channels carrying individual Cys substitutions at these positions. MTSET-modification of D345C channels shifted the pH of half-maximal activation (pH50) from 6.23±0.06 to 4.48±0.08, and increased the desensitization rate. A similar effect was observed with other positively-charged MTS reagents such as MTSMT, MTSPT and MTSEA, as well as with MTSES, a negatively-charged reagent. Moreover, Lys substitution at position 345 shifted pH50 toward more acidic values. Remarkably, Lys substitution at positions 219, 237 or 238 did not significantly change the pH50, suggesting that acid residues at these positions do not substantially contribute to proton activation. Our results indicate that the introduction of a bulky moiety at position 345 restricts conformational movements required for proton-gating, suggesting that the thumb-β-ball interface undergoes dynamic rearrangements during activation.

Dissecting a regulatory calcium-binding site of CLC-K kidney chloride channels

The kidney and inner ear CLC-K chloride channels, which are involved in salt absorption and endolymph production, are regulated by extracellular Ca2+ in the millimolar concentration range. Recently, Gradogna et al. (2010. J. Gen. Physiol. http://dx.doi.org/10.1085/jgp.201010455) identified a pair of acidic residues (E261 and D278) located in the loop between helices I and J as forming a putative intersubunit Ca2+-binding site in hClC-Ka. In this study, we sought to explore the properties of the binding site in more detail. First, we verified that the site is conserved in hClC-Kb and rClC-K1. In addition, we could confer Ca2+ sensitivity to the Torpedo marmorata ClC-0 channel by exchanging its I–J loop with that from ClC-Ka, demonstrating a direct role of the loop in Ca2+ binding. Based on a structure of a bacterial CLC and a new sequence alignment, we built homology models of ClC-Ka. The models suggested additional amino acids involved in Ca2+ binding. Testing mutants of these residues, we could restrict the range of plausible models and positively identify two more residues (E259 and E281) involved in Ca2+ coordination. To investigate cation specificity, we applied extracellular Zn2+, Mg2+, Ba2+, Sr2+, and Mn2+. Zn2+ blocks ClC-Ka as well as its Ca2+-insensitive mutant, suggesting that Zn2+ binds to a different site. Mg2+ does not activate CLC-Ks, but the channels are activated by Ba2+, Sr2+, and Mn2+ with a rank order of potency of Ca2+ > Ba2+ > Sr2+ = Mn2+ for the human CLC-Ks. Dose–response analysis indicates that the less potent Ba2+ has a lower affinity rather than a lower efficacy. Interestingly, rClC-K1 shows an altered rank order (Ca2+ > Sr2+ >> Ba2+), but homology models suggest that residues outside the I–J loop are responsible for this difference. Our detailed characterization of the regulatory Ca2+-binding site provides a solid basis for the understanding of the physiological modulation of CLC-K channel function in the kidney and inner ear.

Identification of Key Residues Essential for the Structural Fold and Receptor Selectivity within the A-chain of Human Gene-2 (H2) Relaxin

Human gene-2 (H2) relaxin is currently in Phase III clinical trials for the treatment of acute heart failure. It is a 53-amino acid insulin-like peptide comprising two chains and three disulfide bonds. It interacts with two of the relaxin family peptide (RXFP) receptors. Although its cognate receptor is RXFP1, it is also able to cross-react with RXFP2, the native receptor for a related peptide, insulin-like peptide 3. In order to understand the basis of this cross-reactivity, it is important to elucidate both binding and activation mechanisms of this peptide. The primary binding mechanism of this hormone has been extensively studied and well defined. H2 relaxin binds to the leucine-rich repeats of RXFP1 and RXFP2 using B-chain-specific residues. However, little is known about the secondary interaction that involves the A-chain of H2 relaxin and transmembrane exoloops of the receptors. We demonstrate here through extensive mutation of the A-chain that the secondary interaction between H2 relaxin and RXFP1 is not driven by any single amino acid, although residues Tyr-3, Leu-20, and Phe-23 appear to contribute. Interestingly, these same three residues are important drivers of the affinity and activity of H2 relaxin for RXFP2 with additional minor contributions from Lys-9, His-12, Lys-17, Arg-18, and Arg-22. Our results provide new insights into the mechanism of secondary activation interaction of RXFP1 and RXFP2 by H2 relaxin, leading to a potent and RXFP1-selective analog, H2:A(4-24)(F23A), which was tested in vitro and in vivo and found to significantly inhibit collagen deposition similar to native H2 relaxin.