Structure-activity relationships for a new family of sulfonylurea herbicides

Acetohydroxyacid synthase (AHAS; EC 2.2.1.6) catalyzes the first common step in branched-chain amino acid biosynthesis. The enzyme is inhibited by several chemical classes of compounds and this inhibition is the basis of action of the sulfonylurea and imidazolinone herbicides. The commercial sulfonylureas contain a pyrimidine or a triazine ring that is substituted at both meta positions, thus obeying the initial rules proposed by Levitt. Here we assess the activity of 69 monosubstituted sulfonylurea analogs and related compounds as inhibitors of pure recombinant Arabidopsis thaliana AHAS and show that disubstitution is not absolutely essential as exemplified by our novel herbicide, monosulfuron (2-nitro-N-(4'-methyl-pyrimidin-2'-yl) phenyl-sulfonylurea), which has a pyrimidine ring with a single meta substituent. A subset of these compounds was tested for herbicidal activity and it was shown that their effect in vivo correlates well with their potency in vitro as AHAS inhibitors. Three-dimensional quantitative structure-activity relationships were developed using comparative molecular field analysis and comparative molecular similarity indices analysis. For the latter, the best result was obtained when steric, electrostatic, hydrophobic and H-bond acceptor factors were taken into consideration. The resulting fields were mapped on to the published crystal structure of the yeast enzyme and it was shown that the steric and hydrophobic fields are in good agreement with sulfonylurea-AHAS interaction geometry.

Herbicide-binding Sites Revealed in the Structure of Plant Acetohydroxyacid Synthase

The sulfonylureas and imidazolinones are potent commercial herbicide families. They are among the most popular choices for farmers worldwide, because they are nontoxic to animals and highly selective. These herbicides inhibit branched-chain amino acid biosynthesis in plants by targeting acetohydroxyacid synthase (AHAS, EC 2.2.1.6). This report describes the 3D structure of Arabidopsis thaliana AHAS in complex with five sulfonylureas (to 2.5 angstrom resolution) and with the imidazolinone, imazaquin (IQ; 2.8 angstrom). Neither class of molecule has a structure that mimics the substrates for the enzyme, but both inhibit by blocking a channel through which access to the active site is gained. The sulfonylureas approach within 5 angstrom of the catalytic center, which is the C2 atom of the cofactor thiamin diphosphate, whereas IQ is at least 7 angstrom from this atom. Ten of the amino acid residues that bind the sulfonylureas also bind IQ. Six additional residues interact only with the sulfonylureas, whereas there are two residues that bind IQ but not the sulfonylureas. Thus, the two classes of inhibitor occupy partially overlapping sites but adopt different modes of binding. The increasing emergence of resistant weeds due to the appearance of mutations that interfere with the inhibition of AHAS is now a worldwide problem. The structures described here provide a rational molecular basis for understanding these mutations, thus allowing more sophisticated AHAS inhibitors to be developed. There is no previously described structure for any plant protein in complex with a commercial herbicide.

Antimicrobial peptides in action

Molecular dynamics simulations of the magainin MG-H2 peptide interacting with a model phospholipid membrane have been used to investigate the mechanism by which antimicrobial peptides act. Multiple copies of the peptide were randomly placed in solution close to the membrane. The peptide readily bound to the membrane, and above a certain concentration, the peptide was observed to cooperatively induce the formation of a nanometer- sized, toroidally shaped pore in the bilayer. In sharp contrast with the commonly accepted model of a toroidal pore, only one peptide was typically found near the center of the pore. The remaining peptides lay close to the edge of the pore, maintaining a predominantly parallel orientation with respect to the membrane.

Crystallization and preliminary X-ray analysis of a Kunitz-type inhibitor, textilinin-1 from Pseudonaja textilis textilis

Textilinin-1 (Txln-1), a Kunitz-type serine protease inhibitor, is a 59-amino-acid polypeptide isolated from the venom of the Australian Common Brown snake Pseudonaja textilis textilis. This molecule has been suggested as an alternative to aprotinin, also a Kunitz-type serine protease inhibitor, for use as an anti-bleeding agent in surgical procedures. Txln-1 shares only 47% amino-acid identity to aprotinin; however, six cysteine residues in the two peptides are in conserved locations. It is therefore expected that the overall fold of these molecules is similar but that they have contrasting surface features. Here, the crystallization of recombinant textilinin-1 (rTxln-1) as the free molecule and in complex with bovine trypsin (229 amino acids) is reported. Two organic solvents, phenol and 1,4-butanediol, were used as additives to facilitate the crystallization of free rTxln-1. Crystals of the rTxln-1-bovine trypsin complex diffracted to 2.0 angstrom resolution, while crystals of free rTxln-1 diffracted to 1.63 angstrom resolution.

Identification and molecular modeling of a novel, plant-like, human purple acid phosphatase

Purple acid phosphatases are a family of binuclear metallohydrolases that have been identified in plants, animals and fungi. Only one isoform of similar to 35 kDa has been isolated from animals, where it is associated with bone resorption and microbial killing through its phosphatase activity, and hydroxyl radical production, respectively. Using the sensitive PSI-BLAST search method, sequences representing new purple acid phosphatase-like proteins have been identified in mammals, insects and nematodes. These new putative isoforms are closely related to the similar to 55 kDa purple acid phosphatase characterized from plants. Secondary structure prediction of the new human isoform further confirms its similarity to a purple acid phosphatase from the red kidney bean. A structural model for the human enzyme was constructed based on the red kidney bean purple acid phosphatase structure. This model shows that the catalytic centre observed in other purple acid phosphatases is also present in this new isoform. These observations suggest that the sequences identified in this study represent a novel subfamily of plant-like purple acid phosphatases in animals and humans. (c) 2006 Elsevier B.V. All rights reserved.

Molecular dynamics simulations of the hydrophobin SC3 at a hydrophobic/hydrophilic interface

Hydrophobins are small (similar to 100 aa) proteins that have an important role in the growth and development of mycelial fungi. They are surface active and, after secretion by the fungi, self-assemble into amphipathic membranes at hydrophobic/hydrophilic interfaces, reversing the hydrophobicity of the surface. In this study, molecular dynamics simulation techniques have been used to model the process by which a specific class I hydrophobin, SC3, binds to a range of hydrophobic/ hydrophilic interfaces. The structure of SC3 used in this investigation was modeled based on the crystal structure of the class II hydrophobin HFBII using the assumption that the disulfide pairings of the eight conserved cysteine residues are maintained. The proposed model for SC3 in aqueous solution is compact and globular containing primarily P-strand and coil structures. The behavior of this model of SC3 was investigated at an air/water, an oil/water, and a hydrophobic solid/water interface. It was found that SC3 preferentially binds to the interfaces via the loop region between the third and fourth cysteine residues and that binding is associated with an increase in a-helix formation in qualitative agreement with experiment. Based on a combination of the available experiment data and the current simulation studies, we propose a possible model for SC3 self-assembly on a hydrophobic solid/water interface.

Determination of the orientational distribution and orientation factor for transfer between membrane-bound fluorophores using a confocal microscope

Orientational fluorophores have been a useful tool in physical chemistry, biochemistry, and more recently structural biology due to the polarized nature of the light they emit and that fact that energy can be transferred between them. We present a practical scheme in which measurements of the intensity of emitted fluorescence can be used to determine limits on the mean and distribution of orientation of the absorption transition moment of membrane-bound. uorophores. We demonstrate how information about the orientation of. uorophores can be used to calculate the orientation factor k(2) required for use in FRET spectroscopy. We illustrate the method using images of AlexaFluor probes bound to MscL mechanosensitive transmembrane channel proteins in spherical liposomes.

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.

Peptide binding prediction for the human class II MHC allele HLA-DP2:a molecular docking approach

MHC class II proteins bind oligopeptide fragments derived from proteolysis of pathogen antigens, presenting them at the cell surface for recognition by CD4+ T cells. Human MHC class II alleles are grouped into three loci: HLA-DP, HLA-DQ and HLA-DR. In contrast to HLA-DR and HLA-DQ, HLA-DP proteins have not been studied extensively, as they have been viewed as less important in immune responses than DRs and DQs. However, it is now known that HLA-DP alleles are associated with many autoimmune diseases. Quite recently, the X-ray structure of the HLA-DP2 molecule (DPA*0103, DPB1*0201) in complex with a self-peptide derived from the HLA-DR a-chain has been determined. In the present study, we applied a validated molecular docking protocol to a library of 247 modelled peptide-DP2 complexes, seeking to assess the contribution made by each of the 20 naturally occurred amino acids at each of the nine binding core peptide positions and the four flanking residues (two on both sides).

Relieving the first bottleneck in the drug discovery pipeline:using array technologies to rationalize membrane protein production

The slow down in the drug discovery pipeline is, in part, owing to a lack of structural and functional information available for new drug targets. Membrane proteins, the targets of well over 50% of marketed pharmaceuticals, present a particular challenge. As they are not naturally abundant, they must be produced recombinantly for the structural biology that is a prerequisite to structure-based drug design. Unfortunately, however, obtaining high yields of functional, recombinant membrane proteins remains a major bottleneck in contemporary bioscience. While repeated rounds of trial-and-error optimization have not (and cannot) reveal mechanistic details of the biology of recombinant protein production, examination of the host response has provided new insights. To this end, we published an early transcriptome analysis that identified genes implicated in high-yielding yeast cell factories, which has enabled the engineering of improved production strains. These advances offer hope that the bottleneck of membrane protein production can be relieved rationally.

Functional recombinant protein is present in the pre-induction phases of Pichia pastoris cultures when grown in bioreactors, but not shake-flasks

Background - Pichia pastoris is a widely-used host for recombinant protein production; expression is typically driven by methanol-inducible alcohol oxidase (AOX) promoters. Recently this system has become an important source of recombinant G protein-coupled receptors (GPCRs) for structural biology and drug discovery. The influence of diverse culture parameters (such as pH, dissolved oxygen concentration, medium composition, antifoam concentration and culture temperature) on productivity has been investigated for a wide range of recombinant proteins in P. pastoris. In contrast, the impact of the pre-induction phases on yield has not been as closely studied. In this study, we examined the pre-induction phases of P. pastoris bioreactor cultivations producing three different recombinant proteins: the GPCR, human A2a adenosine receptor (hA2aR), green fluorescent protein (GFP) and human calcitonin gene-related peptide receptor component protein (as a GFP fusion protein; hCGRP-RCP-GFP). Results - Functional hA2aR was detected in the pre-induction phases of a 1 L bioreactor cultivation of glycerol-grown P. pastoris. In a separate experiment, a glycerol-grown P. pastoris strain secreted soluble GFP prior to methanol addition. When glucose, which has been shown to repress AOX expression, was the pre-induction carbon source, hA2aR and GFP were still produced in the pre-induction phases. Both hA2aR and GFP were also produced in methanol-free cultivations; functional protein yields were maintained or increased after depletion of the carbon source. Analysis of the pre-induction phases of 10 L pilot scale cultivations also demonstrated that pre-induction yields were at least maintained after methanol induction, even in the presence of cytotoxic concentrations of methanol. Additional bioreactor data for hCGRP-RCP-GFP and shake-flask data for GFP, horseradish peroxidase (HRP), the human tetraspanins hCD81 and CD82, and the tight-junction protein human claudin-1, demonstrated that bioreactor but not shake flask cultivations exhibit recombinant protein production in the pre-induction phases of P. pastoris cultures. Conclusions - The production of recombinant hA2aR, GFP and hCGRP-RCP-GFP can be detected in bioreactor cultivations prior to methanol induction, while this is not the case for shake-flask cultivations of GFP, HRP, hCD81, hCD82 and human claudin-1. This confirms earlier suggestions of leaky expression from AOX promoters, which we report here for both glycerol- and glucose-grown cells in bioreactor cultivations. These findings suggest that the productivity of AOX-dependent bioprocesses is not solely dependent on induction by methanol. We conclude that in order to maximize total yields, pre-induction phase cultivation conditions should be optimized, and that increased specific productivity may result in decreased biomass yields.

Statistical deconvolution of enthalpic energetic contributions to MHC-peptide binding affinity

Background - MHC Class I molecules present antigenic peptides to cytotoxic T cells, which forms an integral part of the adaptive immune response. Peptides are bound within a groove formed by the MHC heavy chain. Previous approaches to MHC Class I-peptide binding prediction have largely concentrated on the peptide anchor residues located at the P2 and C-terminus positions. Results - A large dataset comprising MHC-peptide structural complexes was created by re-modelling pre-determined x-ray crystallographic structures. Static energetic analysis, following energy minimisation, was performed on the dataset in order to characterise interactions between bound peptides and the MHC Class I molecule, partitioning the interactions within the groove into van der Waals, electrostatic and total non-bonded energy contributions. Conclusion - The QSAR techniques of Genetic Function Approximation (GFA) and Genetic Partial Least Squares (G/PLS) algorithms were used to identify key interactions between the two molecules by comparing the calculated energy values with experimentally-determined BL50 data. Although the peptide termini binding interactions help ensure the stability of the MHC Class I-peptide complex, the central region of the peptide is also important in defining the specificity of the interaction. As thermodynamic studies indicate that peptide association and dissociation may be driven entropically, it may be necessary to incorporate entropic contributions into future calculations.

Improving T cell-induced response to subunit vaccines:opportunities for a proteomic systems approach

Prophylactic vaccines are an effective strategy to prevent development of many infectious diseases. With new and re-emerging infections posing increasing risks to food stocks and the health of the population in general, there is a need to improve the rationale of vaccine development. One key challenge lies in development of an effective T cell-induced response to subunit vaccines at specific sites and in different populations. Objectives: In this review, we consider how a proteomic systems-based approach can be used to identify putative novel vaccine targets, may be adopted to characterise subunit vaccines and adjuvants fully. Key findings: Despite the extensive potential for proteomics to aid our understanding of subunit vaccine nature, little work has been reported on identifying MHC 1-binding peptides for subunit vaccines generating T cell responses in the literature to date. Summary: In combination with predictive and structural biology approaches to mapping antigen presentation, proteomics offers a powerful and as yet un-tapped addition to the armoury of vaccine discovery to predict T-cell subset responses and improve vaccine design strategies.

Peptide binding to HLA-DP proteins at pH 5.0 and pH 7.0:a quantitative molecular docking study

Background: HLA-DPs are class II MHC proteins mediating immune responses to many diseases. Peptides bind MHC class II proteins in the acidic environment within endosomes. Acidic pH markedly elevates association rate constants but dissociation rates are almost unchanged in the pH range 5.0 - 7.0. This pH-driven effect can be explained by the protonation/deprotonation states of Histidine, whose imidazole has a pKa of 6.0. At pH 5.0, imidazole ring is protonated, making Histidine positively charged and very hydrophilic, while at pH 7.0 imidazole is unprotonated, making Histidine less hydrophilic. We develop here a method to predict peptide binding to the four most frequent HLA-DP proteins: DP1, DP41, DP42 and DP5, using a molecular docking protocol. Dockings to virtual combinatorial peptide libraries were performed at pH 5.0 and pH 7.0. Results: The X-ray structure of the peptide - HLA-DP2 protein complex was used as a starting template to model by homology the structure of the four DP proteins. The resulting models were used to produce virtual combinatorial peptide libraries constructed using the single amino acid substitution (SAAS) principle. Peptides were docked into the DP binding site using AutoDock at pH 5.0 and pH 7.0. The resulting scores were normalized and used to generate Docking Score-based Quantitative Matrices (DS-QMs). The predictive ability of these QMs was tested using an external test set of 484 known DP binders. They were also compared to existing servers for DP binding prediction. The models derived at pH 5.0 predict better than those derived at pH 7.0 and showed significantly improved predictions for three of the four DP proteins, when compared to the existing servers. They are able to recognize 50% of the known binders in the top 5% of predicted peptides. Conclusions: The higher predictive ability of DS-QMs derived at pH 5.0 may be rationalised by the additional hydrogen bond formed between the backbone carbonyl oxygen belonging to the peptide position before p1 (p-1) and the protonated ε-nitrogen of His 79β. Additionally, protonated His residues are well accepted at most of the peptide binding core positions which is in a good agreement with the overall negatively charged peptide binding site of most MHC proteins. © 2012 Patronov et al.; licensee BioMed Central Ltd.

The effects of manipulation of sedimentary iron and organic matter on sediment biogeochemistry and seagrasses in a subtropical carbonate environment

The microbial metabolism of organic matter (OM) in seagrass beds can create sulfidic conditions detrimental to seagrass growth; iron (Fe) potentially has ameliorating effects through titration of the sulfides and the precipitation of iron-sulfide minerals into the sediment. In this study, the biogeochemical effects of Fe availability and its interplay with sulfur and OM on sulfide toxicity, phosphorous (P) availability, seagrass growth and community structure were tested. The availability of Fe and OM was manipulated in a 2 × 2 factorial experiment arranged in a Latin square, with four replicates per treatment. The treatments included the addition of Fe, the addition of OM, the addition of both Fe and OM as well as no addition. The experiment was conducted in an oligotrophic, iron-deficient seagrass bed. Fe had an 84.5% retention efficiency in the sediments with the concentration of Fe increasing in the seagrass leaves over the course of the experiment. Porewater chemistry was significantly altered with a dramatic decrease in sulfide levels in Fe addition plots while sulfide levels increased in the OM addition treatments. Phosphorus increased in seagrass leaves collected in the Fe addition plots. Decreased sulfide stress was evidenced by heavier δ34S in leaves and rhizomes from plots to which Fe was added. The OM addition negatively affected seagrass growth but increased P availability; the reduced sulfide stress in Fe added plots resulted in elevated productivity. Fe availability may be an important determinant of the impact that OM has on seagrass vitality in carbonate sediments vegetated with seagrasses.