Direct protein interaction underlies gene-for-gene specificity and coevolution of the flax resistance genes and flax rust avirulence genes

Plant resistance proteins (R proteins) recognize corresponding pathogen avirulence (Avr) proteins either indirectly through detection of changes in their host protein targets or through direct R-Avr protein interaction. Although indirect recognition imposes selection against Avr effector function, pathogen effector molecules recognized through direct interaction may overcome resistance through sequence diversification rather than loss of function. Here we show that the flax rust fungus AvrLS67 genes, whose products are recognized by the L5, L6, and L7 R proteins of flax, are highly diverse, with 12 sequence variants identified from six rust strains. Seven AvrL567 variants derived from Avr alleles induce necrotic responses when expressed in flax plants containing corresponding resistance genes (R genes), whereas five variants from avr alleles do not. Differences in recognition specificity between AvA567 variants and evidence for diversifying selection acting on these genes suggest they have been involved in a gene-specific arms race with the corresponding flax R genes. Yeast two-hybrid assays indicate that recognition is based on direct R-Avr protein interaction and recapitulate the interaction specificity observed in planta. Biochemical analysis of Escherichia coli-produced AvrL567 proteins shows that variants that escape recognition nevertheless maintain a conserved structure and stability, suggesting that the amino acid sequence differences directly affect the R-Avr protein interaction. We suggest that direct recognition associated with high genetic diversity at corresponding R and Avr gene loci represents an alternative outcome of plant-pathogen coevolution to indirect recognition associated with simple balanced polymorphisms for functional and nonfunctional R and Avr genes.

Characterization of the structure of RAMP1 by mutagenesis and molecular modeling

Receptor activity modifying proteins (RAMPs) are a family of single-pass transmembrane proteins that dimerize with G-protein-coupled receptors. They may alter the ligand recognition properties of the receptors (particularly for the calcitonin receptor-like receptor, CLR). Very little structural information is available about RAMPs. Here, an ab initio model has been generated for the extracellular domain of RAMP1. The disulfide bond arrangement (Cys 27-Cys82, Cys40-Cys72, and Cys 57-Cys104) was determined by site-directed mutagenesis. The secondary structure (a-helices from residues 29-51, 60-80, and 87-100) was established from a consensus of predictive routines. Using these constraints, an assemblage of 25,000 structures was constructed and these were ranked using an all-atom statistical potential. The best 1000 conformations were energy minimized. The lowest scoring model was refined by molecular dynamics simulation. To validate our strategy, the same methods were applied to three proteins of known structure; PDB:1HP8, PDB:1V54 chain H (residues 21-85), and PDB:1T0P. When compared to the crystal structures, the models had root mean-square deviations of 3.8 Å, 4.1 Å, and 4.0 Å, respectively. The model of RAMP1 suggested that Phe93, Tyr 100, and Phe101 form a binding interface for CLR, whereas Trp74 and Phe92 may interact with ligands that bind to the CLR/RAMP1 heterodimer. © 2006 by the Biophysical Society.

Protein deposition at polymer surfaces:the bulk and surface structure-property effects of hydrogels

The primary objective of this research has been to investigate the interfacial phenomenon of protein adsorption in relation to the bulk and surface structure-property effect s of hydrogel polymers. In order to achieve this it was first necessary to characterise the bulk and surface properties of the hydrogels, with regard to the structural chemistry of their component monomers. The bulk properties of the hydrogels were established using equilibrium water content measurements, together with water-binding studies by differential scanning calorimetry (D.S.C.). Hamilton and captive air bubble-contact angle techniques were employed to characterise the hydrogel-water interface and from which by a mathematical derivation, the interfacial free energy (ðsw) and the surface free energy components (ð psv, ðdsv, ðsv) were obtained. From the adsorption studies using the radio labelled iodinated (125I) proteins of human serum albumin (H.S.A.) and human fibrinogen (H.Fb.), it was Found that multi-layered adsorption was occurring and that the rate and type of this adsorption was dependent on the physico-chemical behaviour of the adsorbing protein (and its bulk concentration in solution), together with the surface energetics of the adsorbent polymer. A potential method for the invitro evaluation of a material's 'biocompatibility' was also investigated, based on an empirically observed relationship between the adsorption of albumin and fibrinogen and the 'biocompatibility' of polymeric materials. Furthermore, some consideration was also given to the biocompatibility problem of proteinaceous deposit formation on hydrophilic soft' contact lenses and in addition a number of potential continual wear contact lens formulations now undergoing clinical trials,were characterised by the above techniques.

Novel peptide antagonists of adrenomedullin and calcitonin gene-related peptide receptors:identification, pharmacological characterization, and interactions with position 74 in receptor activity-modifying protein 1/3

Human adrenomedullin (AM) is a 52-amino acid peptide belonging to the calcitonin peptide family, which also includes calcitonin gene-related peptide (CGRP) and AM2. The two AM receptors, AM(1) and AM(2), are calcitonin receptor-like receptor (CL)/receptor activity-modifying protein (RAMP) (RAMP2 and RAMP3, respectively) heterodimers. CGRP receptors comprise CL/RAMP1. The only human AM receptor antagonist (AM(22-52)) is a truncated form of AM; it has low affinity and is only weakly selective for AM(1) over AM(2) receptors. To develop novel AM receptor antagonists, we explored the importance of different regions of AM in interactions with AM(1), AM(2), and CGRP receptors. AM(22-52) was the framework for generating further AM fragments (AM(26-52) and AM(30-52)), novel AM/alphaCGRP chimeras (C1-C5 and C9), and AM/AM(2) chimeras (C6-C8). cAMP assays were used to screen the antagonists at all receptors to determine their affinity and selectivity. Circular dichroism spectroscopy was used to investigate the secondary structures of AM and its related peptides. The data indicate that the structures of AM, AM2, and alphaCGRP differ from one another. Our chimeric approach enabled the identification of two nonselective high-affinity antagonists of AM(1), AM(2), and CGRP receptors (C2 and C6), one high-affinity antagonist of AM(2) receptors (C7), and a weak antagonist selective for the CGRP receptor (C5). By use of receptor mutagenesis, we also determined that the C-terminal nine amino acids of AM seem to be responsible for its interaction with Glu74 of RAMP3. We provide new information on the structure-activity relationship of AM, alphaCGRP, and AM2 and how AM interacts with CGRP and AM(2) receptors.

G-protein-coupled receptor solubilization and purification for biophysical analysis and functional studies, in the total absence of detergent

G-protein coupled receptors (GPCRs) constitute the largest class of membrane proteins and are a major drug target. A serious obstacle to studying GPCR structure/function characteristics is the requirement to extract the receptors from their native environment in the plasma membrane, coupled with the inherent instability of GPCRs in the detergents required for their solubilization. In the present study, we report the first solubilization and purification of a functional GPCR [human adenosine A2A receptor (A2AR)], in the total absence of detergent at any stage, by exploiting spontaneous encapsulation by styrene maleic acid (SMA) co-polymer direct from the membrane into a nanoscale SMA lipid particle (SMALP). Furthermore, the A2AR-SMALP, generated from yeast (Pichia pastoris) or mammalian cells, exhibited increased thermostability (∼5°C) compared with detergent [DDM (n-dodecyl-β-D-maltopyranoside)]-solubilized A2AR controls. The A2AR-SMALP was also stable when stored for prolonged periods at 4°C and was resistant to multiple freeze-thaw cycles, in marked contrast with the detergent-solubilized receptor. These properties establish the potential for using GPCR-SMALP in receptor-based drug discovery assays. Moreover, in contrast with nanodiscs stabilized by scaffold proteins, the non-proteinaceous nature of the SMA polymer allowed unobscured biophysical characterization of the embedded receptor. Consequently, CD spectroscopy was used to relate changes in secondary structure to loss of ligand binding ([³H]ZM241385) capability. SMALP-solubilization of GPCRs, retaining the annular lipid environment, will enable a wide range of therapeutic targets to be prepared in native-like state to aid drug discovery and understanding of GPCR molecular mechanisms.

The redoxomics of PTEN: walking a fine line between damage and signaling:mass spectrometry-based approaches to study the effect of oxidation on PTEN function, structure, and protein-protein interactions

The research described in this PhD thesis focuses on proteomics approaches to study the effect of oxidation on the modification status and protein-protein interactions of PTEN, a redox-sensitive phosphatase involved in a number of cellular processes including metabolism, apoptosis, cell proliferation, and survival. While direct evidence of a redox regulation of PTEN and its downstream signaling has been reported, the effect of cellular oxidative stress or direct PTEN oxidation on PTEN structure and interactome is still poorly defined. In a first study, GST-tagged PTEN was directly oxidized over a range of hypochlorous acid (HOCl) concentration, assayed for phosphatase activity, and oxidative post-translational modifications (oxPTMs) were quantified using LC-MS/MS-based label-free methods. In a second study, GSTtagged PTEN was prepared in a reduced and reversibly H2O2-oxidized form, immobilized on a resin support and incubated with HCT116 cell lysate to capture PTEN interacting proteins, which were analyzed by LC-MS/MS and comparatively quantified using label-free methods. In parallel experiments, HCT116 cells transfected with a GFP-tagged PTEN were treated with H2O2 and PTENinteracting proteins immunoprecipitated using standard methods. Several high abundance HOCl-induced oxPTMs were mapped, including those taking place at amino acids known to be important for PTEN phosphatase activity and protein-protein interactions, such as Met35, Tyr155, Tyr240 and Tyr315. A PTEN redox interactome was also characterized, which identified a number of PTEN-interacting proteins that vary with the reversible inactivation of PTEN caused by H2O2 oxidation. These included new PTEN interactors as well as the redox proteins peroxiredoxin-1 (Prdx1) and thioredoxin (Trx), which are known to be involved in the recycling of PTEN active site following H2O2-induced reversible inactivation. The results suggest that the oxidative modification of PTEN causes functional alterations in PTEN structure and interactome, with fundamental implications for the PTEN signaling role in many cellular processes, such as those involved in the pathophysiology of disease and ageing.

Membrane protein extraction and purification using styrene-maleic acid (SMA) co-polymer:effect of variations in polymer structure

The use of styrene maleic acid (SMA) co-polymers to extract and purify transmembrane proteins, whilst retaining their native bilayer environment, overcomes many of the disadvantages associated with conventional detergent based procedures. This approach has huge potential for the future of membrane protein structural and functional studies. In this investigation we have systematically tested a range of commercially available SMA polymers, varying in both the ratio of styrene to maleic acid and in total size, for the ability to extract, purify and stabilise transmembrane proteins. Three different membrane proteins (BmrA, LeuT and ZipA) which vary in size and shape were used. Our results show that several polymers can be used to extract membrane proteins comparably to conventional detergents. A styrene:maleic acid ratio of either 2:1 or 3:1, combined with a relatively small average molecular weight (7.5-10 kDa) is optimal for membrane extraction, and this appears to be independent of the protein size, shape or expression system. A subset of polymers were taken forward for purification, functional and stability tests. Following a one-step affinity purification SMA 2000 was found to be the best choice for yield, purity and function. However the other polymers offer subtle differences in size and sensitivity to divalent cations that may be useful for a variety of downstream applications.

Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

Analysis of Molecular Dynamics Simulations of Protein Folding

Microsecond long Molecular Dynamics (MD) trajectories of biomolecular processes are now possible due to advances in computer technology. Soon, trajectories long enough to probe dynamics over many milliseconds will become available. Since these timescales match the physiological timescales over which many small proteins fold, all atom MD simulations of protein folding are now becoming popular. To distill features of such large folding trajectories, we must develop methods that can both compress trajectory data to enable visualization, and that can yield themselves to further analysis, such as the finding of collective coordinates and reduction of the dynamics. Conventionally, clustering has been the most popular MD trajectory analysis technique, followed by principal component analysis (PCA). Simple clustering used in MD trajectory analysis suffers from various serious drawbacks, namely, (i) it is not data driven, (ii) it is unstable to noise and change in cutoff parameters, and (iii) since it does not take into account interrelationships amongst data points, the separation of data into clusters can often be artificial. Usually, partitions generated by clustering techniques are validated visually, but such validation is not possible for MD trajectories of protein folding, as the underlying structural transitions are not well understood. Rigorous cluster validation techniques may be adapted, but it is more crucial to reduce the dimensions in which MD trajectories reside, while still preserving their salient features. PCA has often been used for dimension reduction and while it is computationally inexpensive, being a linear method, it does not achieve good data compression. In this thesis, I propose a different method, a nonmetric multidimensional scaling (nMDS) technique, which achieves superior data compression by virtue of being nonlinear, and also provides a clear insight into the structural processes underlying MD trajectories. I illustrate the capabilities of nMDS by analyzing three complete villin headpiece folding and six norleucine mutant (NLE) folding trajectories simulated by Freddolino and Schulten [1]. Using these trajectories, I make comparisons between nMDS, PCA and clustering to demonstrate the superiority of nMDS. The three villin headpiece trajectories showed great structural heterogeneity. Apart from a few trivial features like early formation of secondary structure, no commonalities between trajectories were found. There were no units of residues or atoms found moving in concert across the trajectories. A flipping transition, corresponding to the flipping of helix 1 relative to the plane formed by helices 2 and 3 was observed towards the end of the folding process in all trajectories, when nearly all native contacts had been formed. However, the transition occurred through a different series of steps in all trajectories, indicating that it may not be a common transition in villin folding. The trajectories showed competition between local structure formation/hydrophobic collapse and global structure formation in all trajectories. Our analysis on the NLE trajectories confirms the notion that a tight hydrophobic core inhibits correct 3-D rearrangement. Only one of the six NLE trajectories folded, and it showed no flipping transition. All the other trajectories get trapped in hydrophobically collapsed states. The NLE residues were found to be buried deeply into the core, compared to the corresponding lysines in the villin headpiece, thereby making the core tighter and harder to undo for 3-D rearrangement. Our results suggest that the NLE may not be a fast folder as experiments suggest. The tightness of the hydrophobic core may be a very important factor in the folding of larger proteins. It is likely that chaperones like GroEL act to undo the tight hydrophobic core of proteins, after most secondary structure elements have been formed, so that global rearrangement is easier. I conclude by presenting facts about chaperone-protein complexes and propose further directions for the study of protein folding.

Prediction of Phosphorylation Sites from Protein Sequences

Phosphorylation is amongst the most crucial and well-studied post-translational modifications. It is involved in multiple cellular processes which makes phosphorylation prediction vital for understanding protein functions. However, wet-lab techniques are labour and time intensive. Thus, computational tools are required for efficiency. This project aims to provide a novel way to predict phosphorylation sites from protein sequences by adding flexibility and Sezerman Grouping amino acid similarity measure to previous methods, as discovering new protein sequences happens at a greater rate than determining protein structures. The predictor – NOPAY - relies on Support Vector Machines (SVMs) for classification. The features include amino acid encoding, amino acid grouping, predicted secondary structure, predicted protein disorder, predicted protein flexibility, solvent accessibility, hydrophobicity and volume. As a result, we have managed to improve phosphorylation prediction accuracy for Homo sapiens by 3% and 6.1% for Mus musculus. Sensitivity at 99% specificity was also increased by 6% for Homo sapiens and for Mus musculus by 5% on independent test sets. In this study, we have managed to increase phosphorylation prediction accuracy for Homo sapiens and Mus musculus. When there is enough data, future versions of the software may also be able to predict other organisms.

A biochemical analysis into the co-translational folding of a G protein-coupled receptor; GPR35

The folding and targeting of membrane proteins poses a major challenge to the cell, as they must remain insertion competent while their highly hydrophobic transmembrane (TM) domains are transferred from the ribosome, through the aqueous cytosol and into the lipid bilayer. The biogenesis of a mature membrane protein takes place through the insertion and integration into the lipid bilayer. A number of TM proteins have been shown to gain some degree of secondary structure within the ribosome tunnel and to retain this conformation throughout maturation. Although studies into the folding and targeting of a number of membrane proteins have been carried out to date, there is little information on one of the largest class of eukaryotic membrane proteins; the G-protein-coupled receptors (GPCRs). This project studies the early folding events of the human ortholog of GPR35. To analyse the structure of the 1st TM domain, intermediates were generated and assessed by the biochemical method of pegylation (PEG-MAL). A structurally-similar microbial opsin (Bacterioopsin) was also used to investigate the differences in the early protein folding within eukaryotic and prokaryotic translation systems. Results showed that neither the 1st TM domain of GPR35 nor Bacterioopsin were capable of compacting in the ribosome tunnel before their N-terminus reached the ribosome exit point. The results for this assay remained consistent whether the proteins were translated in a eukaryotic or prokaryotic translation system. To examine the communication mechanism between the ribosome, the nascent chain and the protein targeting pathway, crosslinking experiments were carried out using the homobifunctional lysine cross-linker BS3. Specifically, the data generated here show that the nascent chain of GPR35 reaches the ribosomal protein uL23 in an extended conformation and interacts with the SRP protein as it exits the ribosome tunnel. This confirms the role of SRP in the co-translational targeting of GPR35. Using these methods insights into the early folding of GPCRs has been obtained. Further experiments using site-directed mutagenesis to reduce hydrophobicity in the 1st TM domain of GPR35, highlighted the mechanisms by which GPCRs are targeted to the endoplasmic reticulum. Confirming that hydrophobicity within the signal anchor sequence is essential of SRP-dependent targeting. Following the successful interaction of the nascent GPR35 and SRP, GPR35 is successfully targeted to ER membranes, shown here as dog pancreas microsomes (DPMs). Glycosylation of the GPR35 N-terminus was used to determine nascent chain structure as it is inserted into the ER membrane. These glycosylation experiments confirm that TM1 has obtained its compacted state whilst residing in the translocon. Finally, a site-specific cross-linking approach using the homobifunctional cysteine cross-linker, BMH, was used to study the lateral integration of GPR35 into the ER. Cross-linking of GPR35 TM1 and TM2 could be detected adjacent to a protein of ~45kDa, believed to be Sec61α. The loss of this adduct, as the nascent chain extends, showed the lateral movement of GPR35 TM1 from the translocon was dependent on the subsequent synthesis of TM2.

STRUCTURE AND FUNCTION OF THE GROUP III CHAPERONINS, A UNIQUE CLADE OF PROTEIN FOLDING NANOMACHINES

The survival and descent of cells is universally dependent on maintaining their proteins in a properly folded condition. It is widely accepted that the information for the folding of the nascent polypeptide chain into a native protein is encrypted in the amino acid sequence, and the Nobel Laureate Christian Anfinsen was the first to demonstrate that a protein could spontaneously refold after complete unfolding. However, it became clear that the observed folding rates for many proteins were much slower than rates estimated in vivo. This led to the recognition of required protein-protein interactions that promote proper folding. A unique group of proteins, the molecular chaperones, are responsible for maintaining protein homeostasis during normal growth as well as stress conditions. Chaperonins (CPNs) are ubiquitous and essential chaperones. They form ATP-dependent, hollow complexes that encapsulate polypeptides in two back-to-back stacked multisubunit rings, facilitating protein folding through highly cooperative allosteric articulation. CPNs are usually classified into Group I and Group II. Here, I report the characterization of a novel CPN belonging to a third Group, recently discovered in bacteria. Group III CPNs have close phylogenetic association to the Group II CPNs found in Archaea and Eukarya, and may be a relic of the Last Common Ancestor of the CPN family. The gene encoding the Group III CPN from Carboxydothermus hydrogenoformans and Candidatus Desulforudis audaxviator was cloned in E. coli and overexpressed in order to both characterize the protein and to demonstrate its ability to function as an ATPase chaperone. The opening and closing cycle of the Chy chaperonin was examined via site-directed mutations affecting the ATP binding site at R155. To relate the mutational analysis to the structure of the CPN, the crystal structure of both the AMP-PNP (an ATP analogue) and ADP bound forms were obtained in collaboration with Sun-Shin Cha in Seoul, South Korea. The ADP and ATP binding site substitutions resulted in frozen forms of the structures in open and closed conformations. From this, mutants were designed to validate hypotheses regarding key ATP interacting sites as well as important stabilizing interactions, and to observe the physical properties of the resulting complexes by calorimetry.

Umbelliferone induces changes in the structure and pharmacological activities of Bn IV, a phospholipase A(2) isoform isolated from Bothrops neuwiedi

In this paper was demonstrated that umbelliferone induces changes in structure and pharmacological activities of Bn IV, a lysine 49 secretory phospholipase A(2) (sPLA2) from Both tops neuwiedi. Incubation of Bn IV with umbelliferone virtually abolished platelet aggregation, edema, and myotoxicity induced by native Bn IV. The amino acid sequence of Bn IV showed high sequence similarities with other Lys49 sPLA2s from B. jararacussu (BthTx-I), B. pirajai (PrTx-I), and B. neuwiedi pauloensis (Bn SP6 and Bn SP7). This sPLA2 also has a highly conserved C-terminal amino acid sequence, which has been shown as important for the pharmacological activities of Lys49 sPLA2. Sequencing of Bn IV previously treated with umbelliferone revealed modification of S(1) and S(20). Fluorescent spectral analysis and circular dichroism (CD) studies showed that umbelliferone modified the secondary structure of this protein. Moreover, the pharmacological activity of Bn IV is driven by synergism of the C-terminal region with the a-helix motifs, which are involved in substrate binding of the Asp49 and Lys49 residues of 5PLA2 and have a direct effect on the Ca2+-independent membrane damage of some secretory snake venom PLA2. For Bn IV, these interactions are potentially important for triggering the pharmacological activity of this 5PLA2. (C) 2011 Elsevier Ltd. All rights reserved.