Conformational Changes of Blood ACE in Chronic Uremia

Background

The pattern of binding of monoclonal antibodies (mAbs) to 16 epitopes on human angiotensin I-converting enzyme (ACE) comprise a conformational ACE fingerprint and is a sensitive marker of subtle protein conformational changes.

Hypothesis

Toxic substances in the blood of patients with uremia due to End Stage Renal Disease (ESRD) can induce local conformational changes in the ACE protein globule and alter the efficacy of ACE inhibitors.

Methodology/Principal Findings

The recognition of ACE by 16 mAbs to the epitopes on the N and C domains of ACE was estimated using an immune-capture enzymatic plate precipitation assay. The precipitation pattern of blood ACE by a set of mAbs was substantially influenced by the presence of ACE inhibitors with the most dramatic local conformational change noted in the N-domain region recognized by mAb 1G12. The “short” ACE inhibitor enalaprilat (tripeptide analog) and “long” inhibitor teprotide (nonapeptide) produced strikingly different mAb 1G12 binding with enalaprilat strongly increasing mAb 1G12 binding and teprotide decreasing binding. Reduction in S-S bonds via glutathione and dithiothreitol treatment increased 1G12 binding to blood ACE in a manner comparable to enalaprilat. Some patients with uremia due to ESRD exhibited significantly increased mAb 1G12 binding to blood ACE and increased ACE activity towards angiotensin I accompanied by reduced ACE inhibition by inhibitory mAbs and ACE inhibitors.

Conclusions/Significance

The estimation of relative mAb 1G12 binding to blood ACE detects a subpopulation of ESRD patients with conformationally changed ACE, which activity is less suppressible by ACE inhibitors. This parameter may potentially serve as a biomarker for those patients who may need higher concentrations of ACE inhibitors upon anti-hypertensive therapy.     

Evolution of the class C GPCR Venus flytrap modules involved positive selected functional divergence

Background  

Class C G protein-coupled receptors (GPCRs) represent a distinct group of the GPCR family, which structurally possess a characteristically distinct extracellular domain inclusive of the Venus flytrap module (VFTM). The VFTMs of the class C GPCRs is responsible for ligand recognition and binding, and share sequence similarity with bacterial periplasmic amino acid binding proteins (PBPs). An extensive phylogenetic investigation of the VFTMs was conducted by analyzing for functional divergence and testing for positive selection for five typical groups of the class C GPCRs. The altered selective constraints were determined to identify the sites that had undergone functional divergence via positive selection. In order to structurally demonstrate the pattern changes during the evolutionary process, three-dimensional (3D) structures of the GPCR VFTMs were modelled and reconstructed from ancestral VFTMs.     

Statistical analysis and molecular dynamics simulations of ambivalent <Emphasis Type="Italic">α</Emphasis> -helices

Background  

Analysis of known protein structures reveals that identical sequence fragments in proteins can adopt different secondary structure conformations. The extent of this conformational diversity is influenced by various factors like the intrinsic sequence propensity, sequence context and other environmental factors such as pH, site directed mutations or alteration of the binding ligands. Understanding the mechanism by which the environment affects the structural ambivalence of these peptides has potential implications for protein design and reliable local structure prediction algorithms. Identification of the structurally ambivalent sequence fragments and determining the rules which dictate their conformational preferences play an important role in understanding the conformational changes observed in misfolding diseases. However, a systematic classification of their intrinsic sequence patterns or a statistical analysis of their properties and sequence context in relation to the origin of their structural diversity have largely remained unexplored.     

A previously unobserved conformation for the human Pex5p receptor suggests roles for intrinsic flexibility and rigid domain motions in ligand binding

Background  

The C-terminal tetratricopeptide (TPR) repeat domain of Pex5p recognises proteins carrying a peroxisomal targeting signal type 1 (PTS1) tripeptide in their C-terminus. Previously, structural data have been obtained from the TPR domain of Pex5p in both the liganded and unliganded states, indicating a conformational change taking place upon cargo protein binding. Such a conformational change would be expected to play a major role both during PTS1 protein recognition as well as in cargo release into the peroxisomal lumen. However, little information is available on the factors that may regulate such structural changes.     

Integrated prediction of one-dimensional structural features and their relationships with conformational flexibility in helical membrane proteins

Background  

Many structural properties such as solvent accessibility, dihedral angles and helix-helix contacts can be assigned to each residue in a membrane protein. Independent studies exist on the analysis and sequence-based prediction of some of these so-called one-dimensional features. However, there is little explanation of why certain residues are predicted in a wrong structural class or with large errors in the absolute values of these features. On the other hand, membrane proteins undergo conformational changes to allow transport as well as ligand binding. These conformational changes often occur via residues that are inherently flexible and hence, predicting fluctuations in residue positions is of great significance.     

Molecular modeling studies give hint for the existence of a symmetric hβ2R-Gαβγ-homodimer

Several experimental studies suggest that GPCR dimers or oligomers may play an important role in signal transduction. In 2011 the crystal structure of a hβ₂R-Gαβγ-complex was published and crystal structures of GPCR dimers are known. But until now, no crystal structure of a GPCR dimer including the Gαβγ-complex is available. In order to obtain detailed insights into interactions within hβ₂R dimers including the Gαβγ-complex we performed a potential-energy-surface scan in order to identify favored asymmetric and symmetric hβ₂R-Gαβγ-homodimers. This potential energy surface scan suggests, besides the existence of asymmetric dimers, the existence of a symmetric hβ₂R-Gαβγ-homodimer with a TM I/VII-contact. A subsequent 20 ns MD simulation of the symmetric homodimer revealed large asymmetric conformational changes of both hβ₂Rs, especially regarding TM VII and the interaction network between Asp^2.50, Val^7.44, Ser^7.46 and Tyr^7.43. Since similar conformational changes were not observed during the molecular dynamic simulation of the monomeric hβ₂R-Gαβγ-complex, it may be suggested that the conformational changes in the symmetric homodimer are related to the presence of the second hβ₂R-Gαβγ-complex. Due to the limitations of simulation time, conformational changes within a time scale of μs or ms may of course not be observed. However, the detected conformational changes, especially in TM VII, correspond to minima on the potential energy surface and thus, this study gives new insights into GPCR dimers on molecular level and furthermore, gives suggestions for site-directed mutagenesis studies.

Conformational mobility of active and E-64-inhibited actinidin

Background

Actinidin, a protease from kiwifruit, belongs to the C1 family of cysteine proteases. Cysteine proteases were found to be involved in many disease states and are valid therapeutic targets. Actinidin has a wide pH activity range and wide substrate specificity, which makes it a good model system for studying enzyme–substrate interactions.

Methods

The influence of inhibitor (E-64) binding on the conformation of actinidin was examined by 2D PAGE, circular dichroism (CD) spectroscopy, hydrophobic ligand binding assay, and molecular dynamics simulations.

Results

Significant differences were observed in electrophoretic mobility of proteolytically active and E-64-inhibited actinidin. CD spectrometry and hydrophobic ligand binding assay revealed a difference in conformation between active and inhibited actinidin. Molecular dynamics simulations showed that a loop defined by amino-acid residues 88–104 had greater conformational mobility in the inhibited enzyme than in the active one. During MD simulations, the covalently bound inhibitor was found to change its conformation from extended to folded, with the guanidino moiety approaching the carboxylate.

Conclusions

Conformational mobility of actinidin changes upon binding of the inhibitor, leading to a sequence of events that enables water and ions to protrude into a newly formed cavity of the inhibited enzyme. Drastic conformational mobility of E-64, a common inhibitor of cysteine proteases found in many crystal structures stored in PDB, was also observed.

General significance

The analysis of structural changes which occur upon binding of an inhibitor to a cysteine protease provides a valuable starting point for the future design of therapeutic agents.     

A Monomeric Variant of Insulin Degrading Enzyme (IDE) Loses Its Regulatory Properties

Background

Insulin degrading enzyme (IDE) is a key enzyme in the metabolism of both insulin and amyloid beta peptides. IDE is unique in that it is subject to allosteric activation which is hypothesized to occur through an oligomeric structuture.

Methodology/Principal Findings

IDE is known to exist as an equilibrium mixture of monomers, dimers, and higher oligomers, with the dimer being the predominant form. Based on the crystal structure of IDE we deleted the putative dimer interface in the C-terminal region, which resulted in a monomeric variant. Monomeric IDE retained enzymatic activity, however instead of the allosteric behavior seen with wild type enzyme it displayed Michaelis-Menten kinetic behavior. With the substrate Abz-GGFLRKHGQ-EDDnp, monomeric IDE retained ∼25% of the wild type activity. In contrast with the larger peptide substrates β-endorphin and amyloid β peptide 1–40, monomeric IDE retained only 1 to 0.25% of wild type activity. Unlike wild type IDE neither bradykinin nor dynorphin B-9 activated the monomeric variant of the enzyme. Similarly, monomeric IDE was not activated by polyphosphates under conditions in which the activity of wild type enzyme was increased more than 50 fold.

Conclusions/Significance

These findings serve to establish the dimer interface in IDE and demonstrate the requirement for an oligomeric form of the enzyme for its regulatory properties. The data support a mechanism where the binding of activators to oligomeric IDE induces a conformational change that cannot occur in the monomeric variant. Since a conformational change from a closed to a more open structure is likely the rate-determining step in the IDE reaction, the subunit induced conformational change likely shifts the structure of the oligomeric enzyme to a more open conformation.     

Designing coarse grained-and atom based-potentials for protein-protein docking

Background  

Protein-protein docking is a challenging computational problem in functional genomics, particularly when one or both proteins undergo conformational change(s) upon binding. The major challenge is to define a scoring function soft enough to tolerate these changes and specific enough to distinguish between near-native and "misdocked" conformations.     

Cytochalasin B triggers a novel pertussis toxin sensitive pathway in TNF-alpha primed neutrophils

Background  

Cytochalasin B does not directly activate the oxygen-radical-producing NADPH oxidase activity of neutrophils but transfers desensitized G-protein coupled receptors (GPCR) into an active signaling state by uncoupling GCPR from the cytoskeleton. The receptor uncoupling results in respiratory burst activity when signals generated by reactivated formyl peptide receptors trigger the NADPH-oxidase to produce superoxide anions.     

Riboswitch Detection Using Profile Hidden Markov Models

Background  

Riboswitches are a type of noncoding RNA that regulate gene expression by switching from one structural conformation to another on ligand binding. The various classes of riboswitches discovered so far are differentiated by the ligand, which on binding induces a conformational switch. Every class of riboswitch is characterized by an aptamer domain, which provides the site for ligand binding, and an expression platform that undergoes conformational change on ligand binding. The sequence and structure of the aptamer domain is highly conserved in riboswitches belonging to the same class. We propose a method for fast and accurate identification of riboswitches using profile Hidden Markov Models (pHMM). Our method exploits the high degree of sequence conservation that characterizes the aptamer domain.     

Sparse estimation for structural variability

Background  

Proteins are dynamic molecules that exhibit a wide range of motions; often these conformational changes are important for protein function. Determining biologically relevant conformational changes, or true variability, efficiently is challenging due to the noise present in structure data.     

Forced unbinding of GPR17 ligands from wild type and R255I mutant receptor models through a computational approach

Background  

GPR17 is a hybrid G-protein-coupled receptor (GPCR) activated by two unrelated ligand families, extracellular nucleotides and cysteinyl-leukotrienes (cysteinyl-LTs), and involved in brain damage and repair. Its exploitment as a target for novel neuro-reparative strategies depends on the elucidation of the molecular determinants driving binding of purinergic and leukotrienic ligands. Here, we applied docking and molecular dynamics simulations (MD) to analyse the binding and the forced unbinding of two GPR17 ligands (the endogenous purinergic agonist UDP and the leukotriene receptor antagonist pranlukast from both the wild-type (WT) receptor and a mutant model, where a basic residue hypothesized to be crucial for nucleotide binding had been mutated (R255I) to Ile.     

Conserved protein domains are maintained in an average ratio to proteome size

Background  

Conserved domains (CD) in proteins play a crucial role in protein interactions, DNA binding, enzyme activity, and other important cellular processes. We proposed to study ratios of genes containing these domains to ratios of proteome size of different eukaryotes.     

Identifying allosteric fluctuation transitions between different protein conformational states as applied to Cyclin Dependent Kinase 2

Background  

The mechanisms underlying protein function and associated conformational change are dominated by a series of local entropy fluctuations affecting the global structure yet are mediated by only a few key residues. Transitional Dynamic Analysis (TDA) is a new method to detect these changes in local protein flexibility between different conformations arising from, for example, ligand binding. Additionally, Positional Impact Vertex for Entropy Transfer (PIVET) uses TDA to identify important residue contact changes that have a large impact on global fluctuation. We demonstrate the utility of these methods for Cyclin-dependent kinase 2 (CDK2), a system with crystal structures of this protein in multiple functionally relevant conformations and experimental data revealing the importance of local fluctuation changes for protein function.     

Domain motions of Argonaute,the catalytic engine of RNA interference

Background  

The Argonaute protein is the core component of the RNA-induced silencing complex, playing the central role of cleaving the mRNA target. Visual inspection of static crystal structures already has enabled researchers to suggest conformational changes of Argonaute that might occur during RNA interference. We have taken the next step by performing an all-atom normal mode analysis of the Pyrococcus furiosus and Aquifex aeolicus Argonaute crystal structures, allowing us to quantitatively assess the feasibility of these conformational changes. To perform the analysis, we begin with the energy-minimized X-ray structures. Normal modes are then calculated using an all-atom molecular mechanics force field.     

Pep-3D-Search: a method for B-cell epitope prediction based on mimotope analysis

Background  

The prediction of conformational B-cell epitopes is one of the most important goals in immunoinformatics. The solution to this problem, even if approximate, would help in designing experiments to precisely map the residues of interaction between an antigen and an antibody. Consequently, this area of research has received considerable attention from immunologists, structural biologists and computational biologists. Phage-displayed random peptide libraries are powerful tools used to obtain mimotopes that are selected by binding to a given monoclonal antibody (mAb) in a similar way to the native epitope. These mimotopes can be considered as functional epitope mimics. Mimotope analysis based methods can predict not only linear but also conformational epitopes and this has been the focus of much research in recent years. Though some algorithms based on mimotope analysis have been proposed, the precise localization of the interaction site mimicked by the mimotopes is still a challenging task.     

Analysis of protein-protein docking decoys using interaction fingerprints: application to the reconstruction of CaM-ligand complexes

Background  

Protein-protein docking for proteins with large conformational changes was analyzed by using interaction fingerprints, one of the scales for measuring similarities among complex structures, utilized especially for searching near-native protein-ligand or protein-protein complex structures. Here, we have proposed a combined method for analyzing protein-protein docking by taking large conformational changes into consideration. This combined method consists of ensemble soft docking with multiple protein structures, refinement of complexes, and cluster analysis using interaction fingerprints and energy profiles.     

High-throughput kinase assays with protein substrates using fluorescent polymer superquenching

Background  

High-throughput screening is used by the pharmaceutical industry for identifying lead compounds that interact with targets of pharmacological interest. Because of the key role that aberrant regulation of protein phosphorylation plays in diseases such as cancer, diabetes and hypertension, kinases have become one of the main drug targets. With the exception of antibody-based assays, methods to screen for specific kinase activity are generally restricted to the use of small synthetic peptides as substrates. However, the use of natural protein substrates has the advantage that potential inhibitors can be detected that affect enzyme activity by binding to a site other than the catalytic site. We have previously reported a non-radioactive and non-antibody-based fluorescence quench assay for detection of phosphorylation or dephosphorylation using synthetic peptide substrates. The aim of this work is to develop an assay for detection of phosphorylation of chemically unmodified proteins based on this polymer superquenching platform.